大地主题解算-C#

格式：doc
大小：56.50 KB
文档页数：5

下载文档原格式

/ 5

84坐标系大地主题解算

84坐标系大地主题解算
84坐标系是一种常用的大地坐标系，用于地理测量和定位。

大
地主题解算是指通过观测数据和数学模型，计算出地球上某一点的
大地坐标。

大地主题解算涉及到以下几个方面的内容：
1. 大地椭球模型，大地主题解算基于大地椭球模型，即将地球
近似看作一个椭球体。

常用的大地椭球模型有WGS84、CGCS2000等。

这些模型定义了地球的几何形状和参数，如椭球体的长半轴、扁率等。

2. 观测数据，大地主题解算需要借助于观测数据，包括卫星导
航系统（如GPS）的测量数据、测角仪、水准仪等测量设备的观测
结果。

这些数据包括卫星的位置、接收器的观测值等。

3. 大地测量基本理论，大地主题解算基于大地测量的基本理论，如三角测量、水准测量、重力测量等。

这些理论提供了测量数据的
处理方法和数学模型，用于计算点的坐标和高程。

4. 大地主题解算方法，大地主题解算采用不同的方法，如最小二乘法、无约束最小二乘法、最大似然估计等。

这些方法根据具体的观测数据和模型，选择合适的数学模型和计算方法，以求得最优的解。

5. 大地坐标系统，大地主题解算的结果是点的大地坐标，通常包括经度、纬度和大地高。

经度表示点在东西方向上的位置，纬度表示点在南北方向上的位置，大地高表示点相对于参考椭球体的高度。

综上所述，大地主题解算是一种通过观测数据和数学模型计算地球上某一点的大地坐标的方法。

它涉及到大地椭球模型、观测数据、大地测量基本理论、解算方法和大地坐标系统等方面的内容。

通过大地主题解算，可以实现地理测量和定位的精确计算和分析。

大地测量学重要名词解释简答题

大地测量学基础一、名词解释1、大地测量学：是指在一定的时间与空间参考系中，测量和描绘地球形状及其重力场并监测其变化，为人类活动提供关于地球的空间信息的一门学科。

2、天球：是指以地球质心O（或测站）为中心，半径r为任意长度的一个假想的球体。

3、大地基准：指用以描述地球形状的参考椭球的参数，以及参考椭球在空间中的定位及定向，还有在描述这些位置时所采用的单位长度的定义。

4、岁差：地球绕地轴旋转，由于日、月等天体的影响，地球的旋转轴在空间围绕黄级发生缓慢移动。

5、章动：地球旋转轴在岁差的基础上叠加18.6年的短期周圆周运动，振幅为9.21秒，这种现象称为章动。

6、极移：地球自转使地球体自身内部结构的相对位置变化，从而导致极点在地球表面上的位置随时间而变化，这种现象被称为极移。

7、恒星时(ST)：以春分点作为基本参考点，由春分点周日视运动确定的时间，称为恒星时。

8、真太阳时MT：以真太阳作为基本参考点，由其周日视运动确定的时间，称为真太阳时。

一个真太阳日就是真太阳连续两次经过某地的上中天（上子午圈）所经历的时间。

9、大地水准面：假想海洋处于完全静止的平衡状态时海水面延伸到大陆地面以下所形成的闭合曲面，叫大地水准面。

10、正常椭球：与地球质量相等且质量分布均匀的椭球。

11、正常重力加速度：正常椭球对其表面与外部点所产生的重力加速度。

12、正常位水准面：相应于正常重力加速度的重力等位面。

13、理论闭合差：由水准面不平行而引起的水准环线闭合差，称为理论闭合差。

14、正常椭球面：是大地水准面的规则形状（一般指旋转椭球面）。

因此引入正常椭球后，地球重力位被分成正常重力位和扰动位两部分，实际重力也被分成正常重力和重力异常两部分。

15、总的地球椭球：一个和整个大地体最为密合的。

总地球椭球中心和地球质心重合，总的地球椭球的短轴与地球地轴相重合，起始大地子午面和起始天文子午面重合，总地球椭球和大地体最为密合。

16、参考椭球：具有确定参数(长半径 a和扁率α)，经过局部定位和定向，同某一地区大地水准面最佳拟合的地球椭球。

用C语言实现大地主题解算

６
４
６
４
（争寺＋
‘ ）ｃｏ５４Ａ ’
裴连磊ＰＥＩＬｉａｎ — ｌｅｉ
（新疆地矿局测绘大队，乌鲁木齐８３００１７）（ＸｉｎｊｉａｎｇＧｅｏｌｏｇｙａｎｄＭｉｎｅｒａｌＢｕｒｅａｕ，ＳｕｒｖｅｙｉｎｇａｎｄＭａｐｐｉｎｇＢｒｉｇａｄｅ，Ｕｒｕｍｑｉ８３００１７，Ｃｈｉｎａ）
１－５计算经度差改正数
本文的程序实现是采用的白塞尔大地主题解算的方法，根据白塞尔大地主题解算的方法，分析得出适合计算机编程的大地主题正反算的具体实现步骤。１白塞尔法大地主题正算思想已知量：大地线起点的纬度Ｂ，经度Ｌ，大地方位角Ａ１及大地线长度Ｓ。求解量：大地线终点的纬度Ｂ：，经度Ｌ２
ｃ番一番＋．．
ｃ＿ｂ（薷４一６一）
２
＝
ｃ０ｓｕ２＝
Ｌａ１＝ｓｉｎｕｌｓｉｎｕ２Ｉａ２＿ｃｏｓｕｌｃ０ｓｕ２ｂ＝
ＣＯＳＵｌｓｉｎｕ２，ｂ２＝ｓｉｎｕｌＣＯＳＵ２
之值。即下式：
求解量：大地线长度Ｓ及起、终点处的大地方位角Ａ
及Ａ＇。
Ａ：ｂ（１一鲁十击ｋ６＿－
２．１辅助计算ｗ：

大地主题解算

"

S2 2 2 2 2 4 A" A21 A12 S sin Am t m 1 [ cos A ( 2 7 9 t 5 m m m m m) 2 Nm 24N m 2 2 sin 2 Am t m 2 t m 2 m ] 5次项
本节主要内容
• • • • • 大地测量主题解算的一般说明勒让德级数高斯平均引数正算公式高斯平均引数反算公式高斯平均引数正、反算公式的实现
4.7 大地测量主题解算概述
一、大地测量主题解算的一般说明大地元素：椭球面上点的大地经度L，大地纬度B，两点间的大地线长度S，及其正反大地方位角A12，A21。大地主题解算：已知某些大地元素推求另一些大地元素叫做大地主题解算；包括正解和反解。
大地线长度和正反方位角.S sin Am，S cosAm 及A" 计算公式为：
计算出S sin Am，S cosAm及A"后，按下式计算大地线长度和正反方位角：
S sin Am t an Am S cos A m S S sin Am S cos Am sin Am cos Am 1 A12 Am A"， 2 1 A21 Am A&2 2 3 m cos2 Am 1 t m m 4t m m ] 5次项

L" L2 L1
"
S2 2 2 S sec Bm sin Am 1 [ sin A t m m 2 Nm 24N m 2 2 2 cos2 Am 1 m 9t m m ] 5次项
Am是大地方位角，其值在 0－360之间，设b B2 B1 , l L2 L1 , b看成x, l看成y，平面方位角的象限判别一样。

大地主题解算深度干货超精

大地主题解算深度干货超精在我们生活的世界上，大地是我们所处的基本环境。

它是被压迫、被掠夺的，也是我们努力保护和改善的对象。

深度解算是一种精确的测量和计算方法，可以帮助我们更好地了解和利用大地资源。

本文将介绍大地主题解算和其带来的深度干货。

一、什么是大地主题解算？大地主题解算是一种基于卫星观测数据和大地测量学原理的技术，用于精确计算大地物体的三维位置和形态变化。

通过对不同时刻的卫星图像进行分析和计算，可以得到地面物体的精确坐标、高程等信息。

大地主题解算的核心思想是通过计算大地物体的三维位置和形态变化，来揭示地球表面的变化过程。

这种技术可以用于监测地表的沉降、地壳的变形、建筑物的结构和变化等。

它在测绘、地质灾害监测、城市规划等领域具有广泛的应用前景。

二、大地主题解算的应用领域1. 地质灾害监测地质灾害是世界各地普遍存在的问题，如地震、山体滑坡、地裂缝等。

通过大地主题解算技术，可以实时监测和预测地质灾害的发生和演变趋势。

这对于提前采取相应的防灾措施、保护人民的生命财产安全具有重要意义。

2. 城市规划随着城市化进程的不断加快，城市规划也变得愈发重要。

大地主题解算可以为城市规划提供精确的地面信息，包括土地利用、道路交通、建筑物布局等。

这有助于提高城市规划的科学性和有效性，减少资源浪费和环境污染。

3. 地表沉降监测地表沉降是地下水开采、矿井开采等人类活动造成的一个重要问题。

通过大地主题解算技术，可以实时监测地表的沉降情况，并对地下水开采和矿井开采等活动进行合理调整和管理。

这有助于减少地表沉降对城市建设和生态环境的不利影响。

4. 环境保护大地主题解算可以为环境保护提供准确的数据支持。

例如，通过对森林面积、湿地面积等进行监测和计算，可以及时发现和预防环境变化和破坏。

这对于保护生态环境和维护生态平衡具有重要意义。

三、大地主题解算的优势和挑战大地主题解算作为一种先进的测量和计算技术，具有很多优势。

首先，它可以实时获取和分析大地物体的位置和形态变化，具有高精度和高时效性。

白塞尔大地主题解算的基本思想

白塞尔大地主题解算的基本思想
首先，白塞尔大地主题解算的基本思想之一是建立地壳变形的力学模型。

地壳变形是地球表面的一项重要现象，是由于地质作用和地球内部力
学过程的结果。

地壳变形的力学模型是研究地壳变形的重要方法。

常见的
地壳变形的力学模型有弹性模型、弹塑性模型和拟弹性模型等。

这些模型
可以描述地壳的变形过程，通过对地壳的变形过程进行建模，可以更好地
理解地壳变形的机制和动力学过程。

其次，白塞尔大地主题解算的思想之一是研究现今地球表面的动力学
过程。

地球表面的动力学过程包括板块运动、地震活动、火山喷发等。

这
些过程在地球的长期演化中起着重要的作用。

通过对现今地球表面动力学
过程的研究，可以揭示地球内部的结构和动力学机制，进而更好地理解地
壳变形的成因和发展变化的规律。

第三，白塞尔大地主题解算的基本思想之一是研究地壳变形的成因。

地壳变形的成因是地壳运动的基本原因，也是地球科学研究的一个重要问题。

地壳变形的成因包括构造运动、地壳应力状态的改变、地震活动等。

通过研究地壳变形的成因，可以更好地了解地壳变形的机制和规律，进而
为地震预测和地壳运动的控制提供科学依据。

总之，白塞尔大地主题解算的基本思想是通过建立地壳变形的力学模型，研究现今地球表面的动力学过程，探究地壳变形的成因，揭示地壳变
形与地球动力学过程之间的相互关系，以推动地壳运动的理论和应用研究。

这一思想对于研究地壳运动的机制和规律，了解地壳变形的成因和动力学
过程，提高地震预测和地壳运动控制的水平具有重要意义。

大地测量学基础复习资料

大地测量学基础一、填空题：1、时间的计量包括时间原点和度量单位（尺度）两个元素。

坐标的计量包括坐标原点、坐标轴的指向和坐标的尺度三个元素。

2、测量外业工作的基准线是铅垂线，基准面是大地水准面。

在椭球面上进行大地测量计算的基准线是法线，基准面是椭球面。

3、经纬仪十字丝分划板上丝和下丝的作用是测量视距。

4、衡量精度的指标有中误差、极限误差、或然误差、平均误差、相对误差。

5、过椭球面上一点P 的垂线与赤道面的夹角称为大地纬度，椭球面上一点P 与椭球中心的连线与赤道面的夹角称为地心纬度，在过椭球面上一点P 的子午面上，以椭圆中心O 为圆心,以椭球长半径a 为半径做辅助圆，反向延长过P 点并与x 轴垂直的垂线，与辅助圆交于P 1点，则P 1与椭球中心的连线与赤道面的夹角称为归化纬度，符号q= BBN B M 0cos d 表示等量纬度。

6、某直线的方位角为123°20’，该直线的反方位角为303°20’。

已知P 1点坐标（-2，-2），P 2点坐标（-4，-4），则P 1P 2的方位角为225°，P 2P 1的方位角为45°。

【注释】在同一高斯平面直角坐标系内一条直线的正、反坐标方位角相差180°，即：α12=α21±180°。

（详见数字测图课本23页）7、水准路线按布设形式分为闭合水准路线、附合水准路线和支水准路线。

8、高斯投影属于横轴椭圆柱等角投影，保证了投影的角度不变性，图形的相似性，以及在某点方向上的长度比的同一性。

在高斯平面直角坐标系中，中央子午线的投影为坐标x 轴。

9、旋转椭球的形状和大小是由子午椭圆的5个基本几何参数来决定的，他们分别是长半轴a 、短半轴b 、扁率、第一偏心率、第二偏心率。

两个互相垂直的法截弧的曲率半径，在微分几何中统称为主曲率半径，它们是指子午圈曲率半径和卯酉圈曲率半径，椭球面上任意一点的平均曲率半径R 等于该点的子午圈曲率半径和卯酉圈曲率半径的几何平均值。

控制测量复习题以及参考答案

《控制测量学》试题参考答案一、名词解释：1、子午圈：过椭球面上一点的子午面同椭球面相截形成的闭合圈。

2、卯酉圈：过椭球面上一点的一个与该点子午面相垂直的法截面同椭球面相截形成的闭合的圈。

3、椭园偏心率：第一偏心率e = 4。

2-b2第二偏心率e,= V。

2-b24、大地坐标系：以大地经度、大地纬度和大地高来表示点的位置的坐标系。

P35、空间坐标系：以椭球体中心为原点，起始子午面与赤道面交线为X轴，在赤道面上与X轴正交的方向为Y轴，椭球体的旋转轴为Z轴，构成右手坐标系O-XYZ。

P46、法截线：过椭球面上一点的法线所作的法截面与椭球面相截形成圈。

P97、相对法截线：设在椭球面上任意取两点A和B，过A点的法线所作通过B点的法截线和过B点的法线所作通过A点的法截线，称为AB两点的相对法截线。

P158、大地线：椭球面上两点之间的最短线。

9、垂线偏差改正：将以垂线为依据的地面观测的水平方向观测值归算到以法线为依据的方向值应加的改正。

P1810、标高差改正：由于照准点高度而引起的方向偏差改正。

P19 11、截面差改正：将法截弧方向化为大地线方向所加的改正。

P20 12、起始方位角的归算：将天文方位角以测站垂线为依据归算到椭球面以法线为依据的大地方位角。

P2213、勒让德尔定理：如果平面三角形和球面三角形对应边相等，则平面角等于对应球面角减去三分之一球面角超。

P2714、大地元素：椭球面上点的大地经度、大地纬度，两点之间的大地线长度及其正、反大地方位角。

P2815、大地主题解算：如果知道某些大地元素推求另外一些大地元素，这样的计算称为大地主题解算。

P2816、大地主题正算：已知P点的大地坐标，P至P的大地线长及其大地方位角，计算P点的大地坐标和大地线在P 点的反方位角。

1 / 1217、大地主题反算：如果已知两点的大地坐标，计算期间的大地线长度及其正反方位角。

18、地图投影：将椭球面上各个元素（包括坐标、方向和长度）按一定的数学法则投影到平面上。

Bessel大地主题解算程序

// 计算终点大地坐标及方位角
B2,L2,A2
sinu2=sinu1*cos(sigma)+cosu1*cos(A1)*sin(sigma);
B2=atan(sinu2/sqrt((1-e2)*(1-sinu2*sinu2)));
//计算 B2
lambda=atan(sin(A1)*sin(sigma)/(cosu1*cos(sigma)-sinu1*sin(sigma)*cos(A1))); //求 λ
大地主题解算的意义bessel函数bessel卫星坐标解算程序基线解算besseljmatlabbesselgps基线解算原理gps基线解算南方gps解算软件
#include<stdio.h> #include<math.h> #include<stdlib.h>
double trans1() { double B1,B11,B12,B13,B111;
//计算归化纬度
double sinA0,cotsigma1,sin2sigma1,cos2sigma1; //计算辅助三角函数值 sinA0=cosu1*sin(A1); cotsigma1=cosu1*cos(A1)/sinu1; sin2sigma1=2*cotsigma1/(1+cotsigma1*cotsigma1); cos2sigma1=(cotsigma1*cotsigma1-1)/(cotsigma1*cotsigma1+1);
//计算 A,B,C 以及 α，β 的值 //P144(4-265)
//P146(4-284)
（好像有问
double sigma0,sin2sigma1sigma0,cos2sigma1sigma0,sigma; //计算球面长度 σ sigma0=(S-(B+C*cos2sigma1)*sin2sigma1)/A; sin2sigma1sigma0=sin2sigma1*cos(2*sigma0)+cos2sigma1*sin(2*sigma0); cos2sigma1sigma0=cos2sigma1*cos(2*sigma0)-sin2sigma1*sin(2*sigma0); sigma=sigma0+(B+5*C*cos2sigma1sigma0)*sin2sigma1sigma0/A;

大地主题解算方法综述

2) 以勒让德级数为基础的公式中, 最典型的代表应该
是高斯平均引数公式 [2O4] 。高斯 ( G auss)提出把勒让德级数式改化成以 P1、P2 ( P1, P 2 为大地椭球面上的两点 ) 两点的平均纬度和平均方位角为根据。这样勒让德级数式中所
第 4期
周振宇等大地主题解算方法综述
时, 公式还可以得到简化。 3) 用龙格 ) 库塔 ( R ungeOK utta) 法 [ 5] 解大地线微分方
程。它的实质是用若干点的函数值的线性组合, 代替泰勒
级数展开中的导数计算, 又可保持必要的精度。对于短距
离大地主题解算, 可以只取一步; 如果距离较长, 可以分
( L1, B1 ) 和 P 2 ( L2, B2 ) 点间的大地线弧长 S 积分得:
Q B2 - B1 =
p2 cosAdS p1 M
Q L2 - L1 =
p2 p1
N
s inA cosB
dS
( 2)
Q A 2 - A 1 ? 180b=
p2 tanB sinA dS
p1
N
在初等函数中这些积分不能计算, 所以其精确值不能
椭球的过渡。
白塞尔 ( Besse l) 首先提出并解决了投影条件, 使这一解
法得以实现。
这类公式的特点是计算公式展开 e2 或的幂级数, 解算
精度与距离长短无关。因此它既适用于短距离解算, 也适用于长距离解算 [ 5] 。其主要缺点在于: 由 S 求 R、由 L 求 K, 或相反的运算, 需要进行迭代。同时还要预先计算辅助
种, 其中大部分适用于短距离, 一部分适用于中距离, 只有几种适用于长距离。这几十种方法虽然形形色色 , 各有不同, 但就其推导的理论基础来说, 大致可归纳为以下五类。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

大地主题解算-正算-C#
大地主题解算-正算-程序
using System;
using System.Collections.Generic;
using ponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Windows.Forms;
namespace WindowsFormsApplication2
{
public partial class Form1 : Form
{
public Form1()
{
InitializeComponent();
}
private void button1_Click(object sender, EventArgs e) {
int n,du,fen,miao;
double B1, L1, A1, S, dS, dB1, dL1, M1, N1, W1;
double E = 0.006694384999588, a = 6378140, C;
B1 =( Convert.ToDouble(textBox1.Text) + Convert.ToDouble(textBox9.Text)/60 + Convert.ToDouble(textBox10.Text)/3600 )* Math.PI / 180;
L1 = (Convert.ToDouble(textBox2.Text) + Convert.ToDouble(textBox11.Text)/60 + Convert.ToDouble(textBox12.Text)/3600) * Math.PI / 180; ;
A1 = (Convert.ToDouble(textBox3.Text) + Convert.ToDouble(textBox13.Text)/60 + Convert.ToDouble(textBox14.Text)/3600) * Math.PI / 180; ;
S = Convert.ToDouble(textBox5.Text);
n = Convert.ToInt16(textBox4.Text);
double[] L = new double[100000];
double[] B = new double[100000];
double[] dL = new double[100000];
double[] dB = new double[100000];
double[] A = new double[100000];
double[] W = new double[100000];
double[] N = new double[100000];
double[] M = new double[100000];
W1 = Math.Pow(1 - E * Math.Pow(Math.Sin(B1), 2), 0.5);
N1 = a / W1;
M1 = a * (1 - E) / (W1 * W1 * W1);
C = N1 * Math.Sin(A1) * Math.Cos(B1);
dS = S / n;
dB1 = dS * Math.Cos(A1) / M1;
dL1 = dS * Math.Sin(A1) / N1 / Math.Cos(B1);
B[0] = B1;
L[0] = L1;
dB[0] = dB1;
dL[0] = dL1;
A[0] = A1;
M[0] = M1;
N[0] = N1;
W[0] = W1;
int i;
for (i = 1; i < n; i++)
{
dB[i - 1] = dS * Math.Cos(A[i - 1]) / M[i - 1];
B[i] = B[i - 1] + dB[i - 1];
W[i] = Math.Pow(1 - E * Math.Pow(Math.Sin(B[i]), 2), 0.5);
N[i] = a / W[i];
M[i] = a * (1 - E) / (W[i] * W[i] * W[i]);
dL[i - 1] = dS * Math.Sin(A[i - 1]) / (N[i - 1] * Math.Cos(B[i - 1])); L[i] = L[i - 1] + dL[i - 1];
A[i] = Math.Asin(C / (Math.Cos(B[i]) * N[i]));
}
du = Convert.ToInt16(Math.Floor(B[i - 1] * 180 / Math.PI));//du
fen = Convert.ToInt16(Math.Floor((B[i - 1] * 180 / Math.PI - Math.Floor(B[i - 1] * 180 / Math.PI)) * 60)); //fen
miao = Convert.ToInt16(Convert.ToInt16(((B[i - 1] * 180 / Math.PI - Math.Floor(B[i - 1] * 180 / Math.PI)) * 60 - Math.Floor((B[i - 1] * 180 / Math.PI - Math.Floor(B[i - 1] * 180 / Math.PI)) * 60)) * 60));//秒?
textBox6.Text = Convert.ToString(du) + "度è" + Convert.ToString(fen) + "分?" + Convert.ToString(miao) + "秒?";
du = Convert.ToInt16(Math.Floor(L[i - 1] * 180 / Math.PI));//du
fen = Convert.ToInt16(Math.Floor((L[i - 1] * 180 / Math.PI - Math.Floor(L[i - 1] * 180 / Math.PI)) * 60)); //fen
miao = Convert.ToInt16(Convert.ToInt16(((L[i - 1] * 180 / Math.PI - Math.Floor(L[i - 1] * 180 / Math.PI)) * 60 - Math.Floor((L[i - 1] * 180 / Math.PI - Math.Floor(L[i - 1] * 180 / Math.PI)) * 60)) * 60));//秒?
textBox7.Text = Convert.ToString(du) + "度è" + Convert.ToString(fen) + "分?" + Convert.ToString(miao) + "秒?";
du = Convert.ToInt16(Math.Floor((180 + A[i - 1] * 180 / Math.PI)));//du
fen = Convert.ToInt16(Math.Floor(((180 + A[i - 1] * 180 / Math.PI) - Math.Floor((180 + A[i - 1] * 180 / Math.PI))) * 60)); //fen
miao = Convert.ToInt16(Convert.ToInt16((((180 + A[i - 1] * 180 / Math.PI) - Math.Floor((180 + A[i - 1] * 180 / Math.PI))) * 60 - Math.Floor(((180 + A[i - 1] * 180 / Math.PI) - Math.Floor((180 + A[i - 1] * 180 / Math.PI))) * 60)) * 60));//秒?
textBox8.Text = Convert.ToString(du) + "度è" + Convert.ToString(fen) + "分?" + Convert.ToString(miao) + "秒?";
}
private void Form1_Load(object sender, EventArgs e)
{
} }
}。