全数字B超动态滤波器的设计与实现

周 盛李仙琴王晓春何新喜杨 军计建军王延群*

1(中国医学科学院生物医学工程研究所，天津 300192)2(天津迈达医学科技有限公司，天津 300384)

全数字B超动态滤波器的设计与实现

周 盛1李仙琴1王晓春1何新喜2杨 军1计建军1王延群1*

1(中国医学科学院生物医学工程研究所，天津 300192)2(天津迈达医学科技有限公司，天津 300384)

为在全探测深度内获取最佳分辨力的超声回波图像，本研究基于现场可编程门阵列(FPGA)设计了一种动态滤波器，将全数字B超的探测深度平均分为64段，采用64个32阶带通滤波器与之对应，完成人体超声回波信号的动态滤波处理。所设计滤波器应用于128阵元、3.5 MHz的全数字B超设备中，并与使用常值滤波器的图像进行了分析比较。通过对仿组织超声体模的检测显示，在图像的探测深度、远场分辨力和噪声滤除等方面都有了较为明显的提高。

超声;动态滤波器;现场可编程门阵列(FPGA);有限冲激响应(FIR)

Abstract:This study is aimed to obtain the ultrasonic echo images in all detective depth.A dynamic filter was designed based on field programmable gate array(FPGA).The detective depth of the digital B mode ultrasonography was divided into 64 segments，and sixty four 32-order band-pass filters were adopted to complete the process of the dynamic filtering of human ultrasonic echo signals.Then the design was applied to the digital B mode ultrasonography with 128-elements and the center frequency of 3.5 MHz.With the experiments in tissue mimicking ultrasound phantom，the detective depth，the resolution in far field and the denoising ability were all obviously improved compared with the image created using the fixed filter.

Key words:ultrasonic;dynamic filter;FPGA，FIR

引言

超声图像质量对于疾病诊断有重要的意义。由于人体软组织对超声波的衰减与频率大致呈线性关系，因此在近场应主要获取回波信息中的高频成分以提高图像的分辨力和清晰度;而随着探测深度的增加，高频成分更快地衰减使得超声波逐渐向频带的低端偏移，这时需要提取回波信息中的低频成分以获取更丰富的远场图像信息，实现全探测深度内最佳的成像效果［1-2］。全数字B超系统中动态滤波器用以自动选择以上具有诊断价值的频率信号，滤除体表部分以低频为主的强回声信号和深部以高频为主的干扰。动态滤波器设计的好坏直接关系到全数字B超系统远场和近场成像的分辨率，是整个系统的一个关键组成部分［3］。

在模拟B超中，动态滤波电路的实现一般是利用变容二极管的结电容随其反偏电压的增大而减小，从而改变了选频网络的频率。这种方法虽然比较成熟简单，但选频网络的频谱特性不精确，而且改动麻烦，所以现已不被广泛采用。在全数字B超中，动态滤波器都是基于匹配滤波的概念来实现的。在信号处理进程中，匹配滤波器系数随接收深度的变化而动态改变。数字滤波器选用具有线性相位的有限冲激响应 (finite impulse response，FIR)数字滤波器。而传统的FIR滤波器是采用乘累加的硬件结构来实现的，这种结构实现的滤波器需要大量的硬件资源。本研究采用全并行的分布式算法并加入多级流水线，同时使用分割表技术来减少查找表的资源占用［4］。由于动态滤波器中每个滤波器的系数都不同，采用动态查找表的方式，适时调整滤波器的输入系数。

1 方法

1.1 动态滤波器的设计

为了获得全探测深度内的最佳分辨力的回声图像，需要设计的动态滤波器是一个随着时间(深度)的增加，通频带逐渐向下移动的带通滤波器组，该滤波器组由64个滤波器组成。在滤波器的选择上，采用具有线性相位的FIR型滤波器，摒弃了传统的使用乘法器实现乘累加的硬件描述方法，而是基于分布式算法在现场可编程门阵列(field programmable gate array，FPGA)上实现滤波器的硬件结构。本设计利用全数字B超系统的平台实现，采用的超声探头型号为C3.5-128R50S，为128阵元的凸阵探头，中心频率为3.5 MHz，-6 dB带宽为71%。

1.2 滤波器系数的生成

为实现人体超声回波信号的动态滤波处理，将全数字B超的探测深度平均分为64段，采用64个带通滤波器与之对应。超声回波信号采样频率为40 MHz，滤波器阶数为32，其他基本技术参数如下:

式中:Fc1为通带下截止频率，Fc2为通带上截止频率。i的取值范围为1～64，随着深度的增加i增加，即i=1代表体表处，i=64代表体内最深处。

首先，在Matlab中，通过给定的滤波器阶数和带宽，利用 FDATool工具［5］，选取 Hamming窗函数来生成滤波器系数。然后利用FPGA芯片强大的RAM存储功能［6］，在初始化阶段，将滤波器系数量化为8位补码数，存储到FPGA芯片的RAM中。在系统运行阶段，通过精确的时序控制从RAM中得到滤波器的系数，并输入到FIR型动态滤波器的硬件实现部分。以第1个和第32个带通滤波器为例，图1中的(a)和(b)分别是在Matlab仿真下，获得的两个滤波器特性图。从图中可以看出，随着深度i的增加，滤波器的频谱逐渐向频带的低端移动，因此可以很好地匹配超声回波信号的频率成分，滤波体表部分以低频为主的强回声信号和深部以高频为主的干扰，从而获得全探测深度内最佳分辨力的回声图像。

图1 滤波器特性图。(a)第1个;(b)第32个Fig.1 Diagram of the filter characteristic.(a)the first filter;(b)the thirty-second filter

1.3 动态滤波器的结构

为了实现高速信号处理，采用全并行的分布式算法并加入了多级流水线［7-8］，同时使用分割表［9］技术来减少查找表的资源占用。由于动态滤波器中每个滤波器的系数都不同，采用动态查找表的方式，来适时调整滤波器的输入系数。

根据FIR滤波器的对称性，动态滤波器组的每个32阶滤波器均可以对折成16阶线性FIR滤波器。然后根据滤波器的线性特性和分割表原理，将用来实现分布式算法的查找表以每4个抽头一组，分割成4个小型查找表，即4个4阶 FIR滤波器级联成16阶。这种级联结构相对于全并行分布式算法不仅加快了查表速度，还成倍地节约了资源，很容易实现更高阶扩展，同时滤波器的系数还可以适时调整。算法实现的硬件结构如下图2所示。图中LUDT为动态查询表(look-up dynamic table)，根据系统的需要随时通过C0，C1，C2，C3对系数的取值进行调整，从而实现动态的自适应调整［10］。

图3为在Altera公司的Quartus II开发环境下实现的动态滤波器的顶层模块图，主要由ROM模块、动态滤波器模块和输入输出端口组成。图中所示的ROM配置图中，address［5..0］是地址线，共存储64组系数;q［127..0］是对预先存入的滤波器系数的输出，每组共16个数据，每个数据代表8位有符号系数;clock是读数据同步时钟。工作时，系统利用clock作为时钟控制将64组滤波器系数通过q［127..0］输出给动态滤波器模块。在具体的实现中，滤波器系数的值可以通过ROM表预先定制，也可以通过人机交互的方式实时传输，还可以根据系统的输出端的反馈自动地进行相应地调整。当采用后一种方案实现系数调整时，就实现了自适应滤波功能。然而采用这种方法需要实现FFT(快速傅立叶变换)等运算，计算量较大，而且需要较高性能的FPGA，所以采用预先订制的ROM表来传送滤波器系数。动态滤波器模块中，clk_40M是数据输入同步时钟，reset是复位信号，x_in［12..0］和 y_out［12..0］分别为输入输出数据，coeff［127..0］为动态滤波器的系数输入。在Quartus II开发环境下，采用VHDL语言对以动态分布式算法为核心所设计的系数可编程动态FIR滤波器进行描述、仿真和逻辑综合。

图2 动态分布式算法硬件框图Fig.2 Diagram of the dynamic distributed arithmetic

图3 动态滤波器顶层模块图Fig.3 Diagram of the top module of the dynamic filter

1.4 验证方法

为了验证所设计的64组32阶FIR动态滤波器的性能，将这个动态滤波器逻辑下载到课题组研发的全数字B超成像系统中。该系统由发射模块、接收模块、前端控制器、全数字B型成像处理模块和PC组成。系统内部定义了一条机器总线，所有的模块都连在机器总线上，并通过PCI接口与PC交换数据(PC将控制参数下载到各个模块，并从B型成像处理模块读取成像数据)。超声前端由32个物理通道构成，装配了128阵元的凸阵探头，采用40 MHz同步时钟。

分别采用常值滤波器和本研究设计的动态滤波器的全数字B超成像系统对超声体模(中国科学院声学研究所制造的KS107BD(L)仿组织超声体模)进行图像采集，进行对比。

2 结果

所采集的利用不同滤波器的标准人体体模图如图4所示。可以看出，使用常值滤波器的B超成像系统在近场分辨率达到了要求，但在远场的效果较差(见图4(a))。全探测深度为140 mm，小信号探测深度为100 mm，远场噪声较大，使得体模的深度靶点回波信号被噪声掩盖无法显示;而使用本设计的动态滤波器系统在近场和远场上都获得了较好的成像效果(见图 4(b))。全探测深度≥180 mm，侧向分辨力≤4 mm，轴向分辨力≤1 mm，小信号探测深度为120 mm，且远场噪声得到了抑制。图5为对人体肝组织的实际探测图，近场背膜和远场膈肌清晰可见，肝内组织图像细腻，血管与软组织边缘明显。

图4 所采集利用不同滤波器的标准人体体模图。(a)常值滤波器;(b)所设计的动态滤波器Fig.4 The B mode ultrasonic human model image collected by the different filters.(a)the const filter;(b)the dynamic filter

图5 人体肝组织的超声成像Fig.5 The ultrasonic image of human liver

3 讨论和结论

本系统利用Altera公司Cyclone系列的EP1C12Q240C6N型FPGA芯片，设计并实现了一个FIR型动态滤波器。从使用该动态滤波器的全数字B超成像系统采集的图像中可以看出，在远场和近场都达到了较好的分辨力和信噪比。说明本设计的动态滤波器可以根据探测深度自动选择回声信号中有诊断价值的频率成分，并滤除体表以低频为主的强回声信号和远场以高频为主的干扰，从而提高了近场分辨力和远场信噪比，使回声图像的质量得到了改善。

在本设计中，动态滤波器中心频率的变化函数为线性的经验公式，在实际回波信号中未必能很好的匹配实际回波信号的中心频率，会导致超声成像灵敏度的降低。进一步的研究将考虑采用依据实际回波信号的频谱特性自动调整动态滤波器系数的设计方法，提高信噪比，从而提高图像的灵敏度。

［1］Quistgaard JU.Signal acquisition and processing in medical diagnostic ultrasound ［J］.IEEE Signal Processing Magazine，1997，14:67-73.

［2］冯若，刘忠齐，姚锦钟，等.超声诊断设备原理及设计［M］.北京:中国医药科技出版社，1991.

［3］Christopher JH，James MP，Robert JE，et al.Advance in ultrasound［J］.Clinical Radiology，2002，57:157-177.

［4］Patrick L，Ali M.Area-efficient FIR filter design on FPGA using distributed arithmetic［A］.In:IEEE International Symposium on Signal Processing and Information Technology［C］.Vancouver:IEEE，2006，6:248-252.

［5］飞思科技产品研发中心.MATLAB 7辅助信号处理技术与应用［M］.北京:电子工业出版社，2005.

［6］Girard P，Hbron O，Pravossoudovitch S，et al.Delay fault testing of look-up tables in SRAM based FPGA［J］.Journal of Electronic Testing，2005，21(1):43-55.

［7］EDA先锋工作室.Altera FPGA-CPLD设计［M］.北京:人民邮电出版社，2005.

［8］Abed KH，Vivek V，Nerurkar SB.High speed digital filter design using minimal signed digital representation［A］.In:The IEEE Southeast Conference［C］.Lauderdale:IEEE，2005.105-110.

［9］Uwe MB.Digital Signal Processing with Field Programmable Gate Arrays［M］.Berlin:Springer-Verlag，2002.

［10］Punchalard R，Lorsawatsiri A，Loetwassana W，et al.Direct frequency estimation based adaptive algorithm for a second-order adaptive FIR notch filter［J］.Signal Processing，2008，88(2):315-325.

Design and Realization of a Dynamic Filter in Digital B Mode Ultrasonography

ZHOU Sheng1LI Xian-Qin1WANG Xiao-Chun1HE Xin-Xi2YANG Jun1JI Jian-Jun1WANG Yan-Qun1
*1(Institute of Biomedical Engineering，Chinese Academy of Medical Sciences，Tianjin 300192，China)2(MEDA CO.，LTD，Tianjin 300084，China)

R318

A

0258-8021(2010)03-0418-04

10.3969/j.issn.0258-8021.2010.03.016

2010-01-10，

2010-04-09

中央级公益性科研院所基本科研业务专项资助;天津市科技支撑项目(2007CKFSF01400)

*通讯作者。E-mail:wangyanqun-2008@163.com