102 blend hillshade

Creating a shaded relief map by blending DEM and hillshade

Uncomment the following line to install geemap if needed.

In [ ]:

# !pip install geemap

# !pip install geemap

Import libraries

In [ ]:

import ee
import geemap
import geemap.colormaps as cm

import ee import geemap import geemap.colormaps as cm

Create an interactive map

In [ ]:

Map = geemap.Map()
Map

Map = geemap.Map() Map

Add DEM and hillshade to the map

In [ ]:

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

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

In [ ]:

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

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

In [ ]:

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

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)

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

In [ ]:

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

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')

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')

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

In [ ]:

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

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

---

Last update: 2023-03-19
Created: 2021-12-18