Creating a shaded relief map by blending DEM and hillshade

Uncomment the following line to install geemap if needed.

In [ ]:

# !pip install geemap

Import libraries

In [ ]:

import ee
import geemap
import geemap.colormaps as cm

Create an interactive map

In [ ]:

Map = geemap.Map()
Map

Add DEM and hillshade to the map

In [ ]:

dem = ee.Image("CGIAR/SRTM90_V4")
hillsahde = ee.Terrain.hillshade(dem)

In [ ]:

vis = {'min': 0, 'max': 6000, 'palette': cm.palettes.dem}

In [ ]:

Map.addLayer(hillsahde, {}, 'Hillshade')
Map.addLayer(dem, vis, 'DEM')

Create a blended image by blending DEM and hillshade

In [ ]:

blend = geemap.blend(top_layer=dem, top_vis=vis)

In [ ]:

Map.addLayer(blend, {}, 'Blend')

Add NLCD land cover to the map

In [ ]:

nlcd = ee.Image("USGS/NLCD_RELEASES/2019_REL/NLCD/2019").select('landcover')
nlcd_vis = {'bands': ['landcover']}
Map.addLayer(nlcd, nlcd_vis, 'NLCD')

Create a blended image by blending NLCD and DEM.

In [ ]:

result = geemap.blend(nlcd, dem, top_vis=nlcd_vis, expression='a*b')

In [ ]:

Map.addLayer(result, {}, 'Blend NLCD')