st_transform in duckdb creates all-INF geometry points?

[cboettig]

library(dplyr)
library(duckdbfs)
load_spatial()

pad_tract = open_dataset("https://minio.carlboettiger.info/public-biodiversity/pad-us-4/pad-by-tract.parquet")

pad_tract |>
  select(STATE, FIPS, geom) |>
  filter(STATE == State) |> 
  distinct() |>
  mutate(geom = st_transform(geom, "EPSG:4326", "EPSG:3857")) |>
  head() |> 
  to_sf()

all geom points are INF!

But fine with sf:

  pad_tract |>
  select(STATE, FIPS, geom) |>
  filter(STATE == State) |> 
  distinct() |>
  to_sf(crs = "EPSG:4326") |>
  mutate(geom = st_transform(geom, "EPSG:4326", "EPSG:3857")) |>
  mutate(area  = st_area(geom)) |> 
  summarise(total = sum(area))

[cboettig]

Duckdb needs always_xy := true when going to/from 4326, as per docs, https://duckdb.org/docs/extensions/spatial/functions.html#st_transform, but not easy to accomplish this in the R Code translation.

[cboettig]

library(dplyr)
library(duckdbfs)
load_spatial()

pad_tract = open_dataset("https://minio.carlboettiger.info/public-biodiversity/pad-us-4/pad-by-tract.parquet")

x = pad_tract |>
  select(STATE, FIPS, geom) |>
  filter(STATE == "Rhode Island") |> 
  distinct() |>
  mutate(geom = st_transform(geom, 
                             "EPSG:4326", 
                             "EPSG:3857",  
                             sql("always_xy := TRUE") 
  )
  ) |>
  to_sf()