关于maskout剔除界外数据遇到的问题

男紫汗

根据http://bbs.06climate.com/forum.php?mod=viewthread&tid=893我输出了中国的地图数据，分辨率1*1，出图后感觉阴影和边界不太重合啊，美观倒是不在乎主要是计算，欣赏的话还是cnbasemap比较好。可能版上已有相关内容，时间紧迫实在翻不动帖子了。恳请大家指点！谢谢。
'reinit'
'open D:\fortran\station\station\stationgrid.ctl'
'open D:\fortran\station\station\station.ctl'
'open D:\fortran\station\station\mask.ctl'
'set lon 70 140'
'set lat 15 55'
'set grads off'
'set grid off'
'set mpdset hires'
'define a=oacres(g,rain.2)'
'define a1=maskout(a,g-0.5)'
'define aa=smth9(maskout(a1,mask.3))'
'set ylevs 20 25 30 35 40 45 50 55'
'set gxout shaded'
'd aa'
*'cnbasemap aa'
*'run cbar_interp 1 1 0'
'cbar.gs'
'southsea_value -v1 aa'
'printim D:\fortran\station\station\fffff.png white'
*'printim D:\fortran\station\station\sta2_grid2.5.png white'
;

用了@四叶草的mask（6楼），出图如下

river

用cnbasemap试试，在论坛看看

四叶草

楼主你把这个语句'define a1=maskout(a,g-0.5)'删除看看效果，本人平时用maskout时，都是
'define a=oacres(g,rain.2)'
'define aa=smth9(maskout(a,mask.3))'
这两个语句就完事了。

男紫汗

river 发表于 2016-1-15 21:32
用cnbasemap试试，在论坛看看

嗯，画图我是用的cnbasemap，这是要计算的，先看看效果

男紫汗

四叶草 发表于 2016-1-16 10:01
楼主你把这个语句'define a1=maskout(a,g-0.5)'删除看看效果，本人平时用maskout时，都是
'define a=oacre ...

试了，一样的，如果阴影都在界限内，偶尔突出一点应该是正常，我这个东北，西北西南都大片突出。。你用的时候效果咋样？

四叶草

男紫汗 发表于 2016-1-16 10:36
试了，一样的，如果阴影都在界限内，偶尔突出一点应该是正常，我这个东北，西北西南都大片突出。。你用的 ...

我用maskout，阴影全部都是在边界内的，不会有超出的现象。你应该是mask文件制作环节出问题了，返回去仔细检查下，或者你用我的这个mask文件试试看，记得更改路径 mask.rar (449 Bytes, 下载次数: 31)

2016-1-16 10:42 上传

点击文件名下载附件

男紫汗

四叶草 发表于 2016-1-16 10:39
我用maskout，阴影全部都是在边界内的，不会有超出的现象。你应该是mask文件制作环节出问题了，返回去 ...

哈哈都在边界内了。。我那个是jave版meteoinfo做的。。另一个版的出图时有错误提示。另问海南那块有没有值得，还是说有值但没有阴影，还有东南沿海这块正常吗

四叶草

男紫汗 发表于 2016-1-16 11:22
哈哈都在边界内了。。我那个是jave版meteoinfo做的。。另一个版的出图时有错误提示。另问海南那块有没有 ...

海南那块应该是有值，因为精度问题，按默认的插值半径应该显示不出阴影，建议你换成cnbasemap试试

男紫汗

四叶草 发表于 2016-1-16 12:51
海南那块应该是有值，因为精度问题，按默认的插值半径应该显示不出阴影，建议你换成cnbasemap试试

能否再上传个meteoinfo

四叶草

男紫汗 发表于 2016-1-17 09:48
能否再上传个meteoinfo

本想直接给你个链接，但链接现在好像进不去，所以只能直接上传附件了。 MeteoInfo_1.1.3_setup.exe (36.55 MB, 下载次数: 61)

2016-1-17 14:38 上传

点击文件名下载附件