24小时热门版块排行榜    

查看: 5374  |  回复: 31

dream_yaoyao

银虫 (小有名气)

[交流] 交叉学科的人,在arcgis里做空间插值,竟然用经纬度坐标,我该怎么告诉他们(哭)已有12人参与

很多交叉学科的人,觉得用arcgis做空间插值,他们的数据一下子从几个点上的数据变成了面上的数据,而且还有空间规律,那个万般欣喜啊,于是到处滥用空间插值。
其实做空间插值对数据有很多限制条件的,尤其是克里格插值,对数据的要求很严格,可是在arcgis中做插值,不管是什么样子的数据,都可以有一个所谓的“结果”,也就是arcgis里做空间插值并不提供数据阀门,但是并不代表什么数据都可以做插值,就像你做传统的统计分析,不是还要筛选数据看有没有异常值,可是没学过地理信息的人,就随便乱用,应该说是滥用,让我惊讶最多的是很多人在经纬度坐标下做插值。数据插值的算法是以长度单位计算的,怎么能用经纬度呢,可是真的有好多人在用,包括一些非本学科的老师,弄得我都糊了,难道是我错了?请大神们坚定的告诉我,arcgis里的空间插值不能在经纬度下计算!这样是不对的!我需要坚定我的三观!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

Zoning

专家顾问 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
9楼: Originally posted by dream_yaoyao at 2014-05-03 17:01:56
高手,感觉你是理解我的这个问题的,我之前也是这么想的,可是我却相反了,我想成大区域误差会更大,我就想,用经纬度肯定是会有误差的,这个误差是必然的,那么区域越大误差积累越大,现在想想应该是你这么解释的 ...

首先,对楼主这种格物致知的科研精神很有感触;没像其他层主那样注意到您的性别,根据“回复是否详细判断性别”,这可能会是一个好的idea,但是我注定是个反例,因为我回答也够啰嗦的,且我是个直男。
其次,楼主MM的这个命题,本人也深以为然。最近逛小木虫,发现GIS越来越受交叉亦或非相关学科的关注了,起码在制图方面。武大的秦昆教授说过“空间分析是GIS的灵魂”,我想他应该所指的便是地统计分析了。
最后,回到正题,我本人对ArcMap里面的插值也有些嗤之以鼻,这种放任的态度,很容易让部分非GIS专业误入歧途。没有显著性检验的地统计分析就像是脱缰的野马。
PS:最后一点,“R^2”是拟合优度,非显著检验。
盗用星爷《功夫》台词:小鬼,老衲见你根骨奇好,以后共同切磋武艺如何……呵呵 祝好!
我喜欢黑,因为它可以把我隐藏起来……
17楼2014-05-04 11:30:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangok

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
通常情况下确实是应该将经纬度坐标转换为UTM投影坐标,但是对于大区域的数据也可以用经纬度坐标,SAM这个很不错的空间统计分析软件就可以实现,当然最终的空间统计分析结果是需要经过显著性检验的。其实,对于ARCGIS下的空间插值,如果象你所说的滥用得到的结果,看似很漂亮,但是经检验无效,还是空欢喜一场,所以,你可以建议那些人补充一下显著性检验。
2楼2014-05-03 00:39:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dream_yaoyao

银虫 (小有名气)

引用回帖:
4楼: Originally posted by dibazhang at 2014-05-03 16:12:06
人家只要看看趋势分布,有些sci里的图都是用经纬度插出来的

是的,的确,我身边很多人这么做,但是他们都是非本学科的,用经纬先不管对不对,但是比如他们常用克里格插值,可是根本不检验半变异函数,就直接用,也发表了,只不过他们没放那些东西,有的是放上了半变异函数的那些参数,那些审稿人很多也都是交叉学科的,根本不懂,半变异函数的需要达到某种要求才可以往下用的,我上一篇sci,心想自己学这个的,要严格要求自己,就把半变异函数列上了,其中有的拟合程度根本没达到插值的要求,但是我看那些人也都发表就这么投过去了,结果遇到一个搞地理的大牛,把我的那块插值的东西全都指出来不行不行,半变异的函数的拟合不好是不能用的,最后整个一块内容让我换成IDW插值了,这个对数据的要求没那么苛刻。
所以也不能说发表出来的就是对的。
7楼2014-05-03 16:47:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangok

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
6楼: Originally posted by dream_yaoyao at 2014-05-03 16:38:27
可是那个是属性,无论哪个插值的算法,不都是以长度距离单位计算的么,比如克里格插值中的最关键的半变异函数,range之类的不都是长度单位么?比如IDW插值,反距等权,所谓距离不也是以长度单位计算的?难道arcgis ...

经纬度坐标与公里制的单位相差很大,1个纬度大致相当于110km的距离吧,对于小区域采用经纬度坐标肯定是点与点之间因距离太小而导致重叠,根本无法参与运算。所以对于诸如中国全境或者北美这样的大区域,各样点距离相差很大就可以抵消这样的影响,SAM软件中自带的计算鸟类物种多样性的数据文件我打开看过,各样点距离很大,可以这样直接用经纬度坐标参与运算。关于显著性检验,诸如半变异函数的显著性检验,其实就是标准化之后的z检验,在各类地统计学专业软件中都是自带的,具体可参阅软件使用说明。
8楼2014-05-03 16:48:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dibazhang

铜虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
人家只要看看趋势分布,有些sci里的图都是用经纬度插出来的
4楼2014-05-03 16:12:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dream_yaoyao

银虫 (小有名气)

引用回帖:
8楼: Originally posted by zhangok at 2014-05-03 16:48:31
经纬度坐标与公里制的单位相差很大,1个纬度大致相当于110km的距离吧,对于小区域采用经纬度坐标肯定是点与点之间因距离太小而导致重叠,根本无法参与运算。所以对于诸如中国全境或者北美这样的大区域,各样点距离 ...

高手,感觉你是理解我的这个问题的,我之前也是这么想的,可是我却相反了,我想成大区域误差会更大,我就想,用经纬度肯定是会有误差的,这个误差是必然的,那么区域越大误差积累越大,现在想想应该是你这么解释的。
也就是说即使使用的经纬度,在算法运行中也是把它按照类似于1个纬度相当于110km换成成距离计算的是吧?经纬度是没法直接参与计算的是吧,想想即使一个县这个尺度,经纬度也相差很小,所以应该是换算成长度单位计算的是吧?

另外可不可以请教你具体怎么检验显著性,arcgis不用提,只能留一部分点当做验证点去看验证结果。我用GS+,比如半变异函数也只有一个R^2,貌似不提供显著性检验。
9楼2014-05-03 17:01:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

echoone168

禁虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖内容被屏蔽

10楼2014-05-03 17:24:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qt211

至尊木虫 (文坛精英)

我是一个跑龙套的


小木虫: 金币+0.5, 给个红包,谢谢回帖
既然有时间,就弄个图对比呗
同一组数据,用不同坐标方法,结果一比就明了了。
大范围和小范围都做一下试试。
不同插值方法也可以对比。
☆敬天护地爱人☆
23楼2014-05-05 11:22:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skywonder

新虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
9楼: Originally posted by dream_yaoyao at 2014-05-03 17:01:56
高手,感觉你是理解我的这个问题的,我之前也是这么想的,可是我却相反了,我想成大区域误差会更大,我就想,用经纬度肯定是会有误差的,这个误差是必然的,那么区域越大误差积累越大,现在想想应该是你这么解释的 ...

用r做kriging比较专业

[ 发自小木虫客户端 ]
30楼2014-05-08 23:13:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

sapphirechina

铜虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
插值的本质,是根据一些实测点的属性,来推测其他空间位置的属性. 经纬度表示的是地理坐标.所有的投影坐标都是从地理坐标转过来的.既然用投影坐标可以插值,为什么不能用经纬度?
3楼2014-05-03 08:07:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dream_yaoyao

银虫 (小有名气)

引用回帖:
3楼: Originally posted by sapphirechina at 2014-05-03 08:07:45
插值的本质,是根据一些实测点的属性,来推测其他空间位置的属性. 经纬度表示的是地理坐标.所有的投影坐标都是从地理坐标转过来的.既然用投影坐标可以插值,为什么不能用经纬度?

可是那个是属性,无论哪个插值的算法,不都是以长度距离单位计算的么,比如克里格插值中的最关键的半变异函数,range之类的不都是长度单位么?比如IDW插值,反距等权,所谓距离不也是以长度单位计算的?难道arcgis中插值可以自动转换?可是根本就没有投影呀,比如GS+里面就只能用投影坐标,有点不太懂arcgis里面到底是怎么个算法。
5楼2014-05-03 16:35:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dream_yaoyao

银虫 (小有名气)

引用回帖:
2楼: Originally posted by zhangok at 2014-05-03 00:39:34
通常情况下确实是应该将经纬度坐标转换为UTM投影坐标,但是对于大区域的数据也可以用经纬度坐标,SAM这个很不错的空间统计分析软件就可以实现,当然最终的空间统计分析结果是需要经过显著性检验的。其实,对于ARCGI ...

可是那个是属性,无论哪个插值的算法,不都是以长度距离单位计算的么,比如克里格插值中的最关键的半变异函数,range之类的不都是长度单位么?比如IDW插值,反距等权,所谓距离不也是以长度单位计算的?难道arcgis中插值可以自动转换?可是根本就没有投影呀,比如GS+里面就只能用投影坐标,有点不太懂arcgis里面到底是怎么个算法。

此外,你说,对于大区域的数据可以用经纬度坐标,为什么呢?求原理
此外,怎么进行显著性检验,你指的半变异函数的显著性检验?还是什么的?
6楼2014-05-03 16:38:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dream_yaoyao 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见