if (FALSE) {
dd <- data.frame(x = c(-122.649,-100.348,-75.186,-106.675),
y = c(45.085, 35.405,42.403,38.721))
comids <- sc_get_comid(dd, xcoord='x',
ycoord='y', crsys=4269)
dd <- sf::st_point_on_surface(sf::read_sf(system.file("shape/nc.shp", package="sf")))
comids <- sc_get_comid(dd)
comids <- sc_get_comid(dd, xcoord='x',
ycoord='y', crsys=4269)
dd <- sf::read_sf(system.file("shape/nc.shp", package="sf"))
comids <- sc_get_comid(dd)
}
Run the code above in your browser using DataLab