sp包调用CRS函数传入EPSG坐标系代码返回NA报错问题求助

你遇到的报错是sp 1.4-5版本与GDAL 3.2+、PROJ 7.2+版本的兼容性问题导致的，旧版本sp的EPSG数据库读取逻辑没有适配高版本PROJ的搜索路径，所以传入EPSG编码的调用会返回NA报错。

解决方案1：升级核心空间包（优先推荐） #

直接升级sp、rgdal、sf等核心空间依赖包到最新版本，适配高版本GDAL/PROJ：

• 全平台通用升级命令：

install.packages(c("sp", "rgdal", "sf"), type = "source")

• Windows用户需要提前安装对应R版本的Rtools工具链才能完成源码编译，Mac/Linux用户需提前安装系统级的GDAL/PROJ开发依赖。
• 升级完成后重启R会话，测试sp::CRS(SRS_string='EPSG:4326')即可正常返回CRS对象。

解决方案2：临时兼容调用（无需升级包） #

如果暂时无法升级依赖包，可以通过sf包中转构造CRS对象，避免直接调用sp的CRS接口报错：

# 用sf读取EPSG编码后转为sp兼容的CRS对象
library(sf)
crs_4326 <- as(st_crs("EPSG:4326"), "CRS")

自定义投影警告处理 #

你传入proj4字符串触发的警告是高版本PROJ的正常提示：PROJ 6+版本后已经废弃proj4字符串作为权威坐标系的唯一标识，自定义投影没有关联官方datum时就会弹出该警告，不影响实际坐标转换运算。如果需要消除警告，可以改用WKT格式定义坐标系，示例如下：

# 替换为你自定义投影对应的WKT字符串即可
custom_crs <- sp::CRS(WKT = '你的自定义投影WKT文本')

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