一种基于块匹配的运动估计方法、运动估计系统及其应用
【技术领域】
[0001] 本发明涉及图像处理技术领域，尤其涉及一种基于块匹配算法的图像运动估计方 法、运动估计系统及其在细胞动力学分析上的应用。
【背景技术】
[0002] 基于计算机视觉的运动估计是一门有着广泛的应用前景的研究项目，例如在军事 领域的目标跟踪，工业过程的动态监控，商业领域的视频数据压缩、视频数据分析，医学领 域的心脏运动研究以及虚拟现实方面等方面。
[0003] 在随时间变化的视频序列中，帧与帧之间存在着很大的空间冗余，运动估计的目 标在于实现如下目标：如果景物和摄像设备都是静止的，则景物在当前帧中的位置与在下 一帧中的位置应当是相同的。如果在静止景物中还有运动的物体，则对当前帧中运动物体 上某一像素点，在未来时刻的最佳运动位置估计，应为该像素点在下一帧中的位置。
[0004]因此，运动估计可有效地去除冗余，保留帧间的有效信息，这对于图像序列（包括 视频）数据压缩和传输都非常重要。
[0005] 运动估计算法多种多样，大体上可以把它们分成四类：块匹配法、递归估计法、贝 叶斯估计法和光流法。其中块匹配算法最为简单有效，被广泛采用。块匹配算法中较为典 型的算法包括有：全搜索算法(在搜索窗（M+2w)x(N+2w)内计算所有的像素来寻找具有最 小误差的最佳匹配块。对于当前帧一个待匹配块的运动向量的搜索要计算（2w+l)x(2w+l) 次误差值）；三步搜索算法(三步搜索流程为：（1)以w/2为步长，测试以原点为中心的八 点；（2)以最小匹配误差点为中心，步长折半，测试新的八点；（3)重复第2步得到最后的运 动向量。TSS算法对于每一块的测试点为固定的（9 8 8) 25个。当位移大小w〈7时，三步 搜索相对于全搜索算法的加速因子为9)以及二维对数搜索（以跟踪最小均方差所在的方向 为主要思想。初始化计算五点，一点为原点，其他四点为（土w/2, 土w/2);再以相同的步长， 以上一步搜索到的最小点为中心测试点；然后，步长折半重复以上步骤，直到步长大小变为 1停止）等等。
[0006] 但在已有的块匹配算法中，全搜索算法虽然准确度高，但是速度太慢；而其他的快 速搜索方法通过限制搜索位置的数目来减小计算复杂度，但不利于估计小的运动且搜索容 易陷入局部最优；而且使用了矩形的匹配模板，只能进行平移操作，在目标做旋转运动时跟 踪效率过低。
[0007] 递归估计法虽然在运动目标较小时，递归估计算法收敛较快，还有利用过去信息 进行估计的能力。但是，当图像序列前后帧变化比较大的时候，这种方法很难得到正确的 结果，而且点位移的帧差绝对值常有多个最小值，造成递归估计算法的结果有时是局部最 小值而不是全局最小值。
[0008]贝叶斯估计法的计算量太大。而光流场方法的迭代速度虽然比较快，但得到的光 流场只是速度场的一种近似。在图像的亮度突变和运动不连续处，由于该算法所依赖的假 设不成立，因而所得的结果误差也较大。
【发明内容】
[0009] 鉴于上述现有技术的不足之处，本发明的目的在于提供一种基于块匹配的运动估 计方法、运动估计系统及其应用，旨在解决现有块匹配算法对旋转运动追踪效率过低，无法 兼顾运算速度及准确度的问题。
[0010] 为了达到上述目的，本发明采取了以下技术方案： 一种基于块匹配的运动估计方法，其中，所述方法包括如下步骤：获取微小物体的旋转 运动图像序列；通过基于块匹配的自定义块算法对目标微小物体进行平动补偿；对平动补 偿后所述目标微小物体进行像素块相关分析，以计算所述目标微小物体的旋转圈数；所述 基于块匹配的自定义算法使用圆形旋转模板。
[0011] 所述的运动估计方法，其中，所述方法还包括：对所述微小物体的旋转运动图像序 列通过使用高斯低通滤波器进行降噪以及通过直方图均衡化增强目标微小物体与背景的 对比度。
[0012] 所述的运动估计方法，其中，所述通过基于块匹配的自定义块算法对目标微小物 体进行平动补偿的步骤具体包括：生成圆形旋转模板；利用所述圆形旋转模板进行块匹 配，估计所述目标微小物体的运动轨迹。
[0013] 所述的运动估计方法，其中，所述生成圆形旋转模板的步骤具体包括：将目标图像 转换为灰度等级图像；通过自适应阈值算法将所述灰度等级图像转换为二值图像；计算所 述圆形旋转模板的圆心及半径，形成所述圆形旋转模板。
[0014] 所述的运动估计方法，其中，所述块匹配步骤具体包括：在搜索窗口依据预设的匹 配标准计算获得最佳匹配块作为参考块；计算所述参考块到宏块当前位置的运动向量。
[0015] 所述的运动估计方法，其中，所述匹配标准具体由如下算式表示：
其中，機.，s为搜索窗口，Μ为圆形 旋转模板。
[0016] 所述的运动估计方法，其中，与所述运动向量相对应的搜索窗口中绝对差和最小 的运动向量具体由如下算式计算具体由如下算式计算：
[0017] 所述的运动估计方法，其中，所述像素相关分析具体包括：通过算式1计 算模板与顺序图像块之间的相关系数；寻找所述相关系数的局部最大值以追 踪峰点；依据峰点的索引计算获得所述微小物体的旋转圈数；所述算式1为：
，其中，_、为模板，丨为顺序图像块，:_,和 /〇分别是和:_ι的平均数。
[0018] -种细胞运动动力学特性分析方法，其中，应用如上所述的运动估计方法对处于 介电泳力场的细胞进行旋转速度分析。
[0019] -种基于块匹配的运动估计系统，其中，所述系统包括：图像获取模块，用于获取 微小物体的旋转运动图像序列；平动补偿模块，用于通过基于块匹配的自定义块算法对目 标微小物体进行平动补偿；以及相关系数计算模块，用于对平动补偿后所述目标微小物体 进行像素块相关分析，以计算所述目标微小物体的旋转圈数；所述所述基于块匹配的自定 义算法使用圆形旋转模板。
[0020] 有益效果：本发明提供的一种基于块匹配的运动估计方法、运动估计系统及其应 用，针对目标的旋转运动跟踪，使用圆形旋转模板，从而有效的改进了块匹配算法对于目标 物体的旋转运动的追踪效率，在保持准确度的同时大幅度提高了块匹配算法的效率。并据 此提出了基于上述匹配算法应用于细胞旋转速度估计，对于细胞的动力学特性分析具有非 常深远的意义，能够为细胞动力学特性分析提供极大的便利。
【附图说明】
[0021] 图1为本发明具体实施例的光诱导介电泳平台的系统示意图。
[0022] 图2为本发明具体实施例的基于块匹配的运动估计方法的方法流程图。
[0023]图3为本发明具体实施例的运动估计方法中生成圆形旋转模板步骤的示意图。
[0024] 图4为本发明具体实施例的运动估计方法中块匹配步骤的示意图。
[0025]图5为本发明具体实施例的运动估计方法中生成圆形旋转模板步骤的方法流程 图。
【具体实施方式】
[0026] 本发明提供一种基于块匹配的运动估计方法、运动估计系统及其应用。为使本发 明的目的、技术方案及效果更加清楚、明确，以下参照附图并举实施例对本发明进一步详细 说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
[0027] 在本发明具体实施例中，使用处于介电泳力场中进行自转运动的染色细胞图像序 列为例进行陈述。应当理解的是，本发明所述基于块匹配的运动估计方法及其估计系统尤 其适用于其他图像特点相类似的微小物体进行旋转运动的图像序列分析，也可以应用于包 含所述染色细胞图像序列图像特点的其他任何合适的图像序列或者视频流等分析处理中。
[0028] 如图1所不，为获取所述细胞图像序列的系统的具体实施例。
[0029] 所述系统包括：0DEP芯片100、光学显微镜200、高分辨率投影仪300、可编程的 高精度0DEP芯片驱动器（图中未示出）以及计算机500。（亦即一个完整的光诱导介电泳平 台） 其中，所述计算机500可以采用现有技术中任何合适的，具有能够满足需要的运算能 力的电子运算设备或者平台，例如个人电脑、手提电脑、云端主机等等。
[0030] 0DEP是指光诱导介电泳（optically-induceddielectrophoresis)。所述0DEP芯 片具体可以由三层结构组成：其中，基底110是涂有一层1微米厚的氢化非晶硅（a_Si:H) 涂层11的ΙΤ0玻璃，最上层是普通ΙΤ0玻璃120,基底与最上层之间利用PDMS或是双面胶 封装出一个100微米高的微流体通道130。
[0031] 所述计算机通过数据线10与与投影仪连接，投影仪的镜头则与显微镜的入射光 路通过转接装置相连(如图1中A1所示)。0DEP芯片置于显微镜载物台上(其中图中200具 体为显微镜的物镜)。
[0032] 所述系统的使用过程为： 首先、培养细胞，并将细胞和基质放入ODEP芯片中。所述细胞具体可以为任何类型的 细胞，具体依据实际研究的需求而确定，例如肿瘤细胞，免疫活性细胞（ICC)等等。
[0033] 然后、在计算机上生成图案，由投影仪投射到位于显微镜载物台上的0DEP芯片 上，形成如图1中所示的图案（即图1中矩形框B1)。通过所述芯片驱动器400为0DEP芯片 加上交流电，改变交流电信号的频率和大小，并配合对应的投射图案，使0DEP芯片中的细 胞13做旋转运动。
[0034] 最后，通过显微镜的尚速(XD记录细胞运动的图像序列。
[0035] 如图2所示，为使用本发明所述的基于块匹配的运动估计方法对上述显微镜的高 速CCD记录的细胞运动的图像序列进行分析，计算细胞旋转圈数的方法流程图。
[0036] 所述方法包括如下步骤： S1、获取微小物体的旋转运动图像序列。如上所述，此处获取需要进行分析的图像序列 为显微镜的高速CCD记录