Deconvolution with Least Squares

We are given a noisy observation x0 = h ∗ x + noise. We want to recover the true sharp image x by solving the Tikhonov-regularized least-squares problem minx ‖ h ∗ x − x022 + λ ‖ x ‖22. "Run" solves it using the either closed-form (exact) or gradient descent (iterative).

 | 

Target

The sharp image x we want to recover.
Resize: 100%

Observation

x0 = h ∗ x + noise
Blur σ: 0.5px

Deconvolved (Estimated)

x̂ ≈ argminx ‖ h ∗ x − x02 + λ ‖ x ‖2
Regularization λ: 0.020
(click "Run")

Pixel-wise Difference

|x̂ − x| (Target vs. Deconvolved)
mean |xhat - x|: —