基于菲涅尔域的双随机相位编码图像加密方案，利用角谱传播法实现matlab.zip

%基于fresnel变换的无透镜4－f加解密模拟，fresnel变换用角谱传播理论计算 clear; clc; np=256;%采样点 lambda=633e-9;%入射波长(key) step=50e-6; % um 空域采样 (dx<N*lambda/2) range=np*step; % x,y的大小范围 scale=-range/2:step:range/2-step; %物平面采样 ftscale=(np/range^2)*scale; %等价于d_ftscale＝1/range(空间频率采样） [x,y]=meshgrid(scale,scale); %物平面坐标 [fx,fy]=meshgrid(ftscale,ftscale);%空间频率域坐标 %mt=t(x,y);%输入函数 mt=double(imread('baboon.bmp'))/255; %明文 mask1=double(imread('RM1.bmp','bmp'))/255; mask2=double(imread('RM2.bmp','bmp'))/255; % First stage rt1=exp(1i*2*pi*mask1); %第一块相位模板 fai=mt.*rt1; dt=fftshift(fft2(fai)); %先变换到频域 distance1=2; ff=h2(fx,fy,distance1,lambda);%角谱传播函数 ft=ff.*dt; %传播后的频谱 % Second stage rt2=exp(1i*2*pi*mask2); %第二块相位模板 ft=ft.*rt2; %频域滤波 distance2=2; ff=h2(fx,fy,distance2,lambda);%角谱传播函数 ft=ft.*ff; %最后频域密文 %显示密文 temp=abs(ifft2(fftshift(ft))); %密文 figure(1), imshow(mat2gray(temp)); %---------------------------------------------------- % 解密过程 %---------------------------------------------------- %回传2 distance1=2;%可以测试灵敏度 （0.1量级） distance2=2; lambda=633e-9; %(1nm量级） ff=h2(fx,fy,-distance2,lambda); ft=ft.*ff.*conj(rt2); %回传1 ff=h2(fx,fy,-distance1,lambda); ft=ft.*ff; %由频域变换到空域 ft=ifft2(ft); ft=ft.*conj(rt1); temp=mat2gray(abs(ft)); %显示解密结果 imwrite(temp,'n8.bmp','bmp'); figure(2), imshow(temp);