138 draw control

Clip Earth Engine images interactively with the Draw Control

Uncomment the following line to install geemap if needed.

In [ ]:

# %pip install -U geemap

# %pip install -U geemap

In [ ]:

import ee
import geemap
import ipywidgets as widgets

import ee import geemap import ipywidgets as widgets

Add some Earth Engine data to the map, then you can use the draw control and a button widget to clip the Earth Engine data by the drawn polygons.

The drawn geometries can be retrieved as either an ee.Geometry object (m.user_roi) or ee.FeatureCollection (m._user_rois).

In [ ]:

self = geemap.Map(center=[0, 9.31], zoom=3)
image = ee.Image("USGS/SRTMGL1_003")
vis_params = {
    "min": 0,
    "max": 6000,
    "palette": "terrain",
}
self.add_layer(image, vis_params, "SRTM")

clip_btn = widgets.Button(description="Clip Image")
reset_btn = widgets.Button(description="Reset")


def on_clip_btn_clicked(b):
    if self.user_roi is not None:
        try:
            clipped_image = image.clip(self.user_roi)
            self.add_layer(clipped_image, vis_params, "Clipped Image")
            self.find_layer("SRTM").visible = False
        except:
            pass


clip_btn.on_click(on_clip_btn_clicked)


def on_reset_btn_clicked(b):
    self._draw_control.clear()
    self.find_layer("SRTM").visible = True
    self.layers = self.layers[:3]


reset_btn.on_click(on_reset_btn_clicked)

widget = widgets.VBox([clip_btn, reset_btn])
self.add_widget(widget, position="bottomright")
self

self = geemap.Map(center=[0, 9.31], zoom=3) image = ee.Image("USGS/SRTMGL1_003") vis_params = { "min": 0, "max": 6000, "palette": "terrain", } self.add_layer(image, vis_params, "SRTM") clip_btn = widgets.Button(description="Clip Image") reset_btn = widgets.Button(description="Reset") def on_clip_btn_clicked(b): if self.user_roi is not None: try: clipped_image = image.clip(self.user_roi) self.add_layer(clipped_image, vis_params, "Clipped Image") self.find_layer("SRTM").visible = False except: pass clip_btn.on_click(on_clip_btn_clicked) def on_reset_btn_clicked(b): self._draw_control.clear() self.find_layer("SRTM").visible = True self.layers = self.layers[:3] reset_btn.on_click(on_reset_btn_clicked) widget = widgets.VBox([clip_btn, reset_btn]) self.add_widget(widget, position="bottomright") self

Create a solara web app.

In [ ]:

import ee
import geemap
import solara
import ipywidgets as widgets


class Map(geemap.Map):
    def __init__(self, **kwargs):
        super().__init__(**kwargs)

        image = ee.Image("USGS/SRTMGL1_003")
        vis_params = {
            "min": 0,
            "max": 6000,
            "palette": "terrain",
        }
        self.add_layer(image, vis_params, "SRTM")

        clip_btn = widgets.Button(description="Clip Image")
        reset_btn = widgets.Button(description="Reset")

        def on_clip_btn_clicked(b):
            if self.user_roi is not None:
                try:
                    clipped_image = image.clip(self.user_roi)
                    self.add_layer(clipped_image, vis_params, "Clipped Image")
                    self.find_layer("SRTM").visible = False
                except:
                    pass

        clip_btn.on_click(on_clip_btn_clicked)

        def on_reset_btn_clicked(b):
            self._draw_control.clear()
            self.find_layer("SRTM").visible = True
            self.layers = self.layers[:3]

        reset_btn.on_click(on_reset_btn_clicked)

        widget = widgets.VBox([clip_btn, reset_btn])
        self.add_widget(widget, position="bottomright")


@solara.component
def Page():
    with solara.Column(style={"min-width": "500px"}):
        Map.element(
            center=[0, 9.31],
            zoom=3,
            height="600px",
        )


Page()

import ee import geemap import solara import ipywidgets as widgets class Map(geemap.Map): def __init__(self, **kwargs): super().__init__(**kwargs) image = ee.Image("USGS/SRTMGL1_003") vis_params = { "min": 0, "max": 6000, "palette": "terrain", } self.add_layer(image, vis_params, "SRTM") clip_btn = widgets.Button(description="Clip Image") reset_btn = widgets.Button(description="Reset") def on_clip_btn_clicked(b): if self.user_roi is not None: try: clipped_image = image.clip(self.user_roi) self.add_layer(clipped_image, vis_params, "Clipped Image") self.find_layer("SRTM").visible = False except: pass clip_btn.on_click(on_clip_btn_clicked) def on_reset_btn_clicked(b): self._draw_control.clear() self.find_layer("SRTM").visible = True self.layers = self.layers[:3] reset_btn.on_click(on_reset_btn_clicked) widget = widgets.VBox([clip_btn, reset_btn]) self.add_widget(widget, position="bottomright") @solara.component def Page(): with solara.Column(style={"min-width": "500px"}): Map.element( center=[0, 9.31], zoom=3, height="600px", ) Page()