Skip to content

05 load raster

image image image

Loading local raster datasets with leafmap

Uncomment the following line to install leafmap if needed.

1
# !pip install leafmap

To follow this tutorial, you need to install the leafmap and xarray_leaflet Python packages. Use the following conda commands to create a conda env and install packages. Note that xarray_leaflet does not work properly on Windows (source). Also, the add_raster function is only available for the ipyleaflet plotting backend. Therefore, Google Colab is not supported. Use the binder link above instead.

  • conda create -n gee python
  • conda activate gee
  • conda install mamba -c conda-forge
  • mamba install leafmap xarray_leaflet -c conda-forge

Use the ipyleaflet plotting backend. The folium plotting backend does not support the add_raster function.

1
2
import os
import leafmap.leafmap as leafmap

Specify input raster datasets

1
2
landsat = "landsat.tif"
dem = "dem.tif"

Download samples raster datasets

More datasets can be downloaded from https://viewer.nationalmap.gov/basic/

1
2
3
4
url1 = "https://opengeos.org/data/raster/landsat7.tif"
url2 = "https://opengeos.org/data/raster/srtm90.tif"
satellite = leafmap.download_file(url1, "landsat.tif")
dem = leafmap.download_file(url2, "dem.tif")

Create an interactive map

1
m = leafmap.Map()

Add local raster datasets to the map

More colormap can be found at https://matplotlib.org/3.1.0/tutorials/colors/colormaps.html

1
m.add_raster(dem, colormap="terrain", layer_name="DEM")

1
m.add_raster(landsat, bands=[5, 4, 3], layer_name="Landsat")

Display the map

1
m