51 clip image

Uncomment the following line to install leafmap if needed.

In [1]:

# !pip install leafmap

# !pip install leafmap

In [2]:

# !pip install rasterio fiona

# !pip install rasterio fiona

In [3]:

import leafmap

import leafmap

Download a sample raster dataset.

In [4]:

url = 'https://open.gishub.org/data/raster/srtm90.tif'
dem = 'dem.tif'

url = 'https://open.gishub.org/data/raster/srtm90.tif' dem = 'dem.tif'

In [5]:

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

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

Downloading...
From: https://open.gishub.org/data/raster/srtm90.tif
To: /home/runner/work/leafmap/leafmap/docs/notebooks/dem.tif
100%|██████████| 15.2M/15.2M [00:00<00:00, 253MB/s]

Out[5]:

'/home/runner/work/leafmap/leafmap/docs/notebooks/dem.tif'

Create an interactive map.

In [6]:

m = leafmap.Map()
m.add_raster(dem, palette='terrain', layer_name="DEM")
m

m = leafmap.Map() m.add_raster(dem, palette='terrain', layer_name="DEM") m

Out[6]:

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.

In [7]:

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

# 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.

In [8]:

# mask = m.user_roi

# mask = m.user_roi

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

In [9]:

mask = [
    [-119.679565, 37.256566],
    [-119.679565, 38.061067],
    [-118.24585, 38.061067],
    [-118.24585, 37.256566],
    [-119.679565, 37.256566],
]

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.

In [10]:

output = 'clip.tif'

output = 'clip.tif'

Clip image by mask.

In [11]:

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

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

Reading input: /home/runner/work/leafmap/leafmap/docs/notebooks/clip.tif

Adding overviews...
Updating dataset tags...
Writing output to: /home/runner/work/leafmap/leafmap/docs/notebooks/clip.tif

Add the clipped image to the map.

In [12]:

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

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

---

Last update: 2023-09-22
Created: 2023-09-22