search

Rijkswaterstaat / wm-ws-dl

wm-ws-dl documentation
https://rijkswaterstaatdata.nl/waterdata
11 stars 2 forks source link

hoedanigheid (vertical reference level) is defined in multiple metadata fields #42

Closed veenstrajelmer closed 4 months ago

veenstrajelmer commented 4 months ago

When filtering waterlevel extremes, we have to distinguish between NAP/MSL stations in the Groepering column (or include all). The vertical reference level (NAP or MSL) is already defined in the Hoedanigheid column, so it does not make sense to also distinguish between these in the Groepering : GETETM2/GETETBRKD2/GETETMSL2/GETETBRKDMSL2

import ddlpy

locations = ddlpy.locations()

bool_station = locations.index.isin(['HOEKVHLD'])
bool_grootheid_meas = locations['Grootheid.Code'].isin(['WATHTE'])
bool_grootheid_astro = locations['Grootheid.Code'].isin(['WATHTBRKD'])
bool_hoedanigheid = locations['Hoedanigheid.Code'].isin(['NAP'])
bool_groepering_ts = locations['Groepering.Code'].isin(['NVT'])
bool_groepering_ext = locations['Groepering.Code'].isin(['GETETM2','GETETBRKD2','GETETMSL2','GETETBRKDMSL2'])

# get WATHTE locations (ts) >> measured waterlevel
locs_wathte = locations.loc[bool_station & bool_grootheid_meas & bool_hoedanigheid & bool_groepering_ext]
# get WATHTBRKD locations (ts) >> computed astronomical waterlevel
locs_wathtbrkd = locations.loc[bool_station & bool_grootheid_astro & bool_hoedanigheid & bool_groepering_ext]
TvLoon-RWS commented 4 months ago

This will be partly covered with the ProcesType 'meting' 'astro' or 'verwachting' in Wadar WaterWebservices