导入Shapefile到RStudio时PROJ.db未找到，EPSG坐标返回NA求助

问题原因 #

该错误表明sf包无法定位PROJ库的核心配置文件proj.db，导致坐标系解析失效，最终st_crs()无法返回有效EPSG代码。

解决方案 #

1. 设置PROJ_LIB环境变量 #

先找到proj.db所在的目录：

• Windows：通常在C:\Program Files\PROJ\share\proj（独立安装）或C:\OSGeo4W\share\proj（OSGeo4W安装）
• macOS：Homebrew安装路径为/opt/homebrew/share/proj，系统默认路径为/usr/share/proj
• Linux：Ubuntu/Debian为/usr/share/proj，CentOS/RHEL为/usr/share/proj

在R中临时设置环境变量（重启R后失效）：

Sys.setenv(PROJ_LIB = "/path/to/folder/containing/proj.db")

若需永久生效，配置系统级环境变量：

• Windows：右键「此电脑」→ 属性 → 高级系统设置 → 环境变量，新建PROJ_LIB变量并填入路径
• macOS/Linux：在~/.bashrc或~/.zshrc中添加export PROJ_LIB=/path/to/proj/folder，执行source ~/.bashrc生效

2. 验证依赖库版本兼容性 #

运行以下命令查看sf依赖的GDAL、PROJ版本：

sf::sf_extSoftVersion()

若版本不匹配或存在缺失，重新安装对应库后重装sf：

install.packages("sf", type = "source")

（Windows/macOS用户优先尝试CRAN二进制包，安装速度更快）

3. 手动指定Shapefile坐标系 #

若Shapefile的.prj文件损坏或格式不标准，可手动指定坐标系：

# 示例：指定为WGS84（EPSG:4326），请根据实际数据替换
st_crs(espana) <- 4326

# 或使用PROJ字符串指定
st_crs(espana) <- "+proj=longlat +datum=WGS84 +no_defs"

执行st_crs(espana)确认坐标系已生效。

4. 重装完整的GDAL/PROJ工具链 #

若上述方法无效，重新安装完整GIS工具链：

• Windows：使用OSGeo4W安装器，选择「Advanced Install」，勾选gdal、proj相关组件
• macOS：通过Homebrew安装：brew install gdal proj
• Linux：Ubuntu/Debian执行sudo apt-get install libgdal-dev libproj-dev；CentOS/RHEL执行sudo yum install gdal-devel proj-devel

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