me.m

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

clear;close all;clc;
T=double(imread('hEaplpha001t009.bmp'));A=T(1:64,1:64);
figure,imshow(mat2gray(log(1+abs(fftshift(fft2(A))))));
[M,N]=size(A);
% phase=PH;
Grayholo=A;
% 
% if PH==null;
%     x=0
% else
%     x=1;
% end

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

    end
end

Binary=mat2gray(Binaryholo);
figure,imshow(Binary);