amhalftone1.m

来自「关于半色调编码计算全息图的信息隐藏Matlab程序的说明： 算法实现主程」· M 代码 · 共 32 行

function [Binary]=AMhalftone(AM)

[M,N]=size(AM);

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


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

    end
end

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