用ltrend 命令，快速求出数据的趋势线！或者给数据去掉倾向！

zl.87625 · 发表于 2013-7-8 08:38:34

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 zl.87625 于 2013-7-9 20:47 编辑

平时经常画出一个时间序列以后，还需要画出它的趋势线（倾向线），看到很多同学都用程序算，

其实grads里面就有自带的ltrend 命令，可以快速给求出数据的趋势线！只要设定好时间起始。

比如，我的变量名称为海温异常 ssta，设置语句 'ltrend ssta bb'  ，这里的bb就是变量ssta 的趋势。

直接"d  bb" ，就可以画出源数据的趋势（一条斜的直线）。

如果需要画去过倾向的数据，就用原变量减去趋势线bb, 就是去过倾向的数据了，'define aa=ssta-bb'

'd  aa' 这里的变量 aa  就是去除倾向（趋势>数据，很简单，很方便！

_____________________________________________________________________

(我下面的程序是给一个三维数据去的倾向，如果是一串时间序列的数据，只需要 'set x 1' ，'set  y  1' 就可以了！

这个命令，既可以画图时候直接用，也可以处理数据时用。

希望这个方法，对大家有所帮助！第一次发帖子，感谢支持！

（ps：使用这个命令，需要将一个 ltrend.gs的文件，放在grads按装包里面的 lib文件夹下面！

我把这个gs放在附件里面大家可以下载）

'reinit'
'open  f:\sst\ersst\50-09anomaly.ctl'
'set  fwrite  f:\sst\ersst\1947-2009ssta-detrend.grd'
'set gxout fwrite'

'set lon    30 120'
'set lat    -60    0'
'set time Jan1947 Dec2009'
'ltrend  ssta  bb'
'define aa=ssta-bb'
'd  aa'

'reinit'

liuzf406 · 发表于 2014-10-26 16:17:46

有知道如何检验趋势显著性的吗？

zl.87625 · 发表于 2013-10-25 16:50:51

本帖最后由 zl.87625 于 2013-10-25 16:56 编辑

yaoyao 发表于 2013-10-22 14:21

登录/注册后可看大图

按照楼主的说法bb算出来应该是个二维场，为什么我算出来的bb是个三维的= =而且和原始场很像，一点都不像趋势 ...

ltrend 命令是二维场和三维场都适用的！
如果你的原始数据是个二维的，比如时间序列一串数，那么算出来的bb 就是一个直线（其实也是二维）。
如果你的原始数据是三维的，那算出来的bb确实也是三维的，因为每个格点上的数据在这段时间都有一个趋势，而你有很多个格点。所以bb的分布就是一个空间场。
比如我研究印度洋60年海温的年际变化特征，就可以先用ltrend 命令去掉所有格点上60年海温的趋势（算是年代际变化吧），然后再看海温的年际变化，这个时候反映的信息才是真实的。

kongfeng0824 · 发表于 2013-7-8 08:48:19

留言需要5贡献，贵点

zl.87625 · 发表于 2013-7-8 08:50:14

kongfeng0824 发表于 2013-7-8 08:48

登录/注册后可看大图

留言需要5贡献，贵点

啊不好意思！我立马修改一下！新手报到，现在还不太懂这些贡献啊什么的具体是什么意思，sorry

兰溪之水 · 发表于 2013-7-8 12:34:38

使用这个还需要qdims.gsf这个文件，可以到这里下：
http://www.atmos.umd.edu/~bguan/grads/GrADS_Scripts.html

zl.87625 · 发表于 2013-7-8 12:49:33

兰溪之水发表于 2013-7-8 12:34

登录/注册后可看大图

使用这个还需要qdims.gsf这个文件，可以到这里下：
http://www.atmos.umd.edu/~bguan/grads/GrADS_Scripts ...

去了你那个网址，东西好全啊，好实用，谢谢你啊！另外，想请教一下，一个时间序列，我画了趋势线之后，想知道这个趋势系数是多少，有木有鉴别的方法吗？

看一般趋势线都会配合给出一个 y=kx+b 这样的方程

兰溪之水 · 发表于 2013-7-8 12:55:40

zl.87625 发表于 2013-7-8 12:49

登录/注册后可看大图

去了你那个网址，东西好全啊，好实用，谢谢你啊！另外，想请教一下，一个时间序列，我画了趋势线之后，想 ...

你看看那函数里面的计算就知道了~可以找到相关的系数，然后输出就行了~

zl.87625 · 发表于 2013-7-8 12:58:28

恩好的，明白了，这就去看看，谢谢老师指导！

zl.87625 · 发表于 2013-7-8 14:23:05

zl.87625 发表于 2013-7-8 12:58

登录/注册后可看大图

恩好的，明白了，这就去看看，谢谢老师指导！

您好，我仔细的看了那个ltrend.gs 大概明白了，斜率就是slope. 还有均方根方差 rmse ，

但是我读入ctl的数据，运行了ltrend.gs 文件，都木有问题，只是 d slope 或者 d rmse 的时候，就说不认识这个变量名，这是可能什么原因？

风儿飘飘 · 发表于 2013-7-11 14:01:00

非常实用，谢过！

guobch · 发表于 2013-7-12 08:06:44

懂了谢谢~~

		自动登录	找回密码
密码			立即注册

[秀图] 用ltrend 命令，快速求出数据的趋势线！或者给数据去掉倾向！

登录后查看更多精彩内容~

评分

回帖奖励 +1 金钱