 点击上方“溪流之海洋人生”即可订阅哦我国地域幅员辽阔，地形地貌多样，近年来我国正大力开展海洋战略，对海湾地区的开发需求日趋强烈。但海湾地区地势复杂，周围分布的泥潭、水池等对野外测绘作业带来了诸多不便，为了获取海湾地区的地形数据，传统的野外RTK(Real-timeKinematic)测量不仅需要耗费大量的人力物力，还会因为交通不便等原因增加工作难度，效率较低，且获取地面的坐标数据更新周期长，而海湾地区受海水冲击影响，地形多变，数据精度难以保证。 
8 r3 }% T0 D3 ^* Z- |/ A 随着无人机技术的发展成熟，无人机摄影测量具有作业时间短、测量范围广、成果精度高等优势，一次数据采集可以得到影像、点云、DOM等多种数据，在许多地形测量项目中得以应用。为获取某一海湾地区的陆地及近岸海底地形数据，本研究结合无人机摄影测量与水深测量技术，完成了对海湾地区的地形测量方案，并通过实验验证了该方案的可行性和有效性。 0 _7 L5 N+ v! B* n
一、
! c+ z% ~ }- }* A1 z% k; M 技术原理 / z8 ^* Q8 L) s" E3 M+ V$ | Y
⒈测量基本原理 ! B5 \# M' V! Z0 x4 P& o! L
无人机摄影测量基本原理是通过无人机搭载的航摄相机连续曝光获取测区地面影像，通过内业空三解算得到加密点坐标，进而完成对测区的全部测量。通过搭载不同传感器可在短时内获取大面积范围的点云、影像、多光谱等多源数据，具有测量范围广、作业效率高、获取数据类型多样等特点。通过专业数据处理软件可最终获取测区的正射影像、DOM、DEM、点云数据、三维模型等结果，可满足不同行业应用的各种需求。本实验选用智航V200C垂直起降无人机，配套SONY黑卡相机进行航测实验。
: \" Q# E. O& T 水深测量基本原理是通过使用测深杆、水铊、回声测深仪、多波束回声测深系统和海底地貌探测仪等测得瞬时水面下的深度，经测深仪改正和水位改正，归算到由深度基准面起算的深度，最后获得水下待测点的平面和高程的坐标数据。本实验选用HY1200B型号声速剖面仪，HY1601型号测深仪进行单波束水深测量。
' B" N, K/ u$ y( N) \ ⒉数据融合原理
3 V% S" f/ \% R0 D- ^ 无人机摄影测量采集所获得的数据为影像数据，利用专业软件对数据完成空三解算后可生成DSM(DigitalSurfaceModel)，进而用于地形测量，在与水深数据融合前，需要将2组数据的坐标系统进行统一，即进行测量基准的转换。测量基准的具体要求为转换区域椭球面经纬度在3°以内，转换所用的基准控制点周边应有不少于10个稳定的IGS站、国内IGS站及国家级卫星导航定位基准站，并遵循以下原则：①连续性：测站连续观测3a(或以上)；②稳定性：站点坐标变化很小，具有已知的点位速度；③高精度：点位速度值精度优于3mm/a；④多种解：IGS站点具有至少3个不同国际分析中心的速度值，并且残差小于3mm/a；⑤均衡性：站点尽量均匀分布；6)精度一致性：站点间坐标精度应在同一数量级，并且速度值的精度也应在同一数量级。 / z" V, O0 b; [
本实验采用5个C级控制点进行七参数坐标基准转换，通过测量经坐标转换后的DSM中均匀分布的高程点坐标，与水深测量得到的水下高程点一同导入软件CASS中完成数据的融合，进而可对测区进行地形分析。
2 e! E& c2 b8 o+ O! {% x2 S, ~; a5 c 二、 $ T% x7 a+ K+ s9 g w0 C0 E( c/ T
测量技术方案 ( W% a t: N8 ~! ^+ a
本实验区域选择位于山东海阳与即墨交界处的丁字湾为测区，丁字湾为一“丁字形”的海边湾区。由于测区泥滩分布较多，受潮水冲刷影响，大部分地势比较平缓，从水流向外方向，地势缓慢增高，水流区域由于有河道，所以地势很低，测区西北部施工区域地势偏高，较为复杂，西南部分为陆地部分，地势偏高一些，起伏不定。
2 d* b# B5 A7 z3 X 在丁字湾区域开展的地形测量实验采用2种数据采集方式：陆地部分采用无人机摄影测量、水下部分采用单波束测深仪进行测深，最后通过数据融合得到整个丁字湾测区的地形数据。 & a% H4 }" A' C6 n3 ]8 E' I
⒈航测数据采集 _- f& }4 ~: q! e4 `$ K
⑴确定航飞区域 0 @& Z# T2 o; p8 F
使用专业地图软件确定航飞区域，根据航飞区域规划每个架次的飞行区域，根据精度要求规划航线，设置飞行参数，保存每个架次的区域文件与规划好的航线文件。 + ~4 K n) s, F! t/ j/ T8 B8 p
⑵航线规划与飞行参数设置 , {& J* i8 r1 R7 Q, Y A H8 Y, x
根据实验要求设置参数：比例尺为1∶2000，航高为500～600m，旁向重叠率为70%，航向重叠率为80%，地面分辨率为7cm。按照此参数在每架次区域内生成航线，保存待用。 ) ~. b( r( |7 w8 y* x ]1 N
⑶布设像控点
' h3 j3 @- v( p7 i0 u1 U- J 按飞行区域布设像控点，在地图软件内进行标记。由于使用PPK后差分系统，可适当减少像控点布设，所以约1～2km布设一个像控点，中间布设部分检查点(图1)。 
i4 U7 Z6 H8 I5 [4 `7 V) r 图1 像控点位置分布 5 p( W% b3 ^9 r: }1 ]6 v: n
⑷野外航拍 ( X3 Y2 z. C! t/ F. `8 q. o1 a" t
实施野外飞行作业前，选择视野开阔的平地作为起飞地点，观察风向风力，选择合适的起飞时间。起飞前在地面先架设GPS静态基站，调节飞机相机参数并自检完成后进行数据采集作业。飞机按规划航线完成作业后，检查数据与PPK是否对应。
5 v/ e8 |3 b( I h ⒉水深数据采集
" s2 R& ^) y3 Q( `" m5 v 本次水深测量平面控制采用CGCS2000坐标系统，中央子午线为120°E，3度带高斯投影。高程采用国家1985高程系统，施测前选用5个C级控制点进行了静态测量，并进行了基线解算，三维约束平差，获得测区高精度七参数，选择测区附近的控制点进行检核，单点校正，并进行坐标比对，满足测量规范要求，经质检符合精度要求。坐标校对中误差见表1。
2 P: F- w2 O3 ]% n( ]9 b3 E 表1 使用七参数改正后坐标校对中误差分析(mm）  . b0 e( ^3 S/ v: q( f2 ~
⑴确定测深区域布设测线 0 N: w6 T1 p6 x% O2 R2 `6 [. I7 q
按技术要求布设主测线：间隔50m，主测线方向垂直与等深线方向。检查线与主测线垂直，不少于主测线5%区域测线。具体测线布设见图2。 
& }8 Z2 y: {' l" B7 [* e 图2 测线布设 0 u: q1 Q* q# z. V6 |- A1 T2 Q' y' M
⑵水深改正 2 H* V4 `" j9 m! A
声速改正：为提高数据精度，采用由浅至深再由深至浅的方法投放声速剖面仪，将平均声速曲线直接用于声速改正。②吃水改正：精确量取换能器吃水深度，输入测深仪进行吃水改正，换能器吃水测量前后各测量一次。
6 @9 ^/ S% H6 c, r. s' ^, {5 a ⑶导航定位 # ]% Y; t4 b" @. o) e; l8 N
由于近岸测深，GPS可达到实时固定解，所以采用RTK进行定位，设置GPS接收机输出GGA数据，同时采集固定解，不需要验潮。操作员根据显示器上的测线与船型定位来操作船进行测深，为保证数据质量，船速保持在5～7kn。 # y; F" v( G8 m$ y* S# G
⑷水深施测
% z1 ~# `7 A' g+ S$ Q7 v 测深仪在测量前先将声速仪的声速参数输入测深仪，使用人工测深对比板测量水深值，与测深仪所测得水深值进行比较，对声速进行修正。在测深仪的测深界面上，通过增益方案实时调整合适增益值进行水深校准。 \+ q. q$ A4 {! p- ]" m
检测线的测量方法及测量精度与主测线相同，检测线的定位点间距根据测量比例加密至在规定范围内与主测线保证有重合点，每一条检查线与主测线的交叉点都进行比较。
6 [4 m/ f6 y3 E. q% S 三、 7 q4 T \, e6 Z- f$ N8 Y! o& ^
数据处理及结果分析
+ L: h& M. Q3 b ⒈航测内业数据处理
7 e7 W2 W. V% B 数据采集完成后，检查照片数据与PPK数据，照片数据质量是否合格，照片与PPK是否对应。检查无误后使用PPK解算软件对PPK数据进行解算，得到目标坐标系下的每张照片的精确坐标。最后使用Pix4D软件进行数据处理，通过像控点矫正，得到DOM与DEM。将得到的DOM与DEM导入三维测图软件建立DSM进行高程点提取(图3)，5m一个高程点，分布均匀，保证所选点具有代表性，可体现地貌特征。 
1 V! E& H' p4 }3 r* S 图3 软件提取高程点界面 " @, h/ l3 e5 j
所有区域高程点提取完毕后导入CASS进行检核，与检查点进行对比，并且每片区域接边重合区域高程点进行对比，对比重合及相近的点。结果表明，差值全部在15cm之内，确认无误后融合。
0 r4 ]. T. F9 d ⒉海测内业数据处理 ! S) ?% s$ C7 C
使用HaiDa水深处理软件对水深测量数据进行计算与处理。对DGPS测量自动化系统所采集的水深和定位数据进行内业资料后处理。在软件内设置投影参数、比例尺和图幅范围等基本参数后，根据声速改正和吃水改正测量数据对原始数据进行改正计算，在主测线上均匀选取水深点，删除多余的水深点，最后将处理后的水深数据导出。水深数据处理结果见图4。  3 w: I2 O$ C. p6 g# M6 v: l
图4 水深数据处理结果 / }( m8 [7 X3 ` l6 Q6 Y5 D8 t9 z
⒊结果分析
+ F* ]6 K" A0 z 对航测数据和水深数据处理完成后，首先将2组数据进行融合，将处理结束后的水下高程点在CASS中展出，并导入航测经坐标转换后的高程点。2组数据融合后对重合部分进行对比，并与检查线交叉部分的点进行对比，经过对比得出，重合高程点及相近高程点差值均小于25cm，其他区域如有高程突变情况，联合实际地形检查合理性，是否有陡坎或陡坡等地形变化，最后将融合后的高程点导出(图5)。  ' e h0 i( V# o7 |6 b
图5 数据融合后的高程点 / p3 N/ s) i3 m% A6 P
原始数据资料中测区最大高程点高程5.91m，最低高程点高程-15.24m。本实验测得测区最大高程点高程5.83m，最小高程点高程-15.2m，从数据融合后的等高线图来看，整体地势比较平缓，测区西北部地势偏高，地形较为复杂，西南部分地势偏高。结合以上对比数据可以得出，测区地形与实验前已有资料数据基本符合，2组数据在融合前精度均满足测量标准规范要求，融合部分数据高程误差值小于25cm，满足相关的精度要求，经实验验证本方法可用于实际的工程应用中。本次实验共用时8天，其中踏勘布设像控点2天，无人机摄影测量外业用时2天，水深测量外业用时2天，内业数据处理用时2天，相较于传统外业测量节省了大量的人力和作业周期，显著地提高了工作效率。
! {/ _% C# z9 S1 H- P5 M* ] 四、 $ g/ R/ f) n+ t# Q0 b: f, |9 w
结语
- T P! g' J6 |9 C 地形测量是工程应用开展前的必要准备。通过地形测量，可以获取地形地貌空间数据，其具有可量算性，具备空间分析基础，为工程规划、设计、施工等不同阶段提供了基础数据。我国近年来正在大力推进对海湾地区的开发及保护，面对海湾地区水池洼地密布、潮流影响大、交通困难、地势复杂等环境因素，传统的野外RTK测量作业难度大。本文采用无人机摄影测量技术与单波束水深测量技术相结合的方法降低了外业工作量。在高潮时采用单波束水深测量方法进行测量，在低潮时采用无人机航测方法进行测量，确保了重叠区域便于精度验证，很好地解决了潮区测量。通过数据融合，使得陆地与海底的地形数据相统一，可以更好地应用于海岸带保护利用与规划、海湾生态修复保护、海洋工程建设等工作，提高了作业效率，降低了作业成本。该方法在其他海湾地区地形测量中也具有一定的应用参考价值。
% u+ b& Y- o k1 U6 w 1
: g7 B5 z9 y+ Q% f5 M- k) b9 G8 u: H END
# o J8 X! j0 ~3 W/ B 1
* u0 @, k2 k7 k! m7 Z5 B7 e8 |' k( q Z, r 【作者简介】文/禹云亮 侯健 魏明刚 曾淑辉，分别来自中山市海洋与地质环境监测站和青岛中海基业海洋科技有限公司；第一作者禹云亮，男，1978年出生，高级工程师，主要从事海洋测绘方面研究；文章来自《海岸工程》（2020年第3期），参考文献略，用于学习与交流，版权归作者及出版社共同拥有，转载也请备注由“溪流之海洋人生”微信公众平台编辑与整理。   
8 t3 \( y; b5 n+ s 相关阅读推荐 海测论坛▏赵建虎等：海底地形测量技术现状及发展趋势技术交流▏水位改正方法对水道测量精度的影响分析论文专区▏海道测量规范与IHO标准的比较研究海测技术▏GNSS支持下的无验潮测深模式优化经验交流▏RTK三维水深测量数据处理优化方法探讨溪流笔谈▏水深测量外业探测前的准备工作    
7 D7 L) F- J! w+ F9 v+ }0 X 公众号 5 o5 a1 A3 i4 f+ X% W- `- ~
溪流之海洋人生
6 \! Q9 i+ w3 z/ d 微信号▏xiliu92899  / j% Q+ z6 g; v- U! p8 K$ A; B, D
用专业精神创造价值 + U7 {; q: k/ J8 r* W, b
用人文关怀引发共鸣 您的关注就是我们前行的动力
% }4 p- R: |4 s8 O- ~/ ] 投稿邮箱▏12163440@qq.com  
8 @. @2 _' p( W8 o) ~. g) z4 U3 n7 n& z& u. I$ I
8 L7 v4 `7 G& C4 C5 }; r
5 l; n2 g: V: F( Q
4 P( j0 J+ b% K! N0 B5 ?/ c
| 
