ltrend命令：求线性趋势

半头砖

技能get 谢谢分享

jolincai

求帮忙啊  按着一步一步来 怎么图是这样的...？

jolincai

接楼上！ctl和gs如下：
dset E:\data\NCEP-1\1\35nc1tmin2.5.grd
undef -9.99e+33
title 1979-2013tmin
xdef 29 linear 70 2.5
ydef 17 linear 15 2.5
tdef 35 linear 00Z01jan1979 1yr
zdef 1 linear 1 1
vars 1
tmin 1 0  ** surface Min 2m temp since previous post-processing K
ENDVARS



'reinit'
'open E:\data\NCEP-1\1\35nc1tmin2.5.ctl'
'set fwrite E:\data\NCEP-1\1\liner.grd'  
'set gxout fwrite'   
'set lon 73 135'         
'set lat 17 55'
'set t 1 35'
'ltrend tmin slope'
'd slope'         
'disable fwrite'   
;
至此 是生成文件liner.grd


然后配上ctl和gs画图:
dset E:\data\NCEP-1\1\liner.grd
undef -9.99e+33
title 1979-2013tmin
xdef 29 linear 70 2.5
ydef 17 linear 15 2.5
tdef 1 linear 00Z01jan1979 1yr
zdef 1 linear 1 1
vars 1
lin 1 0  ** surface Min 2m temp since previous post-processing K liner
ENDVARS



'reinit'
'open E:\data\NCEP-1\1\liner.ctl'
'set grads off'
'set grid off'
'set mpdset cnriver'
'set map 1 1 5'
'set xlopts 1 5 0.21'
'set ylopts 1 5 0.21'
'set lat 17 55'
'set lon 73 135'
*'set strsiz 0.3'
*'set string 1 c 6'
*'draw string 3 8.1 (a) GRADS'
*设置等值线标记的颜色
'set clopts 1 5 0.18'
*设置样条插值光滑开关
'set cterp on'
*设置是否将网格值重新插值              
'set csmooth on'
'set clskip 2'
'set gxout shaded'
*'colormap blueyellowred'
'cnbasemap smth9(lin)'
'cbarn 1 0'
*等值线用黑色
'set gxout contour'
'set ccolor 1'
'cnbasemap smth9(lin)'
*'cnbasemap std'
*'set mpdset hires'
'gxprint E:\data\NCEP-1\1\gradslin.pdf white'
'gxprint E:\data\NCEP-1\1\gradslin.eps white'
'gxprint E:\data\NCEP-1\1\gradslin.png white'
;

jolincai

接楼上，我的本意是想画如文献里这样的温度趋势图：

mangguojun

谢谢楼主

tthappy

长见识了，实践中，希望有用，谢谢楼主！

dengshulin12531

楼主，请问Itrend在OpenGrDAS中怎么使用啊，非常感谢！

泡zombie

file:///d:/gyw/lhtfl.png请问楼主大神，为什么我画出来的图是这样的？

yangping

jolincai 发表于 2016-3-9 16:37
接楼上，我的本意是想画如文献里这样的温度趋势图：

请问同学，这类图你有画出来吗？可以请教这是怎么画的吗？急求啊感谢

常山赵子龙

请问楼主，matlab当中是如何计算的线性趋势呢？