专利名称：：一种剔除图像稳像系统周期性特征块的方法
技术领域：
：本发明属于图像稳像
技术领域：
，特别涉及一种对图像稳像算法中提取周期性特征块并加以剔除的方法。
背景技术：
：在图像稳像系统的稳像算法中，特征块应该检査是否具有周期性，周期性的特征块将会影响求解的局部运动矢量的精度进而可能影响整体稳像效果。现有技术中没有对特征块是否具有周期性进行判断的方法，目前多采用在图像的固定位置选取几块图像作为特征块，或者是利用最大均值差(或梯度和)的方法求出几块图像作为特征块。现有技术中的一种求解图像特征块的方法如下将图像划分为互相不重叠的R个区域，每个区域的大小为rs，这样在每一个区域中找一个ab的特征模板是可行的。满足下面条件则认为是特征区域S=Max(F》(2)式中，f(x，y)是每一个像素的值，E是该特征模板的均值，S最大的特征块就认为是特征模板。该方法基本可以保证特征模板内至少在某一点的一个领域内有明显变化；但由于没有考虑到特征模板(特征块)的周期性，不能有效地排除图像中特征模板(特征块)周围和该块图像相似的目标图像，在计算局部运动矢量时将会影响到局部运动矢量估计的精度和可靠性，并最终影响到稳像系统的整体稳像效果。
发明内容本发明的目的是针对上述现有技术的缺陷，提出一种可以有效提取图像稳像系统中周期性特征块以加以剔除的方法，以提高图像稳像系统的精度和稳像效果。本发明的技术方案是一种剔除图像稳像系统周期性特征块的方法，包括以下步骤1)、在第k帧图像上选取供第k+l帧图像求解局部运动矢量使用的、大小为RXS的搜索区域图像块S(g)，在搜索区域图像块S(g)中将图像分为MXN个大小为rXs的待选图像块T(g)，r、s为特征块水平和竖直方向的大小，R、S为搜索区域图像块的水平和竖直方向的大小，单位都是像素；2)、统计每个待选图像块T(g)的梯度和Grad(g)，根据梯度和Grad(g)排序、并同时根据一定阈值要求筛选出预选特征图像块T(h)，0〈h〈E，E=2XP，P为最终进行局部运动矢量提取时所用的图像特征块个数；各预选特征图像块T(h)对应的梯度和为Grad(h)，作为特征图像块的特征；所述筛选时的阈值要求为删除Grad(h)〈threshold的图像块；3)、对第h个预选特征图像块T(h)，利用最多近邻点距离一MCD的确定相似性的相关算法计算出预选特征图像块在搜索区域图像内的相关度量值D(A，B)，其计算公式为1S-5Jt-f广"一^其他式中，Ai，j表示预选特征图像块T(h)，Bi，j表示搜索区域图像中和预选特征图像块T(h)进行MCD相关算法运算的图像块；4)、统计预选特征图像块在搜索区域图像内的相关度量值D(A，B)大于设定阈值的位置出现次数PeriodCNT(h)，其统计公式为fPeriodCnt(h)+l当D(A，B)>Gthreshold!_o其他其中阈值G—threshold为用于判断是否与预先特征图像块的相关度量值D(A，B)相似的参考相似度量值；5)、根据PeriodCNT(h)的大小判断所选的预选特征图像块T(h)是不是周期性特征块，其判断公式为…f1当Peri。dCnt(h)XrXs)/.10Lo其他6)、输出判断结果并剔除周期性特征块，同时将保留的特征块提供给第k+l帧求解图像求解局部运动矢量；7)、继续处理第k+l帧图像。所述的步骤4中的阈值G—threshold为设定值，其取值范围为O.7-0.9。本发明充分利用图像稳像系统中包含候选特征块在内的整体图像的信息，有效的剔除了周期性的特征块，本方法不管图像块的周期性如何，既选取了最大特征的图像块，又最大限度地排除了周围和目标相似的目标，既适用于规则的周期性图像块的剔除，又适用于不规则的周期性图像块的剔除；本发明是稳像系统图像稳像算法中利用包含候选特征块和邻域图像块的相关性来剔除周期性的特征块的方法，它基于较大的比较区域进行相似性检验，提高了算法的抗干扰能力，并通过设定阈值来解决运动适量不一致的问题，从而保证了求解局部运动矢量的精度和可靠性；本方法中采用了最多近邻点距离(MCD)的确定相似性的原理，尽可能简化算法，解决了是规则的周期性特征块和不规则的周期性特征块的剔除问题，从而能大大提高了下一帧块匹配求取的局部运动矢量的精度和稳像效果。图1为本发明进行周期性特征块剔除的MCD相似性图像块求解示意图；图2为本发明剔除周期性特征块的程序流程示意图。具体实施例方式下面结合附图对本发明做进一步详细说明。如图l、图2所示，本发明的剔除周期性特征块方法主要根据第k帧图像，选出N个特征块，供第k+l帧图像求解局部运动矢量使用。在第k帧图像上选取供第&+1帧图像求解局部运动矢量使用的、大小为RXS的搜索区域图像块S(g)，在搜索区域图像块S(g)中将图像分为MXN个大小为rXs的待选图像块T(g)，r、s为特征块水平和竖直方向的大小，R、S为搜索区域图像块的水平和竖直方向的大小，单位都是像素；统计每个待选图像块T(g)的梯度和Grad(g)，根据梯度和Grad(g)排序、并同时根据一定阈值要求筛选出预选特征图像块T(h)，0〈h〈E，(E=2XP，P为最终进行局部运动矢量提取时所用的图像特征块个数)；各预选特征图像块T(h)对应的梯度和为Grad(h)，作为特征图像块的特征；所述筛选时的阈值要求为删除Grad(h)〈threshold的图像块；3)对第h个预选特征图像块T(h)，利用最多近邻点距离的确定相似性抓MCD的相关算法计算出预选特征图像块在搜索区域图像内的相关度量值D(A，B)，该算法如图l所示，是在预选特征图像块T(h)的最多近邻点距离选取和预选特征图像块T(h)相似的图像块以进行计算。MCD算法的计算公式如下其中，Ai，j表示特征块图像和搜索区域图像中，Bi，j特征块图像进行MCD相关算法进行运算的图像。4)统计预选特征图像块在搜索区域图像内的相关度量值D(A，B)大于设定阈值的位置出现次数PeriodCNT(h)，其统计公式为<formula>formulaseeoriginaldocumentpage7</formula>〔5)其中阈值G—threshold为用于判断是否与预先特征图像块的相关度量值D(A，B)相似的参考相似度量值；一般取值为0.7—0.9。5)根据PeriodCNT(h)的大小判断所选的是不是周期性特征块。具体方法是当Fe:riodCnt(h)>(Ws)/100其他(5)利用该方法剔除周期性图像块在图像稳像算法中，，实现步骤见图2为1、在第k帧将图像分为MXN个大小为rXs的图像块；2、统计每个图像块的梯度和，并进行排序；3、利用公式3和4针对第h个特征块，利用MCD相关算法在搜索区域图像块内进行相似性计算；4、利用公式5统计其相关度量值大于设定阈值的位置出现次数PeriodCNT(h)，作为当前特征块和邻域图像块相似性度量；5、利用公式6根据PeriodCNT(h)的大小判断所选的是不是周期性特征块。6、输出判断结果并剔除周期性特征块，同时将保留的特征块提供给第k+l帧求解图像求解局部运动矢量；7、继续处理第k+l帧图像。具体实施数据当前图像中存在周期性图像块。以下是对其中连续10帧图像计算后，经过第一步筛选选取出的特征最大的3个预选特征图像块的梯度和Grad(h)，当前预选特征图像块和邻域图像相关度量大于阈值的次数PeriodCNT(h)和最终确定的当前图像块是否具有周期性的结果blsPeriod(h):<table>tableseeoriginaldocumentpage7</column></row><table><table>tableseeoriginaldocumentpage8</column></row><table>从以上数据可以看出本发明的方法有效的剔除了周期性的特征块，大大提高了下一帧块匹配求取的局部运动矢量的精度，稳像效果良好。本发明是稳像系统图像稳像算法中利用包含候选特征块和邻域图像块的相关性来剔除周期性特征块的方法，它基于较大的比较区域进行相似性检验，提高了算法的抗干扰能力，并通过设定阈值来解决运动适量不一致的问题，从而保证了求解局部运动矢量的精度和可靠性。该方法中采用了最多近邻点距离(MCD)的确定相似性的原理，尽可能简化算法，解决了是规则的周期性特征块和不规则的周期性特征块的剔除问题，从而能大大提高了下一帧块匹配求取的局部运动矢量的精度和稳像效果。本发明的技术关键是利用最多近邻点距离MCD的确定相似性的原理。权利要求1.一种剔除图像稳像系统周期性特征块的方法，其特征在于，包括以下步骤1)、在第k帧图像上选取供第k+1帧图像求解局部运动矢量使用的、大小为R×S的搜索区域图像块S(g)，在搜索区域图像块S(g)中将图像分为M×N个大小为r×s的待选图像块T(g)，r、s为特征块水平和竖直方向的大小，R、S为搜索区域图像块的水平和竖直方向的大小，单位都是像素；2)、统计每个待选图像块T(g)的梯度和Grad(g)，根据梯度和Grad(g)排序、并同时根据一定阈值要求筛选出预选特征图像块T(h)，0＜h＜E，E＝2×P，P为最终进行局部运动矢量提取时所用的图像特征块个数；各预选特征图像块T(h)对应的梯度和为Grad(h)，作为特征图像块的特征；所述筛选时的阈值要求为删除Grad(h)＜threshold的图像块；3)、对第h个预选特征图像块T(h)，利用最多近邻点距离朚CD的确定相似性的相关算法计算出预选特征图像块在搜索区域图像内的相关度量值D(A，B)，其计算公式为<mathsid="math0001"num="0001"><math><![CDATA[<mrow><mi>D</mi><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>B</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>r</mi><mo>×</mo><mi>s</mi></mrow></mfrac><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>S</mi><mo>-</mo><mi>s</mi></mrow></munderover><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>R</mi><mo>-</mo><mi>r</mi></mrow></munderover><mi>R</mi><mrow><mo>(</mo><msub><mi>A</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo><msub><mi>B</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow></mrow>]]></math></maths>式中，Ai，j表示预选特征图像块T(h)，Bi，j表示搜索区域图像中和预选特征图像块T(h)进行MCD相关算法运算的图像块；4)、统计预选特征图像块在搜索区域图像内的相关度量值D(A，B)大于设定阈值的位置出现次数PeriodCNT(h)，其统计公式为其中阈值G_threshold为用于判断是否与预先特征图像块的相关度量值D(A，B)相似的参考相似度量值；5)、根据PeriodCNT(h)的大小判断所选的预选特征图像块T(h)是不是周期性特征块，其判断公式为输出判断结果并剔除周期性特征块，同时将保留的特征块提供给第k+1帧求解图像求解局部运动矢量；7)、继续处理第k+1帧图像。2根据权利要求l所述的剔除图像稳像系统周期性特征块的方法，其特征在于，所述的步骤4中的阈值G—threshold为设定值，其取值范围为O.7-0.9。其判断公式为:1当PeriodCnt(li)XrXs)/100其他全文摘要本发明公开了一种剔除图像稳像系统周期性特征块的方法，其将第k帧图像分成M×N个图像块，统计每块的梯度和，并将最终的梯度和排序选出特征明显的图像块，然后针对第h个特征块，统计其邻域内相关度量值大于阈值的次数PeriodCNT(h)，最后根据PeriodCNT(h)的大小判断所选的是不是周期性特征块，并加以剔除。本发明充分利用了稳像系统中包含候选特征块在内的整体图像的信息，不管图像块的周期性如何，既选取了最大特征的图像块，又最大限度地排除了周围和目标相似的目标，既适用于规则的周期性图像块的剔除，又适用于不规则的周期性图像块的剔除，从而能大大提高下一帧块匹配求取的局部运动矢量的精度和稳像效果。文档编号H04N7/26GK101600115SQ20091030382公开日2009年12月9日申请日期2009年6月29日优先权日2009年6月29日发明者吕高杰申请人:中国航空工业集团公司洛阳电光设备研究所