c#高斯正算

本程序是一测绘学科中的高斯坐标正算，将WGS84坐标系下的坐标投影到高斯平面上。

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 a0 a3 a4 a5 a6 Nlxy;
            //double m1 = （32 * 3600 + 24 * 60 + 57.6522） / （（180 * 3600） / Math.PI）;
            //a0 = 32144.5189 - （135.3646 - （0.7034 - 0.0041 * Math.Cos（m1） * Math.Cos（m1）） * Math.Cos（m1） * Math.Cos（m1）） * Math.Cos（m1） * Math.Cos（m1）;
            //a4 = （0.25 + 0.00253 * Math.Cos（m1） * Math.Cos（m1）） * Math.Cos（m1） * Math.Cos（m1） - 0.04167;
            //a6 = （0.167 * Math.Cos（m1） * Math.Cos（m1） - 0.083） * Math.Cos（m1） * Math.Cos（m1）;
            //a3 = （0.3333333 + 0.001123 * Math.Cos（m1） * Math.Cos（m1）） * Math.Cos（m1） * Math.Cos（m1） - 0.16666667;
            //a5 = 0.00878 - （0.1702 - 0.20382 * Math.Cos（m1） * Math.Cos（m1）） * Math.Cos（m1） * Math.Cos（m1）;
            //l = （118*3600+54*60+57.6522-117*3600） / （（180 * 3600） / Math.PI）;
            //N = 6399596.652 - （21565.047 - （109.003 - （0.612 - 0.004 * Math.Cos（m1） * Math.Cos（m1）） * Math.Cos（m1） * Math.Cos（m1）） * Math.Cos（m1） * Math.Cos（m1）） * Math.Cos（m1） * Math.Cos（m1）;
            //x = 6367452.132 * m1 - （a0 - （0.5 + （a4 + a6 * l * l） * l * l） * l * l * N） * Math.Sin（m1） * Math.Cos（m1）;
            //y = N * Math.Cos（m1） * （1 + （a3 + a5 * l * l） * l * l） * l;
            //弧度转化
            int m3 = 0;
            double l;
            double a0 a3 a4 a5 a6 Nxy;
            int[] biaozhun = { 75 81 87 93 99 105 111 117 123 129 135 141 };

            double m1 = （Convert.ToDouble（32） * 3600 + Convert.ToDouble（24.96087） * 60）/（（180 * 3600） / Math.PI）;
            for （int i = 0; i < biaozhun.Length; i++）
            {
                if （Math.Abs（Convert.ToDouble（118） - biaozhun[i]） <= 3）
                {
                    m3 = biaozhun[i];
                }
            }
            a0 = 32144.5189 - （135.3646 - （0.7034 - 0.0041 * Math.Cos（m1） * Math.Cos（m1）） * Math.Cos（m1） * Math.Cos（m1）） * Math.Cos（m1） * Math.Cos（m1）;
            a4 = （0.25 + 0.00253 * Math.Cos（m1） * Math.Cos（m1）） * Math.Cos（m1） * Math.Cos（m1） - 0.04167;
            a6 = （0.167 * Math.Cos（m1） * Math.Cos（m1） - 0.083） * Math.Cos（m1） * Math.Cos（m1）;
            a3 = （0.3333333 + 0.001123 * Math.Cos（m1） * Math.Cos（m1）） * Math.Cos（m1） * Math.Cos（m1） - 0.16666667;
            a5 = 0.00878 - （0.1702 - 0.20382 * Math.Cos（m1） * Math.Cos（m1）） * Math.Cos（m1） * Math.Cos（m1）;
            l = （118 * 3600 + 54.25367666* 60 - m3 * 3600） / （（180 * 3600） / Math.PI）;
            N = 6399596.652 - （21565.047 - （109.003 - （0.612 - 0.004 * Math.Cos（m1） * Math.Cos（m1）） * Math.Cos（m1） * Math.Cos（m1）） * Math.Cos（m1） * Math.Cos（m1）） 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-12-02 21:17  高斯正算\
     目录           0  2012-12-02 21:36  高斯正算\高斯正算\
     文件         926  2012-12-02 21:17  高斯正算\高斯正算.sln
     文件       13312  2012-12-02 21:48  高斯正算\高斯正算.suo
     目录           0  2012-12-02 21:17  高斯正算\高斯正算\bin\
     目录           0  2012-12-02 21:25  高斯正算\高斯正算\bin\Debug\
     文件       11264  2012-12-02 21:36  高斯正算\高斯正算\bin\Debug\高斯正算.exe
     文件       19968  2012-12-02 21:36  高斯正算\高斯正算\bin\Debug\高斯正算.pdb
     文件       14328  2012-12-02 21:39  高斯正算\高斯正算\bin\Debug\高斯正算.vshost.exe
     文件         490  2009-06-11 05:14  高斯正算\高斯正算\bin\Debug\高斯正算.vshost.exe.manifest
     文件        3438  2012-12-02 21:36  高斯正算\高斯正算\Form1.cs
     文件        8362  2012-12-02 21:25  高斯正算\高斯正算\Form1.Designer.cs
     文件        5814  2012-12-02 21:25  高斯正算\高斯正算\Form1.resx
     目录           0  2012-12-02 21:17  高斯正算\高斯正算\obj\
     目录           0  2012-12-02 21:36  高斯正算\高斯正算\obj\Debug\
     目录           0  2012-12-02 21:17  高斯正算\高斯正算\obj\Debug\TempPE\
     文件         591  2012-12-02 21:39  高斯正算\高斯正算\obj\Debug\高斯正算.csproj.FileListAbsolute.txt
     文件         847  2012-12-02 21:25  高斯正算\高斯正算\obj\Debug\高斯正算.csproj.GenerateResource.Cache
     文件       11264  2012-12-02 21:36  高斯正算\高斯正算\obj\Debug\高斯正算.exe
     文件         180  2012-12-02 21:25  高斯正算\高斯正算\obj\Debug\高斯正算.Form1.resources
     文件       19968  2012-12-02 21:36  高斯正算\高斯正算\obj\Debug\高斯正算.pdb
     文件         180  2012-12-02 21:25  高斯正算\高斯正算\obj\Debug\高斯正算.Properties.Resources.resources
     文件         493  2012-12-02 21:17  高斯正算\高斯正算\Program.cs
     目录           0  2012-12-02 21:17  高斯正算\高斯正算\Properties\
     文件        1376  2012-12-02 21:17  高斯正算\高斯正算\Properties\AssemblyInfo.cs
     文件        2874  2012-12-02 21:17  高斯正算\高斯正算\Properties\Resources.Designer.cs
     文件        5612  2012-12-02 21:17  高斯正算\高斯正算\Properties\Resources.resx
     文件        1097  2012-12-02 21:17  高斯正算\高斯正算\Properties\Settings.Designer.cs
     文件         249  2012-12-02 21:17  高斯正算\高斯正算\Properties\Settings.settings
     文件        3735  2012-12-02 21:25  高斯正算\高斯正算\高斯正算.csproj