王霞

(山東廣播電視大學(xué)計算機與通信學(xué)院，山東 濟南 250014)

?

凸陣式超聲數(shù)字掃描變換技術(shù)的軟件實現(xiàn)

王霞

(山東廣播電視大學(xué)計算機與通信學(xué)院，山東 濟南 250014)

摘要：針對B 型超聲系統(tǒng)中采用硬件電路完成數(shù)字掃描變換缺乏靈活性且成本較高的問題，采用R-theta插補法以及幀相關(guān)降噪實現(xiàn)了凸陣式超聲數(shù)字掃描變換。結(jié)果表明，采用本方法能夠快速地顯示圖像，最終圖像顯示幀頻為30 f/s，滿足B型超聲系統(tǒng)的實時顯像要求。

關(guān)鍵詞：數(shù)字掃描變換; 坐標(biāo)變換;插補處理;B超

B超的系統(tǒng)原理是超聲探頭中的換能器按照一定的順序依次發(fā)射和接收超聲波，形成一個平面，從而得到一幅二維截面黑白聲像圖[1]。在這種顯示圖像為二維圖像的診斷設(shè)備中，如果將超聲回波信號直接顯示到屏幕上，會使人眼感到閃爍，尤其在使用凸陣式超聲探頭掃查時，閃爍感會更加明顯，這就需要應(yīng)用一種圖像存儲器，位置在超聲掃描與圖像顯示之間，讓存入圖像數(shù)據(jù)的速度與超聲聲束掃查的速度相同，而令讀出圖像數(shù)據(jù)的速度適當(dāng)?shù)靥嵘?，采用這種方法就可以令圖像穩(wěn)定地顯示。這就是超聲數(shù)字掃描變換技術(shù)[2]。

目前，大多數(shù)數(shù)字掃描變換還是由硬件來完成的。彭勇等[3]以CORDIC 算法和雙線性插值，在 FPGA 中實現(xiàn)。譚震威等[4]同樣以FPGA為系統(tǒng)的控制中心，實現(xiàn)了數(shù)字掃描變換。隨著高速計算機運算和數(shù)字信號處理技術(shù)的飛速發(fā)展，計算機系統(tǒng)的硬、軟件性能已經(jīng)有了大幅度的提高，而且價格越來越低，利用計算機自身的硬件資源結(jié)合軟件開發(fā)技術(shù)，超聲診斷檢測儀的性能將得到大幅度的提高。這樣不僅可以節(jié)省大量硬件成本，而且增加了系統(tǒng)的靈活性，便于之后的改進(jìn)和升級。 本文摒棄了傳統(tǒng)的硬件實現(xiàn)方式，直接利用計算機自身軟硬件編程來實現(xiàn)數(shù)字掃描變換。

1凸陣式數(shù)字掃描變換技術(shù)

關(guān)于凸陣式數(shù)字掃描變換技術(shù)，需要涉及到兩個步驟，即超聲掃查和顯示掃描。其中顯示掃描是本文的主要工作，主要包括坐標(biāo)變換和插補處理兩部分。

1.1坐標(biāo)變換

凸陣式數(shù)字掃描變換，其超聲前端掃查得到的是極坐標(biāo)形式的回波信號，而顯示器是以直角坐標(biāo)的形式將數(shù)據(jù)進(jìn)行顯示的，這就需要將原本的極坐標(biāo)形式的數(shù)據(jù)變換到直角坐標(biāo)系中正確地顯示[5]。

圖1　坐標(biāo)變換方案Fig.1　Coordinate transform scheme

如圖1中所示，512×512個數(shù)據(jù)存儲在矩陣中，坐標(biāo)變換首先是取一個極坐標(biāo)點，即圖中的(Xp，Yp)，然后映射到標(biāo)準(zhǔn)投影S平面，并且一部分點會作為之后的插補像素點對S平面得到的圖像進(jìn)行插補，最終圖像還需要將坐標(biāo)原點進(jìn)行變換才能適應(yīng)Windows默認(rèn)的坐標(biāo)系，最終得到D平面，即為顯示在顯示器上面的圖像。

關(guān)于極坐標(biāo)到直角坐標(biāo)的變換，直接采用的是公式法[6]，即

(1)

式中，x0、y0是平面S相對于平面D的水平和垂直偏移，需要根據(jù)實際凸陣探頭的內(nèi)圓半徑以及抽取數(shù)目決定具體的值；θ是掃描聲束偏轉(zhuǎn)角；r是沿掃描聲束矢徑上的采樣深度。

圖2　圓弧插補法Fig.2　Arc interpolation method

1.2插補處理

隨著探查深度的增加，相鄰掃描線間距變大，以致掃描線之間某些顯示像素沒有給予回波數(shù)據(jù)賦值而形成顯示器上的“黑洞”顯示。許多黑洞云集成“云紋狀”的線條，使得圖像的質(zhì)量下降。所以需要利用周圍的采樣點進(jìn)行插補來填補這些空白像素點。

圓弧插補即通常所說的二維圖像插補，本項目采用的是R-theta插補法[7](圖2)。在這種插補方法中，每一個空缺的顯示平面上直角坐標(biāo)形式的像素都是從相鄰的極坐標(biāo)形式的采樣數(shù)據(jù)按照線性插值方法插補出來的。對于二維的超聲圖像來說，使用弧形插補原則上應(yīng)該是比較適合的。

凸陣掃查的掃查線是按照等角度增量Δθ均勻分布的，掃查線上采樣點按照等間隔Δr均勻分布。設(shè)f(x,y)為直角坐標(biāo)像素數(shù)據(jù)，s(i，j)為極坐標(biāo)采樣數(shù)據(jù)，則有

f(x，y)=z(i)(1-θerr)+z(i+ 1)θerr,

(2)

式中，z(i)是第i條聲束經(jīng)線的中值，即

z(i)=s(i，j)(1 -Rerr)+s(i，j+ 1)Rerr,

(3)

(4)

分別為半徑方向和角度方向的權(quán)重。

將公式(2)和(3)進(jìn)行合并，得到最終的插補公式

f(x，y)=s(i，j)(1 -θerr)+s(i，j+ 1)(1-Rerr)θerr

+s(i+ 1，j)Rerr(1-θerr) +s(i+1，j+ 1)(1-Rerr)(1-θerr)。

(5)

圖3　圓弧插補流程圖Fig.3　Flowchart of arc interpolation

本文的弧形插補即是根據(jù)上述公式實現(xiàn)的。圖3所示為圓弧插補的流程圖，其中flag為是否為空缺像素的標(biāo)志。假如凸陣探頭的掃查角度是不變的，抽取數(shù)目也是不變的話，那么每幅圖像中每個被插補點的位置就是固定的，所以插補的點的位置無需每幅圖像都計算一次，只需要分別用一個矩陣將Rerr和θerr保存，后續(xù)的圖像插補點直接取矩陣中的數(shù)據(jù)進(jìn)行計算即可。直到掃查的角度以及抽取數(shù)目任何一個發(fā)生變化，再重新計算存入矩陣。這樣一來，只有當(dāng)參數(shù)發(fā)生變化時所用時間較長，其余時間會大大縮短，能夠滿足實時顯示的要求。

2實驗結(jié)果

本文通過上述坐標(biāo)變換以及插補處理過程，完成了數(shù)字掃描變換技術(shù)的關(guān)鍵部分，得到以下實驗結(jié)果，并對結(jié)果圖像進(jìn)行了幀相關(guān)降噪。

2.1實驗結(jié)果

圓弧插補后的超聲圖像如圖4所示，抽取數(shù)目為16、36。實際上模擬圖像是為了邊界明顯，直接將兩輪0～255的按行依次排列的數(shù)據(jù)進(jìn)行處理，行為深度，列為角度，而真實的超聲系統(tǒng)上傳上來的數(shù)據(jù)并不是按行排列的，而是按照列進(jìn)行存儲的，即列為深度，行為角度，所以，還需要將矩陣轉(zhuǎn)90度再進(jìn)行坐標(biāo)變換和數(shù)據(jù)插補才能夠看到正的超聲圖像。OpenCV有提供此方法的函數(shù)cvTranspose[8]：

函數(shù)原型：void cvTranspose ( const CvArr * src, CvArr * dst ) ;

src：輸入矩陣，dst：目標(biāo)矩陣。

圖4　弧形插補后的超聲圖像Fig.4　Ultrasonic image after arc interpolation

2.2后期圖像處理

對于相對靜止的人體臟器作超聲掃查時，對所獲得的B超圖像可進(jìn)行多幅圖像疊加平均處理，使圖像上的噪聲得到抑制，如圖5所示，這就是幀相關(guān)技術(shù)的思想。幀相關(guān)處理是一種多圖像平均法，它將同一像素點的前、后幾幅圖像進(jìn)行處理，可以取平均值、最大值或新的像素值，用來達(dá)到減少噪聲干擾的目的[9]。實際項目中取了4幅圖像進(jìn)行平均來減少噪聲，動態(tài)圖效果更加明顯。

圖5　幀相關(guān)對比圖Fig.5　Comparison chart of frame correlation

3結(jié)語

本文實現(xiàn)了軟件化的數(shù)字掃描變換，并且最終圖像顯示幀頻為30 f/s，能夠滿足實際需求。但本文中凸陣診斷的最小深度就是抽取最小時，是512個采樣點達(dá)到的深度，還可以從這512個采樣點再截取小段進(jìn)行擴展設(shè)計，使得超聲系統(tǒng)能夠進(jìn)行更加精確的觀察。

參考文獻(xiàn)：

[1]錢蘊秋，周曉東，張軍，等. 實用超聲診斷手冊[M]. 2版. 北京：人民軍醫(yī)出版社， 2011.

[2]吳水才，楊春蘭,白燕萍. 醫(yī)學(xué)超聲及應(yīng)用[M]. 北京：北京工業(yè)大學(xué)出版社. 2012.

[3]彭勇，陳菲，董萬利. B超中數(shù)字掃描變換器的 FPGA 設(shè)計與實現(xiàn)[J]. 微型機與應(yīng)用，2013，32(12)：24-26.

[4]譚震威，鄭旭生，羅繼裕. 醫(yī)學(xué)超聲圖像數(shù)字變換方法的硬件實現(xiàn)[J].科技風(fēng)，2011(5):61.

[5]馮若，劉忠齊,姚錦鐘,等. 超聲診斷設(shè)備原理與設(shè)計[M]. 北京：中國醫(yī)藥科技出版社. 1993.

[6]陳紅. 基于Matlab的數(shù)字掃描變換器設(shè)計[J]. 現(xiàn)代電子技術(shù)，2015，38(8)：109-113.

[7]尚兆梅，陳波，彭勇. R-Theta 算法的一種快速實現(xiàn)方式[J]. 科學(xué)技術(shù)與工程，2010，10(31)：7803-7806.

[8]OpenCV.OpenCV參考手冊[EB/OL]. [2014-10-05].http://www.opencv.org.cn/index.php/Template:Doc.

[9]郭德全，楊紅雨，劉東權(quán)，等. 基于序列相關(guān)性的超聲圖像自適應(yīng)去噪[J]. 計算機應(yīng)用研究，2014，31(2)：600-604.

Software implementation of convex ultrasonic digital scan conversation

WANG Xia

(School of Computer and Telecommunication, University of Radio and Television of Shandong, Jinan 250014, China)

Abstract∶We employ R-theta interpolation and frame correlation to achieve convex ultrasonic digital scan conversion for bad flexibility and higher cost of hardware digital scan conversion in B-ultrasonic system.Results show that it can quickly display images and satisfy the requirement of real-time image display, final video frame frequency of 30 f/s.

Key words∶digital scan conversation;conversion of coordinates;interpolation;B-ultrasonic system

中圖分類號：TP391

文獻(xiàn)標(biāo)識碼：A

文章編號：1002-4026(2015)04-0107-04

作者簡介：王霞 (1989-)，女，助教，研究方向為醫(yī)學(xué)超聲圖像顯示。

收稿日期：2014-12-29

DOI:10.3976/j.issn.1002-4026.2015.04.020