matlab程序 surf算法,【求大神帮忙，surf算法源代码解析】

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

% Example, Corresponding points

% Load images

I1=imread('1.JPG');

I2=imread('2.JPG');

% Get the Key Points

Options.upright=true;

Options.tresh=0.0001;

Ipts1=OpenSurf(I1,Options);

Ipts2=OpenSurf(I2,Options);

% Put the landmark descriptors in a matrix

D1 = reshape([Ipts1.descriptor],64,[]);

D2 = reshape([Ipts2.descriptor],64,[]);

% Find the best matches

err=zeros(1,length(Ipts1));

cor1=1:length(Ipts1);

cor2=zeros(1,length(Ipts1));

for i=1:length(Ipts1),

distance=sum((D2-repmat(D1(:,i),[1 length(Ipts2)])).^2,1);

[err(i),cor2(i)]=min(distance);

end

% Sort matches on vector distance

[err, ind]=sort(err);

cor1=cor1(ind);

cor2=cor2(ind);

% Show both images

I = zeros([size(I1,1) size(I1,2)*2 size(I1,3)]);

I(:,1:size(I1,2),:)=I1; I(:,size(I1,2)+1:size(I1,2)+size(I2,2),:)=I2;

figure, imshow(I/255); hold on;

% Show the best matches

for i=1:30,

c=rand(1,3);

plot([Ipts1(cor1(i)).x Ipts2(cor2(i)).x+size(I1,2)],[Ipts1(cor1(i)).y Ipts2(cor2(i)).y],'-','Color',c)

plot([Ipts1(cor1(i)).x Ipts2(cor2(i)).x+size(I1,2)],[Ipts1(cor1(i)).y Ipts2(cor2(i)).y],'o','Color',c)

end