Just a quick note for those working in Python who might need to simulate PSFs: I made python bindings for Francois Aguet’s C++ PSFmodels. (he had previously made MATLAB mex bindings, but not python). It includes a vectorial model in addition to a scalar Gibson-Lanni model. Even implemented in C++, it’s not nearly as fast as Jizhou Li’s approach as implemented in MicroscPSF-Py, but there are some interesting differences between the vectorial and scalar approximations for high NA and spherical aberration with depth from the coverslip. Francois’ thesis has a lot of detail on the differences between models.
python wheels available for mac/windows on pypi:
pip install psfmodels.
source is on github: