查看: 11296|回复: 8

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

[复制链接]

小胡同志

TA的每日心情

郁闷
2020-5-8 21:19

签到天数: 8 天

[LV.3]偶尔看看II

主题

回帖

积分

积分: 0

收听TA

发消息

发表于 2020-7-13 11:50 | 显示全部楼层 |阅读模式

想请教下大家，我得到了成像系统的点扩散函数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

TA的每日心情

郁闷
2024-9-29 16:36

签到天数: 153 天

[LV.7]常住居民III

主题

回帖

积分

积分: 3

收听TA

发消息

发表于 2020-7-23 19:58 | 显示全部楼层

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

评分

参与人数 1金币 +2 收起理由

天空很蓝

+ 2

查看全部评分

发帖求助前要善用【论坛搜索】功能，那里可能会有你要找的答案；

使用道具举报

小胡同志

TA的每日心情

郁闷
2020-5-8 21:19

签到天数: 8 天

[LV.3]偶尔看看II

主题

回帖

积分

积分: 0

收听TA

发消息

楼主| 发表于 2020-7-28 12:03 | 显示全部楼层

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

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

发帖求助前要善用【论坛搜索】功能，那里可能会有你要找的答案；

使用道具举报

chenhuanliufeng

TA的每日心情

郁闷
2024-9-29 16:36

签到天数: 153 天

[LV.7]常住居民III

主题

回帖

积分

积分: 3

收听TA

发消息

发表于 2020-7-29 08:45 | 显示全部楼层

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

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

发帖求助前要善用【论坛搜索】功能，那里可能会有你要找的答案；

使用道具举报

chenhuanliufeng

TA的每日心情

郁闷
2024-9-29 16:36

签到天数: 153 天

[LV.7]常住居民III

主题

回帖

积分

积分: 3

收听TA

发消息

发表于 2020-7-29 08:49 | 显示全部楼层本楼为最佳答案

这个最佳答案由 chenhuanliufeng 给出，感谢 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）

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

评分

参与人数 1贡献值 +1 收起理由

天空很蓝

+ 1

查看全部评分

发帖求助前要善用【论坛搜索】功能，那里可能会有你要找的答案；

使用道具举报

十一月的萧邦

TA的每日心情

开心
2024-11-4 10:18

签到天数: 843 天

[LV.10]以坛为家III

主题

回帖

积分

积分: 67

收听TA

发消息

发表于 2022-6-25 13:46 | 显示全部楼层

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

发帖求助前要善用【论坛搜索】功能，那里可能会有你要找的答案；

使用道具举报

十一月的萧邦

TA的每日心情

开心
2024-11-4 10:18

签到天数: 843 天

[LV.10]以坛为家III

主题

回帖

积分

积分: 67

收听TA

发消息

发表于 2022-6-29 08:53 | 显示全部楼层

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

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

发帖求助前要善用【论坛搜索】功能，那里可能会有你要找的答案；

使用道具举报

chenhuanliufeng

TA的每日心情

郁闷
2024-9-29 16:36

签到天数: 153 天

[LV.7]常住居民III

主题

回帖

积分

积分: 3

收听TA

发消息

发表于 2022-7-11 17:23 | 显示全部楼层

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

奈奎斯特频率

发帖求助前要善用【论坛搜索】功能，那里可能会有你要找的答案；

使用道具举报

十一月的萧邦

TA的每日心情

开心
2024-11-4 10:18

签到天数: 843 天

[LV.10]以坛为家III

主题

回帖

积分

积分: 67

收听TA

发消息

发表于 2022-7-12 09:31 | 显示全部楼层

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

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

发帖求助前要善用【论坛搜索】功能，那里可能会有你要找的答案；

使用道具举报

返回列表发新帖