Create a fresh conda env to run this example if needed.
conda create -n cog python=3.9
conda install mamba -c conda-forge
mamba install geemap rio-cogeo -c conda-forge
In [ ]:
# !pip install geemap rio-cogeo
In [ ]:
import geemap
In [ ]:
url = 'https://github.com/giswqs/leafmap/raw/master/examples/data/cog.tif'
in_cog = 'cog.tif'
out_cog = "ndvi.tif"
Download a sample dataset.
In [ ]:
geemap.download_from_url(url, in_cog)
Convert image to numpy array.
In [ ]:
arr = geemap.image_to_numpy(in_cog)
In [ ]:
arr.shape
Computer NDVI.
In [ ]:
ndvi = (arr[3] - arr[0]) / (arr[3] + arr[0])
In [ ]:
ndvi.shape
Convert numpy array to COG.
In [ ]:
geemap.numpy_to_cog(ndvi, out_cog, profile=in_cog)
In [ ]:
m = geemap.Map()
m.add_raster(in_cog, band=[4, 1, 2], layer_name="Color infrared")
m.add_raster(out_cog, palette="Greens", layer_name="NDVI")
m
