Excel在工程测量中的应用

Excel是测量内业数据处理最基础、最灵活的工具,几乎每个测量人员都离不开它。相比专业测量软件,Excel的优势在于:门槛低、功能强大、灵活性高、可自定义程度高。本文系统介绍Excel在导线计算、水准平差、坐标换算三大核心模块中的应用,并附详细公式与操作思路。

一、导线计算

导线测量是控制测量的主要方法之一,Excel可以实现闭合导线和附和导线的全自动计算,包括角度闭合差检核、坐标增量闭合差调整、坐标推算。

1.1 闭合导线计算

已知条件:起始点A坐标(XA, YA)、起始边方位角α0、观测角β1n(左角)、边长D1~Dn

(一)角度闭合差计算

角度闭合差 fβ 按下式计算:

fβ = Σβ – (n – 2) × 180°

其中 n 为导线角个数。

允许角度闭合差(限差)按《工程测量标准》GB/T 50026-2017 规定:

测量等级 方位角闭合差 (″)
一级 ±16√n
二级 ±24√n
三级 ±40√n
图根级 ±60√n

(二)角度改正数计算

当 |fβ| ≤ fβ允 时,按下式计算各角改正数 vβ

vβ = -fβ / n(平均分配)

各改正后角度:β′i = βi + vβ

(三)推算方位角

起始方位角:α1 = α0 + β′1 ± 180°(左角加,右角减)

各边方位角递推:αi+1 = αi + β′i+1 ± 180°

注意:方位角范围应归算至 0°~360°,可用公式:

=MOD(计算值, 360)=IF(计算值<0, 计算值+360, IF(计算值>=360, 计算值-360, 计算值))

(四)坐标增量计算

ΔXi = Di × cos αi

ΔYi = Di × sin αi

Excel公式(假设方位角在C列、边长在D列):

  • ΔX:=D2*COS(C2*PI()/180)
  • ΔY:=D2*SIN(C2*PI()/180)

(五)坐标闭合差调整

坐标增量闭合差:fx = ΣΔX,fy = ΣΔY

导线全长闭合差:f = √(fx² + fy²)

全长相对闭合差:K = f / ΣD

允许相对闭合差(以分母K表示):

测量等级 相对闭合差
一级 1/15000
二级 1/10000
三级 1/5000
图根级 1/2000

坐标增量改正数(按边长比例分配):

vxi = -(fx / ΣD) × Di

vyi = -(fy / ΣD) × Di

Excel公式(设fx所在单元格为F_x,ΣD为S_D,当前边长为D2):

  • vx:=-($F_x/$S$D)*D2
  • vy:=-($F_y/$S$D)*D2

(六)坐标推算

调整后增量:ΔX′i = ΔXi + vxi

调整后坐标:

  • Xi+1 = Xi + ΔX′i
  • Yi+1 = Yi + ΔY′i

1.2 附和导线计算

附和导线两端有已知控制点,计算方法与闭合导线基本相同,区别在于角度闭合差的计算。

角度闭合差:fβ = Σβ – (α – α + n × 180°)

其中 α 为起始边方位角,α 为终止边方位角,n 为观测角个数。

附和导线的坐标增量闭合差计算以附和路线两端的已知坐标为基准:

fx = ΣΔX – (X终_已知 – X始_已知)

fy = ΣΔY – (Y终_已知 – Y始_已知)

其余步骤与闭合导线完全相同。

1.3 Excel导线计算表结构建议

内容 关键公式
A 点号 直接填写
B 观测角 β 直接填写
C 方位角 α =MOD(前角+β±180, 360)
D 边长 D(m) 直接填写
E ΔX =D×COS(C×π/180)
F ΔY =D×SIN(C×π/180)
G 坐标 X 递推公式
H 坐标 Y 递推公式

二、水准平差

水准测量用于测定地面点高程,Excel可实现二、三、四等水准及图根水准的平差计算。

2.1 闭合水准路线平差

已知条件:起始点A高程 HA、观测高差 h1~hn、路线长度 L1~Ln(km)

(一)闭合差计算

闭合差:fh = Σh(往返测高差代数和)

允许闭合差(按测段距离):

  • 二等:±4√L mm
  • 三等:±12√L mm
  • 四等:±20√L mm
  • 图根级:±40√L mm

其中 L 为路线总长度(km)。

(二)改正数计算

各测段改正数(按测段长度比例分配):

vi = -(fh / ΣL) × Li

各点高程平差值:Hi = HA + Σ(h + v)

2.2 附和水准路线平差

已知条件:起点高程 H、终点高程 H、各测段高差 hi、测段长度 Li

(一)闭合差计算

fh = Σh – (H – H)

允许闭合差:同闭合水准。

(二)改正数计算

各测段改正数:

vi = -(fh / ΣL) × Li

各测段平差后高差:h′i = hi + vi

(三)推算待定点高程

H1 = H + h′1

H2 = H1 + h′2

……

2.3 Excel公式示例

设观测高差在B列,测段距离在C列(km),路线总长在$C$10,闭合差在D1:

  • 路线总长:=SUM(C2:C9)
  • 闭合差:=SUM(B2:B9)(闭合路线)
  • 改正数(每行):=-($D$1/$C$10)*C2
  • 平差高差:=B2+D2

三、坐标换算

坐标换算是测量内业的核心操作,包括高斯投影正反算、坐标正算、坐标反算。

3.1 坐标正算(极坐标法)

已知点A坐标(XA, YA)、边长D、方位角α,求点B坐标(XB, YB)。

XB = XA + D × cos α

YB = YA + D × sin α

Excel公式(XA在A2,YA在B2,D在C2,方位角(°)在D2):

  • XB=A2+C2*COS(D2*PI()/180)
  • YB=B2+C2*SIN(D2*PI()/180)

注意:方位角须按”北偏东为正”换算,不同象限需用 ATAN2 函数自动处理象限:

=ATAN2(ΔY, ΔX) × 180/π

3.2 坐标反算

已知两点坐标(XA, YA)、(XB, YB),反算边长D和方位角α。

边长:D = √[(XB-XA)² + (YB-YA)²]

方位角:α = ATAN2(YB-YA, XB-XA)

Excel公式(XB在A3,XA在A2,YB在B3,YA在B2):

  • 边长:=SQRT((A3-A2)^2+(B3-B2)^2)
  • 方位角(°):=IF(A3>=A2, ATAN2(B3-B2,A3-A2)*180/PI, ATAN2(B3-B2,A3-A2)*180/PI+180)

方位角结果需归算到 0°~360° 范围:

=MOD(结果+360, 360)

3.3 经纬度与平面坐标换算(高斯投影)

我国使用高斯-克吕格投影,3°分带或6°分带。以下公式适用于通用椭球参数(可按项目所在坐标系修正)。

(一)正算:经纬度 (L, B) → 平面坐标 (x, y)

设中央子午线经度为 L0(取整),经差 l = L – L0(秒),求平面坐标:

x = X + l²×N×sinB×cosB / (2×ρ′²) + l⁴×N×sinB×cosB³ / (24×ρ′⁴) × (5 – t² + 9η² + 4η⁴) + …

y = l×N×cosB / ρ′ + l³×N×cosB³ / (6×ρ′³) × (1 – t² + η²) + l⁵×N×cosB⁵ / (120×ρ′⁵) × (5 – 18t² + t⁴) + …

其中:

  • B 为纬度(°),l 为经差(″)
  • N = a / √(1 – e²sin²B)(卯酉圈曲率半径)
  • t = tan B,η² = e′²cos²B,e′² ≈ 0.0067395
  • ρ′ = 206265″(弧秒换算常数)

Excel简化计算(设中央子午线 L0 在 E2,L 在 D2,B 在 C2):

// 计算辅助参数
l = (D2-$E$2)*3600                    // 经差(秒)
B_rad = C2*PI()/180                    // 纬度(弧度)
N = 6378137/SQRT(1-0.00669438*SIN(B_rad)^2)  // 卯酉圈曲率半径
t = TAN(B_rad)
eta2 = 0.0067395*COS(B_rad)^2

// x坐标(米)
x = X(纬度) + l^2*N*SIN(B_rad)*COS(B_rad)/(2*206265^2)

// y坐标(米)附加带号
y = l*N*COS(B_rad)/206265 + 带号(39/40等)

(二)反算:平面坐标 (x, y) → 经纬度 (L, B)

反算相对复杂,需要迭代。先计算辅助量,再求纬度 B,再求经差 l,进而求经度 L。

实际工作中建议:使用项目所在地坐标系统的已知控制点坐标进行参数校准,或使用标准化的高斯投影计算工具,避免手工迭代出错。

3.4 坐标换算综合示例

整理一份标准的坐标换算表,包含以下功能区块:

功能 输入 Excel核心函数
坐标正算 XA, YA, D, α COS/SIN + 加法
坐标反算 XA, YA, XB, YB SQRT + ATAN2
方位角归算 任意角度值 MOD + IF
高斯正算 L, B, L0 多项式展开
高斯反算 x, y, L0 迭代法

四、高级技巧与效率提升

4.1 数组公式应用

Excel数组公式可以一次性对整列数据进行计算,大幅提升效率。

计算导线全长闭合差(单行公式):

=SQRT(SUM(E:E)^2+SUM(F:F)^2)

批量计算坐标增量(数组公式):

选中E2:E10,输入公式 =D2:D10*COS(C2:C10*PI()/180),按 Ctrl+Shift+Enter 完成。

4.2 条件格式实现自动检核

利用条件格式,当测量数据超限时自动标红提示:

  1. 选中角度列 → 条件格式 → 新建规则 → 使用公式确定要设置格式的单元格
  2. 公式:=ABS(fβ)>容许值,设置红色填充
  3. 同样对坐标闭合差 fx、fy 设置条件格式

4.3 数据验证防止输入错误

在导线计算表中,对边长列设置数据验证:

  1. 选中边长列 → 数据 → 数据验证
  2. 允许:小数
  3. 数据:介于 0.001 ~ 999999
  4. 出错警告:停止,标题”边长输入错误”,内容”边长必须在合理范围内”

4.4 VLOOKUP 实现控制点坐标调用

在日常工作中,可以建立一个控制点坐标数据库工作表(Sheet2),用 VLOOKUP 自动引用已知点坐标:

=VLOOKUP(A2, 控制点库!$A:$C, 3, FALSE)

这样在输入测站点号后,已知坐标自动填充,减少人工输入错误。

4.5 常用辅助公式集

目的 公式 说明
度分秒转弧度 =度*PI()/180+分*PI()/10800+秒*PI()/648000 角度制转弧度
弧度转度分秒 =度+分/60+秒/3600 弧度制转角度
方位角象限处理 =MOD(ATAN2(dy,dx)*180/PI()+360,360) 0~360° 方位角
容许闭合差 =容差系数*SQRT(SUMPRODUCT(距离列)) 12√L 形式
多边形面积(坐标法) =SUMPRODUCT(坐标阵列)/2 用于土方量估算

五、总结

Excel在工程测量内业中的应用远不止以上三个方向,还可用于横断面超欠填计算、土方量方格网法、变形监测数据分析、测量仪器检校计算等。掌握好Excel的公式编写、数据验证、条件格式和数据透视表四大技能,可以将大部分测量内业工作从繁琐的手工计算中解放出来,大幅提高工作效率和数据准确性。

建议测量人员根据本单位常用的测量方法和记录格式,建立一套标准化的Excel测量计算模板库,长期积累复用,是提升测量内业效率的有效途径。

声明:土木笔记网提供的部分资料来自互联网下载,纯属交流学习,如有侵权,请联系删除,请网友下载后24小时内删除。

联系方式:发送邮件至460147900@qq.com

(0)
土木笔记的头像土木笔记
南方CASS绘制地形图实战
上一篇 2026年3月23日 上午4:28
测量程序计算器fx-5800P编程
下一篇 2026年3月23日 上午4:28

相关推荐

  • 南方CASS绘制地形图实战

    南方CASS是国产主流的测绘成图软件。 一、数据导入 支持全站仪数据导入、GPS数据导入、CAD格式导入。 二、编图操作 绘制地物符号、标注高程点、绘制等高线。 三、图幅整饰 添加图廓、坐标格网、图例、指北针。 四、质量检查 检查高程点注记、地物完整性、等高线合理性。

    2026年3月23日
    2400
  • 测量数据处理与软件应用:数字化测绘技能

    测量数据处理是测绘工作的重要环节,现代测量已全面进入数字化时代。 一、数据处理基础 平差计算 利用最小二乘法原理,对观测值进行平差,消除矛盾,提高精度。 二、常用测量软件 南方CASS:地形图绘制 AutoCAD:工程制图 纬地/海地:道路工程设计 PKPM/迈达斯:结构计算 ArcGIS:GIS分析 三、Excel在测量中的应用 导线计算表 水准测量记录表…

    2026年3月23日
    3300
  • 测量程序计算器fx-5800P编程

    卡西欧fx-5800P是测量人员常用的程序计算器。 一、坐标正算 已知起点坐标、距离、方位角,计算终点坐标。 二、坐标反算 已知两点坐标,计算距离和方位角。 三、路线曲线计算 计算缓和曲线、圆曲线要素和主点坐标。 四、程序分享 常用测量程序代码示例。

    2026年3月23日
    2700

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

扫码关注
扫码关注
邮箱

如需帮助请发邮件。

邮箱账号:460147900@qq.com

分享本页
返回顶部
Index