大地测量 高斯投影正反算

《大地测量 高斯投影正反算》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、大地测量高斯投影正反算程序代码课程：大地测量学基础姓名：林江伟学号：2008301610045班级：0804界面如下：输入数据计算：usingSystem;usingSystem.Collections.Generic;usingSystem.ComponentModel;usingSystem.Data;usingSystem.Drawing;usingSystem.Text;usingSystem.Windows.Forms;namespace大地{publicpartialclassForm1:Form{doubleB,L;do

2、ublex,y;doubleX,Y;intN;doubleL0,l1;doublep=206264.80625;publicForm1(){InitializeComponent();}privatevoidForm1_Load(objectsender,EventArgse){//自动化控件显示初始值radioButton2.Checked=true;radioButton3.Checked=true;textBox1.Focus();textBox1.Text=textBox2.Text=textBox3.Text=textBox

3、4.Text=textBox5.Text=textBox6.Text="0";this.richTextBox2.Text="说明：输入的坐标需为按6°带投影且采用克氏椭球参数所得的国家统一坐标";}privatevoidbutton1_Click(objectsender,EventArgse){//获取输入数据doublebb1=Convert.ToDouble(this.textBox1.Text);doublebb2=Convert.ToDouble(this.textBox2.Text);doublebb3=Convert.

4、ToDouble(this.textBox3.Text);doublell1=Convert.ToDouble(this.textBox4.Text);doublell2=Convert.ToDouble(this.textBox5.Text);doublell3=Convert.ToDouble(this.textBox6.Text);//检查输入格式的正确性if(bb1>=0&&bb1<90&&bb2>=0&&bb2<60&&bb3>=0&&bb3<60){B=bb1*3600+bb2*60+bb3;}else{MessageBo

5、x.Show("纬度输入格式不正确！","警告");return;}if(ll1>=0&&ll1<360&&ll2>=0&&ll2<60&&ll3>=0&&ll3<60){L=ll1*3600+ll2*60+ll3;}else{MessageBox.Show("经度输入格式不正确！","警告");return;}doubleb=B/p;//获取用户选项值inttyd=radioButton1.Checked?3:6;inttq=radioButton3.Checked?1:2;//按带求带号以及中央子午线经度if(tyd==6){N=

6、(int)(L/(6*3600.0))+1;L0=(6*N-3)*3600.0;l1=(L-L0)/p;}if(tyd==3){N=(int)(L/(3*3600.0)+0.5);L0=3*N*3600.0;l1=(L-L0)/p;}//采用克氏椭球参数的计算公式if(tq==1){doublec=Math.Pow(Math.Cos(b),2);doublec1=Math.Sin(b)*Math.Cos(b);doublec2=Math.Cos(b);doublel2=Math.Pow(l1,2);doublen=6399698.90

7、2-(21562.267-(108.973-0.612*c)*c)*c;doublea0=32140.404-(135.3302-(0.7092-0.0040*c)*c)*c;doublea4=(0.25+0.00252*c)*c-0.04166;doublea6=(0.166*c-0.084)*c;doublea3=(0.3333333+0.001123*c)*c-0.1666667;doublea5=0.0083-(0.1667-(0.1968+0.0040*c)*c)*c;x=6367558.4969*b-(a0-(0.5+(a

8、4+a6*l2)*l2)*l2*n)*c1;y=(1+(a3+a5*l2)*l2)*l1*n*c2;X=x;doubley1=y+500000.0;for(inti=1;y1/i>1;i=i*10){Y=N*i*10+y