简介这是中国矿业大学《大地测量学》课程的期末知识点总结面向测绘工程专业本科生、考研复习者以及需要夯实测量基础的从业者帮助快速抓住高频考点、理清概念关系。资源为一份PDF格式文档体积仅356KB整体内容浓缩方便在移动设备上随时翻阅和标注。目前已有365人学习浏览适合考前冲刺、集中记忆和查漏补缺。文档按单元逐条归纳了坐标系统、正常椭球与椭球面、垂线偏差与水准面、高程系统、国家平面控制网与水准网、GPS与工程控制网、椭球定位等核心专题同时具体整理了天球直角坐标系、天文与大地坐标、正高/正常高/大地高、似大地水准面、高程异常、克莱劳方程、垂线偏差改正、标高等差改正等关键概念与公式有助于系统掌握大地测量学的知识框架与重难点。1. 一张期末知识点总结恰好是控制网内业的完整主线外业测完回来内业计算时最容易翻车的不是平差而是坐标系选错、投影带搞混、高程基准对不上。点位偏离几十厘米、闭合差莫名其妙超限往往都能追溯到这几个基本概念。中国矿业大学这份《大地测量学》期末知识点总结表面上是备考提纲实际上把控制测量内业的主线全串起来了从正常椭球、垂线偏差到高斯投影换带再到抵偿高程面的选择都是工程中每天要面对的问题。本文按坐标系与高程基准 → 椭球几何与大地线 → 高斯投影与换带 → 控制网与水准概算 → 坐标转换与局部坐标系这条线展开把其中可直接用于作业的技术点逐一展开并补上可运行的代码和参数说明。2. 坐标系与高程基准外业数据归属的前提2.1 三套坐标系的定义与定位要素资料里把坐标系分成天球直角坐标系、天文坐标系、大地坐标系和空间大地直角坐标系四类。天球直角坐标系的原点在地心Z 轴与地球自转轴重合XY 平面与赤道面重合X 轴指向春分点 γ。天文坐标系用天文经度 λ 和天文纬度 φ 描述点位基准线是铅垂线基准面是大地水准面。大地坐标系用大地经度 L 和大地纬度 B 描述点位基准线是椭球面法线基准面是参考椭球面所以又叫参心坐标系。三者的关系可以这样理解天文坐标是重力方向的测量结果大地坐标是几何模型上的结果。同一点上两者不一致根本原因就是垂线偏差。外业用全站仪测出来的方向和高程本质上基于铅垂线而内业平差和成果提交需要的是基于法线的大地坐标。这个转换不做后续计算全部建立在错误基准上。2.2 垂线偏差与拉普拉斯方程垂线偏差 μ 是同一测站点上铅垂线与椭球面法线的夹角通常分解为子午圈分量 ξ 和卯酉圈分量 ηξ φ - B南北方向分量η (λ - L)cosφ东西方向分量这里 φ、λ 是天文纬度、天文经度B、L 是大地纬度、大地经度。拉普拉斯方程把天文方位角 α 与大地方位角 A 联系起来A α - (λ - L)sinφ这个式子在导线网起算方位角计算中非常关键。实测天文方位角后必须用拉普拉斯方程改化为大地方位角才能在椭球面上参与解算。很多初学者直接把天文方位角当大地方位角用在纬度 40° 地区经度差 1″ 就会带来约 0.64″ 的方位角误差累积到边长 5km 的导线端点上就是 1.5cm 的横向位移。2.3 正高、正常高、大地高三种高程系统的取舍高程系统是外业水准测量成果归属的核心。资料中整理了三个基准面高程系统基准面基准线与大地高的关系正高 H正大地水准面铅垂线H大 H正 N正常高 H正常似大地水准面铅垂线H大 H正常 ζ大地高 H大参考椭球面椭球面法线—其中 N 为大地水准面差距ζ 为高程异常。似大地水准面是根据地面各点正常高沿垂线向下截取得到的连续曲面它与椭球面的差值就是高程异常。GNSS 测出来的是大地高而我国高程成果统一使用正常高系统所以实际作业中必须知道测区高程异常 ζ 才能完成转换。在缺乏精化大地水准面模型的地区直接用大地高代替正常高误差可能达数十米这也是低精度 RTK 测量中高程误差远大于平面误差的主要原因。2.4 从 54 北京到 CGCS2000椭球换代与参心到地心资料列出了我国使用过的几套坐标系参数。54 北京坐标系采用克拉索夫斯基椭球是参心坐标系80 国家大地坐标系采用 IAG-75 椭球大地原点在陕西省泾阳县永乐镇CGCS2000 则采用总地球椭球属于地心坐标系。GPS 使用的 WGS-84 椭球与 CGCS2000 椭球参数非常接近长半轴都是 6378137m扁率略有差异。这里有几个关键参数值得对照椭球长半轴 a (m)扁率 f 分母克拉索夫斯基 (1940)6378245298.3IAG-756378140298.257WGS-846378137298.257223563CGCS20006378137298.257222101实际作业中常见的一个坑手里拿的是 CGCS2000 成果但坐标转换软件里默认参数还是 WGS-84二者在平面上差异很小但在高程异常计算和长距离控制网中会产生厘米级到分米级的偏差。进行跨年代成果对比时务必先确认椭球参数与坐标系定义。3. 椭球几何参数与大地线理解投影变形的基础3.1 卯酉圈、子午圈与平均曲率半径椭球面上不同方向的法截线曲率半径不同这是高斯投影产生长度变形的几何根源。资料明确给出三个关键曲率半径卯酉圈曲率半径 N a / W其中 W √(1 - e²sin²B)W 称为第一基本纬度函数子午圈曲率半径 M a(1 - e²) / W³平均曲率半径 R √(MN) a√(1 - e²) / W²三者之间恒有 N R M只有在极点处 N M R c极点处的子午圈曲率半径。工程中做距离改化时常以测区平均纬度计算 R作为高斯投影变形估算的基本参数。3.2 用 Python 计算指定纬度的 M、N、R下面这段代码直接计算给定大地纬度下三个曲率半径便于在测区踏勘时快速评估投影变形量级。import math def curvature_radii(B_deg, a6378137.0, f1/298.257223563): 计算椭球面上某纬度的子午圈、卯酉圈和平均曲率半径 参数: B_deg: 大地纬度, 单位度 a: 椭球长半轴, 默认 WGS-84/CGCS2000 f: 扁率, 默认 WGS-84/CGCS2000 返回: (M, N, R), 单位米 e2 2*f - f*f # 第一偏心率的平方 B math.radians(B_deg) W math.sqrt(1 - e2 * math.sin(B)**2) M a * (1 - e2) / W**3 # 子午圈曲率半径 N a / W # 卯酉圈曲率半径 R math.sqrt(M * N) # 平均曲率半径 return M, N, R M, N, R curvature_radii(34.5) print(fB34.5°: M{M:.3f}, N{N:.3f}, R{R:.3f})这段代码的关键在 W 的计算e2由扁率导出sin(B)随纬度变化直接影响 W 值。在低纬度地区 M 与 N 相差不大到北纬 45° 附近两者差距超过 20km所以不能用单一曲率半径代表整个测区。实际作业中距离改化公式里的 R 应取测区平均纬度处的平均曲率半径而不是简单取 6371km 的球近似。3.3 大地线的本质与克莱劳方程椭球面上两点间最短程曲线叫大地线。资料里有几个实用结论值得记住相对法截线之间的夹角 δ大地线与相对法截线的夹角约为 Δ/3两点间距离 600km 时大地线与法截线的长度差仅 0.007mm完全可以忽略。因此实际计算中可直接把法截线长度当大地线长度使用。克莱劳方程 r·sinA C 描述的是大地线上任意一点的平行圈半径 r 与该点方位角 A 的正弦乘积为常数。这个常数 C 称为大地线常数。它的实践价值在于已知大地线起点纬度 B1 和方位角 A1可算出 C 值进而判断大地线在椭球面上的走向和最高纬度。子午线和赤道都是特殊的大地线前者满足 A0° 或 180°后者满足 B0°。3.4 地面观测方向归算至椭球面的三项改正地面方向观测值要归算到椭球面上需要依次加入三项改正改正项产生原因主要影响因素量级垂线偏差改正 δ1测站铅垂线与法线不一致观测方向的天顶距一般小于 1″山地可达数秒标高差改正 δ2照准点高程导致法截线偏离照准点高程 H2与 H2 成正比高山地区显著截面差改正 δ3法截线方向与大地线方向不一致测站到照准点距离 S距离 50km 时约 0.001″通常可忽略三项改正中δ1 和 δ2 是必须加的δ3 仅在长边控制网中考虑。垂直角越大垂线偏差改正越显著所以精密三角高程测量在山区受垂线偏差影响很大这也是 GNSS 水准在山区替代传统三角高程的重要动因。4. 高斯投影与坐标换带内业计算的枢纽4.1 高斯投影的三个条件与带号计算高斯投影是横轴椭圆柱等角投影满足三个条件投影后角度不变形正形投影中央子午线投影为直线且长度不变除中央子午线外其他点长度比都大于 1且离中央子午线越远变形越大。中国采用 6° 带和 3° 带两种分带方式。6° 带从零子午线起每 6° 一带带号 N 与中央子午线经度 LN 的关系为 LN 6N - 33° 带从 1.5° 起每 3° 一带带号 n 与中央子午线经度 Ln 的关系为 Ln 3n。两种分带之间的对应关系为 n 2N - 1。例如某点经度 120°落在 6° 带第 21 带中央子午线 123°也落在 3° 带第 40 带中央子午线 120°。同一控制点在不同投影带中的坐标完全不同跨带平差或成果拼接时必须把这些坐标统一到同一投影带。4.2 用 pyproj 完成高斯投影坐标换带坐标换带的规范流程是先将已知投影带坐标反算为椭球面大地坐标 (B, L)再按新投影带中央子午线 L0 计算经差 l L - L0最后正算得到新带坐标。全程等角投影保证换带前后经纬度严格一致精度只取决于正反算公式的截断误差。实际工程中直接用成熟库完成即可下面是 pyproj 实现 CGCS2000 3° 带 40 带与 41 带之间换带的示例。from pyproj import CRS, Transformer # 定义 CGCS2000 / 3-degree Gauss-Kruger zone 40 和 zone 41 crs_40 CRS.from_proj4(projtmerc lat_00 lon_0120 k1 x_040500000 y_00 ellpsGRS80 unitsm no_defs) crs_41 CRS.from_proj4(projtmerc lat_00 lon_0123 k1 x_041500000 y_00 ellpsGRS80 unitsm no_defs) transformer Transformer.from_crs(crs_40, crs_41, always_xyTrue) # 40 带内某点坐标, 注意使用带号对应的加 500km 后的 Y 值 x1, y1 40500350.25, 3412345.67 x2, y2 transformer.transform(x1, y1) print(f40带坐标: ({x1:.3f}, {y1:.3f})) print(f41带坐标: ({x2:.3f}, {y2:.3f}))代码中projtmerc指定横轴墨卡托投影lon_0是中央子午线k1表示中央子午线长度比取 1x_0是伪东加常数带号 40 对应 40500000带号 41 对应 41500000。always_xyTrue保证输入输出的坐标顺序为 (X, Y)。这里要注意如果数据源是标准 CGCS2000 3° 带坐标Y 值已经包含 500km 加常数伪东应写成 40500000 而不是 500000否则换带结果整体偏移 4000km。4.3 方向改正与距离改正椭球面上的方向观测值归算到高斯平面需要加方向改正椭球面上的大地线长度归算到高斯平面要加距离改正。子午线收敛角 γ 是方向改正的基础它满足三个规律γ 是经差 l 的奇函数在北半球点在中子午线以东 γ 为正以西为负纬度越高 γ 越大极点处达到最大。实用中常用近似公式 γ l·sinB (l³/3)·sinB·cos²B(1 3η²) 计算。距离改化公式可近似写为ΔS/S y_m² / (2R²)其中 y_m 为测线两端点高斯平面横坐标的平均值R 为测区平均曲率半径。这个式子直接说明了为什么离开中央子午线越远变形越大变形与 y_m 的平方成正比。例如 R 取 6371kmy_m 150km 时变形约为 1/3600即 1km 边产生约 28cm 的长度变化这种量级的改正完全不能忽略。5. 国家控制网与水准测量概算5.1 平面控制网的四种建立方法与布网原则资料整理出四种建立国家平面控制网的方法三角测量法利用通视条件构成三角形网精确测定一条边和一边一角后推算其他边角控制面积大且精度高精密导线测量把控制点连成折线已知一点坐标和一边方位角推算其余点灵活方便边长精度均匀三边测量精确测三边长度推算其他边角边角同测则边长角度全部实测精度最高但成本也最高。布网原则归纳为四条分级布网、逐级控制保持必要的精度应有一定的密度应有统一的规格。工程控制网同样遵循这些原则只是精度指标根据服务对象调整测图控制网按测图比例尺确定密度施工控制网侧重点位精度变形监测网则强调重复观测的稳定性。5.2 国家水准网的等级划分与精度指标国家水准网按精度分为四等一等水准网是国家高程控制网的骨干二等水准网是全面基础三、四等直接为地形测图和工程建设提供高程控制。精度评定有两个关键指标偶然中误差 MΔ ±√([ΔΔ/R] / (4n))其中 Δ 为测段往返测高差不符值R 为测段长度kmn 为测段数全中误差 MW ±√([WW/F] / N)其中 W 为经过正常位水准面不平行改正后的环线闭合差F 为环线周长kmN 为环数这两个指标是水准网质量评定的核心。偶然中误差反映的是观测本身的精度全中误差则反映了整个网的内部符合程度。实际平差中如果 MW 超限首先检查环线闭合差较大的区域很可能是标尺尺长误差或水准面不平行改正项遗漏。5.3 水准测量概算的完整流程水准外业结束后不能直接拿观测高差去平差必须经过概算。资料中列出了三个必要步骤水准标尺每米长度改正。当一对标尺每米长度的平均误差大于 ±0.02mm 时必须对观测高差施加改正正常位水准面不平行改正。三等以上水准测量按下式计算每个测段的改正数 εi它与测段平均纬度 φm 和测段始末点近似高程平均值 Hi 有关εi -Hi · α · Δφ其中 α 为正常重力垂直梯度相关常数水准路线闭合差计算与概略高程推求。先求各测段观测高差之和与不平行改正之和再计算环闭合差并分配概略高程的计算按 h_i h_i εi v_i 进行其中 v_i 是闭合差分配改正数。这里有个实操细节正常位水准面不平行改正与测段高差本身无关只与纬度和概略高程有关所以可以在水准网平差之前独立完成再作为固定改正加入观测值。6. 坐标转换与局部坐标系的建立6.1 空间直角坐标与大地坐标互转空间直角坐标 (X, Y, Z) 与大地坐标 (B, L, H) 之间的正算是直接的代入 N a / W按 X (N H)cosBcosLY (N H)cosBsinLZ (N(1 - e²) H)sinB 计算即可。反算需要迭代或用直接公式。下面给出一个不使用迭代的实用反算代码。import math def ecef_to_geodetic(X, Y, Z, a6378137.0, f1/298.257223563): 空间直角坐标转大地坐标 (WGS-84/CGCS2000) 参数: X, Y, Z: 地心直角坐标, 单位米 返回: B, L, H: 纬度(度), 经度(度), 大地高(米) e2 2*f - f*f b a * (1 - f) ep2 (a*a - b*b) / (b*b) # 第二偏心率平方 L math.atan2(Y, X) p math.hypot(X, Y) # 平面距离 theta math.atan2(Z * a, p * b) # 辅助角 B math.atan2(Z ep2 * b * math.sin(theta)**3, p - e2 * a * math.cos(theta)**3) N a / math.sqrt(1 - e2 * math.sin(B)**2) H p / math.cos(B) - N return math.degrees(B), math.degrees(L), H B, L, H ecef_to_geodetic(-2267743.23, 5009153.62, 3221334.23) print(fB{B:.9f}, L{L:.9f}, H{H:.3f})这个算法的关键是引入辅助角 θ 将大地纬度 B 从迭代循环中解耦出来避免了传统方法的反复迭代。ep2是第二偏心率平方在纬度较高时对精度影响明显。这段代码可直接用于 GNSS 成果从地心坐标到大地坐标的批量转换处理十万点级别数据耗时基本可忽略。6.2 七参数布尔莎模型及实现不同空间直角坐标系之间的转换采用七参数模型包含三个平移参数 (x0, y0, z0)、三个旋转参数 (εx, εy, εz) 和一个尺度比参数 m。转换关系为X₂ X₀ (1 m)R(ε)X₁其中 R(ε) 是三个坐标轴的微小旋转矩阵。布尔莎模型在旋转角很小时可线性化为资料中的矩阵形式。下面是核心实现import numpy as np def helmert7(X, Y, Z, dx, dy, dz, rx, ry, rz, scale): 布尔莎七参数坐标转换 rx, ry, rz: 旋转参数, 单位弧度微小角 scale: 尺度参数, 单位 ppm, 实际乘数 1 scale*1e-6 s 1.0 scale * 1e-6 # 旋转矩阵 (小角度近似) R np.array([ [1, -rz, ry], [rz, 1, -rx], [-ry, rx, 1 ] ]) src np.array([X, Y, Z]) dst np.array([dx, dy, dz]) s * R.dot(src) return dst # 示例: 从某参心系转到 CGCS2000, 参数来自已知重合点解算 x2, y2, z2 helmert7(1000.0, 2000.0, 3000.0, dx12.5, dy-8.3, dz6.7, rx0.0002, ry0.0001, rz-0.0003, scale0.8) print(f转换后坐标: {x2:.3f}, {y2:.3f}, {z2:.3f})代码中的旋转矩阵采用了小角度线性化形式适用于旋转角小于几十角秒的应用场景。实际操作中七参数应基于至少 3 个高精度重合点解算且重合点应均匀覆盖测区。解算残差能直观反映两套坐标系在该区域的吻合程度残差过大时说明存在局部变形此时七参数模型不适合应考虑格网改正模型。6.3 抵偿高程面的三种选择方法高程投影变形与高斯投影变形方向相反可以利用这一特性选择合适的高程面和中央子午线使两者抵消。资料给出的判断不等式为 H ≤ y² / (2R)其中 H 为测区平均高程y 为测区中心到中央子午线的距离。第一种方法选择抵偿高程面作为投影面。将原椭球面高程抬高至 H抵使高程归化改正与高斯投影距离改化大致抵消。实用做法是直接改装椭球长半轴令 a_new a H抵仍用 3° 带中央子午线。第二种方法保持椭球面不变选择任意投影带通过移动中央子午线减小 y 值。中央子午线移动量 l″ (H / (N·cosB)) · ρ″ 的平方根关系实际可按 l″ 7362″ · √(H / (N·cosB)) 估算。这种方法适合东西跨度较大但高程变化较小的测区。第三种方法同时改变投影面和高程面的中央子午线适用于既有较大高程又有较大偏离中央子午线的测区。本文还有配套的精品资源点击获取
