dynamic CRS for xy requests

owp-spatial · Feb 20, 2025 · 66a9250 · 66a9250
1 parent 86887f3
commit 66a9250
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 3 deletions.
diff --git a/R/find_origin.R b/R/find_origin.R
@@ -113,10 +113,15 @@ find_origin_query.xy <- function(xy, network, src) {
     src = src
   }
 
-  bb <- sf::st_point(xy) |> 
+  crs <- as_ogr(src, "divides") |> 
+    head(1) |> 
+    sf::st_as_sf() |> 
+    st_crs()
+
+  bb <- sf::st_point(unclass(xy)) |> 
     sf::st_sfc(crs = 4326) |> 
     sf::st_as_sf() |> 
-    sf::st_transform(5070) |> 
+    sf::st_transform(crs$wkt) |> 
     sf::st_geometry() |> 
     sf::st_as_text()
 

diff --git a/R/query_subset.R b/R/query_subset.R
@@ -10,7 +10,7 @@ query_subset <- function(query) {
 
   origin <- find_origin(
     network = query_source_layer(query$source, "network"),
-    src = query$source, 
+    src = query$source[1], 
     id = identifier,
     type = class(identifier)
   )