大地主题正反算

大地主题正反算程序源码，适用于测绘工程大地测量

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;

namespace 大地主题正反算
{
    public partial class Form1 : Form
    {
        public Form1（）
        {
            InitializeComponent（）;
        }
        
        private void button1_Click（object sender EventArgs e）
        {
            double dub = Convert.ToDouble（textBox1.Text）; 
            double fenb = Convert.ToDouble（textBox2.Text）;
            double miaob = Convert.ToDouble（textBox3.Text）;
            double dul = Convert.ToDouble（textBox4.Text）;
            double fenl = Convert.ToDouble（textBox5.Text）;
            double miaol = Convert.ToDouble（textBox6.Text）;
            double dua = Convert.ToDouble（textBox7.Text）;
            double fena = Convert.ToDouble（textBox8.Text）;
            double miaoa = Convert.ToDouble（textBox9.Text）;
            double s = Convert.ToDouble（textBox10.Text）;
            double b1 = （dub+fenb/60+miaob/3600） / 180 * Math.PI;
            double l1 = （dul+fenl/60+miaol/3600） / 180 * Math.PI;
            double a1 = （dua+fena/60+miaoa/3600） / 180 * Math.PI;
            double dB dL dA;
            double a = 6378245; double e2 = 0.006693421622966; double e12 = e2 / （1 - e2）; double B = b1 A = a1 L;
            double N = a * Math.Pow（（1 - e2 * Math.Sin（B） * Math.Sin（B）） 0.5）;
            double M = a * （1 - e2） * Math.Pow（（1 - e2 * Math.Sin（B） * Math.Sin（B）） -1.5）;
            double t = Math.Tan（B）;
            double yita = e12 * Math.Cos（B） * Math.Cos（B）;
            double delta_b = s * Math.Cos（A） / M;
            double delta_l = s * Math.Sin（A） / N / Math.Cos（B）;
            double delta_a = delta_l * Math.Sin（B）;
            B = b1 + 0.5 * delta_b;
            L = l1 + 0.5 * delta_l;
            A = a1 + 0.5 * delta_a;
            for （int i = 0; ; ）
            {
                dB = delta_b;
                dL = delta_l;
                dA = delta_a;
                delta_b = 1/M * s * Math.Cos（A） * （1 + s * s / 24 / N / N * （Math.Sin（A） * Math.Sin（A） * （2 + 3 * t * t + 2 * yita ） + 3 * Math.Cos（A） * Math.Cos（A） * yita *（t * t - 1 - yita - 4 * yita * t * t）））;
                delta_l = 1/N / Math.Cos（B） * s * Math.Sin（A） * （1 + s * s / 24 / N / N * （t * t * Math.Sin（A） * Math.Sin（A） - Math.Cos（A） * Math.Cos（A） * （1 + yita  - 9 * yita *  t * t）））;
                delta_a = 1/N * s * Math.Sin（A） * t * （1 + s * s / 24 / N / N * （Math.Cos（A） * Math.Cos（A） * （2 + 7 * yita + 9 * yita * t * t + 5 * Math.Pow（yita 2）） + Math.Sin（A） * Math.Sin（A） * （2 + t * t + 3 *  yita）））;
                B = b1 + 0.5 * delta_b;
                L = l1 + 0.5 * delta_l;
                A = a1 + 0.5 * delta_a;
                N = a * Math.Pow（（1 - e2 * Math.Sin（B） * Math.Sin（B）） -0.5）;
                M = a * （1 - e2） * Math.Pow（（1 - e2 * Math.Sin（B） * Math.Sin（B）） -1.5）;
      

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      30208  2012-05-05 13:12  大地主题正反算\大地主题正反算\bin\Debug\大地主题正反算.exe

     文件      38400  2012-05-05 13:12  大地主题正反算\大地主题正反算\bin\Debug\大地主题正反算.pdb

     文件      11600  2012-05-05 13:19  大地主题正反算\大地主题正反算\bin\Debug\大地主题正反算.vshost.exe

     文件        490  2010-03-17 22:39  大地主题正反算\大地主题正反算\bin\Debug\大地主题正反算.vshost.exe.manifest

     文件       8632  2012-05-05 13:12  大地主题正反算\大地主题正反算\Form1.cs

     文件      69886  2012-05-05 13:12  大地主题正反算\大地主题正反算\Form1.Designer.cs

     文件       5817  2012-05-05 13:12  大地主题正反算\大地主题正反算\Form1.resx

     文件       4440  2012-04-28 16:51  大地主题正反算\大地主题正反算\obj\x86\Debug\DesignTimeResolveAssemblyReferences.cache

     文件       6525  2012-05-05 13:12  大地主题正反算\大地主题正反算\obj\x86\Debug\DesignTimeResolveAssemblyReferencesInput.cache

     文件        344  2012-05-05 13:12  大地主题正反算\大地主题正反算\obj\x86\Debug\GenerateResource.read.1.tlog

     文件        786  2012-05-05 13:12  大地主题正反算\大地主题正反算\obj\x86\Debug\GenerateResource.write.1.tlog

     文件       7897  2012-04-28 17:03  大地主题正反算\大地主题正反算\obj\x86\Debug\ResolveAssemblyReference.cache

     文件       1321  2012-05-05 13:19  大地主题正反算\大地主题正反算\obj\x86\Debug\大地主题正反算.csproj.FileListAbsolute.txt

     文件      30208  2012-05-05 13:12  大地主题正反算\大地主题正反算\obj\x86\Debug\大地主题正反算.exe

     文件        180  2012-05-05 13:12  大地主题正反算\大地主题正反算\obj\x86\Debug\大地主题正反算.Form1.resources

     文件      38400  2012-05-05 13:12  大地主题正反算\大地主题正反算\obj\x86\Debug\大地主题正反算.pdb

     文件        180  2012-05-01 17:11  大地主题正反算\大地主题正反算\obj\x86\Debug\大地主题正反算.Properties.Resources.resources

     文件        502  2012-04-26 16:42  大地主题正反算\大地主题正反算\Program.cs

     文件       1398  2012-04-26 16:42  大地主题正反算\大地主题正反算\Properties\AssemblyInfo.cs

     文件       2892  2012-04-26 16:42  大地主题正反算\大地主题正反算\Properties\Resources.Designer.cs

     文件       5612  2012-04-26 16:42  大地主题正反算\大地主题正反算\Properties\Resources.resx

     文件       1106  2012-04-26 16:42  大地主题正反算\大地主题正反算\Properties\Settings.Designer.cs

     文件        249  2012-04-26 16:42  大地主题正反算\大地主题正反算\Properties\Settings.settings

     文件       3752  2012-05-01 17:11  大地主题正反算\大地主题正反算\大地主题正反算.csproj

     文件        905  2012-04-26 16:42  大地主题正反算\大地主题正反算.sln

    ..A..H.     20992  2012-05-05 13:22  大地主题正反算\大地主题正反算.suo

     目录          0  2012-04-26 16:42  大地主题正反算\大地主题正反算\obj\x86\Debug\TempPE

     目录          0  2012-05-05 13:12  大地主题正反算\大地主题正反算\obj\x86\Debug

     目录          0  2012-04-26 16:49  大地主题正反算\大地主题正反算\bin\Debug

     目录          0  2012-04-26 16:42  大地主题正反算\大地主题正反算\obj\x86

............此处省略8个文件信息