发明名称 一种基于区域特征的多源异类图像配准方法
摘要 本发明涉及图像处理领域中的图像配准方法,具体为一种基于区域特征的多源异类图像配准方法,对基准图像和待配准图像中进行聚类分割,得到两幅图像的分割区域;利用数字形态学从分割区域提取两幅图像的封闭轮廓;封闭轮廓删除,采用轮廓不变矩理论,引入矩变量距离均值、方差约束机制和一致性检查的匹配策略,获取最佳匹配轮廓对;以最佳封闭轮廓的质点为控制点,计算配准参数,实现两幅图像的配准。本发明克服了已有算法计算复杂、容易造成误匹配,导致配准精度低,甚至配准失败的缺陷,可应用于安全监测、机载多源图像融合、图像目标识别和跟踪等民用或军用领域。
申请公布号 CN103578105B 申请公布日期 2016.04.27
申请号 CN201310533490.7 申请日期 2013.11.01
申请人 中北大学 发明人 杨风暴;王志社;彭智浩;冯裴裴;牛涛
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 太原科卫专利事务所(普通合伙) 14100 代理人 朱源
主权项 一种基于区域特征的多源异类图像配准方法,其特征在于包括以下步骤:S1:对于基准图像和待配准图像分别进行k‑均值聚类处理,分割出基准图像和待配准图像的封闭区域,基准图像和待配准图像分别为可见光图像和SAR图像或SAR图像和可见光图像;S2:对基准图像和待配准图像分割后得到的封闭区域分别运用数学形态学处理,进行数学形态学先腐蚀后膨胀的开运算,腐蚀消除封闭区域中比腐蚀因子小的散点和噪声,膨胀把封闭区域中比膨胀因子小的缺口或孔洞填充上,把断裂的区域连接成封闭轮廓,获得目标和背景各区域的封闭边界,提取出基准图像和待配准图像的封闭轮廓;S3:删除基准图像和待配准图像的封闭轮廓中轮廓长轴长度小于设定阈值的封闭轮廓;S4:在S3提取到的基准图像和待配准图像的封闭轮廓中,采用轮廓不变矩理论,利用封闭轮廓的矩不变量将封闭轮廓进行匹配,获得基准图像和待配准图像的初始匹配轮廓对,轮廓不变矩理论中,一个封闭轮廓的二维区域G,用7个不变矩来表示,二维区域G中心矩定义为<img file="FDA0000879818840000011.GIF" wi="599" he="119" />式中,x,y表示二维区域G内的坐标,p、q分别代表x轴、y轴方向的指数,<img file="FDA0000879818840000012.GIF" wi="123" he="63" />为二维区域G的质心坐标,利用中心矩,构造的7个不变矩表示如下:<img file="FDA0000879818840000013.GIF" wi="389" he="79" /><img file="FDA0000879818840000014.GIF" wi="1333" he="109" /><img file="FDA0000879818840000015.GIF" wi="701" he="110" /><img file="FDA0000879818840000021.GIF" wi="1141" he="215" /><img file="FDA0000879818840000022.GIF" wi="870" he="182" /><img file="FDA0000879818840000023.GIF" wi="1150" he="214" />对于每一个封闭轮廓的二维区域G,可以把它的7个不变矩看作是七维特征空间的一个点,两个二维区域G的相似性根据这两个点的距离来度量,距离D<sub>ij</sub>可以表示为:<img file="FDA0000879818840000024.GIF" wi="573" he="156" />式中,i=1,2,…,m,j=1,2,…,n分别表示基准图像和待配准图像中的二维区域G的序号,对于基准图像中的m个二维区域和待配准图像中的n个二维区域,它们的关系可以用m×n的距离矩阵D来判定,距离矩阵D表示为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>D</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>D</mi><mn>11</mn></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>D</mi><mrow><mi>m</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>D</mi><mrow><mn>1</mn><mi>n</mi></mrow></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>D</mi><mrow><mi>m</mi><mi>n</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000879818840000025.GIF" wi="452" he="230" /></maths>选取距离矩阵D中每行元素最小值对应的匹配轮廓对,并设定一阈值D<sub>T</sub>,当选取的每行的匹配轮廓对的距离≤D<sub>T</sub>时,可以认为该匹配轮廓对满足匹配要求,得到基准图像和待配准图像的初始匹配轮廓对;S5:对基准图像和待配准图像的初始匹配轮廓对利用矩不变量的距离建立距离矩阵,计算矩不变量的距离方差矩阵,把距离方差矩阵中的每一行距离方差最小的两个轮廓作为求取配准参数的优化匹配轮廓对;S6:在基准图像和待配准图像的优化匹配轮廓对中,以封闭轮廓的质心作为控制点,求取配准参数为s、θ、△X、<img file="FDA0000879818840000026.GIF" wi="109" he="54" />s、θ、△X、<img file="FDA0000879818840000027.GIF" wi="78" he="55" />为两图像之间的缩放系数、旋转量、x轴平移量和y轴平移量,基准图像和待配准图像的坐标(X,Y)与(X',Y')之间满足如下关系:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><mi>X</mi></mtd></mtr><mtr><mtd><mi>Y</mi></mtd></mtr></mtable></mfenced><mo>=</mo><mi>s</mi><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mi>c</mi><mi>o</mi><mi>s</mi><mi>&theta;</mi></mrow></mtd><mtd><mrow><mi>s</mi><mi>i</mi><mi>n</mi><mi>&theta;</mi></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mi>s</mi><mi>i</mi><mi>n</mi><mi>&theta;</mi></mrow></mtd><mtd><mrow><mi>cos</mi><mi>&theta;</mi></mrow></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><msup><mi>X</mi><mo>&prime;</mo></msup></mtd></mtr><mtr><mtd><msup><mi>Y</mi><mo>&prime;</mo></msup></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mi>&Delta;</mi><mi>X</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>&Delta;</mi><mi>Y</mi></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000879818840000031.GIF" wi="785" he="161" /></maths>完成基准图像和待配准图像配准。
地址 030051 山西省太原市尖草坪区学院路3号