grid['node']['Elevation'] = elev ids = grid.node_at_cell # Burroughs slope_Burrough, aspect_Burrough = \ grid.calculate_slope_aspect_at_nodes_Burrough(ids, vals='Elevation') # Horn slope_horn, aspect_horn = \ grid.calculate_slope_aspect_at_nodes_horn(ids, vals='Elevation') # BestFitPlane slope_BFP, aspect_BFP = grid.calculate_slope_aspect_at_nodes_bestFitPlane( ids, elev) # node_slopes_using_patches slope_NSP = grid.node_slopes_using_patches(elevs='Elevation', unit='degrees') pic = 0 plt.figure(pic) imshow_field(grid, 'Elevation', values_at='node', grid_units=('m', 'm')) plt.title('Elevation in m') # plt.savefig('Elevation_NS') pic += 1 plt.figure(pic) imshow_grid(grid, np.degrees(slope_Burrough), values_at='cell', grid_units=('m', 'm')) plt.title('Slope in degrees - Burrough 1998') # plt.savefig('Slope_Burrough_NS') pic += 1
grid = RasterModelGrid(53, 67, 10.0) elev = np.load("elevation_NS.npy") grid["node"]["Elevation"] = elev ids = grid.node_at_cell # Burroughs slope_burrough, aspect_burrough = grid.calculate_slope_aspect_at_nodes_burrough(ids, vals="Elevation") # Horn slope_horn, aspect_horn = grid.calculate_slope_aspect_at_nodes_horn(ids, vals="Elevation") # BestFitPlane slope_bfp, aspect_bfp = grid.calculate_slope_aspect_at_nodes_best_fit_plane(ids, elev) # node_slopes_using_patches slope_NSP = grid.node_slopes_using_patches(elevs="Elevation", unit="degrees") pic = 0 plt.figure(pic) imshow_field(grid, "Elevation", values_at="node", grid_units=("m", "m")) plt.title("Elevation in m") # plt.savefig('Elevation_NS') pic += 1 plt.figure(pic) imshow_grid(grid, np.degrees(slope_burrough), values_at="cell", grid_units=("m", "m")) plt.title("Slope in degrees - Burrough 1998") # plt.savefig('Slope_burrough_NS') pic += 1 plt.figure(pic)