基于R语言tmap创建可切换图层与指标的交互式GIS地图问题咨询

需求目标 #

• 支持选择单个指标进行查看
• 可切换查看地理区域（整合护理委员会ICB）及其下属子群体（临床委托小组CCG）的对应指标数据
  示例场景：先查看牛津周边CCG的糖尿病患病率，再切换至所属ICB区域的同指标数据；或先查看约克郡ICB的平均手术等待时长，再切换到下属CCG的对应数据

当前实现与问题 #

基于自定义Shapefile制作英格兰医疗地理图层（ICB和CCG），搭配两个虚拟指标数据，当前代码如下：

tm_shape(england_ICB_maps) + 
  tm_polygons(col=c('Metric 1', 'Metric 2')) +  
  tm_shape(england_CCG_maps) + 
  tm_polygons(col=c('Metric 1', 'Metric 2')) + 
  tm_facets(as.layers = TRUE)

存在的问题：

• 两个指标分别对应ICB、CCG两个图层，选择指标时会同时显示两个层级的图层，无法单独选择查看某一层级（例如只看ICB的Metric 1，或只看CCG的Metric 2）

tm_view参数配置困惑 #

地图基于Leaflet构建，需要配置底图（如默认设置为OpenStreetMap并移除其他选项），但tm_view的leaflet.options参数文档说明不够直观：

通过leafletOptions传递给leaflet.js地图创建的其他选项（参见leaflet文档，Map>Creation）。命名列表，名称对应变量名。提示：使用zoomSnap和zoomDelta实现分数缩放。

尝试添加以下代码时出现报错：

+ tm_view(leaflet.options = addLayersControl(baseGroups = "OpenStreetMap"))

错误信息：

Error in addLayersControl(baseGroups = "OpenStreetMap") :
could not find function "addLayersControl"

解决方案 #

1. 实现层级与指标的独立选择 #

要实现ICB/CCG层级和指标的独立控制，需要调整tmap的图层分组逻辑，为每个层级的指标单独命名，再通过配置实现Leaflet图层控件的独立切换。示例代码思路：

# 定义ICB层级的两个指标图层，单独命名
icb_m1 <- tm_shape(england_ICB_maps) +
  tm_polygons(col = "Metric 1", title = "ICB - 指标1")
icb_m2 <- tm_shape(england_ICB_maps) +
  tm_polygons(col = "Metric 2", title = "ICB - 指标2")

# 定义CCG层级的两个指标图层，单独命名
ccg_m1 <- tm_shape(england_CCG_maps) +
  tm_polygons(col = "Metric 1", title = "CCG - 指标1")
ccg_m2 <- tm_shape(england_CCG_maps) +
  tm_polygons(col = "Metric 2", title = "CCG - 指标2")

# 切换到交互式视图模式
tmap_mode("view")

# 组合图层并配置控件
icb_m1 + icb_m2 + ccg_m1 + ccg_m2 +
  tm_basemap("OpenStreetMap") +  # 设置默认底图
  tm_view(
    # 配置Leaflet地图初始化选项
    leaflet.options = list(
      zoomSnap = 0.5,  # 允许分数缩放的示例配置
      zoomDelta = 0.5
    )
  )

这样配置后，Leaflet控件会显示所有四个独立图层，你可以自由选择查看任意层级的任意指标。

2. tm_view的leaflet.options正确用法示例 #

leaflet.options接受命名列表，对应Leaflet.js地图初始化的底层配置选项，不能直接调用Leaflet R包的函数（比如addLayersControl）。常见配置示例：

• 开启分数缩放：

tm_view(leaflet.options = list(zoomSnap = 0.5, zoomDelta = 0.5))

• 禁用地图缩放控件：

tm_view(leaflet.options = list(zoomControl = FALSE))

• 设置地图最大缩放级别：

tm_view(leaflet.options = list(maxZoom = 12))

• 禁用地图拖拽：

tm_view(leaflet.options = list(dragging = FALSE))

3. 底图配置补充 #

如果需要仅保留OpenStreetMap作为底图选项，可使用tm_basemap指定唯一底图，替代默认的多底图选择：

tm_basemap("OpenStreetMap") +
  # 后续图层代码...

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