transforming error before cropping

[COverton-CWA]

This may be a system issue on my end involving a recent QGIS install as my default PROJ_LIB and GDAL_DATA environment variables originally pointed to a postgres install rather than rtools but...

I can't transform vectors to the abundance raster crs:

> median_abund<-load_raster("Abert's Towhee", product = "abundance", metric = "median", resolution = "9km") 
> northamerica <- st_read(system.file("shapes/world.gpkg", package="spData")) %>% filter(continent == "North America")
> northamerica<-northamerica %>% st_transform(st_crs(median_abund))

results in

Warning in CPL_transform(x, crs, aoi, pipeline, reverse, desired_accuracy,  :
  GDAL Error 6: Cannot find coordinate operations from `{
  "$schema": "https://proj.org/schemas/v0.7/projjson.schema.json",
  "type": "ProjectedCRS",
  "name": "unknown",
  "base_crs": {
    "type": "GeographicCRS",
    "name": "unknown",
    "datum": {
      "type": "GeodeticReferenceFrame",
      "name": "World Geodetic System 1984",
      "ellipsoid": {
        "name": "WGS 84",
        "semi_major_axis": 6378137,
        "inverse_flattening": 298.257223563
      },
      "id": {
        "authority": "EPSG",
        "code": 6326
      }
    },
    "coordinate_system": {
      "subtype": "ellipsoidal",
      "axis": [
        {
          "name": "Latitude",
          "abbreviation": "lat",
          "direction": "north",
          "unit": "degree"
        },
        {
          "name": "Longitude",
          "abbreviation": "lon",
          "direction": "east",
          "unit": "degree"
        }
      ]
    }
  },
  "conversion": {
    "name": "unknown",
    "method": {

Error: OGRCreateCoordinateTransformation(): transformation not available

as does

northamerica<-exact_extract(median_abund, northamerica, 'sum')

[mstrimas]

That's strange, this code works fine for me, but I am on a Mac and I'm guessing you're on windows. Can you try running this and tell me what you get:

> terra::gdal()
[1] "3.8.5"
> st_crs(median_abund)$proj4string
[1] "+proj=eqearth +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs"
> st_crs(median_abund)$epsg
[1] 8857

[COverton-CWA]

Definitely, something amiss.

> terra::gdal()
[1] "3.11.4"
> st_crs(median_abund)$proj4string
[1] NA
> st_crs(median_abund)$epsg
[1] NA
> st_crs(median_abund)
Coordinate Reference System:
  User input: WGS 84 / Equal Earth Greenwich 
  wkt:
ENGCRS["WGS 84 / Equal Earth Greenwich",
    EDATUM["Unknown engineering datum"],
    CS[Cartesian,2],
        AXIS["(E)",east,
            ORDER[1],
            LENGTHUNIT["metre",1,
                ID["EPSG",9001]]],
        AXIS["(N)",north,
            ORDER[2],
            LENGTHUNIT["metre",1,
                ID["EPSG",9001]]]]

[mstrimas]

Sorry, I forgot to respond to this. I looks like it's something on the GDAL side, which is tricky for me to troubleshoot since it's likely totally specific to your system and I'm not too familiar with the details of GDAL. If you have access to a different computer, you might try checking if your code runs there, or maybe uninstall and re-install GDAL? Sorry I don't have any better solutions...