Skip to content

51 clip image

image image image

Uncomment the following line to install leafmap if needed.

1
# !pip install leafmap

1
# !pip install rasterio fiona

1
import leafmap

Download a sample raster dataset.

1
2
url = "https://opengeos.org/data/raster/srtm90.tif"
dem = "dem.tif"

1
leafmap.download_file(url, dem, overwrite=True)

Create an interactive map.

1
2
3
m = leafmap.Map()
m.add_raster(dem, palette="terrain", layer_name="DEM")
m

Define a mask to extract the image. The mask can be a string representing a file path to a vector dataset (e.g., geojson, shp), or a list of coordinates (e.g., [[lon,lat], [lon,lat]]), or a dictionary representing a feature (e.g., m.user_roi).

For example, the mask can be a filepath to a vector dataset.

1
# mask = 'https://raw.githubusercontent.com/opengeos/leafmap/master/examples/data/mask.geojson'

Or you can draw a polygon on the map, then use m.user_roi as the mask.

1
# mask = m.user_roi

Or specify a list of coordinates [lon, lat] as the mask.

1
2
3
4
5
6
7
mask = [
    [-119.679565, 37.256566],
    [-119.679565, 38.061067],
    [-118.24585, 38.061067],
    [-118.24585, 37.256566],
    [-119.679565, 37.256566],
]

Specify the output filename.

1
output = "clip.tif"

Clip image by mask.

1
2
3
4
try:
    leafmap.clip_image(dem, mask, output)
except Exception as e:
    print(e)

Add the clipped image to the map.

1
2
3
4
try:
    m.add_raster(output, palette="gist_earth", layer_name="Clip Image")
except Exception as e:
    print(e)