R/segvert2vector.R

In WeiquanLuo/rivertopo: What the Package Does (One Line, Title Case)

#' convert the segvert.set to vector.set
#' @importFrom magrittr %>%
#' @importFrom sf st_drop_geometry
#' @importFrom dplyr select rename left_join
#' @export
#' @param site_sf a site sf (or dataframe) object that must contain column: id, X, Y, and the snapped point information from the riverdist::xy2segvert() columns: seg, vert
#' @param segvert.set a segvert.set compute by site2segvert() containing column: seg, vert, seg1, vert1, onnected.
#' Variable seg and vert is the snaped point from site to tth river, whereas column seg1, vert1 are the downstream snaped point


# segvert.set to vector.set
segvert2vector <- function(segvert.set, site_sf){
  vector.set <- site_sf %>%
    st_drop_geometry() %>%
    select(id, seg, vert, X, Y, longitude, latitude) %>%
    left_join(segvert.set, by = c("seg", "vert")) %>%
    rename(seg0=seg, vert0 =vert, from = id, x0 = X, y0 = Y,  lon0 = longitude, lat0 = latitude) %>%
    select(-connected); vector.set

  complete_vector.set <- vector.set %>%
    left_join(site_sf %>%
                st_drop_geometry() %>%
                select(id, seg, vert, X, Y, longitude, latitude) %>%
                rename(to=id, seg1=seg, vert1=vert, x1=X, y1=Y, lon1 = longitude, lat1 = latitude),
              by = c("seg1", "vert1")) %>%
    select(from, seg0, vert0, x0, y0, lon0, lat0, to, seg1, vert1, x1, y1, lon1, lat1) %>% 
    distinct; complete_vector.set
  return(complete_vector.set)
}