updating nexus on parquet file #817
Conversation
AminTorabi-NOAA
requested review from
kumdonoaa,
shorvath-noaa and
JurgenZach-NOAA
There was a problem hiding this comment.
I think we need to be careful with averaging the water depth. Look at flowpath ID 2427759. This flowpath lies just downstream of a waterbody. Waterbodies tend to always have much greater depths than flowpaths, so when you average depths of this ID and the waterbody above it, the depth is artificially inflated. Depth before averaging, ((287, d) 1.462409):
But after doing the conversion from flowpaths to nexus points (22.146202m):
|
|
||
| def crosswalk_nex_flowpath_poi(self): | ||
| self._nexus_dict = dict() | ||
| for key, values in self._connections.items(): | ||
| for value in values: | ||
| new_key = f"nex-{value}" | ||
| new_value = f"wb-{key}" | ||
| if new_key not in self.nexus_dict: | ||
| self._nexus_dict[new_key] = [] | ||
| self._nexus_dict[new_key].append(new_value) |
There was a problem hiding this comment.
NHD network doesn't have nexus points, so I'm not sure this is needed here. There should be some method for NHD just to prevent errors, but NHD output will always be indexed by flowpath IDs because that's all it has.
There was a problem hiding this comment.
That's True, and already the output of NHD is by flowpath IDs, I just created this artificially dictionary to use the same function (updated_flowveldepth) I created for analyzing cross-walk between nexus point and flowpath ids in hyfeature. (To avoid repetition of a same function)
AminTorabi-NOAA
force-pushed
the
parquet_update
branch
from
daf4867 to
e7bc8af
Compare
|
@jameshalgren Would you look at this PR, I think this PR and your PR (#842) are related |
| if nexus_dict: | ||
| location_ids = prefix_ids + '-' + df['location_id'].astype(str) | ||
| else: | ||
| location_ids = df['location_id'].astype(str) |
There was a problem hiding this comment.
I think this is a breaking change. For example if your configuration looked like:
output_parameters:
parquet_output:
parquet_output_folder: output/
prefix_ids: wbThe location_id column would no longer be prefixed with wb- (e.g. wb-111 would now be 111).
There was a problem hiding this comment.
Can you elaborate more? I think it still have the prefix if we change the prefix_ids to wb
| parquet_output_folder: output/ | ||
| configuration: short_range | ||
| prefix_ids: nex | ||
| # parquet_output: |
There was a problem hiding this comment.
Can you elaborate on why this was commented out?
There was a problem hiding this comment.
Just to simplify the yaml file. So if user need that, can uncomment it
This PR updates the parquet file output if its nexus point is requested, It creates a crosswalk between upstream and downstream (as nex-, wb-) and if there are more than one segment going to a nexus points it add up the flowrate and average the waterdepth.
Additions
Removals
Changes
Testing
Screenshots
Notes
Todos
Checklist
Testing checklist
Target Environment support
Accessibility
Other