0. 码字不易，点赞加关注（公众号：WZZHHH，后期代码将在公众号可以下载），使用请注明出处（根据我的研究方向，我会不断更新代码）。 1. 读取gamit解算的met文件，根据前缀met_读取相应文件下met文件，met文件的列分别为1为year、2为doy（年积日）、3为小时、4为分钟、5为秒、6为ZTD、7为ZWD、8不知道、9为PWV、11为气压、12为温度、13为ZHD。 2

0. 码字不易，点赞加关注（公众号：WZZHHH，部分资料在公众号可以下载） Part1 垂直补偿的原因 我们在获取GNSS站点位置的PWV有很多方法，其中一个是利用ERA5数据进行积分得到。由于ERA5格网位置和GNSS不并址，需要进行水平方向和垂直方向的补偿的补偿获得。这里讲的是垂直方向补偿，水平方向补偿已经写过反距离加权水平内插，附matlab代码(ERA5和GNSS站点不并址的处理