基于摆动式单波束微地形探测的高程异常值剔除算法研究.doc

第7期 基于摆动式单波束微地形探测的高程异常值剔除算法研究 · 21 ·
更多电子资料请登录赛微电子网
基于摆动式单波束微地形探测的高程异常
值剔除算法研究*
罗柏文1 夏毅敏2 卜英勇2 周知进1
(1. 湖南科技大学机电工程学院, 湘潭 411201; 2. 中南大学机电工程学院, 长沙 410083)
摘 要: 针对摆动式单波束探测水下微地形出现的高程异常值, 提出一种剔除异常值的新算法。首先根据摆动式单波束的探测原理以及与多波束探测的等价性采用数据加窗法, 其中Z方向窗口尺寸通过3s 准则确定, X、Y方向窗口尺寸依据地形变化坡度小于45o的经验准则。然后利用反射回波与散射回波在能量上的差异确定起始加窗点; 加窗点采用实测点, 并要求起始加窗点尽可能靠近一条测线上实测点的高程均值。在摆动式单波束探测装置探测随机微地形中, 其结果表明该算法能成功剔除高程异常值。
关键词: 单波束；微地形；高程异常值；加窗；能量
中图分类号: TB559　　　文献标识码: A　　　国家标准学科分类代码: 420.5050
Study on algorithm of eliminating exceptional elevation value based on
micro-terrain detection with pendulum single beam
Luo Bowen1 Xia Yimin2 Bu Yingyong2 Zhou Zhijin1
(1. College of Mechanical and Electrical Engineering, Hunan University of Science and Technology, Xiangtan 411201, China;
2. College of Mechanical and Electrical Engineering, Central South University, Changsha 410083, China)
Abstract: Aiming at the pendulum single beam bathymeter used to explore underwater micro-terrain appears exceptional elevation value, a new algorithm is proposed. First of all, the data windowing method is used under single-beam pendulum detection principle and multi-beam detection of equivalence, in which the window size of Z direction is determined through the 3s criteria, the window size of X and Y direction is given based on the experience criteria of terrain slope less than 45º. Then, using the difference between reflection echo and scattering echo in the energy the starting windowed point is determined. Actual detection data are used as the windowed point, and near mean value of actual detected points as far as possible is used as the starting windowed point. The detection results of random micro-terrain in the pendulum single beam bathymeter show that the algorithm can successfully remove the exceptional elevation value.
Keywords: single beam; micro-terrain; exceptional elevation value; windowed; energy

1 引 言
在研究深海钴结壳微地形探测方法中, 采用摆动式单波束探测和振幅相关检测法进行信号处理, 能够达到钴结壳微地形探测实时性要求。但是超声信号因掠射角、底质特性、外界噪声等干扰因素的影响, 其反射回波会远离水声换能器的接收位置, 只能接收到探测点的散射回波。当探测点的散射回波隐藏在干扰因素造成的回波中时, 利用振幅相关法检测到的最大相关值有可能不是探测点真实回波产生的, 这样不可避免地导致高程异常值的出现[1-4]。高程异常值是造成严重“伪”地形的首要因素, 因此有必要对其剔除。
吴英姿博士处理H/HCS-017型条带测深仪获得的水深数据, 依据地形变化的坡度小于45º的经验准则[5]。摆动式单波束探测与多波束探测的工作原理不同, 但它通过水声换能器的摆动使之达到横向大扇区、全覆盖的窄带探测, 并且水声换能器本身的高指向性等价于多波束的预成波束, 因此该经验准则也适合于摆动式单波束探测[2]。本文依据该准则, 根据摆动式单波束探测原理和振幅相关检测法工作特点, 提出了一种新的高程异常值剔除算法。由实验证明, 该算法能成功剔除高程异常值。
2 摆动式单波束探测简述
图1是摆动式单波束探测装置的XOZ截面示意图。由收发一体水声换能器发射出频率为500 kHz、波束角b =1.5º的超声波照射到水底面, 当碰到底面时产生反射或散射回波信号。信号接收系统在超声波发射的同时进入工作状态, 接收发射信号和回波信号。当接收完一个探测点的回波信号后, 由步进电机驱动水声换能器以1º的等间隔绕摆动中心转动角度, 进行下一次探测。并在水声换能器工作的同时, 探测装置一直沿着导轨前进。最后由发射信号与回波信号进行振幅相关性分析获得探测点高程数据。

图1 摆动式单波束探测装置示意图
Fig.1 Sketch map of pendulum single beam bathymeter

振幅相关检测算法如下[6-8]:
假设发射信号时序函数为, 回波信号时序函数为, D为采样间隔, h为发射信号采样点数, x 为回波信号采样点数。根据广义平稳随机序列的互相关函数可得
(1)
式中: s为时延数, 当s=DT时, p(t)与q(t)具有最大相关性。那么p(t)与q(t)的互相关函数无偏估计表达式为
(2)
式中: 。设采样容量N=h, 时延数, , 则式(2)离散序列表达式为
(3)
式中: , 由式(3)获取互相关函数最大值所对应的信号延时来确定DT。再由DT计算出水声换能器到探测点的距离:
(4)
式中: V为水中声速。最后在XOZ截面内, 以摆动中心为坐标原点, 第k个探测点对应的X坐标和Z坐标分别为
(5)
式中: , H为水声换能器板面到摆动中心的距离, q 为水声换能器绕摆动中心摆动的角度。假设微地形高程参考值为, 则第k个探测点的高程值为
(6)
另外, 在不增加任何硬件情况下, 利用振幅相关检测法计算高程值的同时, 可获得相应探测点的回波能量[7-8]:
(7)
式中: g就是p(t)与q(t)具有最大相关性时s的时延数。
3 高程异常值剔除算法
3.1 高程参考值确定
1) 将X方向测点的探测距离代入式(5), 计算出第j 条测线上测点的X、Z坐标值。设摆动中心为坐标原点, 则zi,j都为负值。i表示X方向的测点序列, j表示Y方向的测点序列。
2) Z坐标值的均值mz与均方差sz:
(8)
(9)
理论上应采用Min(zi,j)作为,但Min(zi,j)有可能是异常值。根据3s 准则, 的概率为95.45%。当, 则取Min(zi,j); 否则, 取。
3.2 X方向数据窗尺寸确定
设同一测线上两相邻测点在Z方向的坐标值为zi,j、zi+1,j, 其对应的摆角为qi、qi+1, 依据地形变化的坡度小于45º的经验准则, zi,j、zi+1,j应满足:
(10)
式中: m为冗余量系数, 则
Max(m)
m (11)
令 XCH=m, 则设定X方向数据窗宽为XCH, 窗高为2XCH。
3.3 Y方向数据窗尺寸确定
Y方向代表摆动式单波束探测装置的行进方向。假设行进速度为v, 水声换能器摆动周期为T, 沿Y方向的相邻探测点的坐标值zi,j+1、zi,j应满足:
(12)
式中: m为冗余量系数。令YCH=, 则设定Y方向数据窗宽为YCH。
3.4 回波能量阈值确定
由于导致高程异常值产生以掠射角、底质特性、外界噪声为主, 一条测线所对应的被测地形, 因其范围极小, 其底质特性可认为是均质; 而掠射角受微地形随机起伏影响, 无规律可循, 即使水声换能器摆角为0º 时探测的高程值也会出现异常, 由此在数据加窗判别高程异常值中恒定选取正常的高程值作为起始加窗点是一个难题。
一般来说, 如果水声换能器接收到探测点的回波是反射回波, 则认为其高程值是正常的; 而散射回波容易受外界因素干扰, 需要借助其他方法判断其高程值的正确性。故利用回波能量这一特点参与到数据加窗中实现起始加窗点的选取工作。但是反射回波与散射回波没有明确的界限。为提高数据加窗判别的准确性, 计算每条测线上的回波能量的平均值mW:
(13)
式中: n为测线上的测点数, W为测点的回波能量值; 将mW设为阈值, 认为回波能量低于该阈值时, 其对应测点的高程值可能会出现异常。
3.5 异常值判别步骤
1) 将坐标值zi,j排序, 以接近均值mz为原则, 由此将最接近均值mz的坐标值作为起始加窗点;
2) 引入回波能量, 若起始加窗点的回波能量不低于平均值mW, 则以该点为中心向其测线两边的坐标值依次加窗判别; 反之, 则将仅次于它接近均值mz的坐标值作为起始加窗点, 重新引入回波能量进行判别, 依次类推;
3) 以加窗点zi,j为窗口中心, 当zi-1,j与zi+1,j都不在窗口内, 则认为zi,j为异常值, 并剔除; 当zi-1,j或zi+1,j在窗口内, 则借助Y方向数据窗口进行判断, 若不满足式(12)则认为zi,j为异常值, 并剔除。其他情况可认为zi,j是正常值;
4) 加窗处理后的坐标值zi,j由式(6)转化为高 程值。
4 应 用
利用摆动式单波束探测装置对随机微地形进行探测, 表1和表2分别为6条测线上的探测距离和回

表1 随机微地形探测距离/mm
Table 1 Detecting distance of random micro-Terrain/mm
序号
1
2
3
4
5
6
1
1 585.3
1 223.9
1 522.8
1 520.7
1 570.7
1 565.6
2
1 585.3
1 522.3
1 522.7
1 573.0
1 569.3
1 564.1
3
1 577.4
1 515.2
1 515.4
1 509.5
1 572.3
1 568.6
4
1 574.5
1 516.6
1 513.6
1 512.6
1 564.1
1 562.1
5
1 564.7
1 507.5
1 509.1
1 555.6
1 566.7
1 530.1
6
1 564.9
1 507.5
1 514.1
1 552.4
1 553.3
1 539.6
7
1 555.3
1 501.0
1 501.0
1 548.4
1 552.0
1 548.2
8
1 563.2
1 499.4
1 497.9
1 547.3
1 547.5
1 500.4
9
1 565.7
1 499.2
1 499.3
1 550.6
1 552.8
1 498.8
10
1 565.7
1 559.9
1 556.1
1 551.9
1 549.4
1 484.0
11
1 565.2
1 556.9
1 504.2
1 550.3
1 547.9
1 484.0
12
1 563.2
1 506.1
1 506.1
1 504.8
1 544.8
1 486.4
13
1 561.5
1 558.9
1 507.9
1 506.9
1 544.6
1 485.1
14
1 568.5
1 559.4
1 506.5
1 544.3
1 544.8
1 482.9
15
1 565.6
1 558.3
1 555.4
1 551.1
1 554.4
1 487.1
16
1 567.3
1 559.8
1 555.6
1 350.8
1 528.9
1 544.0
17
1 570.7
1 565.4
1 519.7
1 513.8
1 496.1
1 541.1
18
1 574.4
1 568.9
1 567.4
1 510.6
1 496.5
1 523.2
19
1 576.1
1 573.5
1 573.5
1 571.0
1 562.8
1 558.9
20
1 613.8
1 612.3
1 574.9
1 572.5
1 608.3
1 580.2
波能量数据。表3是由上述算法(m=1.2, v=50 mm/s, T=2s)处理的各测线参数值, 处理结果生成图形, 由图2~图7所示。从图中可看出, 高程异常值被成功剔除。

表2 随机微地形回波能量/V2·ms
Table 2 Echo energy of random micro-terrain/V2·ms
序号
1
2
3
4
5
6
1
6.49
4.56
84.35
14.93
86.66
21.18
2
8.75
5.08
83.57
17.30
75.29
15.65
3
33.95
42.25
43.65
6.99
30.71
18.40
4
87.15
59.03
171.22
40.26
104.45
16.84
5
44.56
52.09
54.20
22.59
65.47
13.93
6
24.88
11.11
28.39
9.27
43.46
18.10
7
27.93
366.05
17.99
8.17
89.50
17.87
8
32.29
342.19
137.81
42.88
59.34
52.85
9
97.39
20.84
71.55
47.84
18.47
30.11
10
84.35
30.82
68.66
436.27
180.37
11.59
11
604.70
54.50
78.31
237.28
122.71
30.04
12
666.49
48.09
564.65
45.30
22.33
45.71
13
102.56
21.35
538.18
33.18
57.74
124.06
14
265.98
15.05
30.75
17.01
120.48
15.27
15
189.42
128.58
112.53
30.52
24.29
36.24
16
103.25
24.46
20.17
4.50
20.57
8.05
17
24.45
49.30
18.86
19.50
18.00
36.46
18
23.91
71.42
20.57
28.89
64.45
11.36
19
81.25
40.85
43.08
12.92
13.60
28.45
20
23.40
51.75
24.63
123.57
12.02
13.89

表3 各测线处理后的具体参数
Table 3 Parameter of each disposed detected line
测 线
1
2
3
4
5
6
mW /V2·ms
126.66
71.97
110.66
59.96
61.50
28.30
/mm
1712.3
1723.3
1691.0
1689.2
1719.4
1691.7
XCH /mm
36.9
37.1
36.4
36.4
37.0
36.4
YCH /mm
60.0
60.0
60.0
60.0
60.0
60.0

图2 第1条测线处理前后对比
Fig. 2 Fore-and-aft contrast of the first detected line

图3 第2条测线处理前后对比
Fig. 3 Fore-and-aft contrast of the second detected line

图4 第3条测线处理前后对比
Fig. 4 Fore-and-aft contrast of the third detected line

图5 第4条测线处理前后对比
Fig. 5 Fore-and-aft contrast of the fourth detected line

图6 第5条测线处理前后对比
Fig. 6 Fore-and-aft contrast of the fifth detected line

图7 第6条测线处理前后对比
Fig. 7 Fore-and-aft contrast of the sixth detected line

5 结 论
本文针对摆动式单波束探测出现的高程异常值导致严重“伪”地形进行异常值剔除算法研究。首先从摆动式单波束探测方法和振幅相关检测法的特点分析高程异常值产生的原因, 并根据摆动式单波束探测与多波束探测的等价性提出采用数据加窗法。加窗法中Z方向窗口尺寸通过3s 准则确定, X、Y方向窗口尺寸依据地形变化坡度小于45º的经验准则。然后为了进一步提高异常值剔除的准确性, 利用反射回波与散射回波在能量上的差异确定起始加窗点; 加窗点采用实测点, 并要求起始加窗点尽可能靠近一条测线上实测点的高程均值。最后由随机微地形探测实验证明, 该算法能成功剔除高程异常值。
参考文献:
[1] 赵海鸣,卜英勇,等.摆动式单波束超声波水下微地形探测[J]. 中南大学学报:自然科学版, 2007, 38(5): 932-936.
ZHAO H M, BU Ying-yong, et al. Detection of seabed cobalt-rich crusts micro-topography based on single-beam ultrasonic wave[J]. Journal of Central South University: Natural Sciece Edition, 2007, 38(5): 932- 936.
[2] 刘军.基于分形理论的单探头单波束微地形探测系统研究[D]．长沙: 中南大学, 2006: 11-28.
LIU J. Study of tiny-terrain detection system of single-beam based on fractal theory[D]. Changsha: Central South University, 2006: 11-28.
[3] YAMAZAKI T, SHARMA R. Morphological features of co-rich manganese deposits and their relation to seabed slopes[J]. Marine Georesources & Geotech-nology. 2000, 18(1): 43-76.
[4] 赵海鸣,卜英勇. 一种高精度超声波测距方法的研究[J]. 湖南科技大学学报:自然科学版, 2006, 21(3): 35-38.
ZHAO H M, BU Y Y, et al. Study on the ultrasonic distance measurement method with high precision[J]. Journal of Hunan University of Science & Technology: Natural Science Edition, 2006, 21(3): 35-38.
[5] 吴英姿.多波束测深系统地形跟踪与数据处理技术研究[D]. 哈尔滨: 哈尔滨工程大学, 2001: 96-101.
WU Y Z. Study on technology of multi-beam bathymetric system terrain following and data disposal[D]. Harbin: Harbin Engineering University, 2001: 96-101.
[6] 王纪婵. 基于互相关理论的高精度水下超声波测距系统研究[D]. 长沙: 中南大学, 2006: 24-31.
WANG J CH. Study on underwater ultrasonic distance measurement system with high precision based on cross- correlation theory[D]. Changsha: Central South University, 2006: 24-31.
[7] 罗柏文,卜英勇. 基于能量相关搜索法的海底钴结壳微地形探测方法研究[J]. 矿冶工程, 2007, 1: 21-24.
LUO B W, BU Y Y. Exploration method for tiny terrain of seabed cobalt crust based on energy relativity search[J]. Mining and Metallurgical Engineering, 2007, 1: 21-24.
[8] 张晓东, 高波. 互相关函数法在声学测温技术中的应用研究[J]. 中国电机工程学报, 2003, 23(4) : 185-188.
ZHANG X D, GAO B. The research of acoustic measuring of gas temperature employing cross-correlation algorithm[J]. Proceedings of the Csee. 2003, 23 (4): 185-188.
作者简介:
罗柏文: 男, 1975年出生, 湖南茶陵县人, 2003年于中南大学获得硕士学位, 2008年于中南大学获得博士学位, 现为湖南科技大学讲师, 主要研究方向为深海采矿装备技术及理论。
Email: luobow@
Luo Bowen: male, born in 1975. He received MS and PhD both from Central South University in 2003 and 2008, respectively. Now he is a lecturer at Hunan University of Science and Technology. His main research directions are the deep-sea mining equipment technology and theory.

第3期 汤清虎 等: 非晶态Mn-Ce-O催化芒香醇选择氧化 23

安捷伦科技推出突破1瓦输出功率大关的微波信号发生器

安捷伦科技公司日前宣布,推出业内首款具有超高输出功率、在业界率先突破1瓦输出功率大关的微波信号发生器。安捷伦PSG E8257D 信号发生器选件 521 可在 250 MHz ~20 GHz(低至 10 MHz 仍可使用)的频率范围内提供 +24 至 +28 dBm 的额定输出功率, 无需再使用放大器、耦合器和检波器等其他 硬件。
安捷伦科技副总裁兼信号源业务部总经理Andy Botka表示：“继推出业内首款频率高达 44 GHz 的合成矢量信号发生器和首款频率高达 67 GHz 的模拟信号发生器之后, 安捷伦现在又推出了新型微波信号发生器。它秉承了安捷伦在射频和微波测试应用方面加速创新的传统, 可满足航空航天、国防和通信行业的需求。”
几乎所有的微波信号发生器都有一个共同的局限性, 那就是额定输出功率通常都小于 +25 dBm。对于那些需要 +25 dBm 以上的额定输出功率的应用, 就需要使用放大器、耦合器和检波器等外置硬件, 例如 TWT 放大器测试和 ATE/天线测试配置(由于开关路由和电缆较长, 它们的信号插入损耗通常很高)。而带选件 521的PSG 无需使用外置硬件, 因 此降低了用户的外置硬件成本。
该信号发生器在大部分频率范围内典型的最大输出功率均超过1瓦(大于+30 dBm), 信号发生器内部可借助新型的射频输出功率可调保护电路, 避免功率过高的问题, 从而保护用户的被测设备(UUT)。这个可调功率保护电路可通过内置和外置的电平配置, 在 PSG 输出端设置为 +15 dBm 至 +33 dBm。 对于高功率、窄频率范围的应用, +33 dBm 的最大射频输出保护电路配置提供了获得选件 521 的最大输出功率的能力。

关于安捷伦科技公司

安捷伦科技是全球领先的测量公司, 是通信、电子、生命科学和化学分析领域的技术领导者, 公司的19,000名员工在110多个国家为客户服务。在2008财政年度, 安捷伦的业务净收入为58亿美元。要了解安捷伦科技的信息, 请访问：www.agilent. 。