判断某个坐标点是否在研究区域范围内

前言

       这里主要对采用Meteoinfo这款软件来进行操作和判断，通过fortran程序，夹角之和检验法也可以实现，这里不做过多的介绍。

正文

        对这款软件的了解源自气象论坛（氛围极好的一个论坛），此判断方法在N年前本人运用的非常熟练，但现在基本上忘记的也差不多了，只是这段时间，帮忙几个本科生完成毕业设计，进行相关数据处理，这里大概记录一下，下次再用的时候直接照着这个来就可能会方便很多。

        这里主要针对Meteoinfo软件下的"MIScript.exe"可执行程序来实现，"C:\Program Files (x86)\MeteoInfo\Script"，脚本通常在该文件夹下。既然我们是要判断某个点是否在研究区域范围内，或者换个说法，我们需要提取出某个边界范围内的所有坐标点。那么第一步，我们需要知道研究区域的边界文件，对于如何获取可以通过Surfer的数字化功能来实现。此外还需要简单提醒的是，在Script文件夹下，脚本通常以".py"后缀名结束，".wmp"通常为研究区域的边界文件，只需要是边界即可，文件格式如下：

 

        我们通常按照此格式手动修改此文件即可。

        比如我们希望研究区域为华北地区，研究区域如图：

      如图红色区域，是由319个点组成，我们生成northchina.wmp文件。接下来需要建立.py文件，或者在原来文件基础上进行修改，具体内容我们这里不表。

      那么，接下来最重要的就是就是生成研究区域范围内的所有点，比如我们需要红色线条框起来的所有经纬度，在Meteoinfo文件夹下有个小工具"MIScript.exe"，我们主要通过其来实现。在此目录下打开dos窗口，通过输入cmd切入正常模式。

       我们通过如下命令：MIScript.exe *.py > name.txt  即可实现研究区域坐标点的提取，最终效果如下图：

    至此，我们目的已经达到。通过已知边界来提取边界范围内点位的坐标。