The sparse estimation methods that utilize the l_p-norm, with p being between 0 and 1, have shown better utility in providing optimal solutions to the inverse problem in diffuse optical tomography. These l_p-norm-based regularizations make the optimization function nonconvex, and algorithms that implement l_p-norm minimization utilize approximations to the original l_p-norm function. In this work, three such typical methods for implementing the l_p-norm were considered, namely, iteratively reweighted l_1-minimization (IRL1), iteratively reweighted least squares (IRLS), and the iteratively thresholding method (ITM). These methods were deployed for performing diffuse optical tomographic image reconstruction, and a systematic comparison with the help of three numerical and gelatin phantom cases was executed. The results indicate that these three methods in the implementation of l_p-minimization yields similar results, with IRL1 fairing marginally in cases considered here in terms of shape recovery and quantitative accuracy of the reconstructed diffuse optical tomographic images.
展开▼
