图像的平均梯度AG，信息熵IE，峰值信噪比PSNR，平均值，SSIM（matlab）

1，平均梯度

I10=imread('F:\picture_AIEM\10_aiem.jpg');
x2=I10;
x2=x2(:,:,1);
x2=double(x2);
[m,n]=size(x2);
agvalue=0;
for i=1:m-1
for j=1:n-1
Ix=x2(i+1,j)-x2(i,j);
Iy=x2(i,j+1)-x2(i,j);
agvalue=agvalue+sqrt((Ix.^2+Iy.^2)/2);
end
end
AG1=agvalue/(m*n);
disp(['输入数据的AG为:',num2str(AG1)]);

2，信息熵

I1 = imread('F:\picture_AIEM\10_aiem.jpg');
e_1 = entropy(I1);

3，峰值信噪比

P1 = imread('F:\use\1.jpg');
I1 = imread('F:\picture_AIEM\1_aiem.jpg');
p1 = psnr(P1, I1);%P1和I1必须大小相同。P1是原图，I1是处理后的

4，平均值

子函数

function [a,average,deviation]=mean_deviation(I)
[m,n,q]=size(I);
I1=rgb2gray(I);
[height1,width1,channel1] = size(I1);
average=sum(I1(:))/(height1*width1);
a=sum(I(:))/(m*n*q);
I1 = double(I1);
deviation = 0;
for i = 1:height1
    for j = 1:width1
        deviation = deviation + (I1(i,j) - average)*(I1(i,j) - average);
    end
end
deviation = deviation/(height1*width1 - 1);    %图一的方差
deviation=sqrt(deviation);
end

主函数

I1=imread('F:\picture_AIEM\1_aiem.jpg');
[a1,~,~]=mean_deviation(I1);

5，ssim

结构相似性，比较亮度、对比度和结构三个指标

SSIM = ssim(srcImage,dstImage);%范围0-1，值越大越好，当两张图一样时为1