CN101485565B

CN101485565B - 一种激光散斑血流成像分析方法

Info

Publication number: CN101485565B
Application number: CN2009100607475A
Authority: CN
Inventors: 李鹏程; 骆清铭; 邱建军
Original assignee: Huazhong University of Science and Technology
Current assignee: Huazhong University of Science and Technology
Priority date: 2009-02-13
Filing date: 2009-02-13
Publication date: 2011-04-13
Anticipated expiration: 2029-02-13
Also published as: CN101485565A

Abstract

本发明公开了一种激光散斑血流成像分析方法。将激光光束照射到被测对象上，以相同的曝光时间和帧间隔时间连续采集N帧被测对象反射的激光散斑图像；取出各帧图像中相同位置处对应N_s×N_s大小空间邻域内的像素，组成像素集，计算该空间邻域内的衬比K_st，再由衬比与血流速度的关系式计算对应的生物组织的血流速度V(i，j)，以每个像素对应的血流速度值为灰度，构建二维的血流速度分布图。本发明基于激光散斑血流成像时空联合分析，对生物组织血流进行高时间和空间分辨率成像，实现了对生物组织二维血流分布和血管形态，及血流动力学变化的实时、动态、高时间、空间分辨率的监测。本发明适用于脑功能成像、神经生理学、疾病病理学和药效评价的研究。

Description

一种激光散斑血流成像分析方法

技术领域

本发明属于生物组织血流成像方法领域，具体为一种激光散斑血流成像分析方法，它尤其适合于以高时空分辨率研究生理和病理状态下的生物组织局部血流分布和血管形态等。

背景技术

激光散斑血流成像广泛应用于研究生理和疾病状态下生物组织局部二维血流变化。但现有的激光散斑血流成像分析方法不能同时兼顾图像的时间分辨率和空间分辨率。如专利文献“利用激光散斑成像仪监测肠系膜上微循环血流时空响应特性的方法”(其公开号为CN1391869，公开日为2002.07.12)，其衬比计算方法见附图1(a))，需要计算空间邻域内的衬比值，因而降低了图像的空间分辨率；而“一种经颅脑血流高分辨成像方法及其装置”(其公开号为CN1792323，公开日为2006.06.28)，其衬比计算方法见附图1(b))，则需要对时间序列上的若干帧原始散斑图像进行统计分析，时间分辨率较差，不能满足实时成像的要求。

发明内容

本发明的目的在于提供一种激光散斑血流成像分析方法，该方法能够解决现有激光散斑成像方法不能兼顾空间分辨率和时间分辨率的局限性，实现了对生物组织血流二维分布的实时高分辨成像。

本发明提供的激光散斑血流成像分析方法，其步骤包括：

(1)将激光光束照射到被测对象上；

(2)以相同的曝光时间和帧间隔时间连续采集N帧被测对象反射的激光散斑图像；每帧图像采集的曝光时间为1ms～100ms，帧间隔时间20～1000ms，采集的帧数N为：1≤N≤100；

(3)对采集所得N帧图像，取出各帧图像中相同位置处对应N_s×N_s大小空间邻域内的像素，组成大小为N×N_s×N_s个像素的像素集，2≤N_s≤7，利用该像素集内各像素的灰度值I_p，利用公式(I)计算该空间邻域内的衬比K_st，

K_{st} = \sqrt{\frac{1}{N \cdot N_{s} \cdot N_{s} - 1} Σ_{p = 1}^{N \cdot N_{s} \cdot N_{s}} {(I_{p} - \overset{&OverBar;}{I})}^{2}} / \overset{&OverBar;}{I} - - - (I)

其中N为所采集的图像帧数，N_s代表空间邻域宽度，I_p代表N帧图像中同一位置处对应的空间邻域内共N×N_s×N_s个像素中第p个像素的灰度值，I为这N×N_s×N_s个像素灰度的平均值；

(4)利用所得空间衬比K_st计算该空间邻域处对应的生物组织的血流速度V(i，j)，i，j分别为该空间邻域位置处对应的横坐标和纵坐标值；

(5)按步骤(2)～(4)遍历图像中所有的像素，获得所有像素对应的血流速度V(i，j)；

(6)分别以每个像素对应的血流速度值为灰度，构建二维的血流速度分布图。

如果需要进行动态血流分布监测，则包括步骤(7)，在每一个需要对血流分布进行成像的时间点重复(2)～(6)的步骤，获得不同时刻的二维生物组织血流分布图像。

本发明基于激光散斑血流成像时空联合分析，对生物组织血流进行高时间和空间分辨率成像。与其它现有的激光散斑血流成像方法相比，本发明所提供的激光散斑血流成像时空联合分析方法的优点在于：结合了空间散斑衬比分析方法与时间散斑衬比分析方法的优点，可实现高时间分辨率、高空间分辨率的激光散斑血流成像，用于对生物组织二维血流分布和血管形态，及血流动力学变化的实时、动态、高时间、空间分辨率的监测。其应用范围广泛，可用于研究生理和病理状态下的大鼠、小鼠、兔、猫、猴等实验动物，以及人体的皮肤、眼底、脑皮层局部血流分布，以及神经活动、脑疾病引起的脑皮层血流变化。本发明适用于脑功能成像、神经生理学、疾病病理学和药效评价的研究。

附图说明

图1是已有的激光散斑空间衬比成像方法(图1(a))，激光散斑时间衬比成像方法(图1(b))和本发明提出的激光散斑成像的时空联合分析方法(图1(c))示意图。

图2是本发明公开的方法与激光散斑空间衬比分析方法在大鼠脑皮层血流成像中的实验结果的比较。

图3本发明公开的方法与激光散斑时间衬比分析方法在大鼠脑皮层血流成像中的实验结果的比较。

具体实施方式

生物组织血流分布图像的重建需要利用在每一个需要测量血流的空间邻域所采集的若干帧激光散斑图像，对所采集的激光散斑图像序列进行时间域与空间域上的联合统计特性分析，计算激光散斑图像中每个空间邻域内对应的时间序列上所有像素光强(即图像灰度)的统计量，以此统计量反映该像素所对应生物组织处的血流速度；如此遍历图像中所有像素，即可获得高分辨的二维生物组织血流分布图像。如附图1(c)所示，其计算衬比值所使用的空间邻域小于空间衬比算法所使用的空间邻域，所需时间序列图像帧数小于时间衬比算法所需帧数。因此本发明提出的衬比计算方法为综合考虑统计精度，时间分辨率和空间分辨率的平衡提供了另一种选择。

动物实验：

实验对象为Wistar大鼠，固定于工作台3上，以780nm的半导体激光器为光源，准直扩束后的激光光束1经第一线偏振片2后斜入射到去皮后的大鼠头部。大鼠头部反射形成的激光散斑经第二线偏振片4后由带宏视镜头的电荷耦合器件相机构成的光电成像系统5成像，光学系统放大倍数设为0.5倍。曝光时间5ms，帧间隔时间25ms，连续采集40帧激光散斑图像。利用所采集的前10帧激光散斑图像，按本发明所述方法(见公式(I))计算激光散斑衬比，再利用公式(II)计算空间邻域处对应的生物组织的血流速度V(i，j)，进一步重建二维脑血流分布图像，其中校正系数c取值为1。

V (i, j) = \frac{c}{K_{st}^{2} (i, j)}

(II)其中c为校正系数。

为与其它方法比较，选取原始散斑图中某一位置处5×10个像素的血管区域，分别用空间衬比，时间衬比和本发明所提出的衬比计算方法计算该区域内的衬比图，并计算衬比的均值和标准偏差。图2为采用不同空间邻域计算所得空间衬比值与采用本发明提出的衬比计算方法所得衬比值的比较。可以看出，采用本发明提出的方法得到的衬比值随空间邻域变化趋势并不明显，这说明即使使用很小的空间邻域，仍然能够保证足够统计量；而使用空间衬比方法得到的衬比值则随空间邻域的增大而增大，当使用较小的空间邻域时，由于统计量太少，其衬比值偏离准确值较大。图3(a)中给出了2×2和3×3的空间邻域下，本发明所提出的方法得到的衬比值与采用时间衬比计算方法得到的衬比值随统计帧数的变化。可以看到，当统计帧数较少时，时间衬比计算方法得到的衬比值稍小于时空联合分析方法得到的衬比值，从而说明时本发明提出的衬比计算方法即使使用较少的帧数(例如10帧)，仍然可以得到比时间衬比计算方法所的结果更为准确的衬比值。另外可以看到，采用时本发明提出的衬比计算方法所得结果中，图像噪声要小于采用时间衬比计算方法所的结果的噪声。因此，通过上述比较证明了本发明提出的衬比计算方法可以同时兼顾时间分辨率和空间分辨率，为生物组织血流测量中统计精度，时间衬比，空间衬比三者的平衡提供了更多的选择。

以上所述为本发明的较佳实施例而已，但本发明不应该局限于该实施例和附图所公开的内容。所以凡是不脱离本发明所公开的精神下完成的等效或修改，都落入本发明保护的范围。

Claims

1.一种激光散斑血流成像分析方法，其步骤为：

(1)将激光光束照射到被测对象上；

K_{st} = \sqrt{\frac{1}{N \cdot N_{s} \cdot N_{s} - 1} Σ_{p = 1}^{N \cdot N_{s} \cdot N_{s}} {(I_{p} - \overset{&OverBar;}{I})}^{2}} / \overset{&OverBar;}{I} - - - (I)

(4)利用所得空间衬比K_st计算该空间邻域处对应的生物组织的血流速度V(i，j)；

(5)按步骤(2)～(4)遍历图像中所有的像素，获得所有像素对应的血流速度V(i，j)，i，j分别为该空间邻域位置处对应的横坐标和纵坐标值；

2.根据权利要求1所述的激光散斑血流成像分析方法，其特征在于：它还包括步骤(7)，在每一个需要对血流分布进行成像的时间点重复步骤(2)～(6)，获得不同时刻的二维生物组织血流分布图像。

3.根据权利要求1或2所述的激光散斑血流成像分析方法，其特征在于：步骤(4)中利用公式(II)计算空间邻域处对应的生物组织的血流速度V(i，j)，

V (i, j) = \frac{c}{K_{st}^{2} (i, j)} - - - (II)

其中c为校正系数。