replace duckspatial with terra + sf

duckspatial was causing CI/CD problems.

Author: ashiklom

docker/depends/pecan_package_dependencies.csv

@@ -66,7 +66,6 @@
 "dplyr","*","modules/uncertainty","Imports",FALSE
 "dplyr",">= 0.8.1","modules/data.atmosphere","Imports",FALSE
 "dplyr",">= 1.1.2","base/db","Imports",FALSE
-"duckspatial","*","modules/data.land","Imports",FALSE
 "ecmwfr",">= 2.0.0","modules/data.atmosphere","Suggests",FALSE
 "ellipse","*","modules/assim.batch","Imports",FALSE
 "exactextractr","*","modules/assim.sequential","Suggests",FALSE

modules/data.land/DESCRIPTION

@@ -29,8 +29,6 @@ Imports:
     curl,
     doSNOW,
     dplyr,
-    duckdb,
-    duckspatial,
     dplR,
     foreach,
     fs,

modules/data.land/R/get-landiq-parcels.R

@@ -7,35 +7,15 @@
 #' @return `design_points` `data.frame` with harmonized LandIQ parcel_IDs
 #' @export
 get_landiq_parcel_ids <- function(design_points, parcels_file) {
-  parcel_crs <- sf::st_layers(parcels_file)[["crs"]][[1]]
+  parcels_vect <- terra::vect(parcels_file)
   pts_sf <- design_points |>
     dplyr::select("id", "lat", "lon") |>
     sf::st_as_sf(coords = c("lon", "lat"), crs = 4326) |>
-    sf::st_transform(parcel_crs)
-  conn <- duckspatial::ddbs_create_conn()
-  duckspatial::ddbs_write_vector(
-    conn = conn,
-    data = pts_sf,
-    name = "design_points"
-  )
-  duckdb::dbSendQuery(
-    conn,
-    glue::glue(
-      "
-      CREATE TABLE merged AS
-      SELECT dp.*, p.parcel_id,
-      FROM design_points dp
-      LEFT JOIN ST_Read('{parcels_file}', layer='parcels') p
-      ON ST_Within(dp.geometry, p.geom)
-      "
-    )
-  )
-  dp_with_parcels <- duckspatial::ddbs_read_vector(
-    conn = conn,
-    name = "merged"
-  ) |>
-    sf::st_drop_geometry(dp_parcels) |>
-    dplyr::right_join(design_points, by = "id")
-
+    sf::st_transform(sf::st_crs(parcels_vect))
+  dp_vect <- terra::vect(pts_sf)
+  matched <- terra::intersect(dp_vect, parcels_vect)
+  matched_sf <- sf::st_as_sf(matched)
+  dp_with_parcels <- design_points |>
+    dplyr::left_join(matched_sf, by = "id")
   dp_with_parcels
 }