renyiqujian.m

来自「这么代码是用Matlab编写的」· M 代码 · 共 63 行

clear all;

Ximage=double(imread('D:\work\Lena.jpg'));%待隐藏图像
Yimage=double(imread('D:\work\cameraman.jpg'));%公开图像
% Ih=size(Ximage,1);%图像的长度
% Iw=size(Ximage,2);%图像的宽度
Image=double(zeros(size(Ximage)));
Image=Ximage;
a=2;b=10;
t=1;

                                 %                             隐                            藏
% key1=4;
key2=[3,4,5];
s=size(key2,2);

 for(t=1:s)
figure(t)
subplot(2,2,1)
imshow(Image,[]);
title('待隐藏图像');
figure(t)
subplot(2,2,2)
imshow(Yimage,[]);
title('模板图像');
Image=(key2(t)-a)*Image/(b-a)+Yimage*(b-key2(t))/(b-a);
figure(t)
subplot(2,2,3)
imshow(Image,[]);
title('结果图像');
Image=uint8(Image);
figure(t)
subplot(2,2,4)
imhist(Image);
title('结果图像的直方图');
Image=double(Image);
 end
 PSNRy=PSNR(Image,Yimage);
 RMSEy=RMSE(Image,Yimage);
%                                                        恢                            复
for(t=1:s)
figure(t+s)
subplot(2,2,1)
imshow(Image,[]);
title('结果图像');
figure(t+s)
subplot(2,2,2)
imshow(Yimage,[]);
title('模板图像');
Image=(b-a)*Image/(key2(t)-a)-Yimage*(b-key2(t))/(key2(t)-a);
figure(t+s)
subplot(2,2,3)
imshow(Image,[]);
title('恢复图像');
Image=uint8(Image);
figure(t+s)
subplot(2,2,4)
imhist(Image);
title('恢复图像的直方图');
Image=double(Image);
end
PSNRh=PSNR(Image,Ximage);
 RMSEh=RMSE(Image,Ximage);