24小时热门版块排行榜    

查看: 5826  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +3 sbdksD 2026-03-19 3/150 2026-03-19 23:21 by fmesaito
[考研] 288求调剂 +15 于海海海海 2026-03-19 15/750 2026-03-19 22:41 by 学员8dgXkO
[考研] 296求调剂 +3 www_q 2026-03-18 6/300 2026-03-19 22:28 by zhq0425
[考研] 能源材料化学课题组招收硕士研究生8-10名 +4 脱颖而出 2026-03-16 12/600 2026-03-19 16:17 by 脱颖而出
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 10/500 2026-03-19 16:05 by 余麟余
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 工科材料085601 279求调剂 +6 困于星晨 2026-03-17 6/300 2026-03-18 10:21 by kkcoco25
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
信息提示
请填处理意见