作者a6902835 (无赖正义)
看板MATLAB
标题[讨论] 全像片还原
时间Fri Apr 13 10:44:29 2012
最近在研究全像片还原的数学模型
看到paper上有列举一些公式
想说跟着做做看
却一直 error
附上code 麻烦高手解答
clc
clear all;
clc
%==========================================================================
h=imread('C:\Users\Batista\Desktop\Figure\dice.bmp');
figure(1):imshow(h)
%==============Parameter===================================================
lamda=632.8*10^-9;
dis=1.054;
N=1024;
m=0:1023;
n=0:1023;
del_x=6.8*10^-6;
del_y=6.8*10^-6;
R=200;
%==========================================================================
g1=(lamda.^2*(n+((N.^2*del_x.^2)/(2*dis*lamda))).^2)/N.^2*del_x.^2;
g2=(lamda.^2*(m+((N.^2*del_y.^2)/(2*dis*lamda))).^2)/N.^2*del_y.^2;
G=exp(((-j*2*pi*dis)/lamda)*sqrt(1-g1-g2));
F=fft(h.*R);
Re=ifft(F.*G);
figure(2):imshow(Re)
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.116.155.66
1F：→ ppt1:倒数第二行，F跟G矩阵不一致.. 04/13 12:23
2F：→ ppt1:漏打size 04/13 12:24
3F：→ a6902835:漏打 size 是 ?? 04/18 15:16