劉效勇,盧佩*,曹海賓,田敏
(1石河子大學(xué)理學(xué)院,新疆 石河子832003;2石河子大學(xué)信息科學(xué)與技術(shù)學(xué)院,新疆 石河子832003)
基于壓縮感知的圖像重構(gòu)關(guān)鍵技術(shù)研究
劉效勇1,盧佩2*,曹海賓1,田敏2
(1石河子大學(xué)理學(xué)院,新疆 石河子832003;2石河子大學(xué)信息科學(xué)與技術(shù)學(xué)院,新疆 石河子832003)
圖像在信息領(lǐng)域應(yīng)用的重要性使人們對快速成像技術(shù)提出了更高需求,壓縮感知理論為圖像處理提供了一種新途徑。本文研究了壓縮感知理論應(yīng)用于圖像重構(gòu)的關(guān)鍵技術(shù),重點(diǎn)分析對比了各種因素對重構(gòu)圖像質(zhì)量的影響。實(shí)驗(yàn)結(jié)果表明:通過小波變換基對信號稀疏,以哈達(dá)瑪矩陣作為測量矩陣,使用正交匹配追蹤算法時(shí)僅需50%的采樣率即可獲得較高質(zhì)量的圖像,有效減少了傳統(tǒng)方法中冗余數(shù)據(jù)過多問題,在重構(gòu)圖像視覺效果和PSNR值上均有所提高,同時(shí)提高了恢復(fù)效率。本研究為應(yīng)用壓縮感知尋求最優(yōu)化圖像恢復(fù)方法在理論和技術(shù)上提供了有力依據(jù)和支撐。關(guān)鍵詞:壓縮感知;圖像重構(gòu);稀疏信號;測量矩陣;重構(gòu)算法
現(xiàn)實(shí)世界的虛擬仿真以及對處理信號數(shù)字化,決定了從模擬信號源獲取數(shù)字信息時(shí)必須通過信號采樣。但是在數(shù)字圖像以及語音視頻的獲取過程中,按照傳統(tǒng)的奈奎斯特采樣定理進(jìn)行采樣容易導(dǎo)致采樣數(shù)據(jù)過多;另外,隨著人們對信息需求數(shù)量的不斷增長,要求信息載體的帶寬愈來愈寬,采集速率以及信號處理速度也會愈來愈高,這給信號的采集、處理、傳輸以及存儲等帶來了很大的挑戰(zhàn)。
學(xué)者們[1-3]提出了一種新穎的信息獲取理論,即壓縮感知(compressed sensing,CS)[1-3],該理論指出:當(dāng)信號具有稀疏性或可壓縮性時(shí),采集少量數(shù)據(jù)就可精確或近似精確地重構(gòu)原始信號。壓縮感知突破了傳統(tǒng)的信號采樣瓶頸,能夠以遠(yuǎn)低于奈奎斯特采樣頻率采樣和恢復(fù)信號,因而可以更加有效存儲和傳輸信息。該理論一經(jīng)提出就受到國內(nèi)外科研人員的高度關(guān)注。目前,壓縮感知技術(shù)已被廣泛應(yīng)用在計(jì)算機(jī)圖形學(xué)、天文學(xué)、醫(yī)學(xué)成像、高光譜成像、遠(yuǎn)距離遙感、模擬-信息轉(zhuǎn)換、地球物理學(xué)數(shù)據(jù)分析、生物傳感、分布式壓縮感知等領(lǐng)域中[4-6]。此外,多種壓縮感知硬件相繼出現(xiàn),如美國麻省理工大學(xué)L.Wald等人研制的MRJ RF脈沖設(shè)備,W.Freeman等人研制的編碼孔徑相機(jī)、美國萊斯大學(xué)研制的單像素照相機(jī)、美國伊利諾伊州立大學(xué)Milenkovic等人研制的DNA微陣列傳感器、美國耶魯大學(xué)研制的超光譜成像儀、中科院李廉林等人研制的CS濾波器和混沌腔等[7-11]。這些壓縮感知硬件儀器的實(shí)現(xiàn)使壓縮感知理論向?qū)嵱没斑M(jìn)了一大步。
應(yīng)用壓縮感知,可以有效地提高應(yīng)用系統(tǒng)的性能、縮短數(shù)據(jù)獲取的時(shí)間、降低數(shù)據(jù)獲取的耗能、減少數(shù)據(jù)的存儲空間、提高數(shù)據(jù)的傳輸速度等。壓縮感知與信號的稀疏重構(gòu)尤其是在圖像的壓縮與重建方面有著非常大的應(yīng)用價(jià)值。然而,由于缺乏有效的壓縮傳感矩陣判別理論,重構(gòu)算法需要迭代運(yùn)算,造成復(fù)雜度高速度慢,算法過程不可控,參數(shù)敏感等問題。另外,實(shí)際應(yīng)用時(shí)也容易出現(xiàn)壓縮比太大,恢復(fù)出的圖像效果嚴(yán)重模糊,無法使用。為此,本文從多個(gè)方面對影響壓縮感知信號重構(gòu)的因素做了初步探索研究,以便從中找到壓縮感知重構(gòu)圖像的最優(yōu)化方法,從而提高重構(gòu)信號的效率和質(zhì)量。
壓縮傳感理論主要包括信號的稀疏表示、編碼測量和重構(gòu)算法等三方面[12-15]。
考慮長度為N的一維實(shí)值離散時(shí)間信號x,可看作RN空間N×1維的列向量,元素為x[n],n=1,2,...,如果是圖像或高維數(shù)據(jù)矢量,則轉(zhuǎn)化成一個(gè)長的一維向量。RN空間的任何信號都能用N×1維的正交基的線性組合表示。把向量"作為列向量形成N×N的基矩陣ψ=ψ1,ψ2,…ψN,則任意信號x都可以表示為:
其中α為N×1列向量的加權(quán)系數(shù),T表示轉(zhuǎn)置操作,x和α是同一個(gè)信號的等價(jià)表示,x是信號在時(shí)域的表示,α則是信號在ψ域的表示。如果式(1)中α只有 K(KN)個(gè)非零(或絕對值較大)的系數(shù),x稱為K稀疏信號。信號的可稀疏性表示是壓縮感知的先驗(yàn)條件,在已知信號是可壓縮的前提下,壓縮感知過程可分為兩步:
(1)設(shè)計(jì)1個(gè)與變換基不相關(guān)的M×N(M N)維測量矩陣對信號進(jìn)行觀測,得到M維的測量向量。
(2)由M維的測量向量重構(gòu)信號。
用1個(gè)與變換矩陣不相關(guān)的測量矩陣φ對信號進(jìn)行線性投影,得到線性測量值y:
其中φ是個(gè)M×N的隨機(jī)測量矩陣,=φψ。構(gòu)造合理有效的測量矩陣是壓縮感知研究的核心問題。Candès等[2]指出,高斯隨機(jī)變量形成測量矩陣可以作為普適的CS觀測矩陣。
因?yàn)閥的維數(shù)小于x的維數(shù),所以公式(2)有無窮解,重構(gòu)原始信號比較困難。既然x為K稀疏,可以利用下面的優(yōu)化問題求解x:
當(dāng)M≥cKlog(N)時(shí),φ具有有限等距性(RIP),x[n]能被重構(gòu),此處c是個(gè)小常系數(shù)。
上面的問題可以通過貪婪迭代算法解決,其中最常用的是使用正交匹配追蹤算法(OMP)。
最后,當(dāng)稀疏系數(shù)α被精確重構(gòu)后,通過ψ的逆矩陣重構(gòu)信號x:
壓縮感知指出:當(dāng)信號在某個(gè)變換域是稀疏的或可壓縮的,可以利用與變換矩陣非相干的測量矩陣將變換系數(shù)線性投影為低維觀測向量,同時(shí)這種投影保持了重建信號所需的信息,通過進(jìn)一步求解稀疏最優(yōu)化問題就能夠從低維觀測向量精確地或高概率精確地重建原始高維信號。在該理論框架下,采樣速率不再取決于信號的帶寬,而在很大程度上取決于2個(gè)基本準(zhǔn)則:稀疏性和非相干性,或者稀疏性和等距約束性。
壓縮感知能夠依據(jù)少量采樣數(shù)據(jù),從低維觀測向量精確地重建原始高維信號。該理論包含3個(gè)關(guān)鍵部分:稀疏變換基、測量矩陣和重構(gòu)算法。因此,通過壓縮感知理論對原始信號進(jìn)行采樣與重構(gòu),在重構(gòu)過程中信號自身的稀疏程度,信號采樣率、選用的測量矩陣和重構(gòu)算法等不可避免地影響重構(gòu)信號的質(zhì)量,為此,本研究從以上幾個(gè)重要方面對影響信號壓縮感知信號重構(gòu)的因素做了初步探索研究。為了更能直觀地判斷各算法的恢復(fù)效果,以重構(gòu)圖像和原始圖像的峰值信噪比(PSNR)、信噪比(SNR)、均方差(MSE)值以及運(yùn)算速度(Run-Time,RT)作為評判重構(gòu)圖像效果的客觀評價(jià)標(biāo)準(zhǔn)。
實(shí)驗(yàn)中選用256*256像素的經(jīng)典灰度圖像“l(fā)ena”進(jìn)行測試,測量矩陣使用高斯隨機(jī)矩陣構(gòu)造,重構(gòu)圖像時(shí)采用正交匹配追蹤算法。圖1為采樣率(Sampling Rate,SR)分別為20%至80%時(shí)圖像的重構(gòu)效果圖。由實(shí)驗(yàn)結(jié)果可見,在測量矩陣和恢復(fù)算法相同的情況下,信號的采樣率越高重構(gòu)圖像的質(zhì)量越好。采樣率大于50% 時(shí)恢復(fù)結(jié)果比較清晰;采樣率為40% 時(shí),仍能較好的重構(gòu)圖像;在采樣率為30% 時(shí),恢復(fù)的效果很差,出現(xiàn)了較多噪聲,圖像大部分細(xì)節(jié)信息丟失。這也驗(yàn)證了壓縮感知的測量值與信號稀疏度的關(guān)系,只有當(dāng)測量值大于某個(gè)固定范圍時(shí),即滿足測量值M≥cKlog(N)時(shí),才能較好的構(gòu)建出原始圖像。當(dāng)然,隨著采樣率的增大,采樣的測量值就越多,運(yùn)算所需時(shí)間就越長。
圖1 采樣率不同時(shí)圖像重構(gòu)效果對比Fig.1 Results contrast under different sampling rate
應(yīng)用壓縮感知理論對信號處理性能分析,必須保證該信號在某個(gè)特定正交空間上是稀疏的。由于1個(gè)信號的稀疏性與所選取的正交空間特性密切相關(guān),因而,選用1個(gè)最合適的稀疏變換域?qū)⑵浔硎境鰜矸浅V匾?。為此,此處對同一個(gè)信號數(shù)據(jù)矩陣稀疏變換基分別采用離散余弦變換 (DCT),快速傅里葉變換(FFT)和小波變換(DWT)[16-18],對比重構(gòu)信號的質(zhì)量和性能。在重構(gòu)過程中,測量矩陣使用高斯矩陣,重構(gòu)算法使用正交匹配追蹤算法。
實(shí)驗(yàn)中的測試圖像仍然選用256*256像素的“l(fā)ena”灰度圖像,原始圖像的采樣率(M/N)分別為40%、50%、60%時(shí),各種情況下重構(gòu)質(zhì)量如表1所示。比較分析結(jié)果(表1)可見,相同的采樣率下,采用小波變換時(shí)重構(gòu)圖像質(zhì)量最好,并且運(yùn)行時(shí)間最快。采用離散余弦變換時(shí)恢復(fù)圖像質(zhì)量效果最差,而采用快速傅里葉變換所需時(shí)間最長。因而,對信號進(jìn)行稀疏時(shí)應(yīng)以使用小波變換為宜。
表1 稀疏變換基不同時(shí)圖像重構(gòu)效果對比Tab.1 Results contrast under different sparse transform basis
接表1
測量矩陣是壓縮感知中最重要的部分,在原始信號抽樣的獲取與信號重建過程中起著關(guān)鍵性作用。因此,測量矩陣的選擇和構(gòu)造是壓縮感知的核心內(nèi)容,在測量矩陣設(shè)計(jì)時(shí)應(yīng)首先考慮到2個(gè)基本準(zhǔn)則:具有有限等距性、測量矩陣和稀疏矩陣間具有不相關(guān)性。
測量矩陣之間性能各有差異,實(shí)際使用時(shí)盡量選用性能盡可能好的測量矩陣。此處選擇了目前經(jīng)常使用的已被發(fā)現(xiàn)且具有較好性質(zhì)的高斯矩陣(Gauss Matrix)、哈達(dá)瑪矩陣 (Hadamard Matrix)、伯努利矩陣 (BernouliMatrix)、部分正交矩陣(Partial OrthogonalMatrix)以及諾斯萊特矩陣(Noiselets Matrix)作為測量矩陣時(shí)重構(gòu)信號的效果[19-21],實(shí)驗(yàn)中使用256*256像素的灰度圖像“flower”做為測試圖像,通過小波變換基對信號稀疏,正交匹配追蹤(OMP)算法進(jìn)行信號重構(gòu)。由實(shí)驗(yàn)結(jié)果(圖2)可見,相同采樣率下,使用高斯矩陣、哈達(dá)瑪矩陣和部分正交矩陣作為測量矩陣時(shí)重構(gòu)圖像運(yùn)行時(shí)間較短,而將伯努利矩陣和諾斯萊特矩陣用作測量矩陣時(shí)恢復(fù)圖像耗時(shí)較長;使用諾斯萊特矩陣作為測量矩陣圖像重構(gòu)效果最為理想,使用哈達(dá)瑪矩陣和部分正交矩陣時(shí)恢復(fù)重構(gòu)效果差別不是很大,采樣率較低時(shí),使用伯努利矩陣重構(gòu)效果最不理想。
圖2 測量矩陣不同時(shí)圖像重構(gòu)效果對比Fig.2 Results contrast under different measurement matrix
圖3 為采樣率取50%時(shí)選用不同測量矩陣的重構(gòu)效果圖,通過圖像視覺和PSNR值也可以比較重構(gòu)圖像質(zhì)量優(yōu)劣。
圖3 相同采樣率時(shí)各測量矩陣圖像重構(gòu)效果Fig.3 Results contrast under different measurement matrix
在相同條件下,壓縮感知不同重構(gòu)算法有不同特點(diǎn),應(yīng)用時(shí)盡量采用效果最優(yōu)化算法,此處對3種常用的基追蹤算法(BP)算法、正交匹配追蹤(OMP)算法和階段正交匹配追蹤(StOMP)算法做了對比研究[22-25]。實(shí)驗(yàn)中,對256*256的灰度圖像“Barbara”進(jìn)行重構(gòu),使用小波對圖像進(jìn)行變換,使用效果比較好速度較快的結(jié)構(gòu)化隨機(jī)矩陣作為測量矩陣。
(1)相同采樣率時(shí)各算法的性能比較。實(shí)驗(yàn)中設(shè)定圖像的采樣率(M/N)恒為80%時(shí),得到的對比結(jié)果如圖4所示。由圖4可見,3種算法采樣率相同時(shí),BP算法重構(gòu)圖像效果最好,OMP重構(gòu)圖像效果較好,而StOMP算法重構(gòu)圖像效果最差。3種算法中,BP算法運(yùn)行時(shí)間最長,StOMP算法運(yùn)行時(shí)間最短。
圖4 相同采樣率不同算法時(shí)的性能比較Fig.4 Results contrast under different algorithms
(2)采用率不同時(shí)各算法的性能比較。為了進(jìn)一步比較各算法的性能,對原始圖像在不同采樣率情況下重構(gòu)結(jié)果進(jìn)行了比較,原始圖像的采樣率(M/N) 分別為 80%、70%、60%、50%、40%、30%時(shí),各算法對應(yīng)的重構(gòu)圖像效果對比如圖5所示。
圖5 不同算法時(shí)圖像重構(gòu)效果Fig.5 Results contrast under different algorithms
由圖5可見,3種算法均是隨著采樣率的增加,重構(gòu)圖像與原始圖像的峰值信噪比PSNR增大,如圖5b。說明信號的采樣率越高,重構(gòu)圖像質(zhì)量越高,效果越明顯。相同采樣率下,BP算法重構(gòu)圖像質(zhì)量最好,StOMP算法重構(gòu)圖像的質(zhì)量最差;當(dāng)采用率低于70%時(shí),StOMP算法重構(gòu)圖像的質(zhì)量基本不變,但效果比較不理想。3種算法下,重構(gòu)圖像和原始圖像的均方差總體趨勢下都是隨采樣率的增加而減??;相同采樣率條件下,BP算法均方差最大,StOMP算法均方差最小,如圖5c。
綜合上面所得的結(jié)果可知,OMP算法同時(shí)能夠兼顧重構(gòu)時(shí)間和重構(gòu)質(zhì)量,是一種相對比較實(shí)用的重構(gòu)算法。
本文研究了應(yīng)用壓縮感知理論進(jìn)行圖像重構(gòu)的關(guān)鍵技術(shù),分析對比了信號稀疏表示所選取的正交變換基、編碼所構(gòu)造的測量矩陣和重構(gòu)算法等不同時(shí)對重構(gòu)圖像質(zhì)量的影響。實(shí)驗(yàn)結(jié)果表明:應(yīng)用壓縮感知理論,僅需30%的采樣率即可有效恢復(fù)圖像,當(dāng)采樣率達(dá)到50%時(shí)能夠獲得較高質(zhì)量的圖像,相較于傳統(tǒng)圖像重構(gòu)方法中要求采樣率至少為原始信號的2倍,大大降低了采樣過程中冗余數(shù)據(jù)過多,信息處理過慢的問題,同時(shí)保證了圖像信號恢復(fù)的質(zhì)量。
一般情況下,以相對較大的采樣率進(jìn)行采樣,使用小波變換基對信號稀疏,以諾斯萊特矩陣作為測量矩陣,以O(shè)MP算法進(jìn)行重構(gòu)能夠有效兼顧圖像恢復(fù)效率和恢復(fù)質(zhì)量。該研究對當(dāng)前應(yīng)用壓縮感知理論進(jìn)行圖像處理時(shí)在大大降低硬件采集系統(tǒng)的復(fù)雜度和耗費(fèi)成本的前提下,進(jìn)一步提高重構(gòu)圖像的速度和質(zhì)量,具有重要的應(yīng)用意義。
[1]Donoho D L.Compressed sensing[J].IEEE Transactions on Information Theory,2006,52(4):289-1306.
[2]Candès E J,Romberg J,Tao T.Robust uncertainty principles:Exact signal reconstruction from highly incomplete frequency information[J].IEEE Transactions on Information Theory,2006,52(2):489-509.
[3]Liu X,Yao X,Wang C,et al.Quantum limit of photoncounting imaging based on compressed sensing[J].Optics Express,2017,25(4):3286-3296.
[4]Lu P,Xu Z,Lu X,et al.Digital image information encryption based on compressive sensing and double randomphase encoding technique[J].Optik-International Journal for Light and Electron Optics,2013,124(16):2514-2518.
[5]Sui X B,Chen Q,GU G H,et al.Infrared super-resolution imaging based on compressed sensing[J].Infrared Physics&Technology,2014,63:119-124.
[6]Feng L,Axel L,Chandarana H,et al.XD‐GRASP:Golden angle radial MRI with reconstruction of extra motion state dimensions using compressed sensing[J].Magnetic Resonance in Medicine,2016,75(2):775-788.
[7]Lustig M,Donoho D,Pauly J M.Sparse MRI:The application of compressed sensing for rapid MR imaging[J].Magnetic Resonance in Medicine,2007,58(6):1182-1195.
[8]Tavakoli A,Pourmohammad A.Image denoising based on compressed sensing[J].International Journal of Computer Theory and Engineering,2012,4(2):266.
[9]Zhou W P,Li Y,Liu Q S,et al.Fast compression and reconstruction of astronomical images based on compressed sensing[J].Research in Astronomy and Astrophysics,2014,14(9):1207-1214.
[10]Cetin M,Stojanovic I,Onhon O,et al.Sparsity-driven synthetic aperture radar imaging:Reconstruction,autofocusing,moving targets,and compressed sensing[J].IEEE Signal Processing Magazine,2014,31(4):27-40.
[11]Radwell N,Mitchell K J,Gibson G M,et al.Single-pixel infrared and visible microscope[J].Optica,2014,1(5):285-289.
[12]Chen X,Du Z,Li J,et al.Compressed sensing based on dictionary learning for extracting impulse components[J].Signal Processing,2014,96:94-109.
[13]Zhang Y,Dong Z,Phillips P,et al.Exponential wavelet iterative shrinkage thresholding algorithm forcompressed sensing magnetic resonance imaging[J].Information Sciences,2015,322:115-132.
[14]Bigot J,Boyer C,Weiss P.An analysis of block sampling strategies in compressed sensing[J].IEEE Transactions on Information Theory,2016,62(4):2125-2139.
[15]Paleo P,Desvignes M,Mirone A.A practical local tomography reconstruction algorithm based on a known subregion[J].Journal of Synchrotron Radiation,2017,24(1):257-268.
[16]Kher R,Patel Y.Medical image compression framework based on compressive sensing,DCT and DWT[J].Biol,2017,2(2):1-4.
[17]Krzakala F,Mézard M,Sausset F,et al.Probabilistic reconstruction in compressed sensing:algorithms,phase diagrams,and threshold achieving matrices[J].Journal of Statistical Mechanics:Theory and Experiment,2012(8):P08009.
[18]盧佩,田敏,張銳敏.基于IS模型的視覺顯著性分析及算法改進(jìn)[J].石河子大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,33(3):385-390.Lu P,Tian M,Zhang R M.Visual saliency analysis and algorithm advance based on IS model[J].Journal of Shihezi University(Natural Science),2015,33(3):385-390.
[19]Zhang T,Zhong S,Zhu Z,et al.Reconstruction algorithm based on compressed sensing[J].Journal of Mechanical&Electrical Engineering,2014,31(6):805-818.
[20]程濤.信號重構(gòu)中測量矩陣性能的判據(jù)[J].探測與控制學(xué)報(bào),2016,38(4):72-76.Cheng T.Performance criteria of measurement matrices in signal reconstruction[J].Journal of Detection& Control,2016,38(4):72-76.
[21]Tropp J A,Gilbert A C.Signal recovery from random measurements via orthogonal matching pursuit[J].IEEE Trans-Actions on Information Theory,2007,53(12):4655-4666.
[22]石光明,劉丹華,高大化,等.壓縮感知理論及其研究進(jìn)展[J].電子學(xué)報(bào),2009,37(5):1070-1081.Shi G M,Liu D H,Gao D H,et al.Advance s in theory and application of compressed sensing[J].Acta Electronica Sinica,2009,37(5):1070-1081.
[23]Tan J,Carmon D,Baron D.Signal estimation with additive error metrics in compressed sensing[J].IEEE Transactions on Information Theory,2014,60(1):150-158.
[24]Liu Y,Song T,Zhuang Y,et al.Compressed sensing signal recovery via creditability-estimation based matching pursuit[J].IEICE Transactions on Fundamentals of Electronics,Communications and Computer Sciences,2015,98(6):1234-1243.
[25]Massa A,Rocca P,Oliveri G.Compressive sensing in electromagnetics-a review[J].IEEE Antennas and Propagation Magazine,2015,57(1):224-238.
Research on key techniques of image reconstruction based on compressed sensing
Liu Xiaoyong1,Lu Pei2*,Cao Haibing1,Tian Min2
(1 College of Science,Shihezi University,Shihezi,Xinjiang 832003,China;2 College of Information Science and Technology,Shihezi University,Shihezi,Xinjiang 832003,China)
The rapid imaging technology put forward higher demand due to the importance of image application in information,and compressed sensing providesanew way forimageprocessing.Thispaperstudiesthekey technologiesof image reconstruction based on compressed sensing,comparative analysis of the influences on the quality of reconstructed image is emphases.The experimental results show that only 50%of the sampling rate we can obtain high quality images using wavelet transform based on signal sparse,Hadamard matrix as the measurement matrix and orthogonal matching pursuit algorithm.The method availably reduces the redundant data and the visual effects,PSNR values and recovery efficiency of reconstructed image are improved.This research provides a strong basis and support on the theory and technology for the application of compressed sensing to seek optimal image restoration method.
compressed sensing;image reconstruction;sparse signal;measurement matrix;reconstruction algorithm
TP391.41
A
10.13880/j.cnki.65-1174/n.2017.05.020
1007-7383(2017)05-0641-07
2017-02-24
國家自然科學(xué)基金項(xiàng)目(61465011);石河子大學(xué)高層次人才項(xiàng)目(RCZX201420、RCZX201436)作者簡介:劉效勇(1976-),男,博士,從事信息安全和三維傳感研究,e-mail:llxxyy1017@shzu.edu.cn。
*通信作者:盧佩(1979-),女,副教授,從事信號與信息處理、目標(biāo)識別研究,e-mail:lupei0@163.com。