% function SSIM()
% SSIM 结构相似性图像处理
clear
close all
K1=0.0001;K2=0.0001;a=1;b=1;c=1; % 参数设置
X=imread('2.jpg');
Y=imread('4.jpg');
% X=rgb2hsi(X);
% Y=rgb2hsi(Y);
% X=double(X)./255; % 灰度0-1之间的图像
% Y=double(Y)./255;
[L1,L2]=size(X);
ux=mean2(X); % 亮度平均
uy=mean2(Y);
deltax=sqrt(sum(sum((X-ux).^2))/(L1*L2-1)); % 对比度
deltay=sqrt(sum(sum((Y-uy).^2))/(L1*L2-1));
% sX=double(X-ux)./deltax;
% sY=double(Y-uy)./deltay;
C1=(K1*255)^2; % C1是引入的参数,K1是一个远小于1的常数
l_XY=(2*ux*uy+C1)/(ux^2+uy^2+C1); % 亮度比较函数
C2=(K2*255)^2;
c_XY=(2*deltax*deltay+C2)/(deltax^2+deltay^2+C2); % 对比度比较函数
deltaxy=sum(sum((X-ux).*(Y-uy)))/(L1*L2-1);
% C3=(K3*255)^2;
C3=C2/2;
s_XY=(deltaxy+C3)/(deltax*deltay+C3); % 结构对比函数
SSIM_XY=l_XY^a * c_XY^b * s_XY^c % 信号X与Y的SSIM索引
评论0