求删除，谢谢

李超宇

karenlk 发表于 2017-4-15 19:16
例如
*30年平均78月
*'define sst2=(ave(sst.1,t=1531,t=1879,12)+ave(sst.1,t=1532,t=1880,12))/2'

非常感谢，你写的这个我看懂了，是30年78月的平均，对我之后的画图非常重要。不过我还想得到每一年区域内的平均，然后画出柱状图，类似下图那样子

karenlk

李超宇 发表于 2017-4-15 20:07
非常感谢，你写的这个我看懂了，是30年78月的平均，对我之后的画图非常重要。不过我还想得到每一年区域内 ...

那你把12个月求出来，平均，再进行区域平均

常山赵子龙

这个资料的缺省值需要注意啊。

李超宇

常山赵子龙 发表于 2017-4-16 22:19
这个资料的缺省值需要注意啊。

要怎么处理？？

李超宇

karenlk 发表于 2017-4-15 22:17
那你把12个月求出来，平均，再进行区域平均

'reinit'
'sdfopen F:\grads\4\sst.mnmean.v4.nc'
'set grid off'
'set grads off'

'set lon 120 165'
'set lat 22 46'
'define sst2=(ave(sst,t=1531,t=1879,12)+ave(sst,t=1532,t=1880,12))/2'
'define sst3=tloop(aave(sst2,lon=120,lon=165,lat=22,lat=46))'
'd sst3'
'draw title sst'
'gxprint F:\grads\4\sst.mnmean.v4.1.png white'

你的意思是这样吗？？可是这样什么也画不出来呀

常山赵子龙

常山赵子龙 发表于 2017-4-16 22:19
这个资料的缺省值需要注意啊。

我没有用过grads处理数据，我记得这个缺省值貌似不是NAN，我用的matlab将这里的缺省值赋值成NAN然后进行的处理，结果就对了。因为之前我也出现过你的情况，就是一条直线。

lgz

缺省值需要注意

karenlk

李超宇 发表于 2017-4-16 23:23
'reinit'
'sdfopen F:\grads\4\sst.mnmean.v4.nc'
'set grid off'

set lon 120 165'
'set lat 22 46'
你前面设置的是这个范围，二维的图
但是你最后出图d sst3的时候，没有再重新设定范围，因为你想出的图是一维的。。。。
要重新设定一下

李超宇

karenlk 发表于 2017-4-17 10:45
set lon 120 165'
'set lat 22 46'
你前面设置的是这个范围，二维的图

'reinit'
'sdfopen F:\grads\4\sst.mnmean.v4.nc'
'set grid off'
'set grads off'

'set lon 120 165'
'set lat 22 46'
'define sst2=(ave(sst,t=1531,t=1879,12)+ave(sst,t=1532,t=1880,12))/2'
'set lon 120'
'set lat 22'
'set lev 0'
'set t 1531 1879'
'define sst3=tloop(aave(sst2,lon=120,lon=165,lat=22,lat=46))'
'd sst3'
'draw title sst'
'gxprint F:\grads\4\sst.mnmean.v4.1.png white'
这样子吗？？
然后出来就成一条直线了

karenlk

李超宇 发表于 2017-4-17 11:23
'reinit'
'sdfopen F:\grads\4\sst.mnmean.v4.nc'
'set grid off'

你这还是不行，t是月值，你最后想画的年的，要提前处理。。。还是先提取个grd计算一下，然后配个ctl再画吧。。。这种图我是用matlab处理的