""" if __name__ == "__main__": from pydem.test_pydem import get_test_data from pydem.dem_processing import DEMProcessor from pydem.test_pydem import make_file_names import sys if len(sys.argv) <= 1: testnum = 30 # Default test number NN = 64 # Resolution of tile # Get test data and create DEMProcessor elev, ang, test_data = get_test_data(testnum, NN) filename = make_file_names(testnum, NN)['fel'] dem_proc = DEMProcessor(filename, dx_dy_from_file=False) # Calculate magnitude and direction of slopes mag, direction = dem_proc.calc_slopes_directions() # Calculate UCA/TWI on entire tile area1 = dem_proc.calc_uca(0) twi1 = dem_proc.calc_twi() # Calculate UCA/TWI on chunks dem_proc.resolve_edges = True # If false, won't do edge resolution dem_proc.chunk_size_uca = 9 # Chunk size dem_proc.chunk_overlap_uca = 3 #Overlap between chunks area2 = dem_proc.calc_uca(0) twi2 = dem_proc.calc_twi()
) parser.add_argument( "--variable-coords", "--vc", help="if set, will create a non-" "uniform grid on the elevation", action="store_true", ) args = parser.parse_args() testnum = args.testnum dx_dy_from_file = args.variable_coords testnum += dx_dy_from_file NN = 511 # The resolution of the grid (N x N) if testnum >= 28: filename = make_file_names(testnum, NN)["elev"] else: filename = make_file_names(testnum, NN)["fel"] elev, ang, test_data = get_test_data(testnum, NN) # Make the filenames of the results filename2 = make_file_names(testnum, NN)["uca"] if os.path.exists(filename2): uca_file = GdalReader(file_name=filename2) uca, = uca_file.raster_layers # Make the pydem Processor object and process the elevation dem_proc = DEMProcessor(filename, dx_dy_from_file=dx_dy_from_file) mag, direction = dem_proc.calc_slopes_directions() area = dem_proc.calc_uca(0)
default=30, nargs='?') parser.add_argument('--variable-coords', '--vc', help='if set, will create a non-' 'uniform grid on the elevation', action="store_true") args = parser.parse_args() testnum = args.testnum dx_dy_from_file = args.variable_coords testnum += dx_dy_from_file NN = 64 # The resolution of the grid (N x N) if testnum >= 28: filename = make_file_names(testnum, NN)['elev'] else: filename = make_file_names(testnum, NN)['fel'] elev, ang, test_data = get_test_data(testnum, NN) # Make the filenames of the results filename2 = make_file_names(testnum, NN)['uca'] if os.path.exists(filename2): uca_file = GdalReader(file_name=filename2) uca, = uca_file.raster_layers # Make the pydem Processor object and process the elevation dem_proc = DEMProcessor(filename, dx_dy_from_file=dx_dy_from_file) mag, direction = dem_proc.calc_slopes_directions() area = dem_proc.calc_uca(0)
limitations under the License. """ if __name__ == "__main__": from pydem.test_pydem import get_test_data from pydem.dem_processing import DEMProcessor from pydem.test_pydem import make_file_names import sys if len(sys.argv) <= 1: testnum = 30 # Default test number NN = 64 # Resolution of tile # Get test data and create DEMProcessor elev, ang, test_data = get_test_data(testnum, NN) filename = make_file_names(testnum, NN)['fel'] dem_proc = DEMProcessor(filename, dx_dy_from_file=False) # Calculate magnitude and direction of slopes mag, direction = dem_proc.calc_slopes_directions() # Calculate UCA/TWI on entire tile area1 = dem_proc.calc_uca(0) twi1 = dem_proc.calc_twi() # Calculate UCA/TWI on chunks dem_proc.resolve_edges = True # If false, won't do edge resolution dem_proc.chunk_size_uca = 9 # Chunk size dem_proc.chunk_overlap_uca = 3 #Overlap between chunks area2 = dem_proc.calc_uca(0) twi2 = dem_proc.calc_twi()
description='Compare pydem output to taudem output') parser.add_argument('testnum', help='the number of the test elevation to ' 'process', type=int, default=30, nargs='?') parser.add_argument('--variable-coords', '--vc', help='if set, will create a non-' 'uniform grid on the elevation', action="store_true") args = parser.parse_args() testnum = args.testnum dx_dy_from_file = args.variable_coords testnum += dx_dy_from_file NN = 64 # The resolution of the grid (N x N) if testnum >= 28: filename = make_file_names(testnum, NN)['elev'] else: filename = make_file_names(testnum, NN)['fel'] elev, ang, test_data = get_test_data(testnum, NN) # Make the filenames of the results filename2 = make_file_names(testnum, NN)['uca'] if os.path.exists(filename2): uca_file = GdalReader(file_name=filename2) uca, = uca_file.raster_layers # Make the pydem Processor object and process the elevation dem_proc = DEMProcessor(filename, dx_dy_from_file=dx_dy_from_file) mag, direction = dem_proc.calc_slopes_directions() area = dem_proc.calc_uca(0)