zhanggc.m

这么代码是用Matlab编写的

clear all;

Ximage=double(imread('D:\MATLAB7\work\Lena.jpg'));
Yimage=double(imread('D:\MATLAB7\work\Monkey.jpg'));
Zimage=double(imread('D:\MATLAB7\work\cameraman.jpg'));
Ih=size(Ximage,1);%图像的长度
Iw=size(Ximage,2);%图像的宽度
Image=double(zeros(size(Ximage)));
Image=Ximage;
s=5;
renewing=0;
hiding=0;

                                           %                隐            藏               program
%key1=4.0;
key2=0.1;

for(t=1:s)
figure(t)
subplot(4,2,1)
imshow(Image,[]);;%显示图像
title('待隐藏图像');
Image=uint8(Image);
figure(t)
subplot(4,2,2)
imhist(Image);%显示图像数据的柱状图
title('待隐藏图像的直方图');
Image=double(Image);

figure(t)
subplot(4,2,3)
imshow(Yimage,[]);
title('模板图像1');
Yimage=uint8(Yimage);
figure(t)
subplot(4,2,4)
imhist(Yimage);
title('模板图像1的直方图');
Yimage=double(Yimage);

figure(t)
subplot(4,2,5)
imshow(Zimage,[])
title('模板图像2');
Zimage=uint8(Zimage);
figure(t)
subplot(4,2,6)
imhist(Zimage);
title('模板图像2的的直方图');
Zimage=double(Zimage);


for(i=1:Ih)
        for(j=1:Iw)
%            key2=key1*key2*(1-key2)
%             key2=1-key1*(key2*key2);
            if( key2< 0.0001)
%                 key2=key2;
                Image(i,j)=mod((key2+0.01)*(Image(i,j)-Yimage(i,j))+Yimage(i,j),256);
     %           key2=(key1*key2*(1-key2))/2;
                Image(i,j)=mod((key2+0.01)*(Image(i,j)-Zimage(i,j))+Zimage(i,j),256);
            else
               Image(i,j)=mod(key2*(Image(i,j)-Zimage(i,j))+Zimage(i,j),256);
               Image(i,j)=mod(key2*(Image(i,j)-Yimage(i,j))+Yimage(i,j),256);
      %          key2=(key1*key2*(1-key2))/2;
                
            end
        end
end

figure(t)
subplot(4,2,7)
imshow(Image,[]);
title('结果图像');
Image=uint8(Image);
figure(t)
subplot(4,2,8)
imhist(Image);
title('结果图像的直方图');
Image=double(Image);
end
%                隐  藏  的  计  算  均  方  根  误  差
for(i=1:Ih)
        for(j=1:Iw)
            if(Image(i,j)~=Zimage(i,j))
                hiding=hiding+(Image(i,j)-Zimage(i,j))^2;
            end
        end
    end
    
   hiding=hiding/(Ih*Iw);
   hiding=hiding^(1/2);
                                          %               恢             复             program
% % key1=4;
% % key2=0.9;
% 
% for(t=1:s)
% figure(t+s)
% subplot(4,2,1)
% imshow(Image,[]);
% title('结果图像');
% Image=uint8(Image);
% figure(t+s)
% subplot(4,2,2)
% imhist(Image);
% title('结果图像的直方图');
% Image=double(Image);
% 
% figure(t+s)
% subplot(4,2,3)
% imshow(Yimage,[]);
% title('模板图像');
% Yimage=uint8(Yimage);
% figure(t+s)
% subplot(4,2,4)
% imhist(Yimage);
% title('模板图像的直方图');
% Yimage=double(Yimage);
% 
% figure(t+s)
% subplot(4,2,5)
% imshow(Yimage,[]);
% title('模板图像');
% Yimage=uint8(Yimage);
% figure(t+s)
% subplot(4,2,6)
% imhist(Yimage);
% title('模板图像的直方图');
% Yimage=double(Yimage);
% 
% key=0; 
% 
% for(i=1:Ih)
%        for(j=1:Iw)
% %             key2=1-key1*(key2*key2);
% %             key2=(key1*key2*(1-key2))/2;
%             if( key2< 0.001)
% %                 key2=key2;
%                 key=1/(key2+0.1);
%                 Image(i,j)=mod(key*(Image(i,j)-Zimage(i,j))+Zimage(i,j),256);
%                 Image(i,j)=mod(key*(Image(i,j)-Yimage(i,j))+Yimage(i,j),256)
%             else
%                 key=1/key2;
%                 Image(i,j)=mod(key*(Image(i,j)-Zimage(i,j))+Zimage(i,j),256);
%                 Image(i,j)=mod(key*(Image(i,j)-Yimage(i,j))+Yimage(i,j),256)
%             end
% 
%         end
% end
% figure(t+s)
% subplot(4,2,7)
% imshow(Image,[]);
% title('恢复图像');
% Image=uint8(Image);
% figure(t+s)
% subplot(4,2,8)
% imhist(Image);
% title('恢复图像的直方图');
% Image=double(Image);
% end
% 
% %             计  算  恢  复  的  均  方  根  误  差
% for(i=1:Ih)
%         for(j=1:Iw)
%             if(Image(i,j)~=Ximage(i,j))
%                 renewing=renewing+(Image(i,j)-Ximage(i,j))^2;
%             end
%         end
%     end
%     
%    renewing=renewing/(Ih*Iw);
%    renewing=renewing^(1/2);