发明名称 基于直线二分法的不规则区域自动匹配方法
摘要 本发明涉及一种数字图像中基于位置二分法的不规则区域自动匹配方法,包括:采集图像并输入计算机;利用已有区域检测方法检测区域;计算区域最大对称位置;计算区域内各点的特征向量;利用梯度幅值极值计算区域关键点;计算区域关键点对应的区域描述向量;通过计算各区域描述向量的均值向量与标准差向量获得区域描述子;利用区域描述子进行区域匹配。相对于已有的区域匹配方法,本发明提供的方法既不需要进行固定形状拟合,又不需要进行分块处理,对图像形变具有更好的稳定性。
申请公布号 CN102663442B 申请公布日期 2014.02.26
申请号 CN201210069697.9 申请日期 2012.03.05
申请人 河南理工大学 发明人 刘红敏;王志衡;侯占伟;夏玉玲;王俊峰
分类号 G06K9/64(2006.01)I;G06T7/00(2006.01)I 主分类号 G06K9/64(2006.01)I
代理机构 代理人
主权项 1.一种数字图像中基于直线二分法的不规则区域自动匹配方法,其特征在于,包括步骤:步骤S1:从不同角度拍摄同一场景两幅不同图像并输入计算机;步骤S2:利用已有区域检测技术进行区域检测;步骤S3:计算每个区域G的最大对称位置P<sub>C</sub>,具体方式为,以区域G内的任一位置P为中心,引出18条直线L<sub>i</sub>将整个圆周等分为36份;记G内直线L<sub>i</sub>两侧的像素数分别为N<sub>L</sub>(i),N<sub>R</sub>(i),定义<img file="FSB0000115314640000011.GIF" wi="346" he="129" />为P处的不对称性;将区域G内不对称性最小的位置P<sub>C</sub>确定为区域G的最大对称位置;步骤S4:计算区域G内各点的特征向量,具体方式为,利用高斯梯度模板计算区域G内各点的梯度向量,记G内点X(x,y)处的高斯梯度为[d<sub>x</sub>,d<sub>y</sub>],G内各点的平均梯度为[V<sub>x</sub>,V<sub>y</sub>],计算点X(x,y)处的特征向量s=[s<sub>1</sub>,s<sub>2</sub>],其中s<sub>1</sub>=d<sub>x</sub>·V<sub>x</sub>+d<sub>y</sub>·V<sub>y</sub>,s<sub>2</sub>=d<sub>x</sub>·V<sub>y</sub>-d<sub>y</sub>·V<sub>x</sub>;步骤S5:利用梯度幅值极值计算区域关键点,具体方式为,记区域G内任一点X(x,y)处的梯度幅值为E(x,y),在阈值T约束下,将在3×3邻域内梯度幅值为极大值的点作为区域G内的关键点,即满足如下条件:E(x,y)>T,E(x,y)>E(x+1,y+1),E(x,y)>E(x-1,y-1),E(x,y)>E(x-1,y),E(x,y)>E(x+1,y),E(x,y)>E(x,y-1),E(x,y)>E(x,y+1),E(x,y)>E(x-1,y+1),E(x,y)>E(x+1,y-1);所述阈值T的具体确定方法为:T=Mean(E)+k·Std(E),Mean(E)与Std(E)分别表示所述区域内各点梯度幅值的均值与标准差,比例系数k的取值范围为2~3;步骤S6:计算区域G内关键点对应的区域描述向量,具体方式为,记步骤S5获得区域G的关键点分别为P<sub>1</sub>,P<sub>2</sub>...P<sub>n</sub>,关键点P<sub>i</sub>与区域G的对称位置P<sub>C</sub>确定的直线将区域G分为两个子区域,对两个子区域内各点的梯度幅值求和,记梯度幅值求和较大、较小的子区域分别为G<sub>B</sub>、G<sub>S</sub>;在两个子区域内区分正负对步骤S4获得的各点特征向量分量进行累加,可获得关键点P<sub>i</sub>对应的8维区域描述向量V<sub>i</sub>=[v<sub>i1</sub>,v<sub>i2</sub>,...,v<sub>i8</sub>],其中<maths num="0001"><![CDATA[<math><mrow><msub><mi>v</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>s</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>></mo><mn>0</mn><mi>andX</mi><mo>&Element;</mo><msub><mi>G</mi><mi>B</mi></msub></mrow></munder><msub><mi>s</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>v</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>s</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn><mi>andX</mi><mo>&Element;</mo><msub><mi>G</mi><mi>B</mi></msub></mrow></munder><msub><mi>s</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>v</mi><mrow><mi>i</mi><mn>3</mn></mrow></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>s</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>></mo><mn>0</mn><mi>andX</mi><mo>&Element;</mo><msub><mi>G</mi><mi>B</mi></msub></mrow></munder><msub><mi>s</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>v</mi><mrow><mi>i</mi><mn>4</mn></mrow></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>s</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn><mi>andX</mi><mo>&Element;</mo><msub><mi>G</mi><mi>B</mi></msub></mrow></munder><msub><mi>s</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>v</mi><mrow><mi>i</mi><mn>5</mn></mrow></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>s</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>></mo><mn>0</mn><mi>andX</mi><mo>&Element;</mo><msub><mi>G</mi><mi>S</mi></msub></mrow></munder><msub><mi>s</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>v</mi><mrow><mi>i</mi><mn>6</mn></mrow></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>s</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn><mi>andX</mi><mo>&Element;</mo><msub><mi>G</mi><mi>S</mi></msub></mrow></munder><msub><mi>s</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>v</mi><mrow><mi>i</mi><mn>7</mn></mrow></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>s</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>></mo><mn>0</mn><mi>andX</mi><mo>&Element;</mo><msub><mi>G</mi><mi>S</mi></msub></mrow></munder><msub><mi>s</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>v</mi><mrow><mi>i</mi><mn>8</mn></mrow></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>s</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn><mi>andX</mi><mo>&Element;</mo><msub><mi>G</mi><mi>S</mi></msub></mrow></munder><msub><mi>s</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>步骤S7:计算各关键点区域描述向量的均值与标准差获得区域描述子,具体方式为,记G内关键点P<sub>1</sub>,P<sub>2</sub>...P<sub>n</sub>确定的区域描述向量分别为V<sub>1</sub>,V<sub>2</sub>,...,V<sub>n</sub>,将V<sub>1</sub>,V<sub>2</sub>,...,V<sub>n</sub>各分量的均值组成一个向量并进行归一化获得G的8维均值向量V<sub>M</sub>=[v<sub>M1</sub>,v<sub>M2</sub>,...,v<sub>M8</sub>]/||v<sub>M1</sub>,v<sub>M2</sub>,...,v<sub>M8</sub>||,其中<img file="FSB0000115314640000029.GIF" wi="506" he="118" />||·||表示向量取模运算;将V<sub>1</sub>,V<sub>2</sub>,...,V<sub>n</sub>各分量的标准差组成一个向量并进行归一化获得区域G的8维标准差向量V<sub>S</sub>=[v<sub>S1</sub>,v<sub>S2</sub>,...,v<sub>S8</sub>]/||v<sub>S1</sub>,v<sub>S2</sub>,...,v<sub>S8</sub>||,其中<img file="FSB00001153146400000210.GIF" wi="704" he="143" />将均值向量V<sub>M</sub>与标准差向量V<sub>S</sub>组成一个向量并进行归一化可获得区域G的16维区域描述子D=[V<sub>M</sub>,V<sub>S</sub>]/||[V<sub>M</sub>,V<sub>S</sub>]||;步骤S8:利用区域描述子进行区域匹配,具体方式为,记第1幅图像中待匹配区域G<sub>11</sub>,G<sub>12</sub>、...,G<sub>1m</sub>的区域描述子分别为D<sub>11</sub>,D<sub>12</sub>、...,D<sub>1m</sub>,第2幅图像中的待匹配区域G<sub>21</sub>,G<sub>22</sub>、...,G<sub>2n</sub>的描述子分别为D<sub>21</sub>,D<sub>22</sub>、...,D<sub>2n</sub>,对于D<sub>11</sub>,D<sub>12</sub>、...,D<sub>1m</sub>中的任一区域描述子D<sub>1i</sub>,找到D<sub>21</sub>,D<sub>22</sub>、...,D<sub>2n</sub>中与D<sub>1i</sub>欧氏距离最小的区域描述子D<sub>2j</sub>,如果D<sub>1i</sub>同时也是D<sub>11</sub>,D<sub>12</sub>、...,D<sub>1m</sub>中与D<sub>2j</sub>欧氏距离最小的区域描述子,则G<sub>1i</sub>与G<sub>2j</sub>为一对匹配区域。
地址 454000 河南省焦作市高新区世纪大道2001号