武大测绘大地测量学基础：大地主题解算及高斯投影（C#exe）

武汉大学测绘学院大地测量学基础编程作业之一。 包含大地主题正反算（白塞尔法和高斯引数法），高斯投影正反算。支持四种椭球（CGCS2000，WGS84，国际1975，克氏椭球）。编写时使用的算例计算的精度较高。 原本是课程的小组作业，写的时候忘了，一个人肝了一个多周。用C#写的窗体，优化了一些交互，反正函数也写出来了，就顺带写了角度转换和经度带计算的小工具。看CSDN上测绘的程序好少，就当做点小贡献吧。 写的时候为了方便检验，程序内自带算例，不过算例数据只能在程序内部改。也可以手动输入内容计算。懒得写批量计算了，快考试了，就这样放上来吧。

using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
using System.Windows.Forms;

namespace 大地主题及高斯投影
{
    public partial class Form1 : Form
    {
        public Form1（）
        {
            InitializeComponent（）;
            Size newSize = new Size（660 370）;
            this.MaximumSize = this.MinimumSize = newSize;
        }

        /*******************************变     量************************************/

        //基础值
        double a b = 0;//椭圆长短轴
        double e = 0;//第一偏心率
        double e2 = 0;
        double edot = 0;
        double edot2 = 0;//第二偏心率的平方
        double c = 0;
        double t = 0;
        double gama2 = 0;
        double M N = 0;
        double p = 180 * 3600 / Math.PI;//206265

        double L1 L2 B1 B2 S A12 A21 = 0;//秒
        double refL1 refL2 refB1 refB2 refA12 refA21 refS = 0;//算例值
        //高斯引数法
        double dL0 dB0 dA0 = 0;//秒
        double dL dB dA = 0;
        double Bm Am = 0;
        //白塞尔法
        double l = 0;//L2-L1
        double W1 W2 = 0;
        double sinu1 sinu2 cosu1 cosu2 sinA0 sin2segama1 cos2segam1 = 0;
        double A B C raf beta = 0;
        double a1 a2 b1 b2 = 0;
        double deta x = 0;
        //高斯投影
        double L0 B0 X0 Y0 = 0;
        double refL0 refB0 refX0 refY0 = 0;//算例值
        double dl = 0;//经度与中央子午线经度之差
        double Bf Nf Mf tf gamaf2 = 0;
        int NUM_1 NUM_2 centerL = 0;//6或3度带，经度带编号。中央子午线经度
        double X = 0;


        /**************************** 算 例 数 据 **********************************/
        /*大地主题解算数据
        L1B1L2B2A1A2S       35.49363300=35 49 36.3300*/
        double[] Ke_54 = new double[7]
        {
                35.49363300
                47.46526470
                36.14450504
                48.04096384
                44.12136609522834
                224.3053547669053
                44797.2817
        };
        double[] Xian_80 = new double[7]
        {
                35.49363300
                47.46526470
                36.14450504
                48.04096384
                44.12137060210981
                224.30535927378682
                44796.5461
        };
        double[] WGS_84 = new double[7]
        {
                35.49363300
                47.46526470
                36.14450504
                48.04096384
                44.121370578674686
                224.30535925035167
                44796.5250
        };
        double[] CGCS_2000 = new double[7]
         {
                35.49363300
                47.46526470
                36.14450504
                48.04096384
                44.1213705882785
                224.30535925050498
                44796.5250
         };
        /*高斯投影算例数据
   

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2020-04-29 12:08  大地主题及高斯投影\.vs\
     目录           0  2020-04-29 12:05  大地主题及高斯投影\.vs\大地主题及高斯投影\
     目录           0  2020-04-29 12:08  大地主题及高斯投影\.vs\大地主题及高斯投影\v15\
     目录           0  2020-04-29 12:05  大地主题及高斯投影\.vs\大地主题及高斯投影\v15\Server\
     目录           0  2020-04-29 12:08  大地主题及高斯投影\.vs\大地主题及高斯投影\v15\Server\sqlite3\
     目录           0  2020-04-29 12:08  大地主题及高斯投影\.vs\大地主题解算\
     目录           0  2020-05-08 20:46  大地主题及高斯投影\.vs\大地主题解算\v15\
     文件       74240  2020-05-09 12:35  大地主题及高斯投影\.vs\大地主题解算\v15\.suo
     文件    28004352  2020-05-08 20:46  大地主题及高斯投影\.vs\大地主题解算\v15\Browse.VC.db
     目录           0  2020-05-08 10:41  大地主题及高斯投影\.vs\大地主题解算\v15\ipch\
     目录           0  2020-05-08 10:43  大地主题及高斯投影\.vs\大地主题解算\v15\ipch\AutoPCH\
     目录           0  2020-05-08 10:43  大地主题及高斯投影\.vs\大地主题解算\v15\ipch\AutoPCH\b400595e02be04dc\
     目录           0  2020-04-29 12:08  大地主题及高斯投影\.vs\大地主题解算\v15\Server\
     目录           0  2020-04-29 12:08  大地主题及高斯投影\.vs\大地主题解算\v15\Server\sqlite3\
     文件           0  2020-04-29 12:08  大地主题及高斯投影\.vs\大地主题解算\v15\Server\sqlite3\db.lock
     文件      622592  2020-05-09 12:23  大地主题及高斯投影\.vs\大地主题解算\v15\Server\sqlite3\storage.ide
     文件       32768  2020-05-09 12:35  大地主题及高斯投影\.vs\大地主题解算\v15\Server\sqlite3\storage.ide-shm
     文件     4140632  2020-05-09 12:34  大地主题及高斯投影\.vs\大地主题解算\v15\Server\sqlite3\storage.ide-wal
     目录           0  2020-05-09 12:34  大地主题及高斯投影\大地主题及高斯投影\
     文件         189  2020-04-29 12:05  大地主题及高斯投影\大地主题及高斯投影\App.config
     目录           0  2020-05-08 11:52  大地主题及高斯投影\大地主题及高斯投影\bin\
     目录           0  2020-05-09 11:59  大地主题及高斯投影\大地主题及高斯投影\bin\Debug\
     文件       48640  2020-05-09 12:33  大地主题及高斯投影\大地主题及高斯投影\bin\Debug\大地主题及高斯投影.exe
     文件         189  2020-04-29 12:05  大地主题及高斯投影\大地主题及高斯投影\bin\Debug\大地主题及高斯投影.exe.config
     文件       75264  2020-05-09 12:33  大地主题及高斯投影\大地主题及高斯投影\bin\Debug\大地主题及高斯投影.pdb
     目录           0  2020-05-08 11:52  大地主题及高斯投影\大地主题及高斯投影\bin\Release\
     文件       60150  2020-05-09 12:34  大地主题及高斯投影\大地主题及高斯投影\Form1.cs
     文件       51852  2020-05-09 12:33  大地主题及高斯投影\大地主题及高斯投影\Form1.Designer.cs
     文件        5817  2020-05-09 12:33  大地主题及高斯投影\大地主题及高斯投影\Form1.resx
     目录           0  2020-04-29 12:05  大地主题及高斯投影\大地主题及高斯投影\obj\
     目录           0  2020-05-09 12:33  大地主题及高斯投影\大地主题及高斯投影\obj\Debug\
............此处省略24个文件信息