Reading/importing .tpk 映射到 R 或 QGIS 并用作 shapefile

Question

是否可以在 R[= 中导入 .tpk 地图文件26=] 或 QGIS 并将其用作 shapefile？我需要 .tpk 地图 中特定位置的坐标，可以使用 shapefile 地图提取这些坐标。我无法访问 ArcGIS。

感谢任何指导！

有关 .tpk 的更多信息，请参阅：https://pro.arcgis.com/en/pro-app/latest/help/sharing/overview/tile-package.htm

Answer 1

这是我所做的：

首先使用 tpk 转换实用程序将 tpk 文件转换为 mbtiles 来自 https://github.com/consbio/tpkutils

tpk export mbtiles 04010588800801.tpk 04010588800801.mbtiles

gdalinfo 显示了有关 mbtiles 文件的一些信息：

$ gdalinfo 04010588800801.mbtiles 
Driver: MBTiles/MBTiles
Files: 04010588800801.mbtiles
Size is 23418224, 20166662
...
Metadata:
  ZOOM_LEVEL=17
...
  minzoom=0
  maxzoom=17
...

我们可以将它加载到 R 中，但它的大小太大了，我找不到一种方法来 select 给定的缩放级别和基于有效图块的给定范围。我可以在命令行上使用 gdal_translate 或通过 gdalUtils 包为给定的缩放级别创建 GeoTIFF，使用 USE_BOUNDS=NO 将输出限制在仅存在图块的地方：

命令行：

gdal_translate -oo ZOOM_LEVEL=17 -oo USE_BOUNDS=NO 04010588800801.mbtiles zoom17.tiff

gdalUtils 包裹：

gdal_translate("04010588800801.mbtiles","z17.tiff",
     oo=c("USE_BOUNDS=NO","ZOOM_LEVEL=17"))

读取并绘制此 17 级 RGB 图像即可完成：

> z17 = raster::stack("z17.tiff")
[ignore CRS warnings...]
> plotRGB(z17)

请注意，这是一张非常高分辨率的图像，因此您无法阅读标签，但如果您放大（或加载到 QGIS 并在那里进行交互），您可以阅读标签。这是 QGIS 中 17 级的极端放大显示分辨率的限制：

请记住，这只是光栅图像数据，因此如果您需要这些点的坐标，则必须在 QGIS 中创建一个新图层并在图像上手动创建一个点数据集。如果这是您想要的，那么我强烈建议您尝试从供应商处获取矢量数据，而不必做所有这些！

其他缩放级别可能对您有用，因此请使用上面的 gdal_translate 过程转换它们。当你出去时你会失去细节，而在 13 级以下你只会得到概览图。

16 级：

13 级：

12 级：

更新

从.mbtiles 文件中读取的一种稍微更直接的方法。使用允许您传递 GDAL 选项的 stars 包：

s = stars::read_stars("04010588800801.mbtiles",
    options=c("USE_BOUNDS=NO","ZOOM_LEVEL=17"), proxy=FALSE)
rs = as(s,"Raster")
raster::plotRGB(rs)

需要 proxy=FALSE，否则当使用 as(..,"Raster") 转换时，输出将恢复为完整全局 17 级缩放栅格的 20166662x23418224 维度。某个地方的星星可能存在错误。无论如何，这让你无需在任何地方使用 gdal_translate 就可以得到缩放的光栅。

Answer 2

基于@Spacedmans 的回答并使用他创建的 tpk 文件，我使用 terra 包

显示了一系列缩放级别

library(terra)
plotm <- function(zoom, e=NULL) {
    # read file with GDAL options
    r = rast("04010588800801.mbtiles", opts=c("USE_BOUNDS=NO", paste0("ZOOM_LEVEL=",zoom))); 
    # declare the RGB channels
    RGB(r) <- 1:4
    plot(r, ext=e)
}
par(mfrow=c(3,3))
x <- lapply(6:14, plotm)

这里使用范围放大一点，以便能够更好地看到最高分辨率下的差异

par(mfrow=c(2,2))
e <- ext(3935206, 3935662, 1085835, 1086294)
x <- lapply(14:17, plotm, e=e)

这在 ArcMap 中看起来与原始 tpk 文件相同。

Reading/importing .tpk 映射到 R 或 QGIS 并用作 shapefile

Reading/importing .tpk maps into R or QGIS and use as shapefile

r

arcgis

qgis

sf

更新