我正在尝试实现本文介绍的盲图像反卷积算法http://research.microsoft.com/en-us/um/people/luyuan/paper/dualdeblur_cvpr08.pdf
但是我很难理解算法的某些部分:
B1 * k2是矩阵，对吧？那么如何在矩阵上计算(4)？
任何帮助表示赞赏!谢谢!
最佳答案
问题1:E_d(k_1，k_2)是每个空间位置上rho值的总和(在连续设置中为整数)。因此，最小化的功能是图像每个像素的rho值之和。
对于等式(9)，我给你建议一个小把戏。正交于图像梯度方向的各向异性扩散(此处为nablak)可以实现为中值滤波。因此，我建议您通过几次中值滤波来替换各向异性扩散求解器(即，可能应用10次)。这样可以避免您为此处涉及的数学所困扰。