高斯投影坐标正算C语言代码
#include
#include
#include
main()
{
printf("#####################################################\n");
printf("# 角度输入说明：如26°12′45.2″输入为26,12,45.2 #\n");
printf("#####################################################\n");
double b,B,L,L0,l1;
double d1,f1,m1,d2,f2,m2;
double x,y;
double p=206264.80625;
printf("请输入B：\n");
scanf("%lf,%lf,%lf",&d1,&f1,&m1); //维度数据输入
while (d1 < 0 || d1 >=90 || f1 < 0 || f1>= 60 || m1 < 0 || m1 >= 60)
{
printf("纬度输入格式不正确！\n");
printf("请重新输入B：\n");
scanf("%lf,%lf,%lf",&d1,&f1,&m1);
}
printf("您输入的大地维度为：\nB=%.0lf°%0.lf′%.2lf″\n",d1,f1,m1);//检查输入格式的正确性
printf("请输入L：\n");
scanf("%lf,%lf,%lf",&d2,&f2,&m2); //经度数据输入
while (d2 < 0 || d2 >=180 || f2 < 0 || f2>= 60 || m2 < 0 || m2 >= 60)
{
printf("经度输入格式不正确！\n");
printf("请重新输入L：\n");
scanf("%lf,%lf,%lf",&d2,&f2,&m2);
}
printf("您输入的大地经度为：\nB=%.0lf°%0.lf′%.2lf″\n",d2,f2,m2);//检查输入格式的正确性
B=d1*3600+f1*60+m1;
L=d2*3600+f2*60+m2;
//转化为以″为单位
b=B/p; //把角度转化为弧度制
int i,j,N;
printf("请选择带宽：1. 6°带。\n 2. 3°带。\n");
printf("选择的带宽为：");
scanf("%d",&i);
高斯投影坐标正反算 V B程序 文件编码（GHTU-UITID-GGBKT-POIU-WUUI-8968）
高斯投影坐标正反算 学院： 班级： 学号： 姓名： 课程名称: 指导老师：
实验目的： 1.了解高斯投影坐标正反算的基本思想； 2.学会编写高斯正反算程序，加深了解。 实验原理： 高斯投影正算公式中应满足的三个条件： 1. 中央子午线投影后为直线； 2. 中央子午线投影后长度不变； 3. 投影具有正形性质，即正形投影条件。 高斯投影反算公式中应满足的三个条件： 1. x坐标轴投影成中央子午线，是投影的对称轴； 2. x轴上的长度投影保持不变； 3. 正形投影条件，即高斯面上的角度投影到椭球面上后角度没有 变形，仍然相等。 操作工具： 计算机中的 代码： Dim a As Double, b As Double, x As Double, y As Double, y_#
Dim l_ As Double, b_ As Double, a0#, a2#, a4#, a6#, a8#, m2#, m4#, m6#, m8#, m0#, l0#, e#, e1# Dim deg1 As Double, min1 As Double, sec1 As Double, deg2 As Double, min2 As Double, sec2 As Double Private Sub Command1_Click() Dim x_ As Double, t#, eta#, N#, W#, k1#, k2#, ik1%, ik2%, dh% deg1 = Val min1 = Val sec1 = Val deg2 = Val min2 = Val sec2 = Val l_ = (deg1 * 3600 + min1 * 60 + sec1) / 206265 b_ = (deg2 * 3600 + min2 * 60 + sec2) / 206265 dh = Val k1 = ((l_ * 180 / + 3) / 6) k2 = (l_ * 180 / / 3) ik1 = Round(k1, 0) ik2 = Round(k2, 0) If dh = 6 Then l0 = 6 * ik1 - 3 Else
高斯投影正反算 班级：测绘九班C 姓名：塔娜 学号：2009301610323 指导老师：苏新洲 2011-11-02
高斯投影正反算编程 一、题目： 已知部分数据，根据高斯投影正反算思想进行编程，并采用克氏椭球，按3°或6°带投影。 正算：已知大地坐标B 、L， 二、已知数据： 正算： B=51.38439023 L=111.02131360 反算： x=5724004.723 y=19502559.920 三、计算结果： 正算结果： x=5724004.723 y=19502559.920 反算结果： B=51.38439023 L=111.02131360
（不予画出） 五、源代码： #include"gaosi.h" #include"math.h" #include"stdio.h" #include"tchar.h" #include"stdlib.h" #define pi 3.141592653589793 #define rho 206265 void Calculateellipse2plane(double B,double L); void Calculateplane2ellipse(double x,double y); double Dms2Rad(double Dms); double D2Dms(double D); double Dms2D(double Dms); int main() { double B=0,L=0; double x=0,y=0; int i=0; printf("如使用高斯投影坐标正算，请输入1；反算，请输入2\n"); scanf_s("%d",&i); if(i==1)
高斯投影坐标正反算编程报告 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. 计算模型 高斯投影正算公式 6 4256 4 42234 22)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 ''+-''+ ''++-''+''?''+=ρηηρρ 5 2224255 3 2233 )5814185(cos 120)1(cos 6cos l t t t B N l t B N l B N y ''-++-''+ ''+-''+''?''=ηηρηρρ 高斯投影反算公式 () () ()( ) 2 22425 52 23 36 4254 222232 8624285cos 12021cos 6cos 459061720935242f f f f f f f f f f f f f f f f f f f f f f f f f f f f f t t t B N y t B N y B N y l y t t y N M t y t t N M t y N M t B B ηηηηη+++++++-=++- -+++ -= 3. 程序框图
高斯投影坐标正反算 一、相关概念 大地坐标系由大地基准面和地图投影确定，由地图投影到特定椭圆柱面后在南北两极剪开展开而成，是对地球表面的逼近，各国或地区有各自的大地基准面，我国目前主要采用的基准面为：基准面，为GPS基准面，17届国际大地测量协会上推荐，椭圆柱长半轴a=6378137m，短半轴b=； 2.西安80坐标系，1975年国际大地测量协会上推荐，椭圆柱长半轴a=6378140m，短半轴b=; 3.北京54坐标系，参照前苏联克拉索夫斯基椭球体建立,椭圆柱长半轴a=6378245m, 短半轴b=; 通常所说的高斯投影有三种，即投影后： a)角度不变（正角投影），投影后经线和纬线仍然垂直； b)长度不变； c)面积不变； 大地坐标一般采用高斯正角投影，即在地球球心放一点光源，地图投影到过与中央经线相切的椭圆柱面上而成；可分带投影，按中央经线经度值分带，有每6度一带或每3度一带两种(起始带中央经线经度为均为3度，即：6度带1带位置0-6度，3度带1带位置度)，即所谓的高斯-克吕格投影。
图表11高斯投影和分带 地球某点经度（L）为过该点和地球自转轴的半圆与子午线所在半圆夹角，东半球为东经，西半球为西经；地球某点纬度（B）为所在水平面法线与赤道圆面的线面角。 正算是已知大地坐标（L，B），求解高斯平面坐标（X，Y）,为确保Y值为正，Y增加500公里；反算则是由高斯平面坐标（X，Y）求解大地坐标（L，B）。 二、计算模型： 地球椭球面由椭圆绕地球自转轴旋转180度而成。 图表 1 椭圆 椭圆长半轴a，椭圆短半轴b, 椭圆方程：
(1) 图表2椭球面 椭球面方程： y2 a2+ x2 b2 + z2 a2 =1 /*************************************** 与网上充斥的将函数关系先展开为泰勒级数，再依据投影规则确定各参数不同，本文直接依据空间立体三角函数关系得出结果。 *****/ （一）正算 由图表1，
#include "stdafx.h" #include "iostream.h" #include "math.h" #include "stdio.h" #define P 206264.806247096355 #define PI 3.141592653589793 void GaosZ_fun() { printf("高斯投影的正算\n"); double l,L,B,n2,x,y,N,t,V,c,e2; double i,j,k,n,h,a0,a4,a6,a3,a5,cB2; int m; e2=0.006738525414683; c=6399698.901782711; B=17.33557339*3600/P; L=119.15521159*3600/P; l=L-111*3600/P // l=((m%6)*3600+n*60+h)/P; t=tan(B); n2=e2*cos(B)*cos(B); V=sqrt(1+n2); cB2=pow(cos(B),2); N=6399698.902-(21562.267-(108.973-0.612*cB2)*cB2)*cB2; // N=c/V; a0=32140.404-(135.3302-(0.7092-0.004*cB2)*cB2)*cB2; a4=(0.25+0.00252*cB2)*cB2-0.04166; a6=(0.166*cB2-0.084)*cB2; a3=(0.3333333+0.001123*cB2)*cB2-0.1666667; a5=0.0083-(0.1667-(0.1968+0.0040*cB2)*cB2)*cB2; // x=X+N*sin(B)*cos(B)*l*l/2+N*sin(B)*pow(cos(B),3)*(5-t*t+9*n2+4*n2*n2)*pow(l, 4)/24+N*sin(B)*pow(cos(B),5)*(61-58*t*t+pow(t,4))*pow(l,6)/720; // y=N*cos(B)*l+N*pow(cos(B),3)*(1-t*t+n2)*pow(l,5)/6+N*pow(cos(B),5)*(5-18*t*t +pow(t,4)+14*n2-58*n2*t*t)*pow(l,5)/120; x=6367558.4969*B-(a0-(0.5+(a4+a6*l*l)*l*l)*l*l*N)*sin(B)*cos(B); y=(1+(a3+a5*l*l)*l*l)*l*N*cos(B); printf("x=%f\ny=%f\n",x,y); } void GaosF_fun() { printf("高斯投影的反算\n"); double B,Bf,Nf,b,b2,b3,b4,b5,Z,x,y,L0,l;
高斯投影坐标正反算 学院： 班级： 学号： 姓名： 课程名称: 指导老师：
实验目的： 1.了解高斯投影坐标正反算的基本思想； 2.学会编写高斯正反算程序，加深了解。 实验原理： 高斯投影正算公式中应满足的三个条件： 1. 中央子午线投影后为直线； 2. 中央子午线投影后长度不变； 3. 投影具有正形性质，即正形投影条件。 高斯投影反算公式中应满足的三个条件： 1. x坐标轴投影成中央子午线，是投影的对称轴； 2. x轴上的长度投影保持不变； 3. 正形投影条件，即高斯面上的角度投影到椭球面上后角度没 有变形，仍然相等。 操作工具： 计算机中的VB6.0 代码： Dim a As Double, b As Double, x As Double, y As Double, y_# Dim l_ As Double, b_ As Double, a0#, a2#, a4#, a6#, a8#, m2#, m4#,
m6#, m8#, m0#, l0#, e#, e1# Dim deg1 As Double, min1 As Double, sec1 As Double, deg2 As Double, min2 As Double, sec2 As Double Private Sub Command1_Click() Dim x_ As Double, t#, eta#, N#, W#, k1#, k2#, ik1%, ik2%, dh% deg1 = Val(Text1.Text) min1 = Val(Text2.Text) sec1 = Val(Text3.Text) deg2 = Val(Text4.Text) min2 = Val(Text5.Text) sec2 = Val(Text6.Text) l_ = (deg1 * 3600 + min1 * 60 + sec1) / 206265 b_ = (deg2 * 3600 + min2 * 60 + sec2) / 206265 dh = Val(Text9.Text) k1 = ((l_ * 180 / 3.14159 + 3) / 6) k2 = (l_ * 180 / 3.14159 / 3) ik1 = Round(k1, 0) ik2 = Round(k2, 0) If dh = 6 Then l0 = 6 * ik1 - 3 Else If dh = 3 Then
高斯投影正反算编程一．高斯投影正反算基本公式（1）高斯正算基本公式 （2）高斯反算基本公式
以上主要通过大地测量学基础课程得到，这不进行详细的推导，只是列出基本公式指导编程的进行。 二．编程的基本方法和流程图 （1）编程的基本方法 高斯投影正反算基本上运用了所有的编程基本语句，本文中是利用C++语言进行基本的设计。高斯正算中对椭球参数和带宽的选择主要运用了选择语句。而高斯反算中除了选择语句的应用，在利用迭代算法求底点纬度还应用了循环语句。编程中还应特别注意相关的度分秒和弧度之间的相互转换，这是极其重要的。 （2）相关流程图 1）正算
选择带宽 3/6 度带 计算带号 输入大地坐标 B ，L 和经差 L0 6 度带 3 度带 选择椭球参数 计算带号 计算弧长 计算平面坐标 x,y 打印 x,y 开始 计算平面坐标 x,y 计算弧长 打印 x,y
开始 输入自然值坐标x,y 和经差L0 选择椭球参数 利用迭代算法 求解底点纬度 利用公式计算B 和L 打印B 和L 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;
高斯投影坐标正算公式 高斯投影坐标正反算公式 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式得
高斯投影坐标正反算编程报告 10021班 张鑫 学号：2010302590040 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. 计算模型 高斯投影正算公式 6 4256 4 42234 22)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 ''+-''+ ''++-''+''?''+=ρηηρρ 5 2224255 3 2233)5814185(cos 120)1(cos 6cos l t t t B N l t B N l B N y ''-++-''+''+-''+''?''=ηηρηρρ 高斯投影反算公式 () () ()( ) 222425 52 2336 4254 2222 32 8624285cos 12021cos 6cos 459061720935242f f f f f f f f f f f f f f f f f f f f f f f f f f f f f t t t B N y t B N y B N y l y t t y N M t y t t N M t y N M t B B ηηηηη+++++++-=++- -+++ -=
高斯投影正、反算代码 //高斯投影正、反算 //////6度带宽 54年北京坐标系 //高斯投影由经纬度(Unit:DD)反算大地坐标(含带号，Unit:Metres) void GaussProjCal(double longitude, double latitude, double *X, double *Y) { int ProjNo=0; int ZoneWide; ////带宽 double longitude1,latitude1, longitude0,latitude0, X0,Y0, xval,yval; double a,f, e2,ee, NN, T,C,A, M, iPI; iPI = 0.0174532925199433; ////3.1415926535898/180.0; ZoneWide = 6; ////6度带宽 a=6378245.0; f=1.0/298.3; //54年北京坐标系参数 ////a=6378140.0; f=1/298.257; //80年西安坐标系参数 ProjNo = (int)(longitude / ZoneWide) ; longitude0 = ProjNo * ZoneWide + ZoneWide / 2; longitude0 = longitude0 * iPI ; latitude0=0; longitude1 = longitude * iPI ; //经度转换为弧度 latitude1 = latitude * iPI ; //纬度转换为弧度 e2=2*f-f*f; ee=e2*(1.0-e2); NN=a/sqrt(1.0-e2*sin(latitude1)*sin(latitude1)); T=tan(latitude1)*tan(latitude1); C=ee*cos(latitude1)*cos(latitude1); A=(longitude1-longitude0)*cos(latitude1); M=a*((1-e2/4-3*e2*e2/64-5*e2*e2*e2/256)*latitude1-(3*e2/8+3*e2*e2 /32+45*e2*e2 *e2/1024)*sin(2*latitude1)
高斯投影坐标正反算V B 程序 Jenny was compiled in January 2021
高斯投影坐标正反算 学院： 班级： 学号： 姓名： 课程名称: 指导老师：
实验目的： 1.了解高斯投影坐标正反算的基本思想； 2.学会编写高斯正反算程序，加深了解。 实验原理： 高斯投影正算公式中应满足的三个条件： 1. 中央子午线投影后为直线； 2. 中央子午线投影后长度不变； 3. 投影具有正形性质，即正形投影条件。 高斯投影反算公式中应满足的三个条件： 1. x坐标轴投影成中央子午线，是投影的对称轴； 2. x轴上的长度投影保持不变； 3. 正形投影条件，即高斯面上的角度投影到椭球面上后角度没有 变形，仍然相等。 操作工具： 计算机中的 代码： Dim a As Double, b As Double, x As Double, y As Double, y_#
Dim l_ As Double, b_ As Double, a0#, a2#, a4#, a6#, a8#, m2#, m4#, m6#, m8#, m0#, l0#, e#, e1# Dim deg1 As Double, min1 As Double, sec1 As Double, deg2 As Double, min2 As Double, sec2 As Double Private Sub Command1_Click() Dim x_ As Double, t#, eta#, N#, W#, k1#, k2#, ik1%, ik2%, dh% deg1 = Val min1 = Val sec1 = Val deg2 = Val min2 = Val sec2 = Val l_ = (deg1 * 3600 + min1 * 60 + sec1) / 206265 b_ = (deg2 * 3600 + min2 * 60 + sec2) / 206265 dh = Val k1 = ((l_ * 180 / + 3) / 6) k2 = (l_ * 180 / / 3) ik1 = Round(k1, 0) ik2 = Round(k2, 0) If dh = 6 Then l0 = 6 * ik1 - 3 Else
高斯投影正反算程序设计 一．程序设计流程 本程序的设计思路如下： （1），程序采用VS08版本作为开发平台，并采用C#语言作为开发语言，设计为WindowsForm 窗体程序形式。 （2），程序主要的算法来自于教材。但是本程序为了更加实用，添加了更多的解算基准，包括：WGS-84，国际椭球1975，克氏椭球，和2000国家大地坐标系。 （3），程序为了更方便的读取数据和输出数据，故需要自己定义了固定的数据输入格式和数据输出格式或形式，请老师注意查看。 二．代码 using System; using Syst using https://www.doczj.com/doc/ec14895571.html,ponentModel; using System.Data; using System.Drawing; using System.Text; namespace Gauss { public partial class Form1 : Form { //大地坐标 //Geodetic Coordinate public struct CRDGEODETIC { public double dLongitude; public double dLatitude; public double dHeight; } //笛卡尔坐标
//Cartesian Coordinate public struct CRDCARTESIAN { public double x; public double y; public double z; } public Form1() { InitializeComponent(); } private void button1_Click(object sender, EventArgs e) { double ee = 0; double a = 0; string tt; try { } catch { MessageBox.Show("Gauss Inverse: Choose datum error!"); return; } if (https://www.doczj.com/doc/ec14895571.html,pareTo("克氏椭球")==0) { a = 6378245.00; } if (https://www.doczj.com/doc/ec14895571.html,pareTo("WGS-84") == 0) { a = 6378137.00; } if (https://www.doczj.com/doc/ec14895571.html,pareTo("1975国际椭球") == 0)
高斯投影正反算编程 一、题目： 已知部分数据，根据高斯投影正反算思想进行编程，并采用克氏椭球，按3°或6°带投影。 正算：已知大地坐标B 、L， 二、已知数据： 正算： B=51.38439023 L=111.02131360 反算： x=5724004.723 y=19502559.920 三、计算结果： 正算结果： x=5724004.723 y=19502559.920 反算结果： B=51.38439023 L=111.02131360
五、源代码： #include"gaosi.h" #include"math.h" #include"stdio.h" #include"tchar.h" #include"stdlib.h" #define pi 3.141592653589793 #define rho 206265 void Calculateellipse2plane(double B,double L); void Calculateplane2ellipse(double x,double y); double Dms2Rad(double Dms); double D2Dms(double D); double Dms2D(double Dms);
int main() { double B=0,L=0; double x=0,y=0; int i=0; printf("如使用高斯投影坐标正算，请输入1；反算，请输入2\n"); scanf_s("%d",&i); if(i==1) { printf("请输入大地坐标纬度B(度分秒):\n"); scanf_s("%lf",&B); printf("请输入大地坐标经度L(度分秒):\n"); scanf_s("%lf",&L); Calculateellipse2plane(B,L); } if(i==2) { printf("请输入国家统一坐标x(m):\n"); scanf_s("%lf",&x); printf("请输入国家统一坐标y(m):\n"); scanf_s("%lf",&y); Calculateplane2ellipse(x,y); } return 0; } void Calculateellipse2plane(double B,double L) //高斯投影正算主体{ double l=0,Lo=0,a0=0,a3=0,a4=0,a5=0,a6=0,n=0,c=0; double x=0,y=0; double m=0,p=0,q=0; int N=0,i=0; //带号 printf("如使用6°带请输入1，使用3°带请输入2\n"); scanf_s("%d",&i); if(i==1) //已知a点在6°带的带号和中央子午线经度 { N=int(Dms2D(L)/6); Lo=6*N-3; } if(i==2) //已知a点在3°带的带号和中央子午线经度 { N=int((Dms2D(L)+1.5)/3);
§8.3高斯投影坐标正反算公式 任何一种投影①坐标对应关系是最主要的；②如果是正形投影，除了满足正形投影的条件外（C-R 偏微分方程），还有它本身的特殊条件。 8.3.1高斯投影坐标正算公式： B,l ? x,y 高斯投影必须满足以下三个条件： ①中央子午线投影后为直线；②中央子午线投影后长度不变；③投影具有正形性质，即正形投影条件。 由第一条件知中央子午线东西两侧的投影必然对称于中央子午线，即(8-10)式中，x 为l 的偶函数，y 为l 的奇函数；0330'≤l ，即20/1/≈''''ρl ，如展开为l 的级数，收敛。 +++=++++=553316644220l m l m l m y l m l m l m m x （8-33） 式中 ,,10m m 是待定系数，它们都是纬度B 的函数。 由第三个条件知： q y l x l y q x ??-=????=??, (8-33)式分别对l 和q 求偏导数并代入上式 ----=++++++=+++553315 63424 42204 52 3164253l dq dm l dq dm l dq dm l m l m l m l dq dm l dq dm dq dm l m l m m (8-34) 上两式两边相等，其必要充分条件是同次幂l 前的系数应相等，即
dq dm m dq dm m dq dm m 231 20 13121? =? -== (8-35) (8-35)是一种递推公式，只要确定了 0m 就可依次确定其余各系数。 由第二条件知:位于中央子午线上的点，投影后的纵坐标x 应等于投影前从赤道量至该点的子午线弧长X ，即(8-33)式第一式中，当0=l 时有： 0m X x == (8-36) 顾及(对于中央子午线) B V M r M B N dq dB M dB dX cos cos 2 ==== 得： B V c B N r dq dB dB dX dq dX dq dm m cos cos 01===?===(8-37,38) B B N dq dB dB dm dq dm m cos sin 2 2121112=?-=?-= (8-39) 依次求得6543,,,m m m m 并代入(8-33)式，得到高斯投影正算公式
高斯投影坐标正反算公式 未知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)
高斯投影坐标反算c语言代码 #include
高斯投影正反算c代码 Coca-cola standardization office【ZZ5AB-ZZSYT-ZZ2C-ZZ682T-ZZT18】
高斯投影正反算程序设计 一．程序设计流程 本程序的设计思路如下： （1），程序采用VS08版本作为开发平台，并采用C#语言作为开发语言，设计为WindowsForm窗体程序形式。 （2），程序主要的算法来自于教材。但是本程序为了更加实用，添加了更多的解算基准，包括：WGS-84，国际椭球1975，克氏椭球，和2000国家大地坐标系。 （3），程序为了更方便的读取数据和输出数据，故需要自己定义了固定的数据输入格式和数据输出格式或形式，请老师注意查看。 二．代码 using System; using ; using ; using ; using ; using Gauss { public partial class Form1 : Form { double b = (a * a * (1 - ee * ee)); double c = a * a / b; double epp = ((a * a - b * b) / b / b); CRDGEODETIC pcrdGeo; CRDCARTESIAN pcrdCar; double midlong = 0;
//求X，Y和带号 = ; ytext = ; string temp = (0, 2); num = (temp); ytext = (0, 2); = (ytext) - 500000; try { tt = } catch { ("Choose 3/6 error!"); return; } if ("3度带") == 0) { midlong = num * 3 * pai / 180; } if ("6度带") == 0) { midlong = (6 * num - 3) * pai / 180; } b = (a * a * (1 - ee * ee)); c = a * a / b; epp = (a * a - b * b) / b; double m0, m2, m4, m6, m8; double a0, a2, a4, a6, a8; m0 = a * (1 - ee * ee);
高斯投影是正形投影的一种，同一坐标系中的高斯投影换带计算公式是根据正形投影原理推导出的两个高斯坐标系间的显函数式。在同一大地坐标系中(例如1954北京坐标系或1980西安坐标系)，如果两个高斯坐标系只是主子午线的经度不同，那么显函数式前的系数可以根据坐标系使用的椭球元素和主子午线经度唯一确定。但如果两个高斯坐标系除了主子午线的经度不同以外，还存在其他线性系，则将线性变换公式代入换带计算的显函数式中，仍然可以得到严密的坐标变换公式。此时显函数式前的系数等价于使用两个坐标系主子午线的经度和线性变换参数联合求解得到的，可以唯一确定。 //6度带宽54北京坐标系 //高斯投影由大地坐标(Unit:Metres)反算经纬度(Unit:DD) void GaussProjInvCal(double X, double Y, double *longitude, double *latitude) { int ProjNo; int ZoneWide; ////带宽 double longitude1,latitude1, longitude0,latitude0, X0,Y0, xval,yval; double e1,e2,f,a, ee, NN, T,C, M, D,R,u,fai, iPI; iPI = 0.0174532925199433; ////3.1415926535898/180.0; a = 6378245.0; f = 1.0/298.3; //54年北京坐标系参数 ////a=6378140.0; f=1/298.257; //80年西安坐标系参数 ZoneWide = 6; ////6度带宽 ProjNo = (int)(X/1000000L) ; //查找带号 longitude0 = (ProjNo-1) * ZoneWide + ZoneWide / 2; longitude0 = longitude0 * iPI ; //中央经线 X0 = ProjNo*1000000L+500000L; Y0 = 0; xval = X-X0; yval = Y-Y0; //带内大地坐标
经纬度BL换算到高斯平面直角坐标xy（高斯投影正算） private function bl2xy(byref a2 as double, byref f2 as double, byref e2 as double, _ byref s2 as double, byref t2 as double as boolean 'a2 输入中央子午线，以度.分形式输入，如115度30分则输入115.30; 起算数据l0 'f2 以度小数形式输入经度值, l 'e2 以度小数形式输入纬度值，b 's2 计算结果，横坐标y 't2 计算结果，纵坐标x '投影带号计算n=[l/6]+1 如：测得经度103.xxxx，故n=[103.x/6]+1=17+1=18 '中央经线经度l0 = n*6-3 = [l/6]*6+3 dim b2 as double 'dim g2 as double dim h2 as double dim i2 as double dim j2 as double dim k2 as double dim l2 as double dim m2 as double dim n2 as double dim o2 as double dim p2 as double dim q2 as double
dim r2 as double b2 = int(a2 + (int(a2 * 100 - int(a2 * 100 / 60 + (a2 * 10000 - int(a2 * 100 * 100 / 3600 '把l0化成度(a2 'g2 = f2 - b2 ' l -l0 'h2 = g2 / 57.2957795130823 '化作弧度 h2 = (f2 - b2 / 57.2957795130823 '将经差的单位化为弧度 i2 = tan(e2 / 57.2957795130823 'tan (b j2 = cos(e2 / 57.2957795130823 ' cos (b k2 = 0.006738525415 * j2 * j2 l2 = i2 * i2 m2 = 1 + k2 n2 = 6399698.9018 / sqr(m2 o2 = h2 * h2 * j2 * j2 p2 = i2 * j2 q2 = p2 * p2 r2 = (32005.78006 + q2 * (133.92133 + q2 * 0.7031 s2 = ((((l2 - 18 * l2 - (58 * l2 - 14 * k2 + 5 * o2 / 20 + m2 - l2 * o2 / 6 + 1 * n2 * (h2 * j2 s2 = s2 + 18500000 '在计算的基础上加上了“带号”（18）和“东移”（500km） '计算结果,横坐标y t2 = 6367558.49686 * e2 / 57.29577951308 - p2 * j2 * r2 + ((((l2 - 58 * l2 + 61 * _ o2 / 30 + (4 * k2 + 5 * m2 - l2 * o2 / 12 + 1 * n2 * i2 * o2 / 2 '计算结果，纵坐标x