一种NetCDF数据快速裁剪转TIFF方法

一种NetCDF数据快速裁剪方法

NetCDF即nc数据，经常用于气象、经济、土地覆被等栅格数据的网络传输和共享。一般做空间分析的时候我们需要将nc数据转为TIFF格式，并进行裁剪。

我一般使用R语言进行nc数据的读取、裁剪和转换，代码非常简单，但是R语言在裁剪栅格的时候效率非常低下。QGIS可以直接打开nc数据，但是直接用掩膜裁剪工具裁剪nc会出错，怎么办呢？

以全球公里格网GDP数据为例进行介绍：

R语言nc数据读取裁剪转TIFF

先看看R语言

• 数据源：https://datadryad.org/stash/dataset/doi:10.5061/dryad.dk1j0
• GDP_PPP_30arcsec_v3.nc对这个数据进行裁剪和转换

library(terra)
gdp = rast("D:/GISDATA/GDP/GDP_PPP_30arcsec_v3.nc")
henanSHP = vect("D:/R/CMR/SHP/China.shp")  #研究区SHP
gdp_hn = trim(mask(gdp, henanSHP))   #裁剪
writeRaster(gdp_hn, filename ="./GDP/Henan.tiff", names=TempCMR@ptr[["names"]])

这个GDP数据像元数非常多

漫长的等待……

数据裁剪了一晚上也没结果……

ArcGIS是32位的软件，处理大数据量栅格能力有限，而且不能直接打开nc文件，需要使用ArcToolbox进行转换，太麻烦了，这个方法直接PASS

QGIS裁剪nc数据

QGIS可以直接打开nc数据，但是按掩膜图层裁剪栅格会出问题

按掩膜图层裁剪栅格

报错

怎么办呢？

QGIS快速裁剪NC并转TIFF

方法非常简单粗暴有效，在nc数据图层上右击，导出，另存为

直接另存为，选择地图画布范围

这样就通过另存为的方法，将nc裁剪为地图显示范围，转换为了TIFF文件，用时几秒钟！

数据裁剪结果

公众号后台回复：GDP获取数据下载方法

参考文献

1. https://datadryad.org/stash/dataset/doi:10.5061/dryad.dk1j0