amhalftone2.m

关于半色调编码计算全息图的信息隐藏Matlab程序的说明： 算法实现主程序为CGHwatermark； 位相恢复函数为phaseSL； 用计算全息记录已知振幅和位相的复波前函数为Hol

function [Binary]=AMhalftone(AM)

[M,N]=size(AM);

% phase=PH;
Grayholo=AM;
% 
% if PH==null;
%     x=0
% else
%     x=1;
% end


Binaryholo=double(zeros(2*M,2*N));
for m=1:1:M
    for n=1:1:N
        temp=uint16(Grayholo(m,n)/128);
        
        %phaseoffset=(phase(m,n)+pi)/(pi/4);
        %phaseoffset=uint16(phaseoffset);
        
        cell=ones(temp,1);
      Binaryholo((m-1)*2+1:(m-1)*2+temp,(n-1)*2+2:(n-1)*2+2)=cell;
        

    end
end

 Binary=mat2gray(Binaryholo);
% imshow(Binary);