如何解读R语言Terra包crs()输出并获取EPSG代码

1. 解读输出中的多个EPSG编号 #

输出里的EPSG编号分别对应坐标系的不同组成部分：

• EPSG:4269：是基准地理坐标系（NAD83）的代码
• EPSG:9822：是Albers等面积投影方法的代码
• 8821、8822、8823等：是投影参数（假原点纬度、中央经线、标准纬线等）的代码
• EPSG:9001：是长度单位（米）的代码

你的栅格使用的是基于NAD83基准的Albers等面积投影，参数完全匹配EPSG:5070（北美Albers等面积投影），这就是你要找的完整坐标系EPSG代码。

2. 正确用于投影转换的方法 #

方法一：直接使用EPSG代码 #

直接用"EPSG:5070"作为目标坐标系，完全符合terra的要求：

# 将SpatVector转换到目标坐标系
my_spatvector_proj <- terra::project(my_spatvector, "EPSG:5070")

方法二：提取WKT格式坐标系（推荐） #

如果不想手动查EPSG代码，可以直接从栅格中提取WKT格式的坐标系（terra推荐使用WKT而非PROJ.4）：

# 获取栅格的WKT格式坐标系
target_crs <- terra::crs(ST_CDL, proj=TRUE)
# 执行投影转换
my_spatvector_proj <- terra::project(my_spatvector, target_crs)

proj=TRUE参数会返回WKT格式的坐标系字符串，避免PROJ.4格式的兼容性问题。

方法三：验证坐标系匹配 #

可以用以下命令确认栅格坐标系与EPSG:5070是否一致：

# 检查两者坐标系是否匹配
terra::compareCRS(ST_CDL, "EPSG:5070")

返回TRUE则说明两者完全一致，可放心使用。

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