大地坐标和高斯平面坐标转换

最新推荐文章于 2022-04-23 15:55:52 发布

zscomehuyue

最新推荐文章于 2022-04-23 15:55:52 发布

阅读量4.6k

点赞数 1

文章标签：算法编程 C C++ C#

本文链接：https://blog.csdn.net/zscomehuyue/article/details/83681692

版权

本文详细介绍了大地坐标与高斯平面坐标之间的转换公式和计算方法，包括图幅理论面积计算、椭球面上任意梯形面积计算以及高斯投影反解变换模型。还提供了算法伪代码，涉及的数学公式和参数定义，适用于编程实现坐标转换和面积计算。

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

图幅理论面积与图斑椭球面积计算公式及要求

一、图幅理论面积计算公式
（1）
式中：
a—椭球长半轴(单位：米)，α—椭球扁率，b—椭球短半轴(单位：米)。
е²﹦(a²﹣b²)/a²。
A﹦1﹢（3/6）е²﹢（30/80）е4﹢（35/112）е6﹢（630/2304）е8。
B﹦ （1/6）е²﹢（15/80）е4﹢（21/112）е6﹢（420/2304）е8。
C﹦ （3/80）е4﹢ （7/112）е6﹢（180/2304）е8。
D﹦ （1/112）е6﹢ （45/2304）е8。
E﹦ （5/2304）е8。
ΔL—图幅东西图廓的经差（单位：弧度）。
(B2﹣B1)—图幅南北图廓的纬差（单位：弧度），Bm﹦（B1﹢B2）/2。

二、椭球面上任意梯形面积计算公式
（2）
其中：A,B,C,D,E 为常数，按下式计算:
е²﹦(a²﹣b²)/a²
A﹦1﹢（3/6）е²﹢（30/80）е4﹢（35/112）е6﹢（630/2304）е8
B﹦ （1/6）е²﹢（15/80）е4﹢（21/112）е6﹢（420/2304）е8
C﹦ （3/80）е4﹢ （7/112）е6﹢（180/2304）е8
D﹦ （1/112）е6﹢（45/2304）е8
E﹦ （5/2304）е8
式中：a—椭球长半轴(单位：米)，b—椭球短半轴(单位：米)；
ΔL—图块经差(单位：弧度)； (B2﹣B1)—图块纬差(单位：弧度)
Bm﹦（B1﹢B2）/2。

[color=red]三、高斯投影反解变换（）模型
（若坐标不带带号，则不需减去带号×1000000；）

+中央子午线经度值（孤度）（3）
式中：

公式说明：若坐标为没有带号前缀格式，则不需减去带号×1000000；若坐标为有带号前缀格式，则需减去带号×1000000。

四、计算用到的常数、椭球参数
在计算图幅理论面积