OSA's Digital Library

Biomedical Optics Express

Biomedical Optics Express

  • Editor: Joseph A. Izatt
  • Vol. 3, Iss. 5 — May. 1, 2012
  • pp: 1006–1024
« Show journal navigation

Parametric level set reconstruction methods for hyperspectral diffuse optical tomography

Fridrik Larusson, Sergio Fantini, and Eric L. Miller  »View Author Affiliations


Biomedical Optics Express, Vol. 3, Issue 5, pp. 1006-1024 (2012)
http://dx.doi.org/10.1364/BOE.3.001006


View Full Text Article

Acrobat PDF (1117 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

A parametric level set method (PaLS) is implemented for image reconstruction for hyperspectral diffuse optical tomography (DOT). Chromophore concentrations and diffusion amplitude are recovered using a linearized Born approximation model and employing data from over 100 wavelengths. The images to be recovered are taken to be piecewise constant and a newly introduced, shape-based model is used as the foundation for reconstruction. The PaLS method significantly reduces the number of unknowns relative to more traditional level-set reconstruction methods and has been show to be particularly well suited for ill-posed inverse problems such as the one of interest here. We report on reconstructions for multiple chromophores from simulated and experimental data where the PaLS method provides a more accurate estimation of chromophore concentrations compared to a pixel-based method.

© 2012 OSA

1. Introduction

Near-infrared light has proven to be useful for imaging the human body and providing functional information for applications including breast cancer detection and characterization [1

1. S. Fantini, E. L. Heffer, V. E. Pera, A. Sassaroli, and N. Liu, “Spatial and spectral information in optical mammography,” Technol. Cancer Res. Treat. 4, 471–482 (2005). [PubMed]

5

5. S. Fantini, D. Hueber, M. A. Franceschini, E. Gratton, W. Rosenfeld, P. G. Stubblefield, D. Maulik, and M. R. Stankovic, “Non-invasive optical monitoring of the newborn piglet brain using continuous-wave and frequency-domain spectroscopy,” Phys. Med. Biol. 44, 1543–1563 (1999). [CrossRef] [PubMed]

]. In recent years significant improvements have been made in diffuse optical tomography (DOT) by using data from multiple wavelengths. The use of this kind of hyperspectral information has moved DOT away from the recovery of space and time varying maps of absorption and scattering properties of the breast to the recovery of chromophore concentrations as well as diffusion amplitude. Inverting for multiple chromophores such as hemoglobin, lipids and water has been shown to provide an improved ability to localize tumours or other objects of interest in the breast cancer application [6

6. S. Kukreti, A. E. Cerussi, W. Tanamai, D. Hsiang, B. J. Tromberg, and E. Gratton, “Characterization of metabolic differences between benign and malignant tumors: high-spectral-resolution diffuse optical spectroscopy,” Radiology 254, 277–284 (2010). [CrossRef]

8

8. H. Dehghani, B. W. Pogue, S. P. Poplack, and K. D. Paulsen, “Multiwavelength three-dimensional near-infrared tomography of the breast: initial simulation, phantom, and clinical results,” Appl. Opt. 42, 135–145 (2004). [CrossRef]

]. While these and other advancements have been critical for moving DOT from the lab into the clinic, there remain significant obstacles to be addressed so that the method can be used to stably recover these quantities.

The imaging of chromophore concentration and diffusion amplitude from DOT data requires the solution of an ill-posed and non-linear inverse problem [3

3. D. A. Boas, D. H. Brooks, E. L. Miller, C. A. DiMarzio, M. Kilmer, R. J. Gaudette, and Q. Zhang, “Imaging the body with diffuse optical tomography,” IEEE Signal Process. Mag. 18(6), 57–75 (2001). [CrossRef]

]. The physics associated with the photons travelling through turbid media, such as breast tissue, makes the problem ill-posed and therefore introduces significant challenges for the DOT problem. The ill-posedness make the reconstruction highly sensitive to noise and un-modeled effects which can reduce accuracy in the recovered images [7

7. A. Li, G. Boverman, Y. Zhang, D. Brooks, E. L. Miller, M. E. Kilmer, Q. Zhang, E. M. C. Hillman, and D. A. Boas, “Optimal linear inverse solution with multiple priors in diffuse optical tomography,” Appl. Opt. 44, 1948–1956 (2005). [CrossRef] [PubMed]

]. Additional problems are encountered when recovering multiple chromophores such as crosstalk, e.g. when two spatially disjoint chromophores pollute the images of each other.

The ill-posedness of the DOT problem is traditionally solved by putting the recovery of images in a variational context. To achieve this, numerous regularization techniques [9

9. J. P. Culver, V. Ntziachristos, M. J. Holboke, and A. G. Yodh, “Optimization of optode arrangements for diffuse optical tomography: a singular-value analysis,” Opt. Lett. 26, 701–703 (2001). [CrossRef]

] can be used to reduce artifacts and improve the overall accuracy. Traditionally, methods such as Tikhonov and total variation (TV) regularization have been employed along with L-curve methods to optimally choose regularization parameters [10

10. F. Larusson, S. Fantini, and E. L. Miller, “Hyperspectral image reconstruction for diffuse optical tomography,” Biomed. Opt. Express 2, 947–965 (2011). [CrossRef]

12

12. K. D. Paulsen and H. Jiang, “Enhanced frequency-domain optical image reconstruction in tissues through total-variation minimization,” Appl. Opt. 35, 3447–3458, (1996). [CrossRef] [PubMed]

]. Increasing the amount of data used to solve the DOT problem has also proven effective in generating high accuracy images [8

8. H. Dehghani, B. W. Pogue, S. P. Poplack, and K. D. Paulsen, “Multiwavelength three-dimensional near-infrared tomography of the breast: initial simulation, phantom, and clinical results,” Appl. Opt. 42, 135–145 (2004). [CrossRef]

, 13

13. G. Boverman, Q. Fang, S. A. Carp, E. L. Miller, D. H. Brooks, J. Selb, R. H. Moore, D. B. Kopans, and D. A. Boas, “Spatio-temporal imaging of the hemoglobin in the compressed breast with diffuse optical tomography,” Phys. Med. Biol. 52, 3619–3641 (2007). [CrossRef] [PubMed]

]. To this end hyperspectral information is employed, which along with regularization has been shown to reduce the non-uniqueness of the solution to the image formation problem [8

8. H. Dehghani, B. W. Pogue, S. P. Poplack, and K. D. Paulsen, “Multiwavelength three-dimensional near-infrared tomography of the breast: initial simulation, phantom, and clinical results,” Appl. Opt. 42, 135–145 (2004). [CrossRef]

, 10

10. F. Larusson, S. Fantini, and E. L. Miller, “Hyperspectral image reconstruction for diffuse optical tomography,” Biomed. Opt. Express 2, 947–965 (2011). [CrossRef]

, 11

11. B. W. Pogue, T. O. McBride, J. Prewitt, U. L. Osterberg, and K. D. Paulsen, “Spatially variant regularization improves diffuse optical tomography,” Appl. Opt. 38, 2950–2961 (1999). [CrossRef]

].

The primary contribution of this work is the development of a new approach to the multi-chromophore inverse problem in which the medium is or can be well approximated as piecewise constant. Piecewise constant DOT problems have been considered in the past. In [18

18. O. Dorn, “A shape reconstruction method for diffuse optical tomography using a transport model and level sets,” in 2002 IEEE International Symposium on Biomedical Imaging, 2002. Proceedings (IEEE, 2002), pp. 1015–1018.

] level sets were used in a two-step method for shape estimation assuming that prior information of the absorption parameter was known. Schweiger et al. [19

19. M. Schweiger, O. Dorn, and S. R. Arridge, “3-D shape and contrast reconstruction in optical tomography with level sets,” J. Phys.: Conf. Ser. 124, 012043 (2008). [CrossRef]

] and Kilmer et al. [20

20. M. E. Kilmer, E. L. Miller, A. Barbaro, and David Boas, “Three-dimensional shape-based imaging of absorption perturbation for diffuse optical tomography,” Appl. Opt. 42, 3129–3144 (2003). [CrossRef] [PubMed]

] employed level sets for the DOT problem estimating parameter distributions using a piecewise basis. Arridge et al. investigated shape based methods by estimating level-sets, specifically investigating an explicit method using basis functions and an implicit shape reconstruction to recover absorption and diffusion assuming a known background [21

21. S. R. Arridge, O. Dorn, V. Kolehmainen, M. Schweiger, and A. Zacharopoulos, “Parameter and structure reconstruction in optical tomography,” J. Phys.: Conf. Ser. 135, 012001 (2008). [CrossRef]

]. For a detailed review of the use of level sets in inverse scattering problems we point the reader to [22

22. O. Dorn and D. Lesselier, “Level set methods for inverse scattering,” Inverse Probl. 22, R67–R131 (2006). [CrossRef]

].

The approach we consider in this paper is significantly different from those in [18

18. O. Dorn, “A shape reconstruction method for diffuse optical tomography using a transport model and level sets,” in 2002 IEEE International Symposium on Biomedical Imaging, 2002. Proceedings (IEEE, 2002), pp. 1015–1018.

, 19

19. M. Schweiger, O. Dorn, and S. R. Arridge, “3-D shape and contrast reconstruction in optical tomography with level sets,” J. Phys.: Conf. Ser. 124, 012043 (2008). [CrossRef]

, 21

21. S. R. Arridge, O. Dorn, V. Kolehmainen, M. Schweiger, and A. Zacharopoulos, “Parameter and structure reconstruction in optical tomography,” J. Phys.: Conf. Ser. 135, 012001 (2008). [CrossRef]

]. In addition to the fact that none of these papers have considered the fully hyperspectral case, some of these methods [18

18. O. Dorn, “A shape reconstruction method for diffuse optical tomography using a transport model and level sets,” in 2002 IEEE International Symposium on Biomedical Imaging, 2002. Proceedings (IEEE, 2002), pp. 1015–1018.

20

20. M. E. Kilmer, E. L. Miller, A. Barbaro, and David Boas, “Three-dimensional shape-based imaging of absorption perturbation for diffuse optical tomography,” Appl. Opt. 42, 3129–3144 (2003). [CrossRef] [PubMed]

] require the recovery of unknown quantities defined on a fine scale pixelated discretization of the region of interest. More specifically in [21

21. S. R. Arridge, O. Dorn, V. Kolehmainen, M. Schweiger, and A. Zacharopoulos, “Parameter and structure reconstruction in optical tomography,” J. Phys.: Conf. Ser. 135, 012001 (2008). [CrossRef]

] absorption and scattering is estimated using level sets assuming the background known. With the Born approximation we assume the absorption and diffusion is known in the background but here we estimate chromophore concentrations and diffusion of the object of interest as well as chromophore concentration in the background. TV reconstructions recover images while traditional level set methods work with a level set function defined on a pixel-based grid. In both cases, regularization is required to obtain adequate results and one is faced with the corresponding challenge of choosing regularization parameters [11

11. B. W. Pogue, T. O. McBride, J. Prewitt, U. L. Osterberg, and K. D. Paulsen, “Spatially variant regularization improves diffuse optical tomography,” Appl. Opt. 38, 2950–2961 (1999). [CrossRef]

, 19

19. M. Schweiger, O. Dorn, and S. R. Arridge, “3-D shape and contrast reconstruction in optical tomography with level sets,” J. Phys.: Conf. Ser. 124, 012043 (2008). [CrossRef]

].

In this paper we consider the use of a shape-based approach to the hyperspectral DOT problem based on a newly-developed parametric level set (PaLS) formulation. In [15

15. A. Aghasi, M. Kilmer, and E. L. Miller “Parametric level set methods for inverse problems,” SIAM J. Imaging Sci. 4, 618–650 (2011). [CrossRef]

], a basis function expansion was used to provide a low order representation of the level set function and yielded very strong results for a number of highly ill-posed inverse problems including a restricted form of the DOT problem where a single wavelength was employed to determine only optical absorption. The method in [15

15. A. Aghasi, M. Kilmer, and E. L. Miller “Parametric level set methods for inverse problems,” SIAM J. Imaging Sci. 4, 618–650 (2011). [CrossRef]

] required no explicit regularization and, due to the low-order nature of the model (number of parameters much, much less than number of pixels) was amenable to Newton-type inversion algorithms known to converge more rapidly than gradient-based schemes. Moreover, in [15

15. A. Aghasi, M. Kilmer, and E. L. Miller “Parametric level set methods for inverse problems,” SIAM J. Imaging Sci. 4, 618–650 (2011). [CrossRef]

] our experiments indicated a surprising lack of sensitivity to the initialization of the inversion algorithm.

Although the breast is a highly heterogeneous medium, in the level set method we assume the images to be recovered to be piecewise constant. This approximation is supported in the literature. For example Schweiger et al. assumed anatomical prior information to derive a piecewise constant region basis [4

4. M. Schweiger and S. R. Arridge, “Optical tomographic reconstruction in a complex head model using a priori region boundary information,” Phys. Med. Biol. 44, 2703–2721 (1999). [CrossRef] [PubMed]

]. The co-located structures of different species of chromophores is discussed further in Section 5. In the future, where we move to more complicated simulations and experiments the goal would be to implement multiple characteristic functions that require the estimation of multiple level sets, discussed further in Section 8. To prove the accuracy of our method, we use simulation data to show parametric level-set reconstructions for multiple chromophores and diffusion amplitude. Experimental data are used to show the improvement of the PaLS method over pixel based methods using previously reported data sets [10

10. F. Larusson, S. Fantini, and E. L. Miller, “Hyperspectral image reconstruction for diffuse optical tomography,” Biomed. Opt. Express 2, 947–965 (2011). [CrossRef]

].

The remainder of this paper is structured as follows. In Section 2 we discuss our forward model and derive in detail the Born model used. In Section 3 we discuss the PaLS method and in Section 4 we derive the Levenberg Marquardt algorithm used for the inverse problem. In Section 7.1 we present simulation results with multiple chromophore and diffusion amplitude reconstructions, and in Section 7.2 we present experimental results that show the improvement of the PaLS method over pixel based methods.

2. Forward problem

In this paper we restrict our attention to problems in which the transport physics [2

2. R. J. Gaudette, D. H. Brooks, C. A. DiMarzio, M. E. Kilmer, E. L. Miller, T. Gaudette, and D. A. Boas, “A comparison study of linear reconstruction techniques for diffuse optical tomographic imaging of absorption coefficient,” Phys. Med. Biol 45, 1051–1069 (2000). [CrossRef] [PubMed]

] associated with the propagation of light at wavelength λ through tissue can be adequately approximated using a diffusion model [23

23. M. A. O’Leary, D. A. Boas, B. Chance, and A. G. Yodh, “Experimental images of heterogeneous turbid media by frequency-domain diffusing-photon tomography,” Opt. Lett. 20, 426–428 (1995). [CrossRef]

, 24

24. B. Brooksby, B. W. Pogue, S. Jiang, H. Dehghani, S. Srinivasan, C. Kogel, T. D. Tosteson, J. Weaver, S. P. Poplack, and K. D. Paulsen, “Imaging breast adipose and fibroglandular tissue molecular signatures by using hybrid mri-guided near-infrared spectral tomography,” Proc. Natl. Acad. Sci. U.S.A. 103, 8828–8833 (2006). [CrossRef] [PubMed]

] of the form
D0(r,λ)Φ(r,λ)+vμa0(r,λ)Φ(r,λ)=vS(r,λ)
(1)
where Φ(r,λ) is the photon fluence rate at position r due to light of wavelength λ injected into the medium, v is the electromagnetic propagation velocity in the medium, μa0(r,λ) is the spatially varying absorption coefficient, and S(r,λ) is the photon source with units of optical energy per unit time per unit volume. For the work in this paper the sources are considered to be delta sources in space and can be written as S(r,λ) = S0(λ)δ(rrs) with S0(λ) the source power at wavelength λ. For spatially varying scattering we assume that the diffusion coefficient D0(r,λ) follows Mie scattering theory where a scattering prefactor Ψ depends on the size and density of scatterers while a scattering exponent b depends on the size of the scatterers. Using this, we write the perturbation in the diffusion coefficient as
ΔD(r,λ)=v3ΔΨ(λλ0)b=vΔΨd(λ).
(2)
The arbitrarily chosen reference wavelength λ0 is introduced to achieve a form of the Mie model where Ψ has the units of mm−1 and Ψ′ has units of mm. In this paper, for simplicity we consider the case where the background medium is infinite and homogeneous. Generalization to the more practical case where there are boundaries is straightforward in theory though somewhat more complex in terms of implementation [25

25. D. Boas, M. A. O’Leary, B. Chance, and A. Yodh, “Scattering of diffuse photon density waves by spherical inhomogeneities within turbid media: analytic solution and applications,” Proc. Natl. Acad. Sci. U.S.A. 91, 4887–4891 (1994). [CrossRef] [PubMed]

]. As the primary objective of this work is the demonstration of a new approach to inversion, we prefer to consider the simpler physical problem holding off on the more realistic implementation for the future. Additionally, we only consider continuous wave experiments since this is what our instrument measures, where Eq. (1) is usually written with a /D(λ) term on the right hand side, in our case we consider ω = 0 giving us the form shown in Eq. (1).

We employ the Born approximation by decomposing μa0(r,λ) and D0(r,λ), as the sum of a constant background absorption, μa(λ), and a spatially varying perturbation Δμa(r,λ) as well as constant background diffusion D(λ) and a diffusion perturbation ΔD(r,λ). To obtain a linear relationship between the measurements and the chromophore concentrations, we subtract Eq. (1) from the perturbed version of the diffusion equation which gives
[2+k02(λ)]Φs(r,λ)=Δk2(r,λ)Φ(r,λ)1D(λ)ΔD(r,λ)Φ
(3)
where k02(λ)=vμa(λ)/D(λ) and Δk2(r,λ) = (v/D(λ))Δμa(r,λ). Assuming the availability of a Green’s function, G(r,r′,λ) for the solution of Eq. (3) as is the case for an unbounded medium as well as range of bounded problems [26

26. A. Mandelis, Diffusion-Wave Fields: Mathematical Methods and Green Functions (Springer, 2001).

], we can change Eq. (3) into an integral Eq. [27

27. B. Brendel, R. Ziegler, and T. Nielsen “Algebraic reconstruction techniques for spectral reconstruction in diffuse optical tomography,” Appl. Opt. 47, 6392–6403 (2008). [CrossRef] [PubMed]

]
Φs(r,λ)=vD(λ)(Δμ(r,λ)G(r,r,λ)Φi(r,λ)dr+ΔD(r,λ)G(rd,rj,λ)Φi(rj,rs,λ)dr)
(4)
where rd is the location of the detector and (with a small abuse of notation) Φi(r,rs,λ) is used here to denote the incident field at position r and wavelength λ due to a delta-source located at rs. It should also be noted that we obtain this equation under the assumption that the total fluence rate, Φ, can be approximated as the incident field, Φi, since Φi ≫ Φs [2

2. R. J. Gaudette, D. H. Brooks, C. A. DiMarzio, M. E. Kilmer, E. L. Miller, T. Gaudette, and D. A. Boas, “A comparison study of linear reconstruction techniques for diffuse optical tomographic imaging of absorption coefficient,” Phys. Med. Biol 45, 1051–1069 (2000). [CrossRef] [PubMed]

].

Equation (4) provides a linear relationship between the scattered fluence rate and the absorption perturbation. To relate the scattered fluence rate to concentrations of chromophores, Δμa is decomposed as follows [28

28. A. Corlu, R. Choe, T. Durduran, K. Lee, M. Schweiger, S. R. Arridge, E. M. Hillman, and A. G. Yodh, “Diffuse optical tomography with spectral constraints and wavelength optimization,” Appl. Opt. 44, 2082–2093 (2005). [CrossRef] [PubMed]

]
Δμa(r,λ)=k=1Nsεk(λ)ck(r).
(5)
where Ns is the number of absorbing species for the problem under investigation, εk(λ) is the extinction coefficient for the kth species at wavelength λ, and ck(r) is the concentration of species k at location r. To obtain the fully discrete form of the Born model used in Section 4, we expand each ck(r)
ck(r)=j=1Nvck,jφ(r)
(6)
where ck,j is the value of the concentration for species k in Vj, the jth “voxel”. The φ(r) function is an indicator function where
φ(r)={1,ifrVj0,ifrVj.
(7)
After using Eqs. (5) and (6) in Eq. (4) and performing some algebra we obtain
Φs(λ)=avD(λ)k,j=1Nc,NV(G(rd,rj,λ)Φi(rj,rs,λ)εi(λ)ck,j+G(rd,rj,λ)Φi(rj,rs,λ)ΔDj(λ)).
(8)
We approximate Eq. (4) as the value at the center of each pixel multiplied by the area or volume of each pixel or voxel, so in Eq. (8) we use a as the area of a pixel. This setup is illustrated in Fig. 1(a).

Fig. 1 (a) The setup of sources and detectors for simulation reconstructions. Same orientation of axes is used for experimental data. (b) Definition of domains used for the parametric level-set methods.

Setting up the linear algebraic structure associated with Eq. (8) we define ck ∈ 𝕉Nv as the vector obtained by lexicographically ordering the unknown concentrations associated with the kth chromophore and ck+1 = ΔΨ′ and Φs(λ) is the vector of observed scattered fluence rate associated with all source-detector pairs collecting data at wavelength λ. Now, with Nλ the number of wavelengths used in a given experiment, Eq. (8) is written in matrix-vector notation as
[Φs(λ1)Φs(λ2)Φs(λNλ)]=[ε1(λ1)K1aε2(λ1)K1aεNc(λ1)K1aK1dε1(λ2)K2aε2(λ2)K2aεNc(λ2)K2aK2dε1(λNλ)KNλaε2(λNλ)KNλaεNc(λNλ)KNλaKNλd][c1c2cNcΔΨ]Φs=Kc
(9)
It should be noted in Eq. (9) that the matrix has elements which are also the matrices Kla and Kld. The (m, j)th element of the Kla matrix is given by the product (v/D(λl))G(rm, rj, λli(rj, rm, λl), where m represents the mth source-detector pair and as before j represents the jth voxel. For Kld the matrix elements are given by (v/D(λl))∇G(rm, rj, λ) ·∇Φi(rj, rm, λDj(λ).

To evaluate the forward model for realistically sized problems, we compute the Nλ matrices Ka and Kd and store it along with the Nλ × Nc extinction coefficients. This reduces the amount of memory required for the reconstruction.

3. Parametric level-set method

To counter the ill-posedness of the DOT problem we employ a Parametric Level-Set Method (PaLS). For the purpose of our paper we assume that all chromophore concentrations and diffusion coefficient perturbations are co-located. This choice is supported by reports in literature, where decrease in hemoglobin and water concentration along with scattering power are located at the cancer location and the lipid concentrations increase at the same location [29

29. A. Cerussi, D. Hsiang, N. Shah, R. Mehta, A. Durkin, J. Butler, and B. J. Tromberg “Predicting response to breast cancer neoadjuvant chemotherapy using diffuse optical spectroscopy,” Proc. Nat. Acad. Sci. U.S.A. 104, 4014–4019 (2007). [CrossRef]

31

31. Q. Zhu, P. U. Hegde, A. Ricci, M. Kane, E. B. Cronin, Y. Ardeshirpour, C. Xu, A. Aguirre, S. H. Kurtzman, P. J. Deckers, and S. H. Tannenbaum, “Early-stage invasive breast cancers: potential role of optical tomography with US localization in assisting diagnosis,” Radiology 256, 367–378 (2010). [CrossRef] [PubMed]

]. This means that the geometry of the anomaly in the medium is the same for all chromophores and diffusion amplitude. The domain Ω ⊂ 𝒡 represents the support of the objects of interest, and 𝒡 represents the homogeneous background where the anomaly is located, shown in Fig. 1(b).

Since the same support is used for each object of interest the characteristic function describing the shape is defined as
χ(x,y)={1if(x,y)Ω,0if(x,y)\Ω.
(10)
Then each image to be reconstructed can be written as
ck(x,y)=χ(x,y)cka+[1χ(x,y)]ckb
(11)
where k = 1, 2,...,Nc + 1. In this formulation the unknown values are the constant concentration values of the anomaly and background, cka and ckb respectively.

The characteristic function χ(x,y) is defined to be the zero level set of a Lipschitz continuous object function 𝒪 : 𝒡 → 𝕉 such that 𝒪 > 0 in Ω(x,y), 𝒪 < 0 in Ω\𝒡 and 𝒪(x,y) = 0 in Ω. Using 𝒪(x,y), χ(x,y) is written as
χ(x,y)=H(𝒪(x,y))
(12)
where H is the step function. In practice we use smooth approximations of the step function and the Dirac delta function denoted as Hε and δε respectively where Hε is computed as
Hε(x)={1ifx>ε,0ifx<ε.12[1+xε+1πsin(πxε)]if|x|ε.
and δε is computed as the derivative of Hε [32

32. T. Chan and L. Vese “Active contours without edges,” IEEE Trans. Image Process. 10, 266–277 (2001). [CrossRef]

, 33

33. H. K. Zhao, S. Osher, B. Merriman, and M. Kang, “Implicit, nonparametric shape reconstruction from unorganized points using a variational level set method,” Comput. Vision Image Understanding 80, 295–314 (2000). [CrossRef]

]. As discussed in Section 1 we represent the object function 𝒪(x,y) parametrically, so instead of using a dense collection of pixel or voxel values [34

34. S. Osher and R. Fedkiw, Level Set Methods and Dynamic Implicit Surfaces (Springer, 2002)

], we represent it by using basis functions
𝒪(x,y)=i=1Laipi(x,y)
(13)
where ai’s are the weight coefficients whereas pi(x,y) are the functions which belong to the basis set of 𝒫 = {p1, p2,..., pl}. Possible choices for the 𝒫 basis set include polynomial or radial basis functions. For the purpose of this paper we use Gaussian basis function. The width and number of the Gaussians determines how coarse or fine the reconstruction will be, where a choice of few basis functions will, on the one hand, result in a reduced number of unknowns, it will on the other hand, give a coarser estimation of the shape, which can be a problem for imaging finer more complex structures. In the DOT case, where the physics in the forward model will only allow for a coarse reconstruction of the underlying structure, the Gaussian approach is sufficient, especially for the relatively simple geometries and concentrations presented in this paper. When we will move on to a more complicated non-linear model, the choice of the number and type of basis functions will be more important and should be based on a rigid mathematical framework. This is discussed further in Section 8.

In the PaLS formulation, all the parameters of the model are gathered in one vector θT=[ca1,ca2,,cb1,aT] where a = [a1,..., aL]T. Now our forward model in Eq. (9) can be expressed as
Φs=K(θ)=Kc(θ)
(14)
The forward model has now been parametrized with a vector containing all of the unknowns, which are far fewer then what a pixel based method would attempt to estimate.

4. Inverse problem

The inverse problem, that of using Φs to recover the value of c, is solved as a Levenberg-Marquardt optimization problem of the form
c^=argmincW(K(θ)Φs)22
(15)
The W matrix reflects the structure of the noise corrupting the data [2

2. R. J. Gaudette, D. H. Brooks, C. A. DiMarzio, M. E. Kilmer, E. L. Miller, T. Gaudette, and D. A. Boas, “A comparison study of linear reconstruction techniques for diffuse optical tomographic imaging of absorption coefficient,” Phys. Med. Biol 45, 1051–1069 (2000). [CrossRef] [PubMed]

]. While a Poisson model is technically the most appropriate for DOT data, as is frequently done [35

35. M. Guven, B. Yazici, X. Intes, and B. Chance, “Diffuse optical tomography with a priori anatomical information,” Phys. Med. Biol. 50, 2837–2858 (2005). [CrossRef] [PubMed]

] we employ a Gaussian approximation in which independent, zero mean Gaussian noise is assumed to corrupt each datum. The reason for this is that with a sufficiently large number of detected photons, the Poisson statistics can be approximated by a Gaussian distribution [11

11. B. W. Pogue, T. O. McBride, J. Prewitt, U. L. Osterberg, and K. D. Paulsen, “Spatially variant regularization improves diffuse optical tomography,” Appl. Opt. 38, 2950–2961 (1999). [CrossRef]

]. Letting σm2 denote the variance of the noise corrupting the mth elements of Φ, W is constructed as a diagonal matrix with 1/σm the mth element along the diagonal. For the experimental and simulated data the variance is calculated from
σm2=Ω(m)10SNRm10.
(16)
where Ω(m) corresponds to the photon count for each source-detector pair. The SNR for each element of Φ is then calculated from
SNRm=10log10(Ω(m)/Ω(m)).
(17)
In experimental data Ω(m) is the standard deviation of the Poisson noise distribution. The minimization of the cost function is then achieved by the Levenberg-Marquardt algorithm. For that purpose an error vector is introduced
ε=W(K(θ)Φ).
(18)
which can be used to write the cost function in term of ε as
M(θ)=εTε
(19)

When estimating the parametric vector, we employ a cyclic coordinate decent strategy [38

38. T. T. Wu and K. Lange, “Coordinate descent algorithms for lasso penalized regression,” Ann. Appl. Stat. 2, 224–244 (2008). [CrossRef]

] Essentially this is equivalent to estimating the shape only at even iterations and the concentration values at odd iterations. This is repeated until stopping criteria is reached. This process is expressed in pseudo-code in Algorithm 1, where Jv and Js denote the Jacobian strictly for the concentration value and shape, respectively, and τi represents a tolerance for the stopping critera.

Algorithm 1. Matlab-like code for estimating shape and concentration value

table-icon
View This Table
| View All Tables

5. Simulation analysis

The arrangement of sources and detectors with respect to the simulated turbid medium is displayed in Fig. 1(a). This arrangement of sources and detectors is chosen to represent a common setup in optical mammography where the breast is compressed between two planes containing sources and detectors [28

28. A. Corlu, R. Choe, T. Durduran, K. Lee, M. Schweiger, S. R. Arridge, E. M. Hillman, and A. G. Yodh, “Diffuse optical tomography with spectral constraints and wavelength optimization,” Appl. Opt. 44, 2082–2093 (2005). [CrossRef] [PubMed]

]. The source-detector separations were set to 5 cm as is shown in Fig. 1(a). In simulations, we reconstruct concentration images of oxygenated and deoxygenated hemoglobin, HbO2 and HbR respectively, along with lipid and water concentration and diffusion amplitude. These chromophores are chosen since they mainly cause near-infrared absorption in breast tissue [39

39. B. Brendel and T. Nielsen, “Selection of optimal wavelengths for spectral reconstruction in diffuse optical tomography,” J. Biomed. Opt. 14, 034041 (2009). [CrossRef] [PubMed]

], and breast cancer tumours have been found to have higher HbO2 and HbR concentrations than normal tissue [40

40. S. Srinivasan, B. W. Pogue, S. Jiang, H. Dehghani, C. Kogel, S. Soho, J. J. Gibson, T. D. Tosteson, S. P. Poplack, and K. D. Paulsen, “Interpreting hemoglobin and water concentration, oxygen saturation, and scattering measured in vivo by near-infrared breast tomography,” Proc. Natl. Acad. Sci. U.S.A. 100, 12349–12354 (2003). [CrossRef] [PubMed]

].

In simulations, values for Ψ and b are obtained from [41

41. P. Taroni, A. Pifferi, A. Torricelli, D. Comelli, and R. Cubeddu, “In vivo absorption and scattering spectroscopy of biological tissues,” Photochem. Photobiol. Sci. 2, 124–129 (2003). [CrossRef] [PubMed]

] for the female breast. Values for μa are calculated from extinction coefficients, which are in the units cm−1/mM and are obtained from data tabulated by Scott Prahl [42

42. S. Prahl, “Tabulated molar extinction coefficient for hemoglobin in water” (Oregon Medical Laser Center, 2007), http://omlc.ogi.edu/spectra/hemoglobin/summary.html.

]. The concentration in the simulated images are defined in units of millimolars or millimoles per liter, mM for HbO2 and HbR. For water and lipid the concentrations are in percent by volume and the diffusion amplitude is measured in units of millimeter. The background has HbR concentration of 0.01 mM, HbO2 concentration of 0.01 mM, lipid concentration of 32%, water concentration of 13% and Ψ′ is set to 1.6 mm. The target concentration of the object of interest is set to 0.015 mM, 0.012 mM, 50 %, 20 % and 0.25 mm for Hb02, HbR, lipid, water and ΔΨ′, respectively.

The simulation set is created with all chromophore concentrations and diffusion perturbations in the same location with different target values. The co-located geometry is chosen since it is likely for chromophore concentrations to have similarly located geometries in the real breast [29

29. A. Cerussi, D. Hsiang, N. Shah, R. Mehta, A. Durkin, J. Butler, and B. J. Tromberg “Predicting response to breast cancer neoadjuvant chemotherapy using diffuse optical spectroscopy,” Proc. Nat. Acad. Sci. U.S.A. 104, 4014–4019 (2007). [CrossRef]

, 30

30. H. Soliman, A. Gunasekara, M. Rycroft, J. Zubovits, R. Dent, J. Spayne, M. J. Yaffe, and G. Czarnota, “Functional imaging using diffuse optical spectroscopy of neoadjuvant chemotherapy response in women with locally advanced breast cancer,” Clin. Cancer Res. 16, 2605–2614 (2010). [CrossRef] [PubMed]

]. The ground truth images for simulations are shown in Fig. 6. Reconstruction is done for these images to explore effects of adding hyperspectral information to the problem, i.e. the improvement in quantitative accuracy and the reduction of crosstalk where a concentration of one chromophore creates a false concentration in an image for another chromophore as well as the performance of the shape based approach. The process is initialized with 21 Gaussian basis functions with width of approximately 8 pixels placed uniformly on a grid over the whole medium to be imaged. A representative image of the order of the basis functions is shown in Fig. 2(a). For all experiments presented in this paper, the ai’s weight coefficients are initialized to 1.

Fig. 6 Reconstruction using the PaLS method. Middle column of images are generated with 8 wavelengths and rightmost images are generated with 176 wavelengths. From top to bottom the rows show HbO2, HbR, lipid, water and diffusion amplitude, respectively. Concentration units are in mM.
Fig. 2 (a) Order of the basis functions used for the PaLS method. (b) Absorption spectra of the ink and dye solutions chromophores used in experimental measurements. Specifically chosen wavelengths are marked with an asterisk.

To best understand the utility of a hyperspectral data set, we employ the Born model to generate simulated data. Though this may not be realistic, it allows us to avoid the confounding factor of model mismatch in evaluating the inversion method being considered in this paper. Moreover, the shortcomings of this approach are mitigated in Section 7.2, where we consider the processing of experimental data which, obviously, are not the product of the Born model. In a bit more detail, the data we use for our simulation analysis are computed as
Φ=Kc+n
(22)
where c are the simulated concentration images for all chromophores and diffusion amplitude, whereas n represents additive noise. Specifically, as in [2

2. R. J. Gaudette, D. H. Brooks, C. A. DiMarzio, M. E. Kilmer, E. L. Miller, T. Gaudette, and D. A. Boas, “A comparison study of linear reconstruction techniques for diffuse optical tomographic imaging of absorption coefficient,” Phys. Med. Biol 45, 1051–1069 (2000). [CrossRef] [PubMed]

] n is a vector of zero mean, independent Gaussian random variables with variances σm2, defined in Eq. (16), chosen such that a pre-determined signal to noise ratio (SNR) is achieved. This SNR is calculated from Eq. (17).

The reconstructed images are evaluated in three ways: through visual inspection, using mean square error (MSE) as a measure of overall quantitative accuracy for each chromophore, and examining the Dice coefficient to judge how well the concentrations are located. For the kth chromophore, the mean square error is computed by using the following Eq.
MSEk=ckc^k2ck2
(23)

Along with using MSE to judge the accuracy of the reconstruction, we also employ a Dice Coefficient to quantitatively analyze how well the shape based approach localizes the reconstruction [43

43. A. A. Joshi, A. J. Chaudhari, D. W. Shattuck, J. Dutta, R. M. Leahy, and A. W. Toga, “Posture matching and elastic registration of a mouse atlas to surface topography range data,” in IEEE International Symposium on Biomedical Imaging: from Nano to Macro, 2009. ISBI ’09 (IEEE, 2009), pp. 366–369 (2009).

, 44

44. J. D. Vylder and W. Philips, “A computational efficient external energy for active contour segmentation using edge propagation,” in IEEE 2100 International Conference on Image Processing (ICIP 2010) (IEEE, 2010), pp. 661–664.

]. If S is the reconstructed image and G is the ground truth created for each set, the Dice coefficient between S and G can be defined as
D(S,G)=2|SG||S|+|G|
(24)
|SG| contains all pixels that belong to the detected segment as well as the ground truth segment, so that if S and G are equal the Dice coefficient is equal to one, indicating an accurate reconstruction. To compute the D(S,G) we use the characteristic function, χ, which essentially works as a binary map of the reconstructed anomaly where the object of interest is represented by 1’s.

6. Experimental analysis

Physical measurements were performed in order to validate the simulation results, where the data used in this paper has previously been used in [10

10. F. Larusson, S. Fantini, and E. L. Miller, “Hyperspectral image reconstruction for diffuse optical tomography,” Biomed. Opt. Express 2, 947–965 (2011). [CrossRef]

]. We use the same data again since we are using a completely new image formation process which results in a significant improvement in reconstructions. The background medium consists of water and milk in the ratio of 2:1, respectively. Milk, with 2% fat, is used due to the similarities of the optical properties to human breast tissue. Black India ink and blue food dye were added to mimic tissue chromophores. The ink and dye are mixed into the background of milk and water to achieve μa = 0.029 cm−1, at 600 nm, which is in the range of optical absorption of the female breast [45

45. T. Durduran, R. Choe, J. P. culver, L. Zubkov, M. J. Holboke, J. Giammarco, B. Chance, and A. G. Yodh, “Bulk optical properties of healthy female breast tissue,” Phys. Med. Biol. 47, 2847–2861 (2002). [CrossRef] [PubMed]

, 46

46. L. Spinelli, A. Torricelli, A. Pifferi, P. Taroni, G. M. Danesini, and R. Cubeddu, “Bulk optical properties and tissue components in the female breast from multiwavelength time-resolved optical mammography,” J. Biomed. Opt. 9, 1137–1142 (2004). [CrossRef] [PubMed]

].

The absorption spectra for the ink and dye inclusions, shown in Fig. 2(b), have the most significant effect in the 450–700 nm range. These chromphores are chosen because the spectral shapes of their absorption are similar to those of HbO2 and HbR and have been widely used in literature [47

47. J. P. Culver, R. Choe, M. J. Holboke, L. Zubkov, T. Durduran, A. Slemp, V. Ntziachristos, B. Chance, and A. G. Yodh, “Three-dimensional diffuse optical tomography in the parallel plane transmission geometry: evaluation of a hybrid frequency domain/continuous wave clinical system for breast imaging,” Med. Phys. 30, 235–247 (2003). [CrossRef] [PubMed]

, 48

48. N. Liu, A. Sassaroli, and S. Fantini, “Paired-wavelength spectral approach to measuring the relative concentrations of two localized chromophores in turbid media: an experimental study,” J. Biomed. Opt. 12, 051602 (2007). [CrossRef] [PubMed]

]. Therefore, the 6 specifically chosen wavelengths were selected from this region where these chromophores had the most significant absorption.

In order to obtain multi- and hyperspectral reconstruction values for μa and D(λ) the background has to be known. In the experimental measurements we assume constant scattering, therefore we are not trying to estimate the perturbation, ΔD(r,λ), as was done in simulations. Therefore we have the unperturbed representaion of the reduced scattering coefficient, μs, is given by
μs=Ψ(λλ0)b.
(25)
The diffusion coefficient relates to the reduced scattering coefficient through, D(λ) = v/3μs. Phase, amplitude and average intensity data are obtained at two wavelengths using a frequency-domain tissue spectrometer to estimate the Ψ and b parameters in Eq. (25) as Ψ = 6.5 cm−1 and b = 0.4. This allows us to extrapolate values for μs at any wavelength [49

49. E. Gratton, S. Fantini, M. A. Franceschini, C. Gratton, and M. Fabiani, “Measurements of scattering and absorption changes in muscle and brain,” Philos. Trans. R. Soc. Lond. B. Biol. Sci. 352, 727–735 (1997). [CrossRef] [PubMed]

]. Determining spectrally extrapolated values for the absorption coefficient is harder. Since μa does not follow a law like μs, values are estimated using extinction coefficient data for ink, dye, milk and water. These extinction coefficient are measured in a standard spectrophotometer.

In our experiments, two phantom inclusions, named set 1 and set 2, are created for different absorption contrasts relative to the background in the range of 3:1 to 1:1. The inclusion in set 1 contains 10% ink and 90% dye and the inclusion for set 2 contains 70% dye and 30% ink. This contrast range is comparable to traditional tumor contrasts reported in literature, which have been close to 3:1 and lower [50

50. B. W. Pogue, S. Jiang, H. Dehghani, C. Kogel, S. Soho, S. Srinivasan, X. Song, T. D. Tosteson, S. P. Poplack, and K. D. Paulsen, “Characterization of hemoglobin, water, and nir scattering in breast tissue: analysis of intersubject variability and menstrual cycles changes,” J. Biomed. Opt. 9, 541–552 (2004). [CrossRef] [PubMed]

]. Reconstructions are done for 126 wavelengths equally spaced over the whole spectrum and 6 specifically chosen wavelengths as λ = [480, 550, 610, 630, 650, 690] nm. The wavelengths are chosen around the isosbestic point, where the contrast between the chromophores is the highest and where each chromophore has highest absorption. The absorption spectra and the contrast over the spectrum for set 1 and set 2 are shown in Fig. 3 and Fig. 4, respectively.

Fig. 3 (a) Absorption spectra for the background, μa, and the inclusion, μa + Δμa, in experimental set 1, containing 10% ink and 90% dye. (b) Contrast between the background and the inclusion for experimental set 1.
Fig. 4 (a) Absorption spectra for the background, μa, and the inclusion, μa + Δμa, in experimental set 2, containing 70% ink and 30% dye. (b) Contrast between the background and the inclusion for experimental set 2.

For both experimental sets 1 and 2, a single 25 cm long transparent cylindrical inclusion containing ink and dye solutions is placed in the background medium. Optical properties are assumed constant along the z-axis. For the light source, an arc lamp is used which delivers light with an average illumination power of 280 mW, that translates into a power density of 3.96 W/cm2. For three different locations, at x ± 1 cm, a 5 mm diameter collection optical glass fiber bundle is placed on the opposite side of the inclusions (at a y-axis separation of 5 cm). Experiments are made with the light source placed in succession at 8 positions with 1 cm increments for total of 24 source-detector pairs. To ensure that approximately the same amount of photons are collected for both hyperspectral and multispectral reconstructions, two exposure times are used for the CCD camera, a longer one of 10 s for the 6 wavelength case and 500 ms for the 126 wavelength case. Since the goal of this paper is to demonstrate the improvement of including hyperspectral information, we present an ideal case where the signal to noise ratio is large, thereby providing a best-case scenario for the few wavelength reconstruction against which we compare our approach as well as using realistic absorption contrasts for the inclusions. Further details on the experimental setup can be found at [10

10. F. Larusson, S. Fantini, and E. L. Miller, “Hyperspectral image reconstruction for diffuse optical tomography,” Biomed. Opt. Express 2, 947–965 (2011). [CrossRef]

].

In our experimental setup, we measure the incident field before the perturbation is put into the medium. Then the scattered field is computed as a dataset that has the unperturbed dataset subtracted from it. For in vivo measurements it is possible to use a priori structural information from other modalities, e.g. MRI, to estimate the incident field by determining the optical properties of the assumed piecewise constant chromophore distribution over these segments [51

51. G. Boverman, E. L. Miller, A. Li, Q. Zhang, T. Chaves, D. H. Brooks, and D. A. Boas, “Quantitative spectroscopic diffuse optical tomography of the breast guided by imperfect a priori structural information,” Phys. Med. Biol. 50, 3941–3956 (2005). [CrossRef] [PubMed]

].

A comparison of the absolute concentrations, ci^ and relative concentration, cir^ to target concentration values is done to test the accuracy of the reconstructions. The relative concentrations for ink are calculated as
cinkr^=cink^/(cink^+cdye^)
(26)
and similarly for dye [52

52. S. Fantini, M. A. Franceschini, J. S. Maier, S. A. Walker, B. Barbieri, and E. Gratton, “Frequency-domain multi-channel optical detector for noninvasive tissue spectroscopy and oximetry,” Opt. Eng. 34, 32–42 (1995). [CrossRef]

,53

53. M. A. Franceschini, V. Toronov, M. E. Filiaci, E. Gratton, and S. Fantini, “On-line optical imaging of the human brain with 160-ms temporal resolution,” Opt. Express 6, 49–57 (2000). [CrossRef] [PubMed]

]. The relative concentration is calculated from the peak concentration value in each reconstruction. This allows us to inspect how well our approach manages to separate and estimate each species of chromophores in the process.

7. Results

7.1. Simulations

In Fig. 5 reconstruction results using the pixel based method are shown for 8 wavelengths, λ = [660, 734, 760, 808, 826, 850, 930, 980] nm and hyperspectral reconstruction using 176 wavelengths, which are equally spaced over the 650–1000 nm range. In the 8 wavelength case the 6 first wavelengths are optimally chosen according to [54

54. M. E. Eames, B. W. Pogue, and H. Dehghani, “Wavelength band optimization in spectral near-infrared optical tomography improves accuracy while reducing data acquisition and computational burden,” J. Biomed. Opt. 13, 054037 (2008). [CrossRef] [PubMed]

] with two wavelengths added where water and lipids have peak absorption. Reconstructed images created with the PaLS method are shown in Fig. 6. In simulations the SNR is set to 30 dB, as it is defined by Eq. (16) and Eq. (17). When comparing the pixel based reconstruction in Fig. 5 to the PaLS reconstruction in Fig. 6, it is evident that the PaLS method provides superior reconstructions. Examining the PaLS results, the 8 wavelength case shows reasonable accuracy along the x axis but rather diffuse results in y. We also see noticeable artifacts in the reconstructions. Considering the concentration values, the values for HbO2, HbR and water concentration come close to the actual value. Moving to hyperspectral information, the reconstruction becomes more accurate, estimating the shape close to the ground truth. It should also be noted that the runtime for each reconstruction for the PaLS method is significantly shorter compared to the pixel-based method. A PaLS reconstruction takes around 30 seconds, which is 3–4 times faster than a pixel-based method. Additionally, we do not employ any regularization parameters, freeing us from finding the optimal reconstruction using regularization. This is a major improvement in moving from a pixel-based approach to the PaLS method.

Fig. 5 Reconstruction using a pixel based method. Middle column of images are generated with 8 wavelengths and rightmost images are generated with 176 wavelengths. From top to bottom the rows show HbO2, HbR, lipid, water and diffusion amplitude, respectively. Concentration units are in mM.

The comparison of the Dice coefficient between the PaLS method and pixel-based is tricky, since for the pixel-based method the Dice coefficient is plotted as a function of a threshold. This threshold is required to create a binary map of the location on the anomaly. If the threshold is chosen to only leave extreme peak concentration values in each image, the Dice coefficient would be low due to edge artifacts as in Fig. 7(b). Therefore, in simulations we compare D(S,G) for the pixel based reconstructions using a threshold of 80% to D(S,G) of the PaLS reconstructions. The improvement of the PaLS method is confirmed quantitatively through D(S,G) and MSE displayed in Table 2 and Table 1, respectively. The Dice coefficient, shown in Table 2, gives a clear view of how the shape estimation improves by added wavelengths, where D(S,G) approaches 1 for the hyperspectral case and the PaLS method shows superior performance in the MSE values.

Fig. 7 Pixel reconstruction from both experimental sets, set 1 containing 10% ink and 90% dye and set 2 70% ink and 30% dye.

Table 1. The MSE is compared for each chromphore for multiple wavelength choices. In each case the reconstructions are done with equally spaced wavelengths over the spectrum except for the 8 wavelength case.

table-icon
View This Table
| View All Tables

Table 2. D(S,G) is compared for each chromphore for multiple wavelength choices. In each case the reconstructions are done with equally spaced wavelengths over the spectrum except for the 8 wavelength case. D(S,G) is calculated comparing 80% of the target peak to the reconstruction.

table-icon
View This Table
| View All Tables

7.2. Experimental validation

Pixel based reconstructions of absolute concentrations for both experimental sets are shown in Fig. 7 and PalS approach reconstructions in Fig. 8. As expected, the hyperspectral information provides improved reconstruction for both the pixel-based and PaLS methods. The pixel based results had previously been reported in [10

10. F. Larusson, S. Fantini, and E. L. Miller, “Hyperspectral image reconstruction for diffuse optical tomography,” Biomed. Opt. Express 2, 947–965 (2011). [CrossRef]

]. Focusing on the PaLS methods, it is evident that forming the reconstruction with shape-based constraints yields improved results. The estimation of relative concentrations and MSE of the absolute values are examined in Table 4 and Table 5 for the pixel-based and PaLS method, respectively. The relative concentration values are better estimated in both cases using the PaLS method, although the hyperspectral method does not show significant improvement for experimental set 2, which was also the case for the pixel based method. Examining the images along with the MSE values for experimental set 1, Fig. 78(a) and (c), it is noticeable how the reconstruction does not resolve the structure particurarly well along the x axis. This is somewhat unexpected since in DOT resolving depth information, on the y axis, is usually the more difficult problem. This is noticeable for both the pixel based and PaLS methods, although the PaLS method outperforms the pixel based method, especially in removing edge artifacts. This streaking in the x direction is most likely a combination of how the Gaussian basis are placed within the imaging medium, and measurement error in placing the source and detectors when taking the reference measurement.

Fig. 8 PaLS Reconstruction from both experimental sets, set 1 containing 10% ink and 90% dye and set 2 70% ink and 30% dye.

Table 3. D(S,G) is compared for each chromphore for multiple wavelength choices. In each case the reconstructions are done with equally spaced wavelengths over the spectrum except for the 6 wavelength case where we use optimally chosen wavelengths. D(S,G) is calculated comparing the half maximum of the target peak to the reconstruction.

table-icon
View This Table
| View All Tables

Table 4. Comparison of ci^ and cir^ to target concentration values for experimental results, for the pixel-based method. Best performance is highlighted in bold.

table-icon
View This Table
| View All Tables

Table 5. Comparison of ci^ and cir^ to target concentration values for experimental results, for the PaLS method. Best performance is highlighted in bold.

table-icon
View This Table
| View All Tables

For both experimental sets, the PaLS method resolves the location and the shape of the inclusion more accurately, which is verified by the calculation of the Dice coefficient shown in Table 3, the improvement is notable when compared to the pixel-based reconstruction. As discussed in Section 7.1 a choice of a threshold is needed to compare D(S,G) between pixel based reconstructions and the PaLS method. For the experimental reconstructions we use a threshold of 50% to compare D(S,G) of the PaLS reconstructions. This demonstrates the usefulness of the PaLS method correctly and accurately localizing the anomaly. The PaLS method does very well with eliminating edge artifacts that were severe when doing pixel-based reconstructions for the same data set. These effects are very noticeable in Fig. 7(b) and (d), where, especially in the multispectral case, the edge artifacts were significant. Comparing that to the same data in Fig. 8(b) and (d) it is obvious that the improvement is significant.

8. Conclusion

In this paper, using simulations and experimental measurements we have shown that the PaLS method provides more accurate estimation of chromophore concentrations than a regularized pixel-based inversion scheme. Hyperspectral information results in improved performance in terms of both MSE and spatial localization as measured using the Dice coefficient. The parametric approach is shown to give significant improvements to image reconstruction decreasing run time of the iterative process and increasing the quality of reconstructed images. The PaLS method is also easily expandable to more complicated problems where multiple geometries need to be considered.

Physical measurements were also performed to demonstrate these advantages for actual measurement data. Although exact concentration values were not achieved, there is a notable improvement associated with hyperspectral information in conjunction with the PaLS method. Additionally, improved localization of inclusions was observed for both sets when using hyper-spectral information. This emphasizes the advantage of hyperspectral information when doing reconstructions for more than one chromophore.

Based on the results in this paper, we will be extending the work to address more realistic, clinical conditions. We will first implement a fully 3D model employing boundary conditions and consider an inhomogeneous background more similar to what is found in breast imaging. In that setting, estimating multiple level sets for different geometries could prove useful, especially to estimate different regions of the heterogeneous background such as adipose and fibroglandular tissue. As mentioned in Section 3, it is also important to generate a rigid framework to choose different types, sizes and shapes of basis functions to generate the best reconstruction. To be able to estimate all shapes possible, we aim to increase the number of basis functions to include different types. To avoid over complicating the image reconstruction with a high number of basis functions we aim to pose the image reconstruction in a compressed sensing framework where few optimal basis functions estimate a complex shape. Furthermore, relating to the development of a framework for different basis functions, we will develop the Levenberg Marquardt algorithm to have optimally chosen stopping criteria and step size changes. This should increase robustness of our method and ensure accuracy of estimation for different situations and settings in diffuse optical tomography.

Acknowledgments

The authors thank Prof. Angelo Sassaroli and Alireza Aghasi for useful discussion regarding theory and Yang Yu for his help with experimental measurements. The authors would like to acknowledge financial support provided by the National Institutes of Health grant R01 CA154774.

References and links

1.

S. Fantini, E. L. Heffer, V. E. Pera, A. Sassaroli, and N. Liu, “Spatial and spectral information in optical mammography,” Technol. Cancer Res. Treat. 4, 471–482 (2005). [PubMed]

2.

R. J. Gaudette, D. H. Brooks, C. A. DiMarzio, M. E. Kilmer, E. L. Miller, T. Gaudette, and D. A. Boas, “A comparison study of linear reconstruction techniques for diffuse optical tomographic imaging of absorption coefficient,” Phys. Med. Biol 45, 1051–1069 (2000). [CrossRef] [PubMed]

3.

D. A. Boas, D. H. Brooks, E. L. Miller, C. A. DiMarzio, M. Kilmer, R. J. Gaudette, and Q. Zhang, “Imaging the body with diffuse optical tomography,” IEEE Signal Process. Mag. 18(6), 57–75 (2001). [CrossRef]

4.

M. Schweiger and S. R. Arridge, “Optical tomographic reconstruction in a complex head model using a priori region boundary information,” Phys. Med. Biol. 44, 2703–2721 (1999). [CrossRef] [PubMed]

5.

S. Fantini, D. Hueber, M. A. Franceschini, E. Gratton, W. Rosenfeld, P. G. Stubblefield, D. Maulik, and M. R. Stankovic, “Non-invasive optical monitoring of the newborn piglet brain using continuous-wave and frequency-domain spectroscopy,” Phys. Med. Biol. 44, 1543–1563 (1999). [CrossRef] [PubMed]

6.

S. Kukreti, A. E. Cerussi, W. Tanamai, D. Hsiang, B. J. Tromberg, and E. Gratton, “Characterization of metabolic differences between benign and malignant tumors: high-spectral-resolution diffuse optical spectroscopy,” Radiology 254, 277–284 (2010). [CrossRef]

7.

A. Li, G. Boverman, Y. Zhang, D. Brooks, E. L. Miller, M. E. Kilmer, Q. Zhang, E. M. C. Hillman, and D. A. Boas, “Optimal linear inverse solution with multiple priors in diffuse optical tomography,” Appl. Opt. 44, 1948–1956 (2005). [CrossRef] [PubMed]

8.

H. Dehghani, B. W. Pogue, S. P. Poplack, and K. D. Paulsen, “Multiwavelength three-dimensional near-infrared tomography of the breast: initial simulation, phantom, and clinical results,” Appl. Opt. 42, 135–145 (2004). [CrossRef]

9.

J. P. Culver, V. Ntziachristos, M. J. Holboke, and A. G. Yodh, “Optimization of optode arrangements for diffuse optical tomography: a singular-value analysis,” Opt. Lett. 26, 701–703 (2001). [CrossRef]

10.

F. Larusson, S. Fantini, and E. L. Miller, “Hyperspectral image reconstruction for diffuse optical tomography,” Biomed. Opt. Express 2, 947–965 (2011). [CrossRef]

11.

B. W. Pogue, T. O. McBride, J. Prewitt, U. L. Osterberg, and K. D. Paulsen, “Spatially variant regularization improves diffuse optical tomography,” Appl. Opt. 38, 2950–2961 (1999). [CrossRef]

12.

K. D. Paulsen and H. Jiang, “Enhanced frequency-domain optical image reconstruction in tissues through total-variation minimization,” Appl. Opt. 35, 3447–3458, (1996). [CrossRef] [PubMed]

13.

G. Boverman, Q. Fang, S. A. Carp, E. L. Miller, D. H. Brooks, J. Selb, R. H. Moore, D. B. Kopans, and D. A. Boas, “Spatio-temporal imaging of the hemoglobin in the compressed breast with diffuse optical tomography,” Phys. Med. Biol. 52, 3619–3641 (2007). [CrossRef] [PubMed]

14.

D. A. Boas, “A fundamental limitation of linearized algorithms for diffuse optical tomography,” Opt. Express 1, 404–413 (1997). [CrossRef] [PubMed]

15.

A. Aghasi, M. Kilmer, and E. L. Miller “Parametric level set methods for inverse problems,” SIAM J. Imaging Sci. 4, 618–650 (2011). [CrossRef]

16.

A. Aghasi, E. L. Miller, and L. M. Abriola “Characterization of source zone architecture: a joint electrical and hydrological inversion approach,” presented at 2011 Fall Meeting, AGU, San Francisco, Calif., 5–9 Dec. 2011.

17.

F. Larusson, S. Fantini, and E. L. Miller, “Parametric level-set approach for hyperspectral diffuse optical tomography,” in 2011 IEEE International Symposium on Biomedical Imaging: from Nano to Macro (IEEE, 2011), pp. 949–955.

18.

O. Dorn, “A shape reconstruction method for diffuse optical tomography using a transport model and level sets,” in 2002 IEEE International Symposium on Biomedical Imaging, 2002. Proceedings (IEEE, 2002), pp. 1015–1018.

19.

M. Schweiger, O. Dorn, and S. R. Arridge, “3-D shape and contrast reconstruction in optical tomography with level sets,” J. Phys.: Conf. Ser. 124, 012043 (2008). [CrossRef]

20.

M. E. Kilmer, E. L. Miller, A. Barbaro, and David Boas, “Three-dimensional shape-based imaging of absorption perturbation for diffuse optical tomography,” Appl. Opt. 42, 3129–3144 (2003). [CrossRef] [PubMed]

21.

S. R. Arridge, O. Dorn, V. Kolehmainen, M. Schweiger, and A. Zacharopoulos, “Parameter and structure reconstruction in optical tomography,” J. Phys.: Conf. Ser. 135, 012001 (2008). [CrossRef]

22.

O. Dorn and D. Lesselier, “Level set methods for inverse scattering,” Inverse Probl. 22, R67–R131 (2006). [CrossRef]

23.

M. A. O’Leary, D. A. Boas, B. Chance, and A. G. Yodh, “Experimental images of heterogeneous turbid media by frequency-domain diffusing-photon tomography,” Opt. Lett. 20, 426–428 (1995). [CrossRef]

24.

B. Brooksby, B. W. Pogue, S. Jiang, H. Dehghani, S. Srinivasan, C. Kogel, T. D. Tosteson, J. Weaver, S. P. Poplack, and K. D. Paulsen, “Imaging breast adipose and fibroglandular tissue molecular signatures by using hybrid mri-guided near-infrared spectral tomography,” Proc. Natl. Acad. Sci. U.S.A. 103, 8828–8833 (2006). [CrossRef] [PubMed]

25.

D. Boas, M. A. O’Leary, B. Chance, and A. Yodh, “Scattering of diffuse photon density waves by spherical inhomogeneities within turbid media: analytic solution and applications,” Proc. Natl. Acad. Sci. U.S.A. 91, 4887–4891 (1994). [CrossRef] [PubMed]

26.

A. Mandelis, Diffusion-Wave Fields: Mathematical Methods and Green Functions (Springer, 2001).

27.

B. Brendel, R. Ziegler, and T. Nielsen “Algebraic reconstruction techniques for spectral reconstruction in diffuse optical tomography,” Appl. Opt. 47, 6392–6403 (2008). [CrossRef] [PubMed]

28.

A. Corlu, R. Choe, T. Durduran, K. Lee, M. Schweiger, S. R. Arridge, E. M. Hillman, and A. G. Yodh, “Diffuse optical tomography with spectral constraints and wavelength optimization,” Appl. Opt. 44, 2082–2093 (2005). [CrossRef] [PubMed]

29.

A. Cerussi, D. Hsiang, N. Shah, R. Mehta, A. Durkin, J. Butler, and B. J. Tromberg “Predicting response to breast cancer neoadjuvant chemotherapy using diffuse optical spectroscopy,” Proc. Nat. Acad. Sci. U.S.A. 104, 4014–4019 (2007). [CrossRef]

30.

H. Soliman, A. Gunasekara, M. Rycroft, J. Zubovits, R. Dent, J. Spayne, M. J. Yaffe, and G. Czarnota, “Functional imaging using diffuse optical spectroscopy of neoadjuvant chemotherapy response in women with locally advanced breast cancer,” Clin. Cancer Res. 16, 2605–2614 (2010). [CrossRef] [PubMed]

31.

Q. Zhu, P. U. Hegde, A. Ricci, M. Kane, E. B. Cronin, Y. Ardeshirpour, C. Xu, A. Aguirre, S. H. Kurtzman, P. J. Deckers, and S. H. Tannenbaum, “Early-stage invasive breast cancers: potential role of optical tomography with US localization in assisting diagnosis,” Radiology 256, 367–378 (2010). [CrossRef] [PubMed]

32.

T. Chan and L. Vese “Active contours without edges,” IEEE Trans. Image Process. 10, 266–277 (2001). [CrossRef]

33.

H. K. Zhao, S. Osher, B. Merriman, and M. Kang, “Implicit, nonparametric shape reconstruction from unorganized points using a variational level set method,” Comput. Vision Image Understanding 80, 295–314 (2000). [CrossRef]

34.

S. Osher and R. Fedkiw, Level Set Methods and Dynamic Implicit Surfaces (Springer, 2002)

35.

M. Guven, B. Yazici, X. Intes, and B. Chance, “Diffuse optical tomography with a priori anatomical information,” Phys. Med. Biol. 50, 2837–2858 (2005). [CrossRef] [PubMed]

36.

K. Madsen, H. Bruun, and O. Tingleff “Methods for non-linear least squares problems,” lecture notes (2004).

37.

C. R. Vogel, Computational Methods for Inverse Problems (SIAM, 2002) [CrossRef]

38.

T. T. Wu and K. Lange, “Coordinate descent algorithms for lasso penalized regression,” Ann. Appl. Stat. 2, 224–244 (2008). [CrossRef]

39.

B. Brendel and T. Nielsen, “Selection of optimal wavelengths for spectral reconstruction in diffuse optical tomography,” J. Biomed. Opt. 14, 034041 (2009). [CrossRef] [PubMed]

40.

S. Srinivasan, B. W. Pogue, S. Jiang, H. Dehghani, C. Kogel, S. Soho, J. J. Gibson, T. D. Tosteson, S. P. Poplack, and K. D. Paulsen, “Interpreting hemoglobin and water concentration, oxygen saturation, and scattering measured in vivo by near-infrared breast tomography,” Proc. Natl. Acad. Sci. U.S.A. 100, 12349–12354 (2003). [CrossRef] [PubMed]

41.

P. Taroni, A. Pifferi, A. Torricelli, D. Comelli, and R. Cubeddu, “In vivo absorption and scattering spectroscopy of biological tissues,” Photochem. Photobiol. Sci. 2, 124–129 (2003). [CrossRef] [PubMed]

42.

S. Prahl, “Tabulated molar extinction coefficient for hemoglobin in water” (Oregon Medical Laser Center, 2007), http://omlc.ogi.edu/spectra/hemoglobin/summary.html.

43.

A. A. Joshi, A. J. Chaudhari, D. W. Shattuck, J. Dutta, R. M. Leahy, and A. W. Toga, “Posture matching and elastic registration of a mouse atlas to surface topography range data,” in IEEE International Symposium on Biomedical Imaging: from Nano to Macro, 2009. ISBI ’09 (IEEE, 2009), pp. 366–369 (2009).

44.

J. D. Vylder and W. Philips, “A computational efficient external energy for active contour segmentation using edge propagation,” in IEEE 2100 International Conference on Image Processing (ICIP 2010) (IEEE, 2010), pp. 661–664.

45.

T. Durduran, R. Choe, J. P. culver, L. Zubkov, M. J. Holboke, J. Giammarco, B. Chance, and A. G. Yodh, “Bulk optical properties of healthy female breast tissue,” Phys. Med. Biol. 47, 2847–2861 (2002). [CrossRef] [PubMed]

46.

L. Spinelli, A. Torricelli, A. Pifferi, P. Taroni, G. M. Danesini, and R. Cubeddu, “Bulk optical properties and tissue components in the female breast from multiwavelength time-resolved optical mammography,” J. Biomed. Opt. 9, 1137–1142 (2004). [CrossRef] [PubMed]

47.

J. P. Culver, R. Choe, M. J. Holboke, L. Zubkov, T. Durduran, A. Slemp, V. Ntziachristos, B. Chance, and A. G. Yodh, “Three-dimensional diffuse optical tomography in the parallel plane transmission geometry: evaluation of a hybrid frequency domain/continuous wave clinical system for breast imaging,” Med. Phys. 30, 235–247 (2003). [CrossRef] [PubMed]

48.

N. Liu, A. Sassaroli, and S. Fantini, “Paired-wavelength spectral approach to measuring the relative concentrations of two localized chromophores in turbid media: an experimental study,” J. Biomed. Opt. 12, 051602 (2007). [CrossRef] [PubMed]

49.

E. Gratton, S. Fantini, M. A. Franceschini, C. Gratton, and M. Fabiani, “Measurements of scattering and absorption changes in muscle and brain,” Philos. Trans. R. Soc. Lond. B. Biol. Sci. 352, 727–735 (1997). [CrossRef] [PubMed]

50.

B. W. Pogue, S. Jiang, H. Dehghani, C. Kogel, S. Soho, S. Srinivasan, X. Song, T. D. Tosteson, S. P. Poplack, and K. D. Paulsen, “Characterization of hemoglobin, water, and nir scattering in breast tissue: analysis of intersubject variability and menstrual cycles changes,” J. Biomed. Opt. 9, 541–552 (2004). [CrossRef] [PubMed]

51.

G. Boverman, E. L. Miller, A. Li, Q. Zhang, T. Chaves, D. H. Brooks, and D. A. Boas, “Quantitative spectroscopic diffuse optical tomography of the breast guided by imperfect a priori structural information,” Phys. Med. Biol. 50, 3941–3956 (2005). [CrossRef] [PubMed]

52.

S. Fantini, M. A. Franceschini, J. S. Maier, S. A. Walker, B. Barbieri, and E. Gratton, “Frequency-domain multi-channel optical detector for noninvasive tissue spectroscopy and oximetry,” Opt. Eng. 34, 32–42 (1995). [CrossRef]

53.

M. A. Franceschini, V. Toronov, M. E. Filiaci, E. Gratton, and S. Fantini, “On-line optical imaging of the human brain with 160-ms temporal resolution,” Opt. Express 6, 49–57 (2000). [CrossRef] [PubMed]

54.

M. E. Eames, B. W. Pogue, and H. Dehghani, “Wavelength band optimization in spectral near-infrared optical tomography improves accuracy while reducing data acquisition and computational burden,” J. Biomed. Opt. 13, 054037 (2008). [CrossRef] [PubMed]

OCIS Codes
(100.3190) Image processing : Inverse problems
(170.3010) Medical optics and biotechnology : Image reconstruction techniques
(170.3660) Medical optics and biotechnology : Light propagation in tissues
(170.3830) Medical optics and biotechnology : Mammography
(170.3880) Medical optics and biotechnology : Medical and biological imaging
(170.5280) Medical optics and biotechnology : Photon migration
(170.6960) Medical optics and biotechnology : Tomography
(290.1990) Scattering : Diffusion
(290.7050) Scattering : Turbid media

ToC Category:
Image Reconstruction and Inverse Problems

History
Original Manuscript: January 27, 2012
Manuscript Accepted: March 15, 2012
Published: April 18, 2012

Citation
Fridrik Larusson, Sergio Fantini, and Eric L. Miller, "Parametric level set reconstruction methods for hyperspectral diffuse optical tomography," Biomed. Opt. Express 3, 1006-1024 (2012)
http://www.opticsinfobase.org/boe/abstract.cfm?URI=boe-3-5-1006


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. S. Fantini, E. L. Heffer, V. E. Pera, A. Sassaroli, and N. Liu, “Spatial and spectral information in optical mammography,” Technol. Cancer Res. Treat.4, 471–482 (2005). [PubMed]
  2. R. J. Gaudette, D. H. Brooks, C. A. DiMarzio, M. E. Kilmer, E. L. Miller, T. Gaudette, and D. A. Boas, “A comparison study of linear reconstruction techniques for diffuse optical tomographic imaging of absorption coefficient,” Phys. Med. Biol45, 1051–1069 (2000). [CrossRef] [PubMed]
  3. D. A. Boas, D. H. Brooks, E. L. Miller, C. A. DiMarzio, M. Kilmer, R. J. Gaudette, and Q. Zhang, “Imaging the body with diffuse optical tomography,” IEEE Signal Process. Mag.18(6), 57–75 (2001). [CrossRef]
  4. M. Schweiger and S. R. Arridge, “Optical tomographic reconstruction in a complex head model using a priori region boundary information,” Phys. Med. Biol.44, 2703–2721 (1999). [CrossRef] [PubMed]
  5. S. Fantini, D. Hueber, M. A. Franceschini, E. Gratton, W. Rosenfeld, P. G. Stubblefield, D. Maulik, and M. R. Stankovic, “Non-invasive optical monitoring of the newborn piglet brain using continuous-wave and frequency-domain spectroscopy,” Phys. Med. Biol.44, 1543–1563 (1999). [CrossRef] [PubMed]
  6. S. Kukreti, A. E. Cerussi, W. Tanamai, D. Hsiang, B. J. Tromberg, and E. Gratton, “Characterization of metabolic differences between benign and malignant tumors: high-spectral-resolution diffuse optical spectroscopy,” Radiology254, 277–284 (2010). [CrossRef]
  7. A. Li, G. Boverman, Y. Zhang, D. Brooks, E. L. Miller, M. E. Kilmer, Q. Zhang, E. M. C. Hillman, and D. A. Boas, “Optimal linear inverse solution with multiple priors in diffuse optical tomography,” Appl. Opt.44, 1948–1956 (2005). [CrossRef] [PubMed]
  8. H. Dehghani, B. W. Pogue, S. P. Poplack, and K. D. Paulsen, “Multiwavelength three-dimensional near-infrared tomography of the breast: initial simulation, phantom, and clinical results,” Appl. Opt.42, 135–145 (2004). [CrossRef]
  9. J. P. Culver, V. Ntziachristos, M. J. Holboke, and A. G. Yodh, “Optimization of optode arrangements for diffuse optical tomography: a singular-value analysis,” Opt. Lett.26, 701–703 (2001). [CrossRef]
  10. F. Larusson, S. Fantini, and E. L. Miller, “Hyperspectral image reconstruction for diffuse optical tomography,” Biomed. Opt. Express2, 947–965 (2011). [CrossRef]
  11. B. W. Pogue, T. O. McBride, J. Prewitt, U. L. Osterberg, and K. D. Paulsen, “Spatially variant regularization improves diffuse optical tomography,” Appl. Opt.38, 2950–2961 (1999). [CrossRef]
  12. K. D. Paulsen and H. Jiang, “Enhanced frequency-domain optical image reconstruction in tissues through total-variation minimization,” Appl. Opt.35, 3447–3458, (1996). [CrossRef] [PubMed]
  13. G. Boverman, Q. Fang, S. A. Carp, E. L. Miller, D. H. Brooks, J. Selb, R. H. Moore, D. B. Kopans, and D. A. Boas, “Spatio-temporal imaging of the hemoglobin in the compressed breast with diffuse optical tomography,” Phys. Med. Biol.52, 3619–3641 (2007). [CrossRef] [PubMed]
  14. D. A. Boas, “A fundamental limitation of linearized algorithms for diffuse optical tomography,” Opt. Express1, 404–413 (1997). [CrossRef] [PubMed]
  15. A. Aghasi, M. Kilmer, and E. L. Miller “Parametric level set methods for inverse problems,” SIAM J. Imaging Sci.4, 618–650 (2011). [CrossRef]
  16. A. Aghasi, E. L. Miller, and L. M. Abriola “Characterization of source zone architecture: a joint electrical and hydrological inversion approach,” presented at 2011 Fall Meeting, AGU, San Francisco, Calif., 5–9 Dec. 2011.
  17. F. Larusson, S. Fantini, and E. L. Miller, “Parametric level-set approach for hyperspectral diffuse optical tomography,” in 2011 IEEE International Symposium on Biomedical Imaging: from Nano to Macro (IEEE, 2011), pp. 949–955.
  18. O. Dorn, “A shape reconstruction method for diffuse optical tomography using a transport model and level sets,” in 2002 IEEE International Symposium on Biomedical Imaging, 2002. Proceedings (IEEE, 2002), pp. 1015–1018.
  19. M. Schweiger, O. Dorn, and S. R. Arridge, “3-D shape and contrast reconstruction in optical tomography with level sets,” J. Phys.: Conf. Ser.124, 012043 (2008). [CrossRef]
  20. M. E. Kilmer, E. L. Miller, A. Barbaro, and David Boas, “Three-dimensional shape-based imaging of absorption perturbation for diffuse optical tomography,” Appl. Opt.42, 3129–3144 (2003). [CrossRef] [PubMed]
  21. S. R. Arridge, O. Dorn, V. Kolehmainen, M. Schweiger, and A. Zacharopoulos, “Parameter and structure reconstruction in optical tomography,” J. Phys.: Conf. Ser.135, 012001 (2008). [CrossRef]
  22. O. Dorn and D. Lesselier, “Level set methods for inverse scattering,” Inverse Probl.22, R67–R131 (2006). [CrossRef]
  23. M. A. O’Leary, D. A. Boas, B. Chance, and A. G. Yodh, “Experimental images of heterogeneous turbid media by frequency-domain diffusing-photon tomography,” Opt. Lett.20, 426–428 (1995). [CrossRef]
  24. B. Brooksby, B. W. Pogue, S. Jiang, H. Dehghani, S. Srinivasan, C. Kogel, T. D. Tosteson, J. Weaver, S. P. Poplack, and K. D. Paulsen, “Imaging breast adipose and fibroglandular tissue molecular signatures by using hybrid mri-guided near-infrared spectral tomography,” Proc. Natl. Acad. Sci. U.S.A.103, 8828–8833 (2006). [CrossRef] [PubMed]
  25. D. Boas, M. A. O’Leary, B. Chance, and A. Yodh, “Scattering of diffuse photon density waves by spherical inhomogeneities within turbid media: analytic solution and applications,” Proc. Natl. Acad. Sci. U.S.A.91, 4887–4891 (1994). [CrossRef] [PubMed]
  26. A. Mandelis, Diffusion-Wave Fields: Mathematical Methods and Green Functions (Springer, 2001).
  27. B. Brendel, R. Ziegler, and T. Nielsen “Algebraic reconstruction techniques for spectral reconstruction in diffuse optical tomography,” Appl. Opt.47, 6392–6403 (2008). [CrossRef] [PubMed]
  28. A. Corlu, R. Choe, T. Durduran, K. Lee, M. Schweiger, S. R. Arridge, E. M. Hillman, and A. G. Yodh, “Diffuse optical tomography with spectral constraints and wavelength optimization,” Appl. Opt.44, 2082–2093 (2005). [CrossRef] [PubMed]
  29. A. Cerussi, D. Hsiang, N. Shah, R. Mehta, A. Durkin, J. Butler, and B. J. Tromberg “Predicting response to breast cancer neoadjuvant chemotherapy using diffuse optical spectroscopy,” Proc. Nat. Acad. Sci. U.S.A.104, 4014–4019 (2007). [CrossRef]
  30. H. Soliman, A. Gunasekara, M. Rycroft, J. Zubovits, R. Dent, J. Spayne, M. J. Yaffe, and G. Czarnota, “Functional imaging using diffuse optical spectroscopy of neoadjuvant chemotherapy response in women with locally advanced breast cancer,” Clin. Cancer Res.16, 2605–2614 (2010). [CrossRef] [PubMed]
  31. Q. Zhu, P. U. Hegde, A. Ricci, M. Kane, E. B. Cronin, Y. Ardeshirpour, C. Xu, A. Aguirre, S. H. Kurtzman, P. J. Deckers, and S. H. Tannenbaum, “Early-stage invasive breast cancers: potential role of optical tomography with US localization in assisting diagnosis,” Radiology256, 367–378 (2010). [CrossRef] [PubMed]
  32. T. Chan and L. Vese “Active contours without edges,” IEEE Trans. Image Process.10, 266–277 (2001). [CrossRef]
  33. H. K. Zhao, S. Osher, B. Merriman, and M. Kang, “Implicit, nonparametric shape reconstruction from unorganized points using a variational level set method,” Comput. Vision Image Understanding80, 295–314 (2000). [CrossRef]
  34. S. Osher and R. Fedkiw, Level Set Methods and Dynamic Implicit Surfaces (Springer, 2002)
  35. M. Guven, B. Yazici, X. Intes, and B. Chance, “Diffuse optical tomography with a priori anatomical information,” Phys. Med. Biol.50, 2837–2858 (2005). [CrossRef] [PubMed]
  36. K. Madsen, H. Bruun, and O. Tingleff “Methods for non-linear least squares problems,” lecture notes (2004).
  37. C. R. Vogel, Computational Methods for Inverse Problems (SIAM, 2002) [CrossRef]
  38. T. T. Wu and K. Lange, “Coordinate descent algorithms for lasso penalized regression,” Ann. Appl. Stat.2, 224–244 (2008). [CrossRef]
  39. B. Brendel and T. Nielsen, “Selection of optimal wavelengths for spectral reconstruction in diffuse optical tomography,” J. Biomed. Opt.14, 034041 (2009). [CrossRef] [PubMed]
  40. S. Srinivasan, B. W. Pogue, S. Jiang, H. Dehghani, C. Kogel, S. Soho, J. J. Gibson, T. D. Tosteson, S. P. Poplack, and K. D. Paulsen, “Interpreting hemoglobin and water concentration, oxygen saturation, and scattering measured in vivo by near-infrared breast tomography,” Proc. Natl. Acad. Sci. U.S.A.100, 12349–12354 (2003). [CrossRef] [PubMed]
  41. P. Taroni, A. Pifferi, A. Torricelli, D. Comelli, and R. Cubeddu, “In vivo absorption and scattering spectroscopy of biological tissues,” Photochem. Photobiol. Sci.2, 124–129 (2003). [CrossRef] [PubMed]
  42. S. Prahl, “Tabulated molar extinction coefficient for hemoglobin in water” (Oregon Medical Laser Center, 2007), http://omlc.ogi.edu/spectra/hemoglobin/summary.html .
  43. A. A. Joshi, A. J. Chaudhari, D. W. Shattuck, J. Dutta, R. M. Leahy, and A. W. Toga, “Posture matching and elastic registration of a mouse atlas to surface topography range data,” in IEEE International Symposium on Biomedical Imaging: from Nano to Macro, 2009. ISBI ’09 (IEEE, 2009), pp. 366–369 (2009).
  44. J. D. Vylder and W. Philips, “A computational efficient external energy for active contour segmentation using edge propagation,” in IEEE 2100 International Conference on Image Processing (ICIP 2010) (IEEE, 2010), pp. 661–664.
  45. T. Durduran, R. Choe, J. P. culver, L. Zubkov, M. J. Holboke, J. Giammarco, B. Chance, and A. G. Yodh, “Bulk optical properties of healthy female breast tissue,” Phys. Med. Biol.47, 2847–2861 (2002). [CrossRef] [PubMed]
  46. L. Spinelli, A. Torricelli, A. Pifferi, P. Taroni, G. M. Danesini, and R. Cubeddu, “Bulk optical properties and tissue components in the female breast from multiwavelength time-resolved optical mammography,” J. Biomed. Opt.9, 1137–1142 (2004). [CrossRef] [PubMed]
  47. J. P. Culver, R. Choe, M. J. Holboke, L. Zubkov, T. Durduran, A. Slemp, V. Ntziachristos, B. Chance, and A. G. Yodh, “Three-dimensional diffuse optical tomography in the parallel plane transmission geometry: evaluation of a hybrid frequency domain/continuous wave clinical system for breast imaging,” Med. Phys.30, 235–247 (2003). [CrossRef] [PubMed]
  48. N. Liu, A. Sassaroli, and S. Fantini, “Paired-wavelength spectral approach to measuring the relative concentrations of two localized chromophores in turbid media: an experimental study,” J. Biomed. Opt.12, 051602 (2007). [CrossRef] [PubMed]
  49. E. Gratton, S. Fantini, M. A. Franceschini, C. Gratton, and M. Fabiani, “Measurements of scattering and absorption changes in muscle and brain,” Philos. Trans. R. Soc. Lond. B. Biol. Sci.352, 727–735 (1997). [CrossRef] [PubMed]
  50. B. W. Pogue, S. Jiang, H. Dehghani, C. Kogel, S. Soho, S. Srinivasan, X. Song, T. D. Tosteson, S. P. Poplack, and K. D. Paulsen, “Characterization of hemoglobin, water, and nir scattering in breast tissue: analysis of intersubject variability and menstrual cycles changes,” J. Biomed. Opt.9, 541–552 (2004). [CrossRef] [PubMed]
  51. G. Boverman, E. L. Miller, A. Li, Q. Zhang, T. Chaves, D. H. Brooks, and D. A. Boas, “Quantitative spectroscopic diffuse optical tomography of the breast guided by imperfect a priori structural information,” Phys. Med. Biol.50, 3941–3956 (2005). [CrossRef] [PubMed]
  52. S. Fantini, M. A. Franceschini, J. S. Maier, S. A. Walker, B. Barbieri, and E. Gratton, “Frequency-domain multi-channel optical detector for noninvasive tissue spectroscopy and oximetry,” Opt. Eng.34, 32–42 (1995). [CrossRef]
  53. M. A. Franceschini, V. Toronov, M. E. Filiaci, E. Gratton, and S. Fantini, “On-line optical imaging of the human brain with 160-ms temporal resolution,” Opt. Express6, 49–57 (2000). [CrossRef] [PubMed]
  54. M. E. Eames, B. W. Pogue, and H. Dehghani, “Wavelength band optimization in spectral near-infrared optical tomography improves accuracy while reducing data acquisition and computational burden,” J. Biomed. Opt.13, 054037 (2008). [CrossRef] [PubMed]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited