A collaboration between Open Science Computing and Frontier Geospatial, hence OpenFrontier

• create a VRT from 1000x1000 pixel geotiff files for a single year using a python notebook
• convert the VRT to COG and push to S3

Ran the notebook to create the VRTs locally.

(pangeo) rsignell@OSC:~$ export AWS_ACCESS_KEY_ID=AKIAUxxxxxxxxx (pangeo) rsignell@OSC:~$ export AWS_SECRET_ACCESS_KEY=Wkdzpxxxxxxxxxxxxx (pangeo) rsignell@OSC:~$ export AWS_REGION=us-east-1 (pangeo) rsignell@OSC:~$ export AWS_VIRTUAL_HOSTING=yes (pangeo) rsignell@OSC:~$ export AWS_HTTPS=yes (pangeo) rsignell@OSC:~$ export AWS_DEFAULT_REGION=us-east-1 

Can warp the VRT to Web Mercator projection using:

(pangeo) gdalwarp agb_1990.vrt -s_srs epsg:5070 -t_srs epsg:3857 agb_1990_web.tif 

Convert to COG using rio cogeo:

(pangeo) rsignell@OSC:~$ rio cogeo create agb_1990.vrt agb_1990_cog.tif 

Or convert the warped TIF to a COG using:

(pangeo) rsignell@OSC:~$ rio cogeo create agb_1990_web.tif agb_1990_cog_geo.tif 

Titiler renders COGS as WMTS tiles using AWS lambda. The COGs must be public for titiler, so I moved one to my Open Storage Network s3-compatible endpoint.

Here's one of the COGs viewed through Development Seeds's titiler endpoint:

https://titiler.xyz/cog/viewer?url=https://usgs.osn.mghpcc.org/esip/rsignell/testing/agb_1990_cog.tif

(Need to scale manually from 0 to 300 until I figure out how to specify that on the URL)