大地测量学课程第21讲-高斯投影坐标反算公式

高斯投影正反算

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

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

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

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

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

2）由于高斯投影是换带投影，在每带经差ｌ是不大的，ｌρ是一个微小量，所以可以将X=X （ｌ，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 ⎧====⎪⎪⎨⎪=-=-=-=-⎪⎩经过计算可以得出232244524632235242225sin cos sin cos (594)224 sin cos (6158)720cos cos (1)6cos (5181458)120N N x X B B l B B t l NB B t t l Ny N B l B t l NB t t t l ηηηηη=+⋅+-+++-+=⋅+-++-++-三、高斯投影坐标反算公式推导1.思路：级数展开，应用高斯投影三个条件，待定系数法求解。

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

①高斯投影产生的背景
高斯－克吕格投影也称等角横切椭圆柱投影，它可看作是等角圆柱投影（墨卡托投影，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′）

高斯平面直角坐标系与大地坐标系相互转化

高斯平面直角坐标系与大地坐标系相互转化高斯平面直角坐标系与大地坐标系转换 1. 高斯投影坐标正算公式(1) 高斯投影正算:已知椭球面上某点的大地坐标(L，B)，求该点在高斯投影平面上的直角坐标(x，y)，即(L，B)->(x，y)的坐标变换。

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

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

(4) 计算公式 4 ' ' 2 2 3 4 ' ' 2 ' ' 2 ' ' ) 9 5 ( cos sin 2 sin 2 l t B B N Bl N X x 5 ' ' 4 2 5 5 ' ' 3 ' ' 2 2 3 ' ' ' ' ' ' ) 18 5 ( cos 120 ) 1 ( 6 cos l t t B N l t B N Bl N y 当要求转换精度精确至0.001m时，用下式计算: 6 ' ' 4 2 5 6 ' ' 4 ' ' 4 2 2 3 4 ' ' 2 ' ' 2 ' ' ) 58 61 ( cos sin 720 ) 4 9 5 ( cos sin 24 sin 2 l t t B B N l t B B N Bl N X x5 ' ' 2 2 2 4 2 5 5 ' ' 3 ' ' 2 2 3 3 ' ' ' ' ' ' ) 58 14 18 5 ( cos 720 ) 1( cos 6 cos l t t t B N l t B N Bl N y2. 高斯投影坐标反算公式(1) 高斯投影反算:已知某点的高斯投影平面上直角坐标(x,y)，求该点在椭球面上的大地坐标(L,B)，即(x,y)->(L,B)的坐标变换。

高斯平面直角坐标系与大地坐标系相互转化

高斯平面直角坐标系与大地坐标系相互转化高斯平面直角坐标系与大地坐标系转换 1. 高斯投影坐标正算公式 (1) 高斯投影正算:已知椭球面上某点的大地坐标(L ，B)，求该点在高斯投影平面上的直角坐标(x ，y)，即(L ，B)->(x ，y)的坐标变换。

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

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

(4) 计算公式 4 ' ' 2 2 3 4 ' ' 2 ' ' 2 ' ' ) 9 5 ( cos sin 2 sin 2 l t B B N BlN X x 5 ' ' 4 2 5 5 ' ' 3 ' ' 2 2 3 ' ' ' ' ' ' ) 18 5 ( cos 120 ) 1 ( 6 cos l t t B N l t B NBl N y 当要求转换精度精确至0.001m 时，用下式计算: 6 ' ' 4 2 5 6 ' ' 4 ' ' 4 2 2 3 4 ' ' 2 ' ' 2 ' ' ) 58 61 ( cos sin 720 ) 4 9 5( cos sin 24 sin 2 l t t B B N l t B B N Bl N X x5 ' ' 2 2 2 4 2 5 5 ' ' 3 ' ' 2 2 3 3 ' ' ' ' ' ' ) 58 14 18 5 ( cos 720 ) 1( cos 6 cos l t t t B N l t B N Bl N y2. 高斯投影坐标反算公式 (1) 高斯投影反算:已知某点的高斯投影平面上直角坐标(x,y)，求该点在椭球面上的大地坐标(L,B)，即(x,y)->(L,B)的坐标变换。

高斯平面直角坐标系

大地测量学基础
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/ρ˝不大，是一个微小量，可展成幂级数。
m０，ｍ１，ｍ２，…，是待定系数，它们都是纬度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 ０，ｎ１，ｎ２ … 是待定系数，它们都是纵坐标 x 的函数 ,与y无关。

高斯投影坐标正反算公式

高斯投影坐标正反算公式未知2010-04-03 10:47:15 本站§高斯投影坐标正反算公式任何一种投影①坐标对应关系是最主要的；②如果是正形投影，除了满足正形投影的条件外（ C-R 偏微分方程），还有它本身的特殊条件。

1.1 高斯投影坐标正算公式： B, x,y高斯投影必须满足以下三个条件：①中央子午线投影后为直线；②中央子午线投影后长度不变；③投影具有正形性质，即正形投影条件。

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

（ 8-33 ）式中是待定系数，它们都是纬度 B 的函数。

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

由第二条件知 : 位于中央子午线上的点，投影后的纵坐标 x 应等于投影前从赤道量至该点的子午线弧长 X ，即 (8-33) 式第一式中，当时有：(8-36)顾及 ( 对于中央子午线 )得：(8-37,38)(8-39)依次求得并代入 (8-33) 式，得到高斯投影正算公式(8-42)1.2 高斯投影坐标反算公式x,y B,投影方程：(8-43)满足以下三个条件：①x 坐标轴投影后为中央子午线是投影的对称轴；② x 坐标轴投影后长度不变；③投影具有正形性质，即正形投影条件。

高斯投影坐标反算公式推导要复杂些。

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

(8-45)是待定系数，它们都是 x 的函数 .由第三条件知：，， (8-21)(8-45) 式分别对 x 和 y 求偏导数并代入上式上式相等必要充分条件，是同次幂 y 前的系数相等，第二条件，当 y=0 时，点在中央子午线上，即 x=X ，对应的点称为底点，其纬度为底点纬度，也就是 x=X 时的子午线弧长所对应的纬度，设所对应的等量纬度为。

高斯投影正反算编程一.高斯投影正反算基本公式

高斯投影正反算编程一．高斯投影正反算基本公式（1）高斯正算基本公式（2）高斯反算基本公式以上主要通过大地测量学基础课程得到，这不进行详细的推导，只是列出基本公式指导编程的进行。

二．编程的基本方法和流程图（1）编程的基本方法高斯投影正反算基本上运用了所有的编程基本语句，本文中是利用C++语言进行基本的设计。

高斯正算中对椭球参数和带宽的选择主要运用了选择语句。

而高斯反算中除了选择语句的应用，在利用迭代算法求底点纬度还应用了循环语句。

编程中还应特别注意相关的度分秒和弧度之间的相互转换，这是极其重要的。

（2）相关流程图1）正算2）反算三．编程的相关代码（1）正算# include "stdio.h"# include "stdlib.h"# include "math.h"# include "assert.h"#define pi (4*atan(1.0))int i;struct jin{double B;double L;double L0;};struct jin g[100];main(int argc, double *argv[]){FILE *r=fopen("a.txt","r");assert(r!=NULL);FILE *w=fopen("b.txt","w");assert(r!=NULL);int i=0;while(fscanf(r,"%lf %lf %lf",&g[i].B,&g[i].L,&g[i].L0)!=EOF){double a,b;int zuobiao;printf("\n请输入坐标系：北京54=1，西安80=2，WGS84=3：");scanf("%d",&zuobiao);getchar();if(zuobiao==1){a=6378245;b=6356863.0187730473;}if(zuobiao==2){a=6378140;b=6356755.2881575287;}if(zuobiao==3){a=6378137;b=6356752.3142;} //选择坐标系//double f=(a-b)/a;double e,e2;e=sqrt(2*f-f*f);e2=sqrt((a/b)*(a/b)-1);//求椭球的第一，第二曲率//double m0,m2,m4,m6,m8;double a0,a2,a4,a6,a8;m0=a*(1-e*e);m2=3*e*e*m0/2;m4=5*e*e*m2/4;m6=7*e*e*m4/6;m8=9*e*e*m6/8;a0=m0+m2/2+3*m4/8+5*m6/16+35*m8/128;a2=m2/2+m4/2+15*m6/32+7*m8/16;a4=m4/8+3*m6/16+7*m8/32;a6=m6/32+m8/16;a8=m8/128;double Bmiao,Lmiao, L0miao;Bmiao=(int)(g[i].B)*3600.0+(int)((g[i].B-(int)(g[i].B))*100.0)*60.0+( g[i].B*100-(int)(g[i].B*100))*100.0;Lmiao=(int)(g[i].L)*3600.0+(int)((g[i].L-(int)(g[i].L))*100.0)*60.0+(g [i].L*100-(int)(g[i].L*100))*100.0;L0miao=(int)(g[i].L0)*3600.0+(int)((g[i].L0-(int)(g[i].L0))*100.0)*60. 0+(g[i].L0*100-(int)(g[i].L0*100))*100.0;double db;db=pi/180.0/3600.0;double B1,L1,l;B1=Bmiao*db;L1= Lmiao*db;l=L1-L0miao*db;//角度转化为弧度//double T=tan(B1)*tan(B1);double n=e2*e2*cos(B1)*cos(B1);double A=l*cos(B1);double X,x,y;X=a0*(B1)-a2*sin(2*B1)/2+a4*sin(4*B1)/4-a6*sin(6*B1)/6+a8*sin(8*B1)/8;//求弧长//double N=a/sqrt(1-e*e*sin(B1)*sin(B1));int Zonewide;int Zonenumber;printf("\n请输入带宽：3度带或6度带Zonewide=");scanf("%d",&Zonewide);getchar();if(Zonewide==3){Zonenumber=(int)((g[i].L-Zonewide/2)/Zonewide+1);}else if(Zonewide==6){Zonenumber=(int)g[i].L/Zonewide+1;}else{printf("错误");exit(0);}//选择带宽//doubleFE=Zonenumber*1000000+500000;//改写为国家通用坐标//y=FE+N*A+A*A*A*N*(1-T*T+n*n)/6+A*A*A*A*A*N*(5-18*T*T+T *T*T*T+14*n*n-58*n*n*T*T)/120;x=X+tan(B1)*N*A*A/2+tan(B1)*N*A*A*A*A*(5-T*T+9*n*n+4*n*n *n*n)/24+tan(B1)*N*A*A*A*A*A*A*(61-58*T*T+T*T*T*T)/720;printf("\n所选坐标系的转换结果：x=%lf y=%lf\n",x,y);fprintf(w,"%lf %lf\n",x,y);//输出结果到文本文件//}fclose(r);fclose(w);system("pause");return 0;}（2）反算# include "stdio.h"# include "stdlib.h"# include "math.h"# include "assert.h"#define pi (4*atan(1.0))double X,Y,B1,B2,B3,F,t;double m0,m2,m4,m6,m8;double a0,a2,a4,a6,a8,a1,b1;double BB,LL,Bf;double e,e1;int d,m,s,i,zuobiao;double sort(double,double);struct jin{double x;double y;double L0;};struct jin g[100];//x,y,L0为输入量：x,y坐标和中央子午线经度// main(int argc, double *argv[]){FILE *r=fopen("c.txt","r");assert(r!=NULL);FILE *w=fopen("d.txt","w");assert(r!=NULL);int i=0;while(fscanf(r,"%lf %lf %lf",&g[i].x,&g[i].y,&g[i].L0)!=EOF)//文件为空，无法打开//{double a1=6378245.0000000000;//克拉索夫斯基椭球参数//double b1=6356863.0187730473;double a75=6378140.0000000000;//1975国际椭球参数//double b75=6356755.2881575287;double a84=6378137.0000000000;//WGS-84系椭球参数//double b84=6356752.3142000000;double M,N;//mouyou圈曲率半径，子午圈曲率半径//double t,n;double A,B,C;double BB,LL,Bf,LL0,BB0;double a,b;printf("\n选择参考椭球：1=克拉索夫斯基椭球，2=1975国际椭球，3=WGS-84系椭球:");scanf("%d",&zuobiao);getchar();if(zuobiao==1){a=a1;b=b1;}if(zuobiao==2){a=a75;b=b75;}if(zuobiao==3){a=a84;b=b84;}//选择参考椭球，求解第一偏心率e,第二偏心率e1// Bf=sort(a,b);//调用求解底点纬度的函数//double q=sqrt(1-e*e*sin(Bf)*sin(Bf));double G=cos(Bf);M=a*(1-e*e)/(q*q*q);N=a/q;double H,I;A=g[i].y/N;H=A*A*A;I=A*A*A*A*A;t=tan(Bf);n=e1*cos(Bf);B=t*t;C=n*n;BB0=Bf-g[i].y*t*A/(2*M)+g[i].y*t*H/(24*M)*(5+3*B+C-9*B*C)-g[i] .y*t*I/(720*M)*(61+90*B+45*B*B);LL0=g[i].L0*pi/180.0+A/G-H/(6*G)*(1.0+2*B+C)+I/(120*G)*(5.0+28 *B+24*B*B+6*C+8*B*C);//利用公式求解经纬度//int Bdu,Bfen,Ldu,Lfen;double Bmiao,Lmiao;Ldu=int(LL0/pi*180);Lfen=int((LL0/pi*180)*60-Ldu*60);Lmiao=LL0/pi*180*3600-Ldu*3600-Lfen*60;Bdu=int(BB0/pi*180);Bfen=int((BB0/pi*180)*60-Bdu*60);Bmiao=BB0/pi*180*3600-Bdu*3600-Bfen*60;//将弧度转化为角度//printf("\n所选坐标系的转换结果：%d度%d分%lf秒%d 度%d分%lf秒\n",Bdu,Bfen,Bmiao,Ldu,Lfen,Lmiao);fprintf(w,"%d°%d’%lf”%d°%d’%lf”\n",Bdu,Bfen,Bmiao,Ldu,Lfen,Lmiao);//将结果输出到文本文件//}fclose(r);fclose(w);system("pause");return 0;}double sort(double a,double b){double e,e1;e=sqrt(1-(b/a)*(b/a));e1=sqrt((a/b)*(a/b)-1);double m0,m2,m4,m6,m8;double a0,a2,a4,a6,a8;m0=a*(1-e*e);m2=3*e*e*m0/2;m4=5*e*e*m2/4;m6=7*e*e*m4/6;m8=9*e*e*m6/8;a0=m0+m2/2+3*m4/8+5*m6/16+35*m8/128;a2=m2/2+m4/2+15*m6/32+7*m8/16;a4=m4/8+3*m6/16+7*m8/32;a6=m6/32+m8/16;a8=m8/128;B1=g[i].x/a0;do{F=-a2*sin(2*B1)/2+a4*sin(4*B1)/4-a6*sin(6*B1)/6+a8*sin(8*B1 )/8;B2=(g[i].x-F)/a0;B3=B1;B1=B2;} while(fabs(B3-B2)>10e-10);//利用迭代算法求解底点纬度//return B2; }。

高斯投影坐标反算公式

大地方位角＝坐标方位角－子午线收敛角＋方向改化 A1 2 1 2 1 2
高斯坐标反算实用步骤
1、根据高斯坐标确定带号、计算中央子午线经度 ①计算带号
n int( y / 1000000 )
②计算中央子午线经度六度带 L 6n 3 0 三度带 L0 3n
11 n 24 23 n 49
0.0067 l 2 ]l 2 cos2 B}l sinB )
对于1975国际椭球
{1 (C3 C5l 2 )l 2 cos2 B}l sin B
C3 0.33332 0.00678 cos2 B C5 0.2 cos2 B 0.0667
计算子午线收敛角的意义： 1、用于大地方位角和高斯平面坐标方位角的转换； 2、高斯正反算检核坐标计算的正确性。。

B Bf t3 f 2 4M f N
3 f
tf 2M f N f
y
2
5 3t
2 f
2 9 2 t 2 y 4 f f f
5、中央子午线收敛角和经度
2 y 2 2 y tan B f [1 (1 tan B f f )] 3 Nf 3N f
1 tan sin B l (1 t 2 3 2 2 4 ) sin B cos 2 Bl 3 3 1 (2 4t 2 2t 4 ) sin B cos 4 Bl 5 15

sin B l sin B cos 2 Bl 3 (1 3 2 2 4 )
L L0 l
小结
• • • • 了解反算公式的推导思路；掌握反算公式保符号的意义；用反算公式会进行计算；掌握子午线收敛角的定义及作用；

高斯投影及换带计算

测绘学院《大地测量学基础》课件
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）

1 高斯投影坐标正算公式

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

（2）投影变换必须满足的条件中央子午线投影后为直线；中央子午线投影后长度不变；投影具有正形性质，即正形投影条件。

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

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

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

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

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

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

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

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

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

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

大地测量高斯投影正反算

高斯投影正反算姓名：王义文班级：四班学号：2009301610135程序说明：本程序基于MFC基本对话框，由于MFC程序构建框架的代码冗长，打印出来恐怕影响阅读效率，所以下面的代码是经过剔除之后的一些核心代码，希望老师谅解！界面关联变量说明：纬度值的度、分、秒分别与B_DD,B_MM,B_SS相关联，缺省项自动设定为0。

经度值的度、分、秒分别与L_DD,L_MM,L_SS相关联，同上。

L0_DD则表示为L0的值，单位：度。

x,y分别与坐标值关联，且y为加500KM以后的值。

单位：米。

m_keshi=1是克氏椭球，=0为IAG椭球，且此处不能缺省。

m_zheng=1表示正算，=0为反算，且此处不能缺省。

核心代码说明：void C高斯投影正反算4Dlg::OnBnClickedRadio3() //此处为正算的设定，如果正算，设定输入焦{ 点在B_DD。

x,y为只读项。

m_zheng=1;((CEdit*)GetDlgItem(IDC_EDIT1))->SetFocus();((CEdit*)GetDlgItem(IDC_EDIT1))->SetReadOnly(FALSE);((CEdit*)GetDlgItem(IDC_EDIT2))->SetReadOnly(FALSE);((CEdit*)GetDlgItem(IDC_EDIT3))->SetReadOnly(FALSE);((CEdit*)GetDlgItem(IDC_EDIT4))->SetReadOnly(FALSE);((CEdit*)GetDlgItem(IDC_EDIT5))->SetReadOnly(FALSE);((CEdit*)GetDlgItem(IDC_EDIT6))->SetReadOnly(FALSE);((CEdit*)GetDlgItem(IDC_EDIT7))->SetReadOnly(FALSE);((CEdit*)GetDlgItem(IDC_EDIT8))->SetReadOnly();((CEdit*)GetDlgItem(IDC_EDIT9))->SetReadOnly();}void C高斯投影正反算4Dlg::OnBnClickedRadio4() //反算的设定：设置输入焦点在L0_DD，且B，{ L的角度值为只读。

高斯投影坐标计算

tf
tf
tf
1 1 2 2 3 l y ( 1 2 t ) y f f 3 N o sB 6 N o sB f c f fc f 1 2 2 4 2 2 5 ( 5 2 8 t 6 2 4 t 8 t ) y f f f f f 5 1 2 0 N c o s B f f
(3)距中央子午线愈远的子午线，投影后弯曲愈厉害，长度变形也愈大。
谢谢!
由恒等式两边对应系数相等，建立求解待定系数的递推公式
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=?
３) 由第二条件（中央子午线投影后长度不变）可知，位于中央子午线上的点，投影后的纵坐标 x 应该等于投影前从赤道量至该点的子午弧长。
高斯投影坐标计算
本节要点提要
1、高斯投影坐标正算公式 2、高斯投影坐标反算公式 3、高斯投影坐标正算的数值公式 4、高斯投影坐标反算的迭代计算公式
地图投影的分类
• 按投影变形性质分类: 等角投影等距投影等积投影
a=b
• 按投影面分类 : 圆锥面正轴投影切投影
a=1 or b=1
圆柱(椭圆柱) 面横轴投影割投影
d m d m d m 2 4 0 2 2 4 4 m 3 m l 5 m l l l 3 5 1 d q d q d q d m d m d m 3 5 3 3 1 2 m l 4 m l l l 5l 2 4 d q d q d q
d m d m d m 2 4 0 2 2 4 4 m 3 m l 5 m l l l 3 5 1 d q d q d q d m d m d m 3 3 3 5 5 1 2 m l 4 m l l l l 2 4 d q d q d q

高斯投影正反算公式

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

二、计算模型：基本椭球参数：椭球长半轴a椭球扁率f椭球短半轴：(1)b a f =-椭球第一偏心率：e a= 椭球第二偏心率：e b'=高斯投影正算公式：此公式换算的精度为0.001m6425644223422)5861(cos sin 720)495(cos 24cos sin 2l t t B B N l t B simB N l B B N X x ''+-''+''++-''+''⋅''+=ρηηρρ 5222425532233)5814185(cos 120)1(cos 6cos l t t t B N l t B N l B N y ''-++-''+''+-''+''⋅''=ηηρηρρ其中：角度都为弧度B 为点的纬度，0l L L ''=-，L 为点的经度，0L 为中央子午线经度； N 为子午圈曲率半径，1222(1sin )N a e B -=-；tan t B =； 222cos e B η'=1803600ρπ''=*其中X 为子午线弧长：2402464661616sin cos ()(2)sin sin 33X a B B B a a a a a B a B ⎡⎤=--++-+⎢⎥⎣⎦02468,,,,a a a a a 为基本常量，按如下公式计算：200468242684468686883535281612815722321637816323216128m a m m m m m m a m m m a m m m m a m a ⎧=++++⎪⎪⎪=+++⎪⎪⎪=++⎨⎪⎪=+⎪⎪⎪=⎪⎩02468,,,,m m m m m 为基本常量，按如下公式计算：22222020426486379(1);;5;;268m a e m e m m e m m e m m e m =-====；高斯投影反算公式：此公式换算的精度为0.0001’’.()()()()2222243246532235242225053922461904572012cos 6cos 5282468120cos f f f f f f f f f f f f f f f f f f f f f ff f f f f f ft t B B y t t yM N M N t y t t yM N y y l t N B N B y t t t N B L l L ηηηηη=-+++--++=-+++++++=+其中： 0L 为中央子午线经度。

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

2
x y , q l
x y l q
柯西-黎曼条件（公式）是
椭球面与平面之间的正形投影的一般条件
考虑到F=0，E=G，长度比公式简化为
x y q q E m2 2 = r r2
2
2 2
x y l l G m2 2 = r r2
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 dm5 l 5 2 4 dq dq dq
§4.9.2 正形投影的一般条件
一、长度比的通用公式推导
dS 2 ( MdB)2 ( N cos Bdl )2
M dB
ds 2 dx 2 dy 2
N cos B d l
长度比平方为：
dx 2 dy 2 ds 2 m 2 2 dS ( MdB) ( N cos Bdl )
？
上式为与方向有关的长度比的通用公式。上式在什么条件下与方向无关？
F 0
E G
柯西.黎曼条件（续）
正形条件：m与A无关，即满足： F 0
E G
2 2 2
x x y y 0 q l q l
y y x q l x l q
x y x y q q l l

高斯投影坐标正反算.ppt

2
4
6
8
高斯投影坐标正算（3）
dm0 = dX
dB =M
N cos B =N cos B
,
dq dB dq
M
c
m1 = N cos B
= V
cos B
子午线曲率半径
等量纬度定义式
m2

N 2
sin B cos B

m3

N b
cos3 B(1 t 2
2)

m4

N sin B cos3 24
y 2 dn4 dx
y4

N
cos M
B
(n1

3n3
y2

5n3
y4

)

2n2 y 4n4 y3
N cos B ( dn1 y dn3 y3 dn5 y5 )
M
dx
dx
dx
由恒等式两边对应系数相等，从而得待定系数的递推公式

n1

M N cos B
• 高斯平面直角坐标系：区分为：自然坐标；国家统一坐标。（掌握两者的换算）
§4.9.2 正形投影的一般条件
一、长度比的通用公式推导
dS 2 (MdB)2 ( N cos Bdl )2
ds2 dx2 dy2
长度比平方为：
m2

ds dS
2

dx2 dy2 (MdB)2 (N cos
dn2 dx
y2 dn4 dx
y4

N
cos M
B
(n1

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

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

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

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

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

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

2）由于高斯投影是换带投影，在每带内经差ｌ是不大的，ｌρ是一个微小量，所以可以将 X=X （ｌ，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.思路：级数展开，应用高斯投影三个条件，待定系数法求解。

「高斯投影坐标正反算公式及适合电算的高斯投影公式」

「高斯投影坐标正反算公式及适合电算的高斯投影公式」高斯投影坐标正反算公式是用于计算高斯投影坐标的数学公式。

高斯投影坐标是一种地理坐标系统，常用于测量和测绘工作中。

高斯投影坐标正算是指已知一个点的经纬度坐标，通过公式计算出该点的高斯投影坐标。

而高斯投影坐标反算是指已知一个点的高斯投影坐标，通过公式计算出该点的经纬度坐标。

一、高斯投影坐标正算公式：已知一个点的经纬度坐标(φ,λ)，其中φ为纬度，λ为经度，以及椭球体参数a、f和中央经线经度L0，可以通过以下步骤计算出该点的高斯投影坐标(X,Y)：1.计算扁率f'：f'=(a-b)/a其中，b=a*(1-f)是椭球体的短半轴。

2.计算黄赤交角ε：ε = atan(b / a)3.计算辅助量t：t = tan(π/4 - φ/2) / [(1 - f' * sin²φ)⁰.⁵ * (1 + e' *sinφ)⁰.⁵]其中，e'=f'*(2-f')是椭球体的第一偏心率。

4.计算辅助量η：η = e'^2 * cos²φ5.计算系数A、B、C和D：A = (L - L0) * cosφC = (L - L0) * cos⁵φ * (5 - tan²φ + 9e'^² + 4e'^⁴ - 24e'^² * tan²φ - 45e'^⁴ * tan²φ)D = (L - L0) * cos⁷φ * (61 - 58tan²φ + tan⁴φ + 270e'^² - 330e'^² * tan²φ)6.计算高斯坐标X和Y：X=k0*a*(A+B/2+C/4+D/6)Y=k0*a*(C/2+D/8)其中，k0是比例系数，一般情况下取1二、高斯投影坐标反算公式：已知一个点的高斯投影坐标(X,Y)，以及椭球体参数a、f、中央经线经度L0、比例系数k0和起始经度L1，可以通过以下步骤计算出该点的经纬度坐标(φ,λ)：1.计算扁率f'：f'=(a-b)/a其中，b=a*(1-f)是椭球体的短半轴。

高斯投影计算内容

高斯投影计算内容地面观测归算到椭球面上以后，可以通过两种途径获得各点的高斯平面直角坐标1：按第五章所述方法解算到球面三角形，推算个边大地线长度和大地方位角，解算出各点的大地坐标，然后利用高斯投影坐标计算公式解算出各点的高斯平面直角坐标。

2：将椭球面上的起算元素和观测元素归算到高斯投影平面上，然后解算平面三角形，推算各边的边长和坐标方位角，在平面上进行平差计算，求解各点的高斯平面直角坐标系。

两种方法比较两种方法解算出的结果是一样的，只是第一种计算量大。

所以通常采用第二种方法，将控制网直接归化到高斯投影平面上，在平面上进行平差和各种计算工作。

以三角网为例由于高斯投影为等角投影，所以椭球面上大地线之间的夹角与高斯平面上相应的投影曲线之间的夹角相等，但是各大地线的长度与它们在投影平面上的投影曲线长度并不相等，因为投影存在长度变形。

为了在平面上进行三角网的平差和计算必须把椭球面上以大地线构成的三角网，换算到高斯投影平面上以直线边构成的三角网。

为此应进行以下计算：（1）将起算点的大地坐标（B，L）换算成高斯投影平面上其投影的平面直角坐标（x,y），称为高斯投影坐标计算。

（2）将起算边的大地方位角A换算成投影平面上的平面坐标方位角T。

若由大地方位角转换平面坐标方位角，必须先计算该大地线起点处的平面子午线收敛角和方向改正。

T12=A12-γ+δ（3）将起算边的大地线长度归算到高斯平面的直线长度，其中，由大地线长度改化为高斯平面上直线长度时加入的改正数，成为距离改正D12=S12+∆S（4）对于椭球面上三角网的各观测方向和边长分别进行方向改正和距离改正，归算至高斯投影平面上的直线方向和直线距离，组成由平面三角形构成的三角网，进行平差计算，解算平面三角形，推求各控制点的平面直角坐标。

综上可知，高斯投影坐标计算、平面子午线收敛角计算、方向改正计算、距离改正计算，将此统称为高斯投影计算。

高斯投影正反算编程一.高斯投影正反算基本公式

高斯投影正反算编程一．高斯投影正反算基本公式（1）高斯正算基本公式（2）高斯反算基本公式以上主要通过大地测量学基础课程得到，这不进行详细的推导，只是列出基本公式指导编程的进行。

二．编程的基本方法和流程图（1）编程的基本方法高斯投影正反算基本上运用了所有的编程基本语句，本文中是利用C++语言进行基本的设计。

高斯正算中对椭球参数和带宽的选择主要运用了选择语句。

而高斯反算中除了选择语句的应用，在利用迭代算法求底点纬度还应用了循环语句。

编程中还应特别注意相关的度分秒和弧度之间的相互转换，这是极其重要的。

（2）相关流程图1）正算2）反算三．编程的相关代码（1）正算# include "stdio.h"# include "stdlib.h"# include "math.h"# include "assert.h"#define pi (4*atan(1.0))int i;struct jin{double B;double L;double L0;};struct jin g[100];main(int argc, double *argv[]){FILE *r=fopen("a.txt","r");assert(r!=NULL);FILE *w=fopen("b.txt","w");assert(r!=NULL);int i=0;while(fscanf(r,"%lf %lf %lf",&g[i].B,&g[i].L,&g[i].L0)!=EOF){double a,b;int zuobiao;printf("\n请输入坐标系：北京54=1，西安80=2，WGS84=3：");scanf("%d",&zuobiao);getchar();if(zuobiao==1){a=6378245;b=6356863.0187730473;}if(zuobiao==2){a=6378140;b=6356755.2881575287;}if(zuobiao==3){a=6378137;b=6356752.3142;} //选择坐标系//double f=(a-b)/a;double e,e2;e=sqrt(2*f-f*f);e2=sqrt((a/b)*(a/b)-1);//求椭球的第一，第二曲率//double m0,m2,m4,m6,m8;double a0,a2,a4,a6,a8;m0=a*(1-e*e);m2=3*e*e*m0/2;m4=5*e*e*m2/4;m6=7*e*e*m4/6;m8=9*e*e*m6/8;a0=m0+m2/2+3*m4/8+5*m6/16+35*m8/128;a2=m2/2+m4/2+15*m6/32+7*m8/16;a4=m4/8+3*m6/16+7*m8/32;a6=m6/32+m8/16;a8=m8/128;double Bmiao,Lmiao, L0miao;Bmiao=(int)(g[i].B)*3600.0+(int)((g[i].B-(int)(g[i].B))*100.0)*60.0+( g[i].B*100-(int)(g[i].B*100))*100.0;Lmiao=(int)(g[i].L)*3600.0+(int)((g[i].L-(int)(g[i].L))*100.0)*60.0+(g [i].L*100-(int)(g[i].L*100))*100.0;L0miao=(int)(g[i].L0)*3600.0+(int)((g[i].L0-(int)(g[i].L0))*100.0)*60. 0+(g[i].L0*100-(int)(g[i].L0*100))*100.0;double db;db=pi/180.0/3600.0;double B1,L1,l;B1=Bmiao*db;L1= Lmiao*db;l=L1-L0miao*db;//角度转化为弧度//double T=tan(B1)*tan(B1);double n=e2*e2*cos(B1)*cos(B1);double A=l*cos(B1);double X,x,y;X=a0*(B1)-a2*sin(2*B1)/2+a4*sin(4*B1)/4-a6*sin(6*B1)/6+a8*sin(8*B1)/8;//求弧长//double N=a/sqrt(1-e*e*sin(B1)*sin(B1));int Zonewide;int Zonenumber;printf("\n请输入带宽：3度带或6度带Zonewide=");scanf("%d",&Zonewide);getchar();if(Zonewide==3){Zonenumber=(int)((g[i].L-Zonewide/2)/Zonewide+1);}else if(Zonewide==6){Zonenumber=(int)g[i].L/Zonewide+1;}else{printf("错误");exit(0);}//选择带宽//doubleFE=Zonenumber*1000000+500000;//改写为国家通用坐标//y=FE+N*A+A*A*A*N*(1-T*T+n*n)/6+A*A*A*A*A*N*(5-18*T*T+T *T*T*T+14*n*n-58*n*n*T*T)/120;x=X+tan(B1)*N*A*A/2+tan(B1)*N*A*A*A*A*(5-T*T+9*n*n+4*n*n *n*n)/24+tan(B1)*N*A*A*A*A*A*A*(61-58*T*T+T*T*T*T)/720;printf("\n所选坐标系的转换结果：x=%lf y=%lf\n",x,y);fprintf(w,"%lf %lf\n",x,y);//输出结果到文本文件//}fclose(r);fclose(w);system("pause");return 0;}（2）反算# include "stdio.h"# include "stdlib.h"# include "math.h"# include "assert.h"#define pi (4*atan(1.0))double X,Y,B1,B2,B3,F,t;double m0,m2,m4,m6,m8;double a0,a2,a4,a6,a8,a1,b1;double BB,LL,Bf;double e,e1;int d,m,s,i,zuobiao;double sort(double,double);struct jin{double x;double y;double L0;};struct jin g[100];//x,y,L0为输入量：x,y坐标和中央子午线经度// main(int argc, double *argv[]){FILE *r=fopen("c.txt","r");assert(r!=NULL);FILE *w=fopen("d.txt","w");assert(r!=NULL);int i=0;while(fscanf(r,"%lf %lf %lf",&g[i].x,&g[i].y,&g[i].L0)!=EOF)//文件为空，无法打开//{double a1=6378245.0000000000;//克拉索夫斯基椭球参数//double b1=6356863.0187730473;double a75=6378140.0000000000;//1975国际椭球参数//double b75=6356755.2881575287;double a84=6378137.0000000000;//WGS-84系椭球参数//double b84=6356752.3142000000;double M,N;//mouyou圈曲率半径，子午圈曲率半径//double t,n;double A,B,C;double BB,LL,Bf,LL0,BB0;double a,b;printf("\n选择参考椭球：1=克拉索夫斯基椭球，2=1975国际椭球，3=WGS-84系椭球:");scanf("%d",&zuobiao);getchar();if(zuobiao==1){a=a1;b=b1;}if(zuobiao==2){a=a75;b=b75;}if(zuobiao==3){a=a84;b=b84;}//选择参考椭球，求解第一偏心率e,第二偏心率e1// Bf=sort(a,b);//调用求解底点纬度的函数//double q=sqrt(1-e*e*sin(Bf)*sin(Bf));double G=cos(Bf);M=a*(1-e*e)/(q*q*q);N=a/q;double H,I;A=g[i].y/N;H=A*A*A;I=A*A*A*A*A;t=tan(Bf);n=e1*cos(Bf);B=t*t;C=n*n;BB0=Bf-g[i].y*t*A/(2*M)+g[i].y*t*H/(24*M)*(5+3*B+C-9*B*C)-g[i] .y*t*I/(720*M)*(61+90*B+45*B*B);LL0=g[i].L0*pi/180.0+A/G-H/(6*G)*(1.0+2*B+C)+I/(120*G)*(5.0+28 *B+24*B*B+6*C+8*B*C);//利用公式求解经纬度//int Bdu,Bfen,Ldu,Lfen;double Bmiao,Lmiao;Ldu=int(LL0/pi*180);Lfen=int((LL0/pi*180)*60-Ldu*60);Lmiao=LL0/pi*180*3600-Ldu*3600-Lfen*60;Bdu=int(BB0/pi*180);Bfen=int((BB0/pi*180)*60-Bdu*60);Bmiao=BB0/pi*180*3600-Bdu*3600-Bfen*60;//将弧度转化为角度//printf("\n所选坐标系的转换结果：%d度%d分%lf秒%d 度%d分%lf秒\n",Bdu,Bfen,Bmiao,Ldu,Lfen,Lmiao);fprintf(w,"%d°%d’%lf”%d°%d’%lf”\n",Bdu,Bfen,Bmiao,Ldu,Lfen,Lmiao);//将结果输出到文本文件//}fclose(r);fclose(w);system("pause");return 0;}double sort(double a,double b){double e,e1;e=sqrt(1-(b/a)*(b/a));e1=sqrt((a/b)*(a/b)-1);double m0,m2,m4,m6,m8;double a0,a2,a4,a6,a8;m0=a*(1-e*e);m2=3*e*e*m0/2;m4=5*e*e*m2/4;m6=7*e*e*m4/6;m8=9*e*e*m6/8;a0=m0+m2/2+3*m4/8+5*m6/16+35*m8/128;a2=m2/2+m4/2+15*m6/32+7*m8/16;a4=m4/8+3*m6/16+7*m8/32;a6=m6/32+m8/16;a8=m8/128;B1=g[i].x/a0;do{F=-a2*sin(2*B1)/2+a4*sin(4*B1)/4-a6*sin(6*B1)/6+a8*sin(8*B1 )/8;B2=(g[i].x-F)/a0;B3=B1;B1=B2;} while(fabs(B3-B2)>10e-10);//利用迭代算法求解底点纬度//return B2; }。

[GB][21][高斯投影正反算实习]

4

UTM投影变形的特点：
UTM投影的中央经线长度比为0.999 6，这是为了使得
B＝00，l＝30处的最大变形值小于0.001而选择的数值。两条割线(在赤道上，它们位于离中央子午线大约 ±180km(约± 10 40’)处)上没有长度变形；离开这两条割线愈远变形愈大；在两条割线以内长度变形为负值；
q q q0 ln 0
1
B B B0 t1 ' q t2 ' q2 t3 ' q3 t4 ' q4 t5 ' q5 ,
23

方向改化及距离改化的简化公式：
1.2 " "
6R
2 0
( y2 y1 )(2 x1 x2 )
18

因为 d ( N cos B) dr dB M sin B N cos B dq dB dq M
d dq
sin B0
0 0
0 N0 cot B0 Ke sin B q 将上述两式代入微分方程得:
则有:
（q 0e
） 0 q
Ke
9
§4.11 兰勃脱投影概述 4.11.1兰勃脱投影基本概念兰勃脱(Lambert)投影是正形正轴圆锥投影。设想用一个圆锥套在地球椭球面上，使圆锥轴与椭球自转轴相一致，使圆锥面与椭球面一条纬线相切，将椭球面上的纬线投影到圆锥面上成为同心圆，经线投影圆锥面上成为从圆心发出的辐射直线，然后沿圆锥面某条母线(一般为中央经线L０)，将圆锥面切开而展成平面，从而实现了兰勃脱切圆锥投影。
25
26

兰勃脱投影是正形正轴圆锥投影，它的长度变形(m-1)与经度无关，但随纬差Δ Ｂ，即纵坐标x的增大而迅速增大，为限制长度变形，采用按纬度的分带投影，因此，这种投影适宜南北狭窄，东西延伸的国家和地区。这些国家根据本国实际情况，采用相应的分带方法和统一的坐标系统。但与高斯投影相比较，这种投影子午线收敛角有时过大，精密的方向改化和距离改化公式也较高斯投影要复杂，故目前国际上还是建议采用高斯投影。