-
热声成像信噪比往往较低, 为提高信噪比, 通常需要通过对热声信号进行多次取平均. 然而, 通过取平均提高信噪比的方式会降低热声成像时间分辨率, 阻碍快速热声成像技术的发展. 本文提出一种基于多路放大器加法电路的低成本快速热声成像技术, 将超声换能器接收到的热声信号分成四路同时进入四个放大器, 放大后再将该四路信号经过加法电路累加输入到采集系统中, 实现硬件层面的取平均去噪. 仿体实验结果表明: 通过基于多路放大器加法电路的低成本快速热声成像技术, 成像时间分辨率提高了5倍, 信噪比由~6 dB提升到了~12 dB. 本文通过低成本的简单加法电路原理, 既提高了热声成像技术时间分辨率, 又提升了成像信噪比, 为快速热声成像技术的发展提供了一种崭新的技术, 有助于推动热声成像技术的发展和临床应用.
Thermoacoustic imaging (TAI) is an emerging biomedical imaging method in which microwave is used as an excitation source to generate acoustic signals. The TAI possesses the advantages of high contrast of microwave imaging and high resolution of ultrasound imaging, which is also noninvasive. While the signal-to-noise ratio (SNR) of TAI is often very low. It is usually required by averaging the thermoacoustic signal many times to improve the SNR. However, averaging the signal to improve the SNR can significantly reduce the TAI’s time resolution, which hinders the development of rapid TAI. Here in this paper, we propose to reduce the cost and improve the time resolution of TAI based on multi-channel amplifier and additive circuit. The received thermoacoustic signals are divided into 4 channels and then entered into 4 amplifiers simultaneously. After being amplified, the signals are added and collected by the data acquisition system for reconstructing the image. The phantom results indicate that the time resolution of TAI increases 5 times and the SNR rises from 6 dB to 12 dB, with the multi-channel amplifier and additive circuit adopted. The method proposed in this paper is helpful in promoting the development and clinical application of TAI, especially it has a great significance for developing the ultra-fast TAI. -
Keywords:
- thermoacoustic imaging /
- low cost /
- additive circuit /
- time resolution
1. 引 言
早在1880年, Bell 就发现了光声效应; 1981年, Bowen教授[1]基于光声效应首次提出微波热声成像技术. 直到最近几十年, 该技术才被广泛用于生物医学成像[2-8]. 微波热声成像技术以微波源作为激励源, 通过脉冲微波照射生物组织, 生物组织吸收微波能量使自身温度升高; 由于温度升高, 生物组织出现热弹性膨胀, 这种热弹性膨胀会激发产生超声波并向外传播, 超声波携带着生物组织对微波吸收特性的差异, 采集超声波信号进行图像重建, 就可以实现生物组织内部微波吸收差异的热声成像.
热声成像与光声成像技术的主要区别在于激励源不同, 热声成像和光声成像分别通过脉冲微波和激光激发产生超声波. 然而, 由于激光技术和生物医学光子学的发展相对较快; 因此, 热声成像技术发展相对缓慢, 国内外从事热声成像技术的研究组远少于从事光声成像技术研究组. 其中: 国际上具有代表性的研究组包括: 印第安纳大学Kruger研究组[4,9]、佛罗里达大学 Jiang研究组[7,10]、加州理工学院的Wang研究组[3,11]、电子科技大学的赵志钦研究组[12,13]、华南师范大学的邢达研究组[5,14]. 经过这些科研组的不懈努力, 微波热声成像技术在基本理论[15]、实验方案[16,17]、成像算法[18,19]以及应用[18,20]等方面都得到了充分的完善和补充. 目前报道的热声成像系统放大器基底白噪声是热声信号噪声的主要来源, 一般需要多次采集信号, 通过对信号多次取平均的方法降低白噪声以提高信噪比; 但这会降低热声成像技术的时间分辨率, 极大限制了热声成像的时间分辨率.
为提高热声成像的时间分辨率, 本文提出了一种基于多路放大器加法电路的低成本、快速热声成像技术, 将超声换能器接收到的热声信号分成四路同时进入四个放大器, 放大后再将该四路信号经过加法电路累加输入到采集系统中, 实现硬件层面的取平均去噪. 仿体实验结果表明: 通过基于多路放大器加法电路的低成本快速热声成像技术, 成像时间分辨率提高了5倍, 信噪比由~6 dB提升到了~12 dB. 因此, 这种基于多路放大器的去噪方法, 在快速热声成像技术领域具有潜在的应用价值.
2. 理论依据
在数据采集时, 信号经过设备均会引入白噪声. 在热声成像系统中, 如果不考虑探头和放大器等的带通特性, 热声信号可以表示为
st=s0+nt+na+nA, (1) 其中
st 为采集到的热声信号;s0 为原始热声信号;nt ,na 和nA 分别为超声换能器、放大器和采集卡引入的高斯白噪声. 一般来说, 根据高斯白噪声的随机特点, 可以通过平均的方法降低, 所谓平均指的是将同一个信号采集多次, 然后取其平均值. 如果采集N次取平均, 也就是st=N∑i=1(s0+nt,i+na,i+nA,i)N=s0+N∑i=1(nt,i+na,i+nA,i)N, (2) 其中
nt,i ,na,i 和nA,i 分别代表第i次采集到的nt ,na 和nA . 然而, 取平均的方法需要多次采集数据, 会增加热声数据的采集时间, 降低系统时间分辨率.进一步分析
nt ,na 和nA 三种高斯白噪声可知, 由于放大器的放大特点, 其高斯白噪声要明显强于超声换能器和采集卡. 所以对(2)式做修改:st=s0+nt,i+n∑i=1na,iN+nA,i. (3) (3)式表示对放大器上的高斯白噪声进行平均, 忽略超声换能器和采集卡上高斯白噪声的影响. 本文所提基于多路放大器加法电路的快速热声成像技术, 主要是对放大器上的白噪声进行去噪处理. 该方法在一定程度上既可以提高信噪比又可以减少采集时间, 提高热声成像时间分辨率.
3. 实验方案及系统
3.1 实验方案
传统热声成像的数据采集方式如图1(a)所示. 具体流程为: 超声换能器接收的热声信号首先进入放大器, 经放大器放大后被采集系统获取; 由于信号经过设备会引入白噪声, 采集系统得到的热声信号是低信噪比信号. 通常通过多次采集进行信号平均的方法提高信噪比, 但该方法会降低热声成像时间分辨率, 阻碍了快速热声成像技术的发展. 针对该问题, 为快速、甚至超快热声成像技术提供可能性, 本文提出了一种采用多路放大器加法电路来消除白噪声的方法. 如图1(b)所示, 该方法将热声信号分成四路信分别同时进入四个相同设计的放大器, 然后通过一个加法电路将放大后的热声信号累加后再输入到采集系统.
基于此原理, 本文设计了如图2所示的四路放大和加法电路, 热声信号被分成四路, 超声换能器接收的热声信号分别由Input 1—4输入至四个相同参数的放大器中, 放大后的信号由输出端Output 1—4输入至加法电路中进行累加后, 由Output端输出至采集系统.
3.2 实验系统
为了验证该方案的可行性, 搭建了一套可以完成单通道和四通道切换的热声成像实验系统. 如图3所示, 该系统采用中心频率为3 GHz, 脉宽为750 ns, 峰值功率为70 kW的自制微波源用于热声信号激发; 通过一个偶极子天线和耦合剂(变压器油)将微波能量辐射到样品中; 样品吸收微波能量后激励起的热声信号被一个中心频率为2.25 MHz的非聚焦超声换能器(V323, Olympus Inc. USA)接收, 经四个放大倍数为58 dB 的自制放大器(放大器 –3 dB带宽为260 kHz—2.5 MHz, 输入阻抗为860 Ω)放大后通过加法电路累加被一张50 MPS 采样率的采集卡(PCI4732, Vidts Dynamic, CHINA)采集存于计算机中用于后续图像重建. 超声换能器由旋转电机(RSA100, 卓立汉光, 中国)带动旋转, 以2°间隔进行旋转扫描, 转动180次完成360°热声信号采集, 整个数据采集和电机转动由自行编写的Labview (National Instruments, USA)程序控制. 热声图像通过基于Matlab (MathWorks, USA)编写的延迟叠加图像重建算法重建得到[21].
4. 实验结果与讨论
实验对象为浓度2%的NaCl 溶液, 溶液分别装在5个直径为3 mm的薄壁塑料管中作为吸收体. 首先进行如图1(b)所示的采集方案, 热声信号分成四路进入四通道放大器放大, 放大后的信号输入到加法电路进行累加后进行采集, 完成四通道加法电路数据采集后, 保持吸收体不动, 将图1(b)所示采集方式切换至图1(a)完成单通道的采集.
在图4中展示了两种采集方法在平均50次时采集的热声信号. 定性比较两种采集方法, 四通道信号的信噪比要明显优于单通道. 在图5中绘制了四通道加法电路和单通道的成像结果.
图 5 仿体重建热声图 (a) 四通道加法电路不平均; (b) 四通道加法电路平均 25 次; (c) 四通道加法电路平均 50次; (d) 单通道不平均; (e) 单通道平均25次; (f) 单通道平均50次Fig. 5. Recovered TA images from phantom experimental data using the four-channel amplifier and additive circuit without average (a), the four-channel amplifier and additive circuit with 25 times average (b), the four-channel amplifier and additive circuit with 50 times average (c), the single channel without average (d), the single channel with 25 times average (e), and the single channel with 50 times average (f).其中图5(a)—(c)分别为四通道加法电路不平均、25次平均、50次平均的延迟累加重建热声图像; 其中图5(d)—(f)分别为单通道不平均、25次平均、50次平均的延迟累加重建热声图像. 从图5可知, 若不去平均, 单通道几乎无法实现吸收体的成像, 而四通道加法电路则能实现吸收体的重建. 选定图5中左下角吸收体进行定量分析发现: 图5(a)—(c)中该吸收体的归一化信号幅值分别约为0.82, 0.91和1; 图5(d)—(f)中该吸收体的信号幅值分别约为0.6, 0.87和0.95; 可见, 四通道不平均和取25次平均时该吸收体的归一化信号幅值, 分别与单通道取25次和50次平均时相当. 总体而言: 横向比较, 随着平均次数的增多, 图像质量逐渐提高; 纵向比较, 在相同平均次数下, 四通道加法电路重建结果要优于单通道.
最后, 为了定量比较四通道加法电路和单通道成像的差异, 在图6中绘制了依据频域能量谱估计得到的信噪比随平均次数变化的曲线. 从图6可知, 信噪比先随平均次数增加而提高, 且初期提高明显, 平均次数增多以后变缓. 另外, 单通道采集方案的信噪比则低于四通道加法电路系统. 当两种方案平均次数均为50次时, 单通道的信噪比为6.443 dB, 四通道加法电路信噪比为11.99 dB, 四通道加法电路方案比单通道方案信噪比整整提高了约1.9倍; 在微波源重复频率为50 Hz且信噪比为6.443 dB时, 两种方案分别需要采集50次和10次信号进行平均才可以达到同信噪比, 单通道方案和四通道方案采集信号所需要的时间分别为1 s和0.2 s, 四通道方案的时间分辨率整整提高了5倍. 可见搭建的四通道加法电路成像系统达到了预期的目的.
5. 结 论
本文提出了一种基于多路放大器加法电路的低成本快速热声成像技术. 通过仿体实验证明了该技术在相同时间内实现信噪比的提升, 或者得到相同信噪比信号所需时间更短. 而该系统相比原系统仅仅增加了放大器数量和一个加法电路, 其成本相对整个系统来说相对较低. 不过, 虽然理论上可以通过增加放大器数量来达到多次平均去噪; 但在实际应用时需考虑放大器输入阻抗等问题, 否则过多放大器构成的并联电路会降低整个放大器系统的输入阻抗, 从而影响放大器增益. 虽然本文通过增加单个通道的输入阻抗, 使得并联后的输入阻抗得以增加; 以及通过降低电路的物理温度、选择合适的工作频带和选用低热噪声类型的电阻暂时解决了输入电阻增加带来的热噪声问题; 但在未来设计更多通道的加法电路时还需要对这一问题进行深入的理论研究.
相比提高微波源重复频率、多次取平均降噪的方法, 本文所提基于多路放大器加法电路的低成本快速热声成像方法为提高热声成像技术的时间分辨率提供了一种可行的思路, 是将来实时、快速、甚至超快热声成像技术的一种潜在技术手段.
[1] Bowen T 1981 Ultrasonics Symposium Chicago, USA, Oct. 14−16, 1981 p817
[2] Ku G, Wang L V 2000 Med. Phys. 27 1195
Google Scholar
[3] Xu M, Xu Y, Wang L V 2003 IEEE Trans. Biomed. Eng. 50 1086
Google Scholar
[4] Kruger R A, Reinecke D R, Kruger G A 1999 Med. Phys. 26 1832
Google Scholar
[5] Qin H, Cui Y S, Wu Z J, Chen Q, Xing D 2020 IEEE Photonics J. 99 1
Google Scholar
[6] Yao L, Guo G, Jiang H 2010 Med. Phys. 37 3752
Google Scholar
[7] Yuan Z, Jiang H 2007 Med. Phys. 34 538
Google Scholar
[8] Gao F, Zheng Y, Feng X, Ohl C D 2013 Appl. Phys. Lett. 102 063702
Google Scholar
[9] Kruger R A, Kiser W L, Reinecke D R, Kruger G A 2003 Med. Phys. 30 856
Google Scholar
[10] Huang L, Yao L, Liu L, Rong J, Jiang H 2012 Appl. Phys. Lett. 101 244106
Google Scholar
[11] Wang L V, Zhao X, Sun H, Ku G 1999 Rev. Sci. Instrum. 70 3744
Google Scholar
[12] Chen G P, Yu W B, Zhao Z Q, Nie Z P, Liu Q H 2008 J. Electromagnet. Wave. 22 1565
Google Scholar
[13] 陈国平, 赵志钦, 龚伟, 聂在平, 柳清伙 2009 科学通报 54 1786
Google Scholar
Chen G P, Zhao Z Q, Gong W, Nie Z P, Liu Q H 2009 Chin. Sci. Bull. 54 1786
Google Scholar
[14] Qin H, Yang S, Xing D 2012 Appl. Phys. Lett. 100 033701
Google Scholar
[15] Li C, Pramanik M, Ku G, Wang L V 2008 Phys. Rev. E: Stat. Nonlinear Soft Matter Phys. 77 031923
Google Scholar
[16] Qin H, Qin B, Yuan C, Chen Q, Xing D 2020 Theranostics 10 9172
Google Scholar
[17] Razansky D, Kellnberger S, Ntziachristos V 2010 Med. Phys. 37 4602
Google Scholar
[18] Wang X, Qin T, Qin Y, Abdelrahman A H, Witte R S, Xin H 2019 IEEE T. Antenn. Propag. 67 4803
Google Scholar
[19] Zheng Z, Jiang H 2019 Quant. Imag. Med. Surg. 9 625
Google Scholar
[20] Kellnberger S, Omar M, Sergiadis G, Ntziachristos V 2013 Appl. Phys. Lett. 103 153706
Google Scholar
[21] Hoelen C G A, de Mul F F M 2000 Appl. Opt. 39 5872
Google Scholar
-
图 5 仿体重建热声图 (a) 四通道加法电路不平均; (b) 四通道加法电路平均 25 次; (c) 四通道加法电路平均 50次; (d) 单通道不平均; (e) 单通道平均25次; (f) 单通道平均50次
Fig. 5. Recovered TA images from phantom experimental data using the four-channel amplifier and additive circuit without average (a), the four-channel amplifier and additive circuit with 25 times average (b), the four-channel amplifier and additive circuit with 50 times average (c), the single channel without average (d), the single channel with 25 times average (e), and the single channel with 50 times average (f).
-
[1] Bowen T 1981 Ultrasonics Symposium Chicago, USA, Oct. 14−16, 1981 p817
[2] Ku G, Wang L V 2000 Med. Phys. 27 1195
Google Scholar
[3] Xu M, Xu Y, Wang L V 2003 IEEE Trans. Biomed. Eng. 50 1086
Google Scholar
[4] Kruger R A, Reinecke D R, Kruger G A 1999 Med. Phys. 26 1832
Google Scholar
[5] Qin H, Cui Y S, Wu Z J, Chen Q, Xing D 2020 IEEE Photonics J. 99 1
Google Scholar
[6] Yao L, Guo G, Jiang H 2010 Med. Phys. 37 3752
Google Scholar
[7] Yuan Z, Jiang H 2007 Med. Phys. 34 538
Google Scholar
[8] Gao F, Zheng Y, Feng X, Ohl C D 2013 Appl. Phys. Lett. 102 063702
Google Scholar
[9] Kruger R A, Kiser W L, Reinecke D R, Kruger G A 2003 Med. Phys. 30 856
Google Scholar
[10] Huang L, Yao L, Liu L, Rong J, Jiang H 2012 Appl. Phys. Lett. 101 244106
Google Scholar
[11] Wang L V, Zhao X, Sun H, Ku G 1999 Rev. Sci. Instrum. 70 3744
Google Scholar
[12] Chen G P, Yu W B, Zhao Z Q, Nie Z P, Liu Q H 2008 J. Electromagnet. Wave. 22 1565
Google Scholar
[13] 陈国平, 赵志钦, 龚伟, 聂在平, 柳清伙 2009 科学通报 54 1786
Google Scholar
Chen G P, Zhao Z Q, Gong W, Nie Z P, Liu Q H 2009 Chin. Sci. Bull. 54 1786
Google Scholar
[14] Qin H, Yang S, Xing D 2012 Appl. Phys. Lett. 100 033701
Google Scholar
[15] Li C, Pramanik M, Ku G, Wang L V 2008 Phys. Rev. E: Stat. Nonlinear Soft Matter Phys. 77 031923
Google Scholar
[16] Qin H, Qin B, Yuan C, Chen Q, Xing D 2020 Theranostics 10 9172
Google Scholar
[17] Razansky D, Kellnberger S, Ntziachristos V 2010 Med. Phys. 37 4602
Google Scholar
[18] Wang X, Qin T, Qin Y, Abdelrahman A H, Witte R S, Xin H 2019 IEEE T. Antenn. Propag. 67 4803
Google Scholar
[19] Zheng Z, Jiang H 2019 Quant. Imag. Med. Surg. 9 625
Google Scholar
[20] Kellnberger S, Omar M, Sergiadis G, Ntziachristos V 2013 Appl. Phys. Lett. 103 153706
Google Scholar
[21] Hoelen C G A, de Mul F F M 2000 Appl. Opt. 39 5872
Google Scholar
计量
- 文章访问数: 6749
- PDF下载量: 173