作者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