专利名称：一种激光散斑血流成像分析方法
技术领域：
本发明属于生物组织血流成像方法领域，具体为一种激光散斑血流成 像分析方法，它尤其适合于以高时空分辨率研究生理和病理状态下的生物 组织局部血流分布和血管形态等。
背景技术：
激光散斑血流成像广泛应用于研究生理和疾病状态下生物组织局部二 维血流变化。但现有的激光散斑血流成像分析方法不能同时兼顾图像的时 间分辨率和空间分辨率。如专利文献"利用激光散斑成像仪监测肠系膜上
微循环血流时空响应特性的方法"(其公开号为CN1391869，
公开日为 2002.07.12)，其衬比计算方法见附图1 (a))，需要计算空间邻域内的衬比 值，因而降低了图像的空间分辨率；而"一种经颅脑血流高分辨成像方法 及其装置"(其公开号为CN1792323，
公开日为2006.06.28)，其衬比计算方 法见附图1 (b))，则需要对时间序列上的若干帧原始散斑图像进行统计分 析，时间分辨率较差，不能满足实时成像的要求。
发明内容
本发明的目的在于提供一种激光散斑血流成像分析方法，该方法能够 解决现有激光散斑成像方法不能兼顾空间分辨率和时间分辨率的局限性， 实现了对生物组织血流二维分布的实时高分辨成像。
本发明提供的激光散斑血流成像分析方法，其步骤包括
(1) 将激光光束照射到被测对象上；
(2) 以相同的曝光时间和帧间隔时间连续采集AA帧被测对象反射的激 光散斑图像；每帧图像采集的曝光时间为lms 100ms，帧间隔时间20
41000 ms，采集的帧数7V为,1^V5100;
(3)对采集所得iV帧图像，取出各帧图像中相同位置处对应A^M大
小空间邻域内的像素，组成大小为A^A^x^个像素的像素集，2^V^7，利
用该像素集内各像素的灰度值/p，利用公式(I)计算该空间邻域内的衬比《
S〃
=
、r
其中々为所采集的图像帧数，M代表空间邻域宽度，々代表w帧图像
中同一位置处对应的空间邻域内共7VxMxiV,个像素中第p个像素的灰度值， 7为这7Vxi^xA^个像素灰度的平均值；
(4) 利用所得空间衬比^,计算该空间邻域处对应的生物组织的血流速
度f^',力，/， 分别为该空间邻域位置处对应的横坐标和纵坐标值；
(5) 按步骤(2) (4)遍历图像中所有的像素，获得所有像素对应 的血流速度F(/,力；
(6) 分,别以每个像素对应的血流速度值为灰度，构建二维的血流速度 分布图。
如果需要进行动态血流分布监测，则包括步骤(7)，在每一个需要对 血流分布进行成像的时间点重复(2) (6)的步骤，获得不同时刻的二 维生物组织血流分布图像。
本发明基于激光散斑血流成像时空联合分析，对生物组织血流进行高 时间和空间分辨率成像。与其它现有的激光散斑血流成像方法相比，本发 明所提供的激光散斑血流成像时空联合分析方法的优点在于结合了空间 散斑衬比分才斤方法与时间散斑衬比分析方法的优点，可实现高时间分辨率、 高空间分辨率的激光散斑血流成像，用于对生物组织二维血流分布和血管 形态，及血流动力学变化的实时、动态、高时间、空间分辨率的监测。其 应用范围广泛，可用于研究生理和病理状态下的大鼠、小鼠、兔、猫、猴
5等实验动物，以及人体的皮肤、眼底、脑皮层局部血流分布，以及神经活 动、脑疾病引起的脑皮层血流变化。本发明适用于脑功能成像、神经生理 学、疾病病理学和药效评价的研究。
图1是已有的激光散斑空间衬比成像方法(图1 (a))，激光散斑时间 衬比成像方法(图1 (b))和本发明提出的激光散斑成像的时空联合分析方 法(图1 (c))示意图。
图2是本发明公开的方法与激光散斑空间衬比分析方法在大鼠脑皮层 血流成像中的实验结果的比较。
图3本发明公开的方法与激光散斑时间衬比分析方法在大鼠脑皮层血 流成像中的实验结果的比较。
具体实施例方式
生物组织血流分布图像的重建需要利用在每一个需要测量血流的空间 邻域所采集的若干帧激光散斑图像，对所采集的激光散斑图像序列进行时 间域与空间域上的联合统计特性分析，计算激光散斑图像中每个空间邻域 内对应的时间序列上所有像素光强(即图像灰度)的统计量，以此统计量 反映该像素所对应生物组织处的血流速度；如此遍历图像中所有像素，即 可获得高分辨的二维生物组织血流分布图像。如附图l(c)所示，其计算衬比 值所使用的空间邻域小于空间衬比算法所使用的空间邻域，所需时间序列 图像帧数小于时间衬比算法所需帧数。因此本发明提出的衬比计算方法为 综合考虑统计精度，时间分辨率和空间分辨率的平衡提供了另一种选择。
动物实验
实验对象为Wistar大鼠，固定于工作台3上，以780nm的半导体激光 器为光源，jt直扩束后的激光光束1经第一线偏振片2后斜入射到去皮后 的大鼠头部。大鼠头部反射形成的激光散斑经第二线偏振片4后由带宏视 镜头的电荷耦合器件相机构成的光电成像系统5成像，光学系统放大倍数
6设为0.5倍p曝光时间5ms，帧间隔时间25ms，连续采集40帧激光散斑 图像。利用所采集的前10帧激光散斑图像，按本发明所述方法(见公式(I))
计算激光散斑衬比，再利用公式(n)计算空间邻域处对应的生物组织的血流
速度F(/J)，进一步重建二维脑血流分布图像，其中校正系数c取值为l。
W，, : .、 (II) 其中c为校正系数。
K"力
为与其它方法比较，选取原始散斑图中某一位置处5xl0个像素的血管
区域，分别用空间衬比，时间衬比和本发明所提出的衬比计算方法计算该
区域内的衬比图，并计算衬比的均值和标准偏差。图2为采用不同空间邻
域计算所得空间衬比值与采用本发明提出的衬比计算方法所得衬比值的比
较。可以看出，采用本发明提出的方法得到的衬比值随空间邻域变化趋势
并不明显，这说明即使使用很小的^1WW，仍然能够保证足够统计i
而使用空间衬比方法得到的衬比值则随空间邻域的增大而增大，当使用较
小的空间邻域时，由于统计量太少，其衬比值偏离准确值较大。图3 (a)
中给出了 2x2和3x3的空间邻域下，本发明所提出的方法得到的衬比值与
采用时间衬比计算方法得到的衬比值随统计帧数的变化。可以看到，当统
计帧数较少时，时间衬比计算方法得到的衬比值稍小于时空联合分析方法
得到的衬比值，从而说明时本发明提出的衬比计算方法即使使用较少的帧
数(例如10帧)，仍然可以得到比时间衬比计算方法所的结果更为准确的
衬比值。另外可以看到，采用时本发明提出的衬比计算方法所得结果中， 图像噪声要小于采用时间衬比计算方法所的结果的噪声。因此，通过上述
比较证明了本发明提出的衬比计算方法可以同时兼顾时间分辨率和空间分 辨率，为生物组织血流测量中统计精度，时间衬比，空间衬比三者的平衡 提供了更多的选择。
以上所述为本发明的较佳实施例而己，但本发明不应该局限于该实施 例和附图所公开的内容。所以凡是不脱离本发明所公开的精神下完成的等 效或修改，.都落入本发明保护的范围。
权利要求
1、一种激光散斑血流成像分析方法，其步骤为(1)将激光光束照射到被测对象上；(2)以相同的曝光时间和帧间隔时间连续采集N帧被测对象反射的激光散斑图像；每帧图像采集的曝光时间为1ms～100ms，帧间隔时间20～1000ms，采集的帧数N为1≤N≤100；(3)对采集所得N帧图像，取出各帧图像中相同位置处对应Ns×Ns大小空间邻域内的像素，组成大小为N×Ns×Ns个像素的像素集，2≤Ns≤7，利用该像素集内各像素的灰度值Ip，利用公式(I)计算该空间邻域内的衬比Kst，其中N为所采集的图像帧数，Ns代表空间邻域宽度，Ip代表N帧图像中同一位置处对应的空间邻域内共N×Ns×Ns个像素中第p个像素的灰度值，<overscore>I</overscore>为这N×Ns×Ns个像素灰度的平均值；(4)利用所得空间衬比Kst计算该空间邻域处对应的生物组织的血流速度V(i，j)；(5)按步骤(2)～(4)遍历图像中所有的像素，获得所有像素对应的血流速度V(i，j)，i，j分别为该空间邻域位置处对应的横坐标和纵坐标值；(6)分别以每个像素对应的血流速度值为灰度，构建二维的血流速度分布图。
2、根据权利要求1所述的激光散斑血流成像分析方法，其特征在于 它还包括步骤(7)，在每一个需要对血流分布进行成像的时间点重复步骤 (2) (6)，获得不同时刻的二维生物组织血流分布图像。
3、根据权利要求1或2所述的激光散斑血流成像分析方法，其特征在 于步骤(4)中利用公式(II)计算空间邻域处对应的生物组织的血流速度WO'),其中c为校正系数。
全文摘要
本发明公开了一种激光散斑血流成像分析方法。将激光光束照射到被测对象上，以相同的曝光时间和帧间隔时间连续采集N帧被测对象反射的激光散斑图像；取出各帧图像中相同位置处对应N<sub>s</sub>×N<sub>s</sub>大小空间邻域内的像素，组成像素集，计算该空间邻域内的衬比K<sub>st</sub>，再由衬比与血流速度的关系式计算对应的生物组织的血流速度V(i，j)，以每个像素对应的血流速度值为灰度，构建二维的血流速度分布图。本发明基于激光散斑血流成像时空联合分析，对生物组织血流进行高时间和空间分辨率成像，实现了对生物组织二维血流分布和血管形态，及血流动力学变化的实时、动态、高时间、空间分辨率的监测。本发明适用于脑功能成像、神经生理学、疾病病理学和药效评价的研究。
文档编号A61B5/02GK101485565SQ20091006074
公开日2009年7月22日 申请日期2009年2月13日 优先权日2009年2月13日
发明者李鹏程, 邱建军, 骆清铭 申请人:华中科技大学