高斯投影正反算及换带程序

已于 2023-08-05 21:18:41 修改
阅读量2.2k 收藏 39
点赞数 8
分类专栏： 大地测量学 文章标签： 算法 c#
于 2023-06-18 17:36:12 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/yitian_00/article/details/131273821
版权

源项目地址：高斯投影: Gaussian_projection_coordinates (gitee.com)

git clone：

git clone git@gitee.com:yitian00/gaussian-projection.git

擅长或喜欢c++qt编程的可以前往：

高斯投影正反算及换带_进击中的小龙的博客-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 Gaussian_projection_coordinates
{
    public partial class Form2 : Form
    {
        public Form2()
        {
            InitializeComponent();
        }

        private void button1_Click(object sender, EventArgs e)
        {
            if (textBox1.Text != string.Empty && textBox2.Text != string.Empty&& textBox3.Text != string.Empty &&
                textBox4.Text != string.Empty && textBox5.Text != string.Empty && textBox6.Text != string.Empty &&
                textBox7.Text != string.Empty )
            {
                double a = 0, b = 0, e_2 = 0, e_dot_2 = 0, rou = 206264.806247096355;
                int pp = 0;
                if (radioButton1.Checked)
                {
                    a = 6378245;//椭圆长半轴
                    b = 6356863.0187730473;//椭圆短半轴
                    e_2 = 0.006693421622966;//椭圆第一偏心率e^2
                    e_dot_2 = 0.006738525414683;//椭圆第二偏心率e'^2 
                    pp = 1;
                }
                else if (radioButton2.Checked)//
                {
                    a = 6378140;//椭圆长半轴
                    b = 6356755.2881575287;//椭圆短半轴
                    e_2 = 0.006694384999588;//椭圆第一偏心率e^2
                    e_dot_2 = 0.006739501819473;//椭圆第二偏心率e'^2 
                    pp = 1;
                }
                else if (radioButton3.Checked)//
                {
                    a = 6378137;//椭圆长半轴
                    b = 6356752.3142;//椭圆短半轴
                    e_2 = 0.00669437999013;//椭圆第一偏心率e^2
                    e_dot_2 = 0.00673949674227;//椭圆第二偏心率e'^2 
                    pp = 1;
                }
                else if (radioButton4.Checked)//
                {
                    a = 6378137;//椭圆长半轴
                    b = 6356752.3141;//椭圆短半轴
                    e_2 = 0.00669438002290;//椭圆第一偏心率e^2
                    e_dot_2 = 0.00673949677548;//椭圆第二偏心率e'^2 
                    pp = 1;
                }
                else if (pp == 0)
                {
                    MessageBox.Show("请选择椭球体", "提示", MessageBoxButtons.OK, MessageBoxIcon.Error);
                }
                double c = a * a / b;
                //转化为弧度制
                double B = (Convert.ToDouble(textBox1.Text) + Convert.ToDouble(textBox2.Text) / 60 + Convert.ToDouble(textBox3.Text) / 3600) / 180 * Math.PI;
                double L = (Convert.ToDouble(textBox4.Text) + Convert.ToDouble(textBox5.Text) / 60 + Convert.ToDouble(textBox6.Text) / 3600) / 180 * Math.PI;
                double L0 = (Convert.ToDouble(textBox7.Text)) / 180 * Math.PI;

                double W = Math.Sqrt(1 - e_2 * Math.Sin(B) * Math.Sin(B));
                double V = Math.Sqrt(1 + e_dot_2 * Math.Cos(B) * Math.Cos(B));
                double M = a * (1 - e_2) / (W * W * W);//子午圈曲率半径
                double N = a / W;//卯酉圈曲率半径

                double t= Math.Tan(B), n = e_dot_2 * Math.Cos(B) * Math.Cos(B);
                double l = (L - L0) * rou;//单位弧度秒
                //计算子午圈弧长
           
                double A1 = 1 + 3 / 4 * e_2 + 45 / 64 * e_2 * e_2 + 175 / 256 * Math.Pow(e_2, 3) + 11025 / 16384 * Math.Pow(e_2, 4) + 43659 / 65536 * Math.Pow(e_2, 5);
                double A2= 3 / 4 * e_2 + 15 / 16* e_2 * e_2 + 525/512 * Math.Pow(e_2, 3) + 2205/2048 * Math.Pow(e_2, 4) +72765/65536 * Math.Pow(e_2, 5);
                double A3= 15 / 64 * e_2 * e_2 + 105/256 * Math.Pow(e_2, 3) + 2205 /4096 * Math.Pow(e_2, 4) + 10395/16384 * Math.Pow(e_2, 5);
                double A4= 35/512 * Math.Pow(e_2, 3) + 315/2048 * Math.Pow(e_2, 4) + 31185/131072 * Math.Pow(e_2, 5);
                double A5 = 315 / 16384 * Math.Pow(e_2, 4) + 3465 / 65536 * Math.Pow(e_2, 5);
                double A6 = 693 / 131072 * Math.Pow(e_2, 5);
                double a0 = A1 * a * (1 - e_2), a2 = -A2 / 2 * a * (1 - e_2), a4 = A3 / 4 * a * (1 - e_2), a6 = -A4 * a / 6 * (1 - e_2), a8 = A5 * a / 8 * (1 - e_2), a10 = -A6 * a / 10 * (1 - e_2);
                double X = a0*B+a2*Math.Sin(2*B)+ a4 * Math.Sin(4* B)+ a6* Math.Sin(6 * B)+ a8  * Math.Sin(8 * B)+a10*Math.Sin(10*B);

                double x = X + N  * Math.Sin(B) * Math.Cos(B) * l * l / 2 / rou / rou +
                    N / 24 * Math.Sin(B) * Math.Pow(Math.Cos(B), 3) * (5 - t * t + 9 * n * n+4*Math.Pow(n,4)) * Math.Pow(l, 4) / Math.Pow(rou, 4)+
                    N*Math.Sin(B)*Math.Pow(Math.Cos(B),5)*(61-58*t*t+Math.Pow(t,4))*Math.Pow(l,6)/720/Math.Pow(rou,6) ;
                double y = N * Math.Cos(B) * l / rou + N * Math.Pow(Math.Cos(B), 3) * (1 - t * t + n * n) * Math.Pow(l, 3) / 6 / Math.Pow(rou, 3)
                    + N * Math.Pow(Math.Cos(B), 5) * (5 - 18 * t * t + Math.Pow(t, 4)) * Math.Pow(l, 5) / 120 / Math.Pow(rou, 5);

                double gamma = Math.Sin(B) * l/rou + Math.Sin(B) * Math.Cos(B) * Math.Cos(B) * Math.Pow(l/rou, 3) / 3 * (1 + 3 * n * n + 2 * Math.Pow(n, 4)) +
                    Math.Sin(B) * Math.Pow(Math.Cos(B), 4) * Math.Pow(l/rou, 5) / 15 * (2 - t * t);//弧度制

                gamma = gamma * 180 / Math.PI;//角度制
                //打印
                textBox8.Text = Convert.ToString(x);
                textBox9.Text = Convert.ToString(y);

                //子午线收敛角
                textBox10.Text = Convert.ToString((int)(gamma));
                textBox11.Text = Convert.ToString((int)((gamma -(int)(gamma)) * 60));
                textBox12.Text = Convert.ToString(((gamma - (int)(gamma)) * 60 - (int)((gamma - (int)(gamma)) * 60)) * 60);



            }
            else
            {
                MessageBox.Show("请输入完整参数", "提示", MessageBoxButtons.OK, MessageBoxIcon.Error);
            }

        }

        private void button2_Click(object sender, EventArgs e)
        {
            textBox1.Text = string.Empty;
            textBox2.Text = string.Empty;
            textBox3.Text = string.Empty;
            textBox4.Text = string.Empty;
            textBox5.Text = string.Empty;
            textBox6.Text = string.Empty;
            textBox7.Text = string.Empty;
            textBox8.Text = string.Empty;
            textBox9.Text = string.Empty;
            textBox10.Text = string.Empty;
            textBox11.Text = string.Empty;
            textBox12.Text = string.Empty;
           
        }
    }
}

 反算及换带窗口

 反算算法：

 子午线弧长反算公式：

 领带换算算法

 

 代码：

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 Gaussian_projection_coordinates
{
    public partial class Form3 : Form
    {
        public Form3()
        {
            InitializeComponent();
        }
        //反算
        private void button1_Click(object sender, EventArgs e)
        {
            if (textBox1.Text != string.Empty && textBox2.Text != string.Empty && textBox3.Text != string.Empty)
            {
                double a = 0, b = 0, e_2 = 0, e_dot_2 = 0, rou = 206264.806247096355;
                int pp = 0;
                if (radioButton1.Checked)
                {
                    a = 6378245;//椭圆长半轴
                    b = 6356863.0187730473;//椭圆短半轴
                    e_2 = 0.006693421622966;//椭圆第一偏心率e^2
                    e_dot_2 = 0.006738525414683;//椭圆第二偏心率e'^2 
                    pp = 1;
                }
                else if (radioButton2.Checked)//
                {
                    a = 6378140;//椭圆长半轴
                    b = 6356755.2881575287;//椭圆短半轴
                    e_2 = 0.006694384999588;//椭圆第一偏心率e^2
                    e_dot_2 = 0.006739501819473;//椭圆第二偏心率e'^2 
                    pp = 1;
                }
                else if (radioButton3.Checked)//
                {
                    a = 6378137;//椭圆长半轴
                    b = 6356752.3142;//椭圆短半轴
                    e_2 = 0.00669437999013;//椭圆第一偏心率e^2
                    e_dot_2 = 0.00673949674227;//椭圆第二偏心率e'^2 
                    pp = 1;
                }
                else if (radioButton4.Checked)//
                {
                    a = 6378137;//椭圆长半轴
                    b = 6356752.3141;//椭圆短半轴
                    e_2 = 0.00669438002290;//椭圆第一偏心率e^2
                    e_dot_2 = 0.00673949677548;//椭圆第二偏心率e'^2 
                    pp = 1;
                }
                else if (pp == 0)
                {
                    MessageBox.Show("请选择椭球体", "提示", MessageBoxButtons.OK, MessageBoxIcon.Error);
                }
                double c = a * a / b;
                //转化为弧度制
                double x = Convert.ToDouble(textBox1.Text);
                double y = Convert.ToDouble(textBox2.Text);
                double L0 = Convert.ToDouble(textBox3.Text);
                //子午线弧长求大地纬度(迭代法)
                double Bf = x /111134.8611,cha;//单位度
                double F = -16036.4803 * Math.Sin(2 * Bf/180*Math.PI) + 16.8281 * Math.Sin(4 * Bf / 180 * Math.PI) - 0.022 * Math.Sin(6 * Bf / 180 * Math.PI);
                do
                {
                    cha = Bf;
                    F = -16036.4803 * Math.Sin(2 * Bf / 180 * Math.PI) + 16.8281 * Math.Sin(4 * Bf / 180 * Math.PI) - 0.022 * Math.Sin(6 * Bf / 180 * Math.PI);
                    Bf = (x - F) /111134.8611;
                } while (Math.Abs(cha - Bf) > 1e-8);
                Bf = Bf / 180 * Math.PI;//单位弧度

                double Wf = Math.Sqrt(1 - e_2 * Math.Sin(Bf) * Math.Sin(Bf));
                double Vf = Math.Sqrt(1 + e_dot_2 * Math.Cos(Bf) * Math.Cos(Bf));
                double Mf = a * (1 - e_2) / (Wf * Wf * Wf);
                double Nf = a / Wf;

                double tf = Math.Tan(Bf), nf = e_dot_2 * Math.Cos(Bf) * Math.Cos(Bf);

                double gamma = rou * y * tf / Nf - rou * y * y * y * tf / 3 / Math.Pow(Nf, 3) * (1 + tf * tf - nf * nf) +
                    rou * Math.Pow(y, 5) * tf / 15 / Math.Pow(Nf, 5) * (2 + 5 * tf * tf + 3 * Math.Pow(tf, 4));//单位弧度秒
                gamma = gamma * 180 / Math.PI/rou;//单位度
                double B = Bf - tf * y * y / 2 / Mf / Math.Pow(Nf, 3) + tf * (5 + 3 * tf * tf + nf * nf - 9 * nf * nf * tf * tf) * Math.Pow(y, 4) / 24 / Mf / Math.Pow(Nf, 3) -
                    tf * (61 + 90 * tf * tf + 45 * Math.Pow(tf, 4)) * Math.Pow(y, 6) / 720 / Mf / Math.Pow(Nf, 5);
                double l = y / Nf / Math.Cos(Bf) - (1 + 2 * tf * tf + nf * nf) * Math.Pow(y, 3) / 6 / Math.Cos(Bf) / Math.Pow(Nf, 3) +
                    (5 + 28 * tf * tf + 24 * Math.Pow(tf, 4) + 6 * nf * nf + 8 * tf * tf * nf * nf) * Math.Pow(y, 5) / 120 / Math.Cos(Bf) / Math.Pow(Nf, 5);
                double L = l*180/Math.PI + L0;//单位度
                 B = B * 180 / Math.PI ;
                //打印纬度
                textBox9.Text = Convert.ToString(Math.Floor(B));
                textBox8.Text = Convert.ToString(Math.Floor((B - Math.Floor(B)) * 60));
                textBox7.Text = Convert.ToString(((B - Math.Floor(B)) * 60 - Math.Floor((B - Math.Floor(B)) * 60)) * 60);
                //打印经度
                textBox4.Text = Convert.ToString(Math.Floor(L));
                textBox5.Text = Convert.ToString(Math.Floor((L- Math.Floor(L)) * 60));
                textBox6.Text = Convert.ToString(((L - Math.Floor(L)) * 60 - Math.Floor((L - Math.Floor(L)) * 60)) * 60);
                //打印子午线收敛角
                textBox10.Text = Convert.ToString((int)(gamma));
                textBox11.Text = Convert.ToString((int)((gamma - (int)(gamma)) * 60));
                textBox12.Text = Convert.ToString(((gamma - (int)(gamma)) * 60 - (int)((gamma - (int)(gamma)) * 60)) * 60);

            }
            else
            {
                MessageBox.Show("请输入完整参数", "提示", MessageBoxButtons.OK, MessageBoxIcon.Error);
            }
            }


        private void button2_Click(object sender, EventArgs e)
        {
            textBox1.Text = string.Empty;
            textBox2.Text = string.Empty;
            textBox3.Text = string.Empty;
            textBox4.Text = string.Empty;
            textBox5.Text = string.Empty;
            textBox6.Text = string.Empty;
            textBox7.Text = string.Empty;
            textBox8.Text = string.Empty;
            textBox9.Text = string.Empty;


        }


        //换带
        private void button3_Click(object sender, EventArgs e)
        {
            if (textBox1.Text != string.Empty && textBox2.Text != string.Empty && textBox3.Text != string.Empty && textBox13.Text != string.Empty)
            {
                double a = 0, b = 0, e_2 = 0, e_dot_2 = 0, rou = 206264.806247096355;
                int pp = 0;
                if (radioButton1.Checked)
                {
                    a = 6378245;//椭圆长半轴
                    b = 6356863.0187730473;//椭圆短半轴
                    e_2 = 0.006693421622966;//椭圆第一偏心率e^2
                    e_dot_2 = 0.006738525414683;//椭圆第二偏心率e'^2 
                    pp = 1;
                }
                else if (radioButton2.Checked)//
                {
                    a = 6378140;//椭圆长半轴
                    b = 6356755.2881575287;//椭圆短半轴
                    e_2 = 0.006694384999588;//椭圆第一偏心率e^2
                    e_dot_2 = 0.006739501819473;//椭圆第二偏心率e'^2 
                    pp = 1;
                }
                else if (radioButton3.Checked)//
                {
                    a = 6378137;//椭圆长半轴
                    b = 6356752.3142;//椭圆短半轴
                    e_2 = 0.00669437999013;//椭圆第一偏心率e^2
                    e_dot_2 = 0.00673949674227;//椭圆第二偏心率e'^2 
                    pp = 1;
                }
                else if (radioButton4.Checked)//
                {
                    a = 6378137;//椭圆长半轴
                    b = 6356752.3141;//椭圆短半轴
                    e_2 = 0.00669438002290;//椭圆第一偏心率e^2
                    e_dot_2 = 0.00673949677548;//椭圆第二偏心率e'^2 
                    pp = 1;
                }
                else if (pp == 0)
                {
                    MessageBox.Show("请选择椭球体", "提示", MessageBoxButtons.OK, MessageBoxIcon.Error);
                }
                double c = a * a / b;
                //转化为弧度制
                double x = Convert.ToDouble(textBox1.Text);
                double y = Convert.ToDouble(textBox2.Text);
                double L0 = Convert.ToDouble(textBox3.Text);
                //子午线弧长求大地纬度(迭代法)
                double Bf = x / 111134.8611, cha;//单位度
                double F = -16036.4803 * Math.Sin(2 * Bf / 180 * Math.PI) + 16.8281 * Math.Sin(4 * Bf / 180 * Math.PI) - 0.022 * Math.Sin(6 * Bf / 180 * Math.PI);
                do
                {
                    cha = Bf;
                    F = -16036.4803 * Math.Sin(2 * Bf / 180 * Math.PI) + 16.8281 * Math.Sin(4 * Bf / 180 * Math.PI) - 0.022 * Math.Sin(6 * Bf / 180 * Math.PI);
                    Bf = (x - F) / 111134.8611;
                } while (Math.Abs(cha - Bf) > 1e-8);
                Bf = Bf / 180 * Math.PI;//单位弧度

                double Wf = Math.Sqrt(1 - e_2 * Math.Sin(Bf) * Math.Sin(Bf));
                double Vf = Math.Sqrt(1 + e_dot_2 * Math.Cos(Bf) * Math.Cos(Bf));
                double Mf = a * (1 - e_2) / (Wf * Wf * Wf);
                double Nf = a / Wf;

                double tf = Math.Tan(Bf), nf = e_dot_2 * Math.Cos(Bf) * Math.Cos(Bf);

                double gamma = rou * y * tf / Nf - rou * y * y * y * tf / 3 / Math.Pow(Nf, 3) * (1 + tf * tf - nf * nf) +
                    rou * Math.Pow(y, 5) * tf / 15 / Math.Pow(Nf, 5) * (2 + 5 * tf * tf + 3 * Math.Pow(tf, 4));//单位弧度秒
                gamma = gamma * 180 / Math.PI / rou;//单位度
                
               double B = Bf - tf * y * y / 2 / Mf / Math.Pow(Nf, 3) + tf * (5 + 3 * tf * tf + nf * nf - 9 * nf * nf * tf * tf) * Math.Pow(y, 4) / 24 / Mf / Math.Pow(Nf, 3) -
                    tf * (61 + 90 * tf * tf + 45 * Math.Pow(tf, 4)) * Math.Pow(y, 6) / 720 / Mf / Math.Pow(Nf, 5);
                double l = y / Nf / Math.Cos(Bf) - (1 + 2 * tf * tf + nf * nf) * Math.Pow(y, 3) / 6 / Math.Cos(Bf) / Math.Pow(Nf, 3) +
                    (5 + 28 * tf * tf + 24 * Math.Pow(tf, 4) + 6 * nf * nf + 8 * tf * tf * nf * nf) * Math.Pow(y, 5) / 120 / Math.Cos(Bf) / Math.Pow(Nf, 5);
                double L = l * 180 / Math.PI + L0;//单位度
                B = B * 180 / Math.PI;
                //打印纬度
                textBox9.Text = Convert.ToString(Math.Floor(B));
                textBox8.Text = Convert.ToString(Math.Floor((B - Math.Floor(B)) * 60));
                textBox7.Text = Convert.ToString(((B - Math.Floor(B)) * 60 - Math.Floor((B - Math.Floor(B)) * 60)) * 60);
                //打印经度
                textBox4.Text = Convert.ToString(Math.Floor(L));
                textBox5.Text = Convert.ToString(Math.Floor((L - Math.Floor(L)) * 60));
                textBox6.Text = Convert.ToString(((L - Math.Floor(L)) * 60 - Math.Floor((L - Math.Floor(L)) * 60)) * 60);
                //打印子午线收敛角
                textBox10.Text = Convert.ToString((int)(gamma));
                textBox11.Text = Convert.ToString((int)((gamma - (int)(gamma)) * 60));
                textBox12.Text = Convert.ToString(((gamma - (int)(gamma)) * 60 - (int)((gamma - (int)(gamma)) * 60)) * 60);

               
               
             //求新坐标的x,y
                //转化为弧度制
                B = B/ 180 * Math.PI;
                L = L / 180 * Math.PI;
                L0 = (Convert.ToDouble(textBox13.Text)) / 180 * Math.PI;

                double W = Math.Sqrt(1 - e_2 * Math.Sin(B) * Math.Sin(B));
                double V = Math.Sqrt(1 + e_dot_2 * Math.Cos(B) * Math.Cos(B));
                double M = a * (1 - e_2) / (W * W * W);//子午圈曲率半径
                double N = a / W;//卯酉圈曲率半径

                double t = Math.Tan(B), n = e_dot_2 * Math.Cos(B) * Math.Cos(B);
               l = (L - L0) * rou;//单位弧度秒
                                  //计算子午圈弧长
                double A1 = 1 + 3 / 4 * e_2 + 45 / 64 * e_2 * e_2 + 175 / 256 * Math.Pow(e_2, 3) + 11025 / 16384 * Math.Pow(e_2, 4) + 43659 / 65536 * Math.Pow(e_2, 5);
                double A2 = 3 / 4 * e_2 + 15 / 16 * e_2 * e_2 + 525 / 512 * Math.Pow(e_2, 3) + 2205 / 2048 * Math.Pow(e_2, 4) + 72765 / 65536 * Math.Pow(e_2, 5);
                double A3 = 15 / 64 * e_2 * e_2 + 105 / 256 * Math.Pow(e_2, 3) + 2205 / 4096 * Math.Pow(e_2, 4) + 10395 / 16384 * Math.Pow(e_2, 5);
                double A4 = 35 / 512 * Math.Pow(e_2, 3) + 315 / 2048 * Math.Pow(e_2, 4) + 31185 / 131072 * Math.Pow(e_2, 5);
                double A5 = 315 / 16384 * Math.Pow(e_2, 4) + 3465 / 65536 * Math.Pow(e_2, 5);
                double A6 = 693 / 131072 * Math.Pow(e_2, 5);
                double a0 = A1 * a * (1 - e_2), a2 = -A2 / 2 * a * (1 - e_2), a4 = A3 / 4 * a * (1 - e_2), a6 = -A4 * a / 6 * (1 - e_2), a8 = A5 * a / 8 * (1 - e_2), a10 = -A6 * a / 10 * (1 - e_2);
                double X = a0 * B + a2 * Math.Sin(2 * B) + a4 * Math.Sin(4 * B) + a6 * Math.Sin(6 * B) + a8 * Math.Sin(8 * B) + a10 * Math.Sin(10 * B);

                x = X + N * Math.Sin(B) * Math.Cos(B) * l * l / 2 / rou / rou +
                    N / 24 * Math.Sin(B) * Math.Pow(Math.Cos(B), 3) * (5 - t * t + 9 * n * n + 4 * Math.Pow(n, 4)) * Math.Pow(l, 4) / Math.Pow(rou, 4) +
                    N * Math.Sin(B) * Math.Pow(Math.Cos(B), 5) * (61 - 58 * t * t + Math.Pow(t, 4)) * Math.Pow(l, 6) / 720 / Math.Pow(rou, 6);
                 y = N * Math.Cos(B) * l / rou + N * Math.Pow(Math.Cos(B), 3) * (1 - t * t + n * n) * Math.Pow(l, 3) / 6 / Math.Pow(rou, 3)
                    + N * Math.Pow(Math.Cos(B), 5) * (5 - 18 * t * t + Math.Pow(t, 4)) * Math.Pow(l, 5) / 120 / Math.Pow(rou, 5);

  
                gamma = Math.Sin(B) * l / rou + Math.Sin(B) * Math.Cos(B) * Math.Cos(B) * Math.Pow(l / rou, 3) / 3 * (1 + 3 * n * n + 2 * Math.Pow(n, 4)) +
                    Math.Sin(B) * Math.Pow(Math.Cos(B), 4) * Math.Pow(l / rou, 5) / 15 * (2 - t * t);//弧度制
                gamma = gamma * 180 / Math.PI ;//单位度

                //打印新坐标
                textBox15.Text = Convert.ToString(x);
                textBox14.Text = Convert.ToString(y);
                //打印新子午线收敛角
                textBox18.Text = Convert.ToString( (int) gamma);
                textBox17.Text = Convert.ToString((int)((gamma - (int)(gamma)) * 60));
                textBox16.Text = Convert.ToString(((gamma - (int)(gamma)) * 60 - (int)((gamma - (int)(gamma)) * 60)) * 60);

            }
            else
            {
                MessageBox.Show("请输入完整参数", "提示", MessageBoxButtons.OK, MessageBoxIcon.Error);
            }


        }
    }
}

进击中的小龙
关注
  • 8
    点赞
  • 踩
  • 39
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
专栏目录
GPS定位（五）-高斯投影正反算C程序
小鱼儿飞丫飞
01-17 6183
基本思想 斯投影正算公式就是由大地坐标（L，B）求解高斯平面坐标（x，y），而高斯投影反算公式则是由高斯平面坐标（x，y）求解大地坐标（L，B）。 计算模型 高斯投影正算 这里输入的经纬度应该转化为度，而不是度分秒或者度分的格式 int LongLat2XY(double longitude,double latitude,double *X,double *Y) { int Pr...
高斯正反算—投影坐标转大地坐标、大地坐标转投影坐标(附有完整代码及测试结果)
热门推荐
热爱数据库的小胖
04-23 1万+
本文意在介绍高斯正反算的基本原理和代码实现！针对高斯正反算网上给出的方法很多，但是我试了之后发现多少都有些问题：或公式原理问题、或精度问题！ 通过参考别人的博客进行了总结与测试(费了很大一番功夫)！原理！代码！测试结果！本文都一一给出，此外本文还对常见坐标系：北京54、西安80、WGS84、CGCS2000等坐标系的高斯转换都做出了实现并使用 QT 进行封装可视化！ 一、常用椭球参数 北京54坐标系 西安80坐标...
高斯坐标正反算c语言,高斯投影坐标正反算编程报告
weixin_35723352的博客
05-16 999
《高斯投影坐标正反算编程报告》由会员分享，可在线阅读，更多相关《高斯投影坐标正反算编程报告(11页珍藏版)》请在人人文库网上搜索。1、高斯投影坐标正反算编程报告1. 编程思想进行高斯投影坐标正反算的编程需要牵涉到大量的公式，为了使程序条理更清楚，各块的数据复用性更强，这里采取了结构化的编程思想。程序由四大块组成。GeodesyHomework.cpp文件用于存放main()函数，是整个程序的入口。...
高斯投影正反算EXE程序
12-27
本压缩包包括高斯投影正反算的代码及界面程序，可以实现较高精度的运算，界面设计简单易懂。
高斯投影坐标正反算算法（python）
weixin_64338372的博客
11-04 3426
高斯坐标正反算算法，在正算的时候，y的第一个数字表示是用几度带换算的，第二个数字和第三个数字表示带号。该程序是根据武汉大学出版社出版的第三版的《大地测量学基础》编写， 用于课程实验学习。
高斯投影法正反算代码MATLAB版本
weixin_44714553的博客
04-17 6514
高斯投影法正反算代码MATLAB版本说明高斯投影正算高斯投影反算 说明 之前一个帖子给出了正反算的C语言代码 链接如下： link 我把它移植到了matlab中 支持向量输入 高斯投影正算 高斯投影正算matlab代码 function [X, Y] = GaussProjCal(Lon, Lat) %% 高斯投影正算： 经纬度 算 直角坐标 % Longitude Latitude 均为角度制 Lon = deg2rad(Lon); Lat = deg2rad(Lat); ZoneWide = de
适用于计算机程序的高斯投影正反算公式及程序编译遇到的问题
qq_38431115的博客
10-01 3660
在网络上查询相关资料时，遇到大量错误公式，且在程序编写中应所使用公式返回单位的差异始终无法得出正结果。期间遇到的“坑”均已注释的形式写明。按本公式直接编译即可计算正确。
高斯投影正反算表达式及代码
最新发布
谢公屐的博客
01-25 1890
高斯投影正反算表达式及代码，大地测量控制点坐标转换技术规范异议点
python实现高斯投影正反算方式
09-18
今天小编就为大家分享一篇python实现高斯投影正反算方式，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
CGCS2000高斯投影正反算
08-24
CGCS2000坐标系高斯投影正反算，基于文件操作，批量处理，精度0.0001。软件基础框架来自氧气女神。
高斯投影的正反算C#代码
05-22
C# 高 斯 投 影 正 反 算，适 用 于 坐 标 转 换。
高斯投影正反算（vb程序）
11-23
高斯投影正反算（vb程序） 大地坐标系与高斯投影坐标系互相转换
高斯投影坐标正反算
07-03
进行高斯投影坐标正反算的编程需要牵涉到大量的公式，为了使程序条理更清楚，各块的数据复用性更强，这里采取了结构化的编程思想。 程序由四大块组成。 GeodesyHomework.cpp文件用于存放main()函数，是整个程序的入口。通过结构化的编程尽力使main()函数变得简单。 MyFunction.h和MyFunction.cpp用于存放计算过程中进行角度弧度换算时所要用到的一些自定的转换函数。 Zhengsuan.h和Zhengsuan.cpp用于存放Zhengsuan类，在Zhengsuan类中声明了高斯投影坐标正算所要用到的所有变量，在类的构造函数中进行成员变量的初始化及正算计算。通过get函数获得相应的正算结果。 Fansuan.h和Fansuan.cpp用于存放Fansuan类，类似于Zhengsuan类，Fansuan类中声明了高斯投影坐标反算所要用到的所有变量，在类的构造函数中进行成员变量的初始化及反算计算。通过get函数获得相应的反算结果。
高斯正反算设计实现!!!
04-23
针对高斯正反算网上给出的方法很多，但是我试了之后发现多少都有些问题：或公式原理问题、或精度问题！通过查找资料对其进行了总结与测试(费了很大一番功夫)！原理！代码！测试结果！本文都一一给出，此外本文还对常见坐标系：北京54、西安80、WGS84、CGCS2000、自定义坐标系等坐标系的高斯转换都做出了实现并使用 QT 进行封装可视化！
高斯正反算（高斯投影计算）
12-02
可用于高斯投影的计算，即已知大地坐标BLH计算投影后的平面坐标，或者由平面坐标反算大地坐标。包含全部源代码，供大地测量初学者使用。
matlab高斯投影正反算代码及文件.rar
06-30
由于提供的压缩包文件名仅提到"matlab高斯投影正反算代码及文件"，没有具体列出代码内容，我们无法提供具体的代码示例。不过，理解高斯投影的基本原理和MATLAB实现的关键步骤是至关重要的。对于那些在文件转换部分...
matlab高斯投影正反算程序
10-26
下面我们将详细探讨MATLAB中的高斯投影正反算程序及其相关知识。 一、高斯投影的基本概念 高斯投影，也称为克吕格投影，由德国数学家和天文学家弗里德里希·威廉·贝塞尔于19世纪提出。这种投影方法以中央经线为...
高斯投影正反算，坐标系转换python源码
09-17
基于高斯投影正反算，坐标系转换，输入坐标为高斯直角坐标，地理坐标输入格式为D.mmss,精度较高，初学编程~多多指教
【测绘程序设计】——高斯投影
C_xxy's blog
12-11 1721
测绘程序设计——高斯投影（C#版和Python版）
python 高斯投影正反算程序
11-23
高斯投影是一种常用的地球表面坐标与平面坐标之间的转换方法。Python语言可以通过一些库来实现高斯投影的正反算，以下是一个简单的实现示例： 正算是将地球表面的经纬度坐标转换为高斯投影坐标。首先，需要确定所使用的高斯投影坐标系的参数，包括中央子午线经度、投影原点纬度、假轴、假东原点等参数。在Python中可以使用库例如`pyproj`来进行坐标转换。首先，需要在代码中导入`pyproj`库。然后，创建一个`pyproj.Proj`对象来定义高斯投影坐标系的参数。最后，通过调用`proj()`方法传入经纬度坐标，即可得到对应的高斯投影坐标。 反算是将高斯投影坐标转换为地球表面的经纬度坐标。同样，需要使用`pyproj`库进行坐标转换。首先创建一个`pyproj.Proj`对象来定义高斯投影坐标系的参数。然后，通过调用`inverse()`方法传入高斯投影坐标，即可得到对应的经纬度坐标。 总结起来，使用Python实现高斯投影的正反算，可以通过导入`pyproj`库创建`pyproj.Proj`对象，并调用相应的方法实现坐标的正反转换。该库提供了很多常用的坐标转换方法，具有较高的精度和稳定性，适合地理信息系统开发等应用领域。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值