用R处理nc文件

1.nc文件概述

NetCDF(network Common Data Form)网络通用数据格式是一种面向数组型并适于网络共享的数据的描述和编码标准。目前，NetCDF广泛应用于大气科学、水文、海洋学、环境模拟、地球物理等诸多领域。用户可以借助多种方式方便地管理和操作 NetCDF 数据集。

NetCDF全称为network Common Data Format，中文译法为“网络通用数据格式”；netcdf文件开始的目的是用于存储气象科学中的数据，现在已经成为许多数据采集软件的生成文件的格式。

•从数学上来说，netcdf存储的数据就是一个多自变量的单值函数。用公式来说就是f(x,y,z,…)=value；

•函数的自变量x,y,z等在netcdf中叫做维(dimension) 或坐标轴(axix)，

•函数值value在netcdf中叫做变量(Variables).

一个Netcdf文件的结构包括以下对象：

•变量(Variables) ：变量对应着真实的物理数据。

•维(dimension)：一个维对应着函数中的某个自变量，或者说函数图象中的一个坐标轴，在线性代数中就是一个N维向量的一个分量。

•属性(Attribute) ：属性对变量值和维的具体物理含义的注释或者说解释。

参考文献：气象数据处理——nc文件

2.R处理nc文件

以气候数据为例。

# 处理npp的数据
tmpfile = nc_open("F:\\geodata\\19气温与降水\\气温\\China_Tmp_1901_2023_from_CRUv4.08\\tmp_2000\\tmp_2000.nc") # 查看nc文件的信息
print(tmpfile) time <- ncvar_get(tmpfile, "time") # 获取时间变量
head(time) # 查看时间变量的前几行
print(time) # 查看时间变量# nc文件的数据存成栅格数据集
tmp = brick("F:\\geodata\\19气温与降水\\气温\\China_Tmp_1901_2023_from_CRUv4.08\\tmp_2000\\tmp_2000.nc",varname='tmp')# 读取nc文件中的变量'time# 对栅格数据集进行栅格计算
tmp_mean = calc(tmp, fun = mean, na.rm = TRUE) * 0.1writeRaster(tmp_mean,'F:\\geodata\\19气温与降水\\气温\\China_Tmp_1901_2023_from_CRUv4.08\\tmp_2000\\tmp2000.tif')