44 cog stac
Uncomment the following line to install geemap if needed.
# !pip install geemap
Working with Cloud Optimized GeoTIFF (COG)¶
A Cloud Optimized GeoTIFF (COG) is a regular GeoTIFF file, aimed at being hosted on a HTTP file server, with an internal organization that enables more efficient workflows on the cloud. It does this by leveraging the ability of clients issuing HTTP GET range requests to ask for just the parts of a file they need.
More information about COG can be found at https://www.cogeo.org/in-depth.html
Some publicly available Cloud Optimized GeoTIFFs:
- https://stacindex.org/
- https://cloud.google.com/storage/docs/public-datasets/landsat
- https://www.digitalglobe.com/ecosystem/open-data
- https://earthexplorer.usgs.gov/
For this demo, we will use data from https://www.maxar.com/open-data/california-colorado-fires for mapping California and Colorado fires. A List of COGs can be found here.
import os
import ee
import geemap
os.environ["TITILER_ENDPOINT"] = "https://titiler.xyz"
Map = geemap.Map()
Map
url = "https://github.com/opengeos/data/releases/download/raster/Libya-2023-07-01.tif"
geemap.cog_bounds(url)
geemap.cog_center(url)
geemap.cog_bands(url)
geemap.cog_tile(url)
Map.add_cog_layer(url, name="Fire (pre-event)")
url2 = "https://github.com/opengeos/data/releases/download/raster/Libya-2023-09-13.tif"
Map.add_cog_layer(url2, name="Fire (post-event)")
Working with SpatioTemporal Asset Catalog (STAC)¶
The SpatioTemporal Asset Catalog (STAC) specification provides a common language to describe a range of geospatial information, so it can more easily be indexed and discovered. A 'spatiotemporal asset' is any file that represents information about the earth captured in a certain space and time. The initial focus is primarily remotely-sensed imagery (from satellites, but also planes, drones, balloons, etc), but the core is designed to be extensible to SAR, full motion video, point clouds, hyperspectral, LiDAR and derived data like NDVI, Digital Elevation Models, mosaics, etc. More information about STAC can be found at https://stacspec.org/
Some publicly available SpatioTemporal Asset Catalog (STAC):
For this demo, we will use STAC assets from https://stacindex.org/catalogs/spot-orthoimages-canada-2005#/?t=catalogs
Map = geemap.Map()
Map
url = "https://canada-spot-ortho.s3.amazonaws.com/canada_spot_orthoimages/canada_spot5_orthoimages/S5_2007/S5_11055_6057_20070622/S5_11055_6057_20070622.json"
geemap.stac_bounds(url)
geemap.stac_center(url)
geemap.stac_bands(url)
geemap.stac_tile(url, bands=["B3", "B2", "B1"])
Map.add_stac_layer(url, bands=["pan"], name="Panchromatic")
Map.add_stac_layer(url, bands=["B3", "B2", "B1"], name="False color")