_基于稀疏轮廓扫描的滑动焦点光声显微成像方法.docVIP

基于稀疏轮廓扫描的滑动焦点光声显微成像方法 摘要： 在基于光声共焦的超声分辨光声显微成像系统中，只有焦点处的吸收体可以获得高分辨率高信噪比的成像，而实际的生物组织往往具有不规则的轮廓，导致成像结果分辨率和信噪比不均匀。本文提出了一种基于手动描绘轮廓的滑动焦点光声显微成像方法，首先在固定的焦点上稀疏扫描轮廓，结合人工选择轮廓点获取轮廓的稀疏矩阵，然后通过线性拟合获取完整的轮廓矩阵，并将轮廓矩阵的信息转化为焦点位置信息，最后进行二次扫描，根据焦点的位置信息在每次扫描之前调整焦点位置，使得在整个扫描过程中，样品始终处于焦点的位置，从而获得分辨率和信噪比均匀的光声显微图像。通过在体皮下血管成像和皮下肿瘤血管再生成像，证明了该方法可以有效提高光声成像的质量。 关键词： 图像处理；光声成像；滑动焦点 引言 光声显微成像（Photoacoustic Microscopy）是光声成像技术的一种新兴模式，通过使用聚焦型超声传感器探测[1]，或者聚焦光入射[2]的办法，目前已实现亚波长的高分辨成像[3]，在生物医学研究和临床应用等领域展示了巨大的潜力。通过暗场反射式照明配合聚焦型超声探头，实现了轴向15 um、侧向45 um的分辨率以及3 mm的最大成像深度，目前已用于小动物的皮下血管成像[4]、肿瘤及周围血管结构与功能成像[5]、颅下血管成像[6]等领域。相对于激光散斑成像、光学相干成像、共聚焦显微成像等高分辨光学成像技术，光声显微成像可以在更深的深度实现高分辨的内源及外源信号成像。 暗场反射式光声显微成像的系统描述和成像原理可在文献[7]中找到，这里简要介绍一下系统的图像重建，通过时间分辨的超声探测，乘以超声在组织中的传输速度，可以获取x-y平面上一个位置的深度信息（A scan），通过沿x轴扫描，可以获取x-y平面上一条线的截面信息（B scan），通过在x-y平面上二维光栅扫描，可以获取三维的结构信息（C scan），通过三维可视化重建，即可得到组织的三维结构图。除此之外，二维的最大值投影图（MAP）使用也很广泛，由于血液的吸收远大于周围组织，所以血管的光声信号远大于周围组织，从而在某个方向上求取最大值，可以保留血管的信号，当然最常用的还是在深度方向上求取最大值投影，这样不但可以获取x-y平面上的血管分布，而且可以抵消部分运动伪迹。本文展示的结果主要采用B型扫描和沿深度方向最大值投影图的形式。 在暗场反射式光声显微成像系统中，为了实现高分辨的光声显微成像，需要使用高中心频率、宽带宽、高数值孔径的球形聚焦超声传感器，而这种结构设计的劣势就是焦深也随之降低[8]，这导致当成像物体离焦时，成像质量会迅速下降，在传统的扫描方式下，探头在x-y平面上以固定的高度扫描，样品的轮廓波动大于焦深时就会导致部分区域的成像质量下降，在重建图像上直观的表现为分辨率下降导致的血管变粗以及信噪比的下降。一维[8]和二维[9]数值孔径合成算法可以在一定的离焦范围内改善这种恶化，但是目前只有一定的方向上进行合成才能获得比较好的效果，而且这种算法并未考虑组织与水的不同声阻抗参数，故而存在一定的误差[10]。因此有必要设计出一种能够使得探头的焦点沿样品轮廓扫描的方法和系统，保证成像的分辨率和信噪比在整幅图像上一致。 本文设计了一种基于轮廓扫描的滑动焦点光声显微成像方法，并在现有的光声显微成像系统上实现了该方法，通过在成像前稀疏扫描样品轮廓，人工选择部分轮廓点结合线性插值的办法获取样品的轮廓矩阵，并将轮廓矩阵信息转化为焦点的位置移动信息，在成像时每移动一个位置，就相应调整焦点的位置，使得在整个扫描过程中样品始终处于焦深范围内，从而获得均一、高分辨、高信噪比的图像。 方法 2.1 轮廓信息的获取 焦点的滑动是按照轮廓信息来执行的，因此在成像之前需要首先获取轮廓信息，这一过程主要分为以下两步： 2.1.1 稀疏轮廓B型图像的获取 稀疏轮廓B型扫描基于传统的固定焦点扫描，样品放置好之后，先将焦点调整到感兴趣区域，然后扫描多个位置的B型图像，位置的选取要充分考虑样品表面轮廓的波动情况，当波动较大时，可以适当的密集扫描B型图像，当波动较小时，可以适当的稀疏扫描B型图像。为了方便后面的插值，一般要求轮廓B型图像的扫描要包括视场的边界。 2.1.2 轮廓矩阵的获取 获取轮廓稀疏B型扫描图像之后，需要人工介入选择每幅B型图像上的轮廓位置，然后通过biharmonic样条曲面插值获取整个视场中的轮廓点位置（Aij），该插值方法的原理简述如下[11]： 对于二维空间中分布的N个控制点，i=1~N，Biharmonic样条曲面插值可以转化为求取方程组（1）的解： （1） 其中为Biharmonic算子，为单位冲击函数，P为二维空间中待插值的位置，为P处插值的结果。方程组（1）的解为：