def fill_and_shade(located_elevations):
"""
Given an array of elevations in a LocatedImage, fill any holes in
the data and add a relief (shadows) to give a realistic 3d appearance.
"""
new_elevations = srtm.fill_gaps(located_elevations.image, max_distance=15)
new_img = srtm.add_shading(new_elevations, azimuth=135, altitude=15)
return LocatedImage(new_img, located_elevations.extent)
def main():
ax = plt.axes(projection=ccrs.PlateCarree())
# Get the 1x1 degree SRTM tile for 12E, 47N
elev, crs, extent = srtm.srtm_composite(12, 47, 1, 1)
# Fill the gaps present in the elevation data
elev_filled = srtm.fill_gaps(elev, 15)
# Add shading simulating the Sun at 10am (South-East)
# and with a low angle (15 degrees above horizon)
shaded = srtm.add_shading(elev_filled, 135.0, 15.0)
# The plot the result :
plt.imshow(shaded, extent=extent, transform=crs,
cmap='Greys', origin='lower')
plt.title("SRTM Shaded Relief Map")
gl = ax.gridlines(draw_labels=True)
gl.xlabels_top = False
gl.ylabels_left = False
plt.show()
def main():
ax = plt.axes(projection=ccrs.PlateCarree())
# Get the 1x1 degree SRTM tile for 12E, 47N
elev, crs, extent = srtm.srtm_composite(12, 47, 1, 1)
# Fill the gaps present in the elevation data
elev_filled = srtm.fill_gaps(elev, 15)
# Add shading simulating the Sun at 10am (South-East)
# and with a low angle (15 degrees above horizon)
shaded = srtm.add_shading(elev_filled, 135.0, 15.0)
# The plot the result :
plt.imshow(shaded, extent=extent, transform=crs,
cmap='Greys', origin='lower')
plt.title("SRTM Shaded Relief Map")
gl = ax.gridlines(draw_labels=True,)
gl.xlabels_top = False
gl.ylabels_left = False
plt.show()
