24小时热门版块排行榜    

CyRhmU.jpeg
查看: 5598  |  回复: 25

大傻meier

新虫 (正式写手)

[求助] 折磨我很久的 地统计学分析,望大神指教已有3人参与

在约52×100m的范围内按地形等差异布置了100个点,取0-10cm土样,做土壤CEC的空间变异性。
经纬度坐标在ARCGIS里投影成平面坐标,数据踢除异常值,做正态性转换。
在GS+里,参数Active lag distance(样点间最大距离的一半)和Lag class distance interval(约等于样点间最小距离)采用软件默认值,未做更改。得到的半变异函数模型为高斯模型,具体参数见表1,半变异函数图1。
在Gis里,做趋势分析,数据南北方向呈倒U型趋势(图2)。半变异云图,系统默认的为图3,资料里说“步长*步长组的值应该在最大的样点对距离的1/2左右,Lag size可由平均最近邻点距离求得(Spatial Statistics toolbox—Analyzing Patterns---Average Nearest Neighbor tool)”,故将Lag Size和number of lags更改为图4中的数值。
在地统计学向导里,插值方法选择Ordinary Kring—Prediction,移除2次趋势。更改Lag Size 和number of lags后,软件生成的半方差模型为图5,线性模型,和GS+里得到的结果很不一样。
我把GS+得到的模型参数输入,得到图6,但是Prediction error里标准均方根预测误差达到12了(图7),明显错了,

折磨我很久了,论文卡在这里,干着急上火,不知道可以问谁。没有学过GIS,不知道我这操作过程有没有什么错误。地统计学分析到底该怎么做嘞,求各位指教,我如果有说的不明白的,我再解释。谢谢,所有金币奉上

折磨我很久的 地统计学分析,望大神指教
表1.PNG


折磨我很久的 地统计学分析,望大神指教-1
图1.PNG


折磨我很久的 地统计学分析,望大神指教-2
图2.PNG


折磨我很久的 地统计学分析,望大神指教-3
图3.PNG


折磨我很久的 地统计学分析,望大神指教-4
图4


折磨我很久的 地统计学分析,望大神指教-5
图5.PNG


折磨我很久的 地统计学分析,望大神指教-6
图6.PNG


折磨我很久的 地统计学分析,望大神指教-7
图7.PNG
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

keguigiser

新虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
大傻meier: 金币+20, 有帮助, 谢谢 2014-11-09 13:21:13
不要修改Lag Size,默认的试一试。

[ 发自手机版 http://muchong.com/3g ]
肯定自己,超越自己,创造自己
2楼2014-10-25 13:14:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

keguigiser

新虫 (著名写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 大傻meier at 2014-10-25 16:50:58
默认的是这个,也是线性模型。
explore data里的半变异云图和Geostatistical Wizard里的这一步,怎么得到的semivariogram图和Lag size也不一样呢

0.PNG
...

是一样的啊,4.17*12就是整个半变异的横坐标啊

[ 发自手机版 http://muchong.com/3g ]
肯定自己,超越自己,创造自己
4楼2014-10-25 18:38:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

truebelief

木虫 (著名写手)

都是去除2次趋势再对残差建模的吗?lag size大了,导致没有一个明显上升的semivariogram。还有不要把sill和nugget都定死了,这样根本就没有办法fit semivariogram。

» 本帖已获得的红花(最新10朵)

5楼2014-10-26 09:19:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大傻meier

新虫 (正式写手)

引用回帖:
8楼: Originally posted by myjing at 2014-11-08 17:48:33
还有ArcGIS中怎么看变程呢?去趋势又是啥啊?我选的是simple kriging,和ordinary kriging具体有什么区别呢?也要选那个去趋势么?

Major range就是变程。如果有明显的趋势,还是要去的吧,插值后差别蛮大的。http://webhelp.esri.com/arcgisde ... tat/Geostat_3_2.htm
操作过程讲解。百度文库中也有很多中文教程
10楼2014-11-09 13:12:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大傻meier

新虫 (正式写手)

引用回帖:
14楼: Originally posted by luckyrocks at 2015-04-22 22:54:35
楼主,我也遇到了同样的问题,请问你解决了吗
我也是在做克里格插值的时候,用着两个软件,发现差别很大,我怀疑GS+和ARCGIS的单位有数量级的差别?
你做出来了吗

我算没做出来吧,中期的时候老师请了个数学模型方面的教授,噼里啪啦说了我一通,大概就是R2太低,RSS太大,数据不成正态性有问题,克里格插值后需要精度验证。帮不到你了
15楼2015-04-23 14:40:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

大傻meier

新虫 (正式写手)

引用回帖:
2楼: Originally posted by keguigiser at 2014-10-25 13:14:38
不要修改Lag Size,默认的试一试。

默认的是这个,也是线性模型。
explore data里的半变异云图和Geostatistical Wizard里的这一步,怎么得到的semivariogram图和Lag size也不一样呢
折磨我很久的 地统计学分析,望大神指教-8
0.PNG

3楼2014-10-25 16:50:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大傻meier

新虫 (正式写手)

送红花一朵
引用回帖:
5楼: Originally posted by truebelief at 2014-10-26 09:19:00
都是去除2次趋势再对残差建模的吗?lag size大了,导致没有一个明显上升的semivariogram。还有不要把sill和nugget都定死了,这样根本就没有办法fit semivariogram。

嗯,去除了2次趋势。Lag size怎么确定呢,软件默认的是4.175,用平均最近邻点距离求出来是2.089606. sill和nugget不要定死了什么意思?资料里说GS+半变异函数拟合厉害,GIS作图厉害,大家不都是GS+里计算出半变异模型参数,然后在GIS里插值作图的么?
6楼2014-10-29 19:03:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

myjing

新虫 (初入文坛)

请问楼主这个问题解决了么?我也是在做这方面的东西,这两个软件也是把我折磨坏了,为什么我用GS+做的给出的模拟最好的模型的r2竟然是0,简直抓狂啊
7楼2014-11-08 17:40:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

myjing

新虫 (初入文坛)

还有ArcGIS中怎么看变程呢?去趋势又是啥啊?我选的是simple kriging,和ordinary kriging具体有什么区别呢?也要选那个去趋势么?
8楼2014-11-08 17:48:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大傻meier

新虫 (正式写手)

引用回帖:
7楼: Originally posted by myjing at 2014-11-08 17:40:08
请问楼主这个问题解决了么?我也是在做这方面的东西,这两个软件也是把我折磨坏了,为什么我用GS+做的给出的模拟最好的模型的r2竟然是0,简直抓狂啊

没解决,不知道怎么搞,做实验有几天没学了
你在ARCGIS半变异云图里,观察半变异特别异常的点对,把它去掉,R2就上去了
自己琢磨出来的,不知道对不对
9楼2014-11-09 13:09:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 大傻meier 的主题更新
信息提示
请填处理意见