Add Google Earth Engine data layers

This notebook demonstrates how to add Google Earth Engine data layers to a map.

Uncomment the following line to install leafmap if needed.

In [1]:

# %pip install "leafmap[maplibre]" geemap

# %pip install "leafmap[maplibre]" geemap

In [2]:

import leafmap.maplibregl as leafmap

import leafmap.maplibregl as leafmap

To run this notebook, you will need an API key from MapTiler. Once you have the API key, you can uncomment the following code block and replace YOUR_API_KEY with your actual API key. Then, run the code block code to set the API key as an environment variable.

In [3]:

# import os
# os.environ["MAPTILER_KEY"] = "YOUR_API_KEY"

# import os # os.environ["MAPTILER_KEY"] = "YOUR_API_KEY"

You can find a list of Earth Engine asset ids from here, which does not require an Earth Engine account.

In [4]:

m = leafmap.Map(
    center=[-120.4482, 38.0399], zoom=13, pitch=60, bearing=30, style="3d-terrain"
)
m.add_ee_layer(asset_id="ESA/WorldCover/v200", opacity=0.5)
m.add_legend(builtin_legend="ESA_WorldCover", title="ESA Landcover")
m.add_layer_control()
m

m = leafmap.Map( center=[-120.4482, 38.0399], zoom=13, pitch=60, bearing=30, style="3d-terrain" ) m.add_ee_layer(asset_id="ESA/WorldCover/v200", opacity=0.5) m.add_legend(builtin_legend="ESA_WorldCover", title="ESA Landcover") m.add_layer_control() m

An API key is required to use the 3D terrain feature.

In [5]:

m.layer_interact()

m.layer_interact()

We can also overlay other data layers on top of Earth Engine data layers.

In [6]:

m = leafmap.Map(
    center=[-74.012998, 40.70414], zoom=15.6, pitch=60, bearing=30, style="3d-terrain"
)
m.add_ee_layer(asset_id="ESA/WorldCover/v200", opacity=0.5)
m.add_legend(builtin_legend="ESA_WorldCover", title="ESA Landcover")
m.add_layer_control()
MAPTILER_KEY = leafmap.get_api_key("MAPTILER_KEY")
source = {
    "url": f"https://api.maptiler.com/tiles/v3/tiles.json?key={MAPTILER_KEY}",
    "type": "vector",
}

layer = {
    "id": "3d-buildings",
    "source": "openmaptiles",
    "source-layer": "building",
    "type": "fill-extrusion",
    "min-zoom": 15,
    "paint": {
        "fill-extrusion-color": [
            "interpolate",
            ["linear"],
            ["get", "render_height"],
            0,
            "lightgray",
            200,
            "royalblue",
            400,
            "lightblue",
        ],
        "fill-extrusion-height": [
            "interpolate",
            ["linear"],
            ["zoom"],
            15,
            0,
            16,
            ["get", "render_height"],
        ],
        "fill-extrusion-base": [
            "case",
            [">=", ["get", "zoom"], 16],
            ["get", "render_min_height"],
            0,
        ],
    },
}
m.add_source("openmaptiles", source)
m.add_layer(layer)
m

m = leafmap.Map( center=[-74.012998, 40.70414], zoom=15.6, pitch=60, bearing=30, style="3d-terrain" ) m.add_ee_layer(asset_id="ESA/WorldCover/v200", opacity=0.5) m.add_legend(builtin_legend="ESA_WorldCover", title="ESA Landcover") m.add_layer_control() MAPTILER_KEY = leafmap.get_api_key("MAPTILER_KEY") source = { "url": f"https://api.maptiler.com/tiles/v3/tiles.json?key={MAPTILER_KEY}", "type": "vector", } layer = { "id": "3d-buildings", "source": "openmaptiles", "source-layer": "building", "type": "fill-extrusion", "min-zoom": 15, "paint": { "fill-extrusion-color": [ "interpolate", ["linear"], ["get", "render_height"], 0, "lightgray", 200, "royalblue", 400, "lightblue", ], "fill-extrusion-height": [ "interpolate", ["linear"], ["zoom"], 15, 0, 16, ["get", "render_height"], ], "fill-extrusion-base": [ "case", [">=", ["get", "zoom"], 16], ["get", "render_min_height"], 0, ], }, } m.add_source("openmaptiles", source) m.add_layer(layer) m

An API key is required to use the 3D terrain feature.

If you have an Earth Engine, you can uncomment the first two code blocks to add any Earth Engine datasets.

In [7]:

# import ee
# ee.Initialize(project="YOUR-PROJECT-ID")

# import ee # ee.Initialize(project="YOUR-PROJECT-ID")

In [8]:

# m = leafmap.Map(center=[-120.4482, 38.03994], zoom=13, pitch=60, bearing=30, style="3d-terrain")
# dataset = ee.ImageCollection("ESA/WorldCover/v200").first()
# vis_params = {"bands": ["Map"]}
# m.add_ee_layer(dataset, vis_params, name="ESA Worldcover", opacity=0.5)
# m.add_legend(builtin_legend="ESA_WorldCover", title="ESA Landcover")
# m.add_layer_control()
# m

# m = leafmap.Map(center=[-120.4482, 38.03994], zoom=13, pitch=60, bearing=30, style="3d-terrain") # dataset = ee.ImageCollection("ESA/WorldCover/v200").first() # vis_params = {"bands": ["Map"]} # m.add_ee_layer(dataset, vis_params, name="ESA Worldcover", opacity=0.5) # m.add_legend(builtin_legend="ESA_WorldCover", title="ESA Landcover") # m.add_layer_control() # m