EPA Water Quality Exchange Plugin

[dschristianson]

Add plugin for EPA Water Quality eXchange Data acquired thru the Water Quality Portal: https://www.waterqualitydata.us/

Type of change

• [x] New feature (non-breaking change which adds functionality)
• [x] This change requires a documentation update

How Has This Been Tested?

• [x] Unit tests ~- [ ] Integration tests~
• [x] Test coverage >= 90%
• [x] Flake8 Tests
• [x] Mypy Tests
• [x] Other - manual testing hitting actual data source see forthcoming comment for queries

Test Configuration

• Python Version: 3.9

PR Self Evaluation

• [x] My code follows the agreed upon best practices
• [x] I have commented my code, particularly in hard-to-understand areas
• [x] I have made corresponding changes to the documentation
• [x] My changes generate no new warnings
• [x] I have added tests or modified existing tests that prove my fix is effective or that my feature works
• [x] Existing unit tests pass locally with my changes ~- [ ] Any dependent changes have been merged and published in the appropriate modules~
• [x] I have performed a self-review of my own code

[dschristianson]

A sample of manual tests. Comment / uncomment the various versions.

Setup

from basin3d.plugins import epa
from basin3d import monitor
from basin3d.core.models import MonitoringFeature

monitor.configure(loggers={"basin3d": {"level": "DEBUG"}})
synthesizer = synthesis.register()

Monitoring Feature queries:

# USGS huc
monitoring_features = synthesizer.monitoring_features(parent_feature='EPA-14020001')
# EPA site ids
# monitoring_features = synthesizer.monitoring_features(monitoring_feature=['EPA-CORIVWCH_WQX-176', 'EPA-11NPSWRD_WQX-BLCA_09128000'])

mf: MonitoringFeature
for mf in monitoring_features:
    lat = mf.coordinates.absolute.horizontal_position[0].latitude
    long = mf.coordinates.absolute.horizontal_position[0].longitude
    alt = 'Nope'
    if mf.coordinates.absolute.vertical_extent:
        alt = mf.coordinates.absolute.vertical_extent[0].value
    depth = 'Nope'
    if mf.coordinates.representative:
        depth = mf.coordinates.representative.vertical_position.value
    mf_details = f'{mf.id}; {mf.name}; LAT {lat}; LON {long}; DEPTH {depth}; ELEV {alt}'
    mf_list.update({mf.id: {'details': mf_details, 'desc': mf.description, 'mf': mf}})
    print(f'{mf.id}-{mf.name}')

MeasurementTVPObservation queries

# 4 objects returned
# measurement_timeseries_tvp_observations = synthesizer.measurement_timeseries_tvp_observations(
#     monitoring_feature=["EPA-11NPSWRD_WQX-BLCA_09128000", "EPA-11NPSWRD_WQX-CURE_09127000", "11NPSWRD_WQX-CURE_38193410713350"], observed_property=['SWL'], start_date='2005-01-01', end_date='2010-01-01', aggregation_duration='NONE')
# 5 objects returned
measurement_timeseries_tvp_observations = synthesizer.measurement_timeseries_tvp_observations(
    monitoring_feature=['EPA-CCWC-COAL-26', 'EPA-CCWC-MM-29 WASH #3'], observed_property=['As', 'WT', 'DO'], start_date='2010-01-01', end_date='2011-01-01', aggregation_duration='NONE')

for mvp in measurement_timeseries_tvp_observations:
    foi: MonitoringFeature = mvp.feature_of_interest
    lat = foi.coordinates.absolute.horizontal_position[0].latitude
    lon = foi.coordinates.absolute.horizontal_position[0].longitude
    depth = 'Nope'
    if foi.coordinates.representative:
        depth = f'{foi.coordinates.representative.vertical_position.type}: {foi.coordinates.representative.vertical_position.value}'
    elev = 'Nope'
    if foi.coordinates.absolute.vertical_extent:
        elev = f'{foi.coordinates.absolute.vertical_extent[0].type}: {foi.coordinates.absolute.vertical_extent[0].value}'
    mvp_ids = f'{mvp.id}; {mvp.observed_property}; {foi.id}; LAT {lat}; LONG {lon}; depth {depth}; elev {elev}'
    print(mvp_ids)

[dschristianson]

@vchendrix: please have a look at the new readthedoc plugins page in the documentation. The build is under this branch.

[dschristianson]

Thanks for the great comments, @vchendrix! Address them all plus the few we found during today's discussion.