Image reconstruction of dense scattering media from cw sources using constrained CGD and a matrix rescaling technique

This study reports on results of our efforts to improve the efficiency and stability of previously developed reconstruction algorithms for optical diffusion tomography. The previous studies, which applied regularization and a priori contraints to iterative methods--POCS, CGD, and SART algorithms--showed that in most cases, good quality reconstructions of simply structured media were achievalbe using a perturbation model. The CGD method, which is the most efficient of the three algorithms, was, however, in some instances not able to produce good quality images because of the difficulty in applying range constraints, which can cause divergence. In this study, a scheme is proposed to detect this gradient vector is reset and the CGD reconstruction is restarted using the previous reconstruction as the initial value. In gradient vector is reset and the CGD reconstruction is restarted using the previous reconstruction as the initial value. In addition, a rescaling technique, which rescaled the weight matrix to make it more uniform and less ill-conditioned, is also used to suppress numerical errors and accelerate convergence. Two criteria, rescaling the maximum of each column to 1 and rescaling the average of each column to 1, were applied and compared to results without rescaling. The results show that, with properly imposed constraints, good quality images can be obtained using the CGD method. The convergence speed is much slower when constraints are imposed, but still comparable to the POCS and SART algorithms, The rescaling technique produces more stable and more accurate reconstructions, and speeds up the reconstruction significantly for all three algorithms.