发明名称 基于Skinner概率自动机的图像误匹配对去除方法
摘要 本发明公开了一种基于Skinner概率自动机的图像误匹配对去除方法,包括:进行图像匹配,得到匹配对;确定匹配结果的Skinner概率自动机数学模型;从匹配对中按照概率随机抽取匹配对,并且求出图像间的基础矩阵,然后将此基础矩阵带入所有匹配对,根据判决函数求出每个匹配对之间的极对距离和平均极对距离;根据权重调节函数、每个匹配对的极对距离和平均极对距离对每个匹配对进行权重的调整,并且计算每个匹配对的提取概率;判断当前基础矩阵是否具有最多正确匹配对数目,并结合三个判断条件控制算法迭代过程。本发明引入仿生学和认知心理学,在缺乏错匹配概率先验知识的情况下提出了三种迭代终止条件,提高了算法的自主性。
申请公布号 CN103700080B 申请公布日期 2016.08.03
申请号 CN201310656943.5 申请日期 2013.12.06
申请人 北京工业大学 发明人 阮晓钢;魏若岩;武璇;于乃功;陈志刚;肖尧;瓦达哈·谢
分类号 G06T5/50(2006.01)I 主分类号 G06T5/50(2006.01)I
代理机构 北京思海天达知识产权代理有限公司 11203 代理人 张慧
主权项 一种基于Skinner概率自动机的图像误匹配对去除方法,其特征在于包括以下步骤:步骤一,进行图像特征点或特征区域的检测与匹配,得到匹配对;步骤二,根据图像匹配结果确定匹配结果的Skinner概率自动机数学模型,其表达式为:Skinner‑Ransac={M,W,P,T,O,Inmost,Mo,N,Stc}M={m<sub>i</sub>|i=1,2,…,n},m<sub>i</sub>={(x<sub>i</sub>,y<sub>i</sub>),(x′<sub>i</sub>,y′<sub>i</sub>)}W={w<sub>i</sub>|i=1,2,…,n}<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>P</mi><mo>=</mo><mo>{</mo><msub><mi>p</mi><mi>i</mi></msub><mo>|</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>n</mi><mo>}</mo><mo>,</mo><msub><mi>p</mi><mi>i</mi></msub><mo>=</mo><msub><mi>w</mi><mi>i</mi></msub><mo>/</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>w</mi><mi>j</mi></msub></mrow>]]></math><img file="FDA0000951020530000011.GIF" wi="789" he="127" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>T</mi><mo>:</mo><msub><mi>d</mi><mi>i</mi></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>c</mi><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>d</mi><mi>i</mi></msub><mo>&GreaterEqual;</mo><mi>c</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>d</mi><mi>i</mi></msub><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>d</mi><mi>i</mi></msub><mo>&lt;</mo><mi>c</mi></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>n</mi></mrow>]]></math><img file="FDA0000951020530000012.GIF" wi="678" he="148" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>O</mi><mo>&RightArrow;</mo><msub><mi>w</mi><mi>i</mi></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>w</mi><mi>i</mi></msub><mo>+</mo><mi>m</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mi>r</mi><mi>o</mi><mi>u</mi><mi>n</mi><mi>d</mi><mo>(</mo><mrow><mover><mi>d</mi><mo>&OverBar;</mo></mover><mo>/</mo><msub><mi>d</mi><mi>i</mi></msub></mrow><mo>)</mo><mo>,</mo><mi>H</mi><mo>)</mo></mrow><mo>,</mo></mrow></mtd><mtd><mrow><mover><mi>d</mi><mo>&OverBar;</mo></mover><mo>/</mo><msub><mi>d</mi><mi>i</mi></msub><mo>&gt;</mo><mn>1</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>max</mi><mrow><mo>(</mo><msub><mi>w</mi><mi>i</mi></msub><mo>-</mo><mn>1</mn><mo>,</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mover><mi>d</mi><mo>&OverBar;</mo></mover><mo>/</mo><msub><mi>d</mi><mi>i</mi></msub><mo>&le;</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>n</mi></mrow>]]></math><img file="FDA0000951020530000013.GIF" wi="1246" he="166" /></maths>其中,M是图像间匹配对集合,m<sub>i</sub>为图像间的第i个匹配对,包括两个图像的对应点坐标,n为总的匹配对数目;W是每个匹配对的权值,初始值w<sub>i</sub>=1;P是每个匹配对的提取概率,p<sub>i</sub>的初始概率为1/n;T是判断每个匹配对在当前基础矩阵下是否为正确匹配对的判决函数,并且对所有错误匹配对之间的极对距离d<sub>i</sub>赋予一个常数c,这个常数在算法执行前进行确定;O是针对匹配对权重的调节函数,根据<img file="FDA0000951020530000014.GIF" wi="108" he="71" />的值对每个匹配对的权值进行调节,<img file="FDA0000951020530000015.GIF" wi="42" he="62" />为平均极对距离,等于d<sub>i</sub>的平均值,round为四舍五入函数,H为最大奖励程度;Inmost是用来记录当前最多正确匹配对数量,初始值为0;Mo用来存储当前最多正确匹配对数所对应的基础矩阵,初始值为3×3全零矩阵;N是算法最大迭代次数,当迭代次数大于N时,将停止迭代,将N设置为1000;Stc是概率自动机中的迭代停止条件,利用匹配对抽取概率的变化情况作为算法的终止条件;步骤三,从匹配对中按照概率随机抽取m个匹配对,并且按照m点法求出图像间的基础矩阵,然后将此基础矩阵带入所有匹配对,根据概率自动机中的判决函数T求出每个匹配对之间的极对距离d<sub>i</sub>和平均极对距离<img file="FDA0000951020530000016.GIF" wi="58" he="62" />步骤四,根据概率自动机中的权重调节函数O、每个匹配对的极对距离d<sub>i</sub>和平均极对距离<img file="FDA0000951020530000017.GIF" wi="41" he="62" />对每个匹配对进行权重的调整,并且根据调整后的权重计算每个匹配对的提取概率P;步骤五,判断当前基础矩阵是否具有最多正确匹配对数目,如果没有,转步骤三;如果有,且满足下面条件之一便结束算法迭代:条件1:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>K</mi><mo>&GreaterEqual;</mo><mfrac><mrow><mi>l</mi><mi>o</mi><mi>g</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>p</mi><mo>)</mo></mrow></mrow><mrow><mi>l</mi><mi>o</mi><mi>g</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mrow><mi>I</mi><mi>n</mi><mi>m</mi><mi>o</mi><mi>s</mi><mi>t</mi><mo>/</mo><mi>n</mi></mrow><mo>)</mo></mrow><mi>m</mi></msup><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000951020530000021.GIF" wi="533" he="135" /></maths>其中,p为置信概率,m为每次抽取匹配对的数目;条件2:概率自动机中的迭代停止条件Stc,以匹配对抽取概率的变化情况作为算法的终止条件,其表达式为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mi>K</mi><mo>-</mo><mi>L</mi><mo>+</mo><mn>1</mn></mrow><mi>K</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>|</mo><mrow><mo>(</mo><msubsup><mi>p</mi><mi>i</mi><mi>j</mi></msubsup><mo>-</mo><msubsup><mi>p</mi><mi>i</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>|</mo></mrow><mi>L</mi></mfrac><mo>&le;</mo><mi>&lambda;</mi></mrow>]]></math><img file="FDA0000951020530000022.GIF" wi="531" he="199" /></maths>其中,K为当前迭代次数;L为步长设定为10;λ为概率变化阈值,为一个正数,设定λ为0.01;条件3:当前迭代次数K等于自动机最大迭代次数N,其表达式为:K=N如果以上条件全不满足,则转步骤三。
地址 100124 北京市朝阳区平乐园100号