高速128通道小动物多光谱光声断层成像系统
李辉;曹宇;刘红波;彭冬;朱玉坤;王坤;田捷
【摘 要】为了实现小动物光声断层信号的高速采集和实时高质量图像的重建,采用了覆盖角度为270°的128阵元弧形聚焦超声换能器、4个32通道的NI公司数据采集模块和可调谐脉冲激光器以及正则化优化的基于模型的光声断层重建算法.结果表明,系统的空间分辨率可以达到180μm;此系统可以在1ms内完成光声断层数据的采集,在40s以内获得高质量的重建图像.该系统可以用于开展小动物在体的多光谱光声断层成像实验研究.%In order to achieve high-speed acquisition and real-time high-quality image reconstruction of photoacoustic tomography signals for small animal, a 128 element arc-shaped focused ultrasound transducer with 270° coverage angle, four 32 channel NI data acquisition modules, a tunable pulsed laser and model-based photoacoustic tomography reconstruction algorithm with regularized optimization were used.The spatial resolution of the system can be up to 180μm.The in vitro and in vivo imaging experiments show that the system can complete photoacoustic tomography data acquisition within 1ms, and obtain high quality reconstructed images within
40s.The system can be used to carry out in vivo multi-spectral photoacoustic tomography of small animals.
【期刊名称】《激光技术》
【年(卷),期】2017(041)005
【总页数】6页(P669-674)
【关键词】生物光学;光声断层成像;多元弧形聚焦换能器;高速多通道;多光谱
【作 者】李辉;曹宇;刘红波;彭冬;朱玉坤;王坤;田捷
【作者单位】哈尔滨理工大学 自动化学院,哈尔滨 150080;中国科学院 自动化研究所 分子影像重点实验室,北京 100190;哈尔滨理工大学 自动化学院,哈尔滨 150080;中国科学院 自动化研究所 分子影像重点实验室,北京 100190;中国科学院 自动化研究所 分子影像重点实验室,北京 100190;中国科学院 自动化研究所 分子影像重点实验室,北京 100190;中国科学院 自动化研究所 分子影像重点实验室,北京 100190;中国科学院 自动化研究所 分子影像重点实验室,北京 100190
【正文语种】中 文
【中图分类】Q631
光声成像是过去20年来新兴的一种生物医学成像模态[1-6]，它是一种光激发的混合成像模式，将光学成像和超声成像的优点结合了起来。一方面，在光声成像中用来重建图像的信号是超声信号，生物组织对超声信号的散射要比光学信号低2～3个数量级，因此相比纯光学成像，光声成像具有更深的成像深度和更高的空间分辨率；另一方面，光声成像根据不同组织对可见光、近红外光的选择性吸收，利用特定波长的激光脉冲对组织进行照射，成像的是在生物组织内被吸收的激光脉冲能量的分布，其一定程度上反映组织的吸收系数分布，这在纯超声成像中是无法做到的，因此相比纯超声成像，光声成像具有更高的光学对比度。
目前光声成像系统逐步向实时[7-9]、多维度[10-11]方向发展，因此采用单通道、单超声换能器系统由于存在系统鲁棒性差，需进行多次调整，单切面扫描时间过长等问题已不能够适应最新的发展趋势，所以必须研发能够高速多通道采集光声信号的高鲁棒性光声成像系统。而国际上已有很多开展在体光声断层成像方面的研究[12-21],这些研究中成像系统的实现方式和目的各有不同，有的采用线性阵列[12,15-17]，有的主要开展离体的成像实验[13-14];参考文
献[18]中主要针对光声显微成像方面;参考文献[19]中采用的是一个垂直放置的弧形的换能器阵列，需要旋转成像对象才能获得3维的图像;参考文献[20]和参考文献[21]中采用弧形换能器阵列开展实时在体实验，但是由于其激光是垂直激发成像对象即正交模式，所以只适合小动物脑部成像，不适合断层成像；参考文献[7]和参考文献[8]中专门设计用于开展临床实验。因此本文中采用背向传输模式，即超声换能器阵列与脉冲激光在同一侧，这样就保证了在成像不同断层时不会对成像对象带来影响，同时利用美国国家仪器公司(National Instruments,NI)的数据采集模块和多通道弧形聚焦换能器来构建多通道采集模块以实现光声断层成像数据的高速采集，最后利用正则化优化的基于模型的重建算法对原始数据进行高质量的图像重建。
1.1 光声成像前向模型
光声信号的产生是由于短时脉冲激光照射生物组织，组织中的吸收体吸收一部分能量使得局部温度升高，导致发生热弹性膨胀，从而产生超声波。在满足热力限制和应力限制的条件下，声压p(r,t)满足的关系[22]:
式中，r为3维空间内的位置坐标;t表示时间;c(r)和ρ(r)分别为组织的声速和密度;Γ为无量纲格鲁内森参量;H(r,t)是热源函数，代表单位时间、单位体积内的热量。假设生物组织的密度是
均匀的，即在成像区域内ρ(r)是常数，并假设声速是均匀的，即c(r)为常数c，另外热源函数还可以表示为H(r,t)=H(r)×H(t)，其中H(r)表示单位体积内沉积的热能量，H(t)表示脉冲激光光强随时间的分布函数。在实际成像过程中脉冲激光脉宽很短，理论上光强函数可假设为一个脉冲函数，即H(t)=δ(t)，因此(1)式可表示为:
(2)式可以等价地表示为一个初始值问题:
初始条件为:
上述初始值问题可以通过求解一个泊松类型的积分[23]而得到一个解析解:
式中，积分的对象是一个半径为=ct的球形表面S′(r,t)，在2维平面即断层平面，所有的光声信号源和测量点位于同一平面，此时积分是沿着半径为=ct的圆周L′(t)开展，忽略(5)式中的常数，光声断层成像的前向模型可表示为:
1.2 基于模型的逆向重建
光声断层重建算法主要分为3种:第1种是基于雷登变换的滤波反投影重建算法;第2种是基于时
间反转方法的重建算法;第3种是基于模型的重建算法。滤波反投影算法虽然实现容易且重建速度很快，但是其重建图像含有条状伪影而影响图像质量。时间反转方法虽然通过反向模拟超声波传播来得到更好质量的重建图像，但是这个过程需要复杂的数值仿真，不适合实时成像的要求。而基于模型的重建算法是在采集的声压信号数据和组织的吸收分布之间建立一种线性映射关系，继而通过最优化方法去最小化采集的声压信号与利用模型计算的声压信号之间的误差。因而具有很强的灵活性，且模型矩阵只与所使用的图像网格和信号采集系统的参量有关，而与实际的成像对象无关。
在基于模型的光声断层重建算法中，第1步是计算模型矩阵，需要用到导数的数值近似表示，因此(6)式可近似表示为:
式中,采用的是导数的中间差分近似。I(r,t)为:
第2步就是计算(8)式，方法有很多[24-25]。最后(8)式和(9)式可以表示为:
式中，p∈Rm为向量化表示的超声换能器阵元采集到的声压信号, x∈Rn为向量化表示的吸收分布，也即初始声压分布，A∈Rm×n是模型矩阵或系统矩阵，表示一个线性算子描述组织的光学吸收分布与换能器探测的声压信号数据之间的关系。
基于模型的重建算法可以分为两类，第1类是通过求解(10)式最小二乘问题的算法，成为朴素算法，其解称为朴素解:
式中，xnaive为最终求解的朴素解。
第2类算法为基于正则化的算法，由于基于模型的光声断层图像重建问题通常是病态的，因此朴素解与精确解之间通常有较大的偏离，为了获得更好的近似解，需要使用各种正则化项，比如Tikhonov正则化、稀疏正则化、全变分正则化等。
2.1 硬件系统
整个系统的硬件组成示意图如图1所示。光声断层成像系统的激发脉冲激光器采用的是德国INNOLAS公司生产的掺钕钇铝石榴石(Nd∶YAG)激光器，该激光器分成两部分:前部分为抽运源，抽运波长为532nm，输出能量为400mJ，重复频率为10Hz，脉冲宽度为4ns～8ns;后部分为光学参量振荡器(optical parametric oscillator,OPO)可调谐激光器，波长调谐范围为680nm～960nm的红光及近红光波段，在750nm时，输出能量超过80mJ。其余还有耦合透镜、耦合光纤和10根光线束。
正则化常数
超声换能器采集硬件局部放大图如图2a所示。超声换能器采用的是法国IMASONIC公司的定制中心频率是5MHz，带宽为80%，覆盖角度为270°的128阵元弧形聚焦换能器，结构图如图2b所示，每个阵元的结构图如图3a所示。聚焦的半径为40mm，阵元的高度为15mm，相邻阵元之间的间隔为0.1mm。在成像水箱中，通过一个具有四肢的夹具把换能器夹住，使之水平放置。而换能器夹具的上下面各有5个可以放置耦合光纤端头的楔形槽，上面5个楔形槽均匀分割整个圆周，下面的5个楔形槽与上面相应位置的楔形槽堆成，且每个楔形槽与水平放置的换能器成一定角度，而每一个耦合光纤端头把耦合进来的激光整形为线光
如图2a所示，成像小鼠通过悬挂臂支撑系统垂直固定，在实际在体的成像试验中，小鼠先涂上超声耦合胶，之后体表面覆盖一层薄膜与水隔离，最后在小鼠嘴部通过麻醉系统设置好麻醉。竖直的悬挂臂被固定在水平横梁上，而水平横梁与机械垂直升降台融为一体通过直线电机如图1所示。这里采用直线电机是为了减少电磁干扰对采集系统的噪声影响。之后机械垂直升降台通过直线电机的上下运动来达到控制成像小鼠采集不同的断层切片。图3b中描绘了此超声换能器仿真的声场图。可以看到,成像平面是一个直径近似为20mm、厚度为0.7mm的圆盘，此圆盘切片保证了系统在进行断层成像时的特异性。
超声换能器采集的128通道超声信号通过基于NI的采集模块进行处理。具体是工控机机箱采用的是NI的配有8个插槽和高带宽背板的PXIe-1082，每个插槽都可以插入PXI Express模块，且含有一个系统定时同步插槽。主机采用的是NI的基于Intel core i7-3610QE处理器的高性能嵌入式控制器PXIe-8135,负责协调整个采集过程，这样就不用在需要额外的计算机。定时和同步模块采用的是NI的PXIe-6674T,用于同步脉冲激光器触发信号和多通道采集模块。采集处理模块采用的是NI FlexRIO硬件，其为NI LabVIEW现场可编程门阵列(field-programmable gate array,FPGA)提供灵活且可定制的I/O。此硬件包含两部分，分别是4个FPGA模块PXIe-7965R和4个提供高性能模拟和数字I/O的适配器模块NI 5752。NI 5752为32通道数字化适配器模块，分辨率为12位，采样频率为50MS/s，另外还含有抗混叠滤波器和可编程时变增益控制，这样就可以在模数转换过程中对微弱的光声信号进行高频滤波和增益放大，以增大信号的信噪比。这两部硬件共同构成一个可重新配置的子系统，通过LabVIEW FPGA软件进行编程。

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。