作者：周盛 王晓春 计建军 杨军 王延群 单位：中国医学科学院生物医学工程研究所
1方法
1．1发射编码
发射编码按其内容分为相位编码和频率编码，也可以按照发射次数分为单次发射码和多次发射码，如表1所示。Chirp码脉冲压缩效果受衰减和波束合成等非线性因素的影响相对较小，而且只需要一次发射就可以完成脉冲压缩，受组织运动影响较小，因此是超声成像系统编码信号的理想选择之一。但是由于Chirp码需要多幅值发射电压，对发射电压和硬件的要求较高。因此在数字化高频眼科超声诊断设备中，基于发射编码、系统复杂程度等考虑，一般采用二相编码信号作为发射编码。Golay互补序列虽然可以完全消除旁瓣达到最佳脉冲压缩效果，但两次发射的机制，容易受到组织运动的影响，引入误差。虽然眼睛是基本静止的生物组织，但由于我们所选用的是机械扇形扫描探头，所以还是会在一定程度上影响成像效果。M序列在码长较短时的脉冲压缩效果较Barker码差;又根据文献［10－11］给出的最优编码序列的概念、判断标准和计算方法，通过穷举搜索可以得到不同码长的最优编码序列［8］，Barker码除了码长N=11时不是最优编码序列，其余都是。而且文献［12］中也对13位Barker码应用于高频超声成像中做了一定的实验室研究。综合比较，最终采用13位Barker码作为本项目数字超声眼科诊断系统中编码激励的发射码型。
1．2解码压缩方式
目前医学超声编码激励系统中常用的脉冲压缩方法分为匹配滤波和非匹配滤波两种。对于Golay码和Chirp码，匹配滤波都能取得很好的脉冲压缩效果，并且设计及实现方法相对简单。而对于Barker码，采用匹配滤波脉冲压缩方法具有最低的距离旁瓣，但主瓣幅度也只是旁瓣幅度的N倍，其中N是Barker码的长度。对于最长的13位Barker码，匹配脉冲压缩得到的信噪比增益为11．1dB，而峰值旁瓣水平(peaksidelobelevel，PSL)也只有(22dB，对距离旁瓣的抑制不够好［9］。为此设计出了逆滤波器［13－14］、维纳滤波器［15］、尖峰滤波器［16］以及失配滤波器等非匹配滤波脉冲压缩方法［12，17］，以损失很小的信噪比增益为代价来换取更好的距离旁瓣脉冲压缩效果。对比各种解码压缩方式的优缺点，本研究对匹配滤波与非匹配滤波分别进行了仿真实验，结果表明。失配滤波方法可以达到最好的脉冲压缩效果。通过此方法，可以使平均距离旁瓣水平(integratedsidelobelevel，ISL)和PSL都尽可能的低。
1．3仿真实验
二进制13位Barker码表示为:［+1，+1，+1，+1，+1，－1，－1，+1，+1，－1，+1，－1，+1］，由于采样频率一般选用信号回波频率的8倍，这里用以下序列来构造每周期8个采样点的正弦波，分别用［0，－7，－10，－7，0，7，10，7］表示－1，用［0，7，10，7，0，－7，－10，－7］表示+1，图1～图3为Barker码的Matlab仿真波形。其中，图1为13位Barker码波形。图2为经过匹配滤波后得到的结果，主瓣峰峰值达到了±312，而旁瓣峰峰值为±24，主副比为13。图3为13位Barker码匹配滤波结果的峰值包络波形，可见旁瓣衰减达到了－45dB。13位Barker码匹配滤波的FPGA实现相对简单，需要开辟一个深度13×8，14bit的存储空间进行卷积解码。将解码滤波系数序列中为+1的系数对应的回波数据保持原数，将解码滤波系数序列中为(1的系数对应的回波数据进行补码取反，然后进行对应数据并行累加，即可得到匹配滤波的结果。图4是用Modelsim仿真得到的结果图，上端波形为13位Barker码编码激励信号，下端波形为经过匹配滤波解码压缩后的信号。可见，Barker码匹配滤波解码压缩后存在较大的旁瓣。图5为13位Barker码经过失配滤波后得到的结果，实验中我们设计的失配滤波器的阶数为P=3N=39。主瓣峰峰值达到了±127，而旁瓣峰峰值为±1，主副比为127，是匹配滤波主副比的9．7倍。由此说明失配滤波虽然降低了主瓣峰值，但有效地降低了旁瓣水平，提高了信噪比。图6为13位Barker码失配滤波结果的峰值包络波形，可见旁瓣衰减达到了－60dB。在FPGA的设计中，具体的实现步骤如图7所示。首先将失配滤波器系数进行量化，乘以27=28，并取整存入FPGA的RAM中，然后构建实时延时回波矩阵，将系数分别与回波矩阵各行相乘，随后将矩阵39行相乘的结果累加，并除以量化倍数，即可得到失配滤波解码压缩的结果。图8是用Modelsim仿真的结果图，上端波形为13位Barker码编码激励信号，下端波形为经过失配滤波解码压缩后的信号。可见，13位Barker码失配滤波解码压缩后有效地压缩了旁瓣，在保证分辨率的前提下可以有效地提升信噪比。
1．4系统方案与设计
编码激励系统和传统的脉冲回波成像系统不同之处在于:(1)发射电路采用编码发射激励，激励脉冲的长度远大于单脉冲;(2)接收电路需要对回波信号进行脉冲压缩，脉冲压缩算法的计算量很大。因此，编码发射的实验系统设计有以下特殊要求:(1)发射电路能够产生各种不同的编码激励;(2)该系统成像帧频是10帧/s，每帧由512条线组成，每条线包含768个点，要求FPGA能够实时地对编码激励的回波进行脉冲压缩并且成像。根据以上系统设计要求，图9为系统的整体框架结构图。主要包括上位机主机、硬件电路板、3个超声探头、LCD液晶显示、脚闸、USB外设等。系统中的上位机选用嵌入式的工控主板，这样可以实现仪器的小型化和便携性。上位机主要用来实现对底层数据的接收和处理、以及软件的实时控制。实验系统采用单阵元机械扇形扫描探头，3个超声探头分别为中心频率为10MHz的B超探头、中心频率为20MHz的B超探头和中心频率为10MHz的A超探头，其分别用来作为眼球全景成像、眼部组织高频成像和眼球生物测量等模式的检测探头。LCD液晶显示用来对图像实时成像。双键脚闸实现探头的启动、冻结。USB外设包括存储器、打印机、鼠标和键盘。USB外设的使用，使得整机更加人性化和提高了可操作性。
硬件电路板是整个系统的核心部件，其结构如图10所示。整个硬件电路以FPGA及其电路为核心进行设计。FPGA采用Altera公司的CycloneIII系列EP3C55型产品。机械扇扫探头的偏转控制、图像显示的同步、USB接口逻辑等系统控制信号均由FPGA产生。FPGA产生的4组编码脉冲序列通过MOSFET触发电路将脉冲幅度由3.3V转化为12V，以激励发射电路，如图11所示。发射电路由两组双极性脉冲发射芯片构成，并带有归零功能。编码序列通过发射电路后，成为双极性脉冲序列，激励超声换能器晶片。回波信号经高压隔离电路后进入前置放大与可变增益控制电路。前置放大电路增益范围为10～15dB，FPGA产生TGC曲线数据，控制可变增益放大器，可以使增益变化50dB左右。实验系统时钟为80MHz，高速A/D芯片按照不同模式工作在80MHz和120MHz，采样精度为14bits，其采样模式由FPGA通过SPI接口完成设定。FPGA完成的超声回波数字信号处理包括:高频数据采集与存储、数字滤波、数字检波、对数变换、二次采样等，如图12所示。数字信号处理之后的回波数据通过SRAM缓存后，经USB2．0接口电路上传至计算机，并由人机界面完成整个系统的操作与控制。#p#分页标题#e#
2结果
为了对比13位Barker码脉冲压缩中的匹配滤波和失配滤波在实际应用中的效果，设计了单反射面回波实验。编码激励信号由FPGA产生，经发射电路激励10MHz超声换能器，换能器在编码信号激励下，发射超声波，超声信号经反射面反射驱动换能器产生回波信号。回波信号经放大后通过A/D采样进入FPGA进行数据处理，脉冲压缩解码，并最终通过D/A输出到示波器上显示。图13显示了发射电压为±80V时，13位Barker码的编码激励波形。图14是编码激励的单反射面回波信号。图15和图16显示了13位Barker码的编码激励回波经过匹配滤波和失配滤波后的脉冲压缩回波。
由于该系统中回波信号经A/D采样后为14位输入，经匹配滤波后输出为18位，经失配滤波后输出为22位，而系统所选用的D/A只有10位，图15和图16是示波器实时冻结的回波图像，只能截取适当位数通过DA输出至示波器显示比对，无法如仿真软件做到精确量化比对值。从图15中可见，13位Barker码激励回波经过匹配滤波解码后单反射面回波的峰峰值为0．6V，而且存在一些距离旁瓣，最大旁瓣峰峰值达到了0．15V，在超声回波信号中，强回波旁瓣幅值过高会淹没弱小回波信号，影响信号的动态范围，因此抑制脉压旁瓣有着重要的现实意义。而图16中13位Barker码激励回波经过失配滤波解码后单反射面回波就对旁瓣做了有效的抑制，回波峰峰值为0．72V，而最大距离旁瓣为0．06V，增大了主副比，并且距离旁瓣随着离主瓣距离的增大而明显减小。因此在最终的B超成像中，所选用的是失配滤波的脉冲压缩算法。图17为10MHz探头眼科B超的眼球全景图，其中(a)图为单脉冲发射接收的眼球B超图;(b)图为13位Barker码编码发射且用失配滤波的脉冲压缩方法实现的眼球B超图。对比两图的二维眼球图像以及下端的A超波形，且实验中用分辨力测试线靶进行检测，采用编码发射和解码压缩之后的眼科B超图在保持纵向分辨率为0．2mm和基底噪声的同时，大大提升了信号的探测深度，提高了图像的信噪比。
3讨论和结论
本研究搭建了数字化眼科高频超声成像实验平台，基于Barker码实现了编码激励和解码压缩技术，最终将实验所得眼科超声图像数据，通过USB接口实时的在上位机上处理和显示，并获得了定量的实验数据以及成像结果。实验结果表明，编码激励的回波信噪比明显优于单脉冲激励。采用编码激励成像的图像分辨率与单脉冲激励系统相当，且图像没有出现明显的距离旁瓣，而图像的探测深度和噪声水平明显优于单脉冲激励的图像。编码激励系统还对微弱信号有明显的提升作用，由此可以检测出眼球玻璃体混浊等病变，在眼科诊断中具有重要意义。
本设计初步在系统平台上实现了编码激励技术。但还有较多有待拓展、深入之处，主要包括:对于不同带宽的超声换能器，应对发射编码进行相应的加窗处理，使得加窗后形成的激励信号频带与换能器的频带接近，此时可充分发挥超声换能器的转换能力;在之后的研究中，将应用基础序列调制技术来提高长编码发射信号的质量;在脉冲压缩仿真过程中，强发射界面还存在一定的旁瓣影响，在后续的研究中将借鉴雷达信号处理中的线性规划法、最小二乘法以及迭代算法等旁瓣抑制方法，以达到更好的脉冲压缩效果;为了进一步降低脉冲平均声强(ISPPA．3)和机械指数(MI)等声输出参数，下一步将降低发射电压，并观察对成像质量的影响，以便选择更合适的编码及解码方式。