[发明专利]基于边缘特征的图像自动配准及快速超分辨率融合方法无效
申请号: | 200810020435.7 | 申请日: | 2008-03-06 |
公开(公告)号: | CN101527039A | 公开(公告)日: | 2009-09-09 |
发明(设计)人: | 徐立中;凌静;石爱业;汤敏;黄凤辰;王慧斌;马贞立 | 申请(专利权)人: | 河海大学 |
主分类号: | G06T7/00 | 分类号: | G06T7/00;G06T5/00 |
代理公司: | 暂无信息 | 代理人: | 暂无信息 |
地址: | 21009*** | 国省代码: | 江苏;32 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 一种基于边缘特征的图像自动配准及快速超分辨率融合方法,涉及图像处理领域,本发明的特征在于具有以下处理步骤:(1)计算机从USB等接口读入待处理的两类目标图像至缓冲区,其中一类是低分辨率多光谱图像,另一类是高分辨率灰度图像;(2)采用小波变换和最小二乘拟合模型求取两类图像的配准参数;(3)将配准的多光谱图像用线性IHS变换获取强度分量图像,然后将高分辨率灰度图像直方图匹配到强度分量;(4)用直方图匹配后的高分辨率灰度图像替换强度分量,线性IHS逆变换实现图像超分辨率融合;(5)输出超分辨多光谱率彩色图像。本发明提供的方法复杂度和执行速度低于传统方法,占用计算机系统资源少,便于在嵌入式设备上执行。 | ||
搜索关键词: | 基于 边缘 特征 图像 自动 快速 分辨率 融合 方法 | ||
【主权项】:
1.一种基于边缘特征的图像自动配准及快速超分辨率融合方法,是由各类摄像设备完成待处理的两类图像,一类为低分辨率多光谱彩色图像,另一类为高分辨率的灰度图像,将目标图像数字化为便于通过USB、红外等接口输入计算机处理,处理后的结果图像输出至缓冲器,可以直接在本地进行存储,或直接输出显示,或通过网络存储设备进行远程传输,其特征在于它还包扩下述步骤:1)采用小波变换结合最小二乘拟合的模型来求取M和P的配准参数实现图像的配准。其中小波变换用来提取图像的边缘特征点,再使用仿射变换建立M图像的边缘特征点与P图像的边缘特征点的对应关系,最后依据最小二乘拟合求取仿射变换参数,从而将M图像配准到P图像,其特征为:a)小波变换提取图像的边缘利用小波变换提取图像边缘的算法概括如下:首先利用二进小波变换将图像分解为J层。针对于j(1≤j≤J)分解层的高频分量,求解小波模极大Wj,并设定一个阈值T,如果Wj>T则认为是边缘点;b)仿射变换考虑的图像配准可以认为是刚体变换,故此可以采用仿射变换来表征特征点之间的关系,仿射变换如下式:x y = a 11 a 12 a 21 a 22 x ′ y ′ + b 1 b 2 ]]> ①式中(x,y)表示参考图像的坐标,(x′,y′)表示待配准图像的坐标,a 11 a 12 a 21 a 22 ]]> 表示图像的旋转和缩放(b1,b2)表示图像在x、y方向的平移量;c)最小二乘拟合记A=(a11 a12 b1)B=(b21 b22 b2)Q = ( x ′ ) j 1 ( x ′ ) j 2 . . . ( x ′ ) jn ( y ′ ) j 1 ( y ′ ) j 2 . . . ( y ′ ) jn 1 1 . . . 1 ]]> ②Mx=(j1(x),j2(x),…,jn(x))My=(j1(y),j2(y),…,jn(y))由此建立最小二乘表达式如下:AQ=Mx ③BQ=My ④由上述最小二乘方程可以求解出A、B;2)图像的直方图匹配记经过配准的多光谱图像为MR,其共有三个波段分别为MR1,MR2,MR3。将图像P的直方图匹配到强度分量图像I,并记经过直方图匹配后的图像分别为P1;3)图像的超分辨率融合图像的超分辨率融合采用线性IHS逆变换,这主要考虑线性IHS变换较常规的IHS变换运算量小、速度快。具体超分辨率融合过程为:设I图像、P图像最大值分别为M1和M2,记M=max(M1,M2)。将I图像、P图像的像素值看作特征s,则在像素位置(m,n)处的信任函数分别记为h(sI,m,n)、h(sI,m,n)。令h(sI,m,n)=I(m,n)/M,h(sP,m,n)=P(m,n)/M依据模糊积分的单调性的约束条件,像素位置为(m,n)处,记s={s1(m,n),s2(m,n)},其中{s1(m,n),s2(m,n)}定义如下:s 1 ( m , n ) = I ( m , n ) , s 2 ( m , n ) = P ( m , n ) ifh ( s I ; m , n ) ≤ h ( s P ; m , n ) s 1 ( m , n ) = P ( m , n ) , s 2 ( m , n ) = I ( m , n ) ifh ( s I ; m , n ) > h ( s P ; m , n ) ]]> ⑤因此,模糊测度g(A1)=g({s1(m,n),s2(m,n)})=1 ⑥g(A2)=g({s2(m,n)}) ⑦只要确定g({s2(m,n)})就可以求出融合后的像素值,而g({s2(m,n)})就是图像I或图像P在(m,n)处像素值的模糊密度。为了表述方便,此处记g({(s2(m,n)})为g(m,n)。记模糊积分的结果为FI(m,n),融合后的像素值为P′(2j;m,n)P ′ ( m , n ) = FI ( m , n ) × Mifh ( s I ; m , n ) ≤ h ( s P ; m , n ) FI ( m , n ) × Mifh ( s I ; m , n ) > h ( s P ; m , n ) ]]>= I ( m , n ) + ( P ( m , n ) - I ( m , n ) ) × g ( m , n ) ifh ( s I ; m , n ) ≤ h ( s P ; m , n ) P ( m , n ) + ( I ( m , n ) - P ( m , n ) ) × g ( m , n ) ifh ( s I ; m , n ) > h ( s P ; m , n ) ]]> ⑧其中g(m,n)的大小可以根据对融合图像的属性要求而定,若关于图像P的g(m,n)越来越小,则对I图像的g(m,n)越来越大,则融合图像的光谱信息越来越好,由此可见,g(m,n)在中间取值范围内可以使得融合图像的光谱和空间细节信息最好;
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于河海大学,未经河海大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/200810020435.7/,转载请声明来源钻瓜专利网。
- 上一篇:一种工件裂纹的焊补工艺方法
- 下一篇:粉煤灰输送管道气流粉碎工艺
- 彩色图像和单色图像的图像处理
- 图像编码/图像解码方法以及图像编码/图像解码装置
- 图像处理装置、图像形成装置、图像读取装置、图像处理方法
- 图像解密方法、图像加密方法、图像解密装置、图像加密装置、图像解密程序以及图像加密程序
- 图像解密方法、图像加密方法、图像解密装置、图像加密装置、图像解密程序以及图像加密程序
- 图像编码方法、图像解码方法、图像编码装置、图像解码装置、图像编码程序以及图像解码程序
- 图像编码方法、图像解码方法、图像编码装置、图像解码装置、图像编码程序、以及图像解码程序
- 图像形成设备、图像形成系统和图像形成方法
- 图像编码装置、图像编码方法、图像编码程序、图像解码装置、图像解码方法及图像解码程序
- 图像编码装置、图像编码方法、图像编码程序、图像解码装置、图像解码方法及图像解码程序