ArcMap Zonal Statistic按区域shp统计功能的使用与error的处理
工具
ArcMap 10.2
材料
1、一份tif 格式 Raster数据
2、中国地理分区shp数据
(最近发现一个很好的网站,是中科院地理科学与资源研究所的一个数据中心,数据非常全面,包括行政区划、地理区划、气象数据、土地利用数据、植被类型数据等等,只要注册账号就可以下载)
下面书归正传
区域统计
很多时候我们得到的栅格数据需要经过一些统计分析才能得到有效的结论,尤其是地理数据,更需要结合shp,ArcMap的Zonal Statistic就是一个非常好的工具。
Zonal 统计的一些工具都位于工具箱 Spatial Analyst Tools下面,包括Tabulate、Fill、Histogram、Statistics等,这里主要讲Zonal Histogram和Zonal Statistics as Table,一个是直方图,一个是按表格统计。
直方图(Zonal Histogram)
这个并不是狭义上的画直方图,而是按一定的取值范围得到统计数据表格,输入包括Zonal data(就是shp),Zone field是指使用的shp属性表里的对应列,也就是按什么去统计(按省、按市还是按区的区别),Input Raster是指输入的栅格文件,output输入输出位置与表名。
计算完毕(右下角弹窗有显示是否成功)就可以在文件列表里看见这个表格了,可以在ArcMap查看,如果没有的话,可以看一下是不是文件Content Table选成别的了
按表格统计区域数据
输入和上面的一样,下面的Statistic type可以选择个别项,也可以选择All
需要注意的是,像Majority、Minority、Veriety、Median这些需要输入的Raster为整型数据,可以在Spatial Analyst Tools–Math–Int中进行转换。
Error处理
在一些比较罕见的情况下,可能出现报错(也就是右下角升上来的弹窗显示Failed),这个时候我们可以点进去看错误原因,比如下面这个错误类型为010423的,我在网上查到可能是输入数据不满足要求,比如存在空值之类的。像其他的error也有因为数据存储位置不正确(没有geodatabase、命名带有中文、命名以数字开头等等)。
但是,其实像geodatabase之类的问题,只要用catalog再建一个就好了,命名的问题也比较容易解决,而我当时遇到的问题,所有的方法都试过了,但是还是不行,尤其是同样的数据,一个可以一个不可以,这个时候可以不要纠结于为什么总是失败,因为有可能并不是你的问题,ArcMap有时候可能就是有bug,纠结了一个小时的问题,可能经过关机、重启等操作,突然就好了,所以有一样情况的朋友,可以不要过于纠结。
类似的事情我也在IDL上经历过,如果所有可能的问题都解决了还是不行,那么一定是ArcMap犯病了,让电脑休息一下,人也休息一下,过会儿再试试就好了,过会儿还不行就过一天再试(被折磨之后得出的结论,这是Esri玄学)