【附证明】用ArcGIS中Band Collection Statistics做相关性分析可能存在错误

您所在的位置:网站首页 相关系数可以计算出两个数吗对吗为什么 【附证明】用ArcGIS中Band Collection Statistics做相关性分析可能存在错误

【附证明】用ArcGIS中Band Collection Statistics做相关性分析可能存在错误

2024-07-12 17:30| 来源: 网络整理| 查看: 265

ArcGIS相关性分析

Spatial Analyst Tools——Multivariate(多元分析)——Band Collection Statistics(波段集统计)。

添加图层,勾选Compute covariance and correlation matrices以输出相关第分析结果,结果保存成txt。

使用的是皮尔逊相关系数(Pearson Correlation Coefficient)。

Spatial Analyst Tools——Multivariate——Band Collection Statistics 添加图层,勾选Compute covariance and correlation matrices以输出相关第分析结果,结果保存成txt。 使用的是皮尔逊相关系数(Pearson Correlation Coefficient)

数据的解释

arcgis10.8的结果

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-idBUNV5q-1661595291357)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220826223504777.png)]

arcgis pro的结果

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-0LbnxuJj-1661595291359)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220827152704314.png)]

从这里可以知道无论是ArcGIS还是pro都是一样的。

创建渔网,进行分析。

这里的渔网,随便找一个矢量确定其位置就好。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-jqMKjdRk-1661595291362)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220827153852394.png)]

新建一个字段。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-mw2j0mLk-1661595291362)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220827154008345.png)]

随机生成1-50的值

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-5g81MNEP-1661595291365)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220827154140803.png)]

import random def rannumb(): res=random.randint(0,50) return res

转成栅格

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-rYM1ZcaY-1661595291366)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220827154303666.png)]

用字段计算器生成两个栅格。但是这样生成的两个栅格有相同值,在属性表中的体现为:一个值的count并非为1.

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-d02f66LX-1661595291367)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220827163422612.png)]

这样也因此没有49个数,为了后期方便验证,我还是把每个值设定为不同的值。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-lTrJvSPr-1661595291367)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220827163520847.png)]

代码什么意思自己感悟,不懂的话可以在下面留言。

import random data = random.sample(range(0, 100), 100) rec=0 def Increment(): global rec global data Start = 1 #初始值 Interval = 1 #间隔值 if (rec == 0): rec = Start else: rec = rec + Interval data_rec=data[rec] return data_rec

直接计算。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-7oy9bOOt-1661595291368)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220827163829863.png)]

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-Klj5mIT7-1661595291368)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220827163844875.png)]

在属性表中统计的和用Band Collection Statistic计算的比较,最大最小值和平均值是一样的,但是标准差是不一样的。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-omxYHwgM-1661595291369)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220827164420823.png)]

再用excel验证一下。即在excel中计算一下标准差。这两个是不同的标准差。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-9rihGCqm-1661595291369)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220827170217692.png)]

再看一下协方差,这里的协方差是计算两个相同的数的。乍一看好像和上面Band Collection Statistic计算的没有什么关系。看来还是要看两个,不同变量的协方差来看关系。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-tIjPx32p-1661595291371)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220827171838004.png)]

之前我忘记弄记录两个栅格的位置对应关系了。所以我根据上面的方法新建两个图层。

位置对应关系。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-GiwjaFFi-1661595291371)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220827172403337.png)]

新建的long3和long4的关系。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-5lBMhdK8-1661595291372)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220827172535515.png)]

于是导出两个表格计算协方差。好像也是不对的,这里说的协方差那么我们就不能确定了。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-f11NDxUJ-1661595291373)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220827172809724.png)]

再算算相关性吧,无容置疑,-1.01315的相关性本身就超出了相关性的值域。这是在excel里面算出的相关性。这个更有可能。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-n7cysfD2-1661595291373)(C:/Users/Administrator/AppData/Roaming/Typora/typora-user-images/image-20220827173037860.png)]

总结:ArcGIS来算相关性是不合适的。

参考文献

论ArcGIS10.2的Band Collection Statistics工具计算相关系数的正确性

支撑数据 数据下载



【本文地址】


今日新闻


推荐新闻


CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3