读取/导入 .tpk 映射到 R 或 QGIS 并用作 shapefile

Question

是否可以在 R 或 中导入 .tpk map 文件>QGIS 并将其用作 shapefile？我需要可以使用 shapefile map 提取的 .tpk map 中特定位置的坐标。我无权访问 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 中，但它的大小非常大，我找不到一种方法来根据有效图 block 选择给定的缩放级别和给定的范围。我可以在命令行上使用 gdal_translate 或通过 gdalUtils 包为给定的缩放级别创建 GeoTIFF，使用 USE_BOUNDS=NO将输出限制在仅存在图 block 的地方:

命令行:

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 即可获得缩放的栅格。