01 前言
先起个头!
我们今天要介绍的是ArcGIS处理多维栅格NC格式的全球土壤数据集,提取不同深度的TC(总碳)。+
介绍这个内容主要是,许多学员朋友在处理多维栅格的时候提示创建NC栅格提示一个或全部两个维度的坐标间距是可变的问题:
官方也给出了解决方案:
说明
x 坐标和/或 y 坐标之间的距离不等。大多数情况下,当纬度之间的距离不相等时会遇到此问题。纬度之间的差异向两极递减。
解决方案
要根据 NetCDF变量创建 NetCDF栅格图层,则各个 x 坐标之间的间距必须相等,同时各个 y 坐标之间的间距也必须相等。如果各坐标的间距不等,则需创建一个netCDF要素图层,然后插值为栅格。
但是,很多朋友看不懂这个解决方案该如何处理,那我们今天就详细介绍一下,也会配上视频介绍。
一、下面先介绍一下前期的推文:
要介绍的是前期推文《ArcGIS处理NetCDF(.nc)的多维科学数据》的续集,我们在《ArcGIS处理NetCDF(.nc)的多维科学数据》与《ArcGIS处理多维栅格nc数据升级技巧》专题介绍了以下7个内容:
1、认识NetCDF格式的数据
2、ArcGIS中添加与显示NetCDF数据
3、在特定位置提取每个月降水记录并制图表
4、模型批量提取每月的降水栅格数据
6、ArcGIS处理多维栅格nc数据升级技巧
5、NetCDF_time_slice_to_Raster.zip(zip,4.4MB)
6、GLDAS下载nc4数据
7、ArcGIS处理nc4(调用及批量提取
二、介绍一下今天的全球土壤数据集:
下载数据的网址如下:
http://globalchange.bnu.edu.cn/research/soilwd.jsp
进来是英文网站,不过应该是中国人做的,好像是中山大学的,我们选择翻译一下网页吧。
对于数据我不做过多介绍,大家自己可以去网站看看。
我这边下载了一个5分钟版的TC总碳的数据:
02 ArcGIS中处理全球土壤数据集nc
在ArcGIS中直接利用创建NetCDF栅格图层来处理数据,由于一个或者两个维度间距可变直接报错。直接拖拉数据时,内容列表也只会显示一个点图层。
官方给出了解决方案:
需创建一个netCDF要素图层,然后插值为栅格
1、 那么,我们就按照方案,直接利用创建NetCDF要素层处理一下数据。(数据是有8个不同的深度的数据,每次选择一个深度的数据进行处理)
2、处理好后就会生成一个将近200万的点图层数据。
3、现在就可以对我们得到的点图层利用IDW等插值方法进行插值处理就可以得到全球的土壤数据。但是你会发现说你插值的结果是覆盖全球的,包括了海洋,我们点本来只是陆地的,毕竟是土壤数据集。
4、所以我们在插值工具的环境中加入全球陆地的矢量数据作为掩膜数据。
加入了掩膜的插值数据(南极洲虽然没有点数据,这边先不处理就放那),下图就是我么获取的结果了,然后其他不同深度的数据,可以同样处理哦。
03 FME处理全球土壤数据
如果我们用FME来处理一下这个数据就十分方便了
运行之后,就可以直接提取出8个深度的土壤数据,存储为8个波段的数据。
04 视频教学
ArcGIS处理NetCDF(.nc) 的多维科学数据——6解决提示一个或全部两个维度的坐标间距可变导致创建nc失败