用C语言实现地图坐标投影的转换.doc

南 昌 工 程 学 院毕 业 设 计（论文） 水利与生态工程学院 系（院） 测绘工程 专业毕业设计（论文）题目 高斯-克吕格投影变换的程序实现 学生姓名 卢 浩 桢 班 级 测 绘 工 程 学 号 2008101405 指导教师 张 红 梅 完成日期 2012 年 6 月1日高斯-克吕格投影的程序实现The conversion of Gauss-Kruger projection in C language总计 毕业设计（论文） 25 页 表 格 2 个插 图 6 幅摘 要论文是对用C语言实现地图投影转换的详细论述，前两章是对C语言和地图投影的相关知识进行简单介绍，并将研究地图投影转换的目的、意义和研究现状进行阐述，让我们对相关知识有了深刻认识与了解；后两章是具体实现的程序代码，程序运行实现的截图以及具体数据进行计算对程序进行验证，说明其可行性与正确性，得出最后的结论。关键词：地图投影 C语言 高斯-克吕格投影 高斯正反算 AbstractThis dissertation is to research the Issue of Map projection conversion achieved with c language in detail. The first two chapters introduce us the C language and the knowledge of map projection briefly, and show the purpose, meaning and situation of study of map projection transformation one by one, which give us a deep understanding of this knowledge; the last two chapters describe the programming source code, and the pictures and specific cases and data in order to make sure its feasibility and correctness, finally , draw a conclusion.Keywords: Map projection; C language; Gauss-Kruger projection; Gauss formula for the positive and negative目 录摘要IABSTRACT第一章 引言11.1地图投影简介11.2研究地图投影转换的目的和意义11.3地图投影研究概括及趋势综述11.4研究思路和技术方法21.5 C语言概述3第二章 高斯-克吕格投影（高斯投影）42.1高斯投影的基本概念42.2高斯投影坐标正反算公式7第三章 c语言程序代码133.1地图投影变换的计算程序流程图133.2 C语言程序源代码153.3程序运行截图19第四章 算例21结 语23参考文献24致 谢25南昌工程学院本科毕业论文第一章 引 言1.1地图投影简介 地图投影,Map Projection.，按照一定的数学法则，把参考椭球面上的点、线投影到可展面上的方法。地图投影就是指建立地球表面（或其他星球表面或天球面）上的点与投影平面（即地图平面）上点之间的一一对应关系的方法。即建立之间的数学转换公式。它将作为一个不可展平的曲面即地球表面投影到一个平面的基本方法，保证了空间信息在区域上的联系与完整。这个投影过程将产生投影变形，而且不同的投影方法具有不同性质和大小的投影变形。1.2研究地图投影转换的目的和意义地图投影变换是指从一种地图投影点的坐标变换为另一种地图投影点的坐标。地图投影总类繁多，形式各异。在实际应用中，地理信息系统的数据来源各不相同，很难保证这些数据的投影方式一致。为了对这些不同来源的地理数据进行有效的处理跟分析，就必须通过地图投影变换把这些投影方式不同的数据定位在同一种坐标系下，使其具有相同的地图投影。同时，在地图编制过程中，也经常会遇到地图资料与新编地图之间投影不一致的情况，必须进过投影转换将该地图资料转绘到新编地图投影坐标网格中。人类一切经济活动都离不开地理空间，各类专业信息都必须以地形基础信息为空间载体，所以必须研究地图数据库中数字化地图数据处理、空间信息定位和变换，以满足各类专业信息系统建设的需要。综上所述，采用C语言作为软件开发平台，实现各种常见地图投影之间的任意变换，也具有很强的现实意义。1.3地图投影研究概括及趋势综述 地图投影及其变换模型已渗透到地理信息系统的各个方面，无论是数据的获取、预处理、信息的存贮记录，还是数据的处理、应用和输出都需要有一个空间定位的框架，即共同的地理坐标和平面坐标系统。由此可见地图投影在地理信息系统中的重要作用。另外，在地图编制过程中也常遇到原始资料与新编地图投影不一致的情况，此时也需要通过投影变换来解决此类问题。传统的投影变换方法均有工序繁、速度慢、1第一章 引言精度低等缺点，已不能适应制图自动化的要求。随着计算机技术和数学手段的高速进步，地图的编制以实现了完全自动化。投影变换计算自动化要求建立两种不同投影方式之间点与点的变换关系式，通过各种关系式编程实现各种投影之间的任意变换，不再需要人工参与变换，这样既提高了变换精度又提高了工作效率。目前国内外很多地图制图和地理信息系统软件，例如国外的Arcgis，MapInfo，MGE等，国内的MapGis，CityStar等，都已不同程度地提供了地图投影变换功能模块，因此，采用C语言工具软件实现各种常见的地图投影之间的任意变换，为地图投影变换方法研究提出一套行之有效的方法。1.4 研究思路和技术方法研究思路：1、首先弄懂地图投影转换原理。地图投影的实质是用数学的方法实现球面上点到平面上点的转换。地图投影的分类很多，每种投影都有一套坐标计算公式。 “高斯-克吕格投影”“UTM投影”“兰勃特等角投影”等投影之间的转换公式应熟练掌握。地图投影变换方法主要有以下两类：（1）传统地图投影变换方法，包括网格转绘法，蓝图拼贴法和纠正仪法等。这些方法都有工序繁杂，速度较慢，精度较低的缺点，难于男足制图自动化的要求。 （2）数字地图投影，包括解析变换法、数字变换法和数字解析变换法。其数字投影变换方法的主要思想就是通过一定的数学手段，寻找两个不同投影性质地图之间点位一一对应关系式。变换法是根据原始地图的投影方程式反解原始地图投影点的经纬度，再代入新编图的地图投影方程式，得到两种投影的平面直角坐标关系；又如建立两种不同地图资料间相应点坐标的直接关系，即在两种地图上量出相应点的平面直角坐标值，带入逼近多项式，分别组成线性方程组，求出系数值，再代入逼近多项式，即可求出新编图投影点的坐标。第二种方法具有普遍性，适用于地图投影转换。假定原始地图资料上点的坐标为（x，y），新编地图上对应点坐标为（X,Y），则对应投影点间坐标转换的基本关系式为X=F1(x,y)Y=F2(x,y) (1.4.1)F1,F2皆为制图区域内单值、连续的函数。对于任意两种不同投影性质的地图资料来说，其投影关系式分别为2南昌工程学院本科毕业论文x=f1(L,B)y=f2(L,B) (1.4.2)X=f3(L,B)Y=f4(L,B) (1.4.3)由（1.4.2）式反解得L=L(x,y) B=B(x,y) （1.4.4）把（1.4.4）式代入（1.4.3）式，便有X=f3L(x,y),B(x,y)Y=f4L(x,y),B(x,y) （1.4.5）上式即为地图投影转换的数学通式，使用于任何两种投影间变换函数f3，f4的具体表达式随投影间转换类型和转换方法的不同而不同。2、选择某种投影方式，构建投影变换公式。3、用编程语言实现地图投影转换。1.5 C语言概述C语言是一种计算机程序设计。它既具有高级语言的特点，又具有汇编语言的特点。它于1972年推出，1978后已先后被移植到大、中、小及微型机上。它既可以作为工作系统设计语言，编写系统应用程序，也可以作为应用程序设计语言，编写不依赖计算机硬件的应用程序。它的应用范围广泛，具备很强的数据处理能力，在软件开发上、各类科研都需要用到C语言，适于编写系统软件，三维，二维图形和动画等。之所以选择使用C语言进行地图投影的转换是因为它有以下特点：简洁紧凑、灵活方便 、运算符丰富 、数据类型丰富 、C是结构式语言、语法限制不太严格，程序设计自由度大、允许直接访问物理地址，对硬件进行操作、生成目标代码质量高，程序执行效率高 、适用范围大，可移植性好等。3第二章 高斯-克吕格投影（高斯投影）第二章 高斯-克吕格投影（高斯投影）2.1高斯投影的基本概念高斯投影又称横轴椭圆柱等角投影，是德国测量学家高斯于18251830年首先提出的。实际上，直到1912年，由德国另一位测量学家克吕格推导出实用的坐标投影公式后，这种投影才得到推广，所以该投影又称高斯-克吕格投影。想象有一椭圆柱面横套在地球椭球体外面，并与某一条子午线（称中央子午线或轴子午线）相切，椭圆柱的中心轴通过椭球体中心，然后用一定的投影方法将中央子午线两侧各一定经差范围内的地区投影到椭圆柱面上，再将此柱面展开即成为投影面。 图2.1 高斯-克吕格投影的几何概念（2）分带投影我国规定按经差和度进行分带投影，为大比例尺测图和工程测量采用带投影。特殊情况下工程测量控制网也用带或任意带。高斯投影带：自子午线起，每隔经差自西向东分带，依次编号1,2,3,。我国中央子午线的经度，由起每隔至，共计12带，带号用n表示，中央子午线的经度用表示，。（带号和经度的关系）4南昌工程学院本科毕业论文图2.2 高斯-克吕格投影分带示意图高斯投影带：是在带的基础上分成的，其中央子午线一部分同带中央子午线重合，一部分则同带分界子午线重合，带号用n/表示，带中央子午线用L表示，关系是：。（3）高斯平面直角坐标系在投影面上，中央子午线和赤道的投影都是直线，并且以中央子午线和赤道的交点O作为坐标原点，以中央子午线的投影为纵坐标轴，以赤道的投影为横坐标轴，这样便形成了高斯平面直角坐标系。在我国坐标均为正，坐标的最大值（在赤道上）约为330KM。为避免出现负的横坐标，可在横坐标上加500KM。此外还应在坐标前面冠以带号，这种坐标称为国家统一坐标。如某点Y=19123456.789m，该点位于19带内，其相对于中央子午线而言的横坐标是：首先去掉带号，再减去500KM，最后得y=-376543.211m。（4）高斯投影的特性与优点高斯-克吕格特性（1）等角投影投影前后的角度相等，但长度和面积有变形；（2）等距投影投影前后的长度相等，但角度和面积有变形；（3）等积投影投影前后的面积相等，但角度和长度有变形。高斯平面投影的特点中央子午线无变形；无角度变形，图形保持相似；离中央子午线越远，变形越大。由此可见，在测量中，如果中央子午线输错了，投影的中央子午线就会编离实地坐标系正5第二章 高斯-克吕格投影（高斯投影）确的中央子午线，变形就越大，最终的结果就使用测量的误差更大。3度带与6度带的划分：1我国采用6度分带和3度分带：12.5万及15万的地形图采用6度分带投影，即经差为6度，从零度子午线开始，自西向东每个经差6度为一投影带，全球共分60个带，用1，2，3，4，5，表示即东经06度为第一带，其中央经线的经度为东经3度，东经612度为第二带，其中央经线的经度为9度。11万的地形图采用3度分带，从东经1.5度的经线开始，每隔3度为一带，用1，2，3，表示，全球共划分120个投影带，即东经1.5 4.5度为第1带，其中央经线经度为东经3度，东经4.57.5度为第2带，其中央经线的经度为东经6度。我省处于东经113度至东经120度之间，共跨第38、39、40，此三个带，其中东经115.5度以西便为第38带，其中央经线为东经114度；东经115.5118.5度为39带，其中央经线为东经117度；东经118.5度以东到山海关为40带，其中央经线为东经120度。地形图上公里网横坐标前2位就是带号，例如：15万地形图上的横坐标为20345486，其中20即为带号，345486为横坐标值。2当地中央经线经度的计算六度带中央经线经度的计算：当地中央经线经度6当地带号3，例如：地形图上的横坐标为20345，其所处的六度带的中央经线经度为：6203117（适用于125万和15万地形图）。三度带中央经线经度的计算：中央经线经度3当地带号（适用于11万地形图）。3、怎样计算当地的中央子午线？ 当地中央子午线取决于当地的直角坐标系统，首先确定您的直角坐标系统是3度带还是6度带投影公式推算:6度带中央子午线计算公式：当地经度/6=N；中央子午线L=6* N （带号）如果没有除尽时，即N有余数时，中央子午线L=6*N - 3 3度带中央子午线计算公式： 当地经度/3=N，中央子午线L=3 *N 我国的经度范围西开始于73东至135，可分成六度带十一个（13号带23号带），各带中央经线依次为（75、81、123、129、135）；三度带二十二 个（24号45号带）。各带中央经线依次为（72、75、132、135）；六度带可用于中小比例尺（如 1：250000）测图，三度带可用于大比例尺（如 1：10000）6南昌工程学院本科毕业论文测图，城建坐标多采用三度带的高斯投影4、如何判断投影坐标是3度带坐标还是6度带坐标如(4231898,21655933)其中21即为带号，同样所定义的东伪偏移值也需要加上带号，如21带的东伪偏移值为21500000米。 假如你的工作区经度在120度至126度范围，则该坐标系为6度带坐标系，该带的中央经度为123度。如(2949320,36353822)其中36即为带号，已知该地点位于贵阳市附近，而从地图上我们看到贵阳大概的经度是东经108度左右，因此可以36*3=108，所以该坐标系为3度带坐标系，该带的中央经度为108度。而不可能为6度带：36*6=2162.2高斯投影坐标正反算公式（1）高斯投影坐标正算公式： B, x,y高斯投影须满足以下三个条件：中央子午线投影后为直线；中央子午线投影后长度不变；投影具有正形性质，即正形投影条件。由第一条件知中央子午线东西两侧的投影必然对称于中央子午线，即(8-10)式中，x为的偶函数，y为的奇函数；，即，如展开为的级数，收敛。 （2-1） 式中是待定系数，它们都是纬度B的函数。由第三个条件知：(2-1)式分别对和q求偏导数并代入上式 上两式两边相等，其必要充分条件是同次幂前的系数应相等，即7第二章 高斯-克吕格投影（高斯投影） (2-2) (2-2)是一种递推公式，只要确定了就可依次确定其余各系数。 由第二条件知:位于中央子午线上的点，投影后的纵坐标x应等于投影前从赤道量至该点的子午线弧长X，即(2-1)式第一式中，当时有： (2-4)顾及(对于中央子午线)得： (2-5) (2-6)依次求得并代入(2-1)式，得到高斯投影正算公式8南昌工程学院本科毕业论文 （2）高斯投影坐标反算公式 x,y B,投影方程： (2-7)满足以下三个条件：x坐标轴投影后为中央子午线是投影的对称轴； x坐标轴投影后长度不变；投影具有正形性质，即正形投影条件。由x求底点纬度(垂足纬度),对应的有底点处的等量纬度，求x,y与的关系式， 由于y和椭球半径相比较小(1/16.37)，可将展开为y的幂级数；又由于是对称投影，q必是y的偶函数，必是y的奇函数。 (2-8)是待定系数，它们都是x的函数.由第三条件知：， ， (2-9)(3-8)式分别对x和y求偏导数并代入上式9第二章 高斯-克吕格投影（高斯投影） 上式相等必要充分条件，是同次幂y前的系数相等， 第二条件，当y=0时，点在中央子午线上，即x=X，对应的点称为底点，其纬度为底点纬度，也就是x=X时的子午线弧长所对应的纬度，设所对应的等量纬度为。也就是在底点展开为y的幂级数。 由(2-8) 依次求得其它各系数(2-10) (2-11) 将代入(2-8)式得 (2-12)10南昌工程学院本科毕业论文 (2-13) 将代入(2-8)式得(2-14)式。(最后表达式)求与的关系。由(8-7)式知： (2-14) (2-15) 按台劳级数在展开(2-16) （2-17）由此可求出各阶导数： (2-18) (2-19) (2-20) 11第二章 高斯-克吕格投影（高斯投影）将式(2-13)1,(2-18),(2-19),(2-20)代入(2-17)并按y幂集合得高斯投影坐标反算公式(2-21), (2-21)归纳由求的基本思想：由点得到底点，将底点f作为过渡，也就是说将坐标原点o移到f点，先求关系式，再将关系式代入关系式得关系式，最后将坐标原点移回到o点,从而求得12南昌工程学院本科毕业论文第三章 c语言程序代码3.1地图投影变换的计算程序流程图 开始输入B,L已知a，b已知a，b由a，b得e,e,e2,N,R由L得L0由L得L1由B得B1由a,e,B1得R由e,B得C由L1,L0,B1得A由e,B1得M由a,e,B1得N由L得FE由FE，KO，N，A,T,C得X由KO，M,N,B1,T,C,得Y 结束 输出X,Y13第三章 C语言程序代码开始输入X,Y已知a，b由a，b得f,e,e2,E由X,Y得L0,X0.Y0由X0，Y0得xval,yval由B得B1由ycal得M由e2,a,M得U由U,e1,e2得Bf由a,Bf得C,T由a,e2,Bf得R由xval,N得D由公式得L1,B1由B1得B由L1得L 输出B,L结束14南昌工程学院本科毕业论文3.2 C语言程序源代码高斯正算：/高斯正算：#includestdio.h#includestdlib.h#includemath.hdouble B,L,X,Y;/B是纬度L是精度；double sort(double,double); main()printf(请输入坐标（经纬度以小数表示）：);printf(纬度B=);/B是纬度scanf(%lf,&B);getchar();printf(经度L=);/L是经度scanf(%lf,&L);getchar();double abeijing=6378245; /北京 54double bbeijing=6356863.0188; /北京 54double axian=6378140; /西安 80double bxian=6356755.2882; /西安 80double awgs=6378137; /WGS 84double bwgs=6356752.3142; /WGS 84sort(abeijing,bbeijing);/北京坐标系的转换 函数调用printf(n 北京坐标系的转换结果：X=%lf Y=%lfn,X,Y);/打印结果sort(axian,bxian);/西安坐标系的转换 函数调用printf(西安坐标系的转换结果：X=%lf Y=%lfn,X,Y);/打印结果 sort(awgs,bwgs);/WGS坐标系的转换 函数调用printf(WGS坐标系的转换结果：X=%lf Y=%lfn,X,Y);/打印结果 double sort(double a0,double b0)/函数，传值进入 double a,b;/接受参数 a,b即长半轴，短半轴a=a0;b=b0;double f=(a-b)/a;/偏率double e,e2;/a,b为参数，e,e2 为第一第二偏心率/偏心率平方e=1-(b/a)*(b/a);e2=(a/b)*(a/b)-1;double dblPI;dblPI=3.1415926/180.0;/ 角度弧度转换参数int ZoneWide;/ 带宽3度带或6度带printf(n 请输入带宽 （3度带或是6度带）ZoneWide=);scanf(%d,&ZoneWide);15第三章 C语言程序代码getchar();int ZoneNumber; /带号double L0;/中央经度if (ZoneWide=3)/3度带ZoneNumber=(int)(L-ZoneWide/2)/ZoneWide+1);else if (6=ZoneWide)/6度带ZoneNumber=(int)(L/ZoneWide+1);else printf(输入错误!将退出!); exit(0);/中央经度（用字母LO表示）L0=(ZoneNumber-1)*ZoneWide+ZoneWide/2;/角度转换成弧度double L1=L*dblPI;double B1=B*dblPI;L0=L0*dblPI;/子午圈曲率半径 doubleR=a*(1-e)/sqrt(1-e*sin(B1)*(1-e*sin(B1)*sin(B1)*(1-e*sin(B1)*sin(B1); double T=tan(B1)*tan(B1); double C=e2*cos(B1)*cos(B1); double A=(L1-L0)*cos(B1); /子午线弧长 doubleM=a*(1-e/4-3*e*e/64-5*e*e*e/256)*B1-(3*e/8+3*e*e/32+45*e*e*e/1024)*sin(2*B1)+(15*e*e/256+45*e*e*e/1024)*sin(4*B1)-(35*e*e*e/3072)*sin(6*B1); /卯酉圈曲率半径 double N=a/sqrt(1-e*sin(B1)*sin(B1); /东纬度偏移 double FE=ZoneNumber*1000000+500000; double k0=1;/Y=FE+k0*N*(A=(1-T+C)*A*A*A/6+(5-18*T*T*T+14*C-58*T*C)*A*A*A*A*A*A/120);X= k0*(M+N*tan(B1)*(A*A/2+(5-T+9*C+4*C*C)*A*A*A*A/24)+(61-58*T*T*T+270*C-330*T*C)*A*A*A*A*A*A/720); /printf (%lf %lf %lf %lf,B,L,a，b）； return X,Y ;16南昌工程学院本科毕业论文高斯反算：/高斯反算：#includestdio.h#includestdlib.h#includemath.hdouble B,L,X,Y;/B是纬度L是精度；double abeijing,bbeijing,axian,bxian,awgs,bwgs;double sort(double ,double );main() printf(请输入直角坐标：);printf(横坐标Y=);scanf(%lf,&Y);getchar();printf(纵坐标X=);scanf(%lf,&X);getchar();double abeijing=6378245; /北京 54double bbeijing=6356863.0188; /北京 54double axian=6378140; /西安 80double bxian=6356755.2882; /西安 80double awgs=6378137; /WGS 84double bwgs=6356752.3142; /WGS 84sort(abeijing,bbeijing);/北京坐标系的转换 函数调用printf(n 直角坐标为54北京坐标系时：B=%lf L=%lfn,B,L);/打印结果sort(axian,bxian);/西安坐标系的转换 函数调用printf(n 直角坐标为80西安坐标系时：B=%lf L=%lfn,B,L);/打印结果 sort(awgs,bwgs);/WGS坐标系的转换 函数调用printf(n 直角坐标为WGS坐标系时：B=%lf L=%lfn,B,L);/打印结果 double sort(double a0,double b0)/函数，传值进入 double a,b;/接受参数 a,b即长半轴，短半轴a=a0;b=b0;double f=(a-b)/a;/偏率double e,ee;/a,b为参数，e,ee 为第一第二偏心率/偏心率平方e=1-(b/a)*(b/a);ee=(a/b)*(a/b)-1;double E=ee/(1-ee);double dblPI=3.1415926/180.0;/ 角度弧度转换参数dblPI=3.1415926/180.0;/ 角度弧度转换参数int ZoneWide;/ 带宽3度带或6度带printf(n 请输入带宽 （3度带或是6度带）ZoneWide=);scanf(%d,&ZoneWide);17第三章 C语言程序代码getchar();int ZoneNumber; /带号if (ZoneWide=3)/3度带ZoneNumber=(int)(Y/500000);else if (6=ZoneWide)/6度带ZoneNumber=(int)(Y/1000000);else printf(输入错误!将退出!); exit(0);/中央经度（用字母LO表示）double L0=(ZoneNumber-1)*ZoneWide+ZoneWide/2;/角度转换成弧度L0=L0*dblPI;double X0;double Y0=ZoneNumber*1000000+500000; X0=0.00;double xval=X-X0;double yval=Y-Y0;/子午线弧长double M=yval; double u=M/(a*(1-ee/4-3*ee*ee/64-5*ee*ee*ee/256); /底点纬度double Bf=u+(3*e/2-27*e*e*e/32)*sin(2*u)+(21*e*e/16-55*e*e*e*e/32)*sin(4*u)+(151*e*e*e/96)*sin(6*u)+(1097*e*e*e*e/512)*sin(8*u); double C=E*cos(Bf)*cos(Bf); double T=tan(Bf)*tan(Bf); /卯酉圈曲率半径 double N=a/sqrt(1-e*sin(Bf)*sin(Bf);/底点所对曲率半径double Nf=a*(1-ee)/sqrt(1-ee*sin(Bf)*sin(Bf)*(1-ee*sin(Bf)*sin(Bf)*(1-ee*sin(Bf)*sin(Bf);double D=yval/N;/计算纬度L,经度B double L1=L0+(D-(1+2*T+C)*D*D*D/6+(5-2*C+28*T-3*C*C+8*ee+24*T*T)*D*D*D*D*D/120)/cos(Bf); double B1=Bf-(N*tan(Bf)/Nf)*(D*D/2-(5+3*T+10*C-4*C*C-9*E)*D*D*D*D/24+(31+90*T+298*C+45*T*T-256*E-3*C*C)*D*D*D*D*D*D/720);/转换为度 L=L1/dblPI; B=B1/dblPI;/printf (%lf %lf %lf %lf,B,L,a，b）; return B,L ;18南昌工程学院本科毕业论文3.3程序运行截图 3.3.1 高斯正算程序运行截图一（1）正算程序代码，进过编译运行后，弹出对话框（2）当分别输入纬度和精度L数据后，选择6度带，回车，出现下图结果（3）同理，按任意键继续，再继续输入纬度B和经度L数据，选择3度带，回车 3.3.2 高斯正算程序运行截图二19第三章 C语言程序代码高斯反算：（4）反算与正算类似，输入完程序代码后，经过编译运行后，弹出对话框 （5）当分别输入横坐标Y和纵坐标X数据之后，选择6度带，回车 （6）按任意键继续，再次输入纬度B和经度L数据，选择3度带，回车 图3.3.3 高斯反算程序运行截图一 图3.3.4 高斯反算程序运行截图二20南昌工程学院本科毕业论文第四章 算例正算算例如下表，与之输入相近数据，验证，程序运行结果正确。计算次序 测站符号D1B29341651122L1062514.86633L01054l12505.58175l”5114.86636l=l”10-80.26161869a2,10-3225411b2,10-75448913c2,10-729708a13220.01614a2l0.59017al+a2l3220.60614b126.916653815b2l1425518b1+ b2l26.918079312c10.493505716c2l77717c1+c2l0.493583421第四章 算例计算次序符号 点号 A1x3275611.1882y180853.8683y=y2*10-103271.0298A1,10-6-1.7619B2,10-6181.6