如何对120个netCDF日均温文件求平均生成单张栅格图

一、R方案（推荐，高效处理缺失值+投影） #

用terra和ncdf4包可以轻松搞定批量NetCDF处理，完美解决缺失值和投影问题：

1. 安装并加载依赖 #

install.packages(c("terra", "ncdf4"))
library(terra)
library(ncdf4)

2. 批量读取所有NetCDF文件 #

把所有mintemp_*.nc放在同一个文件夹里，执行以下代码：

# 替换为你的文件路径
nc_path <- "/Users/xxx/your_nc_files/"
nc_files <- list.files(nc_path, pattern = "mintemp_.*\\.nc", full.names = TRUE)

# 一次性读取所有文件为栅格栈
temp_stack <- rast(nc_files)

3. 计算像素级多年平均（自动忽略缺失值） #

mean()函数默认跳过NA，只算有效数据的平均值，完全不用担心里缺失值拖后腿：

annual_mean <- mean(temp_stack, na.rm = TRUE)

4. 导出带正确投影的TIFF #

导出的文件会保留原始NetCDF的空间参考，ArcGIS直接就能识别：

writeRaster(annual_mean, "australia_120y_mintemp_mean.tif", format = "GTiff", overwrite = TRUE)

二、ArcGIS Pro方案（适配熟悉工具） #

如果习惯用ArcGIS，重点处理缺失值和投影两个问题：

1. 创建镶嵌数据集批量导入NetCDF #

• 在Catalog面板找到所有NetCDF文件，右键创建镶嵌数据集，把所有文件的波段都加进去
• 镶嵌数据集会自动统一空间参考，避免投影混乱

2. 配置缺失值规则 #

• 右键镶嵌数据集 → 属性 → 处理标签
• 找到NoData设置，填入原始数据中缺失值的编码（比如-9999或NaN），勾选计算时忽略NoData

3. 计算平均值 #

• 用单元格统计工具，统计类型选“平均值”，勾选“忽略NoData”
• 或者用栅格计算器，输入公式：MEAN("你的镶嵌数据集名称")

4. 修复投影异常 #

如果之前CDO导出的TIFF投影不对：

• 用定义投影工具给TIFF指定正确坐标系（比如GDA2020或WGS84）
• 再用投影栅格转换到目标坐标系

三、CDO异常问题排查 #

你之前遇到的230℃离谱值，90%是单位错了——原始数据可能是开尔文，被当成摄氏度了。先转单位再计算：

# 把开尔文转成摄氏度
cdo subc,273.15 mintemp_2000.nc mintemp_2000_c.nc
# 批量转所有文件（Mac终端用循环）
for file in mintemp_*.nc; do cdo subc,273.15 "$file" "${file%.nc}_c.nc"; done
# 计算多年平均
cdo ensmean mintemp_*_c.nc annual_mean_c.nc
# 用gdal导出带投影的TIFF
gdal_translate -of GTiff annual_mean_c.nc annual_mean_c.tif

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