求助：怎样用MeteoInfo制作海洋的maskout文件？

dqxiaoding

我想用grads画图只画黄海地区，想要把山东省陆地屏蔽掉。
在grads中用basemap命令的话，海岸线的ASC文件分辨率太低太粗糙，不符合我的要求。
所以我想用MeteoInfo制作海洋的maskout文件。

但我看了本版的帖子，只有介绍制作省份的maskout文件，没有制作海洋部分的。
请问该怎样制作海洋部分的maskout文件？
另外，我下载了最新版本的1.4.3  JAVA版的，试着做山东省的maskout文件，
点了select之后鼠标就变成圆圈一直在转，是我的操作有问题吗？还是出文件就是要等很长时间？

MeteoInfo

参考此贴：http://bbs.06climate.com/forum.php?mod=viewthread&tid=2808

dqxiaoding

MeteoInfo 发表于 2017-5-9 12:48
参考此贴：http://bbs.06climate.com/forum.php?mod=viewthread&tid=2808

谢谢版主。我试试

dqxiaoding

MeteoInfo 发表于 2017-5-9 12:48
参考此贴：http://bbs.06climate.com/forum.php?mod=viewthread&tid=2808

版主好！继续向您求助。我试了您说的通过属性数据选择，我要画的区域是黄海。
开始我用了bou2_4p.shp，name<>山东省。得到的maskout文件画出来后江苏的部分有值，但山东和黄海没有值。
后来我改用country1.shp,colormap<>2,选中的地区是中国以外的【陆地】，选不中海洋。因为我的数据范围就是山东和黄海，这个画出来全是空值。
所以我觉得这个地图是不是只包含陆地的？还是我选择的方法不对？
我的gs文件中的maskout语句是这样的，应该没有问题'd smth9(maskout(p.2,mask.1(t=1)))'

dqxiaoding

这个问题我已经解决了。
MeteoInfo生成的dat文件用grads绘图后就可以发现，dat中只有两种数值，分别是1和-1。1就是我们在画图中显示的部分。-1是不画的部分。
如下图就是我在楼上提到的我做出的mask文件，画图效果是只画出江苏的数值，山东和黄海不画。
我想要只画黄海的话，只需要在MeteoInfo中选中中国地图，输出mask文件（这个mask文件的效果是只画中国陆地），然后用fortran读取dat文件，乘以-1后输出，再按照maskout方法的gs文件画图，就可以只画海洋不画陆地了。
maskout方法其实就是以格点为单位，将不画的格点舍弃或者屏蔽掉。本质上跟这个帖子是相同的，不过制作mask文件的优点是可以很方便的控制mask文件的空间分辨率和经纬度范围，使其与要绘图的文件一致。http://bbs.06climate.com/forum.php?mod=viewthread&tid=4780&extra=&page=1
但maskout与basemap相比，缺点就是maskout是屏蔽掉格点，所以绘图的边缘也是方块形状的、不平滑的，跟地图的省界或者国界不会完全贴合。我看到论坛里也有人贴出自己画的图，线条在靠近省界的部分有留白。
如果能把province-basemap.gs改一下，让province-basemap能像basemap那样有不画陆地的选项，或者把gs改成不画陆地。就可以更好的满足我的需求，治愈我的强迫症了。可惜我精力有限，有时间可以研究下。

dqxiaoding

用上面的方法画出来的效果图是这样的，也是不大理想。