6 高斯投影及其计算

合集下载

大地测量学第六章高斯投影及其计算

大地测量学第六章高斯投影及其计算
应用大地测量学
第六章 高斯投影 及其计算
中国矿业大学环境与测绘学院
第六章 高斯投影及其计算概述
1、椭球面上计算复杂; 2、椭球面上表示点位的经度、纬度大地线长、大地
方位角等对大比例尺测图不适应; 3、为了测绘地形图和计算的方便,需通过地图投影
的方法将椭球面上的元素化算到平面上; 4、本章主要介绍正形投影的特性以及高斯投影建立
应用大地测量学
§6.2.2 高斯投影的长度比和长度变形
1、用大地坐标表示的高斯投影长度比m
式中:
2、用平面坐标表示的高斯投影长度比m
m

1

y2 2R 2
y4 24R4
式中y为投影点的横坐标,R为该点处椭球平均曲率半径。
应用大地测量学
§6.2.2 高斯投影的长度比和长度变形
3、长度变形m-1与横坐标y的关系
5 5′
应用大地测量学
§6.3 高斯投影坐标计算
高斯投影坐标正算——由(B,L)求(x,y) 高斯投影坐标反算——由(x,y)求(B,L)
应用大地测量学
§6.3.1 高斯投影坐标正算公式
(6-26)
式中,X为由赤道至纬度B的子午线弧长, 为计算点P点与中央子午线
的经差。N为卯酉圈曲率半径,t=tanB, η=e′cosB。 L-L0若以度为单位,则ρ=57.295779513; L-L0若以分为单位,则ρ=3437.7467708; L-L0若以秒为单位,则ρ=206264.80625。
平面直角坐标系的方法、观测元素的化算、高斯 投影坐标计算。
第六章 高斯投影及其计算
第一节 地图投影概念和正形投影性质 第二节 高斯投影与国家平面直角坐标系(基础) 第三节 高斯投影坐标计算(重点) 第四节 椭球面上的方向和长度归算至高斯投影平面

高斯投影6度和3度分带计算公式

高斯投影6度和3度分带计算公式

高斯投影6度和3度分带计算公式高斯投影6度和3度分带计算公式什么是高斯投影6度和3度分带?•高斯投影是一种常用于大地测量和地图制图的投影方法。

根据地球的形状和表面特征,我们将地球划分成了若干个分带,每个分带的宽度为6度或3度。

•6度和3度分带指的是每个分带的经度跨度。

例如,6度分带就是每个分带的中央经线与相邻分带的中央经线之间跨越6度。

高斯投影6度和3度分带计算公式6度分带投影计算公式1.计算投影平面与地球经度的差值:L=λ−L02.计算弧长元素:N=a/√1−e2⋅sin2φ3.计算卯酉圈曲率半径:M=N⋅(1−e2)=a⋅(1−e2)/(1−e2⋅sin2φ)4.计算子午线弧长:A=(1+3e2/4+45e4/64+175e6/256+11025e8/16384)⋅N5.计算坐标系原点到点的子午线弧长:S=A−A06.计算纬度差:t=tanφ7.计算坐标Y轴偏移量:y=x⋅cosφ8.计算坐标X、Y(单位:m):X=S−N⋅tanφ2⋅L2−N⋅tanφ24⋅(5−t2+9C2+4C4)⋅L4−N⋅tanφ720⋅(61−58t2+t4−270C2+330C4)⋅L6Y=N⋅L⋅cosφ1+N⋅L3⋅cosφ6⋅(1−t2+C2)+N⋅L5⋅cosφ120⋅(5−18t2+t4+14C2−58C4)3度分带投影计算公式1.计算投影平面与地球经度的差值:L=λ−L02.计算弧长元素:N=a/√1−e2⋅sin2φ3.计算卯酉圈曲率半径:M=N⋅(1−e2)=a⋅(1−e2)/(1−e2⋅sin2φ)4.计算子午线弧长:A=(1+3e2/4+45e4/64+175e6/256+11025e8/16384)⋅N5.计算坐标系原点到点的子午线弧长:S=A−A06.计算纬度差:t=tanφ7.计算坐标Y轴偏移量:y=x⋅cosφ8.计算坐标X、Y(单位:m):X=S−N⋅tanφ2⋅L2+N⋅tanφ24⋅(5+t2+9C2+4C4)⋅L4−N⋅tanφ720⋅(61+90t2+45t4+46C2−252C4−90C6)⋅L6Y=N⋅L⋅cosφ1+N⋅L3⋅cosφ6⋅(1+2t2+C2)+N⋅L5⋅cosφ120⋅(5+28t2+24t4+6C2+8C4)示例解释假设我们需要计算某个点在高斯投影6度分带中的投影坐标。

高斯投影及换带计算

高斯投影及换带计算

测绘学院《大地测量学基础》课件 语言优教资源PPT
22
[知识点及学习要求] 1.高斯投影的基本概念; 2.正形投影的一般条件;
3.高斯平面直角坐标与大地坐标的相互转换
—高斯投影的正算与反算 4.椭球面上观测成果归化到高斯平面上的计算; 5.高斯投影的邻带换算; 6.工程测量投影面与投影带的选择。
[难点]在对本章的学习中,首先要理解和掌握高斯投影的
x F1(L, B) y F2 (L, B)
椭球面是一个凸起的、不可展平的曲面,若将这个曲面上 的元素(比如一段距离、一个角度、一个图形)投影到平 面上,就会和原来的距离、角度、图形呈现差异,这一差 异称作投影的变形
测绘学院《大地测量学基础》课件 语言优教资源PPT
44
长度比:
投影面上的边长与原面上的相应长度之比,称为长度比。
• 3)中国大比例尺地图的投影:多面体投影(北洋军阀时 期)、等角割圆锥投影(兰勃特投影)(解放前)、高斯克吕格投影(解放以后)。
测绘学院《大地测量学基础》课件 语言优教资源PPT
99
4、常用的几种地图投影
从世界范围看,各国大中比例尺地形图所使用 的投影很不统一,据不完全统计有十几种之多,最 常用的有横轴等角椭圆柱投影等。中华人民共和国 成立后,我国大中比例尺地形图一律规定采用以克 拉索夫斯基椭球体元素计算的高斯-克吕格投影。我 国新编1:100万地形图,采用的则是边纬与中纬变 形绝对值相等的正轴等角圆锥投影。
午线为对称轴。投影后有长
度变形。 (3) 赤道线投影后为直线,
但有长度变形。
中央子午线
测绘学院《大地测量学基础》课件 语言优教资源PPT
15 15
x
(4) 除赤道外的其余纬线,投

高斯投影坐标计算

高斯投影坐标计算

B
d B dq
2

dX dq dq

c
(
cos B dV V dB
2
dB dq

sin B dB V dq
2
)
2
d B dq
2
cos B c ( tan B V
2 2
3
V
sin B cos B
)
N sin B cos B
同理得
d X dq
3

N cos B ( 1
3
3


2

0
l

L

L
0

高斯投影坐标正算的函数式:
x y
l 是以弧度为单位的经度差。
F B , l F B , l
1 2

一 高斯投影坐标正算公式计算

如图,椭球面上一点投影 到平面后为d点,椭球面上 该点的平行圈(B或q为一 常数)与中央子午线的交 点为e点,若将上式中的展 开点z0设为e处,则很据高 斯投影条件,中央子午线 的长度比m=1,且纵坐标x 等于从赤道起到该平行圈 间的子午线弧长X。此时 可以写出下列方程:
4 2
二、高斯投影坐标反算公式

最后得到坐标反算的公式为:
B B
f
2M
f
t
f
y N
f
2

t 24 M
2 f
f
f
f
N
4 f
3 f
5 3 t
6
2 f

2 f
9 f t
2
2 f
y
4

t

6高斯投影及其计算

6高斯投影及其计算
第一节 地图投影概念和正形投影性质
应用大地测量学
第一节 地图投影概念和正形投影性质
应用大地测量学
第一节 地图投影概念和正形投影性质
应用大地测量学
(二)投影变形 角度变形、长度变形和面积变形三种。 (三)投影长度比与变形指标 投影长度比——投影面上无限小线段 ds与椭球面上该线段实际长度 dS之比,以m表示,m=ds/dS。长度变形—— v= m-1 变形指标:主方向上投影长度比a和b叫变形指标。 若a=b,则为等角投影,既投影后长度比不随方向而变化。 若ab=1,则为等面积投影。 椭球面上微分圆: 投影平面上对应为微分椭圆:
第一节 地图投影概念和正形投影性质
应用大地测量学
二、正形投影特性 1、任一点上,投影长度比m为一常数,不随方向而变,仅与点位置有关。 2、投影后角度不变形。又叫保角映射或叫正形投影。条件是在微小范围内成立。
第一节 地图投影概念和正形投影性质
应用大地测量学
三、正形投影的一般条件 正形投影必要和充分的条件是满足柯西—黎曼方程:
y/(km)
10
20
30
40
50
100
150
200
250
300
长度变形m-1
1/810000
1/202000
1/90000
1/50000
1/32000
1/8000
1/3500
1/2000
1/1300
1/900
第二节 高斯投影与国家平面直角坐标系
应用大地测量学
三、高斯投影的分带 为限制长度投影变形,投影分带有6度分带和3度分带两种方法。
应用大地测量学
三、距离改正计算 距离改正——椭球面上大地线长S改换为平面上投影曲线两端点间的弦长D,要加距离改正△S。

高斯投影及其计算资料

高斯投影及其计算资料

由图可知
tan 90 A P2P3 MdB N cos Bdq dq P1P3 rdl N cos Bdl dl

dl tan Adq
于是
等量纬度
m2
dx2 dy2
r2[(dq)2 dl
2 ]
应用大地测量学
第一节 地图投影概念和正形投影性质
于是有
x F1B, L 相当于 x xq,l
y
F2
B,
L
y yq,l
dx
x q
dq
x l
dl
dy
y
dq
y
dl
q l
代入
ds2 dx2 dy2
ds2
x q
2
dq 2
2
应用大地测量学
第一节 地图投影概念和正形投影性质
应用大地测量学
第一节 地图投影概念和正形投影性质
应用大地测量学
第一节 地图投影概念和正形投影性质
应用大地测量学
第一节 地图投影概念和正形投影性质
应用大地测量学
第一节 地图投影概念和正形投影性质
应用大地测量学
第一节 地图投影概念和正形投影性质
x q
x l
dq
dl
x l
2
dl
2
y q
2
dq
2
2
y q
y l
dq
dl
y l
2
dl
2
应用大地测量学
第一节 地图投影概念和正形投影性质
ds2
x q
2
y q
2
dq
2
2
x q
x l
y q
y l

高斯投影坐标计算

高斯投影坐标计算

x m0 m 2 l 2 m 4 l 4 y m1l m3 l 3 m5 l 5
式中m0 , m1 , 是待定系数,它们都是纬度B的函数
2) 由第三个条件即正形投影条件可知
y x x y 和 l q l q
分别对l 和q求偏导数并代入上式得
2、高斯投影坐标反算公式
已知高斯平面坐标(x,y),求椭球面上的大地坐标(B,L)的 问题称高斯投影坐标反算。 B 1 ( x, y) 函数式:
l 2 ( x, y)
同正算一样,对投影函数提出三个条件 (1) x (2) x (3) 正形投影条件。
1) 由第一个条件(x 坐标轴投影成中央子午线,是投 影的对称轴)可知
Bf为x值对应的底点纬度, tf ηf Mf Nf 均为底点纬度 的函数。
当l<3.5°时,
上式换算精度达0.0001″
高斯投影反算公式的几何解释
B B f ( n2 y 2 n4 y 4 = Bf高斯投影坐标正算的数值公式 将75国际椭球参数代入前面推导的高斯计算公式, 经过一些简单变化,可得高斯投影正算公式。 高斯投影正算公式:
B 2 2 2 x 6367452 .1328 (a0 (0.5 (a4 a6l )l )l N ) cos B sin B y (1 (a3 a5l 2 )l 2 )lN cos B
实用公式的系数
N 6399596 .652 [21565 .045 (108.996 0.603cos2 B) cos2 B] cos2 B 2 2 2 a 32144 . 5189 [ 135 . 3646 ( 0 . 7034 0 . 0041 cos B ) cos B ] cos B 0 cos2 B) cos2 B 0.04167 a4 (0.25 0.00253 2 2 a ( 0 . 167 cos B 0 . 083 ) cos B 6 0.001123 cos2 B) cos2 B 0.1666667 a3 (0.3333333 a 0.00878 (0.1702 0.20382cos2 B) cos2 B 5

高斯投影及换带计算

高斯投影及换带计算

测绘学院《大地测量学基础》课件
10
6.2 高斯投影概述(重点)
1、控制测量对地图投影的要求
1)等角投影(又称正形投影)
2)长度和面积变形不大,并能用简单公式计算由变形而引起 的改正数。
3)能很方便地按分带进行,并能按高精度的、简单的、同样 的计算公式和用表把各带联成整体 。
测绘学院《大地测量学基础》课件
8
• 3、中国各种地图投影:
1)中国全国地图投影:斜轴等面积方位投影、斜轴等角方 位投影、伪方位投影、正轴等面积割圆锥投影、正轴等角割 圆锥投影。
• 2)中国分省(区)地图的投影:正轴等角割圆锥投影、正 轴等面积割圆锥投影、正轴等角圆柱投影、高斯-克吕格投 影(宽带)。
• 3)中国大比例尺地图的投影:多面体投影(北洋军阀时 期)、等角割圆锥投影(兰勃特投影)(解放前)、高斯克吕格投影(解放以后)。
x F1(L, B) y F2 (L, B)
椭球面是一个凸起的、不可展平的曲面,若将这个曲面上 的元素(比如一段距离、一个角度、一个图形)投影到平 面上,就会和原来的距离、角度、图形呈现差异,这一差 异称作投影的变形
测绘学院《大地测量学基础》课件
4
长度比:
投影面上的边长与原面上的相应长度之比,称为长度比。
(1)该点位于6˚ 带的第几带?
(第19带)
(2)该带中央子午线经度是多少?
(L。=6º×19-3º=111˚)
(3)该点在中央子午线的哪一侧?
(先去掉带号,原来横坐标y=367622.380—500000=-132377.620m,在西侧)
(4)该点距中央子午线和赤道的距离为多少?
(距中央子午线132377.620m,距赤道3102467.280m)

高斯投影及换带计算

高斯投影及换带计算

测绘学院《大地测量学基础》课件
24
高斯平面直角坐标系与数学上的笛卡尔平面直角 坐标系的异同点 :
不同点: 1、 x,y轴互异。 2、 坐标象限不同。 3、表示直线方向的方位角
定义不同。 相同点:
数学计算公式相同。
测绘学院《大地测量学基础》课件
Ⅳx
o

α Ⅰp
D
y

x=Dcosα
y=Dsinα
高斯平面直角坐标系
y3
6N
3 f
cos
Bf
1
2t
2 f
2 f
y5
120N
5 f
cos
Bf
5
28t
2 f
24t
4 f
6
2 f
8
2 f
t
2 f
测绘学院《大地测量学基础》课件
30
3、高斯投影坐 标正反算公式的
几何解释 :
①当B=0时x=X=0,y则随l的变化而变化,这就是说,赤道投影为一直线且 为y轴。当l=0时,则y=0,x=X,这就是说,中央子午线投影亦为直线,且为x轴, 其长度与中央子午线长度相等。两轴的交点为坐标原点。
B B f
tf 2M f N f
y2
tf
24M
f
N
3 f
5
3t
2 f
2 f
9
2 f
t
2 f
y4
过所求点P作中央子午线的垂线,
tf
720M
f
N
5 f
y
61

90t
2 f
45t
4 f
y6
该垂线与中央子午线的交点的纬 度,称垂足纬度。其值由子午线 弧长计算公式反算求得。

高斯投影及其计算

高斯投影及其计算

x y

FF12((BB,,ll))
f (z) f (z0)
f
(z0 )( z
x iy f (q il)
z0 )
f
(z0 ) (z 2!
z0 )2

设 z q il,z0 q ,在高斯投影分带中,经差 l 一般仅为 1.5

3 ,

f (q)
l2 2
d 2 f (q) dq 2

l4 24
d 4 f (q) dq 4

l 6 d 6 f (q)
720
dq6

il
df (q) l3
dq 6
d 3 f (q) l5 d 5 f (q)
dq3
120
dq5

l2 d2X l4 d4X l6 d6X
i
y q

f l

(x iy) l
x i y l l

又:
f f (q il) f q (q il) q (q il) f f (q il) i f

l (q il) l
(三)投影长度比与变形指标
投影长度比——投影面上无限小线段 ds与椭球面上该线段实际 长度 dS之比,以m表示,m=ds/dS。长度变形—— v= m-1
变形指标:主方向上投影长度比a和b叫变形指标。 若a=b,则为等角投影,既投影后长度比不随方向而变化。 若ab=1,则为等面积投影。 椭球面上微分圆: 投影平面上对应为微分椭圆:
大地测量学基础
第一节 地图投影概念和正形投影性质
(四)地图投影的分类

高斯投影坐标计算

高斯投影坐标计算
高斯投影坐标计算
本节要点提要
1、高斯投影坐标正算公式 2、高斯投影坐标反算公式 3、高斯投影坐标正算的数值公式 4、高斯投影坐标反算的迭代计算公式
地图投影的分类
• 按投影变形性质分类: 等角投影 等距投影 等积投影
a=b
• 按投影面分类 : 圆锥面 正轴投影 切投影
a=1 or b=1
圆柱(椭圆柱) 面 横轴投影 割投影
(1)中央子午线投影后为直线; (2)中央子午线投影后长度不变; (3)投影具有正形性质,即正形投影 条件。
高斯投影坐标正算
l =3/ρ=0.052
1) 由第一个条件(中央子午线投影后为直线) 可知,由于地球椭球体是一个旋转椭球体,即 中央子午线东西两侧的投影必然对称于中央子 午线。 x 为 l 的偶函数,而y 则为 l 的奇函数。
由恒等式两边对应系数相等,建立求解待定系数的递推公式
d m d m d m 1 1 0 1 m m m = 2 1 2 3 d q 2 d q 3 d q
m0=?
3) 由第二条件(中央子午线投影后长度不变)可 知,位于中央子午线上的点,投影后的纵坐标 x 应 该等于投影前从赤道量至该点的子午弧长。
Байду номын сангаас
a· b=1
平面投影 斜轴投影
• 按投影的中心轴线: • 按椭球面与投影面的切割情况分:
高斯投影特性(三个): – 中央子午线投影后为一直线,且长度不变; 其它经线为凹向中央子午线的曲线,且长 度改变。 – 投影后,赤道为一直线,但长度改变,其 它纬线呈凸向赤道的曲线。 – 投影后,中央子午线与赤道线正交,经线 与纬度也互相垂直,即高斯投影为等角投 影。
将各系数代入,略去高次项,得高斯投影 坐标正算公式精度为0.001m

高斯6度带计算公式

高斯6度带计算公式

高斯6度带计算公式高斯6度带是指地球表面上横跨6度经度的区域,在地图制图中经常用到。

计算高斯6度带的方法是基于高斯投影法的,需要利用椭球体参数、中央子午线经度和目标点经纬度等多个参数进行计算。

具体来说,计算高斯6度带的公式为:1. 计算目标点所在的带号带号 = int((L+180)/6)+1其中,L为目标点所在经度,int为向下取整函数。

2. 计算目标点距离中央子午线的经度差deltaL = L - (带号-1)*6 - 33. 计算椭球体参数首先需要确定使用的椭球体参数。

常用的有WGS84和北京54两种参数。

根据不同的参数,可以得到对应的长半轴a和第一偏心率e1。

4. 计算子午线弧长m利用椭球体参数和目标点所在纬度,可以计算出子午线弧长m,即:m = a(1-e1^2)/(1-e1^2*sin^2B)^(1/2)其中,B为目标点的纬度。

5. 计算底点纬度B0底点纬度B0是指目标点所在带的纬度值,计算公式为:B0 = 0 (当带号为1时)B0 = 2 (当带号为2时)B0 = 4 (当带号为3时)...B0 = (带号-1)*6-3 (当带号大于等于4时)6. 计算底点子午线弧长m0和底点到目标点的纬度差deltaB利用底点纬度B0和椭球体参数,可以计算出底点子午线弧长m0和底点到目标点的纬度差deltaB,即:m0 = a(1-e1^2)/(1-e1^2*sin^2B0)^(1/2)deltaB = (B-B0)*3600其中,3600是角度和弧度之间的换算系数。

7. 计算高斯投影坐标X和Y最终的高斯投影坐标X和Y可以通过以下公式计算得到:X = k*m*(deltaL*3600)^2/2 +k*m*(1-tan^2B)*(deltaL*3600)^4/24*(5-tan^2B+9e'^2cos^2B+4e' ^4cos^4B)Y = k*m*deltaL*3600 +k*m*tanB*(deltaL*3600)^3/6*(1+2tan^2B+e'^2cos^2B) 其中,k为比例尺因子,e'为第二偏心率,由e'^2 =(a^2-b^2)/b^2计算得到,b为短半轴。

6 高斯投影及其计算

6 高斯投影及其计算
2 2 x y E q q x x y y F q l q l 2 2 x y G l l
其推证步骤为:
1、从长度比表达式出发 系式; ,求出m2与dx2,dy2和dB2,dl2关
2、引入等量纬度q,将x、y表为q、l的函数; 3、对 x=f1(q,l),y=f2(q,l)取全微分,引入符号E、F、G; 4、根据长度比m与方向A无关,得F=0,E=G; 5、由E=G、F=0得主要条件。
应用大地测量学
应用大地测量学
第一节 地图投影概念和正形投影性质
二、正形投影特性
两个基本要求: 1、任一点上,投影长度比m 为一常数,不随方向而变, 仅与点位臵有关。 2、投影后角度不变形。又叫保角映射。条件是在微小范 围内成立。 所以,正形投影的特性是:投影长度比m仅与点的位臵有 关,而与方向无关。
应用大地测量学
第一节 地图投影概念和正形投影性质
四、正形投影一般公式
根据复变函数理论,下列复变函数满足柯西(Cauchy)—黎曼 (Riemann)条件,式中,f代表任意解析函数。
x iy f (q il )
证明:
f x iy x y i q q q q f x iy x y i l l l l
x y x y q q l l
2 2 2 2
x y q l
x y l q
从椭球面到平面投影的柯西-黎曼条件
应用大地测量学
第一节 地图投影概念和正形投影性质
2 2

(完整版)高斯投影正反算

(完整版)高斯投影正反算

高斯投影正反算学院:资源与环境工程工程学院专业:测绘工程学号:X51414012姓名:孙超一、高斯投影概述想象有一个椭圆柱面横套在地球椭球体外面,并与某一条子午线相切,椭圆柱的中心轴通过椭球体的中心,然后用一定投影方法,将中央子午线两侧各一定经差范围内的地区投影到椭圆柱面上,再将此柱面展开即成为投影面。

高斯投影由于是正形投影,故保证了投影的角度不变性,图形的相似性以及在某点各方向上长度比的同一性。

由于采用了同样法则的分带投影,这即限制了长度变形,又保证了在不同投影带中采用相同的简便公式和数表进行变形引起的各项改正的计算,并且带与带间的互相换算也能用相同的公式和方法进行。

高斯投影的这些优点必将使它得到广泛的推广和具有国际意义。

二、高斯投影坐标正算公式1.高斯投影必须满足以下三个条件1)中央子午线投影后为直线2)中央子午线投影后长度不变3)投影具有正形性质,即正形投影条件2.高斯正算公式推导1)由第一个条件可知,由于地球椭球体是一个旋转椭球体,所以高斯投影必然有这样一个性质,即中央子午线东西两侧的投影必然对称于中央子午线。

2)由于高斯投影是换带投影,在每带内经差l是不大的,lρ是一个微小量,所以可以将 X=X (l,q ),Y=Y (l ,q )展开为经差为l 的幂级数,它可写成如下的形式X=m 0+m 2l 2+m 4l 4+…Y=m 1l+m 3l 2+m 5l 5+…式中m 0,m1,m2,…是待定系数,他们都是纬度B 的函数。

3)由第三个条件:∂y ∂l =∂x ∂q 和∂x ∂l =-∂y ∂q ,将上式分别对l 和q 求偏导2340123423401234...........x m m l m l m l m l y n n l n l n l n l =+++++=+++++可得到下式0312123403121234111,,,, 234111,,,,234dm dm dm dm n n n n dq dq dq dq dn dn dn dn m m m m dq dq dq dq ⎧====⎪⎪⎨⎪=-=-=-=-⎪⎩L L 经过计算可以得出232244524632235242225sin cos sin cos (594)224sin cos (6158)720cos cos (1) 6cos (5181458)120N N x X B B l B B t l N B B t t l N y N B l B t l N B t t t l ηηηηη=+⋅+-+++-+=⋅+-++-++-三、高斯投影坐标反算公式推导1.思路:级数展开,应用高斯投影三个条件,待定系数法求解。

常用地图投影之高斯

常用地图投影之高斯

常用地图投影之高斯高斯-克吕格投影(简称高斯投影)的概念从几何概念上分析,它是一种横轴等角切圆柱投影。

我们把地球看成是地球椭球体,假想用一个椭圆筒横套在其上,使筒与地球椭球的某一经线(称为中央经线)相切,椭圆筒的中心轴位于赤道上,按等角条件将地球表面投影到椭圆筒上,然后将椭圆筒展开成平面。

高斯投影的基本条件•中央经线和赤道投影成垂直相交的直线•投影后没有角度变形,那么经纬线投影后仍正交•中央经线没有长度变形投影的变形分析其长度比的基本公式为:μ=1+1/2cos2ϕ(1+η2)λ2+1/6cos4ϕ(2−tg2ϕ)λ4−1/8cos4ϕλ4μ=1+1/2cos2ϕ(1+η2)λ2+1/6cos4ϕ(2−tg2ϕ)λ4−1/8cos4ϕλ4长度变形的规律是:•中央经线(λλ=0)上没有长度变形,即λλ=0,μμ=1•同一条纬线上,离中央经线越远变形越大,即λλ增大,μμ也增大•在同一经线上,纬度越低,变形越大,即ϕϕ越小,μμ越大投影分带的规定在1:2.5万到1:50万时,6060分带在大于1:1万地形图中:3030分带6060分带法从格林尼治零度经线起,自东半球向西半球,每经差60分为一个投影带,即东经0~6,6~12,12~18,….174~180,用阿拉伯数字1,2,3,4….60表示投影带号,全球共分为60个投影带。

东半球中央经线的计算公式为:L0=(6n−3)0L0=(6n−3)0n表示投影带号,n<30西半球中央经线的计算公式为:L0=(6n−3)0−3600L0=(6n−3)0−3600n表示投影带号,n>303030分带法从东经1030′1030′算起,自东半球向西半球每3030为一带,将全球划分为120个投影带,1030′−4030′1030′−4030′,4030′−7030′4030′−7030′….其中央经线的位置为30,60,90,150…1800,−1770…−30,30,60,90,150…1800,−1770…−30,。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

2、投影后角度不变形。又叫保角映射或叫正形投影。条件是 在微小范围内成立。正形投影又叫等角投影。
采用正形投影,在有限范围内,使地形图上的图形与椭球面 上的相应图形保持相似。
第一节 地图投影概念和正形投影性质
三、正形投影的一般条件
正形投影必要和充分的条件是满足柯西—黎曼方程:
推导过程:由长度比的定义顾及正形投影的特性导出。
ds dS
2
dx 2 dy 2 2 2 (MdB) (NcosBdL)
2、引入等量纬度q,将x、y表为q、l的函数;
3、对 x=f1(q,l),y=f2(q,l)取全微分,引入符号E、F、G;
4、根据长度比m与方向A无关,F=0,E=G; 5、由E=G、F=0,得一般条件:
第二节 高斯投影与国家平面直角坐标系
四、高斯投影计算内容
第二种方法的具体推算内容如下:
1、将起算点的大地坐标(B1,L1)换算为高斯平面坐标(x1,y1) 2、将起算边的大地方位角A12改换为平面坐标方位角T12;
T12=A12-γ+δ12
式中,γ为子午线收敛角,δ12为方向改正。 3、将起算边的大地线长度S12归算为高斯平面上的直线长度D12:
一、平面子午线收敛角的计算
1、由大地坐标计算平面子午线收敛角γ
第四节 椭球面上的方向和长度归算至高斯投影平面
一、平面子午线收敛角的计算
2、由平面直角坐标计算平面子午线收敛角γ
第四节 椭球面上的方向和长度归算至高斯投影平面 二、方向改正计算
x C
2dδ
y2 - y1= y
ε
2
δ
2dδ dξ
x
21
x2 - x1= x
5
5′
第三节 高斯投影坐标计算
高斯投影坐标正算——由(B,L)求(x,y)
高斯投影坐标反算——由(x,y)求(B,L)
第三节 高斯投影坐标计算
一、由(B,L)计算(x,y)--正算
6 式中,X为由赤道至纬度B的子午线弧长, 为计算点P点与中央子午线
的经差。N为卯酉圈曲率半径,t=tanB,
η=e′cosB。 L-L0若以度为单位,则ρ=57.295779513;
第六章 高斯投影 及其计算
中国矿业大学(北京)地测学院
第六章 高斯投影及其计算
第一节 地图投影概念和正形投影性质 第二节 高斯投影与国家平面直角坐标系(概念) 第三节 高斯投影坐标计算(重点) 第四节 椭球面上的方向和长度归算至高斯投影平面(重点) 第五节 高斯投影坐标换带计算(重点) 第六节 通用横轴墨卡托投影和兰勃特投影简介
2 3
第二节 高斯投影与国家平面直角坐标系
三、高斯投影的分带
6°带带号N和中央子午线经度 LN的关系式:LN=6N-3
3°带带号n和中央子午线经度 Ln的关系式:Ln=3n
6°带与3°带带号之间的关系为:n=2N-1 国家统一坐标的表示方法:X坐标为正,Y坐标加500km后 前面冠以带号。如在39带中Y坐标自然值分别为12345.678 和-12345.678m,国家统一坐标分别为39512345.678和 39487654.322m。但在坐标计算中应去掉带号,减去500km, 恢复坐标自然值。
L-L0若以分为单位,则ρ=3437.7467708;
L-L0若以秒为单位,则ρ=206264.80625。
(B,L)计算(x,y)正算公式中子午弧长 X的计算(见本书150页公式5-41)
3 5 7 X C[β B (β cosB β cos B β cos B β cos B)sinB]; 0 2 4 6 8
第三节 高斯投影坐标计算
一、由(B,L)计算(x,y)--正算
推证过程: 1、高斯投影坐标正算函数式 2、根据正形投影的一般公式 x+iy=f(q+il)以及高斯投影的条件推 导正算公式,可以将一般公式在q处展为il 的台劳级数。 3、根据中央子午线长度比 m=1,有
2
4、由
求各阶导数
5、将各阶导数代入上式得最后正算公式。
f点的坐标为(x+i0)对应大地坐标为(Bf,0) 3、对于f点,F(x)=q,即(F(x))f=qf ,qf为f点等量纬度。 4、求各阶导数并代回上式得:
5、将q、qf变为B-Bf 由B=φ(q)在qf处展为台劳级数,由等量纬度q与B的微分公式求得 各系数,得B的计算公式和l的计算公式。
第四节 椭球面上的方向和长度归算至高斯投影平面
第一节 地图投影概念和正形投影性质
几何投影--垂直投影
第一节 地图投影概念和正形投影性质
几何投影--中心投影
第一节 地图投影概念和正形投影性质
(二)投影变形
投影变形不可避免。有角度变形、长度变形和面积变形三种。根 据实际需要选择某种变形为零或使其减小到某一适当程度。如高斯投 影,保持角度不变形,但长度和面积有变形。
第二节 高斯投影与国家平面直角坐标系
一、高斯投影的基本概念
高斯投影又称横轴椭圆柱等角投影。在高 斯投影平面上,中央子午线和赤道的投影 都是直线,分别为高斯平面直角坐标系的 X轴和Y轴。
N
x N
O
O
y
S
S
第二节 高斯投影与国家平面直角坐标系
一、高斯投影的基本概念
高斯投影的条件 (1)投影后角度不产生变形,满足正形投影要求; (2)中央子午线投影后是一条直线; (3)中央子午线投影后长度不变,其投影长度比恒等 于1。 高斯投影除了在中央子午线上没有长度变形外,不 在中央子午线上的各点,其长度比都大于1,且离开中 央子午线愈远,长度变形愈大。
第三节 高斯投影坐标计算
二、由(x,y)计算(B,L)--反算
式中,带下标fห้องสมุดไป่ตู้各变量为f点 (又叫底点)的有关变量,如
B f ,M f ,N f , tf , ηf
分别为底点纬度,子午圈半径, 卯酉圈半径等。 Bf可由子午弧长公式反算得到。
第三节 高斯投影坐标计算
反算推证过程: 1、将正形投影公式写成反函数形式 2、在f点(Bf,0)将上式展为iy的台劳级数
第二节 高斯投影与国家平面直角坐标系
四、高斯投影计算内容
1 、由椭球面上各点大地坐标( B , L )求解各点高斯 平面坐标( x , y ):先在椭球面上解算球面三角形, 推算各边大地方位角,解算各点大地坐标,然后求解 各点的高斯平面坐标。(计算工作量大) 2 、将椭球面上起算元素和观测元素归算至高斯投影 平面,然后解算平面三角形,推算各边坐标方位角, 在平面上进行平差计算,求解各点的平面直角坐标。
3、长度变形m-1与横坐标y的关系
y/(km) 10
20
30
40
50
100
150
200
250
300
长度变 形m-1
1/810000
1/202000
1/90000
1/50000
1/32000
1/8000
1/3500
1/2000
1/1300
1/900
第二节 高斯投影与国家平面直角坐标系
三、高斯投影的分带
正轴投影——圆柱面中心轴与椭球短轴重合,圆柱面与赤道相切。
横轴投影——圆柱面中心轴与椭球长轴重合,圆柱面与某一子午圈相切。
斜轴投影——圆柱面中心轴与椭球长、短轴都不重合,位于两者之间。
第一节 地图投影概念和正形投影性质
二、正形投影特性
1、任一点上,投影长度比m为一常数,不随方向而变,a=b。 长度比仅与点位臵有关,不同点投影有不同的长度比。
D12=S12+△S
式中△S为距离改正。 4、对于椭球面上三角网的各观测方向和观测边长分别进行方向改正和距离 改正,归算为高斯平面上的直线方向和直线距离。组成平面三角网,平差计 算,推求各控制点的平面直角坐标。
第二节 高斯投影与国家平面直角坐标系
高斯投影坐标计算、平面子午线收敛角计算、方向改正计 算、距离改正计算,统称为高斯投影计算。
3 2 45 4 175 6 11025 8 e e e e , 4 64 256 16384 15 4 175 6 3675 8 β2 β0 1,β e e e , 4 32 384 8192 35 735 8 315 8 β6 e 6 e ,β8 e , 96 2048 1024 C a 2 /b;对应不同的椭球 ,其参数不同,所计算 的X不同。 式中, β0 1 上式可以变换为: X B (β2 cosB β4 cos3 B β6 cos5 B β8 cos7 B)sinB, Cβ 0 依此公式按迭代法可以 由X求B。
D T
dδ A
B dζ
dx
B′ y A′ A
B


δ
12
1
y
η
第四节 椭球面上的方向和长度归算至高斯投影平面
二、方向改正计算
方向改正——正形投影后,椭球面上大地线投影到平面上仍为曲线, 化为直线方向所加的改正δ。 适用于三、四等三角测量的方向改正计算公式
上式的计算精度为0.1″。
第四节 椭球面上的方向和长度归算至高斯投影平面
第一节 地图投影概念和正形投影性质
(四)地图投影的分类
等角投影——投影后角度不变,保持小范围内图形相似。 等面积投影——用于某些专题地图,投影后面积不变。 圆锥面投影——圆锥面与椭球体在某一纬圈相切或某两纬圈相割,按数学 投影。
平面投影——投影平面与椭球面在某一点相切,按数学投影建立函数关系。
圆柱面投影——圆柱面或椭圆柱面与椭球面在赤道或某一子午面上相切, 按数字投影。
三、距离改正计算
第二节 高斯投影与国家平面直角坐标系
二、高斯投影的长度比和长度变形
相关文档
最新文档