R语言逆地理编码：提取乌干达医疗机构层级地理信息

你当前用reverse_geocode仅得到完整地址字符串，要提取子郡、教区、村庄这类层级信息，有两种可行方案：

方案一：挖掘ArcGIS反向地理编码的结构化字段 #

默认调用reverse_geocode只返回精简结果，开启完整结果后能拿到ArcGIS提供的分层字段，再针对性提取：

1. 调整代码获取完整结果

# 开启full_results参数，获取所有结构化地理字段
healthfacilities_reverse_geocode <- HealthFacilties %>% 
  reverse_geocode(lat = Lat, long = Long, method = "arcgis", full_results = TRUE)

2. 查看并提取目标层级
   先运行colnames(healthfacilities_reverse_geocode)查看所有返回字段，乌干达的行政层级通常对应这些字段（实际以返回结果为准）：

• 子郡（Sub-county）：对应Subregion字段
• 村庄（Village）：对应Neighborhood字段
• 教区（Parish）：若没有直接字段，可从address字符串拆分（乌干达地址格式多为「村庄, 教区, 子郡, 区, 乌干达」）

示例提取代码：

# 替换"FacilityName"为你CSV中的机构名称列名
health_facilities_final <- healthfacilities_reverse_geocode %>%
  select(FacilityName, Subregion, Neighborhood, address) %>%
  # 拆分地址字符串获取教区
  separate(address, into = c("Village", "Parish", "Subcounty", "District", "Country"), 
           sep = ", ", fill = "right") %>%
  # 合并可能重复的字段（比如Subregion和拆分出的Subcounty）
  mutate(Subcounty = ifelse(is.na(Subcounty), Subregion, Subcounty),
         Village = ifelse(is.na(Village), Neighborhood, Village)) %>%
  select(FacilityName, Subcounty, Parish, Village)

方案二：空间连接乌干达行政边界数据（精度更高） #

如果反向地理编码的结果不够准确，直接用空间匹配的方法更靠谱：

1. 加载本地行政边界Shapefile
   提前下载乌干达子郡/教区/村庄级别的边界文件（本地存储），用sf包读取：

library(sf)
# 替换为你的边界文件路径
uganda_admin <- st_read("~/myproject/uganda_village_parish_subcounty.shp")

2. 转换医疗机构数据为空间点对象

# 将经纬度转为sf点格式，坐标系用WGS84（EPSG:4326）
health_facilities_sf <- HealthFacilties %>%
  st_as_sf(coords = c("Long", "Lat"), crs = 4326)

3. 空间匹配获取层级信息

# 用st_within判断点落在哪个行政边界内，完成连接
health_facilities_admin <- st_join(health_facilities_sf, uganda_admin, join = st_within)

之后直接从health_facilities_admin中提取子郡、教区、村庄等字段即可。

注意事项 #

• 确保医疗机构经纬度和行政边界数据的坐标系一致（推荐用WGS84）
• 若拆分地址字符串，需先观察几条地址的格式，调整separate函数的分隔符和字段顺序

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