“…This is the basis of image reconstruction algorithms in which an image of Dl a (i.e., the vector x) can be obtained from a set of measurements (i.e., the vector y) through inversion of the matrix A (i.e., the effective pathlengths L i,j ). While several advanced imaging algorithms have been developed-including analytic diffraction tomography approaches (Cheng and Boas, 1998;Li et al, 1997;Matson and Liu, 1999;Schotland, 1997), perturbation approaches (Arridge and Schweiger, 1995;Barbour et al, 1995;O'Leary et al, 1995;Schotland et al, 1993;Yao et al, 1997), the Taylor series expansion approach (Jiang et al, 1996;Paulsen and Jiang, 1995), gradient-based iterative techniques (Arridge and Schweiger, 1998), elliptic systems method (ESM) (Gryazin et al, 1999;Klibanov et al, 1997), and Bayesian conditioning (Barnett et al, 2003;Eppstein et al, 1999) -the most widely used methods for diffuse optical functional brain imaging incorporate a semiinfinite forward model (Kienle and Patterson, 1997a;Patterson et al, 1989) and either backprojection (Colak et al, 1997;Franceschini et al, 2000;Maki et al, 1995, Walker et al, 1997 or perturbation approaches (Arridge, 1999).…”