在R中通过经纬度关联表格与Shapefile的方法求助

1. 安装并加载必要的R包 #

先安装处理空间数据和Excel的工具包，新手直接运行以下代码：

install.packages(c("sf", "readxl", "dplyr"))

加载已安装的包：

library(sf)
library(readxl)
library(dplyr)

2. 导入你的Excel表格 #

假设你的Excel文件名为bike_stations.xlsx，且包含经度和纬度列，用read_excel导入数据：

bike_df <- read_excel("bike_stations.xlsx")

3. 将Excel数据转换为空间点对象 #

Shapefile是多边形要素，需要把经纬度数据转换成空间点要素，且必须和Shapefile的坐标参考系统（CRS）保持一致：

• 先查看Shapefile的CRS信息：

# 替换bg_shape为你实际导入的Shapefile对象名
st_crs(bg_shape)

输出会显示CRS编码（比如EPSG:4269，这是美国常用的地理坐标系）。

• 转换Excel数据为sf点对象：

bike_points <- st_as_sf(bike_df, 
                        coords = c("经度", "纬度"),  # 替换为你Excel里的经纬度列名
                        crs = st_crs(bg_shape))  # 和Shapefile使用相同CRS

4. 空间匹配获取GEOID（FIPS编码） #

用st_join做空间连接，将每个单车停靠点匹配到对应的多边形，提取GEOID：

# 左连接：保留所有停靠点数据，未匹配到的GEOID会显示为NA
bike_with_fips <- st_join(bike_points, bg_shape[, "GEOID"], join = st_within)

这里bg_shape[, "GEOID"]仅保留Shapefile的GEOID列，减少计算量；st_within用于判断点是否在多边形内部，是最常用的匹配规则。

5. 导出结果（可选） #

如果需要把结果转回普通数据框并导出为Excel：

# 移除空间属性，转换为普通数据框
bike_result <- st_drop_geometry(bike_with_fips)

# 导出到Excel文件
write_excel_csv(bike_result, "bike_stations_with_fips.xlsx")

注意事项 #

• 确认Excel中的经纬度是十进制格式（如120.1234，而非度分秒格式）
• 如果匹配结果出现大量NA，检查经纬度列的顺序是否搞反，或CRS是否完全一致
• 如果Shapefile使用的是投影坐标系（如UTM），可先将其转换为地理坐标系（EPSG:4326）再进行匹配

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