[問題] 關於調制轉換函數(MTF)的製作...

看板MATLAB作者 (Sean)時間12年前 (2013/08/05 14:09), 編輯推噓0(003)
留言3則, 1人參與, 最新討論串1/1
大家好,小弟有一個問題想請教 關於調制轉換函數(MTF)是常用於評估影像的解析力好壞指標 而MTF圖橫軸為空間頻率,縱軸為MTF值 我拍攝了一組黑白邊緣(左邊黑右邊白)的影像,想從這樣的邊緣去計算 為簡化說明,我設定影像為[0 0 0 60 150 255 255 255] 這樣的1*8矩陣 計算方式是由灰階值進行梯度的計算得到一個線散函數LSF 再以LSF進行傅立葉轉換得到MTF 程式碼如下: e=imread('A.jpg'); col=length(e(1,:)); esf=double(e); lsf=gradient(esf); f=fft(lsf); mtf=abs(f); plot(fftshift(mtf)); 到這邊結果呈現是一個鐘型的圖形 而MTF圖理論上是一個像這樣的圖形http://www.quickmtf.com/images/mtf_samples.png
到這邊有幾個疑惑: 1.為何圖形形狀不同? 2.x軸的單位是甚麼? 寫了一長串十分抱歉,這個問題困擾我很久了,希望大家可以幫我解答,謝謝 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.113.134.209

08/07 21:22, , 1F
看不太懂你問什麼,但 試著回答你問(1)鐘形應該是fftshift
08/07 21:22, 1F

08/07 21:23, , 2F
的關係,因為FFT左右對稱,經過SHIFT之後變成鐘形(2)X軸是頻
08/07 21:23, 2F

08/07 21:24, , 3F
率(取樣頻率)的NORMALIZE,因為FFT對稱,所以只取到一半0.5
08/07 21:24, 3F
文章代碼(AID): #1H_q6V-W (MATLAB)
文章代碼(AID): #1H_q6V-W (MATLAB)