楼主: zhaowill
18355 106

GeoDA软件的使用 [推广有奖]

21
zhaowill 发表于 2010-3-6 17:33:48 |只看作者 |坛友微信交流群
空间计量分析很多是用Matlab编程来做。
局部G统计的Matlab计算程序如下:
function [GI,ZG]=Getis(Wij,Xij)
% 计算G统计值:G值是不包含自己的G统计值
% 输入:Wij—空间权值矩阵
% 输入:Xij—研究区域的空间属性数据
% 输出:GI— 空间局部自相关指标-G统计值
% 输出:ZG— 对于GIJ的检验Z值的计算结果
% GI=Getis(Wij,Xij,N)

[M,N]=size(Wij);
GI=zeros(N,1);
GIn=0;
Xn=0;
for m=1:N
    for n=1:N
        if n==m
            GIn=GIn;
            Xn=Xn;
        else
            GIn=GIn+Wij(m,n)*Xij(n);
            Xn=Xn+Xij(n);
        end
    end
    GI(m)=GIn/Xn;
    GIn=0;
    Xn=0;
end
save GI GI -ascii;


%计算Z值
EG=zeros(N,1);
VarG=zeros(N,1);
ZG=zeros(N,1);
WI=0;
for m=1:N                                                        %计算E(Gi(d))
    for n=1:N
        if n==m
            WI=WI;
        else
            WI=WI+Wij(m,n);
        end
    end
    EG(m)=WI/(N-1);
    WI=0;
end
WI=0;
Xj=0;
Xj2=0;
%Yi1=zero(N,1);                                                    %计算Yi1
%Yi2=zero(N,1);                                                    %计算Yi2
Wi=zeros(N,1);                                                     %计算Wi
for m=1:N
    for n=1:N
        if n==m
            WI=WI;
            Xj=Xj;
            Xj2=Xj2;
        else
            WI=WI+Wij(m,n);
            Xj=Xj+(Xij(n)*Xij(n));
            Xj2=Xj2+((Xij(n)*Xij(n))/(N-1));
        end
    end
    Wi(m)=WI;
    Yi1=Xj;
    Yi2=Xj2+Yi1*Yi1;
    WI=0;
    Xj=0;
    Xj2=0;
end
for m=1:N                                                           
    VarG(m)=(Wi(m)*(N-1-Wi(m))/((N*N-1)*(N*N-1)*(N-2)))*(Yi2/(Yi1*Yi1));                %计算Var(Gi(d))
    ZG(m)=(GI(m)-EG(m))/(sqrt(VarG(m)));                                                %计算Z(Gi)
end

使用道具

22
zhangtao 发表于 2010-3-6 17:42:03 |只看作者 |坛友微信交流群
以上程序运行后怎么提示:Error: Missing variable or function.那位朋友能告诉为什么?

使用道具

23
zhaowill 发表于 2010-3-6 18:15:39 |只看作者 |坛友微信交流群
对了,geoda可以做局部空间自相关。其操作步骤如下:
1、载入包含shp数据文件;
2、需要一个空间权重文件。如果你还没有制作权重文件,在开始分析前创建一个 rook 权重文件(比如取名为 xxx.GAL)。
3、从菜单选择 Space>Univariate LISA,或点击相应的工具按钮,开始 Local Moran 功能。将会出现一个熟悉的变量选择对话框。选择要分析的变量,然后OK。
4、接着出现一个权重选择对话框。指定权重文件xxx.GAL,继续点击OK。最后一步,出现结果选项对话框,可以全选或部分选择,ok即大功告成了。

使用道具

24
sangfriend 发表于 2010-3-8 10:15:20 |只看作者 |坛友微信交流群
非常感谢!

使用道具

25
buhanying 发表于 2010-3-21 13:15:42 |只看作者 |坛友微信交流群
谢谢哈

使用道具

26
murong2009 发表于 2010-4-1 17:40:39 |只看作者 |坛友微信交流群
楼主是个好人

使用道具

27
阳光小子 发表于 2010-4-2 11:41:22 |只看作者 |坛友微信交流群
谢谢楼主!十分感谢!

使用道具

28
yearend 发表于 2010-4-3 08:32:02 |只看作者 |坛友微信交流群
good PDF for GeoDA
What is the name of this book??

使用道具

29
kevinlxl2006 发表于 2010-4-5 21:20:32 |只看作者 |坛友微信交流群
[url=http://www.pinggthankyou,好好学习以下
本文来自: 人大经济论坛 详细出处参考:http://www.pinggu.org/bbs/viewthread.php?tid=580980&page=1u.org/bbs/redirect.php?goto=findpost&pid=3652931&ptid=580980]1#[/url] zhaowill

使用道具

30
826loving 发表于 2010-4-17 23:54:51 |只看作者 |坛友微信交流群
谢谢啦!楼主好人啊!

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-6-17 07:46