作者简介:王灿进(1987-),男,助理研究员,博士.研究方向:激光主动成像.E-mail:wcjpsh@126.com
为研究激光主动成像中散斑噪声的抑制问题,提出一种基于信号子空间TDC(Time-domain constrained)的散斑去噪方法,并搭建一套基于距离选通ICCD的激光主动照明系统进行实验验证。首先使用同态变换将乘性噪声变为加性噪声,然后利用小波变换估计噪声的协方差;接着对含噪图像进行奇异值分解并估计信号子空间的维数,根据该维数对无噪图像的协方差矩阵进行特征值分解,计算出滤波估计矩阵。将滤波估计矩阵与含噪图像卷积,最后做同态逆变换,得到降噪后的图像。结果证明本文的去噪方法拥有比经典的Lee、Frost和Kuan算法更好的散斑噪声抑制效果,同时计算时间明显缩短。
To investigate the suppression of the speckle noise in laser active imaging, a denoising method based on Time-Domain Constrained (TDC) in signal subspace is proposed and a laser active imaging system based on range-gating ICCD is constructed for experiment. First, homomorphic transformation is performed to convert the multiplicative noise to additive noise. Second, wavelet transformation is performed to estimate the covariance of speckle noise. Third, the noise image is decomposed into signal subspace and noisy subspace, and singular value decomposing is used to estimate the dimension of the signal subspace. The covariance of clean image is decomposed using eigenvalue decomposing, and denoising estimating matrix is computed. Fourth, the denoising estimating matrix is convolved with noisy image. Finally, the inverse homomorphic transform is carried out to get the denoised image. Experiment results indicate that, compared with classical Lee, Frost and Kuan filtering, the proposed method has advanced denoising performance and consts less computation time.
激光主动成像[1]系统具有分辨率高、抗干扰能力强、不受环境光限制等优点, 适用于微光、夜视以及远距离暗目标的探测成像, 目前已广泛应用于远距离监视、自动目标识别等领域[2, 3]。在主动成像过程中, 由于相干光经过目标表面反射或者通过无规则涨落的折射媒介传播, 会形成一种信号相关的乘性噪声, 即散斑噪声, 造成图像质量严重下降。文献[4]通过计算灰度概率分布曲线证明散斑噪声是激光主动成像的主要噪声, 因此研究抑制散斑噪声的算法是大多数激光主动成像系统都必须考虑的内容。
国内外学者提出诸多抑制散斑噪声的经典算法, 包括均值滤波[5]、Lee滤波[6]、Frost滤波[7]、Kuan滤波[8]及其改进算法[9]。这些滤波器使用一定大小的滤波窗口在图像上滑动, 使用不同算法计算出当前像素的替代值, 其中Lee滤波由于在滤除噪声和保留信号完整性之间能够达到较好的平衡而倍受关注。这类窗口算法虽然在同类区域去除散斑噪声效果较好, 但运算时间较长, 对于图像边缘和细节信息的保留能力也不理想。变换域降噪算法中使用较多的是小波软阈值滤波[10, 11], 这类算法缺陷在于阈值的选取较为困难, 阈值过大容易丢失图像细节信息, 阈值过小则不能很好地去除噪声。近年来出现了不少针对散斑噪声的新算法, 如Jarabo-Amores等[12]提出一种基于“ mean-shift” 滤波的去噪方法, 该方法不必事先假设噪声的统计模型而能够达到与统计去噪相当的效果。Hyenkyun等[13]提出基于交互技术(Shifting technology)的降噪方法, 能够快速有效地去除散斑噪声。国内, 文献[14]使用SURE(Stein's unbiased risk estimate)获得近似最优的小波系数的权值, 该方法优势在于计算量很小。
信号子空间方法的基本思想是将原始信号通过K-L变换分解为信号子空间和噪声子空间, 通过滤除噪声子空间噪声分量的同时保留信号子空间的信号部分, 能够恢复出较为纯净的信号, 该方法已经成功运用于语音和雷达信号降噪领域[15, 16]。基于其优异的信噪分离效果, 本文将该方法从一维推广到二维空间, 提出一种基于时域约束(TDC)的激光主动成像散斑降噪方法。将本文的TDC降噪方法同经典的Lee、Frost和Kuan滤波作比较, 结果表明本文算法具有出色的噪声抑制效果, 同时算法也具有很好的实时性。
被激光照射的目标表面相对于激光波长而言凹凸不平, 相干光被不同位置反射后在像面处发生干涉, 即形成散斑图样。散斑噪声的数学模型如式(1)(2)所示:
式中:I表示含噪图像; I0表示无噪图像; Isn表示散斑噪声作用系数, 服从均值为1, 方差为1/L的伽马分布, L是常数。
由式(1)可以看出散斑噪声是一种乘性噪声。为了将信号和噪声分离, 使用同态变换将其变为加性噪声:
记上式为:
根据文献[17], 同态变换后散斑噪声的均值和方差的关系分别如式(5)(6)所示:
信号子空间技术是一种类似于傅里叶变换[18]、小波变换[19]的子空间变换技术。其原理是:带噪信号的空间由信号子空间与噪声子空间构成, 噪声能量均匀分布于整个空间中, 而信号能量则只分布于信号子空间中。将带噪信号分解至信号子空间和噪声子空间, 通过消去噪声子空间的分量, 就可以恢复出近似纯净的信号。该方法对于语音信号具有良好的去相关性, 因而广泛用于语音增强领域。本文接下来将这种方法推广到二维图像空间。
设无噪图像的估计为:
则误差信号表示为:
定义失真信号能量和残余噪声能量分别为:
式中:tr( )表示取对角线上的元素。
为使失真信号能量达到最小, 设置约束条件为:
式中:σ 2是常数。对上述约束最小值问题, 可以采用Kuhn-Tucker条件得到Lagrange多项式:
根据条件
可以计算出滤波矩阵H的估计值:
式中:RX=E(XXT)是无噪图像的协方差矩阵, RN=E(NNT)是噪声的协方差矩阵, μ 是Lagrange算子。将RX进行特征值分解, 表示成信号子空间和噪声子空间两部分:
式中:U1、U2分别是m× r和m× (n-r)的特征值矩阵, r是信号子空间的维数, n-r是噪声子空间的维数;
取RX的信号子空间的分量代入式(15), 可得滤波估计矩阵:
由式(16)(17)可知, 为了得到HTDC, 需要确定以下四个参数:信号子空间的维数r、无噪图像的协方差矩阵RX、噪声的协方差矩阵RN以及常数μ 。
估计信号子空间的维数, 可以用奇异值分解实现。信号Y的奇异值分解表示如下:
式中:S是Y的奇异值矩阵, Sm× m=diag(s1, s2, …, sm), s1≥ s2≥ …≥ sm, 即对角阵上的值采用降序排列。
信号某维度的奇异值大小代表该维度的重要程度, 奇异值越大, 对应该维度越重要。取S的前r个分量, 则Y可以近似表示为:
文中选取r的条件为:
式中:Tr是一个接近于1的阈值, 具体大小由实验决定。
由于像素值为0的像素点经过log变换会出现无穷小, 使Y无法进行奇异值分解, 因此同态变换前先将每个像素亮度增加10, 待取指数恢复之后再减去10。
经过同态变换之后乘性噪声变为加性噪声, 如式(4)所示, 则可以得到
式中:RY=E(YYT)可以直接由含噪图像求得。对于RN, 式(6)给出了它和常数L的关系, 但是L难以估计。本文从小波分析的角度出发:在小波域上, 图像信号集中在低频分量, 而噪声信号均匀分布在整个小波空间[16], 因此可以利用高频子带的小波系数估计噪声的强度:
RN=k· (
式中:
μ 的取值会影响滤波后图像的质量。μ 值大, 可以较大程度地滤除散斑噪声, 却会导致较严重的信号失真; μ 值小, 能够保留较多的细节信号, 但会残余更多的噪声。选择μ 的值需要在信号完整性和滤波性能之间权衡。
根据文献[16], 采用以下公式确定μ :
式中:SNR是图像的信噪比; μ 0和s是常数, 由实验确定。
式(23)表示的是:对于高信噪比的图像, 噪声较小, μ 取较小的值, 以保护信号完整性为主要目标; 对于低信噪比的图像, 噪声较大, μ 取较大的值, 以滤波为主要目标。
本文基于信号子空间的散斑去噪方法如下:
(1)对增加亮度后的图像进行同态变换, 将乘性噪声变为加性噪声。
(2)将图像进行小波变换, 并利用高频分量估计噪声方差RN。
(3)对Y进行奇异值分解, 估计信号子空间的维数r。
(4)计算无噪图像的协方差矩阵RX=RY-RN, 并根据式(16)对其进行特征值分解。
(5)根据式(17)计算变换矩阵HTDC, 得到无噪图像的估计值
(6)同态逆变换, 得到最终的去噪图像。
实验搭建一套基于距离选通技术的激光主动照明系统:使用532 nm波长、5 mrad发散角的半导体泵浦固体脉冲激光器发射脉冲光束照射3.4 km处的目标, 用选通ICCD接收反射信号, 通过同步控制器控制选通快门的开闭。其中选通ICCD由选通像增强器和CCD组成, 具有高增益、高灵敏度、成像清晰等优点。同步控制器控制ICCD在纳秒级的时间内进行选通, 可以抑制杂散光, 保证采集到较为清晰的含散斑图像。图1是实验系统框图, 图2是实验系统实物图。
截取ICCD输出图像中被照亮的区域, 最终得到的实验图像大小为150× 150像素。分别采用Lee、Frost、Kuan和本文的TDC滤波方法进行实验, Lee、Frost和Kuan的滤波窗口均选为5× 5。窗口大小选为5× 5可以在滤除噪声和保留图像细节之间得到较好的平衡。实验编程环境为:Windows XP、i7-2600 3.40 GHz CPU、2G内存、VS2008+opencv软件平台。
图3为4种算法的滤波效果对比, ICCD的积分时间为50 ns。
从图3直观上看, Lee、Frost和Kuan结果图中残余较多的散斑噪声, 而TDC滤波后的图像最清晰, 滤波效果最好。观察图像下方中间的门框, 可以发现TDC滤波后的边缘更加清晰, 这是因为其余3种算法均使用邻域像素进行加权操作, 必将导致细节模糊。而TDC算法将图像矩阵当做整体进行运算, 是对图像信号的稀疏表示, 信号子空间的信号部分没有损失, 故边缘保持能力更强。
本文引入峰值信噪比PSNR[20]以评价滤波算法的有效性。改变ICCD的积分时间, 得到不同噪声强度的图像, 分别用4种滤波方法进行实验, 计算其PSNR值, 比较滤波效果, 结果如图4所示。
由图4看出, 本文TDC算法对于不同强度的噪声能够取得比经典的Lee、Frost和Kuan算法更好的滤除效果。观察PSNR曲线, 随着曝光时间的增加, TDC滤波算法的曲线与其余曲线的差值越来越大, 这表示随着噪声的增大, TDC的滤波性能优势相对于其余算法越来越突出。这是因为对于Lee、Frost和Kuan等窗口滤波而言, 使用有限大小的滤波窗口信息, 并不能很好地估计噪声的方差, 随着噪声的增加, 滤波能力并不会显著提升, 滤波后残余的噪声会越来越大。而TDC算法的约束条件是将残余噪声限制在σ 2范围内, 算法在整个图像上进行噪声强度的估计, 并且只保留信号子空间中的非常小的噪声能量, 因此滤波效果更好。
表1统计的是不同算法对于150× 150图像的平均运算时间。可见本文算法比其余3种算法耗时要小得多。这是因为TDC算法对整个图像矩阵同时进行计算, 而其余3种算法均需要滑动滤波窗口遍历整个图像区域。在本文的实验条件下(成像区域大小为150× 150, 帧频为25 帧/s), TDC滤波算法能够满足实时处理的要求。
受语音增强领域的子空间分解方法的启发, 本文在分析散斑噪声数学模型的基础上, 提出一种基于信号子空间的激光主动成像去噪方法, 并搭建了一套基于距离选通技术的激光主动照明系统进行实验验证。本文算法将含噪图像分解为信号子空间和噪声子空间, 通过去除噪声子空间并尽量保留信号子空间的能量, 以取得良好的去噪效果。结果证明本文TDC方法具有比经典的Lee、Frost和Kuan算法更好的去噪性能, 同时缩短了计算时间, 满足激光主动成像系统实时性的要求。
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|