Abstract
Image restoration problems are often solved by finding the minimizer of a suitable objective function. Usually this function consists of a data-fitting term and a regularization term. For the least squares solution, both the data-fitting and the regularization terms are in the ℓ1 norm. In this paper, we consider the least absolute deviation (LAD) solution and the least mixed norm (LMN) solution. For the LAD solution, both the data-fitting and the regularization terms are in the ℓ1 norm. For the LMN solution, the regularization term is in the ℓ1 norm but the data-fitting term is in the ℓ2 norm. Since images often have nonnegative intensity values, the proposed algorithms provide the option of taking into account the nonnegativity constraint. The LMN and LAD solutions are formulated as the solution to a linear or quadratic programming problem which is solved by interior point methods. At each iteration of the interior point method, a structured linear system must be solved. The preconditioned conjugate gradient method with factorized sparse inverse preconditioners is employed to solve such structured inner systems. Experimental results are used to demonstrate the effectiveness of our approach. We also show the quality of the restored images, using the minimization of mixed ℓ1-ℓ2 and ℓ1-ℓ2 norms, is better than that using only the ℓ2 norm.
Original language | English |
---|---|
Pages (from-to) | 1881-1902 |
Number of pages | 22 |
Journal | SIAM Journal on Scientific Computing |
Volume | 27 |
Issue number | 6 |
DOIs | |
Publication status | Published - 25 Jul 2006 |
Scopus Subject Areas
- Computational Mathematics
- Applied Mathematics
User-Defined Keywords
- Interior point method
- Least absolute deviation
- Least squares
- Linear programming
- Preconditioner
- Quadratic programming