chap5 地球椭球与测量计算

地球椭球与测量计算(2)大地测量学基础·第五章Contents课程目录大地测量学核心理论框架，从椭球面几何到大地问题解算的完整知识体系。01椭球面法截线与曲率半径02大地线理论与微分方程03地面观测值归算至椭球面04椭球面大地问题解算方法CHAPTER01椭球面法截线与曲率半径理解椭球面几何性质的基础——卯酉圈、子午圈与任意方向法截线Geodesy·NormalSection法截线的定义与分类法截线是包含椭球面法线的平面与椭球面的交线，通过一点可作无数条法截线。子午圈和卯酉圈是两条最基本的法截线，其曲率半径决定了椭球面的局部几何形态。01法截线定义：包含椭球面上某点法线的平面与椭球面相截所得的曲线，是研究椭球面几何性质的基本元素。02子午圈：包含椭球短轴的法截线，平面经过椭球中心，是椭球面上唯一的平面曲线族。03卯酉圈：与子午圈正交的法截线，平面垂直于子午面，在赤道处为圆、其他纬度为椭圆。04任意方向法截线：介于子午圈和卯酉圈之间，曲率半径由欧拉公式确定。椭球面法截线几何关系·子午圈与卯酉圈示意GEODESY·大地测量学卯酉圈曲率半径卯酉圈曲率半径N=a/W是大地测量的基本参数，它随纬度变化而变化：赤道处最小等于a，极点处最大等于c。N的物理意义是椭球面上某点沿东西方向的弯曲程度度量，是高斯投影计算的核心参数。公式表达N=a/W=a/√(1−e²sin²B)，其中a为椭球长半轴，e为第一偏心率，B为大地纬度。FORMULA纬度变化规律赤道处（B=0°）N=a=6378137m；极点处（B=90°）N≈6399594m，随纬度单调递增。VARIATION几何意义等于该点到椭球短轴的垂直距离在法线方向的投影，反映椭球面东西方向的弯曲程度。GEOMETRY实用价值是高斯投影坐标计算、大地线微分方程推导、距离归算等多项大地测量计算的基础参数。APPLICATIONGEODETICFUNDAMENTALS子午圈曲率半径子午圈曲率半径M=a(1-e²)/W³是椭球面上南北方向弯曲程度的度量，其变化幅度比卯酉圈更大。FORMULA公式定义M=a(1−e²)/W³=a(1−e²)/(1−e²sin²B)3/2，分母为W的三次方，变化比N更剧烈。W³三次方EXTREMES极值特征赤道处M=a(1−e²)≈6335553m为最小值；极点处M=c=a²/b≈6399594m为最大值。Δ64kmCOMPARISON与卯酉圈比较同一纬度处M始终小于N（极点除外），两者差异在赤道最大约43km，反映椭球各向异性。M<NAPPLICATION弧长计算应用子午线弧长微元dX=M·dB，从赤道到纬度B的子午线弧长需对M积分，是大地测量基本计算。dX=M·dBGEODESY·WGS84曲率半径的纬度变化特征WGS84椭球的曲率半径随纬度单调递增，子午圈曲率半径M的变化幅度(约64km)远大于卯酉圈N(约21km)。两者在赤道处差异最大，在极点处汇合为极曲率半径c，这一规律深刻反映了地球椭球的扁率效应。WGS84椭球曲率半径随纬度变化表纬度B子午圈曲率半径M(m)卯酉圈曲率半径N(m)N-M差值(m)0°6,335,5536,378,13742,58415°6,339,6026,379,57439,97230°6,351,3776,383,46132,08445°6,367,3806,388,83921,45960°6,383,4546,394,20310,74975°6,395,2796,398,1662,88790°6,399,5946,399,5940曲率半径M和N均随纬度单调递增，N-M差值从赤道42.6km递减至极点0，体现椭球扁率效应GEODESY·CURVATUREANALYSIS任意方向法截线曲率半径（欧拉公式）欧拉公式RA=1/(cos²A/M+sin²A/N)统一描述了椭球面上任意方位角方向法截线的曲率半径，它是子午圈曲率半径M和卯酉圈曲率半径N按方位角的加权调和平均，是大地线曲率分析的理论基础。欧拉公式表达1/RA=cos²A/M+sin²A/N，其中A为法截线的大地方位角，从子午线北方向顺时针量取。方位角A特殊方向验证当A=0°或180°时RA=M（子午圈方向）；当A=90°或270°时RA=N（卯酉圈方向）。0°/90°曲率半径范围任意方向的RA满足M≤RA≤N，子午圈方向曲率最大（弯曲最厉害），卯酉圈方向曲率最小。M≤RA≤N平均曲率半径R=√(M·N)是子午圈和卯酉圈的几何平均值，常用于球面近似计算中替代椭球面。√(M·N)Geodesy·EllipsoidGeometry平行圈（纬度圈）曲率半径平行圈是椭球面上平行于赤道的纬线圆，其半径r=N·cosB随纬度递增而单调递减。平行圈不是法截线（赤道除外），但其半径在高斯投影坐标计算、子午收敛角推导和经度方向距离换算中扮演关键角色。01平行圈定义—椭球面与平行于赤道面的平面的交线，是一个圆，半径r=N·cosB=a·cosB/√(1−e²sin²B)r=N·cosB02非等位面特性—平行圈不是法截线（赤道除外），因除赤道外平行面不包含椭球面法线，法线与平行面存在夹角非等位面03纬度变化规律—赤道处r=a=6,378,137m（最大），极点处r=0，中纬度约4,000–5,000km，单调递减6378137m→004测量应用价值—平行圈半径r决定经度差对应的东西方向弧长（dL方向位移=r·dL），是高斯投影计算的核心参数dL=r·dλChapter02大地线理论与微分方程椭球面上的最短路径——大地线的定义、性质与克莱劳方程GEODESY大地线的定义与几何特性大地线是椭球面上两点间的最短路径，其主法线与椭球面法线处处重合，是大地测量中距离量测和方向传递的唯一正确基准线。大地测量野外作业现场01微分几何定义：大地线是曲面上主法线与曲面法线处处重合的曲线，即测地曲率为零的曲线。02最短路径性质：椭球面上两点间的大地线是所有连接曲线中长度最短的一条。03相对法截线差异：正法截线与反法截线一般不重合，大地线是唯一确定的最短路径。04实际测量意义：长距离三角与导线测量中，观测方向应归算到大地线方向。GeodesicDifferentialEquations大地线微分方程组大地线微分方程组由三个一阶常微分方程构成，分别描述了大地线上大地纬度B、大地经度L和大地方位角A随弧长S的变化率。纬度微分方程dB/dS=cosA/M沿大地线前进时，纬度变化率取决于方位角的余弦分量与子午圈曲率半径M的比值关系。Latitude·dB/dS经度微分方程dL/dS=sinA/(N·cosB)经度变化率取决于方位角的正弦分量，并受卯酉圈曲率半径N和纬度B的双重调制。Longitude·dL/dS方位角微分方程dA/dS=(sinA·tanB)/N赤道处（A=0时）方位角不变；非赤道处方位角持续偏转，变化率受纬度正切值影响。Azimuth·dA/dS方程组耦合特征NonlinearCoupling三个方程相互耦合、非线性，无法解析求解，实际应用中需采用数值积分或级数展开方法。Numerical·数值求解GEODETICEQUATION克莱劳方程（ClairautEquation）克莱劳方程N·sinA=C揭示了大地线上一个守恒量：卯酉圈曲率半径与方位角正弦的乘积为常数。这一守恒定律决定了大地线的整体走向形态，是推导大地线性质和检验解算正确性的核心工具。01方程表达：N·sinA=C（常数），其中N为卯酉圈曲率半径，A为大地线在该点的大地方位角，C为克莱劳常数N·sinA=C02物理意义解释：大地线从低纬向高纬延伸时N递减，为保持C不变sinA必须递增，大地线逐渐趋向东西方向sinA↑03顶点纬度确定：大地线达到最高纬度（顶点）时A=90°，此时C=N顶·sin90°=N顶，可据此求出大地线能到达的最高纬度A=90°04与球面类比：球面上克莱劳方程退化为r·sinA=常数（即大圆弧性质），椭球面克莱劳方程是其推广形式r·sinA=cGEODESY·COMPARISON大地线与相对法截线的比较大地线与相对法截线在短距离内差异微小，但在长距离测量中差异显著。大地线是唯一确定的最短路径，解决了正反法截线不重合的问题，是建立国家大地控制网的数学基础。国家大地控制网测量标志01正反法截线不重合：从A到B的正法截线与从B到A的反法截线是两条不同的曲线，造成方向定义的二义性02非最短路径：法截线长度一般大于大地线长度，在长距离情况下差值可达米级量级03三角形闭合差：由法截线构成的椭球面三角形三内角之和不等于180°+球面角超，存在额外闭合差相对法截线的局限04唯一性保证：两点间大地线唯一确定，不存在正反方向的歧义，适合作为统一的测量基准05最短路径最优：大地线长度是椭球面上两点间的最短距离，满足距离测量的最优性原则06数学性质完备：大地线满足克莱劳守恒方程，具有完整的微分几何理论支撑，便于精确计算大地线的优越性GEODETICLINEMORPHOLOGY大地线的空间形态特征大地线在椭球面上呈现赤道两侧对称的S形波动形态，具有明确的顶点和赤道交点。其形态完全由克莱劳常数和起始条件决定。01S形波动特征大地线在赤道两侧来回摆动，呈周期性波形。每次与赤道交点处方位角最小（A₀），到顶点处方位角最大（90°）。02顶点性质大地线最高纬度点处A=90°，由克莱劳方程可得cosB顶=C/a·W顶，顶点纬度决定了大地线的"振幅"。03赤道交点特性大地线穿越赤道时方位角A₀最小，此时sinA₀=C/N赤道=C/a，赤道交点间距约为半个周期。04与球面大圆弧对比球面大圆弧是平面曲线，而椭球面大地线是空间曲线——投影到赤道面为S形，这是椭球扁率造成的本质差异。CHAPTER03地面观测值归算至椭球面从物理地面到数学椭球面——方向改正、距离改正与垂线偏差GeodeticReduction地面观测值归算的必要性与总体流程大地测量以参考椭球面为数学基准面，但实际观测在地面进行，两者之间存在高程差和垂线偏差。归算的目的是将地面观测的方向和距离精确转换为椭球面上对应的量，确保观测值能够用于椭球面坐标计算和控制网平差。01基准面差异观测在地面（物理面）进行，计算在椭球面（数学面）进行，两者高差从数百米到数千米不等。这种差异导致观测值无法直接用于椭球面计算，必须通过归算消除系统误差。02方向归算路径地面观测方向→消除垂线偏差→归算到椭球面法截线方向→归算到大地线方向。垂线偏差改正将地面垂线方向转换为椭球面法线方向。03距离归算路径地面斜距→化算为水平距离→加入高程归化改正→得到椭球面上的大地线长度。高程归化改正将水平距离投影到参考椭球面。04精度要求国家一、二等控制网要求归算精度达到毫米级，必须考虑所有系统性改正项。忽略任何改正项都可能导致控制网整体精度下降。GEODESY·VERTICALDEFLECTION垂线偏差的定义与分量垂线偏差是地面点铅垂线与参考椭球法线之间的夹角，分解为子午分量ξ和卯酉分量η。它是连接天文坐标系与大地坐标系的桥梁，也是将地面天文方位角归算为椭球面大地方位角的关键参数。GNSS接收机野外大地测量01定义：地面点实际铅垂线方向与参考椭球面法线方向之间的夹角，通常用θ表示02分量分解：子午圈分量ξ（北正南负）与卯酉圈分量η（东正西负），分别表征南北与东西方向的偏差03成因分析：地球内部密度不均匀与地形起伏共同引起，山区和地质构造复杂地区偏差显著增大04数值范围：平原2″–5″，山区10″–30″，青藏高原边缘可达60″+GEODETICREDUCTION地面方向归算至椭球面——三差改正三差改正是将地面观测方向精确归算到椭球面大地线方向的三项系统性改正：垂线偏差改正消除铅垂线与法线不一致的影响，标高差改正消除目标高程影响，截面差改正将法截线方向转换为大地线方向。垂线偏差改正（δ₁）消除测站点铅垂线与椭球法线不一致对水平方向的影响，与垂线偏差分量ξ、η和目标天顶距有关计算公式涉及ξ·sinA+η·cosA等项，一等三角测量中改正值可达数角秒，必须精确计算需要已知测站点的垂线偏差值，可通过天文测量或重力测量获得，也可由垂线偏差模型内插标高差改正（δ₂）与截面差改正（δ₃）标高差改正：消除照准目标高出椭球面引起的方向偏差，与目标高程和方位角有关，高差大时显著截面差改正：将正反法截线不重合的偏差改正到大地线方向，与距离平方成正比，短距离可忽略等级适用性：一等三项全做，二等可省δ₃，三四等通常仅做δ₁；GNSS时代三差改正概念已扩展GEODESY·DISTANCEREDUCTION地面距离归算至椭球面地面距离归算是将实测斜距或平距转换为椭球面上大地线长度的过程，核心是高程归化改正。归化改正量与测线平均高程成正比，在高海拔地区可达米级量级，是建立高精度控制网时必须严格处理的关键环节。斜距化平距利用两端点高差Δh将斜距S化算为水平距离D=√(S²−Δh²)，这是距离归算的第一步。D=√(S²−Δh²)高程归化改正将地面水平距离D归算到椭球面弧长S₀，改正量ΔS=−H·D/Rm，H为测线平均大地高。ΔS=−H·D/Rm改正量级估算平均高程1000m、距离10km时改正量约1.57m；平原地区（100m）约0.16m，山区影响显著。1.57m/10km精确公式应用长距离精密测量需使用包含e²项的精确归化公式，考虑子午圈和卯酉圈曲率半径的差异。e²精确归化GEODETICCORRECTION高程归化改正的精确计算高程归化改正的精确公式包含椭球扁率项和高程高阶项，在高海拔和长距离条件下与近似公式存在显著差异。掌握精确归化方法对于西部高海拔地区的大地测量尤为重要。不同高程下的距离归化因子（每公里改正量）平均大地高H归化因子(1-H/N)每公里改正量典型地区100m0.99998432-15.7mm/km东部平原500m0.99992160-78.4mm/km丘陵地区1000m0.99984321-156.8mm/km云贵高原2000m0.99968643-313.6mm/km黄土高原4000m0.99937291-627.1mm/km青藏高原5000m0.99921615-783.8mm/km珠峰大本营高程每增加1000m，每公里距离归化改正约增加157mm，青藏高原地区改正量可达6‰以上GEODETICFUNDAMENTALS天文大地归算与坐标联系通过比较天文坐标与大地坐标的差异确定垂线偏差，是连接物理观测与数学模型的桥梁，至今仍是检验坐标系统精度的重要手段。天文坐标与大地坐标天文经纬度(φ,λ)由铅垂线定义，大地经纬度(B,L)由椭球法线定义，两者差值即为垂线偏差。ξ=φ−B垂线偏差公式子午分量ξ=φ−B，卯酉分量η=(λ−L)·cosB，通过天文观测与大地坐标比较可直接求得。η=(λ−L)·cosB拉普拉斯方程天文方位角α与大地方位角A的关系为A=α−(λ−L)·sinB，用于将天文方位角归算为大地方位角。A=α−(λ−L)·sinBGNSS时代的演变现代GNSS可直接测定大地坐标，垂线偏差可通过重力场模型或数字高程模型间接计算。重力场模型·DEMCHAPTER04椭球面大地问题解算方法大地问题正算与反算——从理论公式到实用算法GEODESY·FUNDAMENTALS大地问题正算与反算的定义大地问题正算是由已知点坐标、大地线长度和方位角推算未知点坐标，反算是由两点坐标推算大地线长度和方位角。大地问题正算已知条件起点大地坐标(B₁,L₁)、大地线长度S、起点方位角A₁₂求解目标终点坐标(B₂,L₂)、反方位角A₂₁物理意义从已知点出发，沿给定方向走给定距离，确定到达点坐标——类似"导航"问题ForwardProblem大地问题反算已知条件起点大地坐标(B₁,L₁)和终点大地坐标(B₂,L₂)求解目标大地线长度S、起点方位角A₁₂和反方位角A₂₁物理意义已知两点地理位置，求它们之间的距离和方向——类似"定位"问题InverseProblemGEODETICPROBLEM大地问题解算的三类方法椭球面大地问题的解算方法分为三大类：勒让德级数展开法、辅助球面法和数值积分法。各类方法在计算精度、适用范围和实现复杂度上各有特点，实际应用中需根据测量等级和测线长度选择最优方案。01将ΔB、ΔL、ΔA展开为S的幂级数，各系数含B、L、A对S的各阶导数，由大地线微分方程求定02适用于中短距离（300km以内），级数截断到e²项可达毫米级精度，计算效率较高03长距离时级数收敛慢，需取更多项；系数推导复杂，不同公式形式繁多04贝塞尔方法引入辅助球面，建立椭球面与球面的映射关系，适用于任意长度大地线05数值积分法直接对微分方程进行龙格-库塔等数值积分，实现简单但计算量较大数学公式推导与大地测量解算方法教学场景GEODESY·METHOD贝塞尔大地问题解算方法贝塞尔方法通过引入辅助球面建立椭球面与球面的映射关系，将复杂的椭球面大地问题转化为简洁的球面三角问题。01辅助球面引入构造一个与椭球面在特定条件下对应的辅助球面，使椭球面大地线对应球面大圆弧02归化纬度替换用归化纬度u替代大地纬度B，简化椭球面到球面的映射关系，核心公式tanu=√(1−e²)·tanB03球面求解在辅助球面上利用球面三角公式（纳皮尔公式等）求解大地线元素，计算过程简洁优雅04反算回椭球面将球面上求得的结果通过反变换公式转换回椭球面上的大地坐标和方位角GEODESY·COMPUTATION贝塞尔正算的具体步骤贝塞尔正算通过五个步骤完成从椭球面到辅助球面再回到椭球面的完整计算：归化纬度替换、初始参数计算、球面弧长求解、球面三角计算和反变换回椭球面。每步都有明确的几何对应关系。01归化纬度替换—将起点大地纬度B₁化为归化纬度u₁，满足tanu₁=√(1−e²)·tanB₁，从而建立椭球面与辅助球面之间的纬度对应关系。tanu₁=√(1−e²)·tanB₁02初始参数计算—由大地方位角A₁₂和归化纬度u₁计算sinA₀=sinA₁₂·cosu₁（克莱劳常数的变形），同时求取辅助球面上的起始经度差。sinA₀=sinA₁₂·cosu₁03弧长换算—利用大地线长度S与辅助球面弧长σ的关系式σ=f(S)，通过迭代方法求得精确的球面弧长σ。σ=f(S)04球面三角求解—在辅助球面上用球面三角公式求解u₂和Δλ，再反变换为大地坐标B₂、L₂及反方位角A₂₁。u₂·Δλ→B₂,L₂,A₂₁05反变换回椭球面—将辅助球面上求得的归化纬度u₂反变换为大地纬度B₂，完成从辅助球面回到椭球面的坐标转换，满足tanB₂=tanu₂/√(1−e²)。tanB₂=tanu₂/√(1−e²)GeodeticComputation大地问题解算方法精度对比不同大地问题解算方法在精度、适用范围和计算复杂度上各有优劣。勒让德级数法适合中短距离高效计算，贝塞尔方法适合任意长度高精度解算，数值积分法灵活性最强。实际选择需综合考虑测量等级和测线长度。三种大地问题解算方法的综合对比解算方法适用距离典型精度主要特点勒让德级数法<300km1-10mm计算快速，系数推导复杂，长距离需多阶项贝塞尔方法任意长度<1mm理论严密，精度均匀，长距离首选方法数值积分法(RK4)任意长度取决于步长实现灵活，编程简单，计算量随精度增大Vincenty公式任意长度0.5mm迭代求解，程序化强，GPS领域广泛使用球面近似法<50km1-10m公式最简单，仅用于低精度或初步估算贝塞尔方法和Vincenty公式