齊建，陳益民，丁明躍，尉遲明

華中科技大學(xué)，武漢，430074

0 引言

在臨床診斷領(lǐng)域，超聲成像是目前應(yīng)用最為廣泛的一種成像模式。在過去的幾十年中，隨著科技的發(fā)展，超聲成像技術(shù)取得了長足的發(fā)展，從傳統(tǒng)的A型、B型、M型和多普勒超聲成像模式，進入到目前的三維超聲、彈性成像和靶向性造影劑諧波成像。

傳統(tǒng)的B型超聲成像系統(tǒng)所提供的是人體某一斷面的二維圖像，醫(yī)生必須根據(jù)自己的經(jīng)驗對多幅二維圖像在大腦中進行構(gòu)建以理解其三維解剖結(jié)構(gòu)。三維超聲的概念由Baun和Greewood在1961年提出，他們在采取一系列平行的人體器官二維超聲界面的基礎(chǔ)上，用疊加的方法獲得了人體器官的三維圖像[1]。

與二維超聲相比，三維超聲成像具有以下明顯的優(yōu)勢：(1) 圖像顯示直觀；(2)在醫(yī)學(xué)教學(xué)和手術(shù)規(guī)劃方面有廣泛的應(yīng)用前景；(3)可進行醫(yī)學(xué)診斷參數(shù)的精確測量；(4)可以縮短醫(yī)生診斷需要的時間。鑒于上述原因，三維超聲在疾病的定性、定位及定量診斷中發(fā)揮著舉足輕重的作用[2]。

當(dāng)前國內(nèi)外三維超聲圖像的獲得主要有：基于面陣探頭和基于機械掃描的兩種方式。

基于面陣探頭的三維超聲圖像的成像模式由Duke大學(xué)Vonn Ramm等人提出，圖1為面陣超聲探頭超聲換能器排列示意圖[3]。

圖1 二維面陣超聲探頭換能器陣列示意圖Fig.1 Transducer arrays of 2D planer ultrasound probe

二維面陣探頭也就是二維相控陣探頭，通過控制二維面陣中的陣元發(fā)射超聲束并實現(xiàn)空間的偏轉(zhuǎn)，直接采集三維數(shù)據(jù)。面陣的每一個陣元都配置了相應(yīng)的延遲線，采集數(shù)據(jù)時，只要控制不同陣元的延遲時間就能改變波束的指向。二維相控陣超聲探頭避免了復(fù)雜的機械定位裝置，采集數(shù)據(jù)的速度較快，并且較穩(wěn)定[4]。它能在兩個橫向方向上聚焦，聚焦聲場為針狀或棒狀，檢測的空間分辨率較高，相控陣系統(tǒng)不需要移動探頭即可在特定空間自由地實現(xiàn)偏轉(zhuǎn)和聚焦，極大改善了聲場的可達性和檢測速度。但二維面陣探頭的陣元數(shù)目大，通道數(shù)多，單陣元的尺寸小，加工困難。同時，陣元的靈敏度低，接收到的信號弱且信噪比低。若探頭由64×64的陣元組成，通道數(shù)變成了4096，技術(shù)難度大，不易實現(xiàn)，因此市場上這類探頭只有GE，Philips，Siemens等大公司才能生產(chǎn)，且價格較昂貴[5-7]。

基于機械掃描的三維超聲圖像主要是由一組二維超聲數(shù)據(jù)重建得到的，大致分為圖像采集、三維重建和三維圖像的顯示三個步驟[8]。其中圖像的采集是三維成像的關(guān)鍵步驟，直接影響重建結(jié)果的優(yōu)劣。一般來說，二維超聲探頭被固定在特點裝置內(nèi)，依靠馬達驅(qū)動超聲探頭掃描完成數(shù)據(jù)的采集。對于二維超聲探頭而言，傳統(tǒng)的掃描模式包括平行掃描(又稱為線性掃描，圖2)、基于前端掃描的扇掃掃描(圖3)和旋轉(zhuǎn)掃描(圖4)三種方式[8]。其中前兩種模式適用于腹部探頭，后一種掃描主要適用于腔內(nèi)探頭。

圖2 平行掃描Fig.2 Parallel scanning mode

圖3 傳統(tǒng)扇形掃描Fig.3 Traditional fan-shaped scanning mode

圖4 旋轉(zhuǎn)掃描Fig.4 Rotational scanning mode

平行掃描模式獲得的每一幅二維圖像之間相互平行[8]。掃描時，二維超聲探頭在馬達和機械裝置的控制下沿人體體表平行滑動，完成二維超聲圖像數(shù)據(jù)的采集。二維探頭平行移動的速度和數(shù)據(jù)采集的抽樣率的控制對三維重建結(jié)果有很大的影響?；谇岸藪呙璧纳葤攉@得的圖像相互間成等距夾角。掃描時，二維超聲探頭后端保持固定，前端（超聲發(fā)射端）扇形旋轉(zhuǎn)采集數(shù)據(jù)。但當(dāng)采用上述兩種掃描模式對病人胸部臟器掃描時，發(fā)射的超聲波易被人體肋骨遮擋，影響二維超聲圖像質(zhì)量，并造成重建后得到的三維超聲圖像質(zhì)量下降[10]。

為了避免扇形掃描和平行掃描模式的缺陷，本文設(shè)計了一種后端掃描的扇形掃描（圖5），二維超聲探頭前端固定，后端扇形旋轉(zhuǎn)。在此掃描模式下，如果探頭前端置于人體肋骨之間，由超聲探頭發(fā)射的超聲束將通過人體肋骨的間隙進入人體，成功采集二維數(shù)據(jù)且避免了肋骨的影響，極大提高了三維超聲圖像的質(zhì)量。

圖5 基于后端掃描的三維超聲的成像模式Fig.5 Fan-shaped scanning approach based on back-end scanning

1 系統(tǒng)架構(gòu)

為了驗證基于后端掃描模式的有效性，本文搭建了一個包括三維超聲圖像數(shù)據(jù)的采集、重建和顯示的系統(tǒng)。該系統(tǒng)由機械掃描裝置、超聲機、工作站（配有圖像采集卡的PC機）、單片機驅(qū)動控制芯片等構(gòu)成。工作站和超聲機由視頻輸出線相連接，單片機驅(qū)動控制芯片通過串口線和工作站通訊，系統(tǒng)架構(gòu)示意圖如圖6所示。

圖6 系統(tǒng)架構(gòu)示意圖Fig.6 System structure

以下我們就掃描裝置的機械設(shè)計，步進電機驅(qū)動的硬件設(shè)計，三維數(shù)據(jù)采集和重建結(jié)果分別介紹。

1.1 機械裝置

機械裝置主要包括探頭支架、探頭夾、弧形導(dǎo)航滑塊和弧形齒輪軌道等，如圖7所示。裝置采用鋁質(zhì)材料加工而成，步進電機固定在探頭支架兩側(cè)。選擇形狀、大小、力矩和歩矩角合適的步進電機兩臺，如寶來公司生產(chǎn)的兩相42BY100BX型號的步進電機。探頭夾可以進行伸縮，可固定當(dāng)前市面上的大部分二維超聲探頭。

為了確保二維超聲探頭前端固定、后端掃描，設(shè)計必須保證機械裝置的弧形軌道和弧形齒輪的半徑R與探頭固定的位置到探頭頂端的距離L始終保持相等，如圖8所示。此裝置設(shè)計的掃描范圍最大角度為60o。

圖7 機械裝置的整體示意圖1手柄；2探頭模型；3支架；4步進電機；5弧形導(dǎo)航滑塊；6弧形齒輪軌道；7探頭夾；Fig.7 Mechanical scanning device: (1)hand shank, (2) mode of the probe,(3)support, (4)stepper moter, (5) the arc navigation slider,(6) the arc gear track, (7) the probe clip

圖8 探頭與探頭夾位置關(guān)系Fig.8 The location relations between probe and probe clip.

1.2 步進電機的硬件控制和二維超聲數(shù)據(jù)的采集

二維超聲探頭由步進電機驅(qū)動掃描，其掃描速度和掃描角度通過步進電機單片機控制板細分驅(qū)動器和驅(qū)動程序控制。為了保證實驗過程中數(shù)據(jù)采集的穩(wěn)定性，步進電機的控制板由單片機和電機驅(qū)動電路組成。上位機PC工作站通過串口將命令傳給單片機，單片機將接收到的命令轉(zhuǎn)化為對應(yīng)的驅(qū)動信號輸出，該信號進入驅(qū)動電路部分，輸出步進電機所需要的脈沖驅(qū)動信號。

電機控制電路是硬件部分的核心，完成對電機的轉(zhuǎn)動控制，從而對探頭的掃描進行控制，此外，控制電路還負責(zé)與上位機進行通訊，為上位機提供電機轉(zhuǎn)動的狀態(tài)。本系統(tǒng)使用的電機采用步進電機，因此采用脈沖方式實現(xiàn)電機的轉(zhuǎn)動控制，通過控制發(fā)送脈沖的個數(shù)來控制電機旋轉(zhuǎn)的角度，通過改變脈沖的頻率來控制電機的速度和加速度。該硬件控制電路主要由電源、單片機系統(tǒng)、光電隔離、電機驅(qū)動電路、過流保護電路等幾部分構(gòu)成。

電機驅(qū)動電路采用A3977 步進電動機驅(qū)動器。此種驅(qū)動器以DMOS 為基礎(chǔ)的PWM控制全橋驅(qū)動器，內(nèi)置譯碼器，不需要外接微處理器，可操作雙極性步進電機的整步、半步、1/4 步和1/8步模式。

超聲機和PC機的圖像采集卡通過視頻傳輸線連接，通過編程控制圖像采集卡獲得超聲機的圖像并在PC機工作站上顯示。PC機和控制板通過串口線通訊。給控制板上電，在串口通訊助手將指令通過串口線傳送給單片機，通過驅(qū)動部分將驅(qū)動信號發(fā)送給步進電機，實現(xiàn)二維超聲探頭的數(shù)據(jù)掃描和采集。三維超聲圖像的采集和重建流程如圖9所示。首先初始化圖像采集卡參數(shù)，然后發(fā)送掃描開始信號給串口?？刂齐娐方邮盏介_始指令后，電機轉(zhuǎn)動帶動探頭掃描?？刂齐娐返挠行呙栝_始信號，通過串口回傳給工作站程序，工作站程序控制數(shù)據(jù)采集卡開始進行循環(huán)采集數(shù)據(jù)到內(nèi)存，直到從串口接收到有效采集停止信號為止。記錄下采集的二維圖像數(shù)據(jù)序列的幀數(shù)并置于內(nèi)存中，用于三維重建。本系統(tǒng)的三維重建采用本實驗室自主開發(fā)的三維重建庫和DLL完成的。

圖9 二維超聲數(shù)據(jù)的采集和重建流程圖Fig.9 The flow chart of 2D data’s acquisition and reconstruction

2 實驗結(jié)果

本文中，應(yīng)用基于后端掃描的三維超聲成像系統(tǒng)，分別以胎兒體模和真人為實驗對象，做掃描實驗，采集二維數(shù)據(jù)并進行三維重建。步進電機開始和結(jié)束工作時，因電機轉(zhuǎn)速等原因工作不穩(wěn)定，因此對此階段超聲探頭掃描得到的二維數(shù)據(jù)不進行采集。經(jīng)計算，得到的有效掃描角度為58.8o，所以三維圖像的重建角度亦定為58.8o。

圖10體模中胎兒的大小為36周，圖10a所示為胎兒面部三維重建后的結(jié)果，圖10b為對胎兒整體掃描得到的重建結(jié)果。

圖10 以胎兒體模為對象做掃描實驗的重建結(jié)果Fig.10 The reconstruction results of fetus

圖11 以人體腎臟和肝臟為對象做掃描實驗的重建結(jié)果Fig.11 The results of human’s liver and kidney

實驗對象為一成年男子，身體健康，年齡28周歲，身高170 cm，體重65 Kg。圖11的(a)、(b)、(c)為該男子左腎的掃描重建圖像，圖11 d為肝臟掃描重建結(jié)果。掃描時，二維超聲探頭前端固定于人體兩肋骨之間。

3 結(jié)論和展望

本文從機械掃描裝置的設(shè)計到步進電機的控制，搭建了一套基于前端固定后端掃描的掃描模式的三維超聲成像系統(tǒng)，先后對胎兒體模和成人做了掃描重建實驗，成功的采集了一系列的二維超聲圖像數(shù)據(jù)，并各自實現(xiàn)了三維重建和顯示。

根據(jù)二維數(shù)據(jù)的重建結(jié)果可以看出，實驗結(jié)果與實驗對象的現(xiàn)實形狀參數(shù)尚有一定的差異（主要是重建面被拉伸變形），這與數(shù)據(jù)采集過程的穩(wěn)定性和精確性有關(guān)。數(shù)據(jù)采集過程中，需要探頭支架保持平穩(wěn)不動，但操作人員人很難做到理想的相對靜止，此外，采集對象身體的運動（如呼吸）等也會影響數(shù)據(jù)采集的精確性。

為了完善此三維超聲成像系統(tǒng)，得到更精準(zhǔn)的三維超聲圖像，后續(xù)工作將從提高電機控制穩(wěn)定性，減小機械掃描裝置體積，增大二維超聲探頭的掃描范圍，以及三維顯示的參數(shù)設(shè)置優(yōu)化等方面進行完善和改進。

[1]Detmer, Lipscomb K, Blomqvust CG, et al.3D Ultrasound Image Feature Localization Based on Magnetic Scanhead Tracking In vitro Calibration and Validation[J].Ultrasound in Med & Biol,1994, 20(9): 923-924.

[2]郝曉輝, 高上凱, 高小榕, 等.一種基于不規(guī)則采樣平面的超聲三維圖像重建方法, 清華大學(xué)學(xué)報(自然科學(xué)版), 1999, 39(9): 9-13.

[3]Smith Stephen E, Ramm Vonn.High-speed ultrasound volumetric imaging system-part i: Transducer design and beam steering[J].IEEE Trans of Ultrason Ferrolec Frep contr, 1991, 38(2): 100-115.

[4]高上凱.醫(yī)學(xué)成像系統(tǒng)[M].北京: 清華大學(xué)出版社, 2000

[5]Lopes R T, Valent e C M , De Jesus E F O, et al.Detection of paraffin deposition inside a draining tubulation by the compton scattering technique[J].Appl Radiat Isot, 1997, 48(10-12) : 1443-1450.

[6]Cao ZJ, Fery E C, T sui B M W.A scatter model parallel and converging beam SPECT based on the KleinNishina formula[J].IEEE Transaction on Nuclear Science, 1994, 41(4) : 1594- 1599.

[7]Hanke R, Bobel F.Determination of material flaw sizeby intensity evaluation of polychromatic Xray transmission[J].NDT & E International, 1992, 25( 2) : 87-88.

[8]Aaron Fenster, Dónal B Downey, H Neale Cardinal.Threedimensional ultrasound imaging[J], Phys Med Biol.46 2001,46:67–99.

[9]D.H.Pretorius, T.R.Nelson.Three-dimensional ultrasound[J].Ultrasound in Obstetrics and Gynecology, 1995, 5: 219-221.

[10]M.S.Kavic, Three-dimensional ultrasound, New technology[J].Surg Endosc, 1996, 10: 74-76.