郭禹姬，宋文愛

(中北大學(xué) 信息與通信工程學(xué)院 太原 030051)

0 引言

虛擬相控陣超聲散射CT成像方法，不需要反投影重構(gòu)，在一定程度上彌補(bǔ)了線性超聲CT方法由于線性假設(shè)或弱散射假設(shè)帶來的不足，提高了超聲檢測(cè)的空間分辨率力。但是由于該方法需要綜合考慮幅值和相位信息，所以需要對(duì)虛擬相控陣中的每一個(gè)陣元發(fā)射、其余陣元接收（包括發(fā)射陣元本身接收）到的所有信號(hào)進(jìn)行處理，信息量較大，所以影響了成像的速度。為了快速而且準(zhǔn)確進(jìn)行超聲CT成像，如何快速查找掃描區(qū)域（檢測(cè)區(qū)域）內(nèi)存在的散射源，是提高成像速度和成像質(zhì)量的關(guān)鍵[1]。

本文通過時(shí)間壓縮的方法，對(duì)信號(hào)進(jìn)行壓縮，通過相移得到的信號(hào)簇一維壓縮，能準(zhǔn)確找到散射界面在時(shí)間軸上的位置，并通過頻域信號(hào)的觀察對(duì)界面的空間位置達(dá)到相位定位。

1 相控陣成像技術(shù)介紹

在探測(cè)區(qū)域相互垂直的兩個(gè)邊分別排列一個(gè)m元的線性超聲陣列探頭（根據(jù)情況也可以扇形超聲陣列探頭）每個(gè)探頭陣元均可以作為發(fā)射和接收陣元，當(dāng)水平放置的陣列的第i個(gè)單元發(fā)射，垂直放置的陣元j（j=1，2，…，）接收時(shí)，其發(fā)射和接收示意圖如圖1所示[2]。

圖1 水平單陣元發(fā)射，垂直各陣元接收示意圖

采集到的信號(hào)是圖4試塊的時(shí)間序列示意圖如圖2所示。

圖2 24號(hào)陣元發(fā)射其他所有陣元接收原信號(hào)簇

圖3 掃描區(qū)域內(nèi)任意散射點(diǎn)P與其距離最近接收陣元關(guān)系示意圖

2 時(shí)間壓縮

當(dāng)檢測(cè)或掃描區(qū)域內(nèi)存在散射源時(shí)，接收信號(hào)的幅值增大，在相應(yīng)的接收陣元上出現(xiàn)局部最大值點(diǎn)，最大值點(diǎn)出現(xiàn)以前的接收序列開始段或反方向椎大值點(diǎn)出現(xiàn)以后的結(jié)尾段，認(rèn)為是冗余信號(hào)，可以進(jìn)行壓縮。

由于接收到的散射信號(hào)存在一定的相位差，為減小計(jì)算誤差、時(shí)間壓縮后不丟失有用的信息，所以nmin和nmax分別向左、向右移動(dòng)2。

這樣對(duì)于所獲得的所有波列的時(shí)間序列的時(shí)間值均壓縮到以n= nmin為起點(diǎn)，n= nmax為終點(diǎn)的范圍內(nèi)，并以該范圍內(nèi)的每個(gè)時(shí)間采樣點(diǎn)作為灰度圖像的坐標(biāo)y。圖2信號(hào)壓縮后信號(hào)如圖5(a)。

3 實(shí)驗(yàn)

圖4所示的鈦合金CSK-IA試塊為本文所研究的檢測(cè)對(duì)象。圖4中標(biāo)出了試塊的尺寸以及陣列傳感器的布置區(qū)域。圖4中，陣列傳感器A具有49個(gè)陣元，陣元寬度為1mm。傳感器頻率為6．25MHz,在被檢試塊中的波長(zhǎng)近似為1mm。

圖4中還標(biāo)出了可能引起散射的界面。界面1為半徑100mm的圓弧界面，界面2為深度91mm的水平界面，界面3為高度為6mm的豎直界面，界面4為深度85mm,寬度2mm 的水平界面，界面5為高度15mm的豎直界面，界面6為深度100mm的水平界面。

圖4 超聲散射CT的實(shí)驗(yàn)樣品及傳感器布置

我們對(duì)陣列傳感器A所獲得的信號(hào)進(jìn)行處理，圖5所示為24號(hào)晶片發(fā)射超聲波所有晶片接收到的回波信號(hào)組圖。圖5(a)表示接收到的原始回波信號(hào)，由于界面4在24號(hào)陣元的正下方，而且界面4與探頭具體舉例最小，從圖中可以看到界面4信號(hào)最早反射回來，且24號(hào)陣元最先接收到信號(hào)，圖5(b)所示為原始回波信號(hào)的沒有經(jīng)過移相的所有49個(gè)信號(hào)的累加和，可以看出由于各個(gè)信號(hào)之間存在相位差，疊加后信號(hào)的幅值并不大。圖5(c)表示為根據(jù)協(xié)方差法對(duì)每一條信號(hào)線移相后的結(jié)果，從圖中可以看出，移相后的各個(gè)信號(hào)在相位上基本保持一致。圖5(d)為移相后的所有信號(hào)的累加和，與圖5(b)相比，信號(hào)幅值提高一倍多。

根據(jù)圖5(d)中圖像的投影分布，可以看到散射界面應(yīng)該位于投影分布的各個(gè)峰值處，因此找到投影分布圖中找到對(duì)應(yīng)的峰值，就可以確定是散射界面。具體算法如下：

(1)對(duì)灰度圖進(jìn)行一維投影，即沿圖像進(jìn)行灰度值的累加，得到圖像投影分布圖。

(2)從左向右搜索，通過動(dòng)態(tài)選取閾值法找到不同區(qū)間的極大值，根據(jù)極大值對(duì)應(yīng)的時(shí)間軸位置，即可找到散射界面相應(yīng)位置。

為了對(duì)散射界面的位置更加準(zhǔn)確識(shí)別，對(duì)已經(jīng)找到的散射區(qū)域做傅立葉變換。如圖6所示。

由圖6分析可以得出，49個(gè)陣元發(fā)射所在各個(gè)界面的傅立葉變換結(jié)果可以得出，峰值數(shù)不同，散射界面的個(gè)數(shù)也不同。圖6(a),(d)只有一個(gè)峰值，表示所在界面4，(b)有兩個(gè)峰值，表示界面2與界面3，依次類推。

通過傅里葉變換變換結(jié)果，可以得知這個(gè)界面在空間上的相對(duì)位置。

圖5 24號(hào)晶片發(fā)射所有晶片接收回波信號(hào)組圖

圖6 12號(hào)陣元，36號(hào)陣元三個(gè)散射區(qū)域傅立葉變換結(jié)果

4 結(jié)論

通過時(shí)間壓縮，時(shí)移方法能夠快速找到散射界面在時(shí)間軸上的位置，從而為超聲CT快速成像打下基礎(chǔ)。

[1] 倪文磊．超聲CT理論與方法綜述[J]． CT理論與應(yīng)用研究，2004，13(1)：50-55．

[2] Macovski A．Ultrasonic imaging using arrays [J]．Proceedings of IEEE, 1979, 67(4): 484-495．

[3] 趙霞． 基于相控陣超聲脫粘檢測(cè)系統(tǒng)研究[D]．太原：中北大學(xué),2006．

[4] 孫景鰲，蔡安妮，何擁軍． CT圖象重建的快速算法[J]．北京郵電大學(xué)學(xué)報(bào)，1995,18(3): 5-15．

[5] K A Dines and SA Goss．Computed Ultrasonic Reflcction tomography[J], IEEE Trans． 1987; UF FC-34(3):309～319．

[6] G．T．Herman． Image reconstruction from projections[J]．Academic Press, INC,1980．

[7] H．H．Barrett and W．Swindell． Radiological Imaging[J],vol．2 ,Academic Press, INC,1981．

[8] 劉波,李朝榮． 超聲CT成像方法及應(yīng)用[J]． 中國(guó)儀器儀表,2007(2):28-32．