42 create cog

Uncomment the following line to install leafmap if needed.

In [1]:

# !pip install leafmap

# !pip install leafmap

In [2]:

import leafmap

import leafmap

Provide a dataset path or URL.

In [3]:

url = "https://github.com/opengeos/data/raw/main/raster/srtm90.tif"

url = "https://github.com/opengeos/data/raw/main/raster/srtm90.tif"

Validate COG.

In [4]:

leafmap.cog_validate(url)

leafmap.cog_validate(url)

The following warnings were found:
- The file is greater than 512xH or 512xW, it is recommended to include internal overviews

Out[4]:

(True,
 [],
 ['The file is greater than 512xH or 512xW, it is recommended to include internal overviews'])

In [5]:

leafmap.cog_validate(url, verbose=True)

leafmap.cog_validate(url, verbose=True)

Out[5]:

Info(Path='https://github.com/opengeos/data/raw/main/raster/srtm90.tif', Driver='GTiff', COG=True, Compression='LZW', ColorSpace=None, COG_errors=None, COG_warnings=['The file is greater than 512xH or 512xW, it is recommended to include internal overviews'], Profile=Profile(Bands=1, Width=4269, Height=2465, Tiled=True, Dtype='int16', Interleave='BAND', AlphaBand=False, InternalMask=False, Nodata=None, ColorInterp=('gray',), ColorMap=False, Scales=(1.0,), Offsets=(0.0,)), GEO=Geo(CRS='EPSG:4326', BoundingBox=(-120.75592734926073, 36.63401097629554, -117.30451019614512, 38.6269234341147), Origin=(-120.75592734926073, 38.6269234341147), Resolution=(0.0008084837557075694, -0.0008084837557075694), MinZoom=7, MaxZoom=11), Tags={'Image Metadata': {'AREA_OR_POINT': 'Area'}, 'Image Structure': {'COMPRESSION': 'LZW', 'INTERLEAVE': 'BAND'}}, Band_Metadata={'Band 1': BandMetadata(Description='elevation', ColorInterp='gray', Offset=0.0, Scale=1.0, Metadata={})}, IFD=[IFD(Level=0, Width=4269, Height=2465, Blocksize=(256, 256), Decimation=0)])

Convert the image to tiled COG.

In [6]:

out_cog = "cog.tif"
leafmap.image_to_cog(url, out_cog)

out_cog = "cog.tif" leafmap.image_to_cog(url, out_cog)

Reading input: https://github.com/opengeos/data/raw/main/raster/srtm90.tif

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

Validate COG.

In [7]:

leafmap.cog_validate(out_cog)

leafmap.cog_validate(out_cog)

Out[7]:

(True, [], [])

In [8]:

leafmap.cog_validate(out_cog, verbose=True)

leafmap.cog_validate(out_cog, verbose=True)

Out[8]:

Info(Path='/home/runner/work/leafmap/leafmap/docs/notebooks/cog.tif', Driver='GTiff', COG=True, Compression='DEFLATE', ColorSpace=None, COG_errors=None, COG_warnings=None, Profile=Profile(Bands=1, Width=4269, Height=2465, Tiled=True, Dtype='int16', Interleave='BAND', AlphaBand=False, InternalMask=False, Nodata=None, ColorInterp=('gray',), ColorMap=False, Scales=(1.0,), Offsets=(0.0,)), GEO=Geo(CRS='EPSG:4326', BoundingBox=(-120.75592734926073, 36.63401097629554, -117.30451019614512, 38.6269234341147), Origin=(-120.75592734926073, 38.6269234341147), Resolution=(0.0008084837557075694, -0.0008084837557075694), MinZoom=7, MaxZoom=11), Tags={'Image Metadata': {'AREA_OR_POINT': 'Area', 'OVR_RESAMPLING_ALG': 'NEAREST'}, 'Image Structure': {'COMPRESSION': 'DEFLATE', 'INTERLEAVE': 'BAND', 'LAYOUT': 'COG'}}, Band_Metadata={'Band 1': BandMetadata(Description='elevation', ColorInterp='gray', Offset=0.0, Scale=1.0, Metadata={})}, IFD=[IFD(Level=0, Width=4269, Height=2465, Blocksize=(512, 512), Decimation=0), IFD(Level=1, Width=2135, Height=1233, Blocksize=(512, 512), Decimation=2), IFD(Level=2, Width=1068, Height=617, Blocksize=(512, 512), Decimation=4), IFD(Level=3, Width=534, Height=309, Blocksize=(512, 512), Decimation=8)])

Add COG to map.

In [9]:

m = leafmap.Map()
m.add_raster(out_cog, palette="terrain", layer_name="Local COG")
m.add_cog_layer(url, palette="gist_earth", name="Remote COG")
m

m = leafmap.Map() m.add_raster(out_cog, palette="terrain", layer_name="Local COG") m.add_cog_layer(url, palette="gist_earth", name="Remote COG") m

Out[9]:

---

Last update: 2023-04-05
Created: 2021-06-10