ch5地面观测元素归算至椭球面

地面观测元素归算至椭球面大地测量学基础·第五章边角观测元素归算理论与方法Contents课程目录大地测量归算方法系统课程，涵盖从基本概念到综合精度分析的完整知识体系。01归算的基本概念与意义02水平观测方向归算至椭球面03观测天顶距的归算处理04地面观测边长归算至椭球面05综合应用与精度分析Chapter01归算的基本概念与意义理解地面观测与椭球面计算之间的几何差异及其数学处理逻辑GEODETICREDUCTION归算的定义与本质归算是将基于铅垂线基准的地面实测元素，通过系统改正转化为以参考椭球面法线为基准的椭球面元素的过程。其本质是统一观测基准与计算基准，确保大地测量数据在数学曲面上的严密性和一致性。全站仪野外控制点观测场景01地面观测以铅垂线为基准全站仪、经纬仪等仪器通过水准管或电子补偿器按铅垂线方向整平，观测数据天然关联物理地球表面02椭球面计算以法线为基准参考椭球是规则数学曲面，所有坐标计算、大地线推算均沿法线方向进行，与铅垂线存在系统性偏差03归算的核心任务对地面观测元素施加适当改正数，消除铅垂线与法线不一致带来的系统误差，使之转化为椭球面上相应的严格数学元素04精度铁律归算过程不得损害野外观测的精度，改正数计算精度须高于观测精度至少一个数量级GEODETICREDUCTION归算的两种经典方法归算存在两条经典路径：皮齐蒂法通过大地水准面中转实现"两步投影"，海姆特法借助椭球法线直接"一步到位"。两者在理论上等价，但因20世纪中叶前大地水准面参数难以精确获取，实际应用长期受限。皮齐蒂法（两步投影）01借助垂线方向，将地面观测数据投影到大地水准面上，消除测站高程对观测的影响02借助椭球法线，将大地水准面上的数据进一步投影到参考椭球面上，完成从物理面到数学面的转换03物理意义清晰，每一步都有明确的地球物理对应面，适合教学和理论分析海姆特法（一步投影）01直接借助椭球法线方向，将地面观测数据一步投影到参考椭球面上，无需大地水准面作为中间过渡面02计算过程更为简洁，减少了中间环节的误差累积，在现代高精度大地测量中更具实用优势03两种方法在理论上严格等价，实际差异取决于大地水准面差距N和垂线偏差分量的确定精度GEODETICREDUCTION归算的总体流程与精度要求归算遵循"方向值→法截线→大地线"三步递进的严密流程，各等三角测量的归算精度要求依次递减：一等0.001″、二等0.01″、三四等0.1″，确保归算过程不损害原始观测精度。STEP01垂线偏差改正将地面观测方向值化为以椭球法线为准的方向值，需施加垂线偏差改正δ₁，消除测站铅垂线与法线不一致的影响δ₁STEP02标高差改正将法线方向值化为椭球面上的法截线方向值，需施加标高差改正δ₂，消除照准点觇标高引起的投影偏差δ₂STEP03截面差改正将法截线方向值化为椭球面上的大地线方向值，需施加截面差改正δ₃，消除法截线与大地线之间的差异δ₃PRECISION精度等级要求一等三角测量归算至0.001″，二等归算至0.01″，三四等归算至0.1″，改正数计算精度须匹配相应等级0.001″GeodeticReference参考椭球的选择与几何参数归算计算依赖参考椭球的几何参数（长半轴、扁率、曲率半径等），不同坐标系对应不同椭球。我国现行CGCS2000坐标系采用的椭球长半轴为6378137m，归算前必须明确椭球基准以确保计算一致性。常用参考椭球几何参数对比椭球名称对应坐标系长半轴a(m)扁率1/f克拉索夫斯基椭球1954北京坐标系6,378,245298.3IUGG-75椭球1980西安坐标系6,378,140298.257WGS-84椭球GPS/全球定位6,378,137298.257223563CGCS2000椭球国家大地坐标系6,378,137298.257222101CGCS2000与WGS-84椭球长半轴相同但扁率微差，归算时须注意坐标系一致性Geodesy·GeometricFoundations归算中的核心几何概念铅垂线（物理基准）、法线（数学基准）和垂线偏差（二者差异的度量）构成归算的几何基础，是连接物理地球与数学椭球的关键桥梁。大地测量控制点标石—野外测量的物理基准标志01铅垂线：由测站重力矢量决定的方向线，仪器水准管整平后竖轴沿此方向，是野外观测的物理基准物理基准02法线：椭球面上过某点且垂直于椭球切平面的直线，所有大地坐标计算均沿法线进行，是数学基准数学基准03垂线偏差(θ)：铅垂线与法线的空间夹角，分解为子午圈分量ξ（南北）和卯酉圈分量η（东西）ξ·η04量级：平原2-5″，丘陵10-20″，山区极端可超过30″，不可忽略>30″CHAPTER02水平观测方向归算至椭球面垂线偏差改正、标高差改正与截面差改正的数学推导与物理意义GEODETICCORRECTION垂线偏差改正δ₁垂线偏差改正δ1消除测站铅垂线与法线不一致对水平方向的影响，将观测方向从铅垂线基准转换到椭球法线基准。01改正公式δ₁=−u·sin(A−θ)·cot(z₁)，分量形式为δ₁=−(ξ·sinA+η·cosA)·tanα₁，其中α₁为垂直角02物理含义测站铅垂线与法线的夹角导致竖轴倾斜，使水平度盘读数产生系统偏差，改正量与天顶距余切成正比03为零条件照准点位于天顶方向（z₁=0或α₁=0）时改正为零；照准方位角与垂线偏差方位角一致时改正亦为零04量级评估垂线偏差10″、天顶距45°条件下，改正量可达10″量级，一等测量中绝不可忽略GeodeticCorrection垂线偏差改正的影响因素与量级分析垂线偏差改正量受三大因素共同控制：垂线偏差分量(ξ,η)的大小、观测方位与垂线偏差方位的夹角关系、以及天顶距的大小。在山区高精度测量中，改正量可达10″以上，远超一等测量的0.001″精度要求，必须精确计算。FACTORS影响因素分析垂线偏差分量(ξ,η)20–30″由测站区域的地球内部质量分布决定，山区可达20–30″，平原通常2–5″方位角差(A−θ)正交最大观测方位角与垂线偏差方位角的差值决定改正的投影分量，正交时最大、平行时为零天顶距z₁cot(z₁)↑cot(z₁)随天顶距增大而增大，低仰角观测的改正量显著大于高仰角观测PRECISION量级与精度控制01一等三角测量0.001″几乎所有方向的垂线偏差改正都必须精确计算并施加02二等三角测量0.01″垂线偏差大于5″时必须施加改正，小于2″时部分方向可酌情忽略03三四等测量0.1″垂线偏差较小时可简化处理，但山区仍需完整计算GEODETICCORRECTION标高差改正δ2标高差改正δ2消除因照准点觇标高度引起的投影偏差——当照准点高于椭球面时，该偏差随照准点高程H2和觇标高度增大而增大。δ2=ρ″·(H2²·sin2A1)/(2·R·N·S)H2为照准点大地高，A1为方位角，S为边长，R、N为曲率半径几何原理：觇标顶部沿铅垂线投影到椭球面时，与B点沿法线投影的位置存在偏差，水平方向分量即为标高差与H2²成正比：高程翻倍时改正量变为四倍，山区测量中影响尤为显著4×为零条件：A1=0°或180°（沿子午线方向）时sin2A1=0；H2=0时改正亦为零大地测量觇标——照准点高程H2直接影响标高差改正量级GEODETICANALYSIS标高差改正的量级评估与计算要点标高差改正量与照准点大地高的平方成正比，对高程极为敏感：H₂=1000m时改正可达0.02″量级，H₂=100m时仅约0.0002″。01典型量级H₂=1000m、A₁=45°、S=10km时δ₂≈0.02″，一等测量必须计算；H₂=100m时δ₂≈0.0002″，接近可忽略阈值0.02″02方位角调制效应sin2A₁在A₁=45°或135°时取最大值1，在A₁=0°或90°时为零，不同观测方向的改正量差异显著sin2A₁03高程基准注意H₂应为照准点的大地高（沿法线到椭球面的距离），而非正常高或海拔高，高程异常差异可能引入数米级误差大地高≠正常高04曲率半径取值R取测站处的平均曲率半径，N取卯酉圈曲率半径，两者在低纬度地区差异较小但高精度计算中仍需区分R·NGEODETICCORRECTION截面差改正δ3截面差改正δ3消除椭球面上法截线与大地线之间的方向偏差。大地线是两点间最短路径（空间曲线），法截线是含法线竖直平面截椭球面的平面曲线，二者差异与边长平方成正比，长边三角测量中不可忽略。FORMULA改正公式δ₃=−(e′²·S²·cos²Bm·sin2A₁)/(12·N²·ρ″)，其中e′为第二偏心率，Bm为中点纬度，S为边长δ₃GEOMETRY几何含义过A点的法截线与过B点的法截线不共面（正反法截线不重合），大地线位于两条法截线之间，改正使方向值指向大地线SCALELAW与S²成正比边长翻倍则改正量变为四倍。S=10km时δ₃通常小于0.001″可忽略，S=50km时可达0.02″以上必须计算S²LATITUDE纬度因子cos²Bm使低纬度地区的截面差大于高纬度地区，赤道附近影响最大，两极趋近于零cos²BmGeodeticCorrections三差改正综合对比三差改正分别解决基准线转换（δ1）、照准点高程投影（δ2）和曲线类型转换（δ3）三个不同层面的问题。各项改正的主要关系量、为零条件和精度要求各不相同，一等测量三项均须精确施加，三四等测量可酌情简化。改正项改正目的主要关系量为零条件各等要求

δ1垂线偏差改正

铅垂线基准→法线基准ξ、η、天顶距z₁α₁=0或A=θ一二等必算，三四等视情

δ2标高差改正

觇标顶部→椭球面投影H₂、方位角A₁、SH₂=0或A₁=0°/180°一二等必算，三四等可忽略

δ3截面差改正

法截线→大地线方向S²、纬度Bₘ、A₁S=0或A₁=0°/90°一等必算，二等长边必算三差改正分别对应基准线、高程和曲线三个独立维度的系统偏差修正CALCULATIONEXAMPLE三差改正计算实例通过一等三角测量典型算例演示三差改正的完整计算流程：在垂线偏差ξ=+3.2″/η=-1.8″、H2=800m、S=28.5km的条件下，三项改正量分别达到秒级、亚秒级和亚秒级，合计改正量对一等测量精度至关重要。已知条件01测站A：B₁=35°00′，L₁=115°00′，H₁=500m，ξ=+3.2″，η=−1.8″02照准点B：H₂=800m，观测方位角A₁=62°15′30″，天顶距z₁=87°22′10″，边长S=28.5km一等三角测量·方向观测归算计算结果垂线偏差改正δ₁：−(3.2×sin62.26°+(−1.8)×cos62.26°)×tan2.63°≈−0.089″标高差改正δ₂：代入H₂=800m、A₁=62.26°，量级较小但一等测量仍须施加≈+0.013″截面差改正δ₃：代入S=28.5km、Bm=35°，边长较长使得该项不可忽略≈−0.007″总改正量：δ=δ₁+δ₂+δ₃，归算后方向值为观测方向值加总改正量≈−0.083″CHAPTER03观测天顶距的归算处理天顶距观测中垂线偏差影响与大气折光差修正的理论方法GEODESY·VERTICALDEFLECTION天顶距归算的基本原理天顶距观测以铅垂线反方向为天顶基准，但椭球面计算需以法线反方向为基准。垂线偏差导致两种天顶方向不一致，改正量Δz=ξ·cosA+η·sinA，等于垂线偏差在观测方位平面上的投影分量，与天顶距大小本身无关。天顶距定义差异野外观测天顶距zobs以铅垂线反方向起算，椭球面天顶距zell以法线反方向起算，二者差值即为垂线偏差改正。改正公式Δz=ξ·cosA+η·sinA，其中ξ为子午圈垂线偏差分量，η为卯酉圈分量，A为观测方位角。方位依赖性观测方位与子午圈一致（A=0°）时改正量等于ξ；与卯酉圈一致（A=90°）时等于η。一阶近似特征该公式为垂线偏差对天顶距影响的一阶近似，偏差小于30″时精度足够，更大值需考虑高阶项。GeodeticSurveying大气折光差对天顶距的影响与修正大气密度梯度使光线路径向地面弯曲，导致观测天顶距偏小，折光差改正Δk=k·S/(2R)·ρ″与边长成正比。精密测量中采用对向观测取平均值的方法有效削弱折光影响。01折光原理大气密度随高度递减，光线在密度梯度中向地面弯曲，使观测天顶距比几何天顶距偏小约k·S/(2R)弧度。这一物理效应是三角高程测量中最主要的系统误差来源。02折光系数k一般取值0.10–0.14，受地面覆盖、温度梯度和观测时段影响显著，日出日落前后k值波动可达0.05以上，需根据气象条件动态估计。03改正公式Δk=k·S/(2R)·ρ″，S=20km、k=0.12时改正量约39″，量级远超观测精度，必须修正。边长越长，折光影响越显著。04对向观测削弱法在A、B两点同时观测天顶距并取平均值，可有效消除折光系数系统偏差，残余误差降至单次观测的1/5以下，是精密三角高程测量的标准作业方法。ZENITHDISTANCEREDUCTION天顶距的完整归算流程天顶距归算遵循"垂线偏差改正→折光差改正→仪器高/目标高改正"的三步流程，将物理观测天顶距逐步转化为椭球面上严格的法线天顶距。1垂线偏差改正施加Δz=ξ·cosA+η·sinA，将铅垂线基准的天顶距zobs转换为法线基准的天顶距znzobs→zn2折光差改正施加Δk=k·S/(2R)·ρ″，消除大气折射导致的光线弯曲影响，得到几何天顶距zgeomzn→zgeom3仪器高/目标高改正根据仪器高i和目标高v，将测站中心的观测天顶距归算到地面标石中心的天顶距i,v校正精度控制一等三角高程测量天顶距归算精度须优于0.5″，折光系数k宜通过对向观测或实测气象数据精确确定≤0.5″CHAPTER04地面观测边长归算至椭球面从地面实测斜距到椭球面大地线长度的系统归化方法DistanceReduction边长归算的总体思路与步骤边长归算遵循"斜距→平距→椭球面距离"的递进流程：首先利用天顶距将实测斜距化为测区高程面上的平距，再通过高程归化改正将平距"压缩"到椭球面上，最终得到椭球面上的大地线长度。每步改正须保持精度不损失。01斜距化平距利用观测天顶距z（或垂直角α），将光电测距仪实测斜距D转化为水平距离d=D·sin(z)=D·cos(α)d=D·sin(z)02高程归化改正将高程面上的平距d归化到椭球面上，改正数Δd_H=−d·Hm/R，其中Hm为测线平均大地高，R为曲率半径Δd_H=−d·Hm/R03投影改正将椭球面距离投影到高斯平面，改正数与离中央子午线的距离平方成正比，本章暂不涉及GaussProjection04精度铁律现代光电测距精度可达1/100万以上，归算各步骤的公式精度和参数精度必须匹配，否则归算成为精度瓶颈1/1,000,000+边长归算·第三项改正高程归化改正的数学推导高程归化改正将高程面上的平距"压缩"到椭球面上，改正数Δd=-d·Hm/R（负号表示椭球面距离更短）。该改正与平均大地高Hm成正比，Hm=500m、d=10km时改正量达-0.78m，量级远超测距精度，是边长归算中影响最大的单项改正。01推导思路高程面半径=R+Hm，椭球面半径=R，由相似比例关系得椭球面距离dell=d×R/(R+Hm)≈d−d·Hm/Rd·Hm/R02改正公式ΔdH=−d·Hm/R，其中Hm=(H₁+H₂)/2为测线两端大地高的平均值，R取测线中点处的平均曲率半径ΔdH03量级分析Hm=500m、d=10km、R≈6370km时，ΔdH≈−0.785m；Hm=100m时ΔdH≈−0.157m，量级始终远大于测距误差−0.785m04曲率半径取值严格计算中R应取测线中点处的平均曲率半径R=√(M·N)，M为子午圈曲率半径，N为卯酉圈曲率半径R=√(M·N)GEODETICREDUCTION精密边长归算的高阶改正与注意事项精密边长归算须在一阶改正基础上增加二阶项d·(Hm/R)²，长边或高海拔条件下该项可达毫米级。高阶改正项Δd₂二阶改正：Δd₂=d·(Hm/R)²，Hm=1000m、d=30km时约+0.7mm，一等测量长边条件下需纳入计算d_ell完整归化公式：d_ell=d×[1-Hm/R+(Hm/R)²-...]，实际应用中取到二阶项即可满足所有精度等级要求实际作业注意事项气象改正优先：光电测距的波速受温度、气压影响，归算前必须先施加气象改正消除波速偏差高程基准统一：归化改正须用大地高H，若仅有正常高h需加高程异常ζ（H=h+ζ），忽略ζ可引入数米级高程误差GPS基线特殊性：GPS直接测定的基线向量已是椭球面弦长，仅需弦长改弧长改正，不做高程归化CALCULATIONEXAMPLE边长归算完整算例通过典型算例演示边长归算全流程：实测斜距15832.456m经天顶距化平距后为15832.142m，再经高程归化改正（Hm=785m）后椭球面距离为15830.195m，总改正量约-2.26m，充分说明高程归化对长边测量的关键影响。01斜距化平距d=D·sin(z)=15832.456×sin(89°12′35″)=15832.456×0.99980≈15832.142m，改正量约-0.314m-0.314m02平均大地高Hm=(H1+H2)/2=(650+920)/2=785m，取R=6371000mHm=785m03高程归化改正Δd_H=-d·Hm/R=-15832.142×785/6371000≈-1.947m-1.947m04椭球面距离d_ell=15832.142-1.947=15830.195m，从斜距到椭球面距离的总改正量约-2.261m-2.261mCHAPTER05综合应用与精度分析不同测量场景下的归算策略选择与系统精度评估方法ReductionStrategy不同测量等级下的归算策略归算策略须与测量等级匹配：一等测量三差改正全项精算、边长归算含二阶项；二等视垂线偏差和边长决定简化范围；三四等在平原短边条件下可忽略δ2和δ3。测量等级δ1垂线偏差δ2标高差δ3截面差边长归算一等（0.001″）全部精算全部精算全部精算含二阶项二等（0.01″）全部计算H2>200m时计算S>20km时计算一阶改正三等（0.1″）ξ,η>5″时计算通常可忽略通常可忽略一阶改正四等（0.1″）酌情简化可忽略可忽略一阶改正测量等级越低，归算简化空间越大，但山区和长边条件仍须严格计算GPS·归算理论GPS时代对传统归算方法的变革与传承GPS直接给出椭球面三维坐标，免去了方向三差改正的繁琐计算，但归算理论并未过时：大地高到正常高的转换、地面控制网加密中的传统观测归算、以及坐标参考框架的理解，都依赖扎实的归算理论基础。GPS带来的变革直接测定椭球面坐标GPS/GNSS直接输出(B,L,H)，无需方向三差改正，外业归算工作量大幅减少基线向量简化处理GPS基线向量为椭球面弦长，仅需弦弧改正，免去高程归化步骤CORS网络实时定位连续运行参考站网络实现实时高精度定位，归算过程由软件自动完成归算理论的持续价值高程转换的核心支撑GPS大地高→工程正常高需精确高程异常模型，本质是物理面到数学面的转换传统场景不可替代精密导线、三角高程测量、工程控制网加密仍依赖地面观测归算理论参考框架理解基础不同坐标系转换、椭球参数变化分析，均以归算理论为根基PRECISIONANALYSIS归算精度的系统分析归算精度须高于观测精度一至两个数量级。垂线偏差改正受ξ、η测定精度制约，标高差改正取决于大地高H2确定精度，截面差改正精度最高。δ1垂线偏差改正精度m(δ1)≈√[m²(ξ)·sin²A+m²(η)·cos²A]·tanα；当m(ξ)=m(η)=0.5″、α=5°时，m(δ1)≈0.04″，满足一等精度要求。该改正项受垂线偏差分量测定精度直接影响。0.04″δ2标高差改正精度m(δ2)/δ2≈2·m(H2)/H2；当m(H2)=0.5m、H2=800m时，相对精度约0.125%，改正量本身的误差极小。标高差改正精度取决于大地高确