基于DSP的超声编码激励发射研究<?xml:namespace prefix = o ns = "urn:schemas-microsoft-com:office:office" />
蔡欣荣
(1.西安交通大学电子与信息工程学院,西安,710049;2.西安交通大学生命学院,西安,710049)
摘要:与传统的脉冲回波成像技术相比,编码激励技术能够提高超声图像的信噪比及降低脉冲峰值功率。文章首先介绍了GOLAY码互补序列对作为超声编码激励的原理和优点,然后采用TI的TMS320F2812 DSP芯片,把预先编制好的GOLAY码脉冲序列存储在DSP芯片的FLASH中,通过DSP的McBSP多通道缓冲串口把编码脉冲发送给电压放大器,驱动MOS管信号放大后直接加载在超声探头上,回波信号增益采集,最后,对采集到的数据进行脉冲压缩及求和。实验数据表明信号信杂比(SCR)约为31dB,已经满足了医学超声成像的要求。
关键字:Golay互补序列对;DSP ;MOS ; McBSP ;超声换能器
中图分类号: TP274+ .2
Pulse Transmitting In Ultrasound Coded Excitation System Based On DSP
CAI Xinrong1,MIAO Xianglin1,WANG Xiaosen1,LI Peng2, BIAN Zhengzhong2
(1.Xi’an Jiaotong University, Xi’an 710049, China;2.Xi’an Jiaotong University, Xi’an 710049, China)
Abstract: Compare with traditional pulse echo ultrasound imaging,coded excitation can get high SNR of ultrasound image and low the high pulse peak power.This paper first introduce the basic theory and excellence of Golay complementary serial pair,then we use tms320f2812 , A type of Dsp which is produced by TI(Texas instrument),we put these coded excitation which programmed in advance into the flash of DSP,using McBSP(multichannel buffered serial port)extracts these pulse to electric voltage enlarger,drive MOS and get enlarger voltage to probe,echo signals gained and collected,at last ,compress the collected signals and add A echo signal and B’s. The results show that the signal to clutter ratio is about 31dB,which already satisfy medical ultrasound system.
Keywords: Golay complementary serial pair; Digital Signal Processing; MOS; Multichannel buffered serial port;Ultrasound transducer
引言
传统医学超声成像系统通常采用单一脉冲波[1,2],为了获取高的信噪比,需要提高发射的脉冲峰值功率,这样就受到安全诊断阈值的限制及超声换能器非线性制约。因此,传统的医学超声系统存在脉冲峰值功率高、信噪比差、穿透力弱等缺陷。
现代超声医学成像系统采用编码激励脉冲序列来替代单一脉冲作为发射信号[3],降低了发射脉冲的峰值。在接收信号时经过相关解码电路,探测到人体深部的微弱回波信号。选择一组二值自相关性好的编码序列(如GOLAY互补序列对),将其作为超声编码激励成像系统的发射编码,来达到提高图像的信噪比和穿透力,且实现动态超声图像的实时处理。
文章以GOLAY码互补序列对为例研究了基于DSP的超声编码激励发射,实验数据表明信号信杂比SCR≈31dB,已经满足医学超声成像的要求。
1.GOLAY 码互补序列对
1.1 采用GOLAY互补序列对模型
GOLAY互补序列对定义[4]:一对由两种元素构成的等长有限序列,且在任何给定间隔下,一个序列中的相同元素对的个数等于另一个序列中的相异元素对的个数。数学语言描述如下:设有一对长度相同的有限二相序列<?xml:namespace prefix = v ns = "urn:schemas-microsoft-com:vml" />和,其非周期自相关函数分别为:
如果
就称序列A和序列B互补,或称A、B为互补序列。二相互补序列长度为N必须是偶数,且为两个完全平方数之和。
1.2 GOLAY互补序列对自相关
我们要发射的是128位GOLAY码互补序列对,设探头的中心频率为5MHz,则GOLAY码脉宽为100,128位需要。这里构造GOLAY互补对分别为A、B(各自长度为128比特),他们的自相关后数据图表示分别为:
(a)A序列码自相关,(b)B序列码自相关,(c)两个自相关结果求和之后为。
<?xml:namespace prefix = w ns = "urn:schemas-microsoft-com:office:word" />
图1双极Golay互补序列对自相关(a)、(b)及自 相关求和结果(c)波形图
由图1可以看出,采用GOLAY码互补序列对作为发射激励码,自相关之后求和,旁瓣已经完全消失了,他们的各自的回波自相关也具有这样的特性,所以,采用GOLAY码互补序列对,其信号穿透力强,图象性噪比高。
2 系统整体设计
2.1 采用DSP理论模型
TMS320F2812是美国德州仪器公司(TI公司)推出的最先进也是功能最强大的32位定点DSP芯片。该芯片处理能力强,运算速度高,具有丰富的片内外设,如内部看门狗、CAN、MCBSP、SCI、ADC、集成Flash等。该处理器芯片主要用于家电产品、工业控制场合等高性能的应用领域。
本文使用TMS320F2812[5]作为发射编码的主控芯片,使用多通道缓冲串口(McBSP[6])来完成连续128位编码的发射,并且具有发射速率可以随时修改,以适应不同频率的超声探头(一般探头有3.5MHz、5MHz等),编码通过DSP的MDXA口发送给超声脉冲电压放大板,进行电压放大,放大的电压送入探头换能器,进行电声转换,超声在组织内的回波被超声探头换能器接收进行声电转换,之后信号进行放大并验证波形。
2。2 系统硬件设计
系统时钟采用30MHz,因为TMS320F2812可以工作在150MHz,所以要想全频率工作,PLL需要5倍频率,,JP3为McBSP多引脚,其中我们用到了MDXA口,其他引脚在该发射编码程序中暂且不用,SPICLKA,SCITXDA,SPISTEA,MDXA这些引脚还具有跳线设置功能,只有当SCITXDA设置为高电平,系统从FLASH开始执行程序。该芯片内部集成看门狗,方便调试、节省资源。XMP/MC一般拉地,那么该芯片就认为微计算机器状态,无须外扩FLASH和RAM。
连续128位的编码激励信号经过DSP的MDXA口发出之后,送入SN75372[8]非门进行电压提升,提升的电压可以驱动IRFU420[7]NMOS管,在MOS管源极下拉电阻并在源极处取电压送入探头。当开始产生回波信号时,NMOS此时截止,回波信号送入运算放大器AD8048[8]同向放大,放大倍数为 。为了改善放大部分的波形,在运放输出端加小电容滤波,去处尖锐沿,图3。
图2 DSP 硬件电路部分
图3采用IRFU420 NMOS管进行电压放大
2.3 软件系统设计
2.3.1程序设计流程描述
DSP编码发射程序是在CCS2.0环境下实现的。整个系统过程包括程序代码的编写,程序调试以及在线把程序下载到片内Flash中。
CCS(Code Composer Studio)是美国德州仪器公司(TI)开发的一个完整的DSP集成开发环境。由于TI 的DSP使用非常广泛,使得CCS也就成为目前使用最为广泛的DSP开发软件之一。整个程序流程图如图4。
判断选择:奇次发A序列,偶数次发B序列 |
初始化通用GPIO口,设置GPIO8~GPIO13为McBSP口。 |
清除所有中断,初始化系统中断向量表、向量控制寄存器。 |
初始化McBSP口,包括串口控制寄存器、发送控制寄存器、引脚控制寄存器等。设置32位字发送模式,128位需要4个字发送,发送速率为5MHz。 |
初始化32位比特字模式,启动发送控制寄存器,开始连续发送编码,并匹配时序。 |
一次发送128位结束,开始发送连续高电平,经过反向之后使得MOS管截止,这期间开始回波采集和运算放大,并且DSP要喂狗,防止程序跑飞。
|
循环发送,奇偶变换 |
初始化系统控制寄存器,包括看门狗寄存器,系统时钟倍频(30MHz*5倍频),外设时钟(McBSP为低速时钟)。 |
图4 程序设计流程图
2.3.2 回波分析与采集
在超声成像系统中,若不考虑超声在介质中的衰减时,信号通道的框图如图6所示[9]。是激励信号,是发射换能器的传递函数,是声场中的反射函数,是接收换能器的传递函数,是接收电路的电子噪声。、、、和分别是他们的傅立叶变换[9]。
激励信号 e(t), E(f)
|
发射换能器的 传递函数 p1(t),P1(f) |
发射函数 u(t),U(f) |
接收换能器的传递函数 p2(t),P2(f) |
电子噪声 n(t), N(f) |
超声回波 r(t),R(f) |
图5超声成像中信号通道框图
所以超声回波可以表示为:
医学超声成像系统的研究中,通常将发射换能器和接收换能器的传递函数视为一致的,即:为了简化分析,只研究声场中单一散射子的反射,不考虑接收电路电子噪声的影响,即: 。因此,(3-2-12)式简化为:
图6 128位GOLAY码互补序列对的各自回波压力场(a),(b),两者求和压缩脉冲如图(c),压缩脉冲包络(d),其分贝表示的压缩脉冲包络为(e)
从图6可以看出,采用编码激励超声脉冲发射方式其信杂比为31dB,完全满足了医学超声成像的要求。
3 结 论
本文采用GOLAY码互补序列对作为发射编码激励脉冲,使用DSP2812实现编码的发射,使用MOS进行电压放大,回波信号因为很小,进行增益补偿放大。整个过程包括DSP板和发射电压放大以及回波信号放大板,软件是在CCS2.0软件下编制和调试,具体的时序调试需要匹配好,达到正确发射,在工程实际中要注意每次编码发射之后的空闲循环时间和帧频、每帧线数搭配好。采集的回波信号在MATLAB中展开,并在MATLAB中进行编码压缩和求和。
参考文献:
[1] D.Nahamoo,A.C.Kak .Ultrasonic echo imaging with pesdu-random and pulsed sources:acomparative study.Ultrasonic Imaging[J]:1981,3:1-36
[2] 关力勋,雷纪胜.超声诊断仪原理和维护[M] 北京:人民卫生出版社,1982
[3] N.AH.K.Rao.Investigation of pulse compression technique for medical ultrasound:A smulation study. med.biol.Eng.comp.[J].1994,32(2):181-188
[4] M.J.Golay.Complementary series IRE trans Inform Theory[J] 1961 Apr IT-7:82-87
[5] TMS320F2812. Data Sheet. <?xml:namespace prefix = st1 ns = "urn:schemas-microsoft-com:office:smarttags" />Texas Instruments 2001
[6] TMS320F28x Multichannel Buffered Serial Port (McBSP) Peripheral Reference Guide 2002
[7] IRFU420. Data Sheet. International Rectifier
[8] SN75372. Data Sheet. Texas Instruments 1986
[9]AD8048 DataSheet Analog device
[9] 彭旗宇.数字化医学超声成像方法的研究(清华大学博士学位论文)[D] 清华大学,2003
用户205677 2009-7-16 22:27
用户1405008 2007-12-3 22:59
请教楼上GOLAY编码具体是怎么实现的不可能是用举例的方法把0,1电平一个一个的输出把,看了生成GOLAY序列的公式 An(k),Bn(k)为两序列,wn为系数,Dn迭代延迟,
A0(k) B0(k)为冲击序列
An(k) = An-1(k) + wnBn-1(k-Dn)
Bn(k) = Bn-1(k) - wnBn-1(k-Dn)
不大明白公式 楼主如果理解解释一下谢谢