Python RasterModelGrid.node_slopes_using_patches Exemples

Exemple #1

Fichier : slope_aspect_routines_comparison.py Projet : emielkater/landlab

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

Exemple #2

Fichier : slope_aspect_routines_comparison.py Projet : jennyknuth/landlab

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)