计算两张图片的相似度-余弦相似度法

理论参考博客：https://blog.csdn.net/zz_dd_yy/article/details/51926305

代码参考博客：https://blog.csdn.net/qq_23100787/article/details/51175775

                         https://blog.csdn.net/qq_23100787/article/details/51175877

%余弦相似度算法求两张图片的相似度
picture1 = imread('C:\Users\Chris\Desktop\test.jpg');
picture2 = imread('C:\Users\Chris\Desktop\test.jpg');
picture1 = rgb2gray(picture1);%图像灰度化处理
picture2 = rgb2gray(picture2);
% figure;
% imshow(picture1);
% figure;
% imshow(picture2);
t1=picture1;
[a1,b1]=size(t1);
t2=picture2;
t2=imresize(t2,[a1 b1],'bicubic');%将两张图片缩放为一致大小
t1=round(t1);
t2=round(t2);
e1=zeros(1,256);
e2=zeros(1,256);
%获取直方图分布/
for i=1:a1
    for j=1:b1
        m1=t1(i,j)+1;
        m2=t2(i,j)+1;
        e1(m1)=e1(m1)+1;;%读入灰度值
        e2(m2)=e2(m2)+1;
    end
end
figure;
imhist(uint8(t1));
figure;
imhist(uint8(t2));
%将直方图分为64个区
m1=zeros(1,64);
m2=zeros(1,64);
for i=0:63
    m1(1,i+1)=e1(4*i+1)+e1(4*i+2)+e1(4*i+3)+e1(4*i+4);
    m2(1,i+1)=e2(4*i+1)+e2(4*i+2)+e2(4*i+3)+e2(4*i+4);
end
%计算余弦相似度
A=sqrt(sum(sum(m1.^2)));
B=sqrt(sum(sum(m2.^2)));
C=sum(sum(m1.*m2));
cos1=C/(A*B);%计算余弦值
cos2=acos(cos1);%弧度
v=cos2*180/pi;%换算成角度
figure;
imshow(uint8([t1,t2]));
title(['余弦值为：',num2str(cos1),'       ','余弦夹角为：',num2str(v),'°']);
%完