cghwatermark.asv

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

clear;close all;clc;
Master=double(imread('lenna.bmp'));
Watermark=double(imread('E.bmp'));

phase=phaseSL(Master,Watermark);

imwrite(mat2gray(phase),'phase.bmp');         %save phase picture
% testphase=phase+pi.*rand(256);

Holo=Hologen(Master,phase,0.01,0.09);   %第一参数为宿主，第二为相位函数，第三参考光空间频率，第四嵌入权重
%Holo=fftshift(Holo);
%DAPM=halftone(Master,phase);
%imwrite(DAPM,'dpm.bmp');
imwrite(mat2gray(Holo),'hEaplpha001t009.bmp');
%imshow(mat2gray(Holo));
figure;imshow(mat2gray(log(1+abs(fftshift(fft2(Holo))))));

%Bholo=AMhalftone(Holo);  
Bholo=AMhalftone1(Holo);
figure;imshow(Bholo);
imwrite(Bholo,'optoexp/BhEaplpha001t009.bmp');   //                 
%Bholo(501:1700,501:1700)=zeros(1200);

% figure(1)
% imshow(mat2gray(Bholo));
% figure(2)
%  imwrite(mat2gray(log(1+abs(fftshift(fft2(Bholo))))),'mark.bmp')
% imshow(mat2gray(log(1+abs(fftshift(fft2(Bholo))))));