% Script belongs to Asad Ali
% GIK Institute of Engineering Sciences & Technology, Pakistan
% Email: asad_82@yahoo.com
function [similarity] = CalculateAffinity(img1, img2, img3)


counter2 = 1;
for j = 1:size(img1,1)
    for i = 1:size(img1,1)
        counter1 = 1;
        for j2= 1:size(img1,1)
            for i2 = 1:size(img1,1)
                dist = sqrt((img1(i,j) - img1(i2,j2))*(img1(i,j) - img1(i2,j2)) ...
                + (img2(i,j) - img2(i2,j2))*(img2(i,j) - img2(i2,j2)) + ...
                (img3(i,j) - img3(i2,j2))*(img3(i,j) - img3(i2,j2)));
                similarity(counter2,counter1) = exp(-dist/(2));
                counter1 = counter1 + 1;
            end
        end
        counter2 = counter2 + 1;
    end
end

Published with MATLAB® R2018a