如何在R或QGIS中导入.tpk地图文件并作为shapefile使用？

核心前提说明 #

.tpk是ArcGIS出品的切片瓦片包，本质是预渲染的栅格瓦片集合，本身不包含矢量要素，因此无法直接转换为Shapefile矢量格式，但你提取特定点位坐标的需求可以通过以下方案实现，无需依赖ArcGIS环境。 #

QGIS操作流程 #

• 解压TPK文件：TPK本质是ZIP压缩包，直接将文件后缀修改为.zip，用常规解压工具解压后，会得到包含瓦片配置文件、分级瓦片文件夹的完整目录，其中tilemapresource.xml为瓦片配置核心文件。
• 加载瓦片到QGIS：打开QGIS后直接将解压目录中的tilemapresource.xml拖拽到画布，即可自动完成瓦片图层的加载、坐标系适配；也可通过「图层-添加图层-添加XYZ瓦片图层」手动配置瓦片路径、缩放层级等参数完成加载。
• 提取点位坐标：开启顶部工具栏的「坐标捕捉」功能，使用「识别工具」点击目标点位，即可在弹出的信息面板中查看该点的实时地理坐标；如需批量提取，可新建空白点Shapefile图层，手动标注所有目标点位后导出属性表，即可得到全量点位坐标数据。

R操作流程 #

• 安装依赖包：需用到tiler（瓦片读取）、sf（矢量处理）、terra（栅格处理）三个工具包，安装命令如下：

  install.packages(c("tiler", "sf", "terra"))
  

• 读取瓦片数据：先按上述方法解压TPK文件，调用tiler::tile_info()读取瓦片配置参数，再通过tiler::tile()函数将指定缩放层级的瓦片合并为单张栅格对象导入R环境。
• 提取点位坐标：如果需要交互式选点，调用sf::st_click()点击画布上的目标点位，即可直接返回对应坐标的sf矢量对象，可直接导出为CSV或Shapefile格式；如果已知点位的像素位置，调用terra::xyFromCell()即可直接转换为对应坐标系的地理坐标。

内容的提问来源于stack exchange，提问作者edo