|
9 M7 k( I& a8 Z j5 p' H / Q* S" a- A1 u! _" _
VQ-840-G
% x1 Y% j6 Q+ [4 C P 无人机激光雷达测深应用潜力剖析 8 f& l4 _( Z6 s* o9 M: i
项目背景. + t$ L3 t6 z0 D# P
随着国家对海洋测绘研究的重视程度的提高及国内外对机载测深激光雷达技术的不断研究，使用无人机搭载测深激光雷达进行水陆联测已成为一种新型的技术手段，并将随着硬件及软件技术的发展，走向大规模的应用。 9 ]! b' @& x$ r j/ Y# p
现有测深技术介绍分析. . {, E5 y6 r6 s' C: N. w \
伴随着国民经济发展及国防建设需要，我国对湖泊、水库及海域的水上、水下地形测量任务需求不断增加，尤其近几年国家所提出的海洋强国战略，重点发展海洋经济，提高海洋资源开发能力以及维护国家主权、安全和发展利益，都将水域或海洋测绘领域研究提升到了战略性高度。而对于近海地区研究更成为重中之重。它对于渔业发展、港口停靠等将更为重要。而目前测量的手段主要包括：
6 O; r Y/ E8 {% Q1 y: r7 p3 P; P 传统的GPS测量手段
0 K3 E4 H( g- F2 e" k, k 技术痛点：需测绘人员坐船到达近海地区或行走于搁浅的浅滩地区进行测量。这种方法效率低、速度慢，可测量的深度有限，无法满足大面积测量的需求。 F# l, m, k- F- Q& j8 K: U: m
水上光学、水下声学的联合探测方法，可通过无人船或者有人船实现
" M: x- g0 x5 p; p( O+ O. y 技术痛点：虽然可解决水上、水下同时测量的问题，但是当靠近近海地区时，多波束将无法满足工作条件，无法获取水下地形；而远离近海地区，虽多波束可满足工作条件，而水上光学由于远离岸边也将无法获取完美的数据，因此给人们带来了工作的难度。 ; D& F: O* z- |4 }# m2 t* T/ k1 ?: j3 H
大型机载激光雷达测深技术手段 , a+ `! U9 B, ]. L. L
技术痛点：虽然可全天候、大面积、通过主动式遥感手段获取水上、水下地形数据，但由于目前主流的机载测深系统均为体积大、安装困难、需要大型飞机才可飞行的要求，导致难以大面积推广使用。
+ V# e% d0 K: B2 Q/ I' z 因此，无人机激光雷达测深系统将是最理想的选择，也将成为未来水陆联测最为重要的技术手段之一。RIEGL VQ-840-G无人机激光雷达测深系统则是一款比较理想de的测深系统，成功引进中国后，在青岛胶东湾海域进行了首次无人机搭载VQ-840-G的测深数据获取。通过完整的数据解析后，实现了与多波束地形的对比与分析。 2 [5 U7 g' s' O2 \. c4 x% W% a& b; }
4 @ G& A+ Y" {3 E: d1 U* G
实验设备. ; j9 g# A- B9 f6 r3 N! k0 K' d
采用哈瓦战斧旋翼无人机，可搭载重量为25公斤级别设备，是VQ-840-G设备理想的无人机搭载平台，单架次可作业30分钟。 , @8 a& z0 v( x# p# ?2 `+ u: e- E5 z
3 g. G3 b. Z, p& O7 J0 g
飞行参数设置. . T( b: V: f3 R) N: C; U
6 c9 L* Y; U0 [/ G- w
/ r" }: g( F& q/ { " u0 X: J) j: b
航间距
& m: A0 t- \; k& j, c' _ 15m / K( ^/ b) f) K4 x- W7 N+ r
平均点密度 7 O0 X% o. w& w5 }! ?( o! [: F
320点/平米 8 h3 N+ E( o, F
采集频率
# N( I4 n; v. q 50KHz ! C2 o& E0 n7 g# j4 Z7 T1 f0 Z
线频率
) ~. f; ~4 J E |8 D! P 100线/秒 & _( w2 h+ r1 f% ]6 C# w( X' W" \1 h
速度 7 s# v+ h t: Z3 ]* Q
5~6m/s
5 K4 J* O: Y1 ` 飞行相对行高
0 z/ ?6 f% L% [: p. O$ y 50m ; l5 v$ z$ S5 }8 V
视场角
6 p3 G# k& Y( Z" t5 n1 K6 A! c 40°
7 h i5 E& t$ C* y 激光发散度 9 `& A7 g$ \4 V* N6 R
0.001rad
! v& p& H* W2 `8 U 水质情况 2 d1 v7 R* i. G* L& R
中等
, Y9 ?$ Z, b! ^/ o% v% J" i: u6 h( v2 [; l 照片曝光间隔 8 M' b/ e5 O9 S$ V# h
2s
/ x: o' r( @7 Q4 H$ V ) p) g7 H' s! ]9 I% K! t
0 o$ ~$ ~1 ~$ I# b1 D
( g1 R# _ Z( Z. ^ 数据处理流程. 原始波形数据解析噪声去除，主要针对点云数据中的噪点的快速删除航带拼接，采集的多条航带间中的匹配工作水面分析，提取出准备的海面数据，用于进行深度计算测深改正，将根据水域折射率计算实际深度最后为DEM、等高线等成果构建
2 A: b3 I1 D! T: S* ]' @ 数据处理技术重点. 波形数据处理
: O g& ~! E7 m' H+ |3 u 水深处理中，波形数据的提取是至关重要的，波形提取将决定水深测量的深度能力。因此本次处理中，我们将使用波形及未使用波形处理所产生的深度差异进行了对比。未使用时的深度约为6米，而使用全波形处理后深度可达7.5~8米。 ) O) u; @" `6 Y+ G- Y# G) V
( U- F+ n ^" [1 h# n$ _4 w ▲未使用波形情况下解析的水下成果
d9 s7 H' Y) }. p8 \  ! L1 h4 X3 \. l! a3 U
▲使用波形情况下解析的水下成果 ) }; E$ q% x- g% I8 |& e$ q) g0 p
2.噪声去除 8 t. y1 o( ~# O; F4 w
对于测深激光雷达系统，由于激光能量强、入水折射等因素导致测深数据，噪声比较大。因此我们需要联合使用多种噪声分离的算法去除噪声数据，准确的分类出有用的水面和水底，才可以得到准备的深度数据成果。
% r% q G) Z4 f9 \ 本次我们使用了三种混合剔除噪声的手段来完成，包括离散点分类法、反射率差异分类法及手动分类法三种方法完成。
$ a* X& p/ R% r% K 
- Y% K& z5 q5 ^ }+ @+ E* f ▲噪声去除前的效果图
& U3 D( n- m u; |/ @5 Z7 u( v 
" Y1 w. ~: x3 B0 ^4 M ▲噪声去除后的效果图
8 i# q% @ V P9 r2 @- V8 L 3.测深改正 ) p) w( w1 @4 V: u0 D- w' ]/ @
由于采集的激光雷达数据，进入水域后会发生折射现象，而这种折射现象则需要根据实际采集水域的情况进行调整。水域的折射因子与水域温度、激光波长、盐度及深度变化而变化的，因此为了计算水域测量的实际精准深度，我们必须先计算出所测地区的折射率。 - F8 S/ n; @: |7 a
使用此折射率对采集的数据进行折射纠正，纠正前的测量深度为7.5米~8米，纠正后为5~5.5米。
8 e0 x' u7 D* s0 G 
: C7 y: ^1 G, w, J9 L ▲改正前的测深深度（7.5-8米） 9 _* D/ j, ? `/ d+ @
 4 W6 a# S8 a& C. i7 u
▲改正后的测深深度（5-5.5米） , G/ H5 m; f4 I) s0 K9 X# |- W/ K
数据成果. 6 {4 \' W z2 J* ]& V
 : V D1 u" a$ P& N
水陆联测点云 % p6 H0 ~7 O/ Y4 C9 X" Q
& Y8 A. Q. G. U `: k5 w/ W5 j ▲不含水面数据的水陆联测点云
: i, W' `; Z5 h1 U+ I; _; g  & e0 ]5 b+ s4 G5 d6 w) e v
▲真彩色水陆联测点云数据 / [+ Q* ?( o+ v
0 c5 y9 H" T. z+ ~( u* ] ▲水陆联测地区所构建的DEM数据成果
; v) a9 X) }) K# z2 z0 b7 Z  * I) w: P" Z: ^6 ~- D7 e
▲海底地形数据 ; @" c' v7 x H( @, s9 f ^
成果精度比对.
- _2 o+ p( w% o" Y 精度对比手段为使用单波束测量成果与激光雷达成果进行对比的方法实现。而此次单波束成果为山东科技大学采集的单波束测量成果。 0 B* m. p, O+ N3 }
由于所使用的为不同时期采集的成果，我们对单波束数据进行相应的潮位修正、基准改正，而后进行了海底数据的精度对比。
& Y2 ^8 N, \' l1 R6 R# ]. D$ p+ w8 S: G 
/ o! N$ E, L# {3 G ▲ALB海底点与单波束点相对分布
/ o$ a% A/ Q2 c. `! x, g 其中红色的点为单波束点，灰色部分数据为激光雷达测深数据。在重叠区域选择18个单波束点与其最临近激光雷达水下点云数据比对，计算其高差，并统计18组高差的RMSE值。
/ ?" \" t; r5 a/ L7 I) u' _ , s0 J& k g& A% g) J, E; d8 v
& g. Z. L# u& c$ m2 Q7 [8 `
9 V4 R: {7 c% i# n. U 点号 - i1 m' U* s' n! [! U- ]' z/ \
与最临近点平面距离(cm) 4 e9 B/ I. ]+ @8 Z) } w: q7 g' ^$ Q
高差(cm) 0 K0 |+ h6 h @9 [. c, I3 h4 W% L
RMSE(cm) 7 S3 E: ~+ ?7 k T8 w
1
" V! k5 f8 Q% n 0.6 , _% v1 {5 q1 T; E5 R3 W& m; ^
16.8
: d/ p+ p) e, K! P N2 N 20.7
3 L" A7 Q, m( Q3 q2 | 2
' d3 B l+ J& ^8 M 4.0 $ G G/ m3 a- H5 v) {% M- G' ?" B
37.1
0 `' F8 X, o' a E( x 3
8 j) H- P+ X" u! W/ M 1.1
& e' c$ B( o: @/ b2 h2 l 10.9
# t( n2 B8 m1 W 4 , |5 a# v. o2 a, [4 Q1 @# M: }
2.3
$ Y& {+ z& X& |3 F* a9 u1 ? 16.4
+ U, ]9 U, Q. G# N5 h: A0 x8 a 5
/ i1 Y8 b8 }1 l: `) |8 v2 ` 2.8 ) o( l1 E l! x
14.8
" L5 M' f3 ^& d. W 6
& U9 m0 H1 B. A1 i 2.8
- s8 ^) @0 M3 X& G5 x# G 14.9
s# B {* g% _6 v3 P3 ]9 S 7 8 H! O8 g* f) f- E+ F
7.2
/ R1 t9 y* {; p3 k% K/ ~- e 7.8 ! b' s4 P5 F/ l+ h: `, W
8 & z, {' s, w" [) N
8.9
3 L, I" z& l$ U5 s/ @ 4.8 8 c% K3 a% M* G4 C8 f- ~* P% s
9 2 g5 Y& G6 ?) @/ a
5.0 , s; z& G) B' N7 l3 _& B1 D- P' k
8.4
% S g! ~! {9 e. B1 N! N* l+ ^8 l 10
0 w& i$ u+ S/ b' y) W6 T2 L' } 4.7 + D) C, C; H# s/ o8 a
8.8 + u2 q6 h7 D7 U& F" H" S+ J1 q5 @
11 ( T$ m7 }/ _0 Y5 W" p D! T
3.9
- y6 \. |4 i% z9 H 9.1
# M$ H) `/ h: h3 G l' Z; H4 ~8 g/ ` 12
" \5 A u/ m9 ~+ y0 W 4.3
7 e. |( G1 z8 @" c 19.0 0 `6 v& l2 V$ v- C
13
+ S" | C- ~& o 1.1 ?2 Y7 B& z+ P- `" t
5.8
) d$ K6 G/ n- ~( |6 w 14
5 ?! I7 k6 p5 t' B 1.1 ! ?8 h/ K r- h0 b3 L: A# M5 X
46.4 r% f3 C3 e) m3 t. o
15
' O, _$ L# U$ G, g 3.9 ) V6 D. S8 M! G* @' S
10.7 + U: U$ p b; J) Q" z5 b, D
16
7 v- z9 a9 o$ R4 T 8.2
5 I( ~% w* _: ~ @ 34.9 : D% D4 P- T/ V* E+ B
17
+ Q9 N: s5 o" d 2.7 7 R& x. [9 ~6 m) }
3.6
" I) p3 n! p0 Q# F 18
% n3 j1 o! n* z% ~! V1 J 8.4 7 k" |8 I8 c/ H w# p2 n( e# W
31.9 ! `- C: x) B) r8 J4 I1 T; ^6 A
. U" H! Y) S" X0 T5 v) I, P 6 ~2 h! j2 |( h! R5 z8 r
6 C4 M$ d: @- Z! }3 d* l
结合以上分析，可以看出RIEGL VQ-840-G无人机测深激光雷达，能够按照中国海洋测绘的要求生产最终的数据成果，同时使用山东科技大学现有单波束数据进行精度比对，虽然由于测量方法等的差异等带来的对比难度和准确度的争议，但整体结果较为满意。
: f7 |3 x9 u. y* e: ? 通过本次测试充分证明了无人机激光雷达测深技术在中国黄海、东海、南海等海域可以很好的完成海洋测绘应用，是一种具有巨大前景、可对浅滩近海域进行周期性数据采集更新、全天候地、快速响应的海洋测绘新手段。
7 e, I' Q, \' I: x: k+ P( w5 [1 e( ]2 x/ |
8 ~, \$ @+ f8 l, J5 `8 ]- `7 Y" q1 Q8 \# j) b
" k+ n' j- e* M$ f' ], X
|
