如何使用proj4string将dataframe转为同名SpatialPointsDataFrame

你可以直接提取California_shp的投影信息复用，注意构造空间点时坐标列要按先经度、后纬度的顺序传入，完整实现代码如下：

# 加载sp依赖包（已加载可跳过该步）
library(sp)

# 提取经纬度坐标列，顺序必须为经度在前、纬度在后
houses_coords <- Houses[, c("longitude", "latitude")]

# 构造SpatialPointsDataFrame，直接复用已有投影
Houses <- SpatialPointsDataFrame(
  coords = houses_coords,
  data = Houses, # 保留原数据框所有属性字段
  proj4string = California_shp@proj4string # 直接复用California_shp的投影参数
)

你可以运行以下代码验证转换结果：

# 查看转换后的Houses对象投影
Houses@proj4string

运行后输出的投影参数会和California_shp完全一致，且Houses会保留原有所有数据字段，同时具备空间点对象的空间运算能力。

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