clear all;

Ximage=double(imread('D:\work\Lena.jpg'));%待隱藏圖像
Yimage=double(imread('D:\work\Monkey.jpg'));%公開圖像
% Ih=size(Ximage,1);%圖像的長度
% Iw=size(Ximage,2);%圖像的寬度
Image=double(zeros(size(Ximage)));
Image=Ximage;
a=2;b=7;
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(i=100:125)%中心切除
     for(j=100:165)
         Image(i,j)=256;
     end
 end
% 增加污點
% wdh=[5,60,127,148,222,175,100,18,78,36,200,198,108];
% wdsh=[200,160,147,55,166,177,100,43,58,25,16,105,236];
% ss=size(wdh,2);
% for(ii=1:ss)
%   for(i=wdh(ii)+1:wdh(ii)+5)
%       for(j=wdsh(ii)+1:wdsh(ii)+5)
%           Image(i,j)=0;
%       end
%   end
% end
  figure(100)
 imshow(Image,[])
%                                                        恢                            復

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);