Restoring Poissonian Images by a Combined First-Order and Second-Order Variation Approach

Algorithm 1

The time-marching gradient descent algorithm for solving (13).

(1) Input 𝑔, 𝐾, 𝑏, MaxIter, 𝛼 and 𝜏. Initialize 𝑢⁰

(2) For k = 1: MaxIter

 compute 𝑢^𝑘 using the explicit scheme (18).

 update the parameter 𝜃 according to (19).

end for