24小时热门版块排行榜    

Znn3bq.jpeg
查看: 5902  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291求调剂 +9 关忆北. 2026-04-14 9/450 2026-04-16 13:49 by asdsyh1111
[考研] 0854调剂 +9 950824he@ 2026-04-09 9/450 2026-04-16 11:11 by noqvsozv
[考研] 344求调剂 +8 丶风雪夜归人丶 2026-04-09 8/400 2026-04-16 10:49 by 圆心602
[考研] 294求调剂 +11 淡然654321 2026-04-15 11/550 2026-04-16 10:23 by beilsong20
[考研] 一志愿华中农业071010,320求调剂 +12 困困困困坤坤 2026-04-14 14/700 2026-04-16 10:18 by beilsong20
[考研] 恳请有学校收留 +6 柯淮然 2026-04-12 6/300 2026-04-16 07:23 by 学员JpLReM
[考研] 求调剂推荐 +8 小聂爱学习 2026-04-14 8/400 2026-04-16 07:22 by 学员JpLReM
[考研] 求调剂学校 +14 不会吃肉 2026-04-13 16/800 2026-04-15 21:59 by noqvsozv
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
[考研] 食品与营养(0955)271求调剂 +15 升格阿达 2026-04-12 16/800 2026-04-14 13:18 by 浮若_安生
[考研] 305求调剂 +8 玛卡巴卡boom 2026-04-11 8/400 2026-04-14 09:04 by pengliang8036
[考研] 300分求调剂 (085501机械专硕,本科扬大) +9 xu@841019 2026-04-11 10/500 2026-04-14 08:48 by 木木mumu~
[考研] 材料复试求调剂 +24 xhhdjdjsjks 2026-04-09 24/1200 2026-04-13 15:49 by 幸免 ..
[考研] 297工科,求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[教师之家] 山东双非院校考核超级无底线,领导幸灾乐祸,教师遭殃恐 +3 qut2026 2026-04-11 7/350 2026-04-12 20:24 by qut2026
[考研] 295分求调剂 +13 ?要上岸? 2026-04-10 13/650 2026-04-12 15:37 by laoshidan
[考研] 291求调剂 +8 关忆北. 2026-04-11 8/400 2026-04-12 09:32 by 逆水乘风
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
[考研] 广东省 085601 329分求调剂 +14 Eddieddd 2026-04-10 14/700 2026-04-11 09:58 by bljnqdcc
[考研] 本科西工大 0856 324求调剂 +10 wysyjs25 2026-04-09 11/550 2026-04-10 08:37 by 5268321
信息提示
请填处理意见