在R语言中利用GADM数据根据坐标查询法国一级行政区

可以通过空间匹配提取该坐标对应的一级行政区，具体实现代码如下：

library(raster)
# 加载法国一级行政区GADM数据（已执行）
fr = getData("GADM", country="FRA", level=1)

# 目标点坐标
latitude = 48.86122650866868
longitude = 2.341541835915652

# 创建空间点对象，确保坐标参考系与GADM数据一致（默认WGS84）
point = SpatialPoints(data.frame(x = longitude, y = latitude), proj4string = crs(fr))

# 执行空间匹配，提取对应行政区名称
target_admin = over(point, fr)$NAME_1
print(target_admin)

执行后会输出：

[1] "Île-de-France"

对应的中文名称为法兰西岛。

关键注意事项 #

• 必须保证点的坐标参考系（CRS）与GADM图层完全一致，getData()获取的GADM数据默认采用WGS84（EPSG:4326），创建点时指定相同CRS即可避免匹配错误。
• over()函数用于实现点与多边形的空间关联，返回点所在多边形的属性数据，这里提取NAME_1字段即为一级行政区名称。

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