基于VisualC++的高斯投影坐标反算许大亮【摘要】文中利用VisualC++语言编写高斯投影坐标反算的程序，采用的数据是西安80椭球参数和高斯坐标；首先根据子午线弧长公式迭代计算求出垂足纬度，其次算出大地纬度和经差，最后可以得到大地坐标。程序计算出的坐标精度为0.0001″。结果表明：用转换程序解算出的坐标存在一定的误差，但是可以满足低精度的工程建设。%UsingVisualC++languagedoestheGaussprojectivecoordinateinversecalculation.ThedataisXian80ellipsoidparametersandgausscoordinates;first,accordingtothemeridianarclengthformulacangetlatitudesoftheperpendicularfootwithiterativecomputation,secondlycalculatethelatitudeandthedifferenceoflongitude,finallythegeodeticcoordinatesareachieved.Theaccuracyofthecoordinatewhichiscalculatedbythecomputerprogramis0.0001″.Theresultsshowthattheconversionprogramhascertainerrorwhenitdoesgaussinversecalculation,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，η=e'cosB.为了简化公式的书写[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.0001″。只适合较低精度的生产建设和工程设计。综上所述，利用计算机编程语言VisualC＋＋编写高斯投影反算程序，计算方便快捷、精度较高，具有很强的实用性。通过改变程序中的大地椭球参数可以计算将高斯坐标反算成为北京54坐标、WGS-84或2000国家大地坐标。【相关文献】[1]张华海,王宝山,赵长胜,等.应用大地测量学[M].徐州:中国矿业大学出版社,2008.[2]王苑,王春青,杨宏海,等.