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

已于 2023-08-05 21:18:41 修改
阅读量2.7k 收藏 46
点赞数 10
分类专栏： 大地测量学 文章标签： 算法 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);
            }


        }
    }
}

进击中的小龙
关注
专栏目录
高斯投影正反算表达式及代码
谢公屐的博客
01-25 2786
高斯投影正反算表达式及代码，大地测量控制点坐标转换技术规范异议点
matlab高斯投影正反算程序
10-26
利用matlab编写高斯投影正反算程序。
C#高斯投影正反算实现教程与代码
最新发布
weixin_42601608的博客
09-19 1458
本文还有配套的精品资源，点击获取 简介：高斯投影是一种在GIS领域常用的方法，用于将地球表面的经纬度转换为平面坐标。本项目通过C#语言实现高斯投影的正算与反算过程，包括公式实现及用户界面设计。项目支持多种椭球参数，能够适应不同的地理环境或历史标准。开发者需要处理精度控制和异常处理，以提供用户友好的交互界面。代码实现包括源代码、测试数据，便于理解和复用。高斯投影的应用包括地图...
高斯投影正反算EXE程序
12-27
本压缩包包括高斯投影正反算的代码及界面程序，可以实现较高精度的运算，界面设计简单易懂。
高斯投影坐标正反算程序
11-20
用程序实现高斯投影坐标正反算，用于已知x，y求B，L，或已知B，L，求x，y坐标
GPS定位（五）-高斯投影正反算C程序
小鱼儿飞丫飞
01-17 6483
基本思想 斯投影正算公式就是由大地坐标（L，B）求解高斯平面坐标（x，y），而高斯投影反算公式则是由高斯平面坐标（x，y）求解大地坐标（L，B）。 计算模型 高斯投影正算 这里输入的经纬度应该转化为度，而不是度分秒或者度分的格式 int LongLat2XY(double longitude,double latitude,double *X,double *Y) { int Pr...
高斯投影正反算（vb程序）
11-23
高斯投影正反算（vb程序） 大地坐标系与高斯投影坐标系互相转换
python实现高斯投影正反算方式
09-18
高斯投影正反算方法在地理信息系统和测绘领域中具有重要的应用价值。它是一种将地球表面的点在平面上进行表示的数学模型。高斯投影将地球的球面转换成平面坐标，使得在地图上可以方便地进行测量和导航。本文中，我们...
CGCS2000高斯投影正反算
08-24
CGCS2000坐标系高斯投影正反算，基于文件操作，批量处理，精度0.0001。软件基础框架来自氧气女神。
matlab高斯投影正反算代码及文件.rar
06-30
由于提供的压缩包文件名仅提到"matlab高斯投影正反算代码及文件"，没有具体列出代码内容，我们无法提供具体的代码示例。不过，理解高斯投影的基本原理和MATLAB实现的关键步骤是至关重要的。对于那些在文件转换部分...
高斯投影正反算，坐标系转换python源码
09-17
基于高斯投影正反算，坐标系转换，输入坐标为高斯直角坐标，地理坐标输入格式为D.mmss,精度较高，初学编程~多多指教
高斯投影正反算c语言程序
06-27
大地测量学高斯投影坐标正反算的c语言代码，简单适合初学者
高斯投影正反算C++源代码
09-02
////////////////////////////////////////////////////////////////////// // Construction/Destruction ////////////////////////////////////////////////////////////////////// const double FN = 0.;//南半球FN ＝ 1000000. const double k0 = 1.;//UTM投影中k0=0.9996 gs::gs() { } gs::~gs() { } //初始化坐标系统，并告诉程序您要使用的是3度带还是6度带 void gs::InitCordSys(enum gs::CordParam cord,enum gs::ProjKind projKind) { switch(cord){ case BeiJing1954: { m_a = 6378245.; m_b = 6356863.0188; } break; case XiAn1980: { m_a = 6378140.; m_b = 6356755.2882; } break; case WGS84: { m_a = 6378137.; m_b = 6356752.3142; } break; }
高斯投影正反算(c#编写）
12-26
大地测量程序必备，本程序使用了C#编写，操作简单，便于实现，如果有也在学习大地测量学的同学可以进行下载参考学习。
高斯投影正反算matlab程序
03-02
高斯-克吕格投影是由德国数学家、物理学家、天文学家高斯于19 世纪20 年代拟定，后经德国大地测量学家克吕格于1912 年对投影公式加以补充，故称为高斯-克吕格投影，又名"等角横切椭圆柱投影”，是地球椭球面和平面间正形投影的一种。 程序绝对可靠
C++高斯投影正反算
11-26
代码内含测试数据；主体程序都包括在内。如果要为自己所用，只须作少许修改~~
高斯投影正反算C语言程序代码,高斯投影正反算 源程序
weixin_32025789的博客
05-17 1494
#include "stdafx.h"#include "iostream.h"#include "math.h"#include "stdio.h"#define P 206264.806247096355#define PI 3.141592653589793void GaosZ_fun(){printf("高斯投影的正算\n");double l,L,B,n2,x,y,N,t,V,c,e2...
高斯投影坐标正反算算法（python）
weixin_64338372的博客
11-04 3737
高斯坐标正反算算法，在正算的时候，y的第一个数字表示是用几度带换算的，第二个数字和第三个数字表示带号。该程序是根据武汉大学出版社出版的第三版的《大地测量学基础》编写， 用于课程实验学习。
C# 大地测量高斯投影正反算公式计算程序
热门推荐
Lynn_whu的博客
09-16 1万+
C# 大地测量高斯投影正反算公式程序 C# 大地测量高斯投影正反算公式计算程序 一、编程作业要求： 编写能够进行高斯投影正反算公式计算的程序。 二、实验效果： 三、部分代码实现 1.选择界面 using System; using System.Collections.Generic; using System.ComponentModel; using S...
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、付费专栏及课程。

余额充值