R语言实现NetCDF气象数据的逐月解析提取

NetCDF格式数据在气象领域应用较多，处理此类数据的思路或软件有：ArcGIS/QGIS/R语言。

以全球降水量数据集“cru_ts4.04.2001.2010.pre.dat.nc”为例。

ArcGIS：利用“多维工具”进行导入，在“属性-NetCDF”一栏选取对应月份，逐一导出。缺点：来回“重复”，容易卡机。

QGIS：查看方便，导出tiff格式后数据容易丢失。

R语言：方便/快捷/数据无丢失。

代码如下：

#打开nc文件，查看nc文件数据结构
library(ncdf4)
setwd("H:/dataset-derived-near-surface")#路径设置
ncdata <- nc_open("cru_ts4.04.2001.2010.pre.dat.nc")
str(ncdata)

#nc转TIFF
library(raster)
wetbr <- brick("cru_ts4.04.2001.2010.pre.dat.nc", varname="pre")
#nc转多波段TIFF
writeRaster(wetbr, "WorldWet.tif", format="GTiff", overwrite = TRUE) 
#nc转单波段TIFF
writeRaster(wetbr, wetbr@data@names, bylayer=T, format="GTiff", overwrite = TRUE)