论文推荐| 赵荻能: CUBE曲面滤波参数联合优选关键技术及应用

[复制链接]
) h' c; o0 b9 @/ K: k

《测绘学报》

1 N% M K8 L, v/ d* _

构建与学术的桥梁 拉近与权威的距离

; M F0 a& @& I. K3 Z# G3 Y: x

CUBE曲面滤波参数联合优选关键技术及应用

- P1 h. W- ~9 f- d% ~. E

赵荻能1,2, 吴自银1,2, 李家彪1,2, 周洁琼1,2, 张田升1,2, 刘洋1,2, 朱超1,2, 余威3

: W& i0 [! s, F7 |2 B

1. 自然资源部第二海洋研究所, 浙江 杭州 310012;

3 j$ s. l+ t6 {) |+ a* M. g

2. 国家海洋局海底科学重点实验室, 浙江 杭州 310012;

, A4 x7 {8 o- C, l

3. 国家海洋局南海环境监测中心, 广东 广州 510300

) X* V! d( c: {; v) f, A

收稿日期:2018-03-04;修回日期:2018-07-04

& v8 R: L+ {5 b

基金项目:国家自然科学基金(41830540;41476049;41606057);卫星海洋环境动力学国家重点实验室自主课题(SOEDZZ1802);科技部科技基础性工作专项(2013FY112900);浙江省公益性项目(2016C31120)

* {! a, y4 L9 Q# e" k# U) N

第一作者简介:赵荻能(1990—), 男,博士,助理研究员,研究方向为海底地形地貌探测。E-mail:zhaodineng@sio.org.cn

! `& k- k) ]9 x2 N2 V3 p$ l) H

通信作者:吴自银, E-mail:zywu@vip.163.com

, Q1 i6 [! i$ r4 H- v

摘要:CUBE(combined uncertainty and bathymetry estimator)算法是国际上主流的多波束测深异常值自动探测与处理算法,在国内外被广泛应用,但对其核心算法和参数知之甚少,不利于该项技术的国产化。本文详细阐述了CUBE算法的基本原理、数学模型、关键参数和处理步骤,进而建立了CUBE曲面滤波参数联合优选方法。通过选取典型地形区、参数试验、对比分析等步骤完成参数的联合优选,并用台湾浅滩实测数据进行了验证。结果表明,优化后的参数可有效提升多波束数据自动处理的精度和效率。本文成果可应用于国产多波束测深处理软件的深化研发以及多波束实测数据处理。

7 b' v- p' O9 z7 @' e0 l

关键词:多波束测深 CUBE算法 参数优选 算法原理 数学模型

& Z% v, X. g# B: }- N' S/ H+ o

The key technology and application of parameter optimization combined CUBE and surface filter

7 q% @4 P3 U7 V. k, M

ZHAO Dineng1,2, WU Ziyin1,2, LI Jiabiao1,2, ZHOU Jieqiong1,2, ZHANG Tiansheng1,2, LIU Yang1,2, ZHU Chao1,2,YU Wei3

2 G0 [" R! I& o

1. Second Institute of Oceanography, Ministry of Natural Resources, Hangzhou 310012, China;

* p1 j: O* l$ G2 ?

2. Key Laboratory of Submarine Geosciences, State Oceanic Administration, Hangzhou 310012, China;

2 [) c+ o) f/ K1 c; F& M+ X. {

3. Environment Monitoring Center of the South China Sea, South China Sea Branch of State Oceanic Administration, Guangzhou 510300, China

- V& t$ t' j4 A |! U8 ^' O: u4 j

Foundation support: The National Natural Science Foundation of China (Nos. 41830540;41476049; 41606057); The Project of State Key Laboratory of Satellite Ocean Environment Dynamics, Second Institute of Oceanography (No. SOEDZZ1802);The Fundamental Project of Science and Technology, Ministry of science and technology of Peoples Republic of China (No. 2013FY112900); The Public Projects of Zhejiang Province (No. 2016C31120)

0 M! R* s! v& |! K6 U5 @& K

First author: ZHAO Dineng(1990—), male, PhD, assistant research fellow, majors in seabed topography detection.E-mail:zhaodineng@sio.org.cn

5 f C- A r+ H8 d" z

Corresponding author: WU Ziyin, E-mail:zywu@vip.163.com

3 _' C5 C& C. A4 w' A1 g

Abstract: CUBE (combined uncertainty and bathymetry estimator) is a widespread algorithm in the world to automatically detect and process multi-beam sounding outliers. While, little is known about its core algorithms and parameters, which is not good for the localization of this technology. In this paper, the basic principle, mathematical model, key parameters and processing steps of CUBE algorithm are elaborated. Then a method of parameter optimization selection combining CUBE and surface filter is established. And an example is presented, in which optimized parameters are selected by choosing typical bathymetric area, parameter test, comparative analysis and etc. And the optimized parameters are verified with the observed data in the Taiwan Banks. The results show that the optimized parameters can effectively improve the accuracy and efficiency of automatic processing of multi-beam data. The technology in this paper can provide reference for deepening research and development of domestic multi-beam sounding processing software and multi-beam data processing.

% F+ s5 Z' ]5 [- ^

Key words: multi-beam bathymetry CUBE algorithm parameter optimization algorithm principle mathematical model

+ _" k/ B, z+ o( D- I6 e

多波束测深技术出现于20世纪60年代[1-2]。这项技术实现了水深测量由点、线到面的飞跃[3-5],提高了人类对于海底地形地貌的认知水平,极大地推动了海洋地学的发展[6-7],尤其是直接推动了海底地貌学由定性向定量化研究的跨越式发展,并且促进了中国在高分辨率海底地貌学研究方面取得重要进展[8-12]。

% I3 _! O; `, W/ F* _

与传统的单波束测深技术相比,多波束测深技术具有高精度、高分辨率和全覆盖探测等显著优点[13-15],但海量的多波束测深数据也大幅增加了后处理的难度和成本[16-17]。第一代商业多波束系统只有16个波束[18],新一代多波束系统的波束数高达上千个[19],尤其在浅水区作业时,每小时可产生上千万甚至上亿数量级的测深点[20]。因此,如何对新一代多波束测深系统探测数据进行高效处理是重点和难点技术问题。

: k! V% V# X: m0 m3 F

针对多波束测深数据自动滤波处理问题,已经开展了大量的研究。如:通过建立深度和波束角阈值门限进行滤波法[21]、通过检测测深点之间的角度和梯度[22-24]、建立局部地形模型[25-29]、模拟人工处理过程[30]和基于统计学[31-33]等方法检测异常值。但这些自动处理方法主要通过一定的准则对测深数据进行分类,处理效果有时并不太理想。

8 a. P: h* E, C& C' a3 B: `

CUBE算法由新罕布什尔大学的B.R.Calder博士于2000年首次提出,之后分别于2001年和2002年引入HGM模型和多水深假设追踪和干预模型,其核心算法得到了较大的改进和提升[34-35]。该算法在2001年的美国海道测量会议上首次公布[34],逐渐引起了世界各国海道测量部门的重视。在美国,NOAA下属的海岸调查办公室将该算法写入标准外业程序手册中[36-37];在英国,UKHO下属的水深数据中心研究得出,通过该算法得到的水深曲面适用于海道测量部门[38],NetSurvey海道测量公司探讨了该算法的水深假设过滤方法[39-40];在智利,海道测量办公室优化了算法中4个关键参数,并将其应用到多波束数据处理中[41];在韩国,海事安全研究中心对比了人工处理和CUBE算法自动处理后的测深数据集,认为该算法处理多波束数据不但高效而且准确[42]。同时,国内不少专家也对CUBE算法的原理和应用做了大量研究工作[42-47]。当前,CUBE算法已被多种国际多波束数据处理软件采用[28],但尚少见国产软件采用该算法。究其原因,或许与该算法过于复杂有关,且其关键参数晦涩难懂,未能依据测区地形特征优选一组最优CUBE参数。

% g% |- o$ y! c7 r

为进一步推动CUBE算法的底层深度研究与国产化应用,本文通过对大量国际文献的系统梳理,结合作者团队的多年实际工作经验总结,详细阐述了CUBE算法的基本原理、数学模型、关键参数和处理步骤,进而提出了一种CUBE、曲面滤波参数联合优选流程,并使用实测数据进行处理验证,表明优化后的参数可有效提升多波束数据自动处理的精度和效率。

3 q# e& D3 i5 s* J* S0 x' w' i

1 CUBE算法基本原理

5 n6 K1 o; \+ v

图 1为CUBE算法的工作流程,可概括为以下8个主要步骤:

k2 |$ a- l1 T+ v* A8 x

图 1 CUBE算法单格网点运算流程Fig. 1 CUBE algorithm single grid operation flow chart

, i, n) q2 C' w

(1) 通过HGM(Hare-Godin-Mayer)模型获得每个测深点的一阶近似测量不确定度,其基本思想[48-50]为:通过深度与原始测量(及其不确定度)之间的关系建立前向传播预测模型,最终得到测深点垂直不确定度,类似的也可得到水平不确定度。垂直和水平不确定度计算的全微分方程参见文献[44]。

& \6 E; x( W- j) i' U- N# {. j

(2) 应用式(1)计算每个测深点的影响半径

+ I! V9 S) W$ r4 o

(1)

6 J2 R* a8 e4 `/ w, s- P/ d2 |

式中,rinflu为测深点的影响半径(m);Gres为格网分辨率(m);σvertmax为用户选择的IHO S-44测量等级所确定的水深不确定度[51];σvert为测深点的垂直不确定度;σhz为测深点的水平不确定度;μdistexp为不确定度随距离增加而增大的比例系数。

3 b3 J/ Z# B7 `' X8 m7 k

(3) 应用式(2)计算格网点的捕捉半径

; k4 I! P$ L9 E/ k9 @, ~, g

(2)

: m. F! w- |% A( W$ S

式中,rcapdist为捕捉半径(m);μcapdistsca为捕捉距离比例系数;Gendepth为格网点期望水深(m);μcapdistmin为最小捕捉距离(m)。μcapdistsca和μcapdistmin共同控制格网节点吸收的测深点数量,值越大,格网节点将吸收更多水深点。每个测深点需同时在以上两个半径内才能参与相应格网点的运算(图 1)。

0 K1 D# s9 _' L7 O# c; |& L# P

(4) 应用式(3)将测深点3个维度上的不确定度(二维的水平不确定度和一维的垂直不确定度)融合并传播到格网点处的一维垂直不确定度

* r3 [7 P. e" U/ J" N7 }9 E0 B/ k

(3)

$ s2 o% C; V& y

式中,σp为测深点传播到格网点后的垂直不确定度(m);Sdist为测深点与格网点的距离;μhes为水平不确定度传播比例系数。

# O% L V7 J/ b t( q2 Q2 ^

(5) 将测深点依次插入到队列中,按水深值大小进行排列。当队列满时,输出队列中间位置的测深点参与下一步运算,并继续向队列插入新的测深点。队列排序可有效延缓异常值输出[52],而最新版本的算法引入了单格网点多水深假设的架构,输出的测深异常值将通过构建新的水深假设与正常值孤立开来[53]。

, ~' S) C- |! Y$ E5 o8 e* O

(6) 寻找当前测深点作用的格网点内最合适水深假设。存在两种情况:①如没有水深假设,则使用当前的测深点建立新的水深假设;②如已有水深假设,则计算格网点内所有水深假设与测深点之间的差值来确定最接近测深点的水深假设,称之为最适合水深假设[35, 53]。运用式(4)计算测深点和最适合水深假设的归一化差值

3 Q' k' M [" L+ k

(4)

8 Z" X" ?& n9 p. ?5 K) e

式中,en为归一化差值;SZ为测深点水深(m);2);σvert2为测深点垂直不确定度方差(m2)。

' H1 K8 }% }8 t! m

应用式(5)计算贝叶斯因子Bn

2 {) n( P. x8 y0 p6 L, ?( F8 z

(5)

- K9 i3 A' t' }2 C7 q5 [/ O

式中,Bn为贝叶斯因子;hestimoff为估计偏移值,用于描述垂直方向偏差的显著性水平。增加hestimoff的值,CUBE算法建立的水深假设将减少[53]。

1 E9 g) C7 G# o( o; c! Y& r8 @

将其与Bn阈值和长度阈值进行比较,如都通过,则运行卡尔曼滤波(Kalman filter)更新水深假设(图 2);如不通过,则创建新的水深假设。

! i1 }: ^$ s. m4 ?+ k& _0 z

图 2 Kalman滤波流程Fig. 2 Flow chart of Kalman filteringQ为过程噪声;x(-)为当前水深估计;P(-)为当前方差估计;σp2为测深点传播到格网点后的不确定度;K为滤波增益;x(+)为更新后水深估计;P(+)为更新后方差估计;z为测深点水深值;σdilute稀释因子

3 \4 y6 I2 v& }1 R2 I

(7) 所有水深假设构建完成后,应用式(6)计算每个水深假设的假设强度σstren

# t9 D4 |- u5 S5 y

(6)

" c1 B5 I. j7 Q6 ?; _% v

式中,σstren为水深假设强度;εcurrent为当前水深假设包含的测深点数量;εnext为序列内下一个水深假设包含的测深点数量。σstren的值越小,则水深假设的假设强度越高,越可靠。

! Q$ }) \9 t- |7 l8 V% L( m7 ?8 P; h( t

(8) 使用消歧算法确定每个格网点的最佳水深假设[53]。

9 E8 h! Z. \* s2 r! E+ u! Z" k- v

2 CUBE、曲面滤波关键参数联合优选方法

0 ~* R( D7 ^5 `3 p

将CUBE算法和曲面滤波相结合是当前主流的多波束数据处理方法。核心关键参数的选取是该方法成功实施的关键,因此需进行CUBE、曲面滤波参数联合优选。

8 J% z, H; U5 ^+ Z9 C, W

2.1 关键参数选取

% ]( P R" e6 E9 z+ g3 Y' E

μcaptdistsca(捕捉距离比例系数)、μcapdistmin(最小捕捉距离)、μhes(水平不确定度传播比例系数)和hestimoff(估计偏移值)是CUBE算法的4个核心参数,通过设置这4个参数可组合形成多组CUBE参数组。由于这些参数之间具有很强的关联性,随意设置可能导致CUBE算法出现严重错误[37, 53],因此使用当前较成熟的4组参数组(表 1):CUBE default为算法的默认参数[53];CUBE deep常用于海底特征地形较少的测区;CUBE shallow常用于海底特征地形较多的测区;CUBE NOAA与其他3组参数的最大区别在于通过减小μcaptdistsca值,并设置

+ G( l m4 i6 _- O u& p9 x

,保证CUBE曲面构建过程中不遗漏任何测深点[37]。

/ J& k; \: y: K1 ]/ V1 I8 y) U

表 1 CUBE算法参数组Tab. 1 CUBE algorithm parameter groups

; h! }+ C# j4 y* f. Y3 s
$ F' V! i" w% n2 F
CUBE参数组类型μcaptdistscaμcapdistminμheshestimoffCUBE default50.52.954CUBE deep2022.953CUBE shallow40.40.52CUBE NOAA0.51.964
! l" j5 u( H7 U1 c" i& i/ D. {- j, Z
) Y/ P4 b7 I ^/ ?8 u" ?4 ~5 H

格网分辨率也是CUBE算法的重要参数,选择原则为:在满足CUBE算法推荐的每个格网点不少于11个测深点的前提下[53],尽量提高分辨率以捕捉小地形特征。而滤波窗口μfilter则是通过设置一定大小的滤波窗口(一般为标准差的倍数),自动剔除滤波窗口外的数据。

4 p P$ X2 ]6 R$ N4 f' n

2.2 关键参数联合优选方法

0 J; ]9 s* w0 a2 i

如图 3所示,CUBE、曲面滤波参数联合优选流程可概括为4个主要步骤:

B! Q$ L4 e' O6 o

图 3 CUBE、曲面滤波参数联合优选流程Fig. 3 Combined optimization flow chart of CUBE and surface filter parameters

) @ a& x! }: n. t" ]

(1) 选取能表征整个调查区地形特征的典型地形区作为参数联合优选试验区;使用测深点密度分析法获得原始多波束数据的水深点密度,得到试验区最优格网分辨率;按表 1所示,设置4个核心参数以形成多组CUBE参数组,生成相应的多组CUBE曲面。

$ m; \( |; j& u

(2) 综合CUBE曲面的多种辅助信息(水深、假设数量、假设强度和水深不确定度等),识别并定位异常水深区域,人工干预编辑算法无法正确判断水深假设的区域,而后更新整个CUBE曲面。

; T& R: b5 Q3 C8 S

(3) 以更新后准确的水深曲面为基准,设置并应用多组滤波参数μfilter,对测深点进行自动滤波,剔除位于滤波窗口外的水深点,最终得到不同CUBE、曲面滤波参数组合处理后的水深点集和水深格网。

; @& m- t+ k; ?/ L

(4) 使用传统的人工编辑法获取试验区清理后的水深数据集,并以此数据集为基准,分别使用测深点集分析法和水深格网分析法,将其与不同CUBE、曲面滤波组合处理后得到的水深数据集进行综合对比和分析,最终得到全区最优CUBE、曲面滤波联合优选参数组。

+ }; Q- ?7 ~- \

3 试验与分析3.1 试验数据

/ f/ S" }; l5 G

2016年7月笔者在台湾浅滩执行地形地貌与水文测量综合航次,使用多波束测深系统和ADCP进行同步观测。多波束测深系统采集设备及其测量精度如表 2所示。

" r+ j) ]7 [% Y6 m, j

表 2 外业调查设备及其精度Tab. 2 Field survey equipments and the accuracy

8 J* f0 ]9 n; ]! [
* u8 r: Y/ _# h
设备类型仪器型号测量精度多波束测深ResonSeabat 7125测深精度:6 mm定位NavCom SF-3050< 0.1 m航向±0.1°升沉IXBlue Octans Ⅲ±5 cm横纵摇±0.01°潮位RBRduoT.D|tide全量程5%表层声速Reson SVP 70±0.05 m/s声速剖面RBR concerto C.T.D±0.5 m/s
' ^" g. o3 l( }2 E
% n8 d2 T% C1 ]: [/ N$ l g

数据在采集过程中已完成安装误差改正和实时声速改正,并在后期进行了潮位改正。图 4为使用人工编辑处理法获得的海底地形图。该区域水深15~40 m,平均水深28 m,其西侧发育有典型的多尺度沙波地貌。测线间距约110 m,由于该区水深变化剧烈,导致沙波波峰处出现由于条带覆盖宽度不够而产生的数据空白区。在该研究区内选取典型地形区作为验证本文所述关键参数联合优选方法的试验区(图 4红框区)。

7 W1 {4 F! S& E$ v0 W

图 4 台湾浅滩研究区测线分布及海底地形Fig. 4 Trackline distribution and seafloor topographic map of Taiwan Banks

2 ~8 g/ m* O" V" r6 X# Y1 G

3.2 参数联合优选结果

% z% c. l& `$ d, y

曲面测深点密度分析表明,格网分辨率为1 m时,99.5%的格网点满足最小水深点数量阈值要求(图 5)。因此,设置试验区最优格网分辨率为1 m。

^# d$ L/ d0 [8 S& [3 n

图 5 网格测深点密度分析Fig. 5 Density analysis of grid sounding points

" X& R6 r5 y/ l4 ^4 b/ z3 ]

表 3总结了不同CUBE、滤波参数组合对典型地形区数据处理效果的影响,测深点集分析(S.A.)和水深格网分析(G.A.)结果具有较好的对应关系。当MACA+MRCR值较大、MRCA+MRCR值和N.D.值较小时(各缩写词含义见表 3注释),表明该组参数处理效果与人工编辑处理效果具有很好的一致性,即该组参数适用于该地形区。表 3中第1—3行显示,使用相同的CUBE参数和不同的滤波参数对最终数据的处理效果具有一定影响:当滤波窗口较小(1.5×STD)时,部分有效测深点被滤波窗口剔除,MACR值增大;而当滤波窗口设置较大(3.5×STD)时,部分测深异常值无法自动过滤,MRCA值增大。经过多次试验,得出该区最优滤波参数为2.5×STD。表 3中第2、4—6行显示,不同的CUBE参数结合相同的滤波参数对最终数据处理效果具有显著影响。综合MACA+MRCR、MRCA+MRCR和N.D.值,得出该区最优参数组为CUBE NOAA和2.5×STD。

) ^4 ~! p0 F4 Y- K

表 3 典型地形区CUBE、曲面滤波参数及处理效果分析Tab. 3 Analysis of CUBE, surface filtering parameters and processing effect in typical terrain area

' r9 Q' Q I0 d
% m- k6 ^ h4 m
序号CUBE参数μfilterS.A./(%)G.A./mN.D.MACAMACRMRCAMRCRMACA+MRCRMRCA+MRCR1default1.5×STD94.720.8314.0630.38395.104.8940.082default2.5×STD95.370.1824.0710.37595.754.2530.063default3.5×STD94.510.1954.9150.38194.895.1100.124deep2.5×STD94.930.7423.9390.38595.324.6810.095shallow2.5×STD95.850.1143.5370.50296.353.6510.036NOAA2.5×STD95.980.1023.3850.53296.513.4870.03注:STD——格网点水深标准差;S.A.——测深点集分析(soundings analysis);MACA——人工编辑保留且CUBE曲面滤波保留(manual accepted and CUBE accepted);MACR——人工编辑保留且CUBE曲面滤波剔除(manual accepted and CUBE rejected);MRCA——人工编辑剔除且CUBE曲面滤波保留(manual rejected and CUBE accepted);MRCR——人工编辑剔除+CUBE——曲面滤波剔除(manual rejected and CUBE rejected);G.A.——水深格网分析(grid analysis);N.D.——在95%置信区间下的格网节点的差值。
, i( b& a! w4 ]9 ? |4 C- v8 B
9 r k, Z: b7 y! v2 C( B: `# D; z

3.3 参数敏感性分析

, @3 w! n6 _9 \3 T

在获取研究区最优参数后,通过固定其中任意4个参数并按一定测试步长调节余下一参数的值(表 4),获取各参数变化对算法数据处理效果的影响。如图 6所示,除μcapdistmin外,其余4个参数的变化都可对算法最终的处理结果(MACA+MRCR值)产生较大的影响。不同参数对算法的影响存在较大差别,如μhes、hestimoff、μfilter与MACA+MRCR值的响应关系呈现两边近似对称的抛物线形态。在最优值附近,MACA+MRCR值达到最大,说明此时算法的数据处理效果最佳。而随着μcaptdistsca值的增大,MACA+MRCR值的变化呈现为先几乎保持不变,随后下降,最后又保持不变的特点。

3 }/ B1 X& g" n' x/ H, z4 D

表 4 CUBE、曲面滤波参数敏感性测试范围与步长Tab. 4 Range and step size of CUBE, surface filter parameter sensitivity test

! g5 i% A& z% W5 G( A/ z2 w
+ O, Q% K* U* @$ i
参数最优参数参数范围测试步长μcaptdistsca0.50.1-200.25μcapdistmin0.710-20.1μhes1.960.05-50.25hestimoff40.1-100.5μfilter(×STD)2.50.5-50.5
' y4 K* Z( s- B' _; n
5 h6 d& e5 J: d; P, B3 d

图 6 CUBE、曲面滤波参数敏感性分析结果Fig. 6 CUBE and surface filter parameter sensitivity analysis results

5 f0 d2 X x8 o0 r5 l# S6 b

3.4 优选参数应用

, {! u7 I7 k3 S/ D; |$ `: ?

如图 7所示,应用优选后的CUBE参数,生成整个研究区的CUBE曲面。虽然原始数据集内含有大量异常值,但算法对绝大部分区域的海底地形构建非常准确。同时也可看到个别区域存在明显的测深异常,经分析是由于连续脉冲状异常值导致数据的“信噪比”较低,需进行人工干预编辑(图 7(a))。

' W8 ~# ~ J0 y3 N" W0 D: [$ B

图 7 最优CUBE水深曲面及辅助信息Fig. 7 Optimal CUBE surface and auxiliary information

8 B: c: S! R: L1 x8 D9 |9 b' v

从图 7(b)可看出大部分地形平坦区只有一个水深假设,而在多尺度沙波发育区,尤其是在沙波波峰处存在多个水深假设,且假设的强度较低(图 7(c))。这是由于算法在不确定度传播过程中,假设测深点附近一定范围内的水深值不变,使用零阶预测进行水深值及不确定度传播[34],而沙波发育区水深变化大,这些差异较大的测深点传播到同一格网点内导致多个水深假设的建立。

. }$ Q% A+ E; r2 x

测区中部存在垂直于航迹线的水深异常区,且按一定间隔出现,经分析为ADCP对多波束测深仪的干扰导致(图 7(b))。由于这些异常点较少,通过消歧算法可有效抑制(图 7(a))。沿航迹线左右舷对称分布,且越远离中央波束越严重的波状假地形为恶劣海况导致安装杆抖动引起的(图 7(b)),此异常目前还没有较好的解决办法[15]。

Q' R( \" X1 B( A3 J

以传统人工编辑生成的水深曲面为基准,对最终水深曲面进行精度评估。如图 8(a)、(c)所示,除由于安装杆抖动导致个别区域的差值较大(10~20 cm)外,大部分区域的水深差值小于3 cm;统计显示全区97.2%和97.8%的区域,两者差值在平均水深的1‰和2‰以内(图 8(b))。表明处理得到的水深曲面与基准水深曲面几乎一致,完全满足处理精度要求。

6 C: j1 E3 P$ _% U) }

图 8 CUBE曲面与手工编辑曲面差异Fig. 8 Differences between CUBE surface and manually edited surface

8 @" Q) F% n6 v" S

如图 9所示,对比统计了最优CUBE曲面法与传统人工编辑法处理该数据集所需时间。饼状图面积与总时间成正比,最优CUBE曲面法总用时109 min,且绝大部分时间用于参数联合优选(61.5%)和CUBE曲面编辑(27.5%);人工编辑法总用时876 min,绝大部分时间用于条带编辑(63.6%)和子区编辑(33.0%)。表明最优CUBE曲面法的处理效率远高于传统人工编辑法。

8 N/ T( r* `' I! K! F4 O: y0 `* ^

图 9 最优CUBE曲面法与人工编辑法时效对比Fig. 9 Comparison between optimal CUBE surface method and manual editing method

, x* V5 N1 l9 R+ I: @

综上所述,将优选后的CUBE、曲面滤波参数应用到数据处理中,不但可保证最终的数据处理效果,且与传统人工处理法相比提升多波束数据处理效率达到8倍以上。

3 b6 a/ D; F4 h: _5 D

4 结论

4 \ p- |* m& I- K4 d% P

(1) 系统阐述了CUBE算法的基本原理、数学模型、关键参数和主要处理步骤,通过测深点不确定度运算、不确定度传播、格网点吸收、水深假设建立与更新等步骤完成格网点水深及其精度的联合运算。

% Q& a) z8 Z2 ~0 \& {7 D

(2) 建立了CUBE、曲面滤波参数联合优选方法。通过选取典型地形区、参数试验、对比分析等关键步骤完成参数的联合优选,并用台湾浅滩实测数据对优选方法进行了验证,对处理过程、关键问题和时效对比进行了详细分析。试验表明,优选的参数可大幅提升多波束数据自动处理的精度和效率,具有实际的工程应用价值。

4 |7 B; v# `* v6 R4 P& V, @% b/ \

【引文格式】赵荻能, 吴自银, 李家彪, 等. CUBE曲面滤波参数联合优选关键技术及应用. 测绘学报,2019,48(2):245-255. DOI:

10.11947/j.AGCS.2019.20180082" P. h: K2 u* M5 o 0 j% P ?2 H9 _/ D9 B1 M

- K* K% {8 B K* G) s' @

, g$ c* b% z; X+ {+ q* E

& q8 D: o& X/ M4 C& w8 U/ S

$ b0 ?, C- [8 H6 `6 ~* j" C3 ?3 I. t

《测绘学报》2019年第2期网刊发布

3 I7 T0 v) G7 w. e6 j4 `* c

2018年中国高校发表SCI论文综合排名报告

+ L- d7 A- R( S

来了!普遍上涨,武大、北大等20所知名高校公布2019年考研复试线

- R( G! [* s$ V; h8 e

397名2019年两院院士拟推荐候选人公示名单汇总!

6 b: f! P [& Z

2019年QS世界大学学科排名出炉

* p) T& V+ l& u

杨必胜、张小红、赵齐乐等测绘信息领域专家入选第四批国家“万人计划”入选人员

s1 r' r! |9 m& f/ M/ \4 o

关于召开“测绘前沿科技大讲堂与科技期刊论文写作理论与方法高级研修班”的(一号)通知

* D* } Y0 H/ _

权威 | 专业 | 学术 | 前沿

, `6 n% C T, T _# J$ ~& y4 Z

微信投稿邮箱 | song_qi_fan@163.com

2 V6 t1 O8 ~1 | ]8 C

微信公众号中搜索「测绘学报」,关注我们,长按上图二维码,关注学术前沿动态。

$ \6 n( p& ]; A/ K1 K

欢迎加入《测绘学报》作者QQ群: 297834524

$ x% I8 C2 H. O3 F8 x$ X! |+ z

进群请备注:姓名+单位+稿件编号

* v8 Z" Y) ?3 ^, ]& ` P+ F 3 u2 o2 Q6 Y2 r- u' t. O# C 8 R b. S' S6 @+ i7 `( ? ( \9 W/ W& Z9 D- l2 |# n " O! r: P* ]) [
回复

举报 使用道具

全部回帖
暂无回帖,快来参与回复吧
懒得打字?点击右侧快捷回复 【吾爱海洋论坛发文有奖】
您需要登录后才可以回帖 登录 | 立即注册
文星雨
活跃在3 天前
快速回复 返回顶部 返回列表