我需要在没有 MATLAB 中的点扩散函数的情况下对荧光珠的图像进行反卷积

I need to deconvolve the image of a fluorescent bead without the Point Spread Function in MATLAB

我知道我可以使用 LUCY, Weiner and Blind 等反卷积算法,但我对如何估计 LUCY 和 Weiner 算法所需的 PSF(点扩散函数)感到困惑。我目前正在使用 Widefield Florescence 显微镜设置。

我在互联网上的很多地方都读到过代表亚分辨率物体的荧光珠,因此代表脉冲函数。虽然我知道 OTF(光学传递函数)由我得到的图像表示,但我如何叠加以使用完全相同的图像对我拥有的图像进行反卷积?

这是我尝试去卷积的图像示例:

我有一个程序可以自动选择珠子并裁剪图像的其余部分。我需要做的就是对此应用反卷积。

我的最终目标是通过对从视频中获得的一堆 2D 切片进行反卷积并用 3D 软件表示它们来实现光学切片。

所以有人可以从理论上(或在 MATLAB 方面,这是首选)解释我需要做什么。


一种可能的解决方案

与几个人交谈后,我注意到 EPFL 有一个 Java resource 可以从 MATLAB 调用。它是一个 PSF 发生器,根据您的显微镜规格而定。此生成的 PSF 可与 LUCY 和 Weiner 一起使用。


备选方案

我还找到了 this 资源,这是一个 MATLAB 库,专用于我目前正在从事的工作。


This(如果有人也尝试回答这个问题,我会很高兴)是我将如何衡量反卷积算法的性能。我相信我的反卷积的最终结果将取决于 PSF 估计(EPFL 资源或 Praveen 的算法)和反卷积算法(Blind、LUCY、Weiner 等)。一旦完成,我将 post 所有结果的 table 放在这里。

正如我在问题中所述,关键是创建 PSF 模型,使用一种理论 PSF 生成器,然后输入 3D 图像堆栈和 3D PSF 堆栈。我还回答了有关在我的特定情况下反卷积算法性能测量的链接问题,所以请不要犹豫让我清除任何内容。