点扩散函数计算

htsama

最近在研究点扩散函数的计算，遇到点问题，点扩散函数定义上是对广义光瞳函数做傅里叶变换后求模。我找到的公式它长这样


问题是，指数函数上的这个波数到底需不要乘。按这个公式乘波数之后生成的PSF图非常奇怪。网上有些代码里是不乘波数的，我试了可以得到比较正常的图，但也只是类似，并不能完全对上。有懂得大佬么？请指点一二。

不带波数计算出来的,对应zernike项的#9，#14还是挺像的，但是图像上的纹路和文献里的对不上。带波数算的就完全不对劲了。

星辰与风

指数部分是虚数，求模不就没了吗

莱卡莱卡

光瞳函数怎么建立的，用matlab计算的吗

htsama

莱卡莱卡 发表于 2025-1-23 14:35
光瞳函数怎么建立的，用matlab计算的吗

光瞳函数就是公式里傅里叶变换括号里的那个式子，光瞳形状乘以相位分布

QWERT

波数是需要的，能否将对应代码贴出，是否参数设置错误

htsama

QWERT 发表于 2025-1-24 10:51
波数是需要的，能否将对应代码贴出，是否参数设置错误

代码你看下

QWERT

htsama 发表于 2025-1-24 13:40
代码你看下

初步看是空间频率设置有问题，    pupilFunc = (pupilSize*np.exp(-1j*k*zWfAber))这一行把负号去掉 ，再试一下，pupilFunc = (pupilSize*np.exp(j*k*zWfAber))

htsama

QWERT 发表于 2025-1-24 16:15
初步看是空间频率设置有问题，    pupilFunc = (pupilSize*np.exp(-1j*k*zWfAber))这一行把负号去掉 ，再 ...

一样的，这个负号只改变生成图像的朝向。而且公式里这里也是带负号的呀？为什么说影响空间频率呢？

htsama

QWERT 发表于 2025-1-24 16:15
初步看是空间频率设置有问题，    pupilFunc = (pupilSize*np.exp(-1j*k*zWfAber))这一行把负号去掉 ，再 ...

我似乎找到了错误的原因，代码里计算光瞳尺寸的时候是转换到了波长单位的，已经除过一个波长了。算波数的时候又除了一遍，感觉是重复了。我网上找到几个代码，光瞳尺寸都是换算到波长单位的，这有什么意义吗？

baojianing

htsama 发表于 2025-1-26 05:01
我似乎找到了错误的原因，代码里计算光瞳尺寸的时候是转换到了波长单位的，已经除过一个波长了。算波数的 ...

我按照类似图片的公式编写也会得到类似第三张图片的样子，非常模糊。想问问您这是某个单一视场选了几项zernike系数构建的吗？