胡勝男，鄭 政

(上海理工大學(xué) 醫(yī)療器械與食品學(xué)院，上海200093)

一種高頻超聲射頻圖像實時處理系統(tǒng)

胡勝男，鄭 政

(上海理工大學(xué) 醫(yī)療器械與食品學(xué)院，上海200093)

為了驗證新的圖像處理算法在實際系統(tǒng)中的處理效果和能力，需要有一個通用的處理平臺。利用FPGA設(shè)計了一種通用高頻超聲射頻圖像實時處理系統(tǒng)。系統(tǒng)同時集成了一片大規(guī)模FPGA，一片250 MHz/14 bit的ADC芯片，一片275 MHz/14 bit的DAC芯片，4片250 MHz，9 MB容量的SSRAM以及USB2.0接口和通用接口等模塊電路。硬件電路完成后，被成功應(yīng)用于常用的高頻超聲射頻圖像處理方法實驗，在系統(tǒng)上運行了帶通濾波、IQ分解、幅度檢測、對數(shù)壓縮、FFT變換、IFFT變換等。

高頻超聲；射頻圖像；實時處理

近年來許多新的超聲圖像處理算法[1-7]被提出，用以提高超聲的成像質(zhì)量和獲取新的信息。比如利用射頻信號中的相位信息進(jìn)行編碼激勵和脈沖壓縮以增強信號[1]；分析超聲圖像斑點從而獲取組織信息[2]；對射頻信號進(jìn)行相關(guān)處理估計組織彈性，提取血流信息[3-4]、諧波成像[5]，對射頻信號進(jìn)行的基于體模點擴散函數(shù)的均衡化斑點噪聲抑制算法[6]；利用超聲導(dǎo)波波譜法來定位和定量管道缺陷[7]等。為了評價這些算法在實際應(yīng)用中的效果，許多實驗室研發(fā)了自己的實驗平臺[8-11]，這些平臺功能各異，但或多或少都存在一些問題。

以此為出發(fā)點設(shè)計了一種通用的高頻超聲射頻圖像實時處理系統(tǒng)，該系統(tǒng)滿足高頻超聲射頻信號的采集、圖像信號處理、數(shù)據(jù)傳輸、數(shù)據(jù)檢驗等要求，并被成功應(yīng)用于常用高頻超聲射頻圖像處理方法的實時實驗。

1 系統(tǒng)設(shè)計

一個通用的高頻超聲射頻圖像實時處理系統(tǒng)應(yīng)該具備3個基本功能，能夠采集超聲掃描信號，提供可編程的器件運行實驗算法，具備檢驗和調(diào)試手段。根據(jù)這些功能，系統(tǒng)主要分為6個模塊：ADC模塊、DAC模塊、FPGA模塊、數(shù)據(jù)存儲模塊、接口模塊和電源?？臁?/p>

1.1 ADC模塊

高頻超聲的頻率在10～50 MHz，根據(jù)奈奎斯特采樣定律，采集模塊的采樣速度必須能達(dá)到100 MHz以上。同時，考慮到精度，避免數(shù)據(jù)位數(shù)不足導(dǎo)致的截斷誤差，設(shè)計選擇250 MHz/14位模數(shù)轉(zhuǎn)換器ADS4149作為ADC模塊電路的核心。

系統(tǒng)上集成了3片DAC芯片，一片高速DAC和兩片低速DAC。高速DAC選擇275 MHz/14 bit的數(shù)模轉(zhuǎn)換器DAC5672A作為電路核心，用于算法調(diào)試和算法運行過程中射頻信號的檢驗。低速DAC選擇125 MHz/8bit的數(shù)模轉(zhuǎn)換器AD9708作為電路核心，一片用于算法運行過程中視頻信號的檢驗，另一片用于超聲回波信號TGC(時間增益補償)的處理。

1.2 FPGA模塊

選用大規(guī)模的FPGA作為信號處理主芯片來實現(xiàn)超聲信號處理算法的實時實驗，型號為Altera公司的Cyclone V E系列的FPGA芯片5CEFA7F31C8N。該芯片具有豐富的片上資源和IO引腳，能滿足復(fù)雜圖像處理算法的資源要求和速度要求。

1.3 數(shù)據(jù)存儲模塊

2維FFT變換時，需對一維變換后的結(jié)果進(jìn)行存儲，待整幅圖像存滿后再作另一維的FFT。此時，F(xiàn)PGA片內(nèi)資源不能滿足一幅圖像的存儲，故需要外部的存儲器。設(shè)計選用同步靜態(tài)隨機存儲器SSRAM，型號為IS61LPS25636A作為數(shù)據(jù)存儲模塊的核心，該芯片速率為250 Msample·s-1，數(shù)據(jù)位數(shù)為36位，總?cè)萘繛? MB。在系統(tǒng)上同時集成了4片這樣RAM，避免了算法復(fù)雜時存儲器分時復(fù)用的復(fù)雜時序設(shè)計，也保證了足夠的資源和處理速度。

1.4 接口模塊

系統(tǒng)同時集成了USB 2.0和GPIO兩種常用的接口，用于數(shù)據(jù)的傳輸以及與外部設(shè)備的連接和通訊，同時也為平臺未來功能的開發(fā)提供可能。

1.5 電源模塊

系統(tǒng)需要4種不同電壓值的工作電源，DC3.3V、DC2.5V、DC1.8V和DC1.1V，并且DC3.3V、DC1.8V還需要同時對數(shù)字域和模擬域進(jìn)行供電。系統(tǒng)采用電池或者3.3 V適配器供電，得到DC3.3V的供電電源，然后將其通過其他電路得到其他工作電源。

2 實驗結(jié)果

系統(tǒng)硬件電路完成后，被成功應(yīng)用于常用的高頻超聲射頻圖像處理方法實驗。搭建的實驗平臺如圖1所示，在平臺上針對10 MHz B模式超聲探頭進(jìn)行了一系列實驗用以驗證其功能。采集了人眼超聲RF圖像，并用系統(tǒng)預(yù)制的模塊對信號進(jìn)行實時處理，將處理后的結(jié)果讀取輸入至PC進(jìn)行觀察。

圖1 超聲圖像處理實驗平臺

2.1 帶通濾波

ADC將采集的超聲回波射頻圖像信號送入FPGA中，總共采集256條線，每條線上為3 072個點。先對圖像信號進(jìn)行實時的帶通濾波處理，采用20階帶通濾波器。圖2顯示了一次超聲脈沖發(fā)射后，即一條線上采集的回波RF信號以及帶通濾波后的結(jié)果，可見低頻噪聲得到了明顯的抑制；圖3對圖2中的虛線區(qū)域進(jìn)行了放大，可見高頻噪聲也得到了有效的抑制，而真實的回波信號得到了完整保留，即說明了帶通濾波器發(fā)揮了較好的濾波效果。

圖2 回波信號帶通濾波前后

圖3 圖2的局部放大

2.2 IQ分解和幅度檢測

對帶通濾波后的圖像數(shù)據(jù)進(jìn)行IQ分解，移頻到基帶后分成IQ兩路。IQ分解后進(jìn)行降3倍采樣，此時數(shù)據(jù)量為256×1 024，即將每條線上的數(shù)據(jù)點數(shù)由原先的3 072點降為1 024點，為后續(xù)的處理減少了計算量，同時也提高了后續(xù)的處理速度。分別采集I、Q兩路數(shù)據(jù)，以及經(jīng)過幅度檢測后的幅度數(shù)據(jù)，用Matlab程序重現(xiàn)，如圖4所示，圖像大小均為1 024×256像素點。

圖4 I、Q兩路以及幅度數(shù)據(jù)顯示

2.3 對數(shù)壓縮

FPGA中的處理模塊對RF圖像信號進(jìn)行帶通濾波、IQ分解、移頻、幅度提取后，經(jīng)過進(jìn)一步1/2降采樣后再進(jìn)行對數(shù)壓縮，采集實驗結(jié)果，用Matlab程序顯示結(jié)果如圖5所示，圖像大小為512×256像素點。

圖5 硬件對數(shù)壓縮后圖像

2.4 2D-FFT和2D-IFFT

對IQ分解降采樣后的數(shù)據(jù)又進(jìn)行了2維的FFT變換，得到頻率域的數(shù)據(jù)，用Matlab程序重現(xiàn)，畫出超聲圖像的二維頻譜圖，如圖7所示。同時，將數(shù)據(jù)與直接在Matlab中做FFT變換的數(shù)據(jù)進(jìn)行對比，檢查系統(tǒng)硬件上FFT變換是否正確。之后又進(jìn)行了2維FFT逆變換，觀察經(jīng)過此處理后數(shù)據(jù)是否與之前沒作FFT變換的數(shù)據(jù)一致。算法流程如圖6所示，先進(jìn)行1 024點的FFT，再進(jìn)行256點的FFT，由于2維的FFT變換需要將一維處理的數(shù)據(jù)進(jìn)行緩存，待整幅圖像存滿后才能進(jìn)行另一維的FFT變換，故將1 024點FFT變換后的數(shù)據(jù)存入一片外部SSRAM中。后續(xù)的2維FFT逆變換也是如此，將一維FFT逆變換后的數(shù)據(jù)存進(jìn)另一片外部SSRAM中，存滿后讀出進(jìn)行另一維FFT逆變換。在此處理過程中，由于系統(tǒng)上同時集成了4片SSRAM，所以不需分時復(fù)用同一片，降低了時序設(shè)計的復(fù)雜性。

圖7 人眼超聲圖像二維頻譜圖

圖6 2D-FFT和2D-IFFT算法流程圖

實驗證明,在系統(tǒng)硬件上進(jìn)行FFT變換后，數(shù)據(jù)與直接在Matlab上進(jìn)行FFT變換的數(shù)據(jù)一致，同時在系統(tǒng)硬件上經(jīng)過FFT逆變換后數(shù)據(jù)與FFT變換前的信號一致，即證明了系統(tǒng)能夠準(zhǔn)確地進(jìn)行FFT變換和FFT逆變換。

3 結(jié)束語

利用FPGA設(shè)計了一種高頻超聲射頻圖像實時處理系統(tǒng)，用于實時實現(xiàn)和驗證新的圖像處理算法。系統(tǒng)設(shè)計完成后，作了實時高頻超聲圖像處理的驗證，驗證了系統(tǒng)的實用性和通用性。然而系統(tǒng)也存在一定的不足之處，系統(tǒng)最后在14層的PCB板上實現(xiàn)，如果將布局進(jìn)行更加合理的規(guī)劃，板層可以減少，從而降低成本。

[1] Misaridis T,Jensen J A. Use of modulated excitation signals inmedical ultrasound，part I: basic concepts and expected benefits[J].IEEE Transactions on Ultrasonics,Ferroelectrics and Frequency Control,2005,52(2): 177-191.

[2] Damerjian V,Tankyevych O,Souag N,et al. Speckle characterization methods in ultrasound images- A review[J]. IRBM,2014,35(4): 202-213.

[3] Parker K J,Doyley M M,Rubens D J. Corrigendum: Imaging the elastic properties of tissue:the 20 year perspective[J].Physics in Medicine and Biology,2011,56(1):R1-R29.

[4] Brands P J,Hoeks A P G,Ledoux L A F,et al.A radio frequency domain complex cross-correlation model to estimate blood flow velocity and tissue motion by means of ultrasound[J].Ultrasound in Medicine & Biology,1997,23(6):911-920.

[5] Duck F A. Nonlinear acoustics in diagnostic ultrasound[J]. Ultrasound in Medicine & Biology,2002,28(1):1 -18.

[6] 李金冬,鄭政.淺表組織超聲圖像的均衡化處理[J].中國生物醫(yī)學(xué)工程學(xué)報,2013,32(2):57-62．

[7] 陳志偉,陳凱歌,蔣俊俊. 基于超聲導(dǎo)波波譜法管道缺陷檢測仿真研究[J].電子科技,2016,29(9): 48-52.

[8] Rohling R,Fung W,Lajevardi P. PUPIL: programmable ultrasound platform and interface library[M].Berlin:Medical Image Computing and Computer-Assisted Intervention - MICCAI 2003,Springer Berlin Heidelberg,2003.

[9] Tortoli P,Bassi L,Boni E,et al. ULA-OP: An advanced open platform for ultrasound research[J]. IEEE Transactions on Ultrasonics,Ferroelectrics and Frequency Control,2009,56(10): 2207-2216.

[10] Qiu Weibao,Yu Yanyan,Tsang F K,et al.An FPGA-based open platformfor ultrasound biomicroscopy[J]. IEEE Transactions on Ultrasonics,Ferroelectrics and Frequency Control,2012,59(7):1432-1442.

[11] 周崢,鄭政.基于FPGA的高頻超聲成像實驗平臺[J].中國生物醫(yī)學(xué)工程學(xué)報,2015,34(5):581-588.

[12] 馮諾.超聲手冊[M]. 南京:南京大學(xué)出版社,1999.

[13] 王威琪,徐智章.醫(yī)學(xué)超聲工程現(xiàn)狀和生物醫(yī)學(xué)超聲概況[J].聲學(xué)技術(shù),1994,13(1):27-35.

[14] 鄒艷榮,鄭政.一種用于高頻超聲探頭的掃描電機及控制系統(tǒng)[J].電子科技,2016,29(10): 111-114.

[15] 王誠. Altera FPGA CPLD設(shè)計:基礎(chǔ)篇[M].北京:人民郵電出版社,2005.

[16] 馮若.超聲診斷設(shè)備原理與設(shè)計[M].北京:中國醫(yī)藥科技出版社,1993.

Research And Implement of Real-time Processing System for High-frequency Ultrasound Radio-frequency Image

HU Shengnan，ZHENG Zheng

(School of Medical Instrumentation and Food Engineering,University of Shanghai for Science and Technology, Shanghai 200093,China)

In order to validate the effect and ability of new image processing algorithms in the actual system,a general processing platform is needed. A general-purpose,real-time image processing system based on FPGA is designed. The system integrates a large-scale FPGA,a 250 MHz / 14 bit ADC chip,a 275 MHz / 14 bit DAC chip,4 pieces of 250 MHz,9 MB SSRAM,USB2.0 interface and the common interface module circuit,etc. After the hardware circuit of the system is finished,the system has been successfully applied in the high-frequency ultrasound RF processing methods,such as band-pass filter,IQ decomposition,amplitude detection,logarithmic compression,FFT transform,IFFT transform etc.

high-frequency ultrasound; RF image; real-time processing

TN911.73;R318.04

A

1007-7820(2017)11-042-04

2017- 01- 02

胡勝男(1991-)，女，碩士研究生。研究方向：醫(yī)學(xué)超聲成像。鄭政(1961-)，男，博士，研究員。研究反向：超聲影像等。

10.16180/j.cnki.issn1007-7820.2017.11.012