资源预览内容
第1页 / 共11页
第2页 / 共11页
第3页 / 共11页
第4页 / 共11页
第5页 / 共11页
第6页 / 共11页
第7页 / 共11页
第8页 / 共11页
第9页 / 共11页
第10页 / 共11页
亲,该文档总共11页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述
高斯投影正反算程序设计一程序设计流程本程序的设计思路如下:(1),程序采用VS08版本作为开发平台,并采用C#语言作为开发语言,设计为WindowsForm窗体程序形式。(2),程序主要的算法来自于教材。但是本程序为了更加实用,添加了更多的解算基准,包括:WGS-84,国际椭球1975,克氏椭球,和2000国家大地坐标系。(3),程序为了更方便的读取数据和输出数据,故需要自己定义了固定的数据输入格式和数据输出格式或形式,请老师注意查看。二代码using System;using System.Collections.Generic;using System.ComponentModel;using System.Data;using System.Drawing;using System.Text;using System.Windows.Forms;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 tt = this.comboBox1.Itemsthis.comboBox1.SelectedIndex.ToString(); catch MessageBox.Show(Gauss Inverse: Choose datum error!); return; if (tt.CompareTo(克氏椭球)=0) a = 6378245.00; ee = Math.Sqrt(0.006693421622); if (tt.CompareTo(WGS-84) = 0) a = 6378137.00; ee = Math.Sqrt(0.00669437999013); if (tt.CompareTo(1975国际椭球) = 0) a = 6378140.00; ee = Math.Sqrt(0.006694384999588); if (tt.CompareTo(2000国家大地坐标系) = 0) a = 6378137.0; ee =Math.Sqrt(0.0066943802290); const double pai = 3.1415926; double b = Math.Sqrt(a * a * (1 - ee * ee); double c = a * a / b; double epp = Math.Sqrt(a * a - b * b) / b / b); CRDGEODETIC pcrdGeo; CRDCARTESIAN pcrdCar; double midlong; /求纬度 string temp; temp = textBox1.Text.Split( ); double tempradius = new double3; for (int i = 0; i 3; i+) tempradiusi = Convert.ToDouble(tempi); pcrdGeo.dLatitude = tempradius0 / 180.0 * pai + tempradius1 / 180.0 / 60.0 * pai + tempradius2 / 180 / 60.0 / 60 * pai; /求经度 temp = textBox2.Text.Split( ); for (int i = 0; i 3; i+) tempradiusi = Convert.ToDouble(tempi); pcrdGeo.dLongitude = tempradius0 / 180.0 * pai + tempradius1 / 180.0 / 60.0 * pai + tempradius2 / 180 / 60.0 / 60 * pai; int deglon = Convert.ToInt32(pcrdGeo.dLongitude * 180 / pai); /求中央经度 int num; /带号 midlong = 0; /默认值,需要制定分带 try tt = this.comboBox3.Itemsthis.comboBox3.SelectedIndex.ToString(); catch MessageBox.Show(Choose 3/6 error!); return; if (tt.CompareTo(3度带) = 0) num = Convert.ToInt32(deglon / 6 + 1); midlong = (6 * num - 3) / 180.0 * pai; if (tt.CompareTo(6度带) = 0) num = Convert.ToInt32(deglon + 1.5) / 3); midlong = num * 3 * pai / 180; double lp=pcrdGeo.dLongitude - midlong; double N = c / Math.Sqrt(1 + epp * epp * Math.Cos(pcrdGeo.dLatitude) * Math.Cos(pcrdGeo.dLatitude); double M = c / Math.Pow(1 + epp * epp * Math.Cos(pcrdGeo.dLatitude) * Math.Cos(pcrdGeo.dLatitude), 1.5); double ita = epp * Math.Cos(pcrdGeo.dLatitude); double t = Math.Tan(pcrdGeo.dLatitude); double Nscnb = N * Math.Sin(pcrdGeo.dLatitude) * Math.Cos(pcrdGeo.dLatitude); double Ncosb = N * Math.Cos(pcrdGeo.dLatitude); double cosb = Math.Cos(pcrdGeo.dLatitude); double X; double m0, m2, m4, m6, m8; double a0, a2, a4, a6, a8; m0 = a * (1 - ee * ee); m2 = 3.0 / 2.0 * m0 * ee * ee; m4 = 5.0 / 4.0 * ee * ee * m2; m6 = 7.0 / 6.0 * ee * ee * m4; m8 = 9.0 / 8.0 * ee * ee * m6; a0 = m0 + m2 / 2.0 + 3.0 / 8.0 * m4 + 5.0 / 16.0 * m6
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号