Skip to content

Instantly share code, notes, and snippets.

View vincentsarago's full-sized avatar

Vincent Sarago vincentsarago

View GitHub Profile
@vincentsarago
vincentsarago / app.py
Created November 17, 2022 16:11
OGC Features and Tiles API
"""OGC Features and Tiles API."""
import re
from fastapi import Depends, FastAPI, HTTPException, Path, Query
from starlette.requests import Request
from starlette.responses import Response, HTMLResponse
from starlette_cramjam.middleware import CompressionMiddleware
titiler.application==0.7.0
server-thread>=0.2.0
folium
httpx
@vincentsarago
vincentsarago / README.md
Created April 19, 2022 15:45
SUM overview + cogger

1. Create input tif

gdal_translate file.tif input.tif -of GTiff -co TILED=TRUE -co COMPRESS=DEFLATE

2. Create Overview

python create_overview.py input.tif

3. Create COG

Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
<!DOCTYPE html>
<html>
<head>
<meta charset='utf-8' />
<title>Simple STAC API Viewer</title>
<meta name='viewport' content='initial-scale=1,maximum-scale=1,user-scalable=no' />
<script src='https://api.tiles.mapbox.com/mapbox-gl-js/v1.6.1/mapbox-gl.js'></script>
<link href='https://api.tiles.mapbox.com/mapbox-gl-js/v1.6.1/mapbox-gl.css' rel='stylesheet' />
This file has been truncated, but you can view the full file.
{"mosaicjson": "0.0.2", "version": "1.0.0", "minzoom": 8, "maxzoom": 12, "quadkey_zoom": 8, "bounds": [-151.74937671440708, -65.08963152892852, 150.96975489094262, 84.53821238222514], "center": [-0.3898109117322264, 9.72429042664831, 8], "tiles": {"00111301": ["http://radarsat-r1-l1-cog.s3.amazonaws.com/1996/10/RS1_X0503117_SCWB_19961023_090415_HH_SCW.tif", "http://radarsat-r1-l1-cog.s3.amazonaws.com/1996/10/RS1_X0503119_SCWB_19961030_075958_HH_SCW.tif"], "00111310": ["http://radarsat-r1-l1-cog.s3.amazonaws.com/1996/10/RS1_X0503117_SCWB_19961023_090415_HH_SCW.tif", "http://radarsat-r1-l1-cog.s3.amazonaws.com/1996/10/RS1_X0503119_SCWB_19961030_075958_HH_SCW.tif", "http://radarsat-r1-l1-cog.s3.amazonaws.com/1996/11/RS1_X0503120_SCWB_19961106_075548_HH_SCW.tif"], "00111311": ["http://radarsat-r1-l1-cog.s3.amazonaws.com/1996/10/RS1_X0503117_SCWB_19961023_090415_HH_SCW.tif", "http://radarsat-r1-l1-cog.s3.amazonaws.com/1996/10/RS1_X0503119_SCWB_19961030_075958_HH_SCW.tif", "http://radarsat-r1-l1-cog.s3.amazonaws.co
{"geometry": {"type": "Polygon", "coordinates": [[[-118.86884165050823, 56.44911916433228], [-118.86884165050823, 52.826289010731195], [-112.57569524761153, 52.826289010731195], [-112.57569524761153, 56.44911916433228], [-118.86884165050823, 56.44911916433228]]]}, "properties": {"path": "http://radarsat-r1-l1-cog.s3.amazonaws.com/1996/10/RS1_C0003310_SCNA_19961004_140742_HH_SCN.tif", "bounds": [-118.86884165050823, 52.826289010731195, -112.57569524761153, 56.44911916433228], "minzoom": 6, "maxzoom": 12, "datatype": "uint8"}, "type": "Feature"}
{"geometry": {"type": "Polygon", "coordinates": [[[-105.06720960654427, 50.710779349317086], [-105.06720960654427, 50.207840157795054], [-104.28207563625668, 50.207840157795054], [-104.28207563625668, 50.710779349317086], [-105.06720960654427, 50.710779349317086]]]}, "properties": {"path": "http://radarsat-r1-l1-cog.s3.amazonaws.com/1996/10/RS1_C0003495_F5_19961013_005804_HH_SGX.tif", "bounds": [-105.06720960654427, 50.207840157795054, -104.28207563625668, 50.7107793493
@attr.s
class SDSTReader(COGReader):
"""Cloud Optimized GeoTIFF Reader."""
subdatasets: Dict[int, str] = attr.ib(init=False)
def __attrs_post_init__(self):
"""Define _kwargs, open dataset and get info."""
super().__attrs_post_init__()
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.