reprex for bug in files with multiple layers and SRS in terra (gpkg example)

test-multi-CRS.R

library(terra)
gdal(lib = "")
data(us_states, package = "spData")

x <- vect(us_states)[1,]
y <- project(x, "EPSG:5070")

writeVector(x, "test.gpkg")
writeVector(y, "test.gpkg", layer = "test2", insert = TRUE)

# CRS is as expected
a <- vect("test.gpkg")
a
plot(a)

# CRS is wrong (matches first layer)
b <- vect("test.gpkg", layer = "test2")
b
plot(b)

# inspect gpkg contents SRS 
con <- DBI::dbConnect(RSQLite::SQLite(), "test.gpkg")
DBI::dbGetQuery(con, "SELECT * FROM gpkg_contents")
DBI::dbGetQuery(con, "SELECT * FROM gpkg_spatial_ref_sys")
DBI::dbDisconnect(con)