![[分享]基于近景摄影测量的基坑土方量计算及精度评价](/skin/img/logo.png){kind=logo}
[分享]基于近景摄影测量的基坑土方量计算及精度评价
摘 要: 土方量的计算是工程建设中的重要环节, 计算精度将对合理选择施工方案有很大影响。针对传统土方量测量方法外业任务繁重、 作业复杂、 不能快速计算等问题, 提出一种以非量测数码相机、 多基线数字近景摄影测量为基础, 通过获取D TM数据来计算土方量的方法。研究表明, 与传统测量方法相比, 计算结果误差小于5%, 满足土方量计算的精度要求。也验证了近景摄影测量用于土方量计算的实用性和可行性, 为快速、 精准计算土方量提供了可靠途径。
关键词: 近景摄影测量; 多基线; 非量测相机; 土方量计算
进行基坑开挖时, 如何较简便、 准确、 快捷计算土方量, 作为安排施工进度计划、 编制预算和进行工程验收、 结算的依据, 是施工中需要解决的一个技术问题。传统的土方量计算方法一般使用全站仪、 全球定位系统( GPS 、 激光雷达(LIDAR) 等设备获取地面点数据导入CASS计算。但这些方法存在作业复杂、 购置费用昂贵、 外业数据海量、 数据处理周期长等问题。因此, 需要根据具体的工程情况将合适的技术应用于土方量计算, 以此来提高计算效率。张新东[1]介绍了利用全站仪进行土地平整过程中的土方量计算的方案, 但该方案测量条件要求苛刻, 测区范围必须通视, 测量目标必须裸露, 需采集大量野外点位数据, 操作繁琐, 任务繁重; J i a -c h o n gD u,Hu n g c h a oT e n g[2]将 G P S和激光扫描技术结合, 用于计算山体滑坡的土方量大小。将两项技术结合之后, 可对不能人员进入的危险区域进行土方量计算。其中激光扫描技术弥补了 G P S不能及时快速获得坐标点云数据的不足; 孟志义将激光雷达技术应用于土方量计算[3], 实验结果表明, 计算速度快, 效率高, 且填挖比基本保持平衡; 张磊,王晏民等人将激光雷达技术用于复杂区域的土方量计算[4], 实验结果表明, 在处理复杂区域时, 相对于传统土方量计算方法, 该方法计算速度快, 精度高, 填挖方误差率保持在2 . 5%左右; F u k a iJ i a,J o n a t h a nL i等人将激光雷达系统应用于道路重建中的土方量计算[5], 与传统计算方法结果对比分析得出, 该系统计算速度得到较大提升, 且计算误差保持在4%以内; F a nQ i a n g,Y a n a nZ h a n g等人将激光雷达应用于露天矿区地表土开挖过程中的土方量计算, 设计了一种基于三维激光扫描技术的露天矿土方量计算程序, 该程序计算精度高, 并且能快速计算土方量[6]; S l a t t e r yK e r r yT, S l a t t e r yD i -a n n eK等人应用激光雷达测量技术设计了一种快速、 高效准确计算土方量的激光扫描道路施工计算模型[7]。但是以上使用激光雷达的方法, 存在使用仪器昂贵, 外业操作繁琐, 不适合进行大规模推广的缺点。目前近景摄影测量用于土方量的计算还没有见到报道, 本文提出一种将近景摄影测量技术应用于土方量计算的方法, 并验证了其可行性。
本文以多基线近景摄影测量理论为基础, 将非量测数码相机SWD C 5 7标定, 用数字近景摄影测量技术处理获取的影像数据, 得到点云数据, 并进一步计算出土方量。实验结果表明, 与传统测量方法相比, 该方法计算结果误差小于5%, 满足土方量计算的精度要求。
1 理论依据
1.1 非量测相机标定
非量测数码相机, 由于其镜头光学畸变和面阵内畸变以及结构的不稳定, 用于近景摄影测量将引入较大误差。而专业的光学量测相机内方位元素固定, 获取的影像质量好, 但价格昂贵, 摄影技术要求过高, 对外业拍摄姿态要求苛刻。而实地地形情况往往较复杂, 不利于量测型相机获取影像数据。因此, 本实验采用非量测数码相机标定的方法, 事先对相机进行稳固和标定。相机标定采用大型检校场的标定方法。该方法需拍摄2 4张影像, 采用9个参数, 标定结果可满足近景摄影测量精度需要。
标定的畸变差改正公式如下[8]:
式中: ( x, y) 为像点坐标; (x0, y0) 为像主点坐标;k1, k2, k3 为径向畸变参数; p1, p2 为偏心畸变参数;b1, b1 为 面 阵 内 畸 变 参 数; r 为 像 点 径 向,
利用标定参数, 对拍摄的检校影像进行重采样, 对采样过后的影像再次标定, 计算畸变量, 其最大值小于0.5μ m, 说明上述标定方法精度高, 可靠性强。
1.2 多基线摄影测量
传统的摄影测量在进行同名点匹配时一般对两张相片进行匹配, 选择不同的相片对目标点进行匹配可能获得多解。若采用多基线影像摄影测量,对多张相片进行同名点匹配, 正确的光线只能交于同一点, 从而能更好地获得匹配点[9]。如图1( a) 所示, 传统摄影测量在对目标点a进行匹配时, 只有两条光线, 可能出现误匹配于( a) 点, 解得错误点。若采用多基线影像进行匹配, 如图1( b) 所示, 正确的光线多于两条, 只能交于同一点, 能够较好地得到正确的物方点。
多基线摄影测量也解决了传统近景摄影测量交会精度低和影像匹配难的矛盾问题。多基线摄影测量采用短基线获取大重叠度的序列影像, 相邻摄站之间的基线长度小于摄影距离的2 0%, 使相邻影像之间的交会角小于10° , 同时由于为多基线, 因此总体交会角较大, 能够确保交会精度。并且由于有多个观测值, 增加了多余观测, 进一步提高了交会精度及影像匹配的可靠性。视场角α与交会角θ关系如图2所示[10-11]。
此外, 近景摄影测量相对于航空摄影测量, 由于摄影高度低, 摄影倾角大, 即使采用多基线立体匹配技术, 影像重叠9 0%也无法使匹配结果满足摄影测量要求, 为减小摄影倾角可抬高相机高度[8]。
1.3 多基线摄影测量的平差
采用多基线摄影方式在确定同名点时, 由于投影光线多于两条, 存在多余的投影光线, 产生多余观测; 同时由于匹配误差的存在, 出现观测值误差所以必须进行平差, 才能求出最佳同名点坐标。由于观测值是同名的像方坐标, 待求的是其物方坐标, 所以本实验采用间接平差模型。
1.4 基于点云的D TM 构建
DTM 数字地面模型,地形表面形态属性信息的数字表达, 是带有空间位置特征和地形属性特征的数字描述, 是土方量计算的可靠方法。建立DTM 可以采用多基线近景摄影测量生成的目标区大量的点云数据, 其中, 点的位置和密度都会影响 DTM 的精度, 进一步插值算法的选择也会影响其精度[12]。CA S S软件中的构网方式有两种: 规则格网和不规则格网。规则格网是将原始数据进行内插, 算出规则形状格网结点的坐标。为了避免内插方格网而牺牲的原始数据精度,不规则格网是经常使用的构网方式。不规则格网将原始坐标位置作为格网结点, 实际应用中主要采用的是不规则三角形格网( T r i a n g l eI r r e g u l a rN e t -w o r k, T I N) 。T I N 每个基本单元的核心是组成不 万方数据规则三角形的3个顶点的三维坐标, 直接利用原始数据, 构造出由邻接三角形组成的格网结构来建立DTM。
2 实验结果与分析
2 . 1 实验与精度分析
河南理工大学双子湖区域地形起伏较大, 四周较高, 中央低为水域, 使用传统测量方法获取高程数据费时费力, 为此采用数字近景摄影测量技术进行实验区数据采集。测区范围约3900m2; 为减少摄影倾角, 抬高相机高度, 将相机安置在高度约为15m的支架上, 距离测区边缘约2 0m, 摄影距离平均为70m。摄影采用的SWDC57相机, 如图3所示, 由哈苏 H3D数码相机经过稳固改造而成, 固定焦距5 0mm, 视场角52 ° , 影像重叠度约为8 0%~9 0%之间, 为提高数据处理速度, 影像重叠度取8 0%, 摄影基线为 8 m。CCD 尺寸: 36.7mm×49mm, 像元宽度6.8μm。采用等倾平行多基线正直摄影方式, 共拍摄8张影像。
2. 1. 1 控制测量
测区内均匀布设3 0个控制点, 控制点的测量,平面坐标采用全站仪前方交会, 高程采用四等水准。对控制点测量结果进行精度评定, 得到点位中误差Mp=5.22mm, 高程中误差 MH=4.3mm, 满足实验要求。
2. 1. 2 数据处理
使用Lensphoto软件对影像数据进行处理, 先进行自由网平差, 平差后的中误差为0.002185mm, 匹配精度达到1 / 3个像元, 能够满足软件限差要求[13],匹配点如图4所示。
进一步进行约束网平差, 30个点在X, Y, Z 方向上的中误差分别为±0.0031m, ±0.0028m,±0.0025m, 在各方向上基本达到1 / 2像元精度,满足基坑土方量计算精度要求。
2. 1. 3 土方量计算
Cass软件中, 土方量计算方法有 T I N法、 方格网法、 断面法及等高线法。T I N法适用于任何地形地貌, 计算精度最高; 方格网法适用于大面积土方量的计算, 尤其是地形起伏较小、 坡度变化较缓的区域; 断面法适用于线性地带; 等高线法在实际工程很少会用到, 一般都是用于估算。结合实验区为大面积片状区域的实际情况, 故分别对近景摄影测量和传统方法获得的点位数据使用 T I N 法和方格网法进行土方量计算, 并对相应结果进行对比分析, 在两种方法中, 都将湖面作为地面看待, 即不考虑湖面以下的填方量, 仅考虑从湖面填至湖边最高点所需的填方量。T I N法计算结果如图5所示, 其具体计算结果如表1所示。
使用方格网法计算土方量时, 为了验证近景摄影测量方法对于不同的方格宽度都具有较高的精度, 根据土方量计算精度要求, 选取方格宽度分别为0.2m, 0.5m和1m进行计算, 具体结果如表2所示。
为了验证近景摄影测量用于土方量计算的可行性, 需要检验其测量精度。以传统方法的土方量计算结果作为真值, 精度评定采用绝对误差和相对误差指标, 计算结果如表3所示。
绝对误差公式： Δ=x-L
式中: Δ为绝对误差; x为测量值; L为真值。
相对误差公式：δ=Δ / L×100%。
式中, δ为相对误差。
由表3可以得出如下结论:
1) 相对于基于传统测量的土方量计算结果, 近景摄影测量方法的精度可达到4.62%~4.68%, 说明该方法计算精度完全满足土方量计算的一般允许误差10%~20%[14]。
2) 相对于 T I N 法, 方格网法精度小于或等于T I N法的精度, 这主要是由于方格网法需要提取4个角点的高程坐标, 而在生成点云数据时, 所获得的点位坐标分布不均匀, 所得的点位不一定会落在角点位置, 需要内插计算, 对于地形起伏较大区域计算精度变低。
3) 针对方格网法, 方格宽度越大, 其结果精度越低。若采用方格网法计算土方量可通过减小方格网宽度来提高精度, 但是方格宽度越低, 工作量会明显增加。故应根据具体施工需要选择合适的方格网宽度。一般在计算土方量时, 当施工面积不超过5000m2, 方格网宽度不超过5m即可满足误差要求[15]。
针对本次实验, 分析误差来源主要有以下几个方面:
1) 实验区周围植被繁多, 导致摄影测量采集数据匹配困难, 使用Lensphoto多基线近景摄影测量系统进行点位坐标解算时, 坐标计算的中误差较大。
2) 摄影距离影响影像分辨率, 因此, 对点位坐标解算精度有较大影响。在满足被摄目标占据图幅2 / 3以上的条件下, 摄影距离越小, 点位坐标解算精度越高[16]。
3) 点云分布不均匀, 分布稀疏区域精度低。
2 . 2 作业流程
利用数字近景摄影测量技术进行土方量计算时, 在确定测区范围之后, 选择合适的摄影距离和摄影倾角, 计算出相机高度并确定摄影基线。利用标定过的数码相机获取已布好控制网的实验区数字立体像对; 然后将采集到的数字影像利用多基线近景摄影测量进行处理, 获得实验区地表的离散点坐标; 继而生成 DTM 计算出实验区填挖方量。整体工作流程如图6所示。
3 结 论
本文介绍了以非量测数码相机的可靠标定为基础, 利用标定后的相机和多基线近景摄影测量技术进行近景摄影测量, 进而获取基坑的密集点云数据, 并据此生成DTM 来计算土方量的原理、 作业流程。通过与传统方法计算结果对比分析, 精度可达到4.62%~4.68%, 证明了利用近景摄影测量方法进行土方量计算数据采集的可行性。近景摄影测量技术具有采集速度快, 费用低, 劳动强度低的特点, 适用于作业员不易到达的危险区域的测量, 将其用于工程土方量的计算具有其独特的优势。
参考文献:
[ 1] 张新东. 土地平整测量计算方法的探讨[ J ] . 测绘通报,2 0 0 2( 增1) : 8 8.
[ 2] DUJC, T ENGHC. 3 Dl a s e rs c a n n i n ga n dG P St e c h -n o l o g yf o rl a n d s l i d ee a r t h w o r kv o l u m ee s t i m a t i o n[ J] .A u t o m a t i o ni nC o n s t r u c t i o n, 2 0 0 7, 1 6( 5) : 6 5 7 - 6 6 3.
[ 3] 孟志义. 激光扫描技术在土方量计算中的应用及精度分析[ J ] . 北京测绘, 2 0 1 2( 4) : 6 4 - 6 6.
[ 4] 张磊, 王晏民, 王国利. 基于地面激光雷达技术计算复杂区域土方量[ J ] . 测绘通报, 2 0 1 4( S 2) : 1 5 5 - 1 5 8.
[ 5] FUKA IJ , J ONATHAN L,CHENG Wa n g, e ta l .E a r t h w o r kV o l u m e sE s t i m a t i o ni nA s p h a l tP a v e m e n t R e c o n s t r u c t i o nU s i n gaMo b i l eL a s e rS c a n n i n gS y s t e -r m[ J ] . I E E EG e o s c i e n c ea n dR e m o t eS e n s i n gS y m p o s i -u m, 2 0 1 4: 9 2 2 - 9 2 5.
[ 6] F AN Q i a n g, ZHANG Y a n a n, L IH a i l o n g. T h eS t u d y o nE a r t h w o r kC a l c u l a t i o no fO p e n - P i tM i n eB a s e do n T h r e e - D i m e n s i o n a lL a s e rS c a n n i n gT e c h n o l o g y[ J ] . I m -a g ea n dS i g n a lP r o c e s s i n g, 2 0 0 9: 1 - 4.
[ 7] S LATT E RY K T, S LATT E RYDK, P e t e r s o nJ a m e s P. R o a dC o n s t r u c t i o nE a r t h w o r k V o l u m eC a l c u l a t i o n U s i n gT h r e e - D i m e n s i o n a lL a s e rS c a n n i n g[ J] . J o u r n a l o fS u r v e y i n gE n g i n e e r i n g - A S C E, 2 0 1 2: 9 6 - 9 9.
[ 8] 李天子, 邹友峰, 刘昌华. 多基线数字近景摄影平面地表测量可行性研究[ J ] . 测绘科学, 2 0 1 3, 3 8( 5) : 2 5 - 2 7.
[ 9] 张剑清, 胡安文. 多基线摄影测量前方交会方法及精度分析[ J ] . 武汉大学学报( 信息科学版) , 2 0 0 7, 3 2( 1 0) :8 4 7 - 8 5 1.
[ 1 0]张祖勋, 杨生春, 张剑清, 等. 多基线 - 数字近景摄影测量[ J ] . 地理空间信息, 2 0 0 7, 5( 1) : 1 - 4.
[ 1 1]柯涛, 张祖勋, 张剑清. 旋转多基线数字近景摄影测量[ J ] . 武汉大学学报( 信息科学版) , 2 0 0 9, 3 4( 1) : 4 4 - 4 7.
[ 1 2]刘先林, 甘济五, 王铁山. CA S S软件中 DTM 的建立[ J ] . 测绘与空间地理信息, 2 0 1 0, 3 3( 4) : 2 1 8 - 2 2 0.
[ 1 3]傅长刚. 民用数码相机可量测化的探讨[ J ] . 地理空间信息, 2 0 0 9, 7( 6) : 1 4 7 - 1 4 9.
[ 1 4]何晗芝. 工程土方量计算方法优化分析及其程序实现的研究[ D] . 长沙: 湖南大学, 2 0 0 9.
[ 1 5]王铁生, 程鹏里, 赵东保, 等. 方格网法土方量计算及误差影响[ J ] . 测绘通报, 2 0 1 2( S 1) : 1 0 9 - 1 1 1.
[ 1 6]李春梅, 张连蓬. 采石场边坡多基线数字近景摄影测量精度的影响因素分析[ J ] . 测绘通报, 2 0 1 1( 1 0) : 2 3 - 2 6.
来源于：勘测联合网
![[分享]美国visual slope隧道设计功能强大的专业软件](/e/data/images/no.jpg)
