王錄濤王 微金 鋼

①(電子科技大學(xué)自動(dòng)化工程學(xué)院 成都 611731)

②(國家空氣動(dòng)力研究與發(fā)展中心 綿陽 621000)

基于Hankel-SVD的非平穩(wěn)超聲血流成像雜波抑制技術(shù)研究

王錄濤*①王 微②金 鋼②

①(電子科技大學(xué)自動(dòng)化工程學(xué)院 成都 611731)

②(國家空氣動(dòng)力研究與發(fā)展中心 綿陽 621000)

有效抑制由血管或血管周圍組織時(shí)變運(yùn)動(dòng)引起的非平穩(wěn)雜波對(duì)于提高診斷超聲彩色血流成像中血流動(dòng)力學(xué)參數(shù)描述的準(zhǔn)確性有著極其重要的意義。該文基于奇異值濾波技術(shù)提出一種改進(jìn)的非平穩(wěn)雜波自適應(yīng)抑制方法。該方法逐次利用單個(gè)慢時(shí)多普勒回波采樣矢量構(gòu)建Hankel矩陣，然后根據(jù)奇異值分解后得到的正交Hankel主成份所代表的頻域內(nèi)容，動(dòng)態(tài)選取高階Hankel主成份重構(gòu)多普勒血流信號(hào)，實(shí)現(xiàn)非平穩(wěn)雜波的有效抑制。為驗(yàn)證算法的有效性，分別對(duì)多普勒回波仿真模型合成數(shù)據(jù)與利用彩色超聲設(shè)備(Sonix RP)采集的頸動(dòng)脈血流基帶回波信號(hào)進(jìn)行濾波處理，然后采用滯一自相關(guān)估計(jì)法計(jì)算血流平均速度與功率并進(jìn)行成像。處理結(jié)果表明，相對(duì)于傳統(tǒng)IIR濾波方法與多項(xiàng)式回歸濾波技術(shù)，利用該文所提算法可對(duì)高強(qiáng)度、非平穩(wěn)雜波進(jìn)行充分抑制，提高血流估計(jì)精度，此外，該算法具有空間自適應(yīng)性，無需人為設(shè)定閾值參數(shù)以估計(jì)雜波空間維數(shù)，與現(xiàn)有基于特征分解的自適應(yīng)濾波方法相比，可以有效提高組織空間高強(qiáng)度時(shí)變運(yùn)動(dòng)時(shí)血流與組織的區(qū)分能力。

彩色血流成像；奇異值分解；自適應(yīng)雜波抑制；非平穩(wěn)雜波

1 引言

在超聲彩色血流成像技術(shù)中，雜波指由血管或血管周圍慢動(dòng)組織反射引起的多普勒回波分量，該分量的強(qiáng)度通常比血流信號(hào)強(qiáng)度要高40～100 dB，即使當(dāng)采樣容積位于血管中心位置處，由于采樣容積的3維特性、合成聲束旁瓣泄露、聲波多重反射以及距離模糊等因素的影響，血流回波信號(hào)中也包含大量的雜波信號(hào)[1,2]。雜波信號(hào)的存在導(dǎo)致血流動(dòng)力學(xué)參數(shù)估計(jì)偏差，嚴(yán)重降低了彩色血流成像系統(tǒng)的血流鑒別能力，因此充分抑制雜波成份對(duì)于提高彩色超聲成像質(zhì)量有著至關(guān)重要的意義[3,4]。

傳統(tǒng)的雜波抑制濾波器是一種高通濾波器，包括有限沖激響應(yīng)(FIR)濾波器、無限沖激響應(yīng)(IIR)濾波器與多項(xiàng)式回歸濾波器等。在彩色超聲血流成像系統(tǒng)中，有限的多普勒回波矢量采樣數(shù)據(jù)導(dǎo)致FIR濾波器的幅頻特性無法滿足高強(qiáng)度雜波抑制需求。與FIR濾波器相比，相同階數(shù)的IIR濾波器具有更加陡峭的過渡帶滾降特性與阻帶抑制特性，但當(dāng)樣本數(shù)據(jù)有限時(shí)，IIR濾波器的幅頻響應(yīng)由暫態(tài)響應(yīng)主導(dǎo)[5]。多項(xiàng)式回歸型濾波器采用正交多項(xiàng)式對(duì)雜波分量進(jìn)行最小二乘擬合以消除雜波成份的影響，為雜波濾波器的設(shè)計(jì)提供了新的途徑[6]。當(dāng)雜波為平穩(wěn)信號(hào)時(shí)，傳統(tǒng)靜態(tài)雜波濾波器可獲得良好的雜波抑制效果，但當(dāng)組織存在加速運(yùn)動(dòng)時(shí)，不充分的雜波抑制將導(dǎo)致血流參數(shù)估計(jì)產(chǎn)生較大偏差[7]。

針對(duì)雜波信號(hào)的非平穩(wěn)特性，文獻(xiàn)[8]提出了基于特征值分解的雜波濾波器(Eigen filter)設(shè)計(jì)方法。該方法采用對(duì)回波多普勒信號(hào)自相關(guān)矩陣進(jìn)行特征分解獲得的特征矢量構(gòu)造雜波濾波器，將特征矢量作為回歸濾波器的正交基矢量可使得期望信號(hào)的擬合具有最小均方誤差性質(zhì)，但特征值分解濾波器假設(shè)整個(gè)觀察范圍內(nèi)的多普勒矢量信號(hào)是平穩(wěn)的，不具有空間自適應(yīng)性，此外，雜波空間維數(shù)需要根據(jù)經(jīng)驗(yàn)值選擇，不恰當(dāng)?shù)倪x擇將導(dǎo)致血流參數(shù)估計(jì)嚴(yán)重偏離真實(shí)值[9?11]。

作為一種數(shù)據(jù)特征提取的有效方法，奇異值分解技術(shù)可以將數(shù)據(jù)主導(dǎo)特征用少數(shù)分解系數(shù)集中體現(xiàn)，具有良好的穩(wěn)定性和不變性，近年來在特征提取、信號(hào)去噪等領(lǐng)域的應(yīng)用引起了廣泛重視[12,13]。本文介紹了一種基于Hankel奇異值分解(Hankel-SVD)的非平穩(wěn)雜波抑制技術(shù)。該技術(shù)采用單個(gè)慢時(shí)回波多普勒矢量構(gòu)建Hankel矩陣，具有空間自適應(yīng)性；根據(jù)Hankel主成份代表的頻率內(nèi)容自適應(yīng)選擇雜波主成份階數(shù)，通過摒棄雜波Hankel主成分，保留血流Hankel主成份并重構(gòu)血流多普勒信號(hào)，實(shí)現(xiàn)非平穩(wěn)雜波的有效抑制。最后采用基于多普勒回波仿真模型合成數(shù)據(jù)的血流參數(shù)估計(jì)結(jié)果與人體頸動(dòng)脈血流成像實(shí)例，對(duì)本文所提算法針對(duì)高強(qiáng)度、非平穩(wěn)雜波抑制的有效性進(jìn)行了驗(yàn)證。

2 Hankel-SVD濾波

超聲血流成像技術(shù)利用超聲聲束對(duì)待成像區(qū)域進(jìn)行掃描，在每一個(gè)掃描方向形成有限個(gè)發(fā)射聲束，通過接收血紅細(xì)胞散射的多普勒信號(hào)估計(jì)血流動(dòng)力學(xué)參數(shù)。Hankel-SVD濾波器采用KL(Karhunen-Loeve)變換將解調(diào)后的回波多普勒矢量信號(hào)表示為一組正交基矢量的線性疊加，通過保留與血流信號(hào)相對(duì)應(yīng)的正交成份并進(jìn)行信號(hào)重構(gòu)，可有效降低雜波對(duì)血流參數(shù)估計(jì)的影響[11]。對(duì)于包含N個(gè)采樣值的慢時(shí)回波多普勒矢量x，其Hankel-SVD逼近可表示為

其中P(P＜N)為用于逼近x的正交基矢量的階數(shù)，T表示轉(zhuǎn)置，γk,?k分別為k階展開系數(shù)與正交基矢量，滿足下述正交條件：

為計(jì)算式(1)中的正交成份，將多普勒矢量x分割成部分重疊的數(shù)據(jù)段并重新組織以構(gòu)造Hankel矩陣：

式(3)的奇異值分解可表示為

其中，Ak表示秩為1的第k階Hankel成份，σk,uk, vk分別為奇異值、左奇異矢量(N維)、右奇異矢量(N?P+1維)。式(4)中的P個(gè)奇異值滿足σ1＞σ2＞…＞σP，因此Ak又可稱為第k階Hankel主成份。

超聲血流多普勒回波信號(hào)由雜波、血流信號(hào)與噪聲3種成份構(gòu)成，其中雜波成份是由心壁、瓣膜、血管壁及靜止的或慢速運(yùn)動(dòng)的肌肉引起的，平均頻率較低但強(qiáng)度較高，在回波信號(hào)中占主導(dǎo)地位，而血流信號(hào)是由血紅細(xì)胞散射的多普勒信號(hào)構(gòu)成，強(qiáng)度較低，但頻率成份高于雜波信號(hào)[1]。文獻(xiàn)[6]對(duì)超聲多普勒回波自相關(guān)矩陣進(jìn)行特征值分解，得到了超聲多普勒信號(hào)特征值譜，并結(jié)合超聲回波多普勒信號(hào)特性對(duì)特征值譜分布與回波信號(hào)不同成份的能量、頻率特性構(gòu)成之間的關(guān)系進(jìn)行了研究，認(rèn)為與包含低頻成份的特征矢量相對(duì)應(yīng)的較大特征值反映了雜波成份構(gòu)成，與包含高頻成份的特征矢量相對(duì)

其中Kc為雜波Hankel主成份階數(shù)。

分析Hankel矩陣的構(gòu)造過程可知，在構(gòu)造Hankel矩陣A時(shí)，輸入信號(hào)矢量中的每個(gè)元素被沿反對(duì)角線方向重復(fù)配置，雖然經(jīng)奇異值分解，反對(duì)角線方向的元素不再相等，但仍可沿Ak反對(duì)角線方向?qū)Υ硗惠斎霐?shù)據(jù)的所有元素求平均值來重建血流多普勒回波信號(hào)所需的正交基矢量γk?k。應(yīng)的數(shù)值較小的高階奇異值則集中體現(xiàn)了血流多普勒成份特征，而噪聲信號(hào)能量則均勻分布在所有奇異值上。因此，當(dāng)雜波分量包含的Hankel主成份階數(shù)已知時(shí)，可由Ak計(jì)算得到的正交基矢量γk?k，進(jìn)而按照式(1)實(shí)現(xiàn)血流多普勒信號(hào)的準(zhǔn)確重構(gòu)：

3 雜波Hankel主成份階數(shù)估計(jì)

為充分抑制雜波成份，必須準(zhǔn)確判斷Hankel主成份Ak是否為雜波成份?；谔卣髦捣纸獾碾s波濾波器設(shè)計(jì)方法根據(jù)雜波成份與血流成份對(duì)奇異值譜的貢獻(xiàn)不同，通過將分解獲得奇異值與預(yù)先設(shè)置的經(jīng)驗(yàn)閾值比較，估計(jì)雜波Hankel主成份階數(shù)[11]。當(dāng)雜波成份在回波信號(hào)中處于主導(dǎo)時(shí)，可實(shí)現(xiàn)雜波子空間維數(shù)的準(zhǔn)確估計(jì)，但對(duì)于高頻超聲血流成像，由于血流信號(hào)的Rayleigh散射特性，雜波與血流信號(hào)的功率比值(CBR)接近0 dB，采用閾值法將導(dǎo)致雜波子空間維數(shù)估計(jì)產(chǎn)生較大偏差[15]。此外，組織空間生理特性差異與組織運(yùn)動(dòng)的空間、時(shí)間不確定性，導(dǎo)致回波信號(hào)隨空間、時(shí)間變化存在巨大差異，因此很難通過一個(gè)特定的閾值準(zhǔn)確實(shí)現(xiàn)雜波信號(hào)與血流多普勒信號(hào)的分離。

由超聲多普勒回波成份構(gòu)成與奇異值分解原理可知，Hankel-SVD的正交基矢量包含了多普勒回波不同成份的頻率信息，雜波成份通常為低頻分量，而血流信號(hào)則由高頻多普勒成份構(gòu)成[6,16]。因此，可根據(jù)?k所反映的頻率信息判定Hankel主成份Ak是否包含于雜波子空間，進(jìn)而估計(jì)雜波Hankel主成份的階數(shù)。首先利用?k計(jì)算第k階平均多普勒頻率：

其中，TPRI為發(fā)射脈沖重復(fù)間隔。fD,1為一階Hankel主成份，可用于區(qū)分當(dāng)前多普勒矢量是否由雜波成份所主導(dǎo)，如為雜波成份主導(dǎo)，則可以根據(jù)當(dāng)前多普勒矢量的瞬時(shí)雜波頻率估計(jì)雜波帶寬：其中，fi,ct為瞬時(shí)雜波頻率，可由當(dāng)前采樣的復(fù)多普勒矢量計(jì)算得出：

其中I與Q分別為復(fù)多普勒信號(hào)采樣值的實(shí)部與虛部。根據(jù)fD,1與fth,ct可確定雜波頻率范圍，進(jìn)而對(duì)Hankel主成份Ak的性質(zhì)進(jìn)行判決：

當(dāng)fD,k(k＞1)滿足式(9)時(shí)，說明Ak仍為雜波成份，此時(shí)，需增加k值進(jìn)行下一步搜索，即將fD,k+1代入式(9)，驗(yàn)證fD,k+1是否位于雜波頻率范圍內(nèi)，雜波Hankel主成份階數(shù)為使式(9)不成立的最小k值，即Kc=k。增加k值使得雜波成份逐次從回波多普勒信號(hào)中抽取出來，當(dāng)k＞Kc時(shí)，則雜波成份被有效剝離，Ak以血流成份為主導(dǎo)，根據(jù)式(5)可準(zhǔn)確重構(gòu)血流多普勒信號(hào)，而累加大于Kc的奇異值并進(jìn)行歸一化，則可得到濾波后的血流平均功率：

4 血流參數(shù)估計(jì)仿真

為了對(duì)本文設(shè)計(jì)的濾波器的雜波抑制性能進(jìn)行量化評(píng)估，采用文獻(xiàn)[16,17]提出的回波合成模型產(chǎn)生仿真多普勒回波信號(hào)，仿真參數(shù)如表1所示，濾波處理后的數(shù)據(jù)經(jīng)滯一自相關(guān)估計(jì)處理以估計(jì)血流平均多普勒頻率與功率。為分析不同流速時(shí)的濾波器雜波抑制性能，設(shè)計(jì)了20組血流信號(hào)，其平均多普勒頻率在0～0.5發(fā)射脈沖重復(fù)頻率(PRF)范圍內(nèi)均勻分布。

表1 壁濾波器性能比較仿真參數(shù)取值表

圖1給出了200次獨(dú)立仿真實(shí)驗(yàn)得到的血流平均多普勒頻率與功率估計(jì)結(jié)果。在圖1中，IIR-Prj為六階投影初始化Chebyshev IIR濾波器，通帶歸一化截止頻率為0.15，為保證相近的阻帶抑制性能，Reg-Pol多項(xiàng)式回歸濾波器采用五階Legendre 多項(xiàng)式重構(gòu)血流信號(hào)，Eigen濾波器雜波空間維數(shù)設(shè)為5。

圖1(a)給出了估計(jì)平均多普勒頻率比較結(jié)果，虛線表示理想估計(jì)平均多普勒頻率。從圖1(a)中可以看出采用Reg-Pol與Eigen得到的平均多普勒頻率估計(jì)精度優(yōu)于IIR-Prj，這是因?yàn)椴捎谜换瘮?shù)重構(gòu)血流信號(hào)可有效提升濾波器的雜波抑制性能，有效降低因雜波殘留導(dǎo)致的頻率估計(jì)偏差。而對(duì)比Reg-Pol與Eigen則可以發(fā)現(xiàn)，Eigen可以獲得更小的逼近誤差，這是由于Eigen將回波信號(hào)作為隨機(jī)過程處理，根據(jù)雜波的統(tǒng)計(jì)特性構(gòu)造雜波子空間的基矢量，當(dāng)雜波子空間維數(shù)小于信號(hào)空間維數(shù)時(shí)，可使得雜波逼近最小均方誤差最小[6]。相對(duì)于其它3種被計(jì)較濾波算法，當(dāng)輸入血流平均頻率從100 Hz增加到200 Hz時(shí)， Hankel-SVD可以顯著提高血流速度估計(jì)精度。

圖1(b)示出了平均功率估計(jì)比較結(jié)果，其中IIR-Prj濾波處理后的數(shù)據(jù)仍保留較多雜波成份，導(dǎo)致功率估計(jì)高于其它3種雜波濾波器。Reg-Pol濾波器具有最強(qiáng)的阻帶抑制性能，但過渡帶較寬，而在更寬的頻率范圍內(nèi)，Hankel-SVD功率估計(jì)值更接近于理想值。

5 頸動(dòng)脈血流成像比較

為驗(yàn)證本文算法在活體診斷應(yīng)用中的有效性，利用Sonix RP (Ultrasonix Inc., Vancouver, Cannada)超聲成像系統(tǒng)采集人體頸動(dòng)脈基帶多普勒回波數(shù)據(jù)，探頭中心頻率為6.6 MHz，脈沖重復(fù)頻率為2.5 kHz，發(fā)射脈沖重復(fù)次數(shù)為10，系統(tǒng)采樣頻率為40 MHz，軸向采樣750點(diǎn)，側(cè)向掃描50線。對(duì)雜波抑制后的采樣數(shù)據(jù)，采用滯一自相關(guān)估計(jì)法計(jì)算每個(gè)空間采樣點(diǎn)的血流速度與功率并進(jìn)行成像。

圖2 給出了一個(gè)完整心動(dòng)周期50幀數(shù)據(jù)中心臟收縮期的一幀數(shù)據(jù)經(jīng)4種雜波濾波器處理后的，未設(shè)置速率門限與功率門限的彩色多普勒血流成像結(jié)果。 其中，IIR-Prj為二階、歸一化截止頻率為0.05的投影初始化IIR濾波器，Reg-Pol為二階Legendre多項(xiàng)式回歸濾波器，構(gòu)建Eigen濾波器自相關(guān)矩陣的慢時(shí)多普勒矢量數(shù)為750，雜波空間維數(shù)固定為2。

在圖2中，不同的顏色代表不同的速度估計(jì)結(jié)果。在遠(yuǎn)離血管區(qū)域雜波成份被充分抑制，噪聲成份使得速度估計(jì)趨向于系統(tǒng)抗混疊速度，在圖像中表現(xiàn)為紅、藍(lán)交錯(cuò)分布。綠色區(qū)域表示未能被充分抑制的雜波成份，在血管壁、血管附近組織以及靠近內(nèi)血管壁的血流區(qū)域均有所分布。對(duì)比圖2可知，圖2(d) 中雜波殘留區(qū)域最小，圖2(c)中最大，圖2 (b)優(yōu)于圖2(a)。這是因?yàn)镽eg-Pol具有更加陡峭的過渡帶，因而雜波性能優(yōu)于IIR-Prj，但在血管壁上緣區(qū)域仍殘留一定數(shù)量的雜波； Eigen濾波器假設(shè)雜波沿聲束軸向是平穩(wěn)的，當(dāng)過多的多普勒矢量用于構(gòu)建自相關(guān)矩陣且濾波器階數(shù)選擇不當(dāng)時(shí)，雜波抑制性能嚴(yán)重惡化，大量殘留的雜波導(dǎo)致血流估計(jì)速度明顯低于其它3種方法；Hankel-SVD根據(jù)慢時(shí)方向單個(gè)多普勒回波矢量的雜波成份所反映的頻率特性，動(dòng)態(tài)調(diào)整Hankel雜波主成份階數(shù)，因而能有效抑制由于組織空間高強(qiáng)度時(shí)變運(yùn)動(dòng)導(dǎo)致的非平穩(wěn)雜波，提高血流與組織區(qū)分能力。

圖3為濾波后功率多普勒成像比較，顯示為歸一化多普勒功率，濾波器參數(shù)設(shè)置同圖2。從圖3中可以看出，經(jīng)4種濾波器進(jìn)行濾波處理后，雜波成份均得到了不同程度的抑制，遠(yuǎn)離血管區(qū)域的雜波抑制最為充分，而在血管壁及其附近區(qū)域，血管的非勻速搏動(dòng)使得在血管壁邊緣仍殘留一定部分的較強(qiáng)功率雜波。圖3(d)為Hankel-SVD濾波處理結(jié)果，可以看出雜波功率抑制最為充分，因而血流區(qū)域平均功率值明顯高于其它3種濾波算法。

圖1 濾波器雜波抑制性能比較

圖2 彩色多普勒血流成像

圖3 功率多普勒成像

6 結(jié)論

針對(duì)超聲彩色血流成像中的非平穩(wěn)雜波，本文提出了一種改進(jìn)的基于奇異值分解的雜波抑制算法。該方法逐次采用慢時(shí)方向單個(gè)回波多普勒矢量構(gòu)建Hankel矩陣，利用Hankel-SVD分解獲得的Hankel主成分作為回歸濾波器的正交基函數(shù)重構(gòu)血流多普勒信號(hào)，可使雜波擬合最小均方誤差最小，因而雜波抑制性能優(yōu)于傳統(tǒng)采用多項(xiàng)式函數(shù)作為正交基函數(shù)的回歸型雜波抑制方法。根據(jù)Hankel主成份包含的頻率內(nèi)容區(qū)分雜波與血流信號(hào)，當(dāng)雜波成份為主導(dǎo)時(shí)，通過估計(jì)雜波瞬時(shí)帶寬與一階Hankel主成分的平均多普勒頻率確定雜波多普勒頻率范圍，并與高階Hankel主成分平均多普勒頻率進(jìn)行比較，動(dòng)態(tài)選擇雜波Hankel主成份階數(shù)，實(shí)現(xiàn)非平穩(wěn)雜波的有效抑制。基于多普勒回波模型合成數(shù)據(jù)的仿真實(shí)驗(yàn)與頸動(dòng)脈血流與功率成像結(jié)果表明該算法在有效抑制由血管加速運(yùn)動(dòng)導(dǎo)致的非平穩(wěn)、強(qiáng)雜波的同時(shí)充分保留血流成份，改善低速血流參數(shù)的估計(jì)精度。由于僅采用單個(gè)慢時(shí)方向的多普勒回波矢量構(gòu)建Hankel矩陣，本文算法具有空間自適應(yīng)性，與Eigen雜波濾波器相比，無需人為設(shè)定閾值參數(shù)以估計(jì)雜波空間維數(shù)，因而可有效提升組織空間高強(qiáng)度時(shí)變運(yùn)動(dòng)時(shí)的血流成像質(zhì)量。

[1] Oeltze S, Lehmann D J, and Kuhn A. Blood flow clustering and applications in virtual stinting of intracranial aneurysms[J]. IEEE Transactions on Visualization and Computer Graphics, 2014, 20(5): 686-701.

[2] Park G Y, Yeo S M, and Yoon C H. New adaptive clutter rejection based on spectral decomposition and tissue acceleration for ultrasound color Doppler imaging[C]. IEEE International Ultrasonics Symposium, Prague, 2013: 1484-1487.

[3] Masuda K, Bossard A, and Sugano Y. Reconstruction and error detection of blood vessel network from ultrasound volume data[C]. IEEE 26th International Symposium on Computer-Based Medical Systems, Porto, 2013: 497-501.

[4] Yap C H, Thiele K, and Wei Q F. Novel method of measuring vascular regurgitation using three-dimensional nonlinear curve fitting of Doppler signals within the flow convergence zone[J]. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 2013, 60(7): 1295-1311.

[5] Kadi A and Loupas T. On the performance of regression and step-initialized iir clutter filters for color Doppler systems in diagnostic medical ultrasound[J]. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 1995, 42(3): 927-937.

[6] Bjaerum S, Torp H, and Kristoffersen K. Clutter filters design for ultrasound color flow imaging[J]. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 2002, 49(2): 693-704.

[7] Yoo Y and Kim Y. New adaptive clutter rejection for ultrasound color Doppler imaging: in vivo study[J]. Ultrasound in Medicine & Biology, 2010, 36(3): 480-487.

[8] Kruse D E and Ferrara K W. A new high resolution color flow system using an eigen decomposition-based adaptive filter for clutter rejection[J]. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 2002, 49(12): 1384-1399.

[9] Yu A C H. Eigen-based signal processing methods for ultrasound color flow imaging[D]. [Ph.D. dissertation], University of Toronto, 2007.

[10] Song F, Zhang D, and Gong X. Performance evaluation of eigendecomposition-based adaptive clutter filter for color flow imaging[J]. Ultrasonics, 2006, 44(12): 67-71.

[11] Yu A C H and Lovstaken L. Eigen-based clutter filter design for ultrasound color flow imaging: a review[J]. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 2010, 57(5): 1096-1111.

[12] 劉松, 張水蓮, 李世凱, 等. 一種低信噪比損失SVD 濾波部分頻帶干擾抑制算法[J]. 電路與系統(tǒng)學(xué)報(bào), 2013, 18(1): 70-75.

Liu Song, Zhang Shui-lian, Li Shi-kai, et al.. A low distortion algorithm based singular value decomposition for partialband interference rejection[J]. Journal of Circuits and Systems, 2013, 18(1): 70-75.

[13] 徐彥凱, 雙凱. 自適應(yīng)奇異值分解瞬變信號(hào)檢測研究[J]. 電子與信息學(xué)報(bào), 2014, 36(3): 583-588.

Xu Yan-kai and Shuang Kai. Detection of transient signal based on adaptive singular value decomposition[J]. Journal of Electronics & Information Technology, 2014, 36(3): 583-588.

[14] Yu A C H and Lovstakken L. Eigen-based clutter filters for clolor flow imaging: single-ensemble vs. multi-ensemble approaches[C]. IEEE International Ultrasonics Symposium, New York, 2007: 1101-1104.

[15] Park G Y, Yeo S M, and Yoon C H. New adaptive clutter rejection based on spectral analysis for ultrasound color Doppler imaging: phantom and in vivo abdominal study[J]. IEEE Transactions on Biomedical Engineering, 2014, 61(1): 55-63.

[16] Yu A C H, Johnston K W, and Cobbold R S C. Frequency based signal processing in ultrasound color flow imaging[J]. Canadian Acoustics, 2007, 35(2): 11-23,

[17] Yu A C H and Cobbold R S C. Single ensemble based eigenprocessing methods for color flow imaging part 1 the Hankel-SVD filter[J]. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 2008, 55(3): 559-573.

王錄濤： 男，1979年生，博士后，研究方向?yàn)殛嚵行盘?hào)處理與醫(yī)學(xué)成像處理.

王 微： 男，1982年生，博士生，研究方向?yàn)槔走_(dá)信號(hào)處理.

金 鋼： 男，1958年生，研究員，博士生導(dǎo)師，從事雷達(dá)信號(hào)處

理、圖像處理與流體力學(xué)方向研究.

Non-stationary Clutter Rejection Based on Hankel-SVD for Ultrasound Color Flow Imaging

Wang Lu-tao①Wang Wei②Jin Gang②

①(School of Automation Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China)

②(China Aerodynamics Research & Development Center, Mianyang 621000, China)

Effective rejection of the time-varying clutter originating from slowly moving vessels and surrounding tissues is very important for depicting hemodynamics in ultrasound color Doppler imaging. In this paper, a new adaptive clutter rejection method based on Hankel Singular Value Decomposition (Hankel-SVD) is presented for suppressing non-stationary clutter. In the proposed method, a Hankel data matrix is created for each slow-time ensemble. Then the orthogonal principle Hankel components can be obtained through the SVD of the Hankel data matrix. It achieves non-stationary clutter suppression by reconstructing the flow signal with only the high order principle Hankel components, which are estimated from the frequency content carried by the principle Hankel components. To assess its efficiency, the proposed Hankel-SVD based method is applied to synthetic slow-time data obtained from a Doppler flow model and carotid arterial complex baseband data acquired by a commercial ultrasound system (Sonix RP). The resulting flow and power images show that the proposed method outperforms the traditional IIR and polynomials regression filter in attenuation of high intense non-stationary clutter signal. It is also adaptive to highly spatially-varying tissue motion and can automatically select the order of the filter, which leads to improved distinguishing between blood and tissue regions compared to other eigen-based filters.

Color flow imaging; Singular Value Decomposition (SVD); Adaptive clutter rejection; Non-stationary clutter

TN911.7; R445.1

： A

：1009-5896(2015)04-0830-06

10.11999/JEIT140893

2014-07-07收到，2014-11-24改回

*通信作者：王錄濤 wltuestc@163.com