[发明专利]一种基于仿射参数估计的图像配准方法无效

专利信息
申请号: 201010263073.1 申请日: 2010-08-25
公开(公告)号: CN101916445A 公开(公告)日: 2010-12-15
发明(设计)人: 何凯;远中文;牟聪翀;卓磊 申请(专利权)人: 天津大学
主分类号: G06T7/00 分类号: G06T7/00
代理公司: 天津市北洋有限责任专利代理事务所 12201 代理人: 程毓英
地址: 300072 天*** 国省代码: 天津;12
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明属于计算机图像处理领域,涉及一种基于仿射参数估计的图像配准方法,包括下列步骤:采用SIFT算法提取图像的特征点,然后对不同图像的特征点进行匹配;根据仿射变换的具体情况,建立数学模型,并估计其中的仿射变换参数;采用估算出的图像间的仿射变换参数,进行运动补偿,得到配准后图像。本发明的图像配准方法,对误匹配具有较强的鲁棒性,实现起来比较简单并且有较好的实时性。
搜索关键词: 一种 基于 参数估计 图像 方法
【主权项】:
1.一种基于仿射参数估计的图像配准方法,包括:第一步:采用SIFT算法提取图像的特征点,然后对不同图像的特征点进行匹配,步骤如下:1)检测尺度空间极值点首先对原图和待匹配的图像分别生成一个尺度空间,然后进行采样,并将每一个采样点和它所有的相邻点比较以得到空间极值点;2)精确定位极值点设定阈值以去除低对比度的极值点和产生不稳定边缘响应的极值点,得到关键点,并通过拟和三维二次函数确定关键点的位置和尺度;3)确定每个关键点方向参数利用关键点邻域像素的梯度方向分布特性为每个关键点指定方向参数;4)关键点描述子的生成对原图和待匹配的图像的每个关键点,首先将坐标轴旋转为关键点的方向,以关键点为中心取8×8的窗口,将8×8窗口分成4个小块,在每个4×4的小块上计算8个方向的梯度方向直方图,绘制每个梯度方向的累加值,即可形成一个种子点,从而生成原图和待匹配的图像各自的SIFT特征向量;设定一个比例阈值,取原图中的某个关键点,并找出其与待配准图像中欧式距离最近的前两个关键点,在这两个关键点中,如果最近的距离除以次近的距离少于某个比例阈值,则接受这一对匹配点,用这种方法即可得到原图中的点列和待配准图像中对应的点列;第二步:根据仿射变换的具体情况,建立数学模型如下:已知待配准图像中的点列(xi,yj)(i=1,2,ΛM,j=1,2.ΛN)和原图像中的点列(Xi,Yj)(i=1,2,ΛM,j=1,2.ΛN)(M,N分别为图像的长和宽)和原图像中初始点X0的坐标(X0,Y0)的值,令它们近似满足:xiyi=f1(Xi,Yi)g1(Xi,Yi)=cosθ-sinθsinθcosθkx00kyXi-X0Yi-Y0+x0y0]]>其中(i=1,2,ΛN),θ∈(-π,π],kx>0,ky>0,估计其中的仿射变换参数θ,kx,ky,x0,y0,使得误差函数)达到最小,其中,θ为图像旋转角度,kx,ky分别为图像在x轴,y轴方向的畸变,x0,y0分别为图像在x轴,y轴方向的位移,记X=(θ,kx,ky,x0,y0),该算法步骤如下:1)选定初始点X0,计算F的梯度令k=0;2)计算3)如果Hk接近奇异,取计算:Xk+1=Xk+tPk;否则,由方程HkPk=-gk解出Pk,进行如下判断:若F(Xk+Pk)<F(Xk),令Xk+1=Xk+Pk,若F(Xk+Pk)≥F(Xk),令计算Xk+1=Xk+tPk;4)计算‖gk+1‖,判断其是否满足迭代终止条件‖gk+1‖<ε,如满足,则Xk+1即为仿射变换参数,否则,k=k+1,重复上述过程直到满足迭代终止条件为止;第三步:采用估算出的图像间的仿射变换参数,进行运动补偿,得到配准后图像。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于天津大学,未经天津大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201010263073.1/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top