OpenCoordV1.2.0新增框架转换和EGM2008查询功能

Posted 2023-01-26 xiaokcehui

OpenCoordV1.2.0已经发布了，欢迎感兴趣的同学可以在xiaok海洋测绘网站下载使用，此版主要新增框架转换和EGM2008查询功能。软件目录下的帮助文档不仅有这两个功能使用的详细描述，而且新增了如何进行历元转换和WGS84与CGCS2000转换方法。

OpenCoordV1.2.0的主界面见下图。

​

1、框架转换的步骤

首先是框架转换的设置。

正在上传…重新上传取消​

然后是选择ITRF2014-ITRF2008框架转换，并将目标历元2010.0改成2011.0。

勾选“框架转换”，输入经纬度坐标，点击“转换”获取转换后的坐标，从转换结果可以看出，ITRF2014与ITRF2008框架坐标只差在秒的第5位小数，即mm位。

2、不同历元的转换

《2000中国大地坐标系_中国大陆速度场_魏子卿》文献给出了历元转换公式和大陆速度场。
速度场是现代大地测量不可或缺的基础性数据。GNSS定位得到的坐标X(t)是当前观测历元t的坐标。为了将观测历元ｔ的坐标归算为参考历元t0的坐标X0，或者为了将不同历元的坐标互相转换，均用速度Vx按下式进行历元变换
X(t)=X0+Vx(t-t0)
中国大陆3°×3°的格网平均速度见下表。不同历元的转换见第3部分。

3、WGS84与CGCS2000的坐标转换

WGS84是动态坐标，每隔一段时间就与新的ITRF框架对准一次，CGCS2000是静态坐标，对准的是ITRF97框架，2000.0历元，由于WGS84和CGCS2000的椭球参数基本一致，且它们的框架差别也不大（cm级），WGS84与CGCS2000的坐标转换主要是历元转换。
下面以2022年的WGS84坐标36°30’00”N，120°30’00”N转成CGCS2000坐标为例说明WGS84与CGCS2000的坐标转换。
查询速度为：-27.5,-7.9,-11.7，单位：mm/a
WGS84转换CGCS2000三参数为：（-27.5*-22*0.01=0.605）,（-27.5*-22*0.01=0.1738）,（-27.5*-22*0.01=0.2574）
转换过程如下：
（1）设置转换参数

（2）布尔沙七参数转换

（3）验证结果

验证的方法就是将经纬度转换空间直角坐标，然后比较差值。

WGS84空间直角坐标为：-2605291.835，4422907.864，3772934.548

CGCS2000空间直角坐标为：-2605291.230，4422908.037，3772934.805

WGS84-CGCS2000差值为：-0.605，-0.174，-0.257，与上面计算的转换参数相同。

显然，转换参数和结果正确。

4、查询EGM2008

首先下载EGM2008模型文件，下载地址为
链接：https://pan.baidu.com/s/1QsLjvwcK5PUVry3yO93c6g?pwd=nxdh
提取码：nxdh
然后将此压缩文件geoids.rar解压到软件目录pulish下，
接着设置环境变量：GEOGRAPHICLIB_DATA，

最后是运行程序。

点击查询按钮，查询水准面大地高（高程异常），计算出水准高和85高。

与EGM2008官方程序查询结果（见下图）相比，差别1mm左右。