# One month of global rainfall; clip to -75,-20,-50,5 with GDAL
curl -fLO "https://data.chc.ucsb.edu/products/CHIRPS-2.0/global_monthly/tifs/chirps-v2.0.2020.01.tif.gz"
gunzip chirps-v2.0.2020.01.tif.gz
gdal_translate -projwin -75 5 -50 -20 chirps-v2.0.2020.01.tif clip.tif
import gzip, shutil, urllib.request
import rasterio
from rasterio.windows import from_bounds
url = "https://data.chc.ucsb.edu/products/CHIRPS-2.0/global_monthly/tifs/chirps-v2.0.2020.01.tif.gz"
urllib.request.urlretrieve(url, "chirps.tif.gz")
with gzip.open("chirps.tif.gz") as src, open("chirps.tif", "wb") as dst:
shutil.copyfileobj(src, dst)
with rasterio.open("chirps.tif") as ds:
rain = ds.read(1, window=from_bounds(-75, -20, -50, 5, ds.transform))
const res = await fetch("https://hydrological.org/api/v1/datasets/chirps/links?bbox=-75,-20,-50,5");
const { groups } = await res.json();
for (const g of groups) g.links.forEach(l => console.log(l.url));