时间序列图求助

gynxz

有四份资料，分别是2003/09/12-2003/10/12；2005/07/03-2005/08/02；2006/08/23-2006/09/22；2010/10/23/2010/11/22
数据时间不连续，且不同月份，都是一个月的长度。数据格式是 nc的，数据是风速资料。
现在想用grads画这四个数据的平均风速的时间序列，忙了两天了也没出来不知道怎么处理。

出图的效果想是这样的：经度固定，横坐标为时间从1-30天，纵坐标是纬度
不知道直接用grads怎么实现，哪位画过类似的图，帮帮我。

'sdfopen e:\data\2003u.nc'
'set lon 130'
'set lat 0 50'
'set t 1 31'
'define a1=uwnd'
'close 1'

'sdfopen e:\data\2005u.nc'
'set lon 130'
'set lat 0 50'
'set t 1 31'
'define a2=uwnd'
'close 1'

'sdfopen e:\data\2006u.nc'
'set lon 130'
'set lat 0 50'
'set t 1 31'
'define a3=uwnd'
'close 1'

'sdfopen e:\data\2010u.nc'
'set lon 130'
'set lat 0 50'
'set t 1 31'
'define a4=uwnd'
'close 1'

接下来想直接求a1,a2,a3,a4平均值，发现时间起点不一致，没法平均。
程序如下：
'define ad=（a1+a2+a3+a4）/4'
'set gxout shaded'
'cbarn'
'd ad'
'printim  e:\data\u.png white'
'print'
'set gxout contour'
'd ad'
'printim   e:\data\u.png white'
'print'
结果报错cannot plot data -all undefined values

试了很多方法，发现四个变量a1,a2,a3,a4不能一起运算，真不知道怎么办？

做个霸气的木头

我有个比较笨的方法你试试，就是把四个文件的数据写到一个grd数据里，写成连续的，就方便处理了

mofangbao

看来你还是自己思考过的。
楼上思路你可以试一试

gynxz

楼上推荐的方法还要用什么软件？Fortran？
只用grads运算，不能实现的吗？

river

你这个gs和你想画的东西一点都不一样啊，而且为什么要把一个时间的序列赋给一个变量呢，这又不是fortran，还能整个数组啊。楼主想画的类似于剖面图啊，可以参考一下有关的帖子吧

gynxz

river 发表于 2013-3-12 09:43

                               
登录/注册后可看大图

你这个gs和你想画的东西一点都不一样啊，而且为什么要把一个时间的序列赋给一个变量呢，这又不是fortran，还 ...

是想赋值完直接运算，我赋值完画过相应的图，赋值完的变量还是时间序列.
目前只会用grads,走进了胡同，自己不知道咋办了

gynxz

一份资料的时间序列可以画出来。
现在是把每个资料的时间都看成是第1天到第30天的，每天都做平均。出图的效果想是这样的：经度固定，横坐标为时间从1-30天，纵坐标是纬度

river

gynxz 发表于 2013-3-12 10:05

                               
登录/注册后可看大图

是想赋值完直接运算，我赋值完画过相应的图，赋值完的变量还是时间序列.
目前只会用grads,走进了胡同，自 ...

你这个要怎么平均啊，没看出来有平均的算法在里面啊。画这种图和剖面图一个道理都得先计算再画图。就是先计算平均，然后再固定经度，让时间和纬度变化。

'sdfopen e:\data\2003u.nc'
'sdfopen e:\data\2006u.nc'
'sdfopen e:\data\2005u.nc'
'sdfopen e:\data\2010u.nc'
'set t 1 31'
'define average=(uwnd.1+uwnd.2+uwnd.3+uwnd.4)/4.'

'set lon 130'
'set lat 0 50'
'set gxout shaded'
'd average'
'cbarn'
'set gxout contour'
'd average'
'printim   e:\data\u.png white'
;
试试吧，应该就是这样

gynxz

river 发表于 2013-3-12 12:50

                               
登录/注册后可看大图

你这个要怎么平均啊，没看出来有平均的算法在里面啊。画这种图和剖面图一个道理都得先计算再画图。就是先 ...

首先谢谢你。
'sdfopen e:\data\2003u.nc'
'sdfopen e:\data\2006u.nc’
'sdfopen e:\data\2005u.nc'
'sdfopen e:\data\2010u.nc'
到这一步时，我先试着画四个变量的图，uwnd.1是对的，画uwnd.2、uwnd.3、uwnd.4时，提示cannot plot data -all undefined values。
是不是grads没有识别出uwnd.2、uwnd.3、uwnd.4这三个变量的值，
我前面也这么处理过，出了问题，所以后来写成赋值给四个变量再处理。

river

gynxz 发表于 2013-3-12 22:25

                               
登录/注册后可看大图

首先谢谢你。
'sdfopen e:\data\2003u.nc'
'sdfopen e:\data\2006u.nc’

你把你的完整gs贴出来看看，如果资料本身没有问题的话，只能是你的gs有问题