using ERDDAP RESTful queries, pandas and hvplot
import pandas as pd
import hvplot.pandas
endpoint = 'https://geoport.usgs.esipfed.org/erddap'
project = 'W_GMAINE'
datasetID, lon, lat
for all datasets¶df = pd.read_csv(f'{endpoint}/tabledap/allDatasets.csv?datasetID%2CminLongitude%2CminLatitude')
ids = pd.read_csv(f'{endpoint}/search/index.csv?page=1&itemsPerPage=10000&searchFor={project}')['Dataset ID'].values
df = df[df['datasetID'].isin(ids)]
df = df.rename(columns={'minLongitude':'lon', 'minLatitude':'lat', 'datasetID':'mooring'})
df = df.set_index('mooring')
df['lon']= df['lon'].astype('float')
df['lat']= df['lat'].astype('float')
df.head()
lon | lat | |
---|---|---|
mooring | ||
10071sc-a | -69.3578 | 43.7148 |
10081sc-a | -68.1087 | 44.1058 |
10091sc-a | -70.4278 | 43.1807 |
10341sc-a | -69.3548 | 43.7158 |
10431sc-a | -70.5657 | 42.5217 |
df.hvplot.points(x='lon', y='lat', geo=True, tiles='OSM', marker='square', color='green')