高斯投影正反算

C=6,399,698.901782710000000e'2=0.006738525414700L=0.000000000000000B=0.000000000000000L0=X=0.000000000000000m=0.000000000000000t=η2=0.006738525414700N=6,378,244.999999950000000l=x=0.000000000000000y=0.000000000000000作者:张传仁(AoXiangZCR@)使用说明:计算表中的"L0","L"及"B"为大地坐标的经纬度值,单位为十进制度.计算出的横坐标不含加常数 椭球参数已知数据项计算参数项待求数据项

0.0000000000000000.0000000000000000.000000000000000

度.计算出的横坐标不含加常数500KM.请在此栏输入数据

高斯投影坐标正反算

- 1 -

高斯投影坐标正反算

学 院：

班 级：

学 号：

姓 名：

课程名称:

指导老师：

高斯投影坐标正反算

- 2 - 实验目的：

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#, 高斯投影坐标正反算

- 3 - 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%

* *

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

（1）高斯正算基本公式

（2）高斯反算基本公式

* *

以上主要通过大地测量学基础课程得到，这不进行详细的推导，只是列出基本公式指导编程的进行。

二．编程的基本方法和流程图

（1）编程的基本方法

高斯投影正反算基本上运用了所有的编程基本语句，本文中是利用C++语言进行基本的设计。高斯正算中对椭球参数和带宽的选择主要运用了选择语句。而高斯反算中除了选择语句的应用，在利用迭代算法求底点纬度还应用了循环语句。编程中还应特别注意相关的度分秒和弧度之间的相互转* *

换，这是极其重要的。

（2）相关流程图

1）正算

* *

输入大地坐标B，L

和经差L0

选择带宽3/6度带

计算带号

计算弧长

计算平面坐标x,y

打印x,y 计算带号

计算弧长

计算平面坐标x,y

打印x,y

开始

选择椭球参数

3度带 6度带 * *

2）反算

开始

输入自然值坐标x,y

和经差L0

利用迭代算法

求解底点纬度

利用公式计算B和L

打印B和L 选择椭球参数 * *

三．编程的相关代码

（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); * *

高斯投影正反算程序设计

一．程序设计流程

本程序的设计思路如下：

（1），程序采用VS08版本作为开发平台，并采用C#语言作为开发语言，设计为WindowsForm窗体程序形式。

（2），程序主要的算法来自于教材。但是本程序为了更加实用，添加了更多的解算基准，包括：WGS-84，国际椭球1975，克氏椭球，和2000国家大地坐标系。

（3），程序为了更方便的读取数据和输出数据，故需要自己定义了固定的数据输入格式和数据输出格式或形式，请老师注意查看。

二．代码

using System;

using Syst

using 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()