R语言合并带crs栅格时projection返回NA报错如何解决

该问题通常由三类原因导致：

• 栅格对象属于terra包的SpatRaster类，调用了raster包的projection()接口，跨包函数不兼容导致返回NA，打印对象展示的是内置元数据而非实际读取到的投影信息
• 赋值的PROJ.4字符串存在语法错误、或包含当前PROJ版本不支持的废弃参数，导致CRS未实际写入对象
• 此前用sapply处理栅格列表时，返回结果为矩阵而非栅格对象列表，导致投影赋值未生效

首先运行class(raster1)确认栅格对象所属的类，再对应执行操作：

场景1：对象为raster包的RasterLayer/RasterStack/RasterBrick类 #

1. 先验证待设置的CRS有效性，优先使用EPSG代码替代PROJ.4字符串提升兼容性

# 替换为你实际需要的CRS，可替换为EPSG代码格式如"EPSG:4326"
target_crs <- "+proj=utm +zone=50 +datum=WGS84 +units=m +no_defs"
# 校验CRS合法性，无报错则为有效
crs(target_crs)

2. 批量为栅格列表设置投影

# 用lapply返回栅格列表，保留栅格对象结构
rasters <- lapply(rasters, function(x) {
  crs(x) <- target_crs
  return(x)
})
# 验证投影是否设置成功
sapply(rasters, function(x) crs(x, asText = TRUE))

3. 正确调用merge函数合并栅格
   merge函数要求首个参数为单个栅格对象，后续传入其他待合并栅格，直接传入列表会触发参数缺失报错，用do.call拆分列表参数即可：

merged_raster <- do.call(raster::merge, rasters)

场景2：对象为terra包的SpatRaster类 #

不要使用raster包的projection()相关接口，使用terra包原生函数操作：

1. 批量设置投影

# 优先用EPSG代码设置CRS
target_crs <- "EPSG:32650"
rasters <- lapply(rasters, function(x) {
  terra::crs(x) <- target_crs
  return(x)
})
# 验证投影
sapply(rasters, terra::crs)

2. 合并栅格

merged_raster <- do.call(terra::merge, rasters)

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