河南理工大学本科毕业论文

图4.2大地坐标向空间直角坐标转换

2．空间直角坐标系之间的七参数转换框图设计

图4.3空间直角坐标系之间的七参数转换框图设计

3. 空间直角坐标向大地坐标转换

21

河南理工大学本科毕业论文

开始 1 输入X、Y、Z 输入椭球参数a、b 取B的迭代初值 计算椭球第一偏心率e B1?arctanZ?X2?Y2? 1 2 3 2 否 ?Z?Ne2sinB?按式B?arctan?? ??22?X?Y?B??B?0.00001? 迭代计算B值 是 3 B?B?输出B、L、H的值 ?Y?L?arctan???X?ZX2?Y22H??N1?e或?NsinBcosB??结束

图4.4空间直角坐标向大地坐标转换

22

河南理工大学本科毕业论文

4.2.2同一椭球参数下大地坐标与高斯平面坐标的转换框图设计 1. 高斯正算

图4.5 高斯正算

2. 高斯反算

图4.6 高斯反算

23

河南理工大学本科毕业论文

5 实例应用

本章从同一参考椭球下不同坐标形式之间的转换和不同参考椭球之间的坐标转换两个方面出发，结合具体工程算例，验证了上一章所述基于MATLAB编程实现坐标转换的可行性和准确度。 5.1大地坐标与空间直角坐标间换算实例

已知珠穆朗玛峰峰顶B?27?59?16.94241??，L?86?55?31.72137??，以及雪面的大地高H?8821.4016m，其坐标系为1954年北京坐标系。先用本文提供的m文件进行大地坐标向空间直角坐标换算，并与中海达坐标转换模块计算的相应值进行了比较，操作步骤如下：

1.在MATLAB系统中运行本文附录中提供的坐标转换的大地坐标转为空间直坐标的dd2kj.m，（输入时以北京54椭球的长轴a、扁率Bl以及珠穆朗玛峰峰顶大地坐标为参数）如图5.1所示；

2.得到大地坐标转为空间直角坐标的计算结果（图5.2）。

图5.1 本文程序运行

图5.2 MATLAB计算结果

3.用中海达软件中的坐标转换模块进行大地坐标向空间直角坐标的换算，界面如图5.3所示。

24

河南理工大学本科毕业论文

图5.3中海达软件大地坐标向空间直角坐标换算界面

将换算结果整理如表5.1所示：

表5.1与中海达软件从原大地坐标向空间直角坐标的换算结果 单位：m 本文换算值 本文换算与中海达软件换算差值 中海达软件换算值（X，Y，Z） （X，Y，Z）302726.854415229 5636102.39017883 2979527.61924339 302726.854415 5636102.390179 2979527.619243 0.000000229 0.00000017 0.00000039 再将计算出的空间直角坐标值用本系统和中海达软件的坐标转换模块分别换算为大地坐标并分别和原大地坐标对比，操作步骤如下；

1.在MATLAB系统中运行本文附录中提供的坐标转换的空间直角坐标转为大地坐标的kj2dd.m，（输入时以北京54椭球的长轴a、扁率Bl，以及上述换算得到的空间直角坐标为参数），得大地坐标值，如表5.2所示；

2.用中海达软件坐标转换模块进行空间直角向大地坐标的换算，界面如图5.4所示。

图5.4中海达软件大地坐标向空间直角坐标换算界面

将换算结果整理在表5.2：

25