《大地高斯正反算-C#.docx》由会员分享，可在线阅读，更多相关《大地高斯正反算-C#.docx（6页珍藏版）》请在三一文库上搜索。
1、大地高斯正反算-C#代码如下using System;using System.Collections.Generic;using System.ComponentModel;using System.Data;using System.Drawing;using System.Linq;using System.Text;using System.Windows.Forms;namespaceWindowsFormsApplication3public partial class Form1 : Form double a, E,E1,bbb;public Form1()Initialize
2、Component();private void menuStrip1_ItemClicked( object sender, ToolStripItemClickedEventArgs e)private void toolStripTextBox1_Click( object sender, EventArgs e) private void toolStripMenuItem1_Click( object sender, EventArgs e)private void button1_Click( object sender, EventArgs e) if (a != 6378245
3、.00 & a != 6378137 & a != 6378137 & a!=6378140) textBox3.Text = t1?选?择？坐？标括？系心”；textBox4.Text = t1?选?择？坐？标括？系心”；elseint du, fen, miao; double B, L, X;double N, x, y, t, n, m, m0, m2, m4, m6, m8,l;B= Convert .ToDouble(textBox1.Text) * Math.PI /180Convert .ToDouble(textBox2.Text)/60 * Math.PI / 180 +
4、Convert .ToDouble(textBox9.Text)/3600 * Math.PI / 180;L= Convert .ToDouble(textBox10.Text) * Math.PI /180Convert .ToDouble(textBox11.Text)/60*Math.PI/180Convert .ToDouble(textBox12.Text)/3600 *Math.PI / 180; ;l = L -Convert .ToInt16(L-3.0) / 6) * 6-3.0 ;l = l *Math.PI / 180;N = a /Math.Pow(1 - E * M
5、ath.Sin(B) * Math.Sin(B), 0.5);m0 = a * (1 - E);m2 = 1.5 * E * m0;m4 = 5 / 4 * E * m2;m6 = 7 / 6 * E * m4;m8 = 9 / 8 * E * m6;X = (m0 + 0.5 * m2 + 3 / 8 * m4 + 5 / 16*m6 +35/ 128 * m8) *B - (0.5* m2 +0.5 * m4 + 15 / 32 * m6 + 7 / 16 * m8) / 2 *Math.Sin(2 *B)+ (m4 / 8 +3 /16*m6+ 7/ 32 * m8) / 4 * Mat
6、h.Sin(4 * B) - (m6/ 32 +m8 /16) / 6 *Math.Sin(6*B)+m8/ 128 / 8 * Math.Sin(8 * B);t =Math.Tan(B);n =Math.Cos(B) *Math.Cos(B) * E1;m =Math.Cos(B) *l;x = X + N / 2 * t *Math.Cos(B) * Math.Cos(B) * l * l + N / 24 * t * (5 -t * t + 9 * n + 4 * n * n) *Math.Pow(Math.Cos(B) * l, 4) + N / 720 * t * (61 - 58
7、 * t *t + t * t * t * t) *Math.Pow( Math.Cos(B) * l, 6);y = N * (1 + (1 / 6 * (1 - t * t + n) + 1 / 120 * (5 - 18 * t * t + t * t * t * t + 14 * n - 58 * n * t * t) * m * m) * m * m) * m);textBox3.Text =Convert .ToString(x);textBox4.Text =Convert .ToString(y);private void button2_Click( object sende
8、r, EventArgs e)int du, fen, miao;if (a != 6378245.00 & a != 6378137 & a != 6378137 & a != 6378140) textBox7.Text = t1?选?择？坐？标括？系心”；textBox8.Text = t1?选?择？坐？标括？系心”； else double Bf, bb, Vf, tf, xf, yf, nf, Nf, B0, l0;xf =Convert .ToDouble(textBox5.Text);yf =Convert .ToDouble(textBox6.Text);bb = xf / b
9、bb;Bf = bb + (50221746 + (293622 + (2350 + 22 *Math.Cos(bb) * Math.Cos(bb)* Math.Cos(bb) * Math.Cos(bb) * Math.Cos(bb) * Math.Cos(bb) * Math.Sin(bb) * Math.Cos(bb) * Math.Pow(10, -10);tf =Math.Tan(Bf);nf =Math.Cos(Bf) *Math.Cos(Bf) * E1;Vf =Math.Pow(1 + E1,0.5);Nf = a /Math.Pow(1 - E * Math.Sin(Bf)
10、* Math.Sin(Bf), 0.5);B0 = Bf * 180 /Math.PI - 0.5 * Vf * Vf * tf * (yf * yf / Nf / Nf - 1 / 12.00* (5 + 3 * tf * tf + nf - 9 * nf * tf * tf) *Math.Pow(yf / Nf, 4) + 1 / 360.00 * (61 + 90* tf * tf + 45 * tf * tf) *Math.Pow(yf / Nf, 6) * 180 / Math.PI;l0= 1 / Math.Cos(Bf) * (yf / Nf - 1 / 6 * (1 + 2 *
11、 tf * tf + nf) * Math.Pow(yf/ Nf, 3) + 1 / 120 * (5 + 28 * tf * tf + 24 * tf * tf + 6 * nf + 8 * nf * tf * tf) * Math.Pow(yf / Nf, 5) * 180 / Math.PI;du =Convert .ToInt16( Math.Floor(B0); /dufen =Convert .ToInt16( Math.Floor(B0 - Math.Floor(B0) * 60); /fenmiao =Convert .ToInt16( Convert .ToInt16( (
12、(B0 - Math.Floor(B0) * 60- Math.Floor(B0 - Math.Floor(B0)* 60) ) * 60);/ 秒？textBox7.Text =Convert .ToString(du) +度 e + Convert .ToString(fen)+分？+Convert .ToString(miao)+ 秒？；du =Convert .ToInt16( Math.Floor(l0); /dufen =Convert .ToInt16( Math.Floor(l0 - Math.Floor(l0) * 60);/fenmiao= Convert .ToInt16
13、( Convert .ToInt16(l0- Math.Floor(l0)Math.Floor(l0 - Math.Floor(l0) * 60) * 60);/ 秒?60 -textBox8.Text =Convert .ToString(du) +分 ? + Convert .ToString(miao) +秒? ;度 e + Convert .ToString(fen) +private void toolStripMenuItem2_Click( a = 6378245;E = 0.00693421622966;E1 = 0.006738525414683;bbb = 6367588.
14、4969;object sender,EventArgs e)private void toolStripMenuItem3_Click( a = 6378137;E = 0.0669438002290;E1 = 0.00673949677548;bbb = 6367452.133;object sender,EventArgs e)private void toolStripMenuItem4_Click( a = 6378137;E = 0.0069437999013;E1 = 0.00673949674227;bbb = 6367452.133;/ 待娜定一object sender,EventArgs e)private void toolStripMenuItem5_Click(a = 6378140;E = 0.006694384999588;E1 = 0.00673852541468;/ 待娜定一object sender,EventArgs e)private void Form1_Load( object sender, EventArgs e)