function czeka=Czekonowsky(imgo,imgc) img1=double(imread(imgo,'bmp')); img2=double(imread(imgc,'bmp')); [sizex sizey dim]=size(img1) Rimg=img1(:,:,1); Gimg=img1(:,:,2); Bimg=img1(:,:,3); RCimg=img2(:,:,1); GCimg=img2(:,:,2); BCimg=img2(:,:,3); %canberra=sum(sum(abs(Rimg-RCimg)./abs(Rimg+RCimg+1)))+sum(sum(abs(Gimg-GCimg)./abs(Gimg+GCimg+1)))+sum(sum(abs(Bimg-BCimg)./abs(Bimg+BCimg+1))); czeka=1-(2*sum(sum((min(Rimg,RCimg)+min(Gimg,GCimg)+min(Bimg,BCimg))))./sum(sum((Rimg+Gimg+Bimg+RCimg+GCimg+BCimg))));