clear all;

Ximage=double(imread('D:\work\Lena.jpg'));%待隱藏圖像
Yimage=double(imread('D:\work\cameraman.jpg'));%公開圖像
Zimage=double(imread('D:\work\Monkey.jpg'));%公開圖像
% Ih=size(Ximage,1);%圖像的長度
% Iw=size(Ximage,2);%圖像的寬度
Image=double(zeros(size(Ximage)));
% Image=Ximage;
a=2;b=150;
t=1;
%                             隱                            藏
% key1=4;
key2=[3,4,5];
s=size(key2,2);

%  for(t=1:s)
figure(1)
% subplot(2,2,1)
imshow(Ximage,[]);
title('待隱藏圖像');
figure(2)
% subplot(2,2,2)
imshow(Yimage,[]);
title('模板圖像');
figure(3)
% subplot(2,2,2)
imshow(Zimage,[]);
title('模板圖像');
Image=(key2(1)-a)*Ximage/(b-a)+Yimage*(b-key2(1))/(b-a);
% figure(4)
% % subplot(2,2,3)
% imshow(Image,[]);
% title('結果圖像');
Image=(key2(2)-a)*Image/(b-a)+Zimage*(b-key2(2))/(b-a);
figure(4)
% 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,Zimage);
 RMSEy=RMSE(Image,Zimage);
% %                          恢                            復
% for(t=1:s)
% figure(6)
% % subplot(2,2,1)
% imshow(Image,[]);
% title('結果圖像');
% figure(7)
% % subplot(2,2,2)
% imshow(Yimage,[]);
% title('模板圖像');
Image=(b-a)*Image/(key2(2)-a)-Zimage*(b-key2(2))/(key2(2)-a);
figure(5)
% subplot(2,2,3)
imshow(Image,[]);
title('恢復圖像');
Image=(b-a)*Image/(key2(1)-a)-Yimage*(b-key2(1))/(key2(1)-a);
figure(6)
% 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);