Error in CPL_get_z_range(obj, 1) : z error - expecting three columns;

[rhijmans]

I see this error

f = "/vsicurl/https://github.com/GEMScienceTools/gem-global-active-faults/raw/master/geopackage/gem_active_faults_harmonized.gpkg"
 
library(terra)
#terra 1.7.6
v <- vect(f)

library(sf)
#Linking to GEOS 3.9.4, GDAL 3.6.0, PROJ 9.1.1; sf_use_s2() is TRUE
s <- read_sf(f)
#Error in CPL_get_z_range(obj, 1) : z error - expecting three columns;

It could be that the file has an issue. "terra" silently drops the z-coordinates; so it may be that it can read the file for that reason.

[edzer]

This has been solved in the dev version, #2046 #1592

It returns a mix of dimensions:

> s <- read_sf(f)
> s
Simple feature collection with 13696 features and 23 fields
Geometry type: LINESTRING
Dimension:     XY, XYZ
Bounding box:  xmin: -180 ymin: -66.163 xmax: 180 ymax: 86.805
z_range:       zmin: NA zmax: NA
Geodetic CRS:  WGS 84