quanxitu.rar_MATLAB 数字全息_全息水印_数字全息_数字全息水印_计算全息

%基于计算全息图的图像隐藏 clear;clf; [f,map]=imread('cameraman.tif');%读入宿主图像 g=imread('testpat1.png');%读入隐藏信息图像，同时作为谱面约束 p=zeros(256); p(1:end)=(1:65536)/65536*2*pi;%初始相位 % %利用相位恢复算法求输入面的相位 f1=double(f).*exp(i*(p)); F1=fft2(f1); mse=erms(abs(F1),double(g));%计算输出振幅与谱面约束的均方根误差 %while mse>=0.05 for num=1:100 aF1=double(g);%重构傅立叶变换的幅值，用谱面约束代替原来的幅值 pF1=angle(F1); G1=aF1.*exp(i*pF1);%重构后的谱面函数 g1=ifft2(G1);%傅立叶逆变换 ag1=double(f);%重构物面函数的幅值，用物面约束来代替 pg1=angle(g1);%相位保持不变 f1=ag1.*exp(i*pg1);%重构后的物面函数 F1=fft2(f1);%对物面函数进行傅立叶变换 mse=erms(abs(F1),double(g));%求均方根误差 pend=pg1;%获得满足条件的相位函数 end figure(1) %subplot(2,3,1) imshow(pend,map) % %计算全息图的生成 tao=0.1;arf=0.5; for x=1:256 for y=1:256 %x，y表示全息图中像素的位置 t(x,y)=double(f(x,y))*((1-tao)+tao*cos(2*pi*arf*x-pend(x,y))); end end figure(2) %subplot(2,3,2) imshow(uint8(t),map); % %半色调编码--振幅调制 dx=8;dy=8;w=4; [m0,n0]=size(t); h=zeros(m0*dy,n0*dx);%定义编码后图像的大小 for x=1:256 for y=1:256 %x，y表示计算全息图中像素的位置 txy=fix(t(x,y)/max(max(t))*dy);%将计算全息图中的灰度数量化 recth=rect(dy,dx,w,txy);%调用矩形函数，生成一个矩形像元 sposx=(x-1)*dx+1;%编码后图像中行的始点 eposx=x*dx;%编码后图像中行的终点 sposy=(y-1)*dy+1;%编码后图像中列的始点 eposy=y*dy;%编码后图像中列的终点 h(sposx:eposx,sposy:eposy)=recth;%给编码后图像赋值,生成半色调编码图像 end end figure(3) %subplot(2,3,3) imshow(uint8(h),map),title('最终结果'); H=fft2(h); figure(4) %subplot(2,3,4) imshow(uint8(abs(H)),map) figure(5) %subplot(2,3,5) imshow(uint8(abs(H(100:512,1:512))),map)