请选择 进入手机版 | 继续访问电脑版

光学社区

 找回密码
 注册
会员须知
会员须知
实用帮助
实用帮助
查看: 290|回复: 4

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

[复制链接]
最佳答案
0 

9

主题

35

帖子

0

积分

小白

Rank: 1

积分
0
发表于 2020-7-13 11:50 | 显示全部楼层 |阅读模式
想请教下大家,我得到了成像系统的点扩散函数PSF之后,通过:OTF=fft2(PSF);  
MTF=abs(OTF)/max(OTF);

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


最佳答案
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)
这个零频率分量在矩阵的中心附近。

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

使用道具 举报

最佳答案
1 

2

主题

13

帖子

5

积分

小白

Rank: 1

积分
5
发表于 2020-7-23 19:58 | 显示全部楼层
先用fftshift将零频分量移动到中心;2D的MTF取X方向的的slice后,横坐标就是-Ny~+Ny

评分

参与人数 1金币 +2 收起 理由
天空很蓝 + 2

查看全部评分

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

使用道具 举报

最佳答案
0 

9

主题

35

帖子

0

积分

小白

Rank: 1

积分
0
 楼主| 发表于 2020-7-28 12:03 | 显示全部楼层
chenhuanliufeng 发表于 2020-7-23 19:58
先用fftshift将零频分量移动到中心;2D的MTF取X方向的的slice后,横坐标就是-Ny~+Ny

感谢回复!感觉你说的很有意思,请问下这个是将什么的零频分量移动到中心啦,不是很明白,期待你的回复。
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
回复

使用道具 举报

最佳答案
1 

2

主题

13

帖子

5

积分

小白

Rank: 1

积分
5
发表于 2020-7-29 08:45 | 显示全部楼层
小胡同志 发表于 2020-7-28 12:03
感谢回复!感觉你说的很有意思,请问下这个是将什么的零频分量移动到中心啦,不是很明白,期待你的回复。 ...

假设你的PSF中心是极大点。OTF=abs(fft2(PSF)),得到的二维图中心为0,四个角是峰值。fftshift就是对OTF的数据数据进行了平移,以使得峰值在图中央。
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
回复

使用道具 举报

最佳答案
1 

2

主题

13

帖子

5

积分

小白

Rank: 1

积分
5
发表于 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)
这个零频率分量在矩阵的中心附近。

评分

参与人数 1贡献值 +1 收起 理由
天空很蓝 + 1

查看全部评分

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

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /2 下一条

联系我们|本论坛只支持PC端注册|手机版|小黑屋|光学社区 ( 粤ICP备15067533号 )

GMT+8, 2020-8-7 02:47 , Processed in 0.218750 second(s), 29 queries .

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表