可移动成像设备的图像运动检测方法及图像降噪方法
【专利摘要】一种可移动成像设备的图像运动检测方法，包括：根据特定样率对前一帧和当前帧图像进行采样，分别获得前一帧尺度图和当前帧尺度图；在前一帧尺度图中确定搜索区，在当前帧尺度图中确定搜索块，所述搜索区覆盖所述搜索块；在所述搜索区搜索与所述搜索块匹配的区域，确定所述匹配的区域为匹配块；根据所述搜索块与所述匹配块的位置关系确定图像的全局运动。
【专利说明】
可移动成像设备的图像运动检测方法及图像降噪方法
技术领域
[0001]本发明涉及图像处理领域，特别涉及图像运动检测方法及图像降噪方法。
【背景技术】
[0002]目前，可移动医学成像设备采集的图像中图像的降噪一般采用基于时阈的递归降噪，对于静止的图像序列，递归系数大，降噪效果好，但对于运动图像，递归降噪的递归系数大就会弓I起较明显的运动伪影。
[0003]以下是现有递归降噪的说明。在采用时域递归滤波方法进行降噪时，先把当前帧图像的像素和前一帧图像的经过降噪处理的像素按一定的比例叠加，其过程如以下公式:
[0004]Pn' (X,y) =Pn(x，y)*k+Pn_l ’(x，y)*(l_k)。
[0005]其中，Pn-Γ是指前一帧降噪后的图像，对于降噪系数k，k值越大，当前帧占的比重就大，递归就浅，降噪不明显;k值越小，则当前帧所占比重越小，递归深，降噪明显。但是由于使用前一帧与当前帧的加权叠加，对于运动的目标就会在降噪后的图中看到前一帧的“残影”，k值取得越小，“残影”越明显。
[0006]为了减少运动伪影，需要根据图像的运动改变递归系数。两帧图像之间相对运动越大，k值就应该取得大，伪影减少的同时降噪效果也相应减小，因此需要在减小运动伪影和消除噪声之间取的比较平衡的效果，k值就要取得合适。K值的选取需要根据图像之间实际的运动情况来进行选取，因此需要运动动情况的检测。
[0007]在现有技术中关于运动情况的检测是基于运动补偿的检测方法:主要应用于图像整体发生平移运动，如设备的移动导致的图像的平移。该方法对当前帧取一块区域，作为搜索块，在前一帧取一搜索区域，用搜索块在搜索区域中找出相似度最大的块区域，作为匹配块，利用这一匹配块的和搜索块的位置估算位移作为整个图像的位移情况。而这种方法计算量大，影响图像的处理效率。
【发明内容】
[0008]本发明为了解决现有技术中提到的问题，提供一种减少计算量的同时又能精确测量图像全局运动的方法及根据所述检测方法对图像进行降噪的图像降噪方法。
[0009]为了解决上述问题，本发明提供一种可移动成像设备的图像运动检测方法，包括:
[0010]根据特定采样率对前一帧和当前帧图像进行采样，分别获得前一帧尺度图和当前帧尺度图；
[0011]在前一帧尺度图中确定搜索区，在当前帧尺度图中确定搜索块，所述搜索区覆盖所述搜索块；
[0012]在所述搜索区搜索与所述搜索块匹配的区域，确定所述匹配的区域为匹配块；
[0013]根据所述搜索块与所述匹配块的位置关系确定图像的全局运动。
[0014]可选的，所述的图像运动检测方法，还包括:
[0015]在采样之前对所述前一帧和当前帧图像进行去无效区的步骤。
[0016]可选的，所述的图像运动检测方法中，
[0017]所述搜索区是在所述前一帧尺度图中对应所述搜索块的区域上向外扩展不小于设备移动距离对应的像素位后形成的区域。
[0018]可选的，所述的图像运动检测方法中，
[0019]所述搜索区是在所述前一帧尺度图中对应所述搜索块的区域上在设备移动方向上扩展不小于设备移动距离对应的像素位后形成的区域。
[0020]本发明还提供一种图像运动检测方法，还包括:
[0021]根据第一采样率对前一帧和当前帧图像进行采样，分别获得前一帧第一尺度图和当前帧第一尺度图；
[0022]在前一帧第一尺度图中确定第一搜索区，在当前帧第一尺度图中确定第一搜索块，所述第一搜索区覆盖所述第一搜索块；
[0023]在所述第一搜索区搜索与所述第一搜索块匹配的区域，确定所述匹配的区域为第一匹配块；
[0024]根据第二采样率，对前一帧和当前帧图像进行采样，分别获得前一帧第二尺度图和当前帧第二尺度图，所述第二采样率要大于所述第一采样率；
[0025]在前一帧第二尺度图中确定对应第一匹配块的区域为第二匹配块；
[0026]将所述第二匹配块所在的区域向外扩展获得第二搜索区；
[0027]在当前帧第二尺度图中确定对应第一搜索块的区域为第二搜索块；
[0028]在所述第二搜索区搜索与所述第二搜索块匹配的区域，确定该匹配的区域为第三匹配块；
[0029]根据所述第二搜索块与所述第三匹配块的位置关系确定图像的全局运动。
[0030]可选的，所述图像运动检测方法中，
[0031 ]所述第一采样率和所述第二采样率均为偶数分之一采样率。
[0032]可选的，所述图像运动检测方法中，
[0033]所述第一采样率为l/2n;所述第二采样率为1/n，其中η为偶数。
[0034]可选的，所述图像运动检测方法中，
[0035]所述第一采样率为1/4，所述第二采样率为1/2。
[0036]可选的，所述图像运动检测方法中，
[0037]所述前一帧第二尺度图中对应第一匹配块的区域向外扩展是指在前一帧第二尺度图中对应第一匹配块顶点的像素中取最外面的像素作为第二匹配块的顶点。
[0038]本发明还提供一种可移动成像设备的图像降噪方法，包括:
[0039]根据所述的方法计算前一帧和当前帧图像的全局运动；
[0040]根据所述全局运动计算前一帧和当前帧图像的重叠区域；
[0041 ]将前一帧和当前帧图像中的重叠区域分别分割为对应的至少两个子块；
[0042 ]对前一帧和当前帧图像中对应子块的灰度值进行比较;及
[0043]根据比较结果调整降噪系数。
[0044]与现有技术相比，本发明综合考虑了图像帧间的运动情况，针对全局运动和局部运动的实时图像提出一种运动检测方法，结合了整体运动位移估计和局部运动检测，更全面的判断图像帧间运动情况，并根据不同的运动情况设定不同的降噪系数，在对图像进行降噪的同时也减少图像的伪影。
[0045]进一步地，通过对图像采样后检测全局运动的方法，在有效的降低计算量的同时不会影响运动检测的结果。
【附图说明】
[0046]图1是本发明的图像全局运动检测方法的流程示意图；
[0047]图2是本发明的根据第一尺度图进行匹配的过程示意图；
[0048]图3是本发明的搜索区的生成示意图；
[0049]图4是本发明的根据第二尺度图进行匹配的过程示意图；
[0050]图5是本发明的根据两次匹配过程的获取的匹配块示意图；
[0051 ]图6是本发明的局部运动检测方法的流程示意图；
[0052]图7是根据本发明的局部运动检测过程示意图。
【具体实施方式】
[0053]为使本发明的上述目的、特征和优点能够更为明显易懂，下面结合附图对本发明的【具体实施方式】做详细的说明。在以下描述中阐述了具体细节以便于充分理解本发明。但是本发明能够以多种不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似推广。因此本发明不受下面公开的【具体实施方式】的限制。
[0054]图1是本发明的图像全局运动检测方法的流程示意图。以下结合附图1对图像全局运动的检测方法进行说明。
[0055]如图1所示，图像全局运动检测包括以下步骤。
[0056]S10，在前一帧和当前帧图像中去掉无效区域，生成有效区域图像。在医学诊断图像中，如X线成像设备生成的图像里，有效区域指与诊断相关的区域，也就是X线穿过人体组织后形成的图像区域，而无效区域是X线未经过人体组织后形成的图像区域，因为不同物质对X线的吸收不同，因此经过人体的图像区域和未经过人体的图像区域在图像中的灰度值会有明显的差异，在这一步骤中首先设置灰度阈值，将图像中的像素灰度值与所述阈值进行比较，以此判断该像素点是否为有效区域的像素点，去掉无效背景区域，只保留感兴趣区域。通过上述步骤可以减少需要计算的像素点数，能有效的提高计算效率。
[0057]Sll，分别生成前一帧和当前帧的第一尺度图。对上一步骤中去掉无效区域后的图像进行采样，所述采样可以是行列各取1/4的下采样方法，是对行或列每隔三行取一行(相当于对原始图像进行1/16采样)。在此也可以其他采样率，如1/3或1/5等，此处不做限定。
[0058]S12，在前一帧第一尺度图中确定第一搜索区。图2是本发明的根据第一尺度图进行匹配的过程示意图，如图2所示，在当前帧第一尺度图120中确定第一搜索块121，在前一帧第一尺度图110中确定第一搜索区111，所述第一搜索区111是在前一帧第一尺度图110的图像中心范围获取，所述第一搜索块121是在当前帧的第一尺度图120的图像中心范围获取。所述第一搜索区111要大于所述第一搜索块121。
[0059]在这一步骤中，所述第一搜索块121不易太大或太小，搜索块的大小的选择一般取决于设备的移动速度，一般对于前后两帧不会出现平移很大的情况，如果出现很大的情况，有可能都是运动区域，在后续计算的过程中会不匹配，算法会终止。若第一搜索块太大对减少计算量没有太大的帮助，若太小则有可能会出现计算误差，所述搜索块的大小可以在当前第一尺度图大小的1/48至1/16之间取，其中1/32为佳。
[0060]图3为搜索区的生成示意图。搜索区的大小根据设备移动情况确定，在本实施例中，若设备移动为L个像素位，则在前一帧第一尺度图中的中心范围对应第一搜索块的区域再向外扩展2L像素位的范围，扩展后的区域定义为第一搜索区。在其他实施例中也可以用其他方式定义搜索区，如，当判断设备的移动仅是在A-B方向上发生，则仅在A-B方向扩展2L像素位的范围，而在C-D方向上只增加IL像素位的范围。也即，在设备移动发生幅度较少或没有发生移动的方向上增加的像素位范围较少或不增加，移动发生幅度较多的方向上增加的范围也多，以全面覆盖搜索块对应的区域为准。上述设备移动可以是设备的主动移动或被动移动。
[0061]S13:在第一搜索区对第一搜索块进行相似度匹配，确定最佳匹配区域，并定义为第一匹配块。结合图2，在确定所述第一匹配块112以后记录第一匹配块112的顶点坐标，如:左下角xl，yl，右上角x2，y2，在此也可以记录右下角和左上角的坐标，在此不做限定。其中，相似度匹配的过程是采用比较块内的灰度差，差值越小，表示越相似。
[0062]S14:分别生成前一帧和当前帧的第二尺度图。所述第二尺度图可以是行或列各取1/2下采样，是对行或列每隔一行取一行。在此也可以是其他方式，但是第二尺度图中采样的像素点要多于第一尺度图。
[0063]S15:在前一帧第二尺度图中确定第一匹配块的位置，定义为第二匹配块。图4是本发明的根据第二尺度图进行匹配的过程示意图。如图4所示，在当前帧第二尺度图220中，确定对应第一搜索块121的区域为第二搜索块221;在前一帧第二尺度图210中，将对应第一匹配块112的区域确定为第二匹配块212，所述第二搜索块221和第二匹配块212所在区域可以通过第一尺度图和第二尺度图的采样比关系确定。
[0064]S16:对第二匹配块212进行扩展，生成第二搜索区211。图5是本发明的根据两次匹配过程获取的匹配块的示意图。如图4、图5所示，在上述步骤中确定的第二匹配块212的两个顶点坐标为左下角(xl ’，yl’)，右上角(x2 ’，y2 ’)，在第二匹配块212周围往外扩一个范围，在前一帧第二尺度图210上形成新的搜索区范围，也即第二搜索区211，新的搜索区211左下角坐标为(xl，_range/2;yl ’-range/2)，右上角坐标为(x2’+range/2，y2’+range/2)。
[0065]所述在第二匹配块212周围向外扩一个范围是指，按特定的像素位向外扩展，如扩10个像素位;或可以是对应像素点位的扩展，如在1/4采样的第一尺度图中的一个顶点像素点对应1/2采样的第二尺度图的四个像素点，在对应的4个像素点中取最靠外的像素点坐标。
[0066]S17:在第二搜索区211对第二搜索块221进行相似度匹配，确定最佳匹配区域，并定义为第三匹配块214。
[0067]S18:在前一帧和当前帧图像中分别确定第三匹配块214和第二搜索块221的位置，并计算整体图像位移。
[0068]在进行全局快速匹配时，获取尺度图时不能使原图像缩小的太多，缩小的多可能导致获取的尺度图中图像的信息损失过多，进而导致在搜索区匹配的位置误差大。在上述实施例中第一尺度图和第二尺度图分别选择1/4和1/2的采样率，使原图信息损失不大，搜索后准确性不降，且达到了减少计算量的目的。
[0069]进一步，在上述方法中，在尺度选择的时候，考虑了以偶数倍进行降采样，这样在缩放的时候图像大小容易控制，方便计算，若采用奇数倍进行降采样(如1/3或1/5等)在计算的过程中有可能会出现除不尽，在这种情况下换算坐标和长度的时候就需要近似换算，因此在采样的过程中最好以偶数分子I的采样较佳。
[0070]图6是本发明的图像局部运动检测方法的流程示意图。以下结合附图6对局部运动检测方法进行说明。
[0071]S20，确定前一帧图像和当前帧图像的重合区域。图7是根据本发明的图像局部运动检测过程示意图，如图7所示，当前帧图像经过全局位移后，某些区域可能和前一帧不重合，这部分区域标识为新增区域311，而前一帧图像中也有部分与当前帧图像不重合，这部分区域标识为移出区域321。
[0072]S21，将前一帧图像和当前帧图像中的重叠区域分别分割为对应的n*m块。其中所述n、m分别为自然数，n、m可以是相同的自然数，本发明中每块大小为32*32个像素点(也可以是其他大小)。
[0073]S22，比对重叠区域对应块的灰度值。对当前帧图像的第I块的灰度值均值与前一帧的第I’块的灰度值均值进行比较，灰度值差异计算对应了前后两帧对应块的匹配程度，作为每个块的运动强度标识。设定两个阈值，灰度差异小于小阈值的为静止或缓慢运动区，差异大于大阈值时为强运动区，小于大阈值、大于小阈值时为中等运动区。
[0074]在根据上述步骤的结果，适当的调整递归系数，如强运动区域采用较小的递归系数;运动缓慢或静止区域采用较大的递归系数;而处于中等运动区的采用相对中等的递归系数。
[0075]根据本发明的技术方法有如下效果，以搜索块的大小定义为128*128，搜索范围的大小定义为384*384为例，在不进行降采样，进行块匹配过程中整个搜索完成的块匹配次数需要257*257 = 66049;在进行降采样为1/4的情况下，搜索块则是128/4，即搜索块的大小为32*32，搜索范围是384/4，也即搜索区的大小为96*96，在1/4降采样图像上，采用由粗到精的过程进行快速搜索匹配，全部搜索完进行的匹配次数是65*65次;在1/4采样图搜索完成后，换算到1/2采样图上去计算，此时搜索块大小变成64*64，且由于1/4采样图上已经找到近似位置，1/2采样图上只需要在近似位置附近往外再扩展误差范围即可，这里假设扩展10个像素误差范围，则搜索范围变为74*74，全部搜索完只需要11*11次;在1/2采样图上搜索完后，换算到原图做最后搜索，搜索块变为128*128，同时由于在1/2采样图上找到近似匹配块，只需要在这个近似位置往外扩展一个误差范围即可，假设也扩展10个像素为误差范围，则搜索范围变为138*138，全部搜索完只需要11*11次。比较的结果为在降采样下全部搜索的次数是65*65+11*11+11*11 =4467。根据上述分析，使用降采样，由粗到细的快速匹配方式使得匹配次数大大减少，算法次数从66049减少到4467次，节省了算法时间；而且在降采样图上的搜索块小，匹配计算也快一些。
[0076]本发明提出了一种基于全局运动估算和局部运动标识相结合的运动检测方法。本算法把全局运动和局部运动结合到一起，综合考虑了运动的情况，全面判断了图像的运动情况，适用于多种运动情况。本算法在进行全局运动估算时用了基于多尺度的快速估算，从低尺度到高尺度的逐级搜索，加快了全局运动的搜索速度。
[0077]本发明虽然已以较佳实施例公开如上，但其并不是用来限定本发明，任何本领域技术人员在不脱离本发明的精神和范围内，都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改，因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰，均属于本发明技术方案的保护范围。
【主权项】
1.一种可移动成像设备的图像运动检测方法，其特征在于，包括: 根据特定采样率对前一帧和当前帧图像进行采样，分别获得前一帧尺度图和当前帧尺度图； 在前一帧尺度图中确定搜索区，在当前帧尺度图中确定搜索块，所述搜索区覆盖所述搜索块； 在所述搜索区搜索与所述搜索块匹配的区域，确定所述匹配的区域为匹配块； 根据所述搜索块与所述匹配块的位置关系确定图像的全局运动。2.根据权利要求1所述的图像运动检测方法，其特征在于，还包括: 在采样之前对所述前一帧和当前帧图像进行去无效区的步骤。3.根据权利要求1所述的图像运动检测方法，其特征在于， 所述搜索区是在所述前一帧尺度图中对应所述搜索块的区域上向外扩展不小于设备移动距离对应的像素位后形成的区域。4.根据权利要求1所述的图像运动检测方法，其特征在于， 所述搜索区是在所述前一帧尺度图中对应所述搜索块的区域上在设备移动方向上扩展不小于设备移动距离对应的像素位后形成的区域。5.一种可移动成像设备的图像运动检测方法，其特征在于， 根据第一采样率对前一帧和当前帧图像进行采样，分别获得前一帧第一尺度图和当前帧第一尺度图； 在前一帧第一尺度图中确定第一搜索区，在当前帧第一尺度图中确定第一搜索块，所述第一搜索区覆盖所述第一搜索块； 在所述第一搜索区搜索与所述第一搜索块匹配的区域，确定所述匹配的区域为第一匹配块； 根据第二采样率，对前一帧和当前帧图像进行采样，分别获得前一帧第二尺度图和当前帧第二尺度图，所述第二采样率要大于所述第一采样率； 在前一帧第二尺度图中确定对应第一匹配块的区域为第二匹配块； 将所述第二匹配块所在的区域向外扩展获得第二搜索区； 在当前帧第二尺度图中确定对应第一搜索块的区域为第二搜索块； 在所述第二搜索区搜索与所述第二搜索块匹配的区域，确定该匹配的区域为第三匹配块； 根据所述第二搜索块与所述第三匹配块的位置关系确定图像的全局运动。6.根据权利要求5所述的图像运动检测方法，其特征在于， 所述第一采样率和所述第二采样率均为偶数分之一。7.根据权利要求6所述的图像运动检测方法，其特征在于， 所述第一采样率为l/2n;所述第二采样率为1/n，其中η为偶数。8.根据权利要求6或7所述的图像运动检测方法，其特征在于， 所述第一采样率为1/4，所述第二采样率为1/2。9.根据权利要求5所述的图像运动检测方法，其特征在于， 所述前一帧第二尺度图中对应第一匹配块的区域向外扩展是指在前一帧第二尺度图中对应第一匹配块顶点的像素中取最外面的像素作为第二匹配块的顶点。10.一种可移动成像设备的图像降噪方法，其特征在于，包括: 根据权利要求1-9任一项所述的方法计算前一帧和当前帧图像的全局运动； 根据所述全局运动计算前一帧和当前帧图像的重叠区域； 将前一帧和当前帧图像中的重叠区域分别分割为对应的至少两个子块； 对前一帧和当前帧图像中对应子块的灰度值进行比较；及根据比较结果调整降噪系数。
【文档编号】H04N5/14GK105872310SQ201610247824
【公开日】2016年8月17日
【申请日】2016年4月20日
【发明人】江春花, 赵书睿, 周海华, 边钺岩
【申请人】上海联影医疗科技有限公司