77 split raster

Uncomment the following line to install leafmap if needed.

In [1]:

# %pip install -U leafmap

# %pip install -U leafmap

In [2]:

import leafmap

import leafmap

Create an interactive map and add a raster layer to it.

In [3]:

m = leafmap.Map()
url = "https://open.gishub.org/data/raster/cog.tif"
m.add_cog_layer(url, name="COG")
m

m = leafmap.Map() url = "https://open.gishub.org/data/raster/cog.tif" m.add_cog_layer(url, name="COG") m

Out[3]:

Split the raster layer into multiple tiles.

In [4]:

leafmap.split_raster(url, out_dir="tiles", tile_size=(100, 100), overlap=0)

leafmap.split_raster(url, out_dir="tiles", tile_size=(100, 100), overlap=0)

cog.tif already exists. Skip downloading. Set overwrite=True to overwrite.

Show the tiles.

In [5]:

leafmap.find_files("tiles", ext=".tif")

leafmap.find_files("tiles", ext=".tif")

Out[5]:

['tiles/tile_0_0.tif',
 'tiles/tile_0_1.tif',
 'tiles/tile_0_2.tif',
 'tiles/tile_1_0.tif',
 'tiles/tile_1_1.tif',
 'tiles/tile_1_2.tif',
 'tiles/tile_2_0.tif',
 'tiles/tile_2_1.tif',
 'tiles/tile_2_2.tif',
 'tiles/tile_3_0.tif',
 'tiles/tile_3_1.tif',
 'tiles/tile_3_2.tif']

Merge the tiles into a single raster dataset.

In [6]:

leafmap.merge_rasters("tiles", output="merged.tif", input_pattern="*.tif")

leafmap.merge_rasters("tiles", output="merged.tif", input_pattern="*.tif")