作者seansho (Sean)
看板MATLAB
标题[问题] 关於调制转换函数(MTF)的制作...
时间Mon Aug 5 14:09:33 2013
大家好，小弟有一个问题想请教
关於调制转换函数(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
1F：→ qeagle:看不太懂你问什麽,但 试着回答你问(1)钟形应该是fftshift 08/07 21:22
2F：→ qeagle:的关系,因为FFT左右对称,经过SHIFT之後变成钟形(2)X轴是频 08/07 21:23
3F：→ qeagle:率(取样频率)的NORMALIZE,因为FFT对称,所以只取到一半0.5 08/07 21:24