采用邻域分析中的邻域丰度 (Neighborhood Enrichment)
作为空间自相关因子 Autocovi,反映土地利用类型转化的邻域关系,将邻域影响纳入地类变化的驱动因子
中,进行二元Logistic回归分析,这样,邻域丰度的引入能够弥补单纯Logistic回归分析
的不足。表达式:http://control.blog.sina.com.cn/admin/article/article_add.php?is_new_editor=1
F(i,k,d)=[n(k,d,i)/n(d,i)]/[N(k)/N]
式中:F(i,k,d)为邻域丰度因子,其中i代表栅格位置,k为土地利用类型,d为邻域半径;n(k,d,i)
为i栅格d半径范围内k土地利用类型的栅格个数;n(d,i)为i栅格d半径内栅格总数量;N(k)为
整个库区内k土地利用类型的栅格总个数;N为整个库区内的总栅格数。
本文在ArcGIS 9.3的Neighborhood
Statistics工具中设置邻域统计大小为10 m×10 m,即邻域距离为10
m。
由于好久没有用9.3的版本,现在用10.3的版本,就以10.3为例展开。
领域分析工具中,有计算重叠邻域和非重叠邻域的统计量。计算非重叠领域的工具,是焦点运算,另一种块统计针对位置不重叠领域。在焦点运算中,只有当前正在处理的像元才会接收邻域中所有像元的计算结果,而在块运算中,处于邻域的最小外接矩形内的所有像元都会接收相同的输出值。
因而利用,焦点统计工具
(FocalStatistics),计算每个地类栅格的总数,输出OutRas文件,获得n(k,d,i),之后利用栅格计算器,除以搜索半径内的栅格总数n(d,i)以及N(k)/N。最后,获得k地类的领域丰富。
欢迎大家指教交流。
参考:http://resources.arcgis.com/zh-cn/help/main/10.1/index.html#/na/009z000000qs000000/