#!/usr/bin/env python
# coding: utf-8

# # Explore Estuarine Bathymetry on https://source.coop
# 30m bathymetry DEM data from NOAA NOS, obtained via bittorrent on May 2, 2025 and synced to source.coop: 
# ```
# $ conda install aria2
# $ aria2c  https://sciop.net/torrents/9840fa77d032a3821180f7f1726df02431e37e53/estuarine-bathymetric-digital-elevation-models.torrent
# $ aws s3 sync s3 ./estuarine-bathymetry-digital-elevation-models/ s3://rsignell/ncei-estuarine-bathymetry/ --profile=sc
# ```
# 

# In[1]:


import fsspec
import xarray as xr


# In[2]:


fs = fsspec.filesystem('s3', anon=True, endpoint_url='https://data.source.coop')


# In[3]:


fs.ls('s3://rsignell/ncei-estuarine-bathymetry/')


# In[4]:


fs.ls('s3://rsignell/ncei-estuarine-bathymetry/Buzzards_Bay')


# In[5]:


ds = xr.open_dataset(fs.open('rsignell/ncei-estuarine-bathymetry/Buzzards_Bay/buzzards_bay_m010_30m.nc'))


# In[6]:


ds


# In[7]:


import hvplot.xarray


# In[8]:


ds['Band1'].hvplot(x='x', y='y', crs='EPSG:26919', rasterize=True, tiles='OSM', cmap='viridis')


# In[ ]: