最近在研究
MTF 的的底層計(jì)算邏輯，將一些心得體會(huì)及實(shí)現(xiàn)方式示例分享；
U for> g{2~G6%;0 1、MTF 計(jì)算的源頭是 波象差數(shù)據(jù)，通過
光線追跡都可以很方便的得到， 這里我們就用 zemax 的原始 Wavefront map data, sampling 取 512 * 512；
u1. 0-Y? 2、MTF 計(jì)算是在 PSF 的基礎(chǔ)上進(jìn)行傅里葉變換，得到頻域的數(shù)據(jù)，所以先需要對(duì) Wavefront 進(jìn)行 傅里葉變換得到 PSF; 這里使用
Matlab 進(jìn)行計(jì)算處理，同時(shí)附上 Zemax 的計(jì)算對(duì)比圖：
: wUi&xw cXYE !( OPD_W= load(
'WF_data.TXT');
w2lO[o~x} l2Rnyb<;; OPD_Cop = exp(1i*2*pi*OPD_W);
x>T+k8[n psfW = fftshift(fft2(OPD_Cop));
z+zEH9.' ]4K4Nh~ 
xjrL@LO# Q=#FvsF#z3 
3、得到如上的 PSF Data 后，使用變換將 PSF 轉(zhuǎn)到 頻域，代碼實(shí)現(xiàn)如下：
Tv~Ho&LS dqFp"Xe"% MTF = abs(fftshift(fft2(PsfData)));
)gAqWbkB MTFA = MTF./max(MTF,[],
'all');
% normalize \,lIPA/L K\mFb LIyb+rH#yg 4、對(duì) MTF Map 從中心點(diǎn)開始，取水平和垂直兩個(gè)方向的的數(shù)據(jù)，使用二維
曲線表示，就是我們經(jīng)�？吹降� MTF curve
k_gl$`A MTFT = interp1(frq,MTFA((Nn+1)/2,:),plotfrq,
'cubic');MTFS = interp1(frq,MTFA(:,(Nn+1)/2),plotfrq,
'cubic');
<B @z>V M<Dvh y[ um2} XI 對(duì)比數(shù)據(jù)看，還是比較吻合的，以上示例大家參考；
