【图像安全性分析】密钥敏感性分析MATLAB代码

%读取更改原图一个像素点后的加密图

image1=imread('改变像素点加密后的lena1.bmp','bmp');

%读取没有更改原图的加密图

image2=imread('加密后的lena.bmp','bmp');

%显示图片

subplot(1,2,1),imshow(image2),title('原图像')

subplot(1,2,2),imshow(image1),title('改变后图像')

%计算相同位置灰度值相等的个数

[M,N]=size(image1);

m=0;

u1=0;

for i=1:M

    for j=1:N

        u1=u1+abs(double(image1(i,j))-double(image2(i,j)))/255;

        if image1(i,j)~=image2(i,j)

            m=m+1;

        end

    end

end

ua=sum(abs(double(image1(:))-double(image2(:)))/255);%简单的式子，一个就搞定

%计算NPCR

NPCR1=(m*100)/(M*N);

NPCR=sum(double(image1(:))~=double(image2(:)))*100/(M*N);%简单的方法

UACI1=(u1*100)/(M*N);

UACI=(ua*100)/(M*N);

disp(['NPCR值为：',num2str(NPCR)]);

disp(['NPCR值为：',num2str(NPCR1)]);

disp(['UACI值为：',num2str(UACI)]);%方法一

disp(['UACI值为：',num2str(UACI1)]);%方法二

（用了两种方法）像素变化率 NPCR（the number of pixels change rate）和归一化平均变化强度UACI（the unified average changing intensity）。其中 NPCR 表示的不同密文图像在相同位置上灰度值互不相同的比率，而 UACI 则表示不同密文图像之间的平均变化密度，通常用于图像加密性能分析