[发明专利]一种改进的微血管血流成像方法

摘要

本发明公开了一种用于血流成像方法，特别是涉及一种改进的光学或超声微血管血流成像方法。其中，用于微血管血流成像的方法包括：对目标进行数据采集；对同一位置重复扫描图像进行匹配；对匹配后图像进行叠加平均得到增强血流信号图像；对匹配后图像进行改进的时间衬比处理得到动态血流位置信息；利用血流位置信息在增强血流信号图像进行标记，得到血流图像。上述技术方案，通过改进的微血管血流成像算法可以提取组织血流信息，从而无损伤、非接触、无需造影剂即可对组织微血管内血流进行成像。

主权项

一种改进的微血管血流成像方法，其特征在于，包括如下步骤：1)进行数据采集，包括：原始信号采集，生成B‑Scan图像；2)图像配准；3)图像灰度叠加平均，是采用如下公式进行图像灰度叠加平均：Is(x,y,z)=1hΣi=1hI(x,yi,z)]]>式中h表示同一位置重复B‑scan的次数，x、y、z为B‑scan空间像素坐标，表示扫描宽度、长度和深度，I(x，yi，z)表示每个位置第i幅B‑scan灰度图像，Is(x，y，z)表示叠加平均灰度图像，包含静态组织信息和动态血流信息。4)改进的时间衬比处理h幅配准后的图像空间上每一个像素都对应一个灰度变化的时间序列，对这一时间序列进行时间衬比统计，则获得该像素点时间衬比，对空间上每一像素点做相同处理，最后获得一组图像统计出的时间衬比图像It，时间衬比公式如下：It(x,y,z)=1hΣi=1h(I(x,yi,z)-I‾(x,y,z))2I‾(x,y,z)]]>式中It(x，y，z)表示时间衬比灰度图像，I(x，yi，z)表示每个位置的第i幅B‑Scan灰度图像，表示每个位置h幅B‑Scan图像的平均灰度图像，连续的h副图像采用时间衬比处理后得到的血流图像有纹路，匹配后相邻图像间存在微小的位移，平均后的信号强度范围增大，与I(x，yi，z)相减时噪声得到增强，因此利用相邻图像做散斑衬比处理，再叠加平均，相邻图像采集时间间隔短，位移量小，计算过程中没有连续的h副图像处理，纹路噪声得到抑制，改进的时间衬比公式采用相邻图像处理再叠加平均，公式如下：Itc(x,y,z)=1h-1Σi=1h-1|I(x,yi+1,z)-I(x,yi,z)|I(x,yi+1,z)+I(x,yi,z)]]>5)图像标记时间衬比图像Itc中包含静态组织和动态组织信息，静态组信号强度变化小，时间衬比值小，动态组织信号强度变化大，时间衬比值大，所以通过改进的时间衬比图像Itc利用下面的公式，对叠加平灰度图像Is进行处理，降低静态组织信息，增强血流信号，得到一幅二维血流灰度图像Iflow：Iflow(x，y，z)＝Is(x，y，z)*Itc(x，y，z)式中Iflow(x，y，z)表示血流灰度图像，Is(x，y，z)表示叠加平均灰度图像，Itc(x，y，z)表示改进的时间衬比灰度图像。对n个位置的h幅B‑Scan图像采取相同的处理，得到三维血流图像。