USGS-R / drb-gw-hw-model-prep

Code repo to prepare groundwater and headwater-related datasets for modeling river temperature in the Delaware River Basin
Other
0 stars 3 forks source link

25 all geospatial switch to projected crs 5070 #35

Closed msleckman closed 2 years ago

msleckman commented 2 years ago

This PR addresses issue #25 by setting the crs projection to 5070 across the main sf objects such as nhd reaches and catchments, that require raster processing. 4326 proj was kept for p1_drb_temp_sites_sf and p1_dendritic_nhd_reaches_along_NHM (for this one, actually specified that the crs = 4326 when read in) and downstream targets tied to the dataRetrieval nwis data pulls.

Targets that have been reprojected and or modified due to reprojection : p1_nhd_reaches p1_nhd_reaches_along_NHM p1_nhd_catchments
p1_nhm_catchments_dissolved --> no longer get s2_as_sf() warnings p1_ref_gages_sf --> not sure if this one should be changed or not.

p2_buffered_nhd_reaches_along_nhm p2_depth_to_bedrock_reaches_along_nhm p2_depth_to_bedrock_catchments_along_nhm_dissolved

Please note that for the raster processing - actually had to keep the proj to 4326 and then transform the dtb output to 5070 (e.g. 2_process.R#L59 and 2_process.R#L69)

Let me know if this runs for you and if there are any unecessary steps. I approached this as transforming more vs. less.