05 load raster
Loading local raster datasets with leafmap
Uncomment the following line to install leafmap if needed.
# !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.
import os
import leafmap.leafmap as leafmap
Specify input raster datasets
landsat = "landsat.tif"
dem = "dem.tif"
Download samples raster datasets
More datasets can be downloaded from https://viewer.nationalmap.gov/basic/
url1 = "https://open.gishub.org/data/raster/landsat7.tif"
url2 = "https://open.gishub.org/data/raster/srtm90.tif"
satellite = leafmap.download_file(url1, "landsat.tif")
dem = leafmap.download_file(url2, "dem.tif")
Downloading... From: https://open.gishub.org/data/raster/landsat7.tif To: /home/runner/work/leafmap/leafmap/docs/notebooks/landsat.tif
0%| | 0.00/4.87M [00:00<?, ?B/s]
100%|██████████| 4.87M/4.87M [00:00<00:00, 83.4MB/s]
Downloading... From: https://open.gishub.org/data/raster/srtm90.tif To: /home/runner/work/leafmap/leafmap/docs/notebooks/dem.tif
0%| | 0.00/16.6M [00:00<?, ?B/s]
100%|██████████| 16.6M/16.6M [00:00<00:00, 184MB/s]
Create an interactive map
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
m.add_raster(dem, colormap="terrain", layer_name="DEM")
m.add_raster(landsat, bands=[5, 4, 3], layer_name="Landsat")
Display the map
m