st_as_sfc.pq_geometry

[SymbolixAU]

If a table is read from a postgis database as a regular data.frame using RPostgres::dbGetQuery the geometry column is returned as class pq_geometry.

This can be converted to sfc by

sf:::st_as_sfc.WKB( df$geom, EWKB = TRUE )

Can we have a sf_as_sfc.pq_geometry method which does this automatically as part of sf?

For example

#' @export
sf_as_sfc.pq_geometry <- function(x, ..., EWKB = FALSE, spatialite = FALSE, pureR = FALSE, crs = NA_crs_) {
  st_as_sfc.WKB(x, ..., EWKB = TRUE, spatiallite = spatialite, pureR = pureR, crs = crs)
}

[edzer]

Makes sense (@etiennebr ?) - doesn't it make more sense to use

#' @export
sf_as_sfc.pq_geometry <- function(x, ..., EWKB = TRUE, spatialite = FALSE, pureR = FALSE, crs = NA_crs_) {
  st_as_sfc.WKB(x, ..., EWKB = EWKB, spatiallite = spatialite, pureR = pureR, crs = crs)
}

or is the column guaranteed to be EWKB?

[dcooley]

@edzer sorry yes, that makes more sense. I missed setting the EWKB argument to TRUE.