如何通过RSQLite读取MBTiles文件tile_data列Blob内的要素ID

核心思路 #

直接解析RSQLite读取到的tile_data二进制Blob（矢量瓦片一般为Google Protobuf格式的pbf文件），无需通过sf全量读取MBTiles，可保留原始xyz瓦片坐标信息，且可按需仅提取要素ID不加载几何数据。

依赖包安装 #

需要安装protolite包用于解析MVT格式的矢量瓦片，无需GDAL依赖：

install.packages("protolite")

完整实现代码 #

library(RSQLite)
library(protolite)
library(dplyr)

# 1. 连接MBTiles并读取瓦片数据
con <- dbConnect(SQLite(), "替换为你的MBTiles文件路径.mbtiles")
# 全量读取可删除limit限制
tiles_raw <- dbGetQuery(con, "select zoom_level, tile_column, tile_row, tile_data from tiles limit 10")
dbDisconnect(con)

# 2. 定义单瓦片解析函数，仅提取要素ID
extract_feature_ids <- function(tile_blob) {
  # 解析MVT格式pbf文件
  mvt_content <- read_mvt(tile_blob)
  # 提取所有图层的全部要素ID，可按需调整为仅读取指定图层
  all_ids <- lapply(mvt_content$layers, function(layer) {
    layer$features$id
  })
  return(unlist(all_ids, use.names = FALSE))
}

# 3. 批量处理所有瓦片，保留xyz关联关系
result <- tiles_raw %>%
  rowwise() %>%
  mutate(feature_id_list = list(extract_feature_ids(tile_data))) %>%
  ungroup()

# 可选：将列表格式的ID展开为一行一个ID的长表格式
# result_long <- result %>% unnest(feature_id_list)

注意事项 #

• 如果你的矢量瓦片要素ID没有存储在pbf原生ID字段，而是存放在要素属性中，可将提取逻辑调整为layer$features$properties$你的ID字段名，根据实际数据结构修改即可。
• MBTiles默认存储的tile_row为TMS规范的坐标，如需转换为通用XYZ规范的行号，可新增列计算：tile_row_xyz = (2^zoom_level - 1) - tile_row。

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