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)