水印评价 NC psnr

以matlab为基础 所编写的水印图像所需的评价标准

clear;
close all;
I=imread（‘lena‘‘bmp‘）;%I为原图像 
% I1=wiener2（I[3 3]）;
% L=medfilt2（I1[3 3]）;
% I2=imadjust（L[0.1 0.9][0 1]1.2）;%I2为变换后图像
I2=imread（‘lenaq0‘‘bmp‘）
%以上I是原图像，I2是变换后的图像，具体变换可以自定义
%下面是均方误差mse峰值信噪比psnr和熵entropy的计算
X=I;
Y=I2;
mse_m=double（zeros（256256））; %均方误差
for i=1:256
for j=1:256
mse_m（ij）=（X（ij）-Y（ij））^2;
end
end
mse=sum（mse_m（1:256））/256*256;
psnr=10*log（double（255*255/mse））/log（10）;%峰值信噪比
% % 下面计算熵
% t（1:256）=0; 
% %给每个灰度值计数
% for i=1:180
% for j=1:180
% t（Y（ij）+1）=t（Y（ij）+1）+1;
% end
% end
% %计算各个灰度出现频率
% for i=1:256
% t（i）=t（i）/（256*256）;
% end
% entropy=0; %计算熵
% for i=1:256
% if t（i）>0
% entropy=entropy-t（i）*log（t（i））/log（2）;
% end
% end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

    I.A....       802  2013-03-31 17:15  水印评价代码\pingjia.m

     文件        884  2011-04-07 07:55  水印评价代码\PSNR&NC.txt

    I.A....       419  2013-03-31 20:42  水印评价代码\PSNR和归一化相关系数NC值.txt

     目录          0  2013-03-31 20:51  水印评价代码

----------- ---------  ---------- -----  ----

                 2105                    4