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

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

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

ArcGIS:利用“多维工具”进行导入,在“属性-NetCDF”一栏选取对应月份,逐一导出。

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

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

 

代码如下:

#打开nc文件,查看nc文件数据结构
library(ncdf4)
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)

 

上一篇:升级到macOS Big Sur后,网络实用工具没有了怎么办?


下一篇:Python安装netCDF4