一种改进的激光散斑衬比血流成像方法技术

本发明专利技术公开了一种用于血流成像方法，特别是涉及一种改进的激光散斑衬比血流成像方法。其中，用于血流成像的方法包括：对目标进行数据采集；对原始散斑图像插值，再对插值后的图像进行去相关与标准差处理，得到血流图像。上述技术方案，通过去相关算法来降低静态散斑对图像的影响，有效降低因动物呼吸心跳所引起抖动噪声，提升了图像分辨率，尤其是对小血管识别度，可以为将来研究新生血管或增生血管提供有效检测方式。

An improved laser speckle contrast imaging method for blood flow

The invention discloses a blood flow imaging method, in particular to an improved laser speckle contrast imaging method. Among them, the methods used for blood flow imaging include data acquisition of the target, interpolation of the original speckle image, and then the correlation and standard deviation processing of the interpolated images, and the blood flow images are obtained. The above scheme can reduce the impact of static speckle on the image by using the correlation algorithm, effectively reduce the noise caused by the animal breathing and heartbeat, improve the image resolution, especially the small blood vessel recognition, which can provide an effective detection method for the study of new blood vessels or proliferating blood vessels in the future.

【技术实现步骤摘要】
一种改进的激光散斑衬比血流成像方法
本专利技术涉及一种血流成像方法。特别是涉及一种改进的激光散斑衬比血流成像方法。
技术介绍
激光散斑衬比成像技术于上世纪80年代被首次提出，其后逐渐应用于生物医学的血流监测中。激光散斑现象产生的本质是激光器照射粗糙物体表面反射或者透射所形成的干涉现象，其中的散斑颗粒含有被测物体的运动信息，并且散斑场会随着被测生物组织内的运动元素如血流或结构状态变化而产生相应的变化。Briers和Webster提出散斑图像的时间和空间变化的统计特性均包含了被测物质的运动信息，而且由于激光散斑的随机性，故需采用概率统计的方法来研究其运动特性。Bonner于1976年提出用激光散斑衬比度K来表征血流信息，该激光散斑衬比度等于光强标准偏差比光强平均值，Fercher提出LSSCA算法，该算法通过计算空间衬比度Ks来近似估计衬比度K,从而得到激光散斑血流图，LSSCA算法虽然成像速度快但牺牲了图像空间分辨率，为此Duun等提出sLSCA(spatialbasedLaserSpeckleContrastAnalysis)算法，该算法对Nt帧由Ns*Ns的空间窗计算所得到空间衬比图做时间平均，提高了血流图像空间分辨率，Li提出了LSTCA算法，该算法不需要像LSSCA算法使用滑动窗口，具有较高的空间分辨力，但在一定程度上增加了计算量。尽管已经提出很多激光散斑衬比分析的方法，但是这些算法成像原理本质上还是相近的，均是利用概率统计方式获取衬比值。考虑到现有散斑成像技术存在着动态范围较小，并受心跳、呼吸等噪声影响使其信噪比降低等问题，本专利技术提出了一种改进的激光散斑衬比血流成像方法，该方法结合插值、去相关与标准差这三种算法，有效提升散斑血流图像的动态范围，并降低因心跳、呼吸等噪声对其影响。
技术实现思路
本专利技术所要解决的技术问题是，提供一种使组织结构及微血管血流反射和散射信号得到增强，提高了信噪比的激光散斑衬比血流成像方法；本专利技术所采用的技术方案是：一种改进的激光散斑衬比血流成像方法，包括如下步骤：1)进行数据采集，包括：原始信号采集，生成N帧激光散斑图像；2)图像插值；3)相邻两帧散斑图像进行去相关运算得到去相关图像Ri(x,y)：其中Ri(x,y)为去相关图像，Ii(x,y)第i帧插值后的散斑图像，Ii+1(x,y)第i+1帧插值后的散斑图像；4)对N-1帧去相关图像Ri(x,y)进行标准差运算得到σ(x,y)其中Ri(x,y)是第i帧去相关图像，是N-1帧Ri(x,y)图像的平均值,σ(x,y)是生成血流图像。附图说明为了更加清楚地说明本专利技术示例性实施例的技术方案，下面对描述实施例中所需要用到的附图做一简单介绍。显然，所介绍的附图只是本专利技术所要描述的一部分实施例的附图，而不是全部的附图，对于本领域普通技术人员，在不付出创造性劳动的前提下，还可以根据这些附图得到其他的附图。图1是本专利技术一种改进的激光散斑衬比血流成像方法的流程图；图2是本专利技术中全频谱去相关算法计算原理图；图3是本专利技术实施例中的雌性裸鼠耳部改进的激光散斑衬比血流图。具体实施方式下面结合实施例和附图对本专利技术的一种微血管血流成像方法做出详细说明。如图1所示，本专利技术的一种改进的激光散斑衬比血流成像方法，包括如下步骤：1)进行数据采集，包括：原始信号采集，生成激光散斑图像；所述的生成激光散斑图像包括：目标位置重复采集N幅激光散斑图像。2)图像插值；插值算法包括：最邻近元法(nearestinterpolation)、双线性插值法(bilinearinterpolation)、双三次插值等。最邻近元法处理速度快，但可能出现插值后图像灰度不连续，灰度变化大的地方会出现锯齿状。双线性插值法虽不具有像最近邻元法灰度不连续的缺点，但它本身具有低通滤波效果，使得图像整体轮廓出现较为模糊的情况。双三次插法能有效的降低因线性插值而引起的数据失真，且能得到比双线性插值更为平滑的图像边缘，能丰富图像的细节信息。本专利技术中所述的图像插值，是采用双三次插值函数，数学表达式如公式1所示，|x|为到插值点的距离，利用公式1处理后获得插值后图像Ii(x,y)。。3)去相关处理在获得插值后的图像Ii(x,y)之后，对该图像进行去相关运算。去相关算法能在相关运算的基础上，有效的去除静态组织信息，保留动态组织信息，去相关度较大的为动态组织，去相关度较小的为静态组织。本专利技术采用的去相关算法是针对多帧散斑图像对应像素点进行去相关分析，根据像素点间的去相关度大小区分静态组织和动态组织，进而获取组织内血流信息，能有效降低静态散斑对动态信息的影响。经典的去相关算法如公式2所示，其中M*N为窗大小，IA和IB为相邻两帧图像，和分别对应IA、IB光强平均值。经典的去相关算法在处理大数据时，其运行时间较长。Jia等人提出了一种运用于光相干断层扫描微血管成像技术的全频谱去相关算法(Full-spectrumde-correlation)。该算法首先对OCT原始谱域图像数据先进行滤波，再对滤波后的数据进行FFT变换以获得组织OCT灰度图像，再对灰度图像数据进行去相关运算，进而得到组织OCT微血流图像。全频谱去相关算法如公式3所示：其中Ri(x,y)为去相关图像，Ii(x,y)第i帧插值后的散斑图像，Ii+1(x,y)第i+1帧插值后的散斑图像；图2为N帧像素大小为R*C的图像进行去相关运算原理示意图。该算法针对同一时间序列点的多帧图像进行数据处理，其优点是能对同一位置不同图像进行数据处理，能较好的保持原始图像的空间分辨率，有效降低静态散斑影响，相较于经典去相关算法其能极大缩短运算时间。综合考虑OCT成像与激光散斑成像特点的差异，本专利技术提出的算法去相关过程中略去OCT全频谱去相关算法中谱域滤波和FFT变换过程，直接将公式3算法概念引入到激光散斑血流成像分析中，最后得到去相关后的图像Ri(x,y)。4)对N-1帧去相关图像Ri(x,y)进行标准差运算得到σ(x,y)由于激光照射在被测物体表面时，被测物体的流动信息会引起散斑亮暗的变化，在一定条件下，散斑亮暗变化过程是一种各态遍历的随机过程。在同一曝光时间内，被测物运动速度越快，模糊程度越大，散斑的光强分布的标准差越小，流动信息图像显示更暗，因而需要对去相关后的图像Ri(x,y)做标准差运算，标准差计算过程如公式4所示：其中Ri(x,y)是第i帧去相关图像，是N-1帧Ri(x,y)图像的平均值,σ(x,y)是生成血流图像实例：雌性裸鼠的耳部改进的激光散斑衬比血流成像激光散斑成像系统对裸鼠耳朵中部和左下位置分别采集80帧图像。选取裸鼠耳朵中部位置的80帧原始散斑图像，用改进的激光散斑衬比血流成像算法处理，受裸鼠自身呼吸心跳的影响，裸鼠耳朵会出现微弱的震动，导致生物组织产生一个低速移动，而血流相对于组织具有较高速度。速度低时，标准差较高，速度高时，标准差较低，故iLSCA算法得到的血流图像的整体背景偏亮，血流部分偏暗。如图3所示：注意，上述仅为本专利技术的较佳实施例及所运用技术原理。本领域技术人员会理解，本专利技术不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本专利技术的保护范围。因此，虽然通过以上实施例对本专利技术进行了较本文档来自技高网...

【技术保护点】
一种改进的激光散斑衬比血流成像方法，其特征在于，包括如下步骤：1)进行数据采集，包括：原始信号采集，生成激光散斑图像N帧；2)图像插值；3)相邻两帧散斑图像进行去相关运算得到去相关图像Ri(x，y)：

【技术特征摘要】
1.一种改进的激光散斑衬比血流成像方法，其特征在于，包括如下步骤：1)进行数据采集，包括：原始信号采集，生成激光散斑图像N帧；2)图像插值；3)相邻两帧散斑图像进行去相关运算得到去相关图像Ri(x，y)：其中Ri(x，y)为去相关图像，Ii(x，y)第i帧插值后的散斑图像，Ii+1(x，y)第i+1帧插值后的散斑图像；4)对N-1帧去相关图像Ri(x，y)进行标准差运算得到σ(x，y)

【专利技术属性】
技术研发人员：王立伟，李跃杰，李易陆，
申请(专利权)人：中国医学科学院生物医学工程研究所，
类型：发明
国别省市：天津,12