请大家帮个忙，看看我的图为什么这样？

漂洋过海

river 发表于 2013-4-9 18:23

                               
登录/注册后可看大图

把画图gs里的'set cmin 2'都去掉试试吧

数据是负的

river

漂洋过海 发表于 2013-4-9 15:52

                               
登录/注册后可看大图

谢谢你的回复，麻烦你看看我画图用的ctl和gs文件。
这是插值后格点文件的ctl:
dset  H:\drought\1981-2 ...

把画图gs里的'set cmin 2'都去掉试试吧

漂洋过海

lqouc 发表于 2013-4-9 16:22

                               
登录/注册后可看大图

其实看着这个画图的gs也没什么问题，但是确实缺的那三个纬度带很奇怪，你有没有看过插值之后的某一时刻的 ...

.恩，你说的对，我画画看看，谢谢啦！

lqouc

漂洋过海 发表于 2013-4-9 16:22

                               
登录/注册后可看大图

谢谢你的回复，看来我对缺测这个数值的理解还有欠缺~

太客气了，我也不擅长grads，祝你好运啦。

lqouc

漂洋过海 发表于 2013-4-9 15:52

                               
登录/注册后可看大图

谢谢你的回复，麻烦你看看我画图用的ctl和gs文件。
这是插值后格点文件的ctl:
dset  H:\drought\1981-2 ...

其实看着这个画图的gs也没什么问题，但是确实缺的那三个纬度带很奇怪，你有没有看过插值之后的某一时刻的空间分布是什么样？要是你的插值场就有缺失的纬度带，那就应该是插值的问题。话说我也没做过降水，你的数据既然没有缺测，那么为什么还要把那些小值去掉呢。

漂洋过海

冰点 发表于 2013-4-9 15:58

                               
登录/注册后可看大图

第一张图标题Prep 110E-120E 是绘图的时候对纬向求了平均，作纬向平均降水随时间和经向的变化图么？
...

谢谢你的回复，看来我对缺测这个数值的理解还有欠缺~

漂洋过海

lqouc 发表于 2013-4-9 14:56

                               
登录/注册后可看大图

你画图的gs呢，不能就断定是插值不对吧，很可能是你画图的问题，或者是差值的数据ctl有问题？就说这个插值的 ...

我把maskout去掉之后画图结果是高纬度的那两处断面连续了，但是25°N下面的还是没有出现.....
非常感谢你的指点，我再继续试试看~！

漂洋过海

冰点 发表于 2013-4-9 15:58

                               
登录/注册后可看大图

第一张图标题Prep 110E-120E 是绘图的时候对纬向求了平均，作纬向平均降水随时间和经向的变化图么？
...

是做径向平均的降水岁时间和纬度变化图。

冰点

      第一张图标题Prep 110E-120E 是绘图的时候对纬向求了平均，作纬向平均降水随时间和经向的变化图么？
一年365天分73候，也就是你的73个时次。每个时次各站点均有数据，要把全部时次的数据绘制到一张图上还真不简单啊！
      你说数据资料没有缺测值，'define b=maskout(a,prc.1-0.5)'这句话不就说当prc.1-0.5<0时对应格点上的a值就赋值为缺测值么，另外1号文件中的变量直接写出来就是了，没必要再加个.1。
     

漂洋过海

lqouc 发表于 2013-4-9 14:56

                               
登录/注册后可看大图

你画图的gs呢，不能就断定是插值不对吧，很可能是你画图的问题，或者是差值的数据ctl有问题？就说这个插值的 ...

谢谢你的回复，麻烦你看看我画图用的ctl和gs文件。
这是插值后格点文件的ctl:
dset  H:\drought\1981-2010pre\pentad\30ave-pentad-grid.grd
undef 3276.6
title sample data
XDEF 141 LINEAR 70 0.5
YDEF 81 LINEAR 15 0.5
zdef 1 linear 1000 1
tdef 73 linear 1jan1981 5dy
vars 1
r 0 99 grid data prepared oacres function
endvars
画图的gs:
'reinit'
'enable print H:\drought\1981-2010pre\pentad\pentad110-120E.gmf'
'open H:\drought\1981-2010pre\pentad\30ave-pentad-grid.ctl'
'set grads off';'set grid off'

'set parea 1 10.5 2 7.5'
'set xlopts 1 6 0.15'
'set ylopts 1 6 0.15'
'set ylint 5'


'set lat 20.5 45'
'set t 1 73'
'set z 1'
'set x 1'
'set xyrev on'
'define p=ave(r,lon=110,lon=120)'
'set gxout shaded'
'set cmin 2'
'set cint 2'

'set clevs 2 4 6 8 10 12 14 16'
'set ccols 0 14 11 3  7 12 8 2 6'
'd p'
'set gxout contour'
'set ccolor 1';'set cthick 5'
'set cmin 2'
'set cint 2'
'd p'

'cbarn 1 0 5.8 1.3'
'set line 0'
*'draw recf 0.71 1.50 1.25 1.71'
'draw title Prep   110E-120E  '

'print'
'disable print'
;