24小时热门版块排行榜    

查看: 976  |  回复: 0
【悬赏金币】回答本帖问题,作者木心的云雀将赠送您 10 个金币

木心的云雀

新虫 (初入文坛)

[求助] 傅里叶光学中广义光瞳函数关于采样点的一点疑问

已知参数 像面640*640个像素,单个像素大小15μm;光瞳大小60mm,焦距100mm 波长1.2μm
我的心路过程:
广义光瞳函数函数做傅里叶变换得到脉冲相应,如果是非相干系统则取模平方,相干系统则不用,再做傅里叶变换得传递函数。
为了方便matlab编程,希望光学传递函数是一个和图像一样大小640*640的矩阵,这样的话在频域里就能直接相乘了。
为了能是两者直接相乘,需要保证图像频谱和光学传递函数频谱在频域内的频率间隔相等。
图像的频谱间隔为fx = 1/(640*15*10^-6);
第一次傅里叶变换时,有关系δx =fx *lambda *di 。考虑无穷远成像,di即为焦距。利用60mm/δx就得到了光瞳面上的采样点间隔。
经计算,光瞳面上需采样4800*4800个采样点,此时频率间隔相等。

疑问:
虽然频率域的频率间隔相等了,但是图像频谱矩阵和传递函数矩阵大小不同,又不能直接相乘了。
请问我的采样点计算方法是否有问题吗?如果没有问题,那么后续应该怎么处理呢?
我最初认为,应当在第二次傅里叶变换后(即由脉冲相应求傅里叶变换得到传递函数后),对传递函数矩阵进行裁剪,即外围高频信息无用。
但是转念又一想,连做两次傅里叶变换则相当于还是广义光瞳函数本身,那么就相当于对广义光瞳函数本身进行裁剪,那么4800*4800的矩阵若仅保留中间640*640矩阵的话,岂不是整个广义光瞳函数外部大量数据都是无用的?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 木心的云雀 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见