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

吾爱光设

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

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

[复制链接]
  • TA的每日心情
    郁闷
    2020-5-8 21:19
  • 签到天数: 8 天

    [LV.3]偶尔看看II

    9

    主题

    28

    回帖

    0

    积分

    小白

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

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

    使用道具 举报

  • TA的每日心情
    擦汗
    2022-1-19 14:17
  • 签到天数: 152 天

    [LV.7]常住居民III

    2

    主题

    17

    回帖

    3

    积分

    小白

    积分
    3
    发表于 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

    9

    主题

    28

    回帖

    0

    积分

    小白

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

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

    使用道具 举报

  • TA的每日心情
    擦汗
    2022-1-19 14:17
  • 签到天数: 152 天

    [LV.7]常住居民III

    2

    主题

    17

    回帖

    3

    积分

    小白

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

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

    使用道具 举报

  • TA的每日心情
    擦汗
    2022-1-19 14:17
  • 签到天数: 152 天

    [LV.7]常住居民III

    2

    主题

    17

    回帖

    3

    积分

    小白

    积分
    3
    发表于 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

    查看全部评分

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

    使用道具 举报

  • TA的每日心情
    开心
    2024-3-28 15:08
  • 签到天数: 809 天

    [LV.10]以坛为家III

    103

    主题

    423

    回帖

    67

    积分

    入门

    积分
    67
    发表于 2022-6-25 13:46 | 显示全部楼层
    请教下,你这里得到成像系统的点扩散函数PSF,这里的PSF指的是二维的矩阵吗?
    发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
    回复

    使用道具 举报

  • TA的每日心情
    开心
    2024-3-28 15:08
  • 签到天数: 809 天

    [LV.10]以坛为家III

    103

    主题

    423

    回帖

    67

    积分

    入门

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

    请问这里的Ny指的是什么?
    发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
    回复

    使用道具 举报

  • TA的每日心情
    擦汗
    2022-1-19 14:17
  • 签到天数: 152 天

    [LV.7]常住居民III

    2

    主题

    17

    回帖

    3

    积分

    小白

    积分
    3
    发表于 2022-7-11 17:23 | 显示全部楼层
    十一月的萧邦 发表于 2022-6-29 08:53
    请问这里的Ny指的是什么?

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

    使用道具 举报

  • TA的每日心情
    开心
    2024-3-28 15:08
  • 签到天数: 809 天

    [LV.10]以坛为家III

    103

    主题

    423

    回帖

    67

    积分

    入门

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

    这个”奈奎斯特频率“值应该是多少呢?我得到的是一个一维矩阵,是归一化的矩阵元素值,那这个矩阵元素值对应的应该是多少奈奎斯特频率?
    发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
    回复

    使用道具 举报

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

    本版积分规则

    关闭

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

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

    GMT+8, 2024-3-28 21:41 , Processed in 0.109375 second(s), 25 queries .

    Powered by Discuz! X3.5

    © 2001-2024 Discuz! Team.

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