基于图像匹配算法提取电动力场中细胞自旋转速度的方法技术

本发明专利技术提出基于图像匹配算法提取电动力场中细胞自旋转速度的方法，包括：用CCD相机获取RGB视频图像；将RGB视频图像灰度化；手动选择感兴趣细胞；跟踪感兴趣细胞；确定参考帧；计算自旋转速度。本发明专利技术除在第一帧中手动选择感兴趣细胞外，其余步骤可使用背景相减法实现对感兴趣细胞的自动跟踪，无需反复观察录像和手动选择参考图像，减少了对实验结果的人为干涉。使用曲线拟合技术，能够准确地确定感兴趣细胞旋转一次所需的时间，提高细胞自旋转速度的精度和稳定性。

【技术实现步骤摘要】
基于图像匹配算法提取电动力场中细胞自旋转速度的方法
本专利技术属于电动力微流控芯片
，具体涉及基于图像匹配算法提取电动力场中细胞自旋转速度的方法。
技术介绍
作为细胞行为的一种，细胞自旋转特性对于探索和阐明细胞的内在特性，如电生理、生物力学和介电性质具有重要意，而这些内在特性又与人类疾病相关，因此研究细胞自旋转特性可以确定细胞的自旋转速度，也可以获取细胞的运动轨迹，还可以提高我们对癌症治疗效果的理解。然而目前对于细胞自旋转速度的提取均缺乏有效的手段，只能粗略地提取细胞的自旋转速度，还需要人工干预。为实现提高确定细胞旋转的精度和自动化的整个计算过程，就需要开发了一种新的基于匹配的检测算法，简化和优化自旋转速的提取过程的新方法。目前，具有基于图像匹配算法准确实时提取电动力场中细胞自旋转速度的方法，尚未报道。
技术实现思路
本专利技术提出基于图像匹配算法提取电动力场中细胞自旋转速度的方法，包括如下步骤：步骤1：用CCD(charge-coupled-device)相机获取RGB视频图像；步骤2：将RGB视频图像灰度化，灰度化公式如公式(1)所示：Gray＝a×R+b×G+c×B(1)其中，Gray是灰度值，R、G、B为获取RBG图像的红色分量值、绿色分量值、蓝色分量值，a、b、c为灰度化参数值；步骤3：在灰度化后第一帧图像中用矩形框手动选择感兴趣细胞，所选择矩形区域被称为感兴趣细胞区域，感兴趣细胞区域可记为Rect[(uR,vR),2W+1,2H+1]表示中心为(uR,vR)、宽度为2W+1、高度为2H+1的矩形区域；步骤4：跟踪感兴趣细胞：包括步骤4.1～步骤4.5：步骤4.1：计算背景图像：用步骤2中灰度化后图像前N帧的均值作为背景图像，如公式(2)所示：其中，fB表示背景图像，fi表示视频的第i帧图像，N是用于计算背景图像的帧数；步骤4.2：获得差分图像：步骤2灰度化后的图像的每一帧图像中减去背景图像，根据公式(3)获得差分图像△fi：△fi＝|fi-fB|(3)其中，fB表示背景图像，fi表示视频的第i帧图像；步骤4.3：二值化差分图像：根据公式(4)将差分图像二值化，其中阈值level是通过Otsu方法确定的，Otsu方法是一种选择阈值以最大化类间方差的经典阈值选择方法；△fi(u,v)为二值化后图像；步骤4.4：计算二值化后图像连通区域中心与最大半径：将二值化差分图像运用形态学开运算去除噪声，运用形态学闭运算连接相邻前景区域，根据公式(5)计算各二值图像各连通区域的中心与最大半径：其中i表示图像编号，上标q表示连通区域编号，表示第i幅图像第q个连通区域的中心，表示第i幅图像第q个连通区域像素的坐标，Miq表示第i幅图像第q个连通区域中的像素数量，Ro(·)和max(·)分别表示四舍五入函数和最大值函数，riq为第i幅图像第q个连通区域的半径，连通区域中心与半径的示例见图4；步骤4.5：寻找感兴趣细胞的运动区域和搜索窗口，具体为：在视频的第一帧中，将最靠近手动选择的矩形区域的中心(uR,vR)的连通区域的中心作为第一个感兴趣细胞的运动区域的中心，以(u1,v1)表示，将该连通区域的半径以r1表示，则以(u1,v1)为中心、宽度和高度均为2r1+1的矩形区域，被认为是第一帧中感兴趣细胞的运动区域，记作Rect[(u1,v1),w,h]；寻找第i帧中与第i-1帧中的感兴趣细胞运动区域中心(ui-1,vi-1)最为接近的连通区域，由其中心和半径可确定第i帧的感兴趣细胞的运动区域，记为Rect[(ui,vi),wi,hi]，检测出所有视频图像的感兴趣细胞运动区域，将Rect[(ui,vi),wi+2W,hi+2H]作为感兴趣的细胞的搜索窗口；步骤5：确定参考帧：采用归一化相关系数Rij(up,vp)表示两个相同大小的图像区域之间的相似性，如公式(6)，除第一帧以外其他图像中感兴趣细胞区域被定义为：与前一帧中的感兴趣细胞区域具有最大相关系数的区域，由公式(8)可确定。其中，m、n为变量，(uRj,vRj)表示第j帧中的感兴趣细胞区域的中心，i∈[1，N]，表示第j帧图像中Rect[(uRj,vRj),2W+1,2H+1]区域的平均灰度值，表示第i帧图像中Rect[(up,vp),2W+1,2H+1]区域的平均灰度值，由公式(7)计算可得。公式(7)表示中心为(up,vp),宽度为2W+1,高度为2H+1区域的平均灰度值，m、n为变量；其中Si代表第i帧图像中感兴趣细胞的运动区域，Rij(up,vp)表示第i帧和第j帧中两个区域之间的相关系数；令j＝i-1，由公式(8)可计算第i帧与第i-1帧图像之间的最大相关系数Cii-1，其中i∈[2，N]，如此便可得到最大相关系数序列[C21,C32,…,CNN-1]，对最大相关系数序列[C21,C32,C43,…,CN-1N-2]和[C32,C43,C54,…,CNN-1]分别按升序进行排列，排序后的序列分别记为FS和BS。将Cii-1在FS中的序号与Ci+1i在BS中的序号之和最小的第i帧图像作为参考帧，用fT表示，相应的，Rect[(uRT,vRT),2W+1,2H+1]表示参考帧中的感兴趣细胞区域；步骤6：计算自旋转速度：将公式(8)中的下标j替换为T，计算出各帧与参考帧之间的最大相关系数CiT，并根据公式(8)确定各帧图像中感兴趣细胞区域的中心(uRi，vRi)，将序列{CiT}中大于指定阈值的局部极大值作为该序列的峰值点，即其中T为阈值，各峰值点及其两相邻点被用来计算抛物线曲线Y＝aX2+bX+c，采用多项式拟合抛物线中参数A、B、C，其中，X表示图像的帧号，Y表示对应的最大相关系数，相邻抛物线的最大点之间的时间间隔即为感兴趣细胞旋转一周所需的时间，即Xi-Xi-1，由公式(9)即可算出细胞的平均自旋转速度：其中，ffps表示视频的帧速率，ffps＝15fps，fps为CCD相机画面每秒传输的帧数，Xi和Xi-1表示相邻最大点的帧号，Xi-Xi-1为相邻抛物线的最大点之间的时间间隔即为感兴趣细胞旋转一周所需的时间，n是平均自旋转速度，单位为：转/分钟。公式(9)推导过程：ffps表示每秒传输视频的帧速率，那么1/ffps为传输一帧视频所需要的时间，为传输(Xi-Xi-1)帧图像所需的时间，单位为秒；n为转速，单位为转/分钟。有益技术效果：本专利技术的优点与积极效果为：(1)本专利技术除在第一帧中手动选择感兴趣细胞外，其余步骤实现细胞的自动跟踪；(2)本专利技术使用背景相减法技术，对感兴趣的细胞进行自动跟踪，大大提高了分析效率；(3)本专利技术无需反复观察录像和手动选择参考图像，减少了对实验结果的人为干涉，实现了更大的稳定性；(4)本专利技术使用曲线拟合技术，能够准确地确定感兴趣细胞旋转一次所需的时间，提高细胞获得自旋转速度的精度。附图说明图1为本专利技术实例的细胞自旋转速度获取系统流程图；图2为本专利技术实例的手动选择第一帧中感兴趣的细胞示意图；图3为本专利技术实例的自动跟踪感兴趣的细胞流程图；图4为本专利技术实例的每个标记连通区域的中心和半径示意图；图5为本专利技术实例的寻找感兴趣细胞的运动区域和搜索窗口的示意图；图6为本专利技术实例的选择参考帧的示意图；图7为本专利技术实例的细胞在不同周期下的本文档来自技高网...

【技术保护点】
1.基于图像匹配算法提取电动力场中细胞自旋转速度的方法，其特征在于，包括如下步骤：步骤1：用CCD相机获取RGB视频图像；步骤2：将RGB视频图像灰度化，灰度化公式如公式(1)所示：Gray＝a×R+b×G+c×B  (1)其中，Gray是灰度值，R、G、B为获取RBG图像的红色分量值、绿色分量值、蓝色分量值，a、b、c为灰度化参数值；步骤3：在灰度化后第一帧图像用矩形框手动选择感兴趣细胞，所选择矩形区域被称为感兴趣细胞区域，感兴趣细胞区域可记为Rect[(uR,vR),2W+1,2H+1]表示中心为(uR,vR)、宽度为2W+1、高度为2H+1的矩形区域；步骤4：跟踪感兴趣细胞：包括步骤4.1～步骤4.5：步骤4.1：计算背景图像：用步骤2中灰度化后图像前N帧的均值作为背景图像，如公式(2)所示：

【技术特征摘要】
1.基于图像匹配算法提取电动力场中细胞自旋转速度的方法，其特征在于，包括如下步骤：步骤1：用CCD相机获取RGB视频图像；步骤2：将RGB视频图像灰度化，灰度化公式如公式(1)所示：Gray＝a×R+b×G+c×B(1)其中，Gray是灰度值，R、G、B为获取RBG图像的红色分量值、绿色分量值、蓝色分量值，a、b、c为灰度化参数值；步骤3：在灰度化后第一帧图像用矩形框手动选择感兴趣细胞，所选择矩形区域被称为感兴趣细胞区域，感兴趣细胞区域可记为Rect[(uR,vR),2W+1,2H+1]表示中心为(uR,vR)、宽度为2W+1、高度为2H+1的矩形区域；步骤4：跟踪感兴趣细胞：包括步骤4.1～步骤4.5：步骤4.1：计算背景图像：用步骤2中灰度化后图像前N帧的均值作为背景图像，如公式(2)所示：其中，fB表示背景图像，fi表示视频的第i帧图像，N是用于计算背景图像的帧数；步骤4.2：获得差分图像：步骤2灰度化后的图像的每一帧图像中减去背景图像，根据公式(3)获得差分图像Δfi：Δfi＝|fi-fB|(3)其中，fB表示背景图像，fi表示视频的第i帧图像；步骤4.3：二值化差分图像：根据公式(4)将差分图像二值化，其中阈值level是通过Otsu方法确定的，Otsu方法是一种选择阈值以最大化类间方差的经典阈值选择方法；Δfi(u,v)为二值化后图像；步骤4.4：计算二值化后图像连通区域中心与最大半径：将二值化差分图像运用形态学开运算去除噪声，运用形态学闭运算连接相邻前景区域，根据公式(5)计算各二值图像各连通区域的中心与最大半径：其中i表示图像编号，上标q表示连通区域编号，表示第i幅图像第q个连通区域的中心，表示第i幅图像第q个连通区域像素的坐标，Miq表示第i幅图像第q个连通区域中的像素数量，Ro(·)和max(·)分别表示四舍五入函数和最大值函数，riq为第i幅图像第q个连通区域的半径；步骤4.5：寻找感兴趣细胞的运动区域和搜索窗口，具体为：在视频的第一帧中，将最靠近手动选择的矩形区域的中心(uR,vR)的连通区域的中心作为第一个感兴趣细胞的运动区域的中心，以(u1,v1)表示，将该连通区域的半径以r1表示，则以(u1,v1)为中心、宽度和高度均为2r1+1的矩形区域，被认为是第一帧中感兴趣细胞的运动区域，记作Rect[(u1,v1),w,h]；寻找第i帧中与第i-1帧中的感兴趣细胞运动区域中心(ui-1,vi-1)最为接近的连通区域，由其中心和半径可确定第i帧的感兴趣细胞的运动区域，记为Re...

【专利技术属性】
技术研发人员：梁文峰，杨谢柳，牛玺辉，邹媛媛，
申请(专利权)人：沈阳建筑大学，
类型：发明
国别省市：辽宁,21