高斯投影坐标正反算编程

高斯投影坐标正反算编程报告

1. 编程思想

进行高斯投影坐标正反算的编程需要牵涉到大量的公式，为了使程序条理更清楚，各块的数据复用性更强，这里采取了结构化的编程思想。

程序由四大块组成。

GeodesyHomework.cpp文件用于存放main()函数，是整个程序的入口。通过结构化的编程尽力使main()函数变得简单。

MyFunction.h和MyFunction.cpp用于存放计算过程中进行角度弧度换算时所要用到的一些自定的转换函数。

Zhengsuan.h和Zhengsuan.cpp用于存放Zhengsuan类，在Zhengsuan类中声明了高斯投影坐标正算所要用到的所有变量，在类的构造函数中进行成员变量的初始化及正算计算。通过get函数获得相应的正算结果。

Fansuan.h和Fansuan.cpp用于存放Fansuan类，类似于Zhengsuan类，Fansuan类中声明了高斯投影坐标反算所要用到的所有变量，在类的构造函数中进行成员变量的初始化及反算计算。通过get函数获得相应的反算结果。

2. 计算模型

高斯投影正算公式

x?X??NN232244????sinBcosB?l?simBcosB(5?t?9??4?)l2???22

§8.3高斯投影坐标正反算公式

任何一种投影①坐标对应关系是最主要的；②如果是正形投影，除了满足正形投影的条件外（C-R偏微分方程），还有它本身的特殊条件。 8.3.1高斯投影坐标正算公式： B,l x,y

高斯投影必须满足以下三个条件：

①中央子午线投影后为直线；②中央子午线投影后长度不变；③投影具有正形性质，即正形投影条件。

由第一条件知中央子午线东西两侧的投影必然对称于中央子午线，即(8-10)式中，x为l的偶函数，y为l的奇函数；l 30

30 ，即l /

1/20，如展开为l的级数，收敛。

x m24

6

0 m2l m4l m6l y m3

5

（8-33）

1l m3l m5l

式中

m0,m1, 是待定系数，它们都是纬度B的函数。

由第三个条件知：

xx y

q

y l, l q

(8-33)式分别对l和q求偏导数并代入上式

m1 3m24

dm0dm22

dm44

3l 5m5l

dq

dql dql 2ml3

6m5

6l

dm1 (8-34)

dq

l

dm355

2l 4m4dq

l3

dmdq

l

上两式两边相等，其必要充分条件是同次幂l前的系数应相等，即

mdm01

dq

m 1dm1

22

dq

(8-35)

m1dm2

3 3

高斯投影坐标正反算

一、基本思想：

高斯投影正算公式就是由大地坐标（L，B）求解高斯平面坐标（x，y），而高斯投影反算公式则是由高斯平面坐标（x，y）求解大地坐标（L，B）。

二、计算模型：

基本椭球参数： 椭球长半轴a 椭球扁率f

椭球短半轴：b?a(1?f)

a2?b2椭球第一偏心率 ：e?

aa2?b2椭球第二偏心率 ：e?? b高斯投影正算公式：此公式换算的精度为0.001m

x?X??NN2??sinBcosB?l?simBcos3B(5?t2?9?2?4?4)l??4242???24???

N5246??sinBcosB(61?58t?t)l720???6y?NN3223??cosB?l???cosB(1?t??)l???6???3N5242225???cosB(5?18t?t?14??58?t)l120???5

其中：角度都为弧度

B为点的纬度，l???L?L0，L为点的经度，L0为中央子午线经度；

N为子午圈曲率半径，N?a(1?esinB)； t?tanB；

22?12?2?e?2cos2B

????180??3600

其中X为子午线弧长：

1616??X?a0B?sinBcosB?(a2?a4?a6)?(2a4?a6)sin

大地测量学编程实习报告

姓名：鲁尼 学号：10 班级：曼联

编程思想:

这个投影是由德国数学家、物理学家、天文学家高斯于19 世纪20 年代拟定，后经德国大地测量学家克吕格于1912 年对投影公式加以补充，故称为高斯-克吕格投影。

即等角横切椭圆柱投影。假想用一个圆柱横切于地球椭球体的某一经线上，这条与圆柱面相切的经线，称中央经线。以中央经线为投影的对称轴，将东西各3°或1°30′的两条子午线所夹经差6°或3°的带状地区按数学法则、投影法则投影到圆柱面上，再展开成平面，即高斯-克吕格投影，简称高斯投影。这个狭长的带状的经纬线网叫做高斯-克吕格投影带。

高斯投影正算公式就是由大地坐标（L，B）求解高斯平面坐标（x，y），而高斯投影反算公式则是由高斯平面坐标（x，y）求解大地坐标（L，B）。现行的高斯投影用表都是采用克拉索夫斯基椭球参数，这次编程计算就是采用这种椭球参数，并采用实用公式按6度分带投影。编程环境是在VC下，采用C++语言编写。程序主要分为两部分，第一部分是高斯正反算函数，第二部分是主函数。 高斯正反算函数，参考书上175页的电算公式，正算时先将度数换算成秒，再定带号n，求中央经线l0，经度差l''，然后根据实用公式计算高斯

高斯投影坐标正反算

一、相关概念

大地坐标系由大地基准面和地图投影确定，由地图投影到特定椭圆柱面后在南北两极剪开展开而成，是对地球表面的逼近，各国或地区有各自的大地基准面，我国目前主要采用的基准面为：

1.WGS84基准面，为GPS基准面，17届国际大地测量协会上推荐，椭圆柱长半轴a=6378137m，短半轴b=6356752.3142451m；

2.西安80坐标系，1975年国际大地测量协会上推荐，椭圆柱长半轴a=6378140m，短半轴b=6356755.2881575m;

3.北京54坐标系，参照前苏联克拉索夫斯基椭球体建立,椭圆柱长半轴a=6378245m, 短半轴b=6356863.018773m;

通常所说的高斯投影有三种，即投影后：

a) 角度不变（正角投影），投影后经线和纬线仍然垂直； b) 长度不变； c) 面积不变；

大地坐标一般采用高斯正角投影，即在地球球心放一点光源，地图投影到过与中央经线相切的椭圆柱面上而成；可分带投影，按中央经线经度值分带，有每6度一带或每3度一带两种(起始带中央经线经度为均为3度，即：6度带1带位置0-6度，3度带1带位置1.5-4.5 度)，即所谓的高斯-克吕格投影。

图表 11高斯投影和分带

§8.3高斯投影坐标正反算公式

任何一种投影①坐标对应关系是最主要的；②如果是正形投影，除了满足正形投影的条件外（C-R偏微分方程），还有它本身的特殊条件。 8.3.1高斯投影坐标正算公式： B,l x,y

高斯投影必须满足以下三个条件：

①中央子午线投影后为直线；②中央子午线投影后长度不变；③投影具有正形性质，即正形投影条件。

由第一条件知中央子午线东西两侧的投影必然对称于中央子午线，即(8-10)式中，x为l的偶函数，y为l的奇函数；l 3030 ，即l / 1/20，如展开为l的级数，收敛。

x m0 m2l2 m4l4 m6l6 y m1l m3l3 m5l5

（8-33）

式中m0,m1, 是待定系数，它们都是纬度B的函数。 由第三个条件知：

x y x y

, q l l q

(8-33)式分别对l和q求偏导数并代入上式

dm0dm22dm44

m1 3m3l 5m5l l l

dqdqdq

(8-34) dm33dm55dm135

2m2l 4m4l 6m6l l l l

dqdqdq

2

4

上两式两边相等，其必要充分条件是同次幂l前的系数应相等，即

61

dm0

m1

dq

1dm1

m2

2dq

1dm2

一、坐标正算与坐标反算 1、坐标正算 已知点的坐标、

边的方位角、

两点间的水平距离，计算待

定点的坐标，称为坐标正算。 如图6-6 所示，点的坐标可由下式计算：

式中

、

为两导线点坐标之差，称为坐标增量，即：

=1000、

=1000、方位角

【例题6-1】已知点A坐标，=35°17＇36.5\，

2、坐标反算 已知

两点的坐标，计算

两点水平距离

=200.416，计算点的坐标？ 35o17＇36.5\ 35o17＇36.5\

两点的水平距离与坐标方位角，

称为坐标反算。如图6-6 可知，由下式计算水平距离与坐标方位角。

（6-3）

（6-4）

式中反正切函数的值域是-90°～+90°，而坐标方位角为0°～360°，因此坐标方位角的值，可根据正切角值换算为坐标方位角。

、

的正负号所在象限，将反

【例题6-2】=3712227.860、

、水平距离

=3712232.528、=523620.436、

=523611.598，计算坐标方位角计算坐标方位角。

=62°0

摘要

本设计主要阐述了高斯投影分带以及高斯投影坐标正、反算的推导公式，从而根据公式来编写基于VB语言基础上的换带及坐标转换程序。作者系统介绍了测量中经常使用的坐标系以及地图投影的概念和高斯投影的具体含义，叙述了换带和临带计算的原因以及它们在运算时的原理、过程，详细叙述了在VB语言中实现的原理基础以及代码的编写设计。

在设计中根据高斯的正反算公式写出了基于VB语言的程序设计，其程序设计任务完成了由地理坐标向54平面坐标系和80平面坐标系转换的功能，以及由54坐标系和80坐标系向地理坐标系转换的功能，同时也有同一平面坐标系不同投影带之间的换带计算和同一平面坐标系相同投影带临带计算等相互转换的功能。

关键词：高斯投影、坐标正反算、换带计算、临带换算、程序设计

5 程序设计

5.1界面设计

本程序要实现的功能是根据所选择的椭球参数和指定的分带情况，将已知地理坐标或高斯投影坐标经正算和反算求得相应的高斯坐标和地理坐标，以及相应的换带计算和临带计算。因此需要用一个框架控件来组织椭球参数、两个框架分别组织分带选择和换算方式选择，两个框架组织地理坐标和高斯坐标，三个命令按钮分别执行投影计算、换带和临带计算。程

1

序设计界面如图5-1[9]

图5-1 高斯

线路中边桩任意里程坐标正、反算改编程序（CASIO fx-4800改CASIO fx-5800）

1. 正算子程序(SUB1)

0.1739274226→Z[3]：0.3260725774→Z[4]：0.0694318442→Z[5]：0.3300094782→Z[6]：1-Z[6]→Z[7]：1-Z[5]→Z[8]：A+Z[10] (Z[3]cos(F +57.2958 O Z[5] Z[10] (1/C+ Z[5] Z[10] Z[9]))+Z[4]cos(F +57.2958 O Z[6] Z[10] (1/C+ Z[6] Z[10] Z[9]))+Z[4]cos(F +57.2958 O Z[7] Z[10] (1/C+ Z[7] Z[10] Z[9]))+ Z[3]cos(F +57.2958 O Z[8] Z[10] (1/C+ Z[8] Z[10] Z[9]))) →Z[11]：B+ Z[10] (Z[3]sin(F +57.2958 O Z[5] Z[10] (1/C+ Z[5] Z[10] Z[9]))+Z[4]sin(F +57.2958 O Z[6] Z[10] (1/C+ Z[6] Z[10] Z[9]))+Z[4

利用常用的办公软件EXCL表格自身的运算函数,通过一系列的编写程序达到坐标的正反算。一切都是自动化输出结果,只要简单的输入基本的测量数据即可。它不仅是一份文件,更是一款实用的软件,特别适合于工程技术人员大量的坐标正反换算,计算坐标和放样数据。

点号 D1 D2 D3 D4 D5 D6 D7 D8 D9 D10 D11 D12 D13 D14 D15 D16 D17 D18 D19 X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X

坐标值(m)

利用常用的办公软件EXCL表格自身的运算函数,通过一系列的编写程序达到坐标的正反算。一切都是自动化输出结果,只要简单的输入基本的测量数据即可。它不仅是一份文件,更是一款实用的软件,特别适合于工程技术人员大量的坐标正反换算,计算坐标和放样数据。

点号 D19 D20 D21 D22 D23 D24 D25 D26 D27 D28 D29 D30 D31 D32 D33 D34

坐标值(m) Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y X Y

利用常用