利用EGM计算高程异常和垂线偏差软件

利用EGM模型,地球重力场参数计算软件,可以单点和多点计算高程异常和垂线偏差软件，采用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;
using System.IO;
using Microsoft.VisualBasic.Devices;
using Microsoft.VisualBasic.FileIO;
using System.Collections;

namespace duodianjisuan
{
    public partial class dandianjisuan : Form
    {
        public dandianjisuan（）
        {
            InitializeComponent（）;
        }
        double[] Cnm = new double[361 361];//完全正规化位系数变量
        double[] Snm = new double[361 361];
        double[] Pnm = new double[361 361];  //缔合勒让德多项式变量
        double[] dPnm = new double[361 361];  //缔合勒让德多项式变量 （导数）
        private void button1_Click（object sender EventArgs e）
        {
            OpenFileDialog openEgm = new OpenFileDialog（）;
            openEgm.InitialDirectory = Environment.GetFolderPath（Environment.SpecialFolder.Personal）;
            openEgm.Filter = “文本文件（*.txt）|*.txt|所有文件（*.*）|*.*“;
            string textline = ““;//读取一行的值
            string FileName = ““;
            if （openEgm.ShowDialog（this） == DialogResult.OK）
            {
                int m1 = 0;//计数变量，记录文件里面有多少行数据
                FileName = openEgm.FileName;
                StreamReader reader = new StreamReader（FileName）;
                #region//读取ＥＧＭ文件里的行数　　共６５３３８行数据
                do
                {
                    textline = reader.ReadLine（）;
                    if （textline != null）
                    {
                        m1++;
                    }
                } while （textline != null）;
                reader.Close（）;
                //共６５３３８行数据
            }
                #endregion
            #region//从ＥＧＭ文本文件中读取位系数
            StreamReader reader1 = new StreamReader（FileName）;
            reader1 = new StreamReader（FileName）;
            int n = 3; int m = 0; int m2 = 0;
            MessageBox.Show（“即将读取EGM96数据，大概需要几秒时间，请稍候...“ “TIPS“
         MessageBoxButtons.OK MessageBoxIcon.Exclamation）;
            do
            {
                for （n = 2; n <= 360; n++）
                {
                    for （m = 0; m <= n; m++）
                    {
                        double nnn = 0;
                        double mmm = 0;
                        textline = reader1.ReadLine（）; m2 = m2 + 1;
                        if （textline != null）
                        {
                            string[] split = textline.Split（new char[] { ‘\t‘ }）;
                            if （split.GetLength（0） > 3）
                            {
                                nnn = double.Parse（split[2]）;
                                mmm = double.Parse（split[3]）;
                                Cnm[n m] = nnn;
                                Snm[n m] = mmm;
                            }

                        }

                    }

                    progressBar1.Value = n;
       

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2015-06-24 23:07  C#地球重力场参数计算软件\
     目录           0  2015-06-24 23:07  C#地球重力场参数计算软件\duodianjisuan\
     文件       10568  2008-05-25 20:20  C#地球重力场参数计算软件\duodianjisuan\Form1.Designer.cs
     文件        1251  2008-05-25 20:20  C#地球重力场参数计算软件\duodianjisuan\Form1.cs
     文件     3007000  2008-05-25 20:20  C#地球重力场参数计算软件\duodianjisuan\Form1.resx
     文件        2330  2008-05-21 16:19  C#地球重力场参数计算软件\duodianjisuan\Form2.Designer.cs
     文件         363  2008-05-21 16:19  C#地球重力场参数计算软件\duodianjisuan\Form2.cs
     文件        5814  2008-05-21 16:19  C#地球重力场参数计算软件\duodianjisuan\Form2.resx
     文件         494  2008-05-10 10:01  C#地球重力场参数计算软件\duodianjisuan\Program.cs
     文件      159406  2008-05-24 15:10  C#地球重力场参数计算软件\duodianjisuan\cqg2000.ico
     文件       24680  2008-05-21 15:48  C#地球重力场参数计算软件\duodianjisuan\dandianjisuan.Designer.cs
     文件       12693  2008-05-25 20:07  C#地球重力场参数计算软件\duodianjisuan\dandianjisuan.cs
     文件        8710  2008-05-21 15:48  C#地球重力场参数计算软件\duodianjisuan\dandianjisuan.resx
     文件       17563  2008-05-21 15:48  C#地球重力场参数计算软件\duodianjisuan\duodianjisuan.Designer.cs
     文件       16769  2008-05-21 15:48  C#地球重力场参数计算软件\duodianjisuan\duodianjisuan.cs
     文件        7680  2008-05-24 15:10  C#地球重力场参数计算软件\duodianjisuan\duodianjisuan.csproj
     文件         486  2008-05-21 16:06  C#地球重力场参数计算软件\duodianjisuan\duodianjisuan.csproj.user
     文件      246783  2008-05-21 15:48  C#地球重力场参数计算软件\duodianjisuan\duodianjisuan.resx
     文件        1668  2008-05-21 16:04  C#地球重力场参数计算软件\duodianjisuan\duodianjisuan_TemporaryKey.pfx
     文件        8835  2008-05-21 14:31  C#地球重力场参数计算软件\duodianjisuan\jianjie.Designer.cs
     文件        1373  2008-05-21 14:31  C#地球重力场参数计算软件\duodianjisuan\jianjie.cs
     文件        5814  2008-05-21 14:31  C#地球重力场参数计算软件\duodianjisuan\jianjie.resx
     文件       25081  2008-05-25 20:20  C#地球重力场参数计算软件\duodianjisuan\quyujisuan.Designer.cs
     文件       16903  2008-05-25 20:20  C#地球重力场参数计算软件\duodianjisuan\quyujisuan.cs
     文件        6617  2008-05-25 20:20  C#地球重力场参数计算软件\duodianjisuan\quyujisuan.resx
     文件         929  2008-05-10 10:02  C#地球重力场参数计算软件\duodianjisuan.sln
     文件       26112  2008-06-01 13:28  C#地球重力场参数计算软件\duodianjisuan.suo