高斯投影坐标正算公式

合集下载

第6章 高斯投影简介

第6章 高斯投影简介

(4)高斯平面投影的特点:
高斯平面直角坐标系与大地坐标系
1、高斯投影坐标正算公式
(1)高斯投影正算:已知椭球面上某点的大地坐标 L , B ,求该点 在高斯投影平面上的直角坐标 x , y ,即 L , B ( x , y ) 的坐标变换。 (2)投影变换必须满足的条件: 中央子午线投影后为直线; 中央子午线投影后长度不变; 投影具有正形性质,即正形投影条件。 (3)投影过程 在椭球面上有对称于中央子午线的两点 P1 和 P 2 ,它们的大地坐标 分别为(L , B)及(l , B ),式中 为椭球面上 P 点的经度与中央子 l 午线 ( L )的经度差: L L 0 , l 点在中央子午线之东, l 为正,在西 P 则为负,则投影后的平面坐标一定 P1( x , y ) 为 P2 ( x , y ) 和。
归算的内容
水平观测方向
地面长度归算
水平观测方向归算到椭球面
三差改正:水平方向归算到椭球面上,
需进行垂线偏差改正、标高差改正和截 面差改正,通常把这三项改正简称为三 差改正
① 垂线偏差改正(δu)
定义:地面上以铅垂线为准观测的水平方向值,归 算为以椭球面法线为准的水平方向值时,顾及测站 点垂线偏差的影响所加的改正
第六章 高斯投影简介
高斯投影概述
1、投影与变形
地图投影:就是将椭球面各元素(包括坐标、方向和长度)按一 定的数学法则投影到平面上。研究这个问题的专门学科叫地图投影 学。可用下面两个方程式(坐标投影公式)表示:
x F1 ( L , B ) y F2 ( L , B )
式中L,B是椭球面上某点的大地坐标,而是x,y该点投影后的平面直角 坐标。
f f
t

高斯投影高斯投影正算公式

高斯投影高斯投影正算公式
①高斯投影产生的背景
高斯-克吕格投影也称等角横切椭圆柱投 影,它可看作是等角圆柱投影(墨卡托投 影,1569)的一种,它由德国科学家高斯 处理三角测量成果时首先提出,后经克吕 格完善(1919) ,我国于1952年起正式采 用高斯-克吕格投影。
四个世纪以来,世界各国都用墨卡托投影作 为海图的数学基础。当代常用于较大比例尺 分幅海图或赤道附近的航空图。
《大地测量学基础》(FOUNDATION OF GEODESY)
高斯-克吕格投影 高斯平面坐标系与大地坐标系
的关系(1)
测绘学院一系大地测量教研室
上节课内容回顾
☺ 长度比? m d s
dS
☺ 椭球面到平面的长度比在什么方 向取极值?
子午方向和卯酉方向 MNcosB
☺ 最大角度变形? sin a b
② 分带的方法
1) 6°带划分 (n为带号 )
6°带中央子午线的经度计算公式 L0 6 n3
已知6°带中央子午线的经度反算带号
n
1 6
(L0
3
)
计算任意经度所在投影带的带号公式
nL的 整 数 商 ( 1有 余 数 时 ) 6
2、高斯投影的分带
Zone-dividing of Gauss Projection
② 分带的方法
2) 3°带划分 (n'为带号 )
3°带中央子午线的经度计算公式 L0 3 n
已知3°带中央子午线的经度反算带号 n L 0 3
计算任意经度所在投影带的带号公式 n L 1.5 1 3
③ UTM分带的方法
UTM的分带是从经度180°起向东每6°为一 带,即与国际百万分之一地形图的划分一致;
(135°02′30″)
南海南沙群岛的曾母 暗沙(3°52′)

高斯平面直角坐标与大地坐标的相互转换——高斯投影的正算与反算.

高斯平面直角坐标与大地坐标的相互转换——高斯投影的正算与反算.

昆明冶金高等专科学校测绘学院 (4)计算公式
3 2 2 2 4 ( 5 3 t 9 t ) y f f f f 2M f N f 2 4M f N 3 f tf 2 4 6 (6 1 9 0t f 4 5t f ) y 7 2 0M f N 5 f 1 1 2 2 3 l y (1 2t f f ) y 3 N f co s B f 6 N f co s B f 1 2 5 (5 2 8t 2 t4 2 2 f 24 f 6 f 8 f t f )y 5 1 2 0N f co s B f B Bf tf y2 tf
式中:

2 e 2 cos2 B
t 2 tan2 B l (L L0) X为B对应子午线弧长 N为卯酉圈曲率半径 20626 5
昆明冶金高等专科学校测绘学院
2
高斯投影坐标反算公式
(1)高斯投影反算:
已知某点 x, y ,求该点 L, B ,即 x, y ( L, B) 的坐标变换。 (2)投影变换必须满足的条件
昆明冶金高等专科学校测绘学院
二、高斯投影坐标正反算得实用公式及算例
1 高斯投影坐标正算公式 (1)高斯投影正算: 已知某点的 L, B ,求该点的 x, y ,即 (2)投影变换必须满足的条件: 中央子午线投影后为直线; 中央子午线投影后长度不变; 投影具有正形性质,即正形投影条件。 (3)投影过程 在椭球面上有对称于中央子午线的两点 P1 和 P2 ,它们的大地坐标 分别为 ( L1 , B1 )或(l1 , B1)及 (L2 , B2)或(l2 , B2 ) 式中 l 为椭球面上点的经 度与中央子午线 ( L0 ) 的经度差:l L L0 ,点在中央子午线之东, l 为正,在西则为负,则投影后的平面坐标一定为P1 ( x1 , y1 ) 和 P2 ( x 2 , y 2 ) 。

高斯投影坐标计算

高斯投影坐标计算

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

高斯平面直角坐标系

高斯平面直角坐标系

大地测量学基础
4.9 高斯平面直角坐标系 三、高斯投影坐标正反算公式
(4)反算公式
当l<3.5°时,上式换算精度达0.0001″。 欲使换算精确至0.01″,可对上式简化成:
大测量学基础
4.9 高斯平面直角坐标系 三、高斯投影坐标正反算公式
平 时 作 业 用编程进行高斯投影正反算。 已知
B 51 3843.9023 L 111 0213.1360
大地测量学基础
4.9 高斯平面直角坐标系 三、高斯投影坐标正反算公式
即有:
在数学上,F1为 l 的偶函数,F2为 l 的奇函数。 因为在每带中,l/ρ˝不大,是一个微小量,可展成幂级 数。
m0,m1,m2,…,是待定系数,它们都是纬度B的函 数。
大地测量学基础
4.9 高斯平面直角坐标系 三、高斯投影坐标正反算公式
大地测量学基础
4.9 高斯平面 直角坐标系
大地测量学基础
4.9 高斯平面直角坐标系 三、高斯投影坐标正反算公式
三、高斯投影坐标正反算公式 1、高斯投影坐标正反算的定义 (1)高斯投影正算: 已知椭球面上某点的大地坐标B、L,求其 该点在高斯平面直角坐标系中的坐标x、y的工作 叫高斯投影正算。 (2)高斯投影反算: 已知椭球面上某点在高斯平面直角坐标系中 的坐标x、y,求其该点的大地坐标B、L的工作 叫高斯投影反算。
大地测量学基础
4.9 高斯平面直角坐标系 三、高斯投影坐标正反算公式
(3)反算公式推导思路: 和正算公式基本一样,也是根据高斯投影的3个条件来 推导的。 ①由对称条件,同样可得: 把B、l 展成y的幂级数,而φ1为y的偶函数, φ2为y的奇 函数。
式中 n 0 ,n 1 ,n 2 … 是待定系数,它们都是纵坐标 x 的函数 ,与y无关。

第四章 7高斯投影坐标正反算

第四章 7高斯投影坐标正反算
第四章 Ⅶ 高斯投影坐标正反算
——正形投影的一般条件 ——高斯投影坐标正算 ——高斯投影坐标反算 ——高斯投影几何解释
提前在黑板上写出四个m2
上一讲应掌握的内容
1、地图(数学)投影:将椭球面上元素(包括坐标,方位和 距离)按一定的数学法则投影到可展平面上。 x F1 ( L, B) 坐标投影公式: y F2 ( L, B) 2、地图投影变形几个概念: 长度比,主方向,变形椭圆 3、四种投影变形: 长度变形,方向变形,角度变形,面积变形
x m0 m 2 l 2 m 4 l 4 y m1l m3 l 3 m5 l 5
分别对l 和q 求偏导数
2) 由第三个条件正形投影条件
y x x y 和 l q l q
dm0 dm2 2 dm4 4 2 4 m1 3m3 l 5m5 l dq dq l dq l 2m l 4m l 3 dm1 l dm3 l 3 )
将各系数代入,略去高次项,得高斯投影坐标正算公式 精度为0.001m
xX N N sin B cos Bl 2 + sin B cos 3 B(5 - t 9 2 4 4 )l 4 + 2 24
N sin B cos 5 B(61 - 58t 2 t 4 )l 6 720
dl tan Adq
2 2 2 2 E ( dq ) 2 F tan A ( dq ) G tan A ( dq ) m2 2 2 2 r2 ( dq ) tan A ( dq )
E 2 F tan A G tan 2 A = r 2 sec 2 A E cos 2 A 2 F sin A cos A G sin 2 A = r2

大地坐标BLH转平面坐标xyh(高斯投影坐标正算)Java版

大地坐标BLH转平面坐标xyh(高斯投影坐标正算)Java版

⼤地坐标BLH转平⾯坐标xyh(⾼斯投影坐标正算)Java版技术背景 做过位置数据处理的⼩伙伴基本上都会遇到坐标转换,⽽基于⾼斯投影原理的⼤地坐标转平⾯坐标就是其中⼀种坐标转换,坐标转换的⽬的就是⽅便后⾯数据的处理⼯作,⼤地坐标转⾼斯平⾯坐标常⽤的有两种,即3°带和6°带,具体采⽤哪种根据实际情况⽽定。

计算原理 6°带带号n与相应的中央⼦午线L0经度的关系为: 3°带带号n’与相应的中央⼦午线L0’经度的关系为: 设参考椭球的长半轴为 a,第⼀偏⼼率为 e,并令: 设中央⼦午线的经度为 L0,再记: 则⾼斯投影正算公式为: 其中: 没学过测绘的同学对以上原理不是很理解,也很正常,⼤家可以查阅相关《⼤地测量学》书籍,具体武⼤版还是矿⼤版的,没什么区别。

具体实现 具体实现平台依然是IBM的Eclipse软件,编程语⾔为Java,下⾯是以3°带为例,进⾏⾼斯投影坐标正算,具体内容请看代码1package package1;23public class BLH_xyh {45public static double a = 6378137;6public static double e = Math.sqrt(0.0066943799013);7public static double scale_wide = 3;8public static doubleπ = 3.14159265358979323846;910public static void main(String[] args) {11// TODO ⾃动⽣成的⽅法存根12 Point3d xyh = the_coordinates_are_counting(36.0307523111,120.184664478,9.7065);13 System.out.println(xyh.getX()+","+xyh.getY()+","+xyh.getZ());14 }15public static Point3d the_coordinates_are_counting(double B,double L,double H){16double A_ = 117 +3*e*e/418 +45*e*e*e*e/6419 +175*e*e*e*e*e*e/25620 +11025*e*e*e*e*e*e*e*e/1638421 +43659*e*e*e*e*e*e*e*e*e*e/65536;22double B_ = 3*e*e/4+15*e*e*e*e/1623 +525*e*e*e*e*e*e/51224 +2205*e*e*e*e*e*e*e*e/204825 +72765*e*e*e*e*e*e*e*e*e*e/65536;26double C_ = 15*e*e*e*e/6427 +105*e*e*e*e*e*e/25628 +2205*e*e*e*e*e*e*e*e/409629 +10395*e*e*e*e*e*e*e*e*e*e/16384;30double D_ = 35*e*e*e*e*e*e/51231 +315*e*e*e*e*e*e*e*e/204832 +31185*e*e*e*e*e*e*e*e*e*e/13072;33/*34 double E_ = 315*e*e*e*e*e*e*e*e/1638435 +3465*e*e*e*e*e*e*e*e*e*e/65536;36 double F_ = 693*e*e*e*e*e*e*e*e*e*e/13072;37 */3839doubleα = A_*a*(1-e*e);40doubleβ = -B_*a*(1-e*e)/2;41doubleγ = C_*a*(1-e*e)/4;42doubleδ = -D_*a*(1-e*e)/6;43/*44 double ε = E_*a*(1-e*e)/8;45 double ζ = -F_*a*(1-e*e)/10;46 */4748double C0 = α;49double C1 = 2*β+4*γ+6*δ;50double C2 = -8*γ-32*δ;51double C3 = 32*δ;5253double x,y,sign;54double scale_number = Math.floor(L/scale_wide);55if(L > (scale_number * scale_wide + scale_wide/2)){56 scale_number =scale_number + 1;57 sign = -1;58 }else{59 sign = 1;60 }6162double L0 = scale_wide*scale_number;63double l = Math.abs(L-L0);64 B = B*π/180;65 l = l*π/180;66double t = Math.tan(B);67double m0 = Math.cos(B)*l;68doubleη = Math.sqrt(e*e*Math.pow(Math.cos(B),2)/(1-e*e));69double N = a/Math.sqrt(1-e*e*Math.pow(Math.sin(B), 2));7071double X0 = C0*B+Math.cos(B)*(C1*Math.sin(B)+C2*Math.pow(Math.sin(B),3)+C3*Math.pow(Math.sin(B), 5));7273 x = X074 +N*t*m0*m0/275 +N*t*m0*m0*m0*m0*(5-t*t+9*η*η+4*η*η*η*η)/2476 +N*t*m0*m0*m0*m0*m0*m0*(61-58*t*t+t*t*t*t)/720;77 y = N*m0+78 N*m0*m0*m0*(1-t*t+η*η)/6+79 N*m0*m0*m0*m0*m0*m0*(5-18*t*t+t*t*t*t+14*η*η-58*η*η*t*t)/120;8081 y = y*sign+500000;8283double h = H;8485 Point3d xyh = new Point3d(x,y,h);86return xyh;87 }88 }其中Point3d的定义如下:1package package1;23public class Point3d {4private double x;5private double y;6private double z;7public Point3d(double x,double y,double z){8this.x=x;9this.y=y;10this.z=z;11 }12public double getX() {13return x;14 }15public void setX(double x) {16this.x = x;17 }18public double getY() {19return y;20 }21public void setY(double y) {22this.y = y;23 }24public double getZ() {25return z;26 }27public void setZ(double z) {28this.z = z;29 }30 }测试数据为36.0307523111,120.184664478,9.7065,运⾏结果如下:⾄此结束致谢 感谢⼭东科技⼤学北⽃星光创客兴趣学习⼩组的王⽼师对于原理⽂档的整理以及郑** C++代码的技术分享!参考⽂档1、⼭东科技⼤学”北⽃星光创客”兴趣学习⼩组GNSS技术⽂档。

1 高斯投影坐标正算公式

1 高斯投影坐标正算公式

1 高斯投影坐标正算公式(1)高斯投影正算:已知椭球面上某点的大地坐标,求该点在高斯投影平面上的直角坐标,即的坐标变换。

(2)投影变换必须满足的条件中央子午线投影后为直线;中央子午线投影后长度不变;投影具有正形性质,即正形投影条件。

(3)投影过程在椭球面上有对称于中央子午线的两点和,它们的大地坐标分别为()及(),式中为椭球面上点的经度与中央子午线的经度差:, 点在中央子午线之东, 为正,在西则为负,则投影后的平面坐标一定为和。

(4)计算公式当要求转换精度精确至0.OOlm时,用下式计算:2 高斯投影坐标反算公式(1)高斯投影反算:已知某点的高斯投影平面上直角坐标,求该点在椭球面上的大地坐标,即的坐标变换。

(2)投影变换必须满足的条件坐标轴投影成中央子午线,是投影的对称轴;轴上的长度投影保持不变;投影具有正形性质,即正形投影条件。

(3)投影过程根据计算纵坐标在椭球面上的投影的底点纬度,接着按计算()及经差,最后得到、。

(4)计算公式当要求转换精度至时,可简化为下式:3高斯投影相邻带的坐标换算(1)产生换带的原因高斯投影为了限制高斯投影的长度变形,以中央子午线进行分带,把投影范围限制在中央子午线东、西两侧一定的范围内。

因而,使得统一的坐标系分割成各带的独立坐标系。

在工程应用中,往往要用到相邻带中的点坐标,有时工程测量中要求采用带、带或任意带,而国家控制点通常只有带坐标,这时就产生了带同带(或带、任意带)之间的相互坐标换算问题,如图所示:(2)应用高斯投影正、反算公式间接进行换带计算计算过程把椭球面上的大地坐标作为过渡坐标。

首先把某投影带(比如Ⅰ带)内有关点的平面坐标,利用高斯投影反算公式换算成椭球面上的大地坐标,进而得到;然后再由大地坐标,利用投影正算公式换算成相邻带的(第Ⅱ带)的平面坐标。

在这一步计算时,要根据第Ⅱ带的中央子午线来计算经差,亦即此时。

算例在中央子午线的Ⅰ带中,有某一点的平面直角坐标,,现要求计算该点在中央子午线的第Ⅱ带的平面直角坐标。

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

高斯投影坐标正算公式
高斯投影坐标正反算公式
2.2.2. 1高斯投影坐标正算公式: B, x,y
高斯投影必须满足以下三个条件:
⑴中央子午线投影后为直线;⑵中央子午线投影后长度不变;⑶投影具有正形性质,即正形投影条件。

由第一条件知中央子午线东西两侧的投影必然对称于中央子午线,即
式中,x为的偶函数,y为的奇函数;,即,
如展开为的级数,收敛。

(2-10)
式中是待定系数,它们都是纬度B的函数。

由第三个条件知:
分别对和q求偏导数并代入上式
(2-11)
上两式两边相等,其必要充分条件是同次幂前的系数应相等,即
(2-12)
(2-12)是一种递推公式,只要确定了就可依次确定其余各系数。

由第二条件知:位于中央子午线上的点,投影后的纵坐标x应等于投影前从赤道量至该点的子午线弧长X,即(2-10)式第一式中,当时有:
(2-13)
顾及(对于中央子午线)
得:
(2-14,15)
(2-16)
依次求得并代入(2-10)式,得到高斯投影正算公式
(2-17)
2.2.2. 2高斯投影坐标反算公式
x,y B,
投影方程:
(2-18)
高斯投影坐标反算公式推导要复杂些。

⑴由x求底点纬度(垂足纬度),对应的有底点处的等量纬度,求x,y与
的关系式,仿照式有,
由于y和椭球半径相比较小(1/16.37),可将展开为y的幂级数;又由于是对称投影,q必是y的偶函数,必是y的奇函数。

(2-19)
是待定系数,它们都是x的函数.
由第三条件知:

, (2-20)
(2-19)式分别对x和y求偏导数并代入上式
上式相等必要充分条件,是同次幂y前的系数相等,
第二条件,当y=0时,点在中央子午线上,即x=X,对应的点称为底点,其纬度为底点纬度,也就是x=X时的子午线弧长所对应的纬度,设所对应的等量纬度为。

也就是在底点展开为y的幂级数。

由(2-19)1式
依次求得其它各系数
(2-21)
(2-21)1
…………
将代入(2-19)1式得
(2-22)1
(2-22)
将代入(2-19)2式得(2-23)2式。

(最后表达式) ⑵求与的关系。

由式知:
(2-23)
(2-24)
按台劳级数在展开
(2-25)
(2-25) 由式可求出各阶导数:
(2-26)
(2-27)1
(2-27)2
…………………
将式(2-22)1(2-22)(2-26)(2-27)入(2-25)并按y幂集合得高斯投影坐标反算公式(2-28)
(2-28)
归纳由求的基本思想:由点得到底点,将底点f作为过渡,也就是说将坐标原点o移到f点,先求关系式,再将
关系式代入关系式得关系式,最后将坐标原点移回到o点,从而求得点。

2.2.2. 3高斯投影坐标正反算公式的几何解释:
(图4:高斯投影坐标正反算公式的几何解释示意图)
⑴当B=0时x=X=0,y则随的变化而变化,这就是说,赤道投影
为一直线且为y轴。

当=0时,则y=0,x=X,这就是说,中央子午线投影亦为直线,且为x轴,其长度与中央子午线长度相等。

两轴的交点为坐标原点。

⑵当=
常数时(经线),随着B值增加,x值增大,y值减小,这就告诉我们,经线是凹向中央子午线的曲线,且收敛于两极。

又因,即当用-B代替B时,y 值不变,而x值数值相等符号相反,这就说明赤道是投影的对称轴。

⑶当B=常数时(纬线),随着的增加,x值和y值都增大,这就是说,纬线是凸向赤道的曲线。

又当用-代替时,x值不变,而y值数值相等符号相反,这就说明,
中央子午线是投影对称轴。

由于满足正形投影条件,所以经线和纬线的投影是互相垂直的。

⑷距中央子午线愈远的子午线,投影后弯曲愈厉害,表明长度变形愈大。

相关文档
最新文档