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

木心的云雀

已知参数 像面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矩阵的话，岂不是整个广义光瞳函数外部大量数据都是无用的？