24小时热门版块排行榜    

Znn3bq.jpeg
查看: 5882  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085400 328分 求调剂 +5 喂你一个大橙子 2026-04-09 5/250 2026-04-09 11:34 by yongzhesheng
[考研] 285求调剂 +9 AZMK 2026-04-07 10/500 2026-04-09 10:29 by 小物理化学
[考研] 考研调剂 +12 冰冰,,, 2026-04-07 12/600 2026-04-09 09:04 by 5268321
[考研] 349学科化学045106求调剂,化学类都可以 +7 保好懂懂 2026-04-08 7/350 2026-04-09 07:37 by jp9609
[考研] 318求调剂 +13 ykyhsa 2026-04-05 15/750 2026-04-08 21:37 by wj165256
[考研] 265求调剂 +3 风说她早忘了 2026-04-07 3/150 2026-04-08 20:20 by 逆水乘风
[考研] 293分求调剂,外语为俄语 +7 加一一九 2026-04-07 10/500 2026-04-08 20:14 by yutian743
[考研] 一志愿生物与医药,296分,求调剂 +14 66鹿 2026-04-03 16/800 2026-04-08 10:38 by tjzhao
[考研] 本科郑州大学,一志愿华东师范大学282求调剂 +21 熊哥xtk 2026-04-07 24/1200 2026-04-08 09:48 by 123考研123
[考研] 326分,一志愿沪9,求生物学调剂 +4 刘墨墨 2026-04-05 4/200 2026-04-08 06:22 by lijunpoly
[考研] 材料调剂 +13 汉123456 2026-04-07 14/700 2026-04-07 22:53 by 来看流星雨10
[考研] 277求调剂 数一104分 +9 瓶子PZ 2026-04-05 14/700 2026-04-07 17:52 by 蓝云思雨
[考研] 277工科求调剂 +10 1915668 2026-04-04 11/550 2026-04-06 23:53 by 心媛意北
[考研] 0857大类环境工程B区求调剂 +3 龚禹铭 2026-04-05 3/150 2026-04-06 10:22 by 蓝云思雨
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 22:23 by hemengdong
[考研] 考研调剂生寻找导师 +3 顾瞻考研啊 2026-04-05 3/150 2026-04-05 18:18 by 啵啵啵0119
[考研] [调剂信息]085408光电信息 求调剂 总分291分数一英一 +3 iz11az 2026-04-02 3/150 2026-04-04 19:09 by 蓝云思雨
[考研] 321求调剂 +13 认真求上学 2026-04-02 13/650 2026-04-04 18:23 by macy2011
[考研] 数二英二348求调剂 +4 hxdzj1 2026-04-03 5/250 2026-04-03 21:25 by zhq0425
[考研] 一志愿大工学硕,求调剂 +4 yub0811 2026-04-02 4/200 2026-04-02 21:36 by 百灵童888
信息提示
请填处理意见