zssd.m

harris角点检测算子

N=1;
[ligne1,col1]= size(posc1);
[ligne2,col2]= size(posc2);

for i=1:ligne1
    somme = 0;
    for l=-N:1:N
        for k=-N:1:N
            somme=somme+double(ori_im1(posr1(i)+l,posc1(i)+k));
        end
    end
    moyenne1(i)=somme/(2*N+1)^2;
end

for j=1:ligne2
    somme =0;
    for l=-N:1:N
        for k=-N:1:N
            somme=somme+double(ori_im2(posr2(j)+l,posc2(j)+k));
        end
    end
    moyenne2(j)=somme/(2*N+1)^2;
end

for i=1:ligne1
    for j=1:ligne2-1
        ZSSD1 = 0;
        for l=-N:1:N
            for k=-N:1:N
                ZSSD1=ZSSD1+((double(ori_im1(posr1(i)+l,posc1(i)+k))-double(moyenne1(i))) - (double(ori_im2(posr2(j)+l,posc2(j)+k))-double(moyenne2(j))))^2;
            end
        end
        ZSSD(i,j)=ZSSD1;
    end
end

for i = 1:ligne1
    pt = find(ZSSD(i,:) == min(ZSSD(i,:)));
    figure(3)
    imshow(ori_im1);
    hold on;
    plot(posc1(i),posr1(i),'r+')
    figure(4)
    imshow(ori_im2);
    hold on;
    plot(posc2(pt),posr2(pt),'r+')
    pause
end