import os
os.environ['GDAL_DISABLE_READDIR_ON_OPEN']='EMPTY_DIR'
import intake
from satsearch import Search
def get_STAC_items(url, collection, dates, bbox):
results = Search.search(url=url,
collections=[collection],
datetime=dates,
bbox=bbox,
sortby=['-properties.datetime'])
items = results.items()
print(f'Found {len(items)} Items')
return intake.open_stac_item_collection(items)
url='https://cmr.earthdata.nasa.gov/stac/ASF'
bbox = [-122.4, 41.3, -122.1, 41.5] # (min lon, min lat, max lon, max lat)
collection = 'C1595422627-ASF' # sentinel-1 beta interferograms
dates = '2020-01-01/2020-12-31'
items_s1 = get_STAC_items(url, collection, dates, bbox)
url='https://earth-search.aws.element84.com/v0'
collection = 'sentinel-s2-l2a-cogs' # sentinel-2 w/ atmospheric corrections
items_s2 = get_STAC_items(url, collection, dates, bbox)
gf = items_s2.to_geopandas()
gf.head()
%%time
item = items_s2['S2A_10TEL_20201130_0_L2A']
da = item.B04.to_dask()
da
%%time
import hvplot.xarray
subset = da.isel(x=slice(6000,10980), y=slice(0,3000))
title = item.name + ' - B04'
subset.hvplot.image(rasterize=True, aspect='equal', cmap='magma', frame_width=500, title=title)