提取中国范围内的格点数据

lysx

众所周知，清风大侠各种编程，软件设计能力都很出色，为气象家园的建设发展做出了不可磨灭的贡献。每一次看清风的帖子，我都能学到不少东西。再次，真诚致谢！！！

因为项目研究需要，必须将中国范围以外的数据设置成undef，然后才能进行下一步工作。于是乎，本人遍历了家园中fortran和grads下的大部分帖子，找到了清风的帖子：fortran提取任意区域内的点，可用来生成GrADS的mask文件。http://bbs.06climate.com/forum.p ... 18693&fromuid=21651


清风的程序中用到了一个子程序，很牛的，用来判断某个站点是否属于多边形区域内。看到这里，我觉得需要有更简单的方法作为替代。由于地形数据需要用MeteInfo生成，然后就拜读了传说中的谁的帖子：maskout方法的使用以及maskout文件制作
http://bbs.06climate.com/forum.php?mod=viewthread&tid=893&fromuid=21651
(出处: 气象家园)，学会了生成地形文件。和传说中的谁一样，我将地形文件保存成为二进制格式。这样，保存的二进制文件中的所有为1的元素表示中国范围内，为-1的元素为中国范围外的。

接下来，我用fortran读入该地形文件以及所研究的数据，分别存放在数组china和array中。然后仅用一条语句就把中国范围以外的数据给屏蔽掉了:
WHERE ( china .EQ. -1 )
            array = undef
ENDWHERE
看看俺的图吧


很简单吧！！！如果您也按照我的方法试成功了，请您点个赞PS :: 我的mask文件china.dat的数据范围是从70-136,15-55°，分辨率为0.25°。如果您的数据分辨率为0.25°的整数倍，不想学习MeteInfo生成地形文件也可以。只要在fortran中度取数据后，每个几个数进行设置为undef就可以了。下面上传我的地形文件。

river

感谢分享，好方法

aa22556688

给力，赞！

topmad

这样的流程介绍的帖子很好 一步一步来

lysx

topmad 发表于 2014-9-12 07:51
这样的流程介绍的帖子很好 一步一步来

感谢您的鼓励！

更大华

好贴

黑马骑士

有点空闲就发帖。。。

言深深

很不错的一份工作！

mofangbao

不错的思路~我之前写的那个Fortran其实只能应对简单的多边形，复杂的就不行了，MeteoInfo应该实现的比较完善了~

lysx

mofangbao 发表于 2014-9-12 09:46
不错的思路~我之前写的那个Fortran其实只能应对简单的多边形，复杂的就不行了，MeteoInfo应该实现的比较完 ...

不过你的子程序太牛了！