[5pt] PR: Update nws_lid dataset with Alaska sites #1159
Merged
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
This PR addresses issue #1132 and produces an updated input dataset to replace
/dev_fim_share/foss_fim/inputs/ahps_sites/nws_lid.gpkg
with/dev_fim_share/foss_fim/outputs/ali_nws_lid/nws_lid.gpkg
The PR includes the following changes on
tools/generate_nws_lid.py
for updatingnws_lid.gpkg
dataset.In this revised version, stations only from these two groups are retrieved:
rfc_forecast_point= True
/data/inputs/ahp_sites/evaluated_ahps_sites.csv
The lid stations in AK (Alaska), HI, and PR, with above two criteria have also been selected, as shown in the map below. In the previous version of the code, all of lid stations in PR and HI (regardless of meeting above two criteria), were also being retrieved. I have updated this version to exclude such stations.
Also, In this revised version, I've eliminated the code sections that previously generated the "is_headwater" and "is_colocated" columns, which are not needed in FIM4. Therefore, in this updated version, these columns are no longer present.
Similar to 'usgs_gages.gpkg' dataset, all lid stations, including those in Alaska, are stored in a single gpkg file (
nws_lid.gpkg
) with EPSG=5070. The Alaska stations can be identified using their HUC8 numbers (beginning with '19').Changes
Testing
The
nws_lid.gpkg
dataset is used only inusgs_gage_unit_setup.py
for generatingusgs_subset_gages.gpkg
for simulated HUC8s. I confirmed the validity of this PR by running FIM for a HUC8 using the updatednws_lid.gpkg-- contents of
usgs_subset_gages.gpkgwere aligned with the data in the updated
nws_lid.gpkg`.This map displays 4381 NWS_LID points compiled using the updated tools/generate_nws_lid.py:
Deployment Plan (For developer use)
How does the changes affect the product?
Issuer Checklist (For developer use)
You may update this checklist before and/or after creating the PR. If you're unsure about any of them, please ask, we're here to help! These items are what we are going to look for before merging your code.
[_pt] PR: <description>
dev
branch (the default branch), you have a descriptive Feature Branch name using the format:dev-<description-of-change>
(e.g.dev-revise-levee-masking
)dev
branchpre-commit
hooks were run locally/foss_fim/
, run:pytest unit_tests/
)4.x.x.x
Merge Checklist (For Technical Lead use only)