发明名称 基于特征点邻域伪Zernike矩的图像拼接方法
摘要 本发明涉及一种基于特征点邻域伪Zernike矩的图像拼接方法。该方法是将图像统计信息(伪Zernike矩)与图像的特征点相联系,首先利用Harris角检测器提取出参考图像与输入图像的兴趣点,并将兴趣点为中心的矩形邻域作为特征匹配的局部特征区域,然后对这些矩形特征区域计算伪Zernike矩作为该特征区域的描述符,并通过比较各个特征区域描述符特征向量的欧氏距离来实现特征点的匹配。这样的匹配中必然会存在少部分不正确的匹配点,通过RANSAC(RANdom Sample Consensus)算法剔除不正确的匹配点对,并计算出正确的匹配关系以实现图像的拼接。本发明能有效地实现具有平移,旋转,小尺度缩放等几何变换关系的图像配准和拼接,可用于图像处理与通信、多媒体技术等领域的图像合成。
申请公布号 CN101556692A 申请公布日期 2009.10.14
申请号 CN200810017909.2 申请日期 2008.04.09
申请人 西安盛泽电子有限公司 发明人 郭宝龙;杨占龙;闫允一
分类号 G06T5/50(2006.01)I;G06T7/00(2006.01)I;G06T3/00(2006.01)I 主分类号 G06T5/50(2006.01)I
代理机构 代理人
主权项 1.一种基于特征点邻域伪Zernike矩的图像拼接方法,包括如下步骤:(1)利用Harris角点检测方法对参考图像和待拼接图像提取特征点,分别记作p={p<sub>1</sub>,p<sub>2</sub>,p<sub>3</sub>,…,p<sub>m</sub>}和p′={p′<sub>1</sub>,p′<sub>2</sub>,p′<sub>3</sub>,…,p′<sub>n</sub>},以特征点为中心取矩形的区域作为该特征点的特征区域;(2)在所有的特征区域中,分别计算该区域的伪Zernike矩,并将生成的特征向量作为该特征点的描述符;(3)通过比较兴趣点之间特征向量S的欧氏距离可以得到两幅图像兴趣点之间的对应关系。即对于参考图像中的每一个兴趣点S<sub>p</sub>,分别计算该点的特征向量与待拼接图像中每一个兴趣点特征向量的欧氏距离,距离越小则说明两点之间的特征向量越相似,取距离最小的点作为参考图像中点S<sub>p</sub>在待拼接图像中的匹配点,记作S′<sub>p</sub>;(4)通过上述第三步,可以计算出参考图像中每一个兴趣点在待拼接图像中的匹配点,当然这样的匹配对中会包含部分伪匹配对,应用RANSAC算法估计出这些匹配点对的几何变换模型T,同时剔除伪匹配对;<maths num="0001"><![CDATA[<math><mrow><mi>T</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>s</mi><mi>cos</mi><mi>&theta;</mi></mtd><mtd><mo>-</mo><mi>s</mi><mi>sin</mi><mi>&theta;</mi></mtd><mtd><msub><mi>t</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><mi>s</mi><mi>sin</mi><mi>&theta;</mi></mtd><mtd><mi>s</mi><mi>cos</mi><mi>&theta;</mi></mtd><mtd><msub><mi>t</mi><mi>y</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中s是尺度变换参数,θ是旋转的角度,t<sub>x</sub>,t<sub>y</sub>分别是x方向和y方向上的平移量。(5)将待拼接图像按照几何变换模型T变换到参考图像的坐标系上;x=Tx′其中x′是待拼接图像的原始坐标向量,T是几何变换矩阵,x是变换到参考图像坐标系上后的新的坐标向量(6)使用加权平均法将两幅图像的重叠区域进行图像融合,实现拼接图像中重叠区域的平滑过渡,<maths num="0002"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>f</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mtd><mtd><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&Element;</mo><msub><mi>f</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>d</mi><mn>1</mn></msub><msub><mi>f</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>d</mi><mn>2</mn></msub><msub><mi>f</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mtd><mtd><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&Element;</mo><mrow><mo>(</mo><msub><mi>f</mi><mn>1</mn></msub><mo>&cap;</mo><msub><mi>f</mi><mn>2</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>f</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mtd><mtd><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&Element;</mo><msub><mi>f</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中d<sub>1</sub>,d<sub>2</sub>表示权重值,与重叠区域的宽度有关,并且d<sub>1</sub>+d<sub>2</sub>=1,0≤d<sub>1</sub>,d<sub>2</sub>≤1。在重叠区域中d<sub>1</sub>由1渐变至0,d<sub>2</sub>由0渐变至1,由此实现在重叠区域中由f<sub>1</sub>平滑过渡到f<sub>2</sub>,f<sub>1</sub>,f<sub>2</sub>是两幅待拼接的图像,f是融合后的图像。d<sub>1</sub>和d<sub>2</sub>的计算如下:假设当前像素的横坐标为x<sub>i</sub>,重叠区域左右边界的横坐标分别为x<sub>l</sub>和x<sub>r</sub>,那么<maths num="0003"><![CDATA[<math><mrow><msub><mi>d</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>x</mi><mi>r</mi></msub><mo>-</mo><msub><mi>x</mi><mi>i</mi></msub></mrow><mrow><msub><mi>x</mi><mi>r</mi></msub><mo>-</mo><msub><mi>x</mi><mi>l</mi></msub></mrow></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mi>d</mi><mn>2</mn></msub><mo>=</mo><mn>1</mn><mo>-</mo><msub><mi>d</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>x</mi><mi>l</mi></msub></mrow><mrow><msub><mi>x</mi><mi>r</mi></msub><mo>-</mo><msub><mi>x</mi><mi>l</mi></msub></mrow></mfrac><mo>.</mo></mrow></math>]]></maths>
地址 710077陕西省西安市高新技术开发区