我正在开发一个处理显微镜图像的管道。每个图像都是tiff平面的堆栈，我想在堆栈上运行反卷积，但我没有关于psf的信息。到目前为止，我一直在使用Matlab函数deconvblind。我深入研究了scikit-image restoration.modules，描述的那个工作正常，但都需要psf。在Matlab中的deconvblind的描述是
- 假设泊松噪声统计，该算法最大化了结果图像与结果PSF卷积时是模糊图像示例的可能性。
有没有一种方法可以在python中实现这个算法，然后估计psf并使用scikit-image提供的其中一个restoration.modules？
1条答案
按热度按时间yv5phkfx1#
从描述中可以清楚地看出，scipy模块实现了Richardson-lucy算法。它很容易实现，也很容易扩展到双盲迭代。例如，Google Biggs在盲解卷积方面的工作。你的数据是否适合双盲反卷积是另一回事。除非数据非常稀疏，否则盲反卷积问题是非常不适定的。
相关原因：别对普通的宽场反褶积期望太高。当涉及到正确地进行反卷积时，魔鬼在细节上非常重要;谁来处理边界条件等等。最好的免费反卷积软件（不幸的是）以imagej插件的形式出现，正如评论中所指出的那样。