R程序改写为matlab程序

ywh19860616 发表于 2012-5-16 20:16
非常感谢epoh老师
哈哈，epoh老师，对于标准化，我的看法还是和您的不同，请不要介意啊。
您看我上传的 ...

呵呵,讨论原本就应有不同想法,
假设依你的想法,标准化的话,
可以这样处理
[Xij_bar,GGI,ZZG,ZZ]=Getis_revised(row_norm(WWij),Xij(p,:));
但是抓不到k
k一直往上爬,到1000,3000都是一样的值,无法辨识
你只要作Spatial Filtering 就知道答案了

epoh 发表于 2012-5-16 20:01
底下这两种结果是相同的
针对75楼copy错误修改过
仔细推敲文献,应该是标准化再代入exp()没错

epoh老师，刚用您79楼的程序试算了一下，结果还是不理想。
其实我有点好奇，为什么计算出来的Xij_bar会和Xij差那么大呢，
而且Xij_bar都比Xij大，就是说过滤后的值比未过滤的更大。
从理论来说，没有规定Xij_bar不能比Xij大，但是从别人文献结果来看，
L=Xij-Xij_bar基本都是有正也有负的。
不明白是不是真是我的数据问题。

ywh19860616 发表于 2012-5-16 21:23
epoh老师，刚用您79楼的程序试算了一下，结果还是不理想。
其实我有点好奇，为什么计算出来的Xij_bar会和 ...

可能是数据转换问题
刚又看了下
这篇文献也是要cartesian coordinates
Spatial Filtering with EViews and MATLAB
2.1 Global Spatial Structures
(theta,w) are the cartesian coordinates

epoh 发表于 2012-5-16 21:42
可能是数据转换问题
刚又看了下
这篇文献也是要cartesian coordinates

是的啊，我看了那么久，怎么都没有注意到这个，非常感谢epoh老师的提醒。
现在的程序调用的是distance_wm,那是根据经纬度计算的 the Great Circle distance，
而不是在cartesian coordinates下的坐标系统。

ywh19860616 发表于 2012-5-16 22:06
是的啊，我看了那么久，怎么都没有注意到这个，非常感谢epoh老师的提醒。
现在的程序调用的是distance_w ...

我记得这个问题答复过你
但是我忘记转换函数
你还记得吗?

epoh 发表于 2012-5-16 22:17
我记得这个问题答复过你
但是我忘记转换函数
你还记得吗?

epoh老师，您是指从经纬度到cartesian coordinates的转换函数吗？这个好像没有看见过，我不知道。
您还记得大致是什么时候告诉我的吗？

ywh19860616 发表于 2012-5-16 22:25
epoh老师，您是指从经纬度到cartesian coordinates的转换函数吗？这个好像没有看见过，我不知道。
您还记 ...

NDSF Utility: LatLong to XY
http://www.whoi.edu/marine/ndsf/ ... om=LatLon&to=XY

epoh 发表于 2012-5-17 07:53
NDSF Utility: LatLong to XY
http://www.whoi.edu/marine/ndsf/cgi-bin/NDSFutility.cgi?form=0&from=L ...

epoh老师，太感谢您了，还有这个功能
不过我试了下，怎么都转换不成功哦。比如：

35° 18′ 34″ N, 106° 33′ 7″ W

35.309444, -106.551944

转换为X和Y，得出的结果是NAN,NAN。
不知道这一栏要求填什么内容？Enter the Origin of X/Y coordinate:

ywh19860616 发表于 2012-5-17 10:08
epoh老师，太感谢您了，还有这个功能
不过我试了下，怎么都转换不成功哦。比如：

你先填0,0试试
另请再看看这个function如何
wgs2utm.rar (1.67 KB) 本附件包括：

• wgs2utm.m

2012-5-17 10:17:10 上传

epoh 发表于 2012-5-17 10:16
你先填0,0试试
另请再看看这个function如何

好的，谢谢您
那个Matlab函数可以运行，结果：
Lat=35.309444;Lon=-106.551944;
[x,y,utmzone] = wgs2utm(Lat,Lon)
x =
3.5891e+005
y =

  3.9085e+006
utmzone =

    13从google资料来看，这里应该为6的。


另外，如果那个网页形式的第二个框都是输入0,0，结果和这个不一样。

我再学习下其他地理软件，确认下结果。

google参考：
http://www.google.com.hk/#hl=zh-CN&newwindow=1&safe=strict&site=&source=hp&q=Bernalillo&btnK=Google+%E6%90%9C%E7%B4%A2&oq=&aq=&aqi=&aql=&gs_l=&bav=on.2,or.r_gc.r_pw.,cf.osb&fp=49354fa1230e3dc5&biw=1440&bih=693