从 GeoPackage 中删除图层

Question

我正在尝试使用 sf 包从 GeoPackage 文件中删除矢量图层。 "delete" 我的意思是永久删除而不是覆盖或更新。我知道 delete_layer 选项，但据我了解，这仅用于在用同名图层替换之前删除图层。

不幸的是，我使用非标准编码将一个名称的图层写入 GeoPackage，这实际上使整个 gpkg 文件在 QGIS 中不可读。因此，我试图找到一种解决方案，通过 R 将其删除。

Answer 1

geopackage 也是一个 SQLite 数据库，因此您可以使用 RSQLite 数据库函数来删除 tables。

设置测试：

> d1 = st_as_sf(data.frame(x=runif(10),y=runif(10),z=1:10), coords=c("x","y"))
> d2 = st_as_sf(data.frame(x=runif(10),y=runif(10),z=1:10), coords=c("x","y"))
> d3 = st_as_sf(data.frame(x=runif(10),y=runif(10),z=1:10), coords=c("x","y"))

将这些写入 GPKG：

> st_write(d1,"deletes.gpkg","d1")
Writing layer `d1' to data source `deletes.gpkg' using driver `GPKG'
features:       10
fields:         1
geometry type:  Point
> st_write(d2,"deletes.gpkg","d2",quiet=TRUE)
> st_write(d3,"deletes.gpkg","d3",quiet=TRUE)

现在要删除，使用RSQLite包（来自CRAN），创建数据库连接：

library(RSQLite)
db = SQLite()
con = dbConnect(db,"./deletes.gpkg")

并删除 table：

dbRemoveTable(con, "d2")

有一个小问题 - 这会删除 table 但不会删除 GPKG 用来说明此包是空间包的元数据 table。因此，您使用 GDAL 工具会收到这样的警告：

$ ogrinfo -so -al deletes.gpkg 
ERROR 1: Table or view 'd2' does not exist
Warning 1: unable to read table definition for 'd2'

QGIS 很高兴地正确读取了剩余的两层。我认为这可以通过在 SQLite 模块旁边加载 Spatialite 模块扩展来在 R 中解决，或者手动删除元数据 tables gpkg_geometry_columns 和 gpkg_ogr_contents 中的行，但似乎没有任何问题那些没有更新的东西很难。

从 GeoPackage 中删除图层

Delete layer from GeoPackage

r

gis

qgis

sf

geopackage