白塞尔大地主题正算

最新推荐文章于 2023-08-05 16:29:33 发布

小白--君

最新推荐文章于 2023-08-05 16:29:33 发布

阅读量3.8k

点赞数 3

分类专栏：测绘文章标签：算法 c++

本文链接：https://blog.csdn.net/m0_57768434/article/details/122469968

版权

本文介绍了白塞尔大地主题正算的算法原理和流程。首先，阐述了基于椭球面的白塞尔投影及其在辅助球面上的解算方法。接着，详细描述了算法的具体步骤，包括输入起点坐标、角度转换、椭球参数计算、辅助函数及系数的求解、球面长度计算、经度改正以及最终的大地坐标和大地方位角的确定。整个过程涉及弧度转换、几何变换等多个数学操作。

摘要由CSDN通过智能技术生成

算法原理：

白塞尔到底主题正算公式推导的基本思想是：将椭球面上的大地元素按照白塞尔投影条件投影到辅助球面上，继而在球面上进行大地主题解算，最后再将球面上的计算结果换算到椭球面上。

算法流程：

首先输入大底线起点纬度B1，经度，大地方位角和大底线长度。其中经纬度与大地方位角角度弧度转换函数转换为弧度再参与计算。而后给定椭球参数（本次使用克氏椭球），计算起点的归化纬度，然后计算辅助函数值，计算系数A,B,C，以及α，β。然后计算球面长度，计算经度改正差数，最后计算终点大地坐标及大地方位角，判断大地方位角符号后输出结果。

#include <iostream>
#include<cmath>
using namespace std;

double angle_to_radian(double degree, double min, double second)
{
    double flag = (degree < 0) ? -1.0 : 1.0;	//判断正负
    double pi = 4 * atan(1);
    if (degree < 0)
    {
        degree = degree * (-1.0);
    }
    double angle