查看原文
其他

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

走天涯徐小洋 走天涯徐小洋地理数据科学 2022-07-17

一种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
  2. 可以白嫖的数据资源!土地覆被、气象、GDP、人口、土壤等等各种空间数据
  3. NetCDF(nc)数据读取与格式转换


您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存