Introduction - If you have any usage issues, please Google them yourself
This is the code that was used in the papers "A Nonnnegatively Constrained Convex Programming Method for Image Reconstruction", "Total Variation-Penalized Poisson Likelihood Estimation for Ill-Posed Problems", "Tikhonov Regularized Poisson Likelihood Estimation: Theoretical Justification and a Computational Method", "An Efficient Computational Method for Total Variation with Poisson Negative-Log Likelihood", "An Analysis of Regularization by Diffusion for Ill-Posed Poisson Likelihood Estimation," "An Iterative Method for Edge-Preserving MAP Estimation when Data-Noise is Poisson", and finally, "Regularization Parameter Selection Methods for Ill-Posed Poisson Maximum Likelihood Estimation". See my publications page for more details. The main algorithm is for nonnegatively constrained, regularized Poisson likelihood estimation. At this point you can choose Tikhonov, total variation regularization, and diffusion regularization. A number of other methods are also implemented. Regularizatio