站点资料画图并屏蔽掉省外的等值线出错了

凝蓝_@薇微笑

站点资料生成了二进制文件并且生成了map文件以及格点ctl,用meteoinfo输出了mask文件的ctl,想画出华中三省的气温填色图，但是图出错了，大家帮忙
看看哪里有问题，多谢了：

站点气温ctl(zhandian.ctl):
dset e:/data/huazhongzhandian/qiwendata.grd
dtype station
stnmap e:/data/huazhongzhandian/1.map
undef -9999.0
title station data sample
tdef 1 linear 01Z01may2011 1hr
vars 1
t  0 99 tem data
endvars
格点ctl(gedian.ctl):
dset e:/data/huazhongzhandian/grid.grd
title grid data sample
undef -9999.0
XDEF 10 LINEAR 108.35 1.0
YDEF 13 LINEAR 24.66 1.0
ZDEF 1 LEVELS 1000
TDEF 1 LINEAR 01Z01may2011 1hr
vars 1
g 0 99 grid data
endvars
华中三省mask文件的ctl（huazhong.ctl）：
DSET e:\data\huazhongzhandian\huazhong.dat
TITLE Mask data
UNDEF -9999.0
XDEF 10 LINEAR 108.35 1
YDEF 13 LINEAR 24.66 1
ZDEF 1 LINEAR 1 1
TDEF 1 LINEAR 01Z01may2011 1hr
VARS 1
mask 0 99  background mask data
ENDVARS
GS文件：
'reinit'
'open e:\data\huazhongzhandian\gedian.ctl'  
'open e:\data\huazhongzhandian\zhandian.ctl'   
'open e:\data\huazhongzhandian\huazhong.ctl'
'set grads off'
'set grid off'
'set frame off'
'set xlab off '
'set ylab off '
'set lon 108.35 117.35'
'set lat 24.66 36.66'
'define a=oacres(g,t.2)'
'define a1=maskout(a,mask.3)'
'set gxout shaded'
'd a1'
'set gxout contour'
'd a1'
'cbarn'
'printim e:\data\huazhongzhandian\sta_grid1.png  white'

男紫汗

似乎有maskout就不需要格点了

river

你先不用mask，直接画看看是个什么情况

凝蓝_@薇微笑

river 发表于 2017-5-24 21:03
你先不用mask，直接画看看是个什么情况

没用mask是这样子的，我想抠出来华中三省的（河南湖北湖南），谢谢帮忙看看哪里出了问题

river

凝蓝_@薇微笑 发表于 2017-5-25 08:39
没用mask是这样子的，我想抠出来华中三省的（河南湖北湖南），谢谢帮忙看看哪里出了问题

这么看其他地方问题不大，可能是你的mask文件制作的有问题