最近在研究
MTF 的的底層計算邏輯，將一些心得體會及實現(xiàn)方式示例分享；
}.&;NgZS 2nC,1%kxhq 1、MTF 計算的源頭是 波象差數(shù)據(jù)，通過
光線追跡都可以很方便的得到， 這里我們就用 zemax 的原始 Wavefront map data, sampling 取 512 * 512；
Jy% ? "wn 2、MTF 計算是在 PSF 的基礎(chǔ)上進行傅里葉變換，得到頻域的數(shù)據(jù)，所以先需要對 Wavefront 進行 傅里葉變換得到 PSF; 這里使用
Matlab 進行計算處理，同時附上 Zemax 的計算對比圖：
E/a2b(,Tg e2N K7 OPD_W= load(
'WF_data.TXT');
J ffaT_"\ 0QW=2rs OPD_Cop = exp(1i*2*pi*OPD_W);
j}",+ H v psfW = fftshift(fft2(OPD_Cop));
;m#4Q6k)V? 9$z$yGj l 
7+hF1eoI TUC)S&bC 
3、得到如上的 PSF Data 后，使用變換將 PSF 轉(zhuǎn)到 頻域，代碼實現(xiàn)如下：
@-Q l6k <>9!oOa MTF = abs(fftshift(fft2(PsfData)));
r>F wJm! MTFA = MTF./max(MTF,[],
'all');
% normalize >8k Xa.)84 &=jPt%7#M l#1#3F 4、對 MTF Map 從中心點開始，取水平和垂直兩個方向的的數(shù)據(jù)，使用二維
曲線表示，就是我們經(jīng)常看到的 MTF curve
>~rlnRX MTFT = interp1(frq,MTFA((Nn+1)/2,:),plotfrq,
'cubic');MTFS = interp1(frq,MTFA(:,(Nn+1)/2),plotfrq,
'cubic');
kt)Et `3i>e <m~ w/ rQ OHV{ 對比數(shù)據(jù)看，還是比較吻合的，以上示例大家參考；
