search

DOI-USGS / mda.streams

backend tools for powstreams
Creative Commons Zero v1.0 Universal
3 stars 8 forks source link

use different CRS in stage_ldas_ts? #279

Open aappling-usgs opened 8 years ago

aappling-usgs commented 8 years ago

@jread-usgs, in light of https://github.com/USGS-CIDA/stream_metab_usa/pull/121#discussion_r70883034

there's a column in get_meta('basic') called coord_datum. The values there are almost always NAD83 (723 of those, 1 of WGS84). does that mean we should be using EPSG:4269 for most sites and EPSG: 4326 just for that one WGS84 site? (i'm getting my info from https://www.nceas.ucsb.edu/~frazier/RSpatialGuides/OverviewCoordinateReferenceSystems.pdf)

should/can we specify the CRS when creating the site stencil in stage_ldas_ts at https://github.com/USGS-R/mda.streams/blob/3e9791c0d6ce407d53b989f0cc60614e20ce89ad/R/stage_ldas_ts.R#L55?

stencil <- simplegeom(lon_lat_df)
jordansread commented 7 years ago

yes, what I think we should do here is do a spatial transform to get all of them into the datum and proj that geoknife uses for the simplegeom class. Good catch.