大地測量學白塞爾大地主題解算

阿新 • • 發佈：2018-11-26

namespace 大地測量學1
{
    public partial class Form1 : Form
    {
        public Form1()
        {
            InitializeComponent();
        }

        private void button1_Click(object sender, EventArgs e)
        {
            double L1 = 130 + (10 * 60 + 12.2676) / 3600.0, B1 = 41 + (1 * 60 + 35.6874 
) / 3600.0, A12 = 1 + (49 * 60 + 43) / 3600.0, S = 80000.0;

            double ρ = 206264.806247096355,
            a = 6378245, b = 6356863.0187730473, c = 6399698.9017827110, α = 1 / 298.3, e2 = 0.006693421622966, eΒ = 0.006738525414683;//先BLA是角度，BmLmAm是弧度直接進行運算，可以直接進行三角函式運算，感覺在最上邊宣告MATH後就不用再一個一個的代用;
            L1 = L1 /180.0 * Math 
.PI;
            B1 = B1 /180.0 * Math.PI;
            A12 = (A12 / 180.0) * Math.PI;

            //轉化為弧度
            double Wm = Math.Sqrt((1 - e2 * (Math.Sin(B1)) * (Math.Sin(B1))));//在此處e2為平方
            double Nm = a / Wm;
            double Mm = a * (1 - e2) / (Wm * Wm * Wm);
            double B02 = ρ / Mm * S * Math 
.Cos(A12);
            double L02 = ρ / Nm * S * Math.Sin(A12) * (Math.Cos(B1) / Math.Sin(B1));
            double A02 = L02 * Math.Sin(B1);
            L1 = 130 + (10 * 60 + 12.2676) / 3600.0;
            B1 = 41 + (1 * 60 + 35.6874) / 3600.0;
            A12 = 1 + (49 * 60 + 43) / 3600.0;
            double Bm = B1 + 1 / 2.0 * B02 / 3600;
            double Am = A12 + 1 / 2.0 * A02 / 3600;//此時為角度

            Bm = Bm/180*Math.PI;
            Am = Am/180*Math.PI;//此時為弧度

            double L2 = ρ / Nm * S * Math.Sin(Am) * (Math.Cos(Bm) / Math.Sin(Bm)) * (1 + A02 * A02 / (24 * ρ * ρ) - B02 * B02 / (24 * ρ * ρ));
            double B2 = ρ / Mm * S * Math.Cos(Am) * (1 + L02 * L02 / 12 / (ρ * ρ) + A02 * A02 / 24 / (ρ * ρ))/*我用了一個連續相除*/;
            double A2 = ρ / Nm * S * Math.Sin(Am) * Math.Tan(Bm) * (1 + B02 * B02 / 12 / ρ / ρ + L02 * L02 / 12 / ρ / ρ * Math.Cos(Bm) * Math.Cos(Bm) + A02 * A02 / 24 / ρ / ρ);
            double i = 1;
            for (i = 1; (Math.Abs(B2 - B02) > 0.0001); i++)
            {
                Bm = B1 + 1 / 2.0 * B2 / 3600;
                Am = A12 + 1 / 2.0 * A2 / 3600;//此時為角度
                Bm = Bm / 180 * Math.PI;
                Am = Am / 180 * Math.PI;//此時為弧度

                B02 = B2;
                L02 = L2;
                A02 = A2;

                L2 = ρ / Nm * S * Math.Sin(Am) * (Math.Cos(Bm) / Math.Sin(Bm)) * (1 + A02 * A02 / (24 * ρ * ρ) - B02 * B02 / (24 * ρ * ρ));
                B2 = ρ / Mm * S * Math.Cos(Am) * (1 + L02 * L02 / 12 / (ρ * ρ) + A02 * A02 / 24 / (ρ * ρ))/*我用了一個連續相除*/;
                A2 = ρ / Nm * S * Math.Sin(Am) * Math.Tan(Bm) * (1 + B02 * B02 / 12 / ρ / ρ + L02 * L02 / 12 / ρ / ρ * Math.Cos(Bm) * Math.Cos(Bm) + A02 * A02 / 24 / ρ / ρ);
            }

            double B = B1 + B2/3600;//角度
            double L = L1 + L2 / 3600;
            double A21 = A12 + A2/3600+180;


            Console.WriteLine("{0} {1}", i , B);
            textBox1.Text = Convert.ToString(B);
            textBox3.Text = Convert.ToString(L);
            textBox5.Text = Convert.ToString(A21); 
         }

        private void Form1_Load(object sender, EventArgs e)
        {

        }

        private void button2_Click(object sender, EventArgs e)
        {
            double L1 = 130 + (10 * 60 + 12.2676) / 3600.0, B1 = 41 + (1 * 60 + 35.6874) / 3600.0;

            double ρ = 206264.806247096355,
            a = 6378245, b = 6356863.0187730473, c = 6399698.9017827110,  e2 = 0.006693421622966, eΒ = 0.006738525414683;//先BLA是角度，BmLmAm是弧度直接進行運算，可以直接進行三角函式運算，感覺在最上邊宣告MATH後就不用再一個一個的代用;
            double L2 = 130.196203938625;
            double B2 = 41.7465640175246;
            double L = L2 - L1;
            L= L / 180.0 * Math.PI;
            double Bm = (B1 + B2) / 2.0;
            double Lm = (L1 + L2) / 2.0;//此時為角度
            L1 = L1 / 180.0 * Math.PI;
            B1 = B1 / 180.0 * Math.PI;
            //轉化為弧度
            L2 = L2 / 180.0 * Math.PI;
            B2 = B2/ 180.0 * Math.PI;
            double W1 = Math.Sqrt(1 - e2 * Math.Sin(B1) * Math.Sin(B1));
            double W2 = Math.Sqrt(1 - e2 * Math.Sin(B2));

            double u1, u2;
            //Math.Sin(u1) = Math.Sin(B1) * Math.Sqrt(1 - e2) / W1;
            //Math.Sin(u2) = Math.Sin(B2) * Math.Sqrt(1 - e2) / W2;
            //Math.Cos(u1) = Math.Cos(B1) / W1;
            //Math.Cos(u2) = Math.Cos(B2) / W2;
            double a1 = Math.Sin(B1) * Math.Sqrt(1 - e2) / W1 * Math.Sin(B2) * Math.Sqrt(1 - e2) / W2;
            double a2 = Math.Cos(B1) / W1 * Math.Cos(B2) / W2;
            double b1 = Math.Cos(B1) / W1 * Math.Sin(B2) * Math.Sqrt(1 - e2) / W2;
            double b2 = Math.Sin(B1) * Math.Sqrt(1 - e2) / W1 * Math.Cos(B2) / W2;

            double δ = 0.0;
            double λ=L+δ;//這時候sitar為弧度
            double p=Math.Cos(B2)/W2*Math.Sin(λ);
            double q = b1 - b2 * Math.Cos(λ);
            double A1 = Math.Atan(p / q);
            A1 = A1 / Math.PI * 180;//此時為角度
            if (p > 0)
            {
                if (q > 0)
                {
                    A1 = Math.Abs(A1);
                }
                else
                {
                    A1 = 180 - Math.Abs(A1);
                }
            }
            else
            {
                if (q < 0)
                {
                    A1 = 180 + Math.Abs(A1);
                }
                else
                {
                    A1 = 360 - Math.Abs(A1);
                }
            }
            //Math.Sin(δ) = p * Math.Sin(A1) + q * Math.Cos(A1);
            //Math.Cos(δ) = a1 + a2 * Math.Cos(λ);
            double o = Math.Atan(p * Math.Sin(A1) + q * Math.Cos(A1) / a1 + a2 * Math.Cos(λ));//此時為弧度
            o = o / Math.PI * 180;//此時為角度
            if ((a1 + a2 * Math.Cos(λ)) > 0)
            {
                o = Math.Abs(o);
            }
            else
            {
                o = 180 - Math.Abs(o);
            }
            //Math.Cos(A0)
            //把A1化為弧度
            A1=A1/Math.PI*180;
            double zy=Math.Cos(B1) / W1*Math.Sin(A1);
            o=o/180*Math.PI;//化弧
            //Math .Sqrt(1-zy*zy);
            //sin a0=cos u1*sin A1;
            double α = (33523299 - (28189 - 70 *(1-zy*zy)) *(1-zy*zy)) * 0.0000000001;
            double Β1 = (28189 - 94 * (1-zy*zy)) * 0.0000000001;
            double x=2*a1-(1-zy*zy)*Math.Cos(δ);
            double δ1 = (α * o - Β1 * x * Math.Sin(o)) * (Math.Cos(B1) / W1 * Math.Sin(A1));
            ///////////////////////////
            ///////////////////////////
            ///////////////////////////
            double i=1;
            for(i=1;Math.Abs(δ1-δ)>0.000001;i++)
            {
                 δ = δ1;
                 λ = L + δ;
                 p = Math.Cos(B2) / W2 * Math.Sin(λ);
                 q = b1 - b2 * Math.Cos(λ);
                 A1 = Math.Atan(p / q);
                A1 = A1 / Math.PI * 180;//此時為角度
                if (p > 0)
                {
                    if (q > 0)
                    {
                        A1 = Math.Abs(A1);
                    }
                    else
                    {
                        A1 = 180 - Math.Abs(A1);
                    }
                }
                else
                {
                    if (q < 0)
                    {
                        A1 = 180 + Math.Abs(A1);
                    }
                    else
                    {
                        A1 = 360 - Math.Abs(A1);
                    }
                }
                //Math.Sin(δ) = p * Math.Sin(A1) + q * Math.Cos(A1);
                //Math.Cos(δ) = a1 + a2 * Math.Cos(λ);
                o = Math.Atan(p * Math.Sin(A1) + q * Math.Cos(A1) / a1 + a2 * Math.Cos(λ));//此時為弧度
                o = o / Math.PI * 180;//此時為角度
                if ((a1 + a2 * Math.Cos(λ)) > 0)
                {
                    o = Math.Abs(o);
                }
                else
                {
                    o = 180 - Math.Abs(o);
                }
                //Math.Cos(A0)
                //把A1化為弧度
                A1 = A1 /180*Math.PI;
                 zy = Math.Cos(B1) / W1 * Math.Sin(A1);
                o = o / 180 * Math.PI;//化弧
                //Math .Sqrt(1-zy*zy); zy=sin A0 平方  cos  (1 - zy * zy)
                //sin a0=cos u1*sin A1;
                 α = (33523299 - (28189 - 70 * (1 - zy * zy)) * (1 - zy * zy)) * 0.0000000001;
                 Β1 = (28189 - 94 * (1 - zy * zy)) * 0.0000000001;
                 x = 2 * a1 - (1 - zy * zy) * Math.Cos(δ);
                 δ1 = (α * o - Β1 * x * Math.Sin(o)) * (Math.Cos(B1) / W1 * Math.Sin(A1));
            }
            double y = ((1 - zy * zy)*(1 - zy * zy) - 2 * x * x) * Math.Cos(o);
            double A = 6356863.020 + (10708.949 - 13.474 * (1 - zy * zy)) * (1 - zy * zy);
            double B22 = 10708.938 - 17.956 * (1 - zy * zy);
            double C22 = 4.487;
            double S = A * o + (B22 * x + C22 * y) * Math.Sin(o);
            double A2 = Math.Atan(Math.Cos(B1) / W1 * Math.Sin(λ) / (b1 * Math.Cos(λ) - b2));

            textBox2.Text = Convert.ToString(S);
            textBox4.Text = Convert.ToString(A2);
            textBox6.Text = Convert.ToString(y); 
        }
    }
}

程式執行介面

這個實驗運用了兩次迭代感覺自己已經基本掌握了迭代的用法，但是不知道反算哪裡出了問題，傷心。。。不過逐漸搞懂了大一學習C語言時老師教的基礎課，需要程式的苦逼測繪人可以聯絡我，學長幫人到底。
謝謝閱讀

大地測量學白塞爾大地主題解算

namespace 大地測量學1 { public partial class Form1 : Form { public Form1() { InitializeComponent(); }

可視化n次貝塞爾曲線及過程動畫演示--大寶劍

ike all AS 2個 pat title pre while todo 先拋一個動畫模擬的一個例子，吊一吊Xing趣（4次）不夠強？再來一個這樣子，滿足你。demo說明 git倉庫地址示例我眼睛花，沒看懂，能暫停不了？可以控制動畫暫停與繼續。（供大家清

Android貝塞爾曲線————波浪效果（大波浪）

Hello大家好，很高興又一次與大家見面，今天是農曆丁酉雞年（大年初四），現在跟大家拜年有點晚，算是拜晚年，祝大家晚年幸福。這麼快大夥都到了晚年了，Android貝塞爾曲線我也準備以一個大波浪來結束。所以今天給大家帶來的就是android貝塞爾曲線製作波浪效果

Unity遊戲中使用貝塞爾曲線

str net 順序復雜讓我創建函數高程 gin 孫廣東 2015.8.15比方在3D rpg遊戲中。我們想設置彈道，不同的軌跡類型！目的:這篇文章的主要目的是要給你關於在遊戲怎樣使用貝塞爾曲線的基本想法。貝塞爾曲線是最主要的曲線，一般用

貝塞爾曲線實現的購物車添加商品動畫效果

right map 繪制開始 enter 監聽 idg 過程 protected 效果圖如下： 1.activity_main.xml <?xml version="1.0" encoding="utf-8"?> <RelativeLayout xm

把商品添加到購物車的動畫效果（貝塞爾曲線）

param from mat 位置 hold pos 開始 onclick border 目錄(?)[+] 如圖：參考： Android補間動畫，屬性動畫實現購物車添加動畫思路：確定動畫的起終點在起終點之間使用二次貝塞爾曲線填充起終點之間的點的軌跡設置屬

iOS 使用貝塞爾曲線繪制路徑

繪制直線多邊形 ini rcc 三個點 memory images arc poi 使用貝塞爾曲線繪制路徑大多數時候，我們在開發中使用的控件的邊框是矩形，或者做一點圓角，是使得矩形的角看起來更加的圓滑。但是如果我們想要一個不規則的圖形怎麽辦？有人說，叫UI

貝塞爾曲線

線性數值 int 優化連續多項式三次繪制 http 在數學的數值分析領域中，貝塞爾曲線（英語：Bézier curve）是電腦圖形學中相當重要的參數曲線。更高維度的廣泛化貝塞爾曲線就稱作貝塞爾曲面，其中貝塞爾三角是一種特殊的實例。貝塞爾曲線於1

貝塞爾曲線與CSS3動畫、SVG和canvas的應用

document cnblogs blank otto style 函數 alt one absolut 簡介貝塞爾曲線是可以做出很多復雜的效果來的，比如彈跳球的復雜動畫效果，首先加速下降，停止，然後彈起時逐漸減速的效果。使用貝塞爾曲線常用的兩個網址如下：緩動函

css貝塞爾曲線模仿餓了麽購物車小球動畫

viewport title output 代碼 put 動畫效果 doc int class 在線觀看貝塞爾曲線值：傳送門在線觀看動畫效果：傳送門代碼： <!DOCTYPE html> <html> <head> <me

canvas 繪制三次貝塞爾曲線

png 繪制 ges nload 代碼 idt head mage src 代碼： <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8">

canvas 繪制二次貝塞爾曲線

logs text lineto quad img utf 技術分享 bsp element 代碼： <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8

貝塞爾曲線原理（轉載）

image cnblogs itblog gif 二次 fcc ike 轉載依次最近在做圖形學的實驗，關於貝塞爾曲線在網上看到一個非常好的資料：以下內容轉載自：https://www.cnblogs.com/hnfxs/p/3148483.html 原理和簡單推導

【轉】貝塞爾曲線介紹

con 系列 ani 線段課本 start http wpa -- 原文鏈接： http://blog.csdn.net/sangxiaonian/article/details/51984013 http://blog.csdn.net/sangxiaonian/art

【轉】: 塞爾達組在GDC2017演講的文字翻譯：技術的智慧

ogre 遊戲 style 標準檢驗復制 bce 好用很多大家好，我是堂田卓宏，在《荒野之息》的制作中我擔任技術總監的職位。我在2003年加入任天堂，並且作為程序員參與了許多不同遊戲的制作。在本次的制作的過程中，我們的程序員團隊也需要打破許多遊戲內技術的常規。從

【轉】: 塞爾達組在GDC2017演講的文字翻譯：顯示的力量

受限聯想例子 content awl load 樣式決定時也　　塞爾達系列推出新作的時候，美術風格都有明顯變化。本作的風格比起寫實，筆觸輕快變化幅度大是其特征。2011年公開的技術演示中，畫面風格要更加寫實。最終版則更接近於卡通。5年裏到底發生了什麽呢

【轉】: 塞爾達組在GDC2017演講的文字翻譯：創新的勇氣

遊戲設計發生當我 cover ima 不能縮小正在 swipe 大家好，我是藤林秀麿，以導演的身份參與《荒野之息》的制作，感謝大家的出席。我曾經作為設計者和導演制作了諸多塞爾達遊戲（大地與時空之章、縮小帽、四支劍、幻影沙漏、天空之劍），回首望去，我已經為塞爾達系列

iOS-電池圖標【結合貝塞爾曲線控制電量顯示】

radius 顯示 ios ner att OS super 1.0 self. 基於UIView類：WKJBatteryView WKJBatteryView.h #import <UIKit/UIKit.h> @interface WKJBatter

5年時間，Java程序員如何從小白晉升為大牛

默認流行工作 jdk 任務不理解高端那是 mybatis 程序員之言?2018-06-25在程序界流行著一種默認的說法叫“黃金5年”，也就是一個程序員從入職的時候算起，前五年的選擇直接影響著整個職業生涯中的職業發展方向和薪資走向。因此如何走好這5年，徹底從一個剛入

使用貝塞爾曲線進行插值一種非常簡單的平滑多邊形的方法

ear loop line coin algorithm tro itl art mov 原文 Interpolation with Bezier Curves A very simple method of smoothing polygons Init

大地測量學白塞爾大地主題解算

相關推薦