test3.m

Go to the documentation of this file.

 
startImg = round(transpose(Phase_Image(Gerchberg_Saxton())).*255)
imshow(startImg./255)
starti=startImg;
        startImg = reshape(startImg,Hologram.pixelX*Hologram.pixelY,1)
        startImg = dec2bin(startImg)
        startImg = permute(startImg, [2,1])
        startImg = reshape(startImg,Hologram.pixelX*Hologram.pixelY*Hologram.ga_imageDepth,1)
        startbit = logical(str2num(startImg))
        startbit = reshape(startbit,Hologram.ga_imageDepth,Hologram.pixelX,Hologram.pixelY)
        startbit = permute(startbit, [2,3,1])
        %startbit = reshape(startbit,Hologram.pixelX,Hologram.pixelY,Hologram.ga_imageDepth)
        %startbit = permute(startbit, [2,3,1]);
 
 
 
img=zeros(Hologram.pixelX,Hologram.pixelY)
        for q=1:size(startbit,3)
            img(:,:) = single(img(:,:) + 2^(8-q)*reshape(startbit(:,:,q),[size(startbit,1),size(startbit,2)]))
        end
 
        figure()
        imshow(img./255)
 
        figure()
        imshow((starti-img)./255)