function?[OTFPSF]?=?Aperture2OTF（im）
%-----------------------------------
%?從孔徑函數轉換為OTF
%?clc;clear;close?all;
%?f=imread（‘aperture512_16.tif‘）;
im=im2double（im）;??%轉換成雙精度型
%-------------------------------------
%?%因為OTF是孔徑函數的歸一化自相關函數，而自相關函數與功率譜互為FT變換
%?ff=fft2（im）;???
%?ff=abs（ff）.^2;?????????%求功率譜
%?%計算OTF
%?fff=fft2（ff）;?????????%求反FT變換
%?
%?ff_max=max（fff（:））;????%歸一化處理
%?OTF=fff/ff_max;
%?%???計算PSF
%?PSF=ifft2（OTF）;
%?PSF=IMRA（PSF）;
%因為OTF是孔徑函數的歸一化自相關函數，而自相關函數與功率譜互為FT變換
ff=fft2（im）;???
ff=abs（ff）.^2;?????????%求功率譜
fff=ifft2（ff）;?????????%求反FT變換
fff=IMRA（fff）;
ff_max=max（fff（:））;????%歸一化處理
OTF=fff/ff_max;
PSF=ifft2（OTF）;
PSF=IMRA（PSF）;
%---------------------------------------------------
%T2=abs（fftshift（OTF））;
%figureimshow（log（1+T2）[]）;
%figuremesh（T

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2170??2011-06-20?10:20??RL.m

?????文件??????10455??2012-04-29?16:16??super_resolution.m

?????文件???????1078??2011-06-20?10:20??Aperture2OTF.m

?????文件???????1845??2011-06-20?10:20??ElipsoidalAperture.m

-----------?---------??----------?-----??----

????????????????15548????????????????????4