基于SIFT特征的图像配准(Matlab源代码).docVIP

基于SIFT特征的图像配准（附Matlab源代码） 本文先给出了采用SIFT方法进行图像配准的实验原图以及实验结果，最后附上Matlab源代码。 实验一： 实验一的图像（见图1.1）是本人自己拍摄的，然后由软件裁剪成400×400像素而成，其中参考图像和待配准图像之间有重叠部分，且具有一定的旋转。这是一般难度的图像配准。 （图1.1 参考图像（左图）和待配准图像（右图）） 经过筛选后的能作为配准控制点的SIFT匹配特征点对见下图1.2所示。 （图1.2 经过筛选后的能作为配准控制点的SIFT匹配特征点对） 配准后的图像见下图1.3所示。 （1.3 配准后的图像） 实验二： 实验二的图像（见图2.1）是本人自己在同一地点不同角度拍摄的，然后由软件裁剪成600×450像素而成，其中参考图像和待配准图像之间有重叠部分，且具有一定的旋转。这是具有一定难度的图像配准。 （图2.1 参考图像（左图）和待配准图像（右图）） 经过筛选后的能作为配准控制点的SIFT匹配特征点对见下图2.2所示。 （图2.2 经过筛选后的能作为配准控制点的SIFT匹配特征点对） 配准后的图像见下图1.3所示。 （2.3 配准后的图像） 实验三： 实验三的图像（见图3.1）是两幅SAR图像，大小为400×400像素，其中参考图像和待配准图像之间有重叠部分，且具有一定的旋转。这是具有难度的图像配准。 （图3.1 参考图像（左图）和待配准图像（右图）） 经过筛选后的能作为配准控制点的SIFT匹配特征点对见下图3.2所示。 （图3.2 经过筛选后的能作为配准控制点的SIFT匹配特征点对） 配准后的图像见下图3.3所示。 （3.3 配准后的图像） 下面是Matlab源代码，共计13个*.m文件和1个*.exe文件，其中siftWin32 .exe文件需要自行下载，下载链接为 /s/1pJL4SB5，需要*.m格式的源代码的，到 /download/destiny0321/9204877进行下载。 Zoo_main.m close all;clear all;clc; im1=imread(image data\image3_1.jpg); im2=imread(image data\image3_2.jpg); gray1=zoo_x2gray(im1); gray2=zoo_x2gray(im2); [des1,loc1]=zoo_sift(gray1); [des2,loc2]=zoo_sift(gray2); figure;zoo_drawPoints(im1,loc1,im2,loc2); Num=3;Thresh=0.85; match=zoo_BidirectionalMatch(des1,des2,Num,Thresh); clear des1 des2 loc1=loc1(match(:,1),:); loc2=loc2(match(:,2),:); figure;zoo_linePoints(im1,loc1,im2,loc2); agl=zoo_getRotAgl(loc1,loc2); figure;zoo_drawRotAglHist(agl); opt=zoo_optIndex(agl); loc1=loc1(opt,:); loc2=loc2(opt,:); figure;zoo_linePoints(im1,loc1,im2,loc2); T=zoo_getTransMat(gray1,loc1,gray2,loc2); im=zoo_imRegist(im1,im2,T); figure,imshow(im); zoo_x2gray.m function gray=zoo_x2gray(im) if length(size(im))==3 gray=rgb2gray(im); else gray=im; end gray=uint8(medfilt2(double(gray))); zoo_sift.m function [des,loc]=zoo_sift(im) [row,col]=size(im); f=fopen(tmp.pgm,w); if f==-1 error(Could not create file tmp.pgm.); end fprintf(f, P5\n%d\n%d\n255\n, col, row); fwrite(f,im,uint8); fclose(f); if isunix command = !./sift ; else command = !siftWin32 ; end comman