[已解决]通过PSF计算几何MTF的问题？

小胡同志

想请教下大家，我得到了成像系统的点扩散函数PSF之后，通过:OTF=fft2(PSF);  
MTF=abs(OTF)/max(OTF);

计算出MTF之后，然后取了X方向的值，现在的问题是不知道这个X方向的MTF对应的横坐标值应该是多少，横坐标单位为lp/mm；感谢大家的回复！

最佳答案

月排行榜 / 总排行榜

chenhuanliufeng

2020-7-29 08:49

小胡同志 发表于 2020-7-28 12:03
感谢回复！感觉你说的很有意思，请问下这个是将什么的零频分量移动到中心啦，不是很明白，期待你的回复。 ...

[size=+3]fftshift

将离散傅立叶变换的零频分量移到频谱中心

[size=+1]句法

• Y = fftshift（X）Y = fftshift（X，dim）
  

[size=+1]描述

Y = fftshift(X) 重新排列的输出fft，fft2和fftn通过零的频率成分移动到阵列的中心。这对于可视化频谱中间具有零频率分量的傅立叶变换很有用。

对于向量，fftshift(X)交换的左右两半X。对于矩阵，fftshift(X)将第一象限与第三象限交换，将第二象限与第四象限交换。

对于高维数组，fftshift(X)交换X沿每个维的“半空间” 。

Y = fftshift(X,dim) fftshift沿维应用操作dim。

[size=+1]例子

对于任何矩阵 X

• Y = fft2（X）
  

有Y(1,1) = sum(sum(X)); 信号的零频率分量位于二维FFT的左上角。对于

• Z = fftshift（Y）
  

这个零频率分量在矩阵的中心附近。

跳转到最佳答案楼层

chenhuanliufeng

先用fftshift将零频分量移动到中心；2D的MTF取X方向的的slice后，横坐标就是-Ny~+Ny

小胡同志

chenhuanliufeng 发表于 2020-7-23 19:58
先用fftshift将零频分量移动到中心；2D的MTF取X方向的的slice后，横坐标就是-Ny~+Ny

感谢回复！感觉你说的很有意思，请问下这个是将什么的零频分量移动到中心啦，不是很明白，期待你的回复。

chenhuanliufeng

小胡同志 发表于 2020-7-28 12:03
感谢回复！感觉你说的很有意思，请问下这个是将什么的零频分量移动到中心啦，不是很明白，期待你的回复。 ...

假设你的PSF中心是极大点。OTF=abs(fft2(PSF))，得到的二维图中心为0，四个角是峰值。fftshift就是对OTF的数据数据进行了平移，以使得峰值在图中央。

chenhuanliufeng

小胡同志 发表于 2020-7-28 12:03
感谢回复！感觉你说的很有意思，请问下这个是将什么的零频分量移动到中心啦，不是很明白，期待你的回复。 ...

[size=+3]fftshift

将离散傅立叶变换的零频分量移到频谱中心

[size=+1]句法

• Y = fftshift（X）Y = fftshift（X，dim）
  

[size=+1]描述

Y = fftshift(X) 重新排列的输出fft，fft2和fftn通过零的频率成分移动到阵列的中心。这对于可视化频谱中间具有零频率分量的傅立叶变换很有用。

对于向量，fftshift(X)交换的左右两半X。对于矩阵，fftshift(X)将第一象限与第三象限交换，将第二象限与第四象限交换。

对于高维数组，fftshift(X)交换X沿每个维的“半空间” 。

Y = fftshift(X,dim) fftshift沿维应用操作dim。

[size=+1]例子

对于任何矩阵 X

• Y = fft2（X）
  

有Y(1,1) = sum(sum(X)); 信号的零频率分量位于二维FFT的左上角。对于

• Z = fftshift（Y）
  

这个零频率分量在矩阵的中心附近。

十一月的萧邦

请教下，你这里得到成像系统的点扩散函数PSF，这里的PSF指的是二维的矩阵吗？

十一月的萧邦

chenhuanliufeng 发表于 2020-7-23 19:58
先用fftshift将零频分量移动到中心；2D的MTF取X方向的的slice后，横坐标就是-Ny~+Ny

请问这里的Ny指的是什么？

chenhuanliufeng

十一月的萧邦 发表于 2022-6-29 08:53
请问这里的Ny指的是什么？

奈奎斯特频率

十一月的萧邦

chenhuanliufeng 发表于 2022-7-11 17:23
奈奎斯特频率

这个”奈奎斯特频率“值应该是多少呢？我得到的是一个一维矩阵，是归一化的矩阵元素值，那这个矩阵元素值对应的应该是多少奈奎斯特频率？