《图像配准ppt课件.ppt》由会员分享,可在线阅读,更多相关《图像配准ppt课件.ppt(51页珍藏版)》请在三一办公上搜索。
1、1,图像配准技术,郑雪梅,2,内容概要,1.图像配准介绍,4.图像配准的主要方法,5.部分配准算法的实验结果,2.图像配准的分类,3.图像配准的一般步骤,6.图像配准的评价,3,图像配准介绍-图像配准的定义,什么是图像配准?图像配准可以视为源图像和目标图像关于空间和灰度的映射关系:其中:f 为二维空间坐标变换(如仿射变换),g 为一维亮度或其他度量值变换。图像配准问题的关键:最佳空间变换。图像配准的实质:不同图像中表征同一位置的物理点一一对应。,4,图像配准是处理多源图像信息融合、目标识别等实际问题中的一个重要步骤。其广泛应用于军事、医学、航空等多个重要科研领域。经过图像配准,可以获得质量更高
2、、清晰度更好、定位更准确的目标信息。,图像配准介绍-图像配准的意义,5,图像配准介绍-图像配准的应用领域,计算机视觉-视频监控-对跟踪的目标区域进行配准-人脸识别,6,图像配准介绍-图像配准的应用领域,医学-不同模态的(CT,MRT)-配准后进行融合可以得到 更多的信息-单模态-同一病人不同时间,不同病人之间的,7,图像配准介绍-图像配准的应用领域,遥感-不同时间、不同视角,不同传感器-信息融合、环境监视、图像拼接、天气预报以及地图 更新等,8,图像配准介绍-图像配准的应用领域,军事-变化检测、运动目标指示、视点稳定、大视场重构、多通道融合、地形定位和导航,9,图像配准的分类,按图像的维数分类
3、-2D/2D:平面图像之间的配准-2D/3D:平面图像与空间图像之间的配准(医学上的体数据与二维 图像数据的配准)-3D/3D:两幅三维空间图像之间的配准 如果在空间维数的基础上再加上时间维,则原来的2D、3D就分别变成了3D、4D,这在临床医学上比较常见,主要是用来观察儿童骨骼发育,跟踪肿瘤变换等。按成像模式分类-单模:待配准的两幅图像是用同一种成像设备获取的-多模:待配准的两幅图像来源于不同的成像设备,10,图像配准的分类,按交互性分类-人工配准:完全由人工凭借经验进行,输入计算机后实现的只 是显示工作,不需要复杂的配准算法。-半自动:由人工给出一定的初始条件,如人工勾画出轮廓、控制优化参
4、数。-全自动:全自动化配准不需人工干预,由计算机自动完成。按对图像信息的利用情况分类-基于灰度:直接利用局部图像的灰度信息建立两幅图像之间的相似性度 量,然后采用某种搜索方法寻找使相似性度量值最大或最小的 变换模型的参数值。-基于特征:提取各类图像中保持不变特征如边缘点,闭区域的中心 等作为两幅图像配准的参考信息。,11,图像配准的分类,按图像配准的应用领域分类-军事,医学,遥感,计算机视觉按配准图像的来源分类-不同视点的图像配准(多视点分析):目的在于获取更大范围的 二维视角或重构被摄场景的三维信息。应用实例:遥感被摄区域图像镶嵌、计算机视觉形状恢复;-不同时间的图像配准(多时段分析):目的
5、在于寻找并度量两幅 不同时间内拍摄的图像中场景的变化。应用实例:遥感区域规划、计算机视觉运动跟踪、医学成像肿瘤病变检测、白天和黑夜的场景以及夏天和冬天的图像特征研究;,12,图像配准的分类,按配准图像的来源分类-不同传感器的图像配准(多模式分析):目的在于融合不同传感器信息,以获取更为丰富细致的场景信息。如可见光和红外图像配准、医学成像CT和MRI、多波段的人脸识别;-场景与场景模型图像配准:场景模型多为场景的数字表示,如地理信息系统(GIS)中的数字高程图等。配准场景与场景模型的目的在于在场景中定位模型,或在模型中定位场景,或者对场景与模型作相应的比较。应用实例:遥感将航片或卫片与地图或GI
6、S相配准、计算机视觉匹配模板图像与实时场景、医学成像将数字解剖图与病人的图片相比照。,13,配准算法的一般步骤,特征提取 特征匹配 估计变换模型 图像重采样及变换,14,配准算法的一般步骤特征提取,特征提取是指分别提取两幅图像中共有的图像特征。这种特征是出现在两幅图像中对比例、旋转、平移等变换保持一致性的特征,如线交叉点、物体边缘角点、虚圆闭区域的中心等可提取的特征。特征包括:点、线和面三类。-点特征是最常采用的一种图像特征,包括物体边缘点、角点、线交叉点等;根据各特征点的兴趣值将特征点分成几个等级。对不同的目的,特征点的提取应有所不同。-线特征是图像中明显的线段特征,如道路河流的边缘,目标的
7、轮廓线等。线特征的提取一般分两步进行:首先采用某种算法提取出图像中明显的线段信息,然后利用限制条件筛选出满足条件的线段作为线特征;-面特征是指利用图像中明显的区域信息作为特征。在实际的应用中最后可能也是利用区域的重心或圆的圆心点等作为特征。,15,配准算法的一般步骤特征提取,点特征提取方法-Harris(Harris Corner Detector)算法 受信号处理中相关函数的启发,给出与自相关函数相联系的矩阵M,M矩阵的特征值是自相关函数的一阶曲率,如果两个曲率值都高,那么就认为该点是角点,此方法对图像旋转、亮度变化、视角变化和噪声的影响具有很好的鲁棒性。,16,配准算法的一般步骤特征提取,
8、点特征提取方法-Susan(Susan Corner Detector)算法 SUSAN算法使用一个圆形的模板在图像上滑动,将位于圆形模板中心的待检测的像素点称为核心点。假设图像为非纹理,核心点的邻域被划分为两个区域:其一为亮度值等于(或相似于)核心点亮度的区域,称为核值相似区(Univalue Segment Assimilating Nucleus,USAN),其二为亮度值不相似于核心点亮度的区域。,17,配准算法的一般步骤特征提取,点特征提取方法-Harris-Laplace Harris算子能最稳定地在图像旋转、光照变化、透视变换条件下提取二维平面特征点,但在三维尺度空间中,Harri
9、s探测子的重复探测性能不好,不同尺度Harris特征点存在位置误差,Harris探测子不具有尺度和仿射不变性。而三维尺度空间中最稳定高效的特征尺度探测算子是归一化的Laplace算子。K.Mikolajczyk和C.Schmid结合了Harris和Laplace算子的优点,提出了Harris-Laplace算子。Harris-Laplace特征点具有尺度和旋转不变的特性,且对光照变换和小范围视角变化具有稳定性。,18,配准算法的一般步骤特征提取,点特征提取方法-SIFT特征点提取 使用Difference of Gaussian(DoG)filter来建立尺度空间。在尺度空间上提取极值点。,1
10、9,配准算法的一般步骤特征提取,点特征提取方法-SURF特征点提取 基于Hessian矩阵,它依靠Hessian矩阵行列式的局部最大值定位兴趣点位置。对于图像I中的某点X在尺度空间上的Hessian矩阵定义为:,其中,表示高斯二阶偏导在X处与图像I的卷积。、具有相似的含义。,20,配准算法的一般步骤特征提取,点特征提取方法-FAST特征点提取 FAST(Features from Accelerated Segment Test)是对SUSAN角点提取算法的简化,通过比较一个圆上16个像素点与中心像素点的灰度对比情况来判断中心像素点是否为特征点。,21,配准算法的一般步骤特征提取,点特征提取方
11、法-Harris-Affine-Hessian-Laplace-Hessian-Affine-Moravec算子-Forstner算子,22,配准算法的一般步骤特征提取,线特征提取方法-Robert-Sobel-Prewitt-Kirsch-Gauss-Laplace-Canny,23,配准算法的一般步骤特征提取,面特征提取方法-Mser 使用不同的阈值对图像进行二值化,这个过程中,所有阈值图像上形成的连接区域都是极小值区域,在阈值图像的变化过程中,形成了一系列嵌套的极值区域组。在每组嵌套区域里,有一类性质较为稳定的区域,这类区域在较大阈值范围内具有较小的变化,被定义为“最稳定极值区域”。,2
12、4,配准算法的一般步骤特征匹配,特征匹配 特征匹配分两步:对特征作描述 现有的主要特征描述子:SIFT特征描述子,SURF特征描述 子,对比度直方图(CCH),DAISY特征描述子,矩方法。利用相似度准则进行特征匹配 常用的相似性测度准则有如欧氏距离、马氏距离、Hausdorff 距离等。,25,配准算法的一般步骤特征匹配,特征描述-SIFT特征描述子 主要思想:一种基于图像梯度分布的特征描述子。特点:抗干扰性好。但维数高,计算复杂度大。,26,配准算法的一般步骤特征匹配,特征描述-SURF特征描述子 主要思想:将特征点的周围区域分成几个子区域,用每个子区域内像素点的X,Y方向的偏导和及其绝对
13、值的和组成特征点的描述子。特点:有较好的抗亮度变化能力。但是该描述子要求使用积分图像,限定了其应用范围。,27,配准算法的一般步骤特征匹配,特征描述-对比度直方图 主要思想:将特征点周围区域的像素点与特征点的对比度形成直方图来描述该特征点。特点:该方法比基于梯度的描述子要快。但描述力比基于梯度的要略微弱一点。,28,配准算法的一般步骤特征匹配,特征描述-DAISY特征描述子 主要思想:受SIFT算法和GLOH算法的启发,将梯度加权和用几个高斯方向偏导滤波器与原图像进行积分代替。特点:该描述子有和SIFT特征算子相似的优点,但是速度比SIFT特征算子要快。,29,配准算法的一般步骤特征匹配,特征
14、描述-矩方法 Hu矩,zernike矩,30,配准算法的一般步骤估计变换模型,空间变换模型是所有配准技术中需要考虑的一个重要因素,各种配准技术都要建立自己的变换模型,变换空间的选取与图像的变形特性有关。常用的空间变换模型有:刚体变换、仿射变换、投影变换、非线性变换。,1.刚体变换模型刚体变换是平移、旋转与缩放的组合,适用于具有相同视角,但拍摄位置不同的来自同一传感器的两幅图像的配准。刚体变换模型下,若点,分别为参考图像和待配准图像中对应的两点,则它们之间满足以下关系:,31,配准算法的一般步骤估计变换模型,2.仿射变换模型如果第一幅图像中的一条直线经过变换后,映射到第二幅图像上仍然为直线,且平
15、行直线仍旧被映射为平行直线,这样的变换称为仿射变换。该变换保持直线间的平行关系,但由于引入了缩放参数,故它不能保持直线段的长度和角度,若点,分别为参考图像和待配准图像中对应的两点,则它们之间满足以下关系:,32,配准算法的一般步骤估计变换模型,3.投影变换模型如果第一幅图像中的一条直线经过变换后,映射到第二幅图像上依然为直线,但平行关系不再保持,则称这样的变换称为投影变换。投影变换具有8个参数,可以将成像设备的运动(如平移、旋转、缩放等)描述得更为全面。若点,分别为参考图像和待配准图像中对应的两点,则它们之间满足以下关系:,33,配准算法的一般步骤估计变换模型,4.非线性变换模型 若第一幅图像
16、中的一条直线经变换后,映射至第二幅图像上不再是直线,我们把这样的变换称为非线性变换。在二维空间中,点 经非线性变换至点 变换公式为:F表示把第一幅图像映射到第二幅图像上的任意一种函数形式。典型的非线性变换如多项式变换,在2D空间中,多项式函数可写成如下形式:,34,配准算法的一般步骤图像重采样及变换,在得到两幅图像间的变换参数后,要将输入图像做相应参数的变换,使之与参考图像处于同一坐标系下,则校正后的输入图像与参考图像可用作后续的图像融合、目标变化检测处理或图像镶嵌;涉及输入图像变换后所得点坐标不一定为整像素数,则应进行插值处理。常用的插值算法有最近邻域法、双线性插值法和立方卷积插值法。,35
17、,配准算法的主要方法,基于灰度的-互相关-互信息-对数极坐标基于特征的-基于点-基于边缘(线)-基于区域,36,配准算法的主要方法,基于点-SIFT算法 SIFT特征匹配算法包括两个阶段,SIFT特征的生成与SIFT特征向量的匹配。SIFT特征向量的生成算法共包括4步:尺度空间极值检测,以初步确定关键点位置和所在尺度;拟和三维二次函数精确确定位置和尺度,同时去除低对比度的关键点和不稳定的边缘响应点;利用关键点邻域像素的梯度方向分布特性为每个关键点指定方向参数,使算子具备旋转不变性;生成SIFT特征向量;SIFT特征向量的匹配 对图像1中的某个关键点,找出其与图像2中欧式距离最近的前两个关键点的
18、距离NN和SCN,如果NN/SCN小于某个比例阈值,则接受这一对匹配点。,37,配准算法的主要方法,基于点-SIFT算法的几种扩展PCA-SIFT PCA-SIFT与标准sift有相同的亚像素位置(sub-pixel),尺度(scale)和主方向(dominant orientations),但在第4步计算描述子的时候,它用特征点周围的4141的像元计算它的主元,并用PCA-SIFT将原来的23939维的向量降成20维,以达到更精确的表示方式。GLOH(Gradient location-orientation histogram)把原来SIFT中44棋盘格的location bins改成仿射
19、状的同心圆的17 location bins来表示,并计算其中的梯度方向直方图(梯度方向分为16种),因此共1617272维,之后再作PCA将其降成128维,因此保有跟SIFT一样精简的表示方法。,38,配准算法的主要方法,基于点-SIFT算法的几种扩展 ASIFT(Affine-SIFT)通过原始图像来模拟得到场景在各个视角下的图像,再对这些得到的图像提起SIFT特征点,然后进行匹配。其仿射性要好于SIFT,具有全仿射不变性。,39,配准算法的主要方法,基于点-SUFR算法积分图像 对于积分图像中某X,X=(x,y),该点的值表示:原始图像中,原点和点X形成的矩形区域里面所有像素值之和。框装
20、滤波器 在surf算法中用框装滤波器来近似高斯二阶导,40,配准算法的主要方法,基于点-SUFR算法 特征点提取计算原图像的积分图像;用不同尺寸的框状滤波器来计算不同阶以及不同层上的每个点图像点的行列式。一般计算4阶4层;在3维(x、y、S)尺度空间中,在每个333的局部区域里,进行非最大值抑制。只有比邻近的26个点的响应值都大的点才被选为兴趣点;特征匹配 对图像1中的某个关键点,找出其与图像2中欧式距离最近的前两个关键点的距离NN和SCN,如果NN/SCN小于某个比例阈值,则接受这一对匹配点。,41,配准算法的主要方法,基于边缘-基本思想用边缘检测算子提取出边缘边缘匹配根据边缘的相似度提取边
21、缘上的控制点,如曲率比较大的点等,然后用这些点来进行匹配。将边缘拟合成直线,然后匹配直线。估计变换参数 用边缘上的控制点,或直线的端点,中点等等。,42,部分配准算法实验结果,特征点提取时间对比,43,部分配准算法实验结果,匹配效果对比,44,部分配准算法实验结果,拼接效果对比,45,部分配准算法实验结果,拼接效果对比,46,部分配准算法实验结果,实验结果总结对于医学图像,常用的方法是互信息,因为其特征不明显。对于模板匹配,常用的方法是互相关。当两幅图像之间只有简单的旋转,平移和缩放,则可以用对数极坐标算法。当待配准的图像特征十分明显时可以用SURF算法当待配准的图像特征信息比较弱一些的话可以
22、用SIFT算法。当待配准的图像之间存在较大的视角变换的话可以用ASIFT算法。,47,配准算法的评价,精度 配准过程中很容易引入各种各样的误差,而且很难区分是由配准算法引起的,还是由图像间的固有差异引起的。在评估配准精度时,主要将误差分为3类:位置误差、匹配误差和对齐误差。位置误差是指由不精确检测引起的控制点坐标偏移。匹配误差则是指在候选控制点之间建立匹配关系时误匹配的控制点对数目。对齐误差是指配准过程中采用的变换模型和图像真实畸变(包括比例缩放、旋转、平移以及传感器影响等)之间的差异。鲁棒性 鲁棒性是指如果让输入图像有一点小的变动,配准算法还能收敛到相同的结果。,48,配准算法的评价,自动化 自动化指配准算法的自动化执行程度,包括人工、半自动、自动3种形式。实时性 只有满足实时性要求,配准结果才具有实际应用价值。实时性研究涉及资源需求和算法的复杂度。可靠性 在实际应用中,要注意满足不同算法的前提约束条件。若给定算法所需要的输入条件,配准结果应该和理论上计算的一样,即算法具有较强的容错性。,49,配准算法的仍存在的问题,精度实时性通用性,50,附录,局部不变特征对比,51,谢谢!,