[原创]VMF1模型使用

    VMF1模型（Vienna Mapping Function）是目前模拟对流层状况精度最高的模型，主要发展了基于测站的VMF1（S-VMF1）和基于格网的VMF1（G-VMF1），我们主要使用后者，因为其没有任何测站限制，只需输入所需参数即可得到全球任何地点的大气参数。

    在VMF1网站中提供了包括算法和格网数据在内的所有资料。（http://ggosatm.hg.tuwien.ac.at/DELAY/）。

    VMF1模型的输入参数包括：ah,aw,dmjd,dlat,ht,zd（干湿参数，简化儒略日，纬度，大地高，天顶距），其中干湿参数可以通过GPT2模型获得也可以通过VMF1模型网站中提供的6h分辨率的格网数据内插得到。作为估计大气状况，进行精密单点定位的重要模型，VMF1模型输出结果为HMF,WMF（干湿映射函数）进而可以计算倾斜路径水汽含量。

    在网站中提供的格网数据中包括了ah,aw,ZHD,ZWD参数，其格网分辨率为2°*2.5°。当我们单纯使用格网内插获得大气参数时就要考虑内插的问题，文献中提出使用双线性内插方法来获得网格内部某点的大气状况参数。本人参考了维基百科中双线性内插方法的使用说明来实现相应的算法（https://en.wikipedia.org/wiki/Bilinear_interpolation）。

    需要注意的是，在内插ZHD,ZWD时，需要将格网高程换算为计算站点的大地高程，在文献（http://link.springer.com/content/pdf/10.1007/s00190-007-0170-0.pdf）中给出了具体的转换方法。此外，网站中提供了格网点高程数据（orography_ell），其数据结构与格网数据向对应，初次学习时没有弄明白其数据排列的规律，经过与格网数据的仔细比较后豁然开朗，现将格网点高程数据与格网数据对比如下：

 1     90.00   -90.00      0.00    360.00      2.00      2.50 
 2        9       9       9       9       9       9       9       9       9       9
 3        9       9       9       9       9       9       9       9       9       9
 4        9       9       9       9       9       9       9       9       9       9
 5        9       9       9       9       9       9       9       9       9       9
 6        9       9       9       9       9       9       9       9       9       9
 7        9       9       9       9       9       9