基于Visual C++的高斯投影坐标反算.pdfVIP

基于VisualC++的高斯投影坐标反算

许大亮

【摘要】文中利用VisualC++语言编写高斯投影坐标反算的程序，采用的数据是

西安80椭球参数和高斯坐标；首先根据子午线弧长公式迭代计算求出垂足纬度，

其次算出大地纬度和经差，最后可以得到大地坐标。程序计算出的坐标精度为

0.0001″。结果表明：用转换程序解算出的坐标存在一定的误差，但是可以满足低

精度的工程建设。%UsingVisualC++languagedoestheGaussprojective

coordinateinversecalculation.ThedataisXian80ellipsoidparametersand

gausscoordinates;first,accordingtothemeridianarclengthformulacan

getlatitudesoftheperpendicularfootwithiterativecomputation,

secondlycalculatethelatitudeandthedifferenceoflongitude,finallythe

geodeticcoordinatesareachieved.Theaccuracyofthecoordinatewhichis

calculatedbythecomputerprogramis0.0001″.Theresultsshowthatthe

conversionprogramhascertainerrorwhenitdoesgaussinverse

calculation,butitissuitableforengineeringconstructionoflowaccuracy.

【期刊名称】《科技创新与生产力》

【年(卷),期】2015(000)008

【总页数】3页(P59-60,63)

【关键词】VisualC++语言;坐标反算;高斯坐标;垂足纬度

【作者】许大亮

【作者单位】安徽理工大学测绘学院，安徽淮南232001

【正文语种】中文

【中图分类】TP311

当前，在我国常用的坐标系有北京54坐标系、西安80坐标系、WGS-84坐标和

2000国家大地坐标系[1-2]。高斯投影坐标转换成为大地坐标（高斯反算），其

中一个重要的问题就是如何准确、快速地计算出垂足纬度。因此，垂足纬度的准确

性决定着换算后坐标的精度。在垂足纬度的计算中，常采用迭代计算的方法。笔者

参考了一些文献、著作，利用VisualC++语言设计了一个程序，实现了高斯投影

坐标转换成西安80大地坐标，计算结果精确可靠，使用方便。

1.1高斯坐标反算公式

式中：Bf为垂足纬度；Mf为子午圈曲率半径[1]。其中Nf，tf，ηf等脚注为“f”

的均为Bf的函数，也就是说，它们都是将垂足纬度Bf的数值带入相应的公式后求

得[2]。

1.2子午圈弧长公式

1.3Mf，Nf，tf，ηf等关于Bf的公式

地球椭球的长半径为a，短半径为b。

第一偏心率

第二偏心率

t=tanB，η=ecosB.

为了简化公式的书写[2]，引入两个大地纬度B的函数符号

极曲率半径

子午圈曲率半径

卯酉圈曲率半径

2.1程序运行环境

该程序是基于VisualC++6.0环境下开发的控制台程序。

2.2程序运行流程

程序运行流程见图1。

2.3主要代码

编写doublegetBf（doublex）函数计算出垂足纬度，代码为

函数体中重复迭代直到

现在工程建设中一般要求（L，B）计算至0.01″和0.0001″，当要求精度为0.01″

时，可将式（5）中不等式右端的值换为10-6，可以按照具体要求进行改动[3]。

大地纬度B和经差的计算函数为

在VisualC++6.0环境下运行，按顺序输入高斯值坐标（x，y），然后单击回车

键运行出结果。程序先计算出垂足纬度，然后算得以垂足纬度为变量的子午圈曲率

半径，卯酉圈曲率半径等数值，将他们带入高斯坐标反算的公式计算出大地经度和

纬度。对于不同的大地椭球，只需要知道椭球参数，即可在高斯反算下进行计算。

由于计算误差的存在，该程序程序精确度为0.000