using?System;
using?System.Collections.Generic;
using?System.Linq;
using?System.Text;

namespace?Solution_of_geodetic_problem
{
????//白塞爾的算法
????class?Bessel
????{??/*********************************************
?????????*?正算方法
?????????*?形參為橢球選擇標數、正算起算數據
?????????*?返回值為需求算數據的數組
?????????*?******************************************/

????????public?double[]?Forward_calculate（int?ellipsoid_choose?double?B1?double?L1?double?A1?double?S）
????????{
????????????//選擇橢球參數（ellipsoid為橢球標志）
????????????double?b?=?0;
????????????double?e1?=?0;
????????????double?e2?=?0;
????????????if?（ellipsoid_choose?==?1）
????????????{
????????????????b?=?6356863.0187730473;
????????????????e1?=?0.006693421622966;
????????????????e2?=?0.006738525414683;
????????????}
????????????else?if?（ellipsoid_choose?==?2）
????????????{
????????????????b?=?6356755.2881575287;
????????????????e1?=?0.006694384999588;
????????????????e2?=?0.006739501819473;
????????????}
????????????else?if?（ellipsoid_choose?==?3）
????????????{
????????????????b?=?6356752.3142;
????????????????e1?=?0.00669437999013;
????????????????e2?=?0.00673949674227;
????????????}

????????????//計算起點的規劃緯度
????????????double?W1?=?Math.Sqrt（1?-?e1?*?Math.Sin（B1）?*?Math.Sin（B1））;
????????????double?sinu1?=?Math.Sin（B1）?*?Math.Sqrt（1?-?e1）?/?W1;
????????????double?cosu1?=?Math.Cos（B1）?/?W1;

????????????//計算輔助函數值
????????????double?sinAo?=?cosu1?*?Math.Sin（A1）;
????????????double?cota1?=?cosu1?*?Math.Cos（A1）?/?sinu1;
????????????double?sin2sgm1?=?2?*?cota1?/?（cota1?*?cota1?+?1）;
????????????double?cos2sgm1?=?（cota1?*?cota1?-?1）?/?（cota1?*?cota1?+?1）;
????????????double?cosAo_2?=?1?-?sinAo?*?sinAo;
????????????double?k2?=?e2?*?cosAo_2;
????????????double?k4?=?k2?*?k2;
????????????double?k6?=?k2?*?k2?*?k2;

????????????//計算系數A、B、C以及alf、beita；這里可能存在精度問題
????????????double?A?=?b?*?（1?+?k2?/?4?-?3?*?k4?/?64?+?5?*?k6?/?256）;
????????????double?B?=?b?*?（k2?/?8?-?k4?/?32?+?15?*?k6?/?1024）;
????????????double?C?=?b?*?（k4?/?128?-?3?*?k6?/?512）;
????????????double?alf?=?（e1?/?2?+?e1?*?e1?/?8?+?e1?*?e1?*?e1?/?16）
????????????????-?（e1?*?e1?/?16?+?e1?*?e1?*?e1?/?16）?*?cosAo_2
????????????????+?（3?*?e1?*?e1?*?e1?/?128）?*?cosAo_2?*?cosAo_2;
????????????double?bta?=?（e1?*?e1?/?32?+?e1?*?e1?*?e1?/?32）?*?cosAo_2
????????????????-?（e1?*?e1?*?e1?/?64）?*?cosAo_2?*?cosAo_2;


????????????//計算球面長度
????????????double?sgm0?=?（S?-?（B?+?C?*?cos2sgm1）?*?sin2sgm1）?/?A;
????????????double?sin2sgm1_sgm0?=?sin2sgm1?*?Math.Cos（2?*?sgm0）?+?cos2sgm1?*?Math.Sin（2?*?sgm0）;
????????????double?cos2sgm1_sgm0?=?cos2sgm1?*?Math.Cos（2?*?sgm0）?-?sin2sgm1?*?Math.Sin（2?*?sgm0）;
????????????double?sgm?=?sgm0?+?（B?+?5?*?C?*?cos2sgm1_sgm0）?*?sin2sgm1_sgm0?/?A;

????????????//計算經差改正數
????????????double?dlta?=?（alf?*?sgm?+?bta?*?（sin2sgm1_sgm0?-?sin2sgm1））?*?sinAo;


????????????//計算終點大地坐標以及大地方位角
????????????double?sinu2?=?sinu1?*?Math.Cos（sgm）?+?cosu1?*?Math.Cos（A1）?*?Math.Sin（sgm）;
?????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2017-12-25?00:22??Solution_of_geodetic_problem\
?????文件???????10038??2017-10-22?17:58??Solution_of_geodetic_problem\Bessel.cs
?????文件???????45168??2017-10-22?18:58??Solution_of_geodetic_problem\Form1.Designer.cs
?????文件????????6274??2017-10-22?18:58??Solution_of_geodetic_problem\Form1.cs
?????文件????????5817??2017-10-22?18:58??Solution_of_geodetic_problem\Form1.resx
?????文件?????????509??2017-10-22?15:55??Solution_of_geodetic_problem\Program.cs
?????目錄???????????0??2017-12-25?00:22??Solution_of_geodetic_problem\Properties\
?????文件????????1388??2017-10-22?15:55??Solution_of_geodetic_problem\Properties\AssemblyInfo.cs
?????文件????????2908??2017-10-22?15:55??Solution_of_geodetic_problem\Properties\Resources.Designer.cs
?????文件????????5612??2017-10-22?15:55??Solution_of_geodetic_problem\Properties\Resources.resx
?????文件????????1115??2017-10-22?15:55??Solution_of_geodetic_problem\Properties\Settings.Designer.cs
?????文件?????????249??2017-10-22?15:55??Solution_of_geodetic_problem\Properties\Settings.settings
?????文件????????3792??2017-10-22?19:03??Solution_of_geodetic_problem\Solution_of_geodetic_problem.csproj
?????文件?????????897??2017-10-22?15:55??Solution_of_geodetic_problem\Solution_of_geodetic_problem.sln
?????文件???????29184??2017-12-18?23:43??Solution_of_geodetic_problem\Solution_of_geodetic_problem.suo
?????目錄???????????0??2017-12-25?00:22??Solution_of_geodetic_problem\bin\
?????目錄???????????0??2017-12-25?00:22??Solution_of_geodetic_problem\bin\Debug\
?????文件???????28160??2017-12-24?23:52??Solution_of_geodetic_problem\bin\Debug\Solution_of_geodetic_problem.exe
?????文件???????42496??2017-12-24?23:52??Solution_of_geodetic_problem\bin\Debug\Solution_of_geodetic_problem.pdb
?????文件???????11600??2017-12-24?23:49??Solution_of_geodetic_problem\bin\Debug\Solution_of_geodetic_problem.vshost.exe
?????文件?????????490??2016-07-16?19:44??Solution_of_geodetic_problem\bin\Debug\Solution_of_geodetic_problem.vshost.exe.manifest
?????目錄???????????0??2017-12-27?20:40??Solution_of_geodetic_problem\bin\Release\
?????目錄???????????0??2017-12-25?00:22??Solution_of_geodetic_problem\obj\
?????目錄???????????0??2017-12-25?00:22??Solution_of_geodetic_problem\obj\x86\
?????目錄???????????0??2017-12-25?00:22??Solution_of_geodetic_problem\obj\x86\Debug\
?????文件????????5420??2017-11-04?17:18??Solution_of_geodetic_problem\obj\x86\Debug\DesignTimeResolveAssemblyReferences.cache
?????文件????????6314??2017-12-24?23:52??Solution_of_geodetic_problem\obj\x86\Debug\DesignTimeResolveAssemblyReferencesInput.cache
?????文件?????????180??2017-12-18?23:34??Solution_of_geodetic_problem\obj\x86\Debug\Solution_of_geodetic_problem.Form1.resources
?????文件?????????180??2017-12-18?23:34??Solution_of_geodetic_problem\obj\x86\Debug\Solution_of_geodetic_problem.Properties.Resources.resources
?????文件?????????810??2017-12-24?23:49??Solution_of_geodetic_problem\obj\x86\Debug\Solution_of_geodetic_problem.csproj.FileListAbsolute.txt
?????文件?????????975??2017-10-22?18:58??Solution_of_geodetic_problem\obj\x86\Debug\Solution_of_geodetic_problem.csproj.GenerateResource.Cache
............此處省略6個文件信息