Abstract

We present a wide-field method for obtaining three-dimensional images of turbid media. By projecting patterns of light of varying spatial frequencies on a sample, we reconstruct quantitative, depth resolved images of absorption contrast. Images are reconstructed using a fast analytic inversion formula and a novel correction to the diffusion approximation for increased accuracy near boundaries. The method provides more accurate quantification of optical absorption and higher resolution than standard diffuse reflectance measurements.

© 2009 OSA

1. Introduction

Optical tomography with multiply scattered light [1] is a rapidly growing field with applications ranging from human breast [2] and brain [3], to small animals [4]. Recently, in an effort to image large fields of view quickly, we and others have begun projecting spatially extended patterns of light on tissue, eliminating the need to use large arrays of optical fibers or raster scan a collimated beam [57]. In particular, it has been experimentally demonstrated that by projecting sinusoidal patterns of light onto tissue, one can determine the tissue’s optical properties by measuring the relative decay of spatial patterns of differing frequencies [5,8]. This technique has been applied to imaging stroke in the rat cerebral cortex [9], visualizing cortical spreading depression [10], and assessing layered structures in skin [11]. The principle advantages of projecting a spatially extended source are that one can image a large field of view quickly, it is non-contact, inexpensive, and it eliminates the need for detectors with a wide dynamic range since all detectors are close to a spatially extended source. This is especially important for reflectance measurements, since the remitted light intensity due to a point source decays by many orders of magnitude within millimeters.

In our previous work using a spatially modulated light source, tissue optical properties were determined by fitting to a model of light propagation that assumed the tissue was a homogeneous or a layered structure. The data acquired by each detector (i.e. each pixel on the CCD array) was analyzed independently of all the other detectors. This approach tacitly assumed that tissue optical properties (within layers) did not vary with depth, and that neighboring volumes of tissue did not affect the light intensity measured at a given pixel. This approach is sufficient for samples in which optical property changes are both gradual and located near the surface. However, we expect images of small buried objects to appear blurred, and to underestimate optical property changes. This is because each detector samples a volume much larger than the absorber, and because each detector is primarily sensitive to the optical properties at the surface. Thus, the optical properties assigned to any region of tissue are in fact the average values for a much larger region, and this average is highly weighted by the optical properties at the surface. The goal of image reconstruction is ameliorate this partial volume effect, by modeling the sensitivity of each measurement to the entire tissue volume.

The theoretical framework for describing the propagation of plane waves of diffuse light and reconstructing images of sub-surface heterogeneities in the spatial Fourier domain have been known for several years [12,13]. However, the original reconstructed images using these methods were acquired by raster scanning a collimated laser beam and numerically taking a Fourier transform of the resulting data [14,15]. Recently, an algorithm was proposed for reconstructing three-dimensional images directly from measurements made by illuminating tissue with sinusoidal patterns [16]. The algorithm takes advantage of the insight that illuminating with sinusoidal functions of different frequencies is theoretically equivalent to raster scanning a point source and then taking a Fourier transform of the data with respect to the source positions.

In this manuscript we experimentally validate the tomographic approach of Ref [16]. using tissue simulating phantoms with buried heterogeneities resembling superficial veins. We also modify this approach by implementing a novel correction to the diffusion approximation which gives increased accuracy in the reconstructed images near the tissue surface [17]. We quantitatively measure the optical absorption of buried objects, determine their relative depths, and image them with higher resolution than is obtained by fitting to a homogeneous model.

2. Methods

2.1 Theory

The propagation of multiply scattered photons in biological tissue is governed by the radiative transport equation (RTE). For a continuous-wave light source at a single wavelength the RTE for light appears as:

s^I(r,s^)+(μa+μs)I(r,s^)μsd2s'A(s^,s^')I(r,s^')=S(r,s^)  .
Here I(r,s^) is the specific intensity, μa is the absorption coefficient, μs is the scattering coefficient, A(s^,s^') is the scattering kernel, and S(r,s^) is the source. The diffusion approximation to the RTE is obtained by expanding I(r,s^) to first order in s^ resulting in an approximation for I(r,s^)that consists of both a photon density and a current term:

I(r,s^)=c4πΦ(r)+34πJs^  .

Near boundaries, where the specific intensity is not isotropic, both the density and current terms are important. Retaining both terms, the Green’s function to the RTE can be approximated in terms of the Greens function for the diffusion equation as [12]:

G(r,s^,r',s^')=c4π(1+s^r)(1s^'r')G(r,r') ,
where c is the speed of light in the medium, =1/(μa+μs'), and G(r,r') is the Greens function for the diffusion equation which obeys the equation
[D+ cμa] G(r,r') = δ(rr') ,
subject to the boundary condition
G(r,r')+n^G(r,r')=0  .
Here D=c/3 is the diffusion coefficient, is the extrapolation length, and n^ is the outward pointing unit normal. As demonstrated in Ref [18], the Greens function for the diffusion equation in a semi-infinite medium can be decomposed into plane waves according to
G(r,r')=1(2π)2d2q g(q,z,z')exp[i(q(ρ'ρ))] ,
where the Fourier component
g(q,z,z')=Dexp[Q(q)|zz'|]Q(q)+1
governs the decay of spatial waves as they propagate to or from the surface of the tissue. The decay rate Q(q)=(cμa/D+q2)1/2 is determined by the optical properties of the medium and the frequency of the spatial wave.

Our goal is to determine the spatial variations in absorption by projecting sinusoidal patterns of differing frequencies and phases on the sample, and measuring the spatially varying intensity of the remitted light at the surface. The change in specific intensity due to a perturbation of the absorption coefficient in the medium obeys to first order the following integral equation:

φ(k,rd,s^d)=d3r d2s I(r,s^) G(r,s^,rd,s^d) δμa(r) ,
where φ(k,rd,s^d) is the experimentally measured data function at detector location rd, in direction s^d, due to a light source with spatial wave number k. δμa(r) is the change in absorption, and I(r,s^) is the specific intensity do to the spatially extended source and obeys the equation
I(r,s^)=d3r'd2s' G(r',s^',r,s^)S(r',s^')  .
In our experiment, the source consists of light which is sinusoidal in space with wave number k, has intensity I0, and is directed into the sample at the surface. It is described as:
S(r,s^)=I0exp(ikρ)δ(z)δ(s^+n^)  .
By substituting Eqs. (3) and (9) into Eq. (8) and using Eqs. (5-7) and (10), Eq. (8) can be put in the form
φ(k,ρd)=c I04π(1+)2d3r K(k,ρd;r) c δμa(r) ,
where
K(k,ρd;r)=1(2π)2d2q κ(k,q;z) exp[i(kq)ρ+iqρd] ,
and

κ(k,q;z)=(D)2[123(kq+Q(k)Q(q))]exp[(Q(k)+Q(q)) z](Q(k)+1)(Q(q)+1)  .

By taking a 2D Fourier transform with respect to the detector locations ρd, Eq. (11) can be put in block diagonal form and inverted by analytical methods which have been previously described [12,19].

Note, the source function specified by Eq. (10) contains a complex exponential, while the source intensity in the actual experiment is real and positive. It has the form

S(k,ρ,θ)=I02(1+Mcos(kρ+θ)) δ(z) δ(s^+n^)  .

Here M and θ are the modulation depth and phase of the source. Since the diffusion equation is linear, we can obtain the data that would be generated by a source of the form of Eq. (10) by using the source described by Eq. (14). To do so, we acquire data for at least three phases, and add a linear combination of the resulting measurements of light intensity I(ρ) such that the total source distribution from the three measurements equals Eq. (10). For example, using three phases (θ=2π/3,0,2π/3), we combine the resulting measurements as

Itot=13(2I2I1I3)+i3(I1I3)  .

Due to the relatively large subsurface absorption perturbation we use the Rytov approximation [20]. In order to calculate the Rytov data function we use the standard rules for complex numbers to write Itot in terms of an amplitude and phase (i.e. Itot=Aeiϕ), where for three phases

A=23{(I1I2)2+(I2I3)2+(I3I1)2}1/2 ,
and
ϕ=tan1{3(I1I3)/(2I2I1I3)}  .
Then the Rytov data function is calculated as:
φRytov=Itot(0){log(A/A(0))+i(ϕϕ(0))}  .
This formulation differs from the usual temporal frequency domain in that the amplitude and phase represent the Fourier components in space, and not in time. The data acquired this way is equivalent to the data acquired by using point sources and then taking a Fourier transform with respect to their spatial positions.

2.2 Instrumentation

A schematic of the instrument is shown in Fig. 1 . Broadband light was generated with a Newport Corporation power source (Newport Corporation, Irvine CA) and a 250 watt tungsten lamp. This light was used to illuminate a DLP Developers Kit 1024 x 768 pixel digital micro-mirror device (Texas Instrument, Dallas TX) which spatially modulated the light to be projected onto the phantom. Frequency patterns for the spatial modulation were generated using customized C# code (Modulated Imaging Inc., Irvine CA) and ALP (Application Program Interface) software. The patterned light was projected directly onto the phantoms approximately 500 mm from the projector system. Eleven frequencies evenly spaced from 0 to 1 cm−1 with four phases each were used. The diffuse reflected light images were acquired at 650 nm with a Nuance Multispectral Imaging System (CRi Inc., Woburn, MA) which consists of a liquid crystal tunable filter capable of acquiring at discrete wavelengths (10nm bandwidth) between 650 nm and 1100 nm, and a 1040 x 1392 pixel, front illuminated charge coupled device (CCD). A 6x6 cm field of view was acquired by the CCD using 4x4 hardware binning and an exposure time of ~100 ms. Cross-polarizers were placed at the source and detector in order to eliminate any specular reflection.

 

Fig. 1 Schematic of the experimental setup. Broadband light from a lamp is projected on the sample using a digital micro-mirror device (DMD) and detected at discrete wavelengths using a liquid crystal tunable filter and CCD. Crossed polarizers eliminate specular reflections.

Download Full Size | PPT Slide | PDF

2.3 Phantoms

Tissue simulating phantoms were made from silicone (P4, Eager Plastics, Chicago, IL), with titanium dioxide (TI-602, Atlantic Equipment Engineers, Bergenfield, NJ) used as a scattering agent, and nigrosin (Sigma, St. Louis, MO) as an absorber [21]. Optical property measurements obtained by fitting the diffuse reflectance at the different measured spatial frequencies to an analytic model for a homogeneous semi-infinite medium [5] were μa=0.004 mm-1 and μs'=0.9 mm-1 for 650 nm light. These values were calculated using standard reflectance phantoms measured previously using “gold standard” Diffuse Optical Spectroscopy (DOS), a multi-frequency, multi-distance method combining frequency-domain photon migration and broadband spectroscopy [22,23]. Phantoms were constructed with cylindrical tube shaped void(s) 1.8 mm in diameter running parallel to the surface centered at depths of 2 and 3 mm. A mixture of a fatty emulsion (Lyposin II, Abbott Laboratories, Chicago, IL), nigrosin, and water was injected into the tubes. The scattering properties of the mixture were matched to that of the silicone phantom, whereas the absorption varied so that the contrast between the tube and background material was varied from 3:1 to 100:1.

3. Results and discussion

An example of a three-dimensional reconstructed image is shown in Fig. 2 . The phantom, depicted in Fig. 2(a), consisted of a single absorbing tube with an absorption contrast of 50:1. This is the absorption contrast in the near-infrared that we expect between superficial veins and the surrounding tissue. The center of the tube was located at a depth of 3 mm. Figure 2(b) is a volume rendering of the three-dimensional reconstructed absorption image. The tube surface in the volume rendering represents an isosurface of 60% of the maximum reconstructed value in the image. Individual slices through the three-dimension image at depths of 1-5 mm are shown in Fig. 2(c). At the surface the tube is not visible. The maximum amount of absorption occurs in the 1 mm slice. Absorption then decreases with depth until at a depth of 5 mm no absorption contrast is present.

 

Fig. 2 Example of a reconstructed image of a single absorbing tube located 3 mm below the surface of a tissue simulating phantom. (a) Schematic of the tissue simulating phantom. (b) Volume rendering of the three-dimensional reconstructed image. (c) Slices through the reconstructed image at depths of d = 1-5 mm.

Download Full Size | PPT Slide | PDF

Image reconstruction leads to a significant increase in resolution compared to images created by fitting for optical properties on pixel by pixel basis while assuming a homogeneous model. Profiles of tomographic images are shown in Fig. 3 . The profiles are from experiments with a single absorbing tube at depths of 2 mm (Fig. 3(a)) and 3 mm (Fig. 3(b)). In both cases the contrast in 50:1. The curves are normalized to a maximum value of one to facilitate a comparison of the two methods. The solid lines correspond to the lateral profiles of the reconstructed images. The full-width-half-maximum (FWHM) for the plots shown are 1.8 mm and 2.2 mm at depths of 2 mm and 3 mm respectively. We also measured the FWHM for tubes with contrasts ranging from 3:1 to 100:1. The FWHM ranged from 1.8 to 2.3 mm at a depth of 2 mm and from 2.2 to 3.2 mm at a depth of 3 mm, with the low contrast images having slightly lower resolution. For comparison, profiles of the images obtained by fitting pixel by pixel to a homogeneous model are shown with dotted lines. The FWHM is approximately three times larger. For the different concentrations, the FWHM determined from homogeneous model fits ranged from 4.4 to 6.8 mm at 2 mm depth, and from 5.5 to 9.5 mm at 3 mm depth.

 

Fig. 3 Line profiles of the reconstructed images of a single absorbing tube at depths of (a) 2mm and (b) 3 mm. Solid line corresponds to tomographic reconstruction, and dotted lines correspond to pixel by pixel fitting to a homogeneous model. Curves are normalized to facilitate comparison.

Download Full Size | PPT Slide | PDF

The improved resolution not only allows one to distinguish between nearby absorbers, but also decreases errors in quantification due to the blurring of nearby objects in the image. As a demonstration, we imaged four absorbing tubes at a depth of 2 mm. Each tube contained the same concentration of nigrosin, with the expected contrast being 50:1. The lateral spacing between the tubes was 3, 4, and 5 mm (see Fig. 4(a) ). In the reconstructed image and line profile (Fig. 4(b)) the tubes are clearly resolved and all have approximately the same reconstructed value of absorption. In contrast, for the image obtained by fitting pixel by pixel (Fig. 4(c)), the center two tubes appear to contain much more contrast than the outer tubes. This is because the image the absorption from nearby tubes, which appears broadened, is added to absorption of any given tube such that the central tubes falsely appear to have more absorption. Note, the maximum intensity also differs greatly, as images created using the homogeneous model underestimate the actual contrast by about a factor ten.

 

Fig. 4 (a) Schematic of tissue simulating phantom with four absorbing tubes located at a depth of z0 = 2 mm, and with lateral separations of d1 = 5 mm, d2 = 4 mm, and d3 = 3 mm. (b) Image and line profile of the reconstructed image using the tomographic method. (c) Image and line profile of the image produced by fitting to a homogeneous model.

Download Full Size | PPT Slide | PDF

In order to quantify the ability of both imaging methods to accurately recover the amount of absorption in a small sub-surface heterogeneity, we performed titration experiments in which the expected contrast in a single tube was varied from 3:1 to 100:1. Figure 5 shows the measured contrast versus the expected contrast using both the tomographic (Fig. 5(a,c)) and pixel by pixel (Fig. 5(b,d)) methods. The depth of the absorbing tube was 2 mm for Fig. 5(a,b), and 3 mm for Fig. 5(c,d). The contrast was determined by selecting a region of interest (ROI) 1 mm wide along the length of the tube and dividing the mean absorption value for the voxels in the ROI by the background value. For the tomographic images we selected all voxels from the depth at which the contrast was greatest. Error bars represent the standard deviation of voxels in the ROI.

 

Fig. 5 Plot of measured (reconstructed) vs. expected (known) contrast for absorbing tubes whose contrast ranged from 3:1 to 100:1 at depths of (a,b) 2 mm and (c,d) 3 mm. The tomographic method was used in (a) and (c), whereas (b) and (d) were calculated by fitting to a homogeneous model. Error bars denote the standard deviation of voxels within a region of interest (see text). The dotted line represents the known value.

Download Full Size | PPT Slide | PDF

Regardless of depth or imaging method, all curves are linear up to a contrast of 30:1, and then begin to saturate as expected due to the non-linearity of the inverse problem for diffuse light. We note that for this experiment, the linear range for tomographic reconstruction was larger than in previous work by a factor of about three [14]. However, non-linearity is expected to be more severe as the size of the inclusion becomes larger. Here the absorbing tube has a thickness of 1.5 , where as in the previous work the thickness of the absorber was 13 . Within the linear range both methods accurately give the relative change in contrast between the absorber and the background. However, the tomographic method comes closer to the actual values. Within the linear range, the percent error in measured absorption (i.e. 100 x (measured – expected)/expected) had a mean value of −2% and −28% at depths 2 mm and 3 mm respectively using the tomographic method.

The slight decrease in resolution at 3 mm suggests that the decrease in contrast measured at 3 mm is in part due to a partial volume effect. That is, the 3 mm tube appears broader and less intense due to the fact that the sensitivity of each detector is sharply peaked at the surface, and becomes both broader and smaller in magnitude with depth. While the reconstruction algorithm accounts for this effect, it does not remove it completely. Unlike the maximum absorption value, the total amount of reconstructed absorption does remain relatively constant with the depth of the tube. By integrating the reconstructed absorption of the entire tube, the total amount of reconstructed contrast is equal for the 2 mm and 3 mm depths to within 10%. In contrast, when fitting pixel by pixel to a homogeneous model the absolute values of absorption are underestimated by up to a factor of about ten. This was expected, since the trajectories of the detected photons are primarily located in background medium, and not in the absorbing tube.

We were able to determine the relative depths of absorbers, but not their absolute depths. Figure 6 show the depth profiles of the tomographic images of the single tubes at both depths. For clarity, only depth profiles for the 5x, 10x, 20x, 30x, and 50x concentrations are shown. The grey area denotes the range of depths at which the maximum absorption was reconstructed. For both the 2 and 3 mm depths, the reconstructed images underestimate the depth of the tube. The maximum absorption occurred from 0.6 to 0.7 mm and 0.9-1 mm for actual depths of 2 mm and 3 mm respectively. We attempted to use spatially varying regularization [24] to overcome this effect. However, for these particular experiments we found that the depth at which the tubes appeared in the image was more determined by our choice of regularization parameter, than by the actual depth of the absorber.

 

Fig. 6 Depth profiles of the three-dimensional images. For clarity only contrasts of 5x (dashed line), 10x (dash-dot line), 20x (dotted line), 30x (solid line with dots), and 50x (solid line) are shown. The grey area corresponds to the range in depth were the maximum values of absorption occur.

Download Full Size | PPT Slide | PDF

In these experiments the absorbers were located at 2-3 from the surface of the scattering medium, where the diffusion approximation to the radiative transport equation (RTE) is known to break down. It may be possible to more accurately determine the depths of buried absorbers by modeling photon transport using the RTE. Reconstruction algorithms such as the one used in this manuscript can be modified to incorporate radiative transport provided that the Fourier components of the Green’s functions (see Eq. (7)) can be calculated for the RTE [25,26]. We are currently exploring the use of both analytic [27] and Monte Carlo methods to solve for these Fourier components. Use of the RTE may also result in more accurate quantification and improved resolution.

5. Conclusion

We used structured illumination to generate three-dimensional images of absorption contrast. The use of a new image reconstruction formula allowed us to obtain improved quantification which was linear up to contrasts of 30:1 for a small absorber similar to a superficial vein. It also led to a reduction in the full-width-half-maximum in the images of approximately a factor of three. We were able to determine the relative depths of the absorbers based on the reconstructed images, but always underestimated the absolute depth.

Acknowledgements

This research was made possible by the Laser Microbeam and Medical Program ‖LAMMP|, an NIH Biomedical Technology Resource, Grant No. P41-RR01192; the Beckman Foundation; and the Military Photomedicine Program, AFOSR Grant No. FA9550-08-1-0384. Soren Konecky was supported by a fellowship from the Hewitt Foundation for Medical Research.

References and links

1. A. P. Gibson, J. C. Hebden, and S. R. Arridge, “Recent advances in diffuse optical imaging,” Phys. Med. Biol. 50(4), R1–R43 (2005). [PubMed]  

2. B. J. Tromberg, B. W. Pogue, K. D. Paulsen, A. G. Yodh, D. A. Boas, and A. E. Cerussi, “Assessing the future of diffuse optical imaging technologies for breast cancer management,” Med. Phys. 35(6), 2443–2451 (2008). [PubMed]  

3. J. C. Hebden, “Advances in optical imaging of the newborn infant brain,” Psychophysiology 40(4), 501–510 (2003). [PubMed]  

4. V. Ntziachristos, J. Ripoll, L. V. Wang, and R. Weissleder, “Looking and listening to light: the evolution of whole-body photonic imaging,” Nat. Biotechnol. 23(3), 313–320 (2005). [PubMed]  

5. D. J. Cuccia, F. Bevilacqua, A. J. Durkin, and B. J. Tromberg, “Modulated imaging: quantitative analysis and tomography of turbid media in the spatial-frequency domain,” Opt. Lett. 30(11), 1354–1356 (2005). [PubMed]  

6. A. Bassi, C. D’Andrea, G. Valentini, R. Cubeddu, and S. Arridge, “Temporal propagation of spatial information in turbid media,” Opt. Lett. 33(23), 2836–2838 (2008). [PubMed]  

7. A. Joshi, W. Bangerth, K. Hwang, J. C. Rasmussen, and E. M. Sevick-Muraca, “Fully adaptive FEM based fluorescence optical tomography from time-dependent measurements with area illumination and detection,” Med. Phys. 33(5), 1299–1310 (2006). [PubMed]  

8. D. J. Cuccia, F. Bevilacqua, A. J. Durkin, F. R. Ayers, and B. J. Tromberg, “Quantitation and mapping of tissue optical properties using modulated imaging,” J. Biomed. Opt. 14(2), 024012 (2009). [PubMed]  

9. D. Abookasis, C. C. Lay, M. S. Mathews, M. E. Linskey, R. D. Frostig, and B. J. Tromberg, “Imaging cortical absorption, scattering, and hemodynamic response during ischemic stroke using spatially modulated near-infrared illumination,” J. Biomed. Opt. 14(2), 024033 (2009). [PubMed]  

10. D. J. Cuccia, D. Abookasis, R. D. Frostig, and B. J. Tromberg, “Quantitative in vivo imaging of tissue absorption, scattering, and hemoglobin concentration in rat cortex using spatially-modulated structured light,” in In Vivo Optical Imaging of Brain Function, 2nd ed., R. D. Frostig, ed. (CRC, 2009).

11. J. R. Weber, D. J. Cuccia, A. J. Durkin, and B. J. Tromberg, “Noncontact imaging of absorption and scattering in layered tissue using spatially modulated structured light,” J. Appl. Phys. in press.

12. V. A. Markel and J. C. Schotland, “Symmetries, inversion formulas, and image reconstruction for optical tomography,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 70(5), 056616 (2004). [PubMed]  

13. J. C. Schotland and V. A. Markel, “Inverse scattering with diffusing waves,” J. Opt. Soc. Am. A 18(11), 2767–2777 (2001).

14. S. D. Konecky, G. Y. Panasyuk, K. Lee, V. Markel, A. G. Yodh, and J. C. Schotland, “Imaging complex structures with diffuse light,” Opt. Express 16(7), 5048–5060 (2008). [PubMed]  

15. Z. M. Wang, G. Y. Panasyuk, V. A. Markel, and J. C. Schotland, “Experimental demonstration of an analytic method for image reconstruction in optical diffusion tomography with large data sets,” Opt. Lett. 30(24), 3338–3340 (2005).

16. V. Lukic, V. A. Markel, and J. C. Schotland, “Optical tomography with structured illumination,” Opt. Lett. 34(7), 983–985 (2009). [PubMed]  

17. G. Y. Panasyuk, V. A. Markel, and J. C. Schotland, “Superresolution and corrections to the diffusion approximation in optical tomography,” Appl. Phys. Lett. 87(10), 101111 (2005).

18. V. A. Markel and J. C. Schotland, “Inverse problem in optical diffusion tomography. II. Role of boundary conditions,” J. Opt. Soc. Am. A 19(3), 558–566 (2002).

19. V. A. Markel, V. Mital, and J. C. Schotland, “Inverse problem in optical diffusion tomography. III. Inversion formulas and singular-value decomposition,” J. Opt. Soc. Am. A 20(5), 890–902 (2003).

20. A. C. Kak, and M. Slaney, Principles of Computerized Imaging (IEEE, 1988).

21. http://www.bli.uci.edu/ntroi/phantoms.php, retrieved April 24th, 2009.

22. F. Bevilacqua, A. J. Berger, A. E. Cerussi, D. Jakubowski, and B. J. Tromberg, “Broadband absorption spectroscopy in turbid media by combined frequency-domain and steady-state methods,” Appl. Opt. 39(34), 6498–6507 (2000).

23. H. P. Tuan, O. Coquoz, J. B. Fishkin, E. Anderson, and B. J. Tromberg, “Broad bandwidth frequency domain instrument for quantitative tissue optical spectroscopy,” Rev. Sci. Instrum. 71(6), 2500–2513 (2000).

24. 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(13), 2950–2961 (1999).

25. J. C. Schotland and V. A. Markel, “Fourier-Laplace structure of the inverse scattering problem for the radiative transport equation,” Inverse Problems and Imaging 1, 147–154 (2007).

26. G. Y. Panasyuk, J. C. Schotland, and V. A. Markel, “Radiative Transport Equation in Rotated Reference Frames,” J. Phys. A 39(1), 115–137 (2006).

27. V. A. Markel, “Modified spherical hamonics method for solving the radiative transport equation,” Waves Random Media 14(1), L13–L19 (2004).

References

  • View by:
  • |
  • |
  • |

  1. A. P. Gibson, J. C. Hebden, and S. R. Arridge, “Recent advances in diffuse optical imaging,” Phys. Med. Biol. 50(4), R1–R43 (2005).
    [PubMed]
  2. B. J. Tromberg, B. W. Pogue, K. D. Paulsen, A. G. Yodh, D. A. Boas, and A. E. Cerussi, “Assessing the future of diffuse optical imaging technologies for breast cancer management,” Med. Phys. 35(6), 2443–2451 (2008).
    [PubMed]
  3. J. C. Hebden, “Advances in optical imaging of the newborn infant brain,” Psychophysiology 40(4), 501–510 (2003).
    [PubMed]
  4. V. Ntziachristos, J. Ripoll, L. V. Wang, and R. Weissleder, “Looking and listening to light: the evolution of whole-body photonic imaging,” Nat. Biotechnol. 23(3), 313–320 (2005).
    [PubMed]
  5. D. J. Cuccia, F. Bevilacqua, A. J. Durkin, and B. J. Tromberg, “Modulated imaging: quantitative analysis and tomography of turbid media in the spatial-frequency domain,” Opt. Lett. 30(11), 1354–1356 (2005).
    [PubMed]
  6. A. Bassi, C. D’Andrea, G. Valentini, R. Cubeddu, and S. Arridge, “Temporal propagation of spatial information in turbid media,” Opt. Lett. 33(23), 2836–2838 (2008).
    [PubMed]
  7. A. Joshi, W. Bangerth, K. Hwang, J. C. Rasmussen, and E. M. Sevick-Muraca, “Fully adaptive FEM based fluorescence optical tomography from time-dependent measurements with area illumination and detection,” Med. Phys. 33(5), 1299–1310 (2006).
    [PubMed]
  8. D. J. Cuccia, F. Bevilacqua, A. J. Durkin, F. R. Ayers, and B. J. Tromberg, “Quantitation and mapping of tissue optical properties using modulated imaging,” J. Biomed. Opt. 14(2), 024012 (2009).
    [PubMed]
  9. D. Abookasis, C. C. Lay, M. S. Mathews, M. E. Linskey, R. D. Frostig, and B. J. Tromberg, “Imaging cortical absorption, scattering, and hemodynamic response during ischemic stroke using spatially modulated near-infrared illumination,” J. Biomed. Opt. 14(2), 024033 (2009).
    [PubMed]
  10. D. J. Cuccia, D. Abookasis, R. D. Frostig, and B. J. Tromberg, “Quantitative in vivo imaging of tissue absorption, scattering, and hemoglobin concentration in rat cortex using spatially-modulated structured light,” in In Vivo Optical Imaging of Brain Function, 2nd ed., R. D. Frostig, ed. (CRC, 2009).
  11. J. R. Weber, D. J. Cuccia, A. J. Durkin, and B. J. Tromberg, “Noncontact imaging of absorption and scattering in layered tissue using spatially modulated structured light,” J. Appl. Phys.in press.
  12. V. A. Markel and J. C. Schotland, “Symmetries, inversion formulas, and image reconstruction for optical tomography,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 70(5), 056616 (2004).
    [PubMed]
  13. J. C. Schotland and V. A. Markel, “Inverse scattering with diffusing waves,” J. Opt. Soc. Am. A 18(11), 2767–2777 (2001).
  14. S. D. Konecky, G. Y. Panasyuk, K. Lee, V. Markel, A. G. Yodh, and J. C. Schotland, “Imaging complex structures with diffuse light,” Opt. Express 16(7), 5048–5060 (2008).
    [PubMed]
  15. Z. M. Wang, G. Y. Panasyuk, V. A. Markel, and J. C. Schotland, “Experimental demonstration of an analytic method for image reconstruction in optical diffusion tomography with large data sets,” Opt. Lett. 30(24), 3338–3340 (2005).
  16. V. Lukic, V. A. Markel, and J. C. Schotland, “Optical tomography with structured illumination,” Opt. Lett. 34(7), 983–985 (2009).
    [PubMed]
  17. G. Y. Panasyuk, V. A. Markel, and J. C. Schotland, “Superresolution and corrections to the diffusion approximation in optical tomography,” Appl. Phys. Lett. 87(10), 101111 (2005).
  18. V. A. Markel and J. C. Schotland, “Inverse problem in optical diffusion tomography. II. Role of boundary conditions,” J. Opt. Soc. Am. A 19(3), 558–566 (2002).
  19. V. A. Markel, V. Mital, and J. C. Schotland, “Inverse problem in optical diffusion tomography. III. Inversion formulas and singular-value decomposition,” J. Opt. Soc. Am. A 20(5), 890–902 (2003).
  20. A. C. Kak, and M. Slaney, Principles of Computerized Imaging (IEEE, 1988).
  21. http://www.bli.uci.edu/ntroi/phantoms.php , retrieved April 24th, 2009.
  22. F. Bevilacqua, A. J. Berger, A. E. Cerussi, D. Jakubowski, and B. J. Tromberg, “Broadband absorption spectroscopy in turbid media by combined frequency-domain and steady-state methods,” Appl. Opt. 39(34), 6498–6507 (2000).
  23. H. P. Tuan, O. Coquoz, J. B. Fishkin, E. Anderson, and B. J. Tromberg, “Broad bandwidth frequency domain instrument for quantitative tissue optical spectroscopy,” Rev. Sci. Instrum. 71(6), 2500–2513 (2000).
  24. 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(13), 2950–2961 (1999).
  25. J. C. Schotland and V. A. Markel, “Fourier-Laplace structure of the inverse scattering problem for the radiative transport equation,” Inverse Problems and Imaging 1, 147–154 (2007).
  26. G. Y. Panasyuk, J. C. Schotland, and V. A. Markel, “Radiative Transport Equation in Rotated Reference Frames,” J. Phys. A 39(1), 115–137 (2006).
  27. V. A. Markel, “Modified spherical hamonics method for solving the radiative transport equation,” Waves Random Media 14(1), L13–L19 (2004).

2009 (3)

D. J. Cuccia, F. Bevilacqua, A. J. Durkin, F. R. Ayers, and B. J. Tromberg, “Quantitation and mapping of tissue optical properties using modulated imaging,” J. Biomed. Opt. 14(2), 024012 (2009).
[PubMed]

D. Abookasis, C. C. Lay, M. S. Mathews, M. E. Linskey, R. D. Frostig, and B. J. Tromberg, “Imaging cortical absorption, scattering, and hemodynamic response during ischemic stroke using spatially modulated near-infrared illumination,” J. Biomed. Opt. 14(2), 024033 (2009).
[PubMed]

V. Lukic, V. A. Markel, and J. C. Schotland, “Optical tomography with structured illumination,” Opt. Lett. 34(7), 983–985 (2009).
[PubMed]

2008 (3)

2007 (1)

J. C. Schotland and V. A. Markel, “Fourier-Laplace structure of the inverse scattering problem for the radiative transport equation,” Inverse Problems and Imaging 1, 147–154 (2007).

2006 (2)

G. Y. Panasyuk, J. C. Schotland, and V. A. Markel, “Radiative Transport Equation in Rotated Reference Frames,” J. Phys. A 39(1), 115–137 (2006).

A. Joshi, W. Bangerth, K. Hwang, J. C. Rasmussen, and E. M. Sevick-Muraca, “Fully adaptive FEM based fluorescence optical tomography from time-dependent measurements with area illumination and detection,” Med. Phys. 33(5), 1299–1310 (2006).
[PubMed]

2005 (5)

A. P. Gibson, J. C. Hebden, and S. R. Arridge, “Recent advances in diffuse optical imaging,” Phys. Med. Biol. 50(4), R1–R43 (2005).
[PubMed]

V. Ntziachristos, J. Ripoll, L. V. Wang, and R. Weissleder, “Looking and listening to light: the evolution of whole-body photonic imaging,” Nat. Biotechnol. 23(3), 313–320 (2005).
[PubMed]

G. Y. Panasyuk, V. A. Markel, and J. C. Schotland, “Superresolution and corrections to the diffusion approximation in optical tomography,” Appl. Phys. Lett. 87(10), 101111 (2005).

D. J. Cuccia, F. Bevilacqua, A. J. Durkin, and B. J. Tromberg, “Modulated imaging: quantitative analysis and tomography of turbid media in the spatial-frequency domain,” Opt. Lett. 30(11), 1354–1356 (2005).
[PubMed]

Z. M. Wang, G. Y. Panasyuk, V. A. Markel, and J. C. Schotland, “Experimental demonstration of an analytic method for image reconstruction in optical diffusion tomography with large data sets,” Opt. Lett. 30(24), 3338–3340 (2005).

2004 (2)

V. A. Markel and J. C. Schotland, “Symmetries, inversion formulas, and image reconstruction for optical tomography,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 70(5), 056616 (2004).
[PubMed]

V. A. Markel, “Modified spherical hamonics method for solving the radiative transport equation,” Waves Random Media 14(1), L13–L19 (2004).

2003 (2)

2002 (1)

2001 (1)

2000 (2)

F. Bevilacqua, A. J. Berger, A. E. Cerussi, D. Jakubowski, and B. J. Tromberg, “Broadband absorption spectroscopy in turbid media by combined frequency-domain and steady-state methods,” Appl. Opt. 39(34), 6498–6507 (2000).

H. P. Tuan, O. Coquoz, J. B. Fishkin, E. Anderson, and B. J. Tromberg, “Broad bandwidth frequency domain instrument for quantitative tissue optical spectroscopy,” Rev. Sci. Instrum. 71(6), 2500–2513 (2000).

1999 (1)

Abookasis, D.

D. Abookasis, C. C. Lay, M. S. Mathews, M. E. Linskey, R. D. Frostig, and B. J. Tromberg, “Imaging cortical absorption, scattering, and hemodynamic response during ischemic stroke using spatially modulated near-infrared illumination,” J. Biomed. Opt. 14(2), 024033 (2009).
[PubMed]

Anderson, E.

H. P. Tuan, O. Coquoz, J. B. Fishkin, E. Anderson, and B. J. Tromberg, “Broad bandwidth frequency domain instrument for quantitative tissue optical spectroscopy,” Rev. Sci. Instrum. 71(6), 2500–2513 (2000).

Arridge, S.

Arridge, S. R.

A. P. Gibson, J. C. Hebden, and S. R. Arridge, “Recent advances in diffuse optical imaging,” Phys. Med. Biol. 50(4), R1–R43 (2005).
[PubMed]

Ayers, F. R.

D. J. Cuccia, F. Bevilacqua, A. J. Durkin, F. R. Ayers, and B. J. Tromberg, “Quantitation and mapping of tissue optical properties using modulated imaging,” J. Biomed. Opt. 14(2), 024012 (2009).
[PubMed]

Bangerth, W.

A. Joshi, W. Bangerth, K. Hwang, J. C. Rasmussen, and E. M. Sevick-Muraca, “Fully adaptive FEM based fluorescence optical tomography from time-dependent measurements with area illumination and detection,” Med. Phys. 33(5), 1299–1310 (2006).
[PubMed]

Bassi, A.

Berger, A. J.

Bevilacqua, F.

Boas, D. A.

B. J. Tromberg, B. W. Pogue, K. D. Paulsen, A. G. Yodh, D. A. Boas, and A. E. Cerussi, “Assessing the future of diffuse optical imaging technologies for breast cancer management,” Med. Phys. 35(6), 2443–2451 (2008).
[PubMed]

Cerussi, A. E.

B. J. Tromberg, B. W. Pogue, K. D. Paulsen, A. G. Yodh, D. A. Boas, and A. E. Cerussi, “Assessing the future of diffuse optical imaging technologies for breast cancer management,” Med. Phys. 35(6), 2443–2451 (2008).
[PubMed]

F. Bevilacqua, A. J. Berger, A. E. Cerussi, D. Jakubowski, and B. J. Tromberg, “Broadband absorption spectroscopy in turbid media by combined frequency-domain and steady-state methods,” Appl. Opt. 39(34), 6498–6507 (2000).

Coquoz, O.

H. P. Tuan, O. Coquoz, J. B. Fishkin, E. Anderson, and B. J. Tromberg, “Broad bandwidth frequency domain instrument for quantitative tissue optical spectroscopy,” Rev. Sci. Instrum. 71(6), 2500–2513 (2000).

Cubeddu, R.

Cuccia, D. J.

D. J. Cuccia, F. Bevilacqua, A. J. Durkin, F. R. Ayers, and B. J. Tromberg, “Quantitation and mapping of tissue optical properties using modulated imaging,” J. Biomed. Opt. 14(2), 024012 (2009).
[PubMed]

D. J. Cuccia, F. Bevilacqua, A. J. Durkin, and B. J. Tromberg, “Modulated imaging: quantitative analysis and tomography of turbid media in the spatial-frequency domain,” Opt. Lett. 30(11), 1354–1356 (2005).
[PubMed]

J. R. Weber, D. J. Cuccia, A. J. Durkin, and B. J. Tromberg, “Noncontact imaging of absorption and scattering in layered tissue using spatially modulated structured light,” J. Appl. Phys.in press.

D’Andrea, C.

Durkin, A. J.

D. J. Cuccia, F. Bevilacqua, A. J. Durkin, F. R. Ayers, and B. J. Tromberg, “Quantitation and mapping of tissue optical properties using modulated imaging,” J. Biomed. Opt. 14(2), 024012 (2009).
[PubMed]

D. J. Cuccia, F. Bevilacqua, A. J. Durkin, and B. J. Tromberg, “Modulated imaging: quantitative analysis and tomography of turbid media in the spatial-frequency domain,” Opt. Lett. 30(11), 1354–1356 (2005).
[PubMed]

J. R. Weber, D. J. Cuccia, A. J. Durkin, and B. J. Tromberg, “Noncontact imaging of absorption and scattering in layered tissue using spatially modulated structured light,” J. Appl. Phys.in press.

Fishkin, J. B.

H. P. Tuan, O. Coquoz, J. B. Fishkin, E. Anderson, and B. J. Tromberg, “Broad bandwidth frequency domain instrument for quantitative tissue optical spectroscopy,” Rev. Sci. Instrum. 71(6), 2500–2513 (2000).

Frostig, R. D.

D. Abookasis, C. C. Lay, M. S. Mathews, M. E. Linskey, R. D. Frostig, and B. J. Tromberg, “Imaging cortical absorption, scattering, and hemodynamic response during ischemic stroke using spatially modulated near-infrared illumination,” J. Biomed. Opt. 14(2), 024033 (2009).
[PubMed]

Gibson, A. P.

A. P. Gibson, J. C. Hebden, and S. R. Arridge, “Recent advances in diffuse optical imaging,” Phys. Med. Biol. 50(4), R1–R43 (2005).
[PubMed]

Hebden, J. C.

A. P. Gibson, J. C. Hebden, and S. R. Arridge, “Recent advances in diffuse optical imaging,” Phys. Med. Biol. 50(4), R1–R43 (2005).
[PubMed]

J. C. Hebden, “Advances in optical imaging of the newborn infant brain,” Psychophysiology 40(4), 501–510 (2003).
[PubMed]

Hwang, K.

A. Joshi, W. Bangerth, K. Hwang, J. C. Rasmussen, and E. M. Sevick-Muraca, “Fully adaptive FEM based fluorescence optical tomography from time-dependent measurements with area illumination and detection,” Med. Phys. 33(5), 1299–1310 (2006).
[PubMed]

Jakubowski, D.

Joshi, A.

A. Joshi, W. Bangerth, K. Hwang, J. C. Rasmussen, and E. M. Sevick-Muraca, “Fully adaptive FEM based fluorescence optical tomography from time-dependent measurements with area illumination and detection,” Med. Phys. 33(5), 1299–1310 (2006).
[PubMed]

Konecky, S. D.

Lay, C. C.

D. Abookasis, C. C. Lay, M. S. Mathews, M. E. Linskey, R. D. Frostig, and B. J. Tromberg, “Imaging cortical absorption, scattering, and hemodynamic response during ischemic stroke using spatially modulated near-infrared illumination,” J. Biomed. Opt. 14(2), 024033 (2009).
[PubMed]

Lee, K.

Linskey, M. E.

D. Abookasis, C. C. Lay, M. S. Mathews, M. E. Linskey, R. D. Frostig, and B. J. Tromberg, “Imaging cortical absorption, scattering, and hemodynamic response during ischemic stroke using spatially modulated near-infrared illumination,” J. Biomed. Opt. 14(2), 024033 (2009).
[PubMed]

Lukic, V.

Markel, V.

Markel, V. A.

V. Lukic, V. A. Markel, and J. C. Schotland, “Optical tomography with structured illumination,” Opt. Lett. 34(7), 983–985 (2009).
[PubMed]

J. C. Schotland and V. A. Markel, “Fourier-Laplace structure of the inverse scattering problem for the radiative transport equation,” Inverse Problems and Imaging 1, 147–154 (2007).

G. Y. Panasyuk, J. C. Schotland, and V. A. Markel, “Radiative Transport Equation in Rotated Reference Frames,” J. Phys. A 39(1), 115–137 (2006).

Z. M. Wang, G. Y. Panasyuk, V. A. Markel, and J. C. Schotland, “Experimental demonstration of an analytic method for image reconstruction in optical diffusion tomography with large data sets,” Opt. Lett. 30(24), 3338–3340 (2005).

G. Y. Panasyuk, V. A. Markel, and J. C. Schotland, “Superresolution and corrections to the diffusion approximation in optical tomography,” Appl. Phys. Lett. 87(10), 101111 (2005).

V. A. Markel, “Modified spherical hamonics method for solving the radiative transport equation,” Waves Random Media 14(1), L13–L19 (2004).

V. A. Markel and J. C. Schotland, “Symmetries, inversion formulas, and image reconstruction for optical tomography,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 70(5), 056616 (2004).
[PubMed]

V. A. Markel, V. Mital, and J. C. Schotland, “Inverse problem in optical diffusion tomography. III. Inversion formulas and singular-value decomposition,” J. Opt. Soc. Am. A 20(5), 890–902 (2003).

V. A. Markel and J. C. Schotland, “Inverse problem in optical diffusion tomography. II. Role of boundary conditions,” J. Opt. Soc. Am. A 19(3), 558–566 (2002).

J. C. Schotland and V. A. Markel, “Inverse scattering with diffusing waves,” J. Opt. Soc. Am. A 18(11), 2767–2777 (2001).

Mathews, M. S.

D. Abookasis, C. C. Lay, M. S. Mathews, M. E. Linskey, R. D. Frostig, and B. J. Tromberg, “Imaging cortical absorption, scattering, and hemodynamic response during ischemic stroke using spatially modulated near-infrared illumination,” J. Biomed. Opt. 14(2), 024033 (2009).
[PubMed]

McBride, T. O.

Mital, V.

Ntziachristos, V.

V. Ntziachristos, J. Ripoll, L. V. Wang, and R. Weissleder, “Looking and listening to light: the evolution of whole-body photonic imaging,” Nat. Biotechnol. 23(3), 313–320 (2005).
[PubMed]

Osterberg, U. L.

Panasyuk, G. Y.

S. D. Konecky, G. Y. Panasyuk, K. Lee, V. Markel, A. G. Yodh, and J. C. Schotland, “Imaging complex structures with diffuse light,” Opt. Express 16(7), 5048–5060 (2008).
[PubMed]

G. Y. Panasyuk, J. C. Schotland, and V. A. Markel, “Radiative Transport Equation in Rotated Reference Frames,” J. Phys. A 39(1), 115–137 (2006).

Z. M. Wang, G. Y. Panasyuk, V. A. Markel, and J. C. Schotland, “Experimental demonstration of an analytic method for image reconstruction in optical diffusion tomography with large data sets,” Opt. Lett. 30(24), 3338–3340 (2005).

G. Y. Panasyuk, V. A. Markel, and J. C. Schotland, “Superresolution and corrections to the diffusion approximation in optical tomography,” Appl. Phys. Lett. 87(10), 101111 (2005).

Paulsen, K. D.

B. J. Tromberg, B. W. Pogue, K. D. Paulsen, A. G. Yodh, D. A. Boas, and A. E. Cerussi, “Assessing the future of diffuse optical imaging technologies for breast cancer management,” Med. Phys. 35(6), 2443–2451 (2008).
[PubMed]

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(13), 2950–2961 (1999).

Pogue, B. W.

B. J. Tromberg, B. W. Pogue, K. D. Paulsen, A. G. Yodh, D. A. Boas, and A. E. Cerussi, “Assessing the future of diffuse optical imaging technologies for breast cancer management,” Med. Phys. 35(6), 2443–2451 (2008).
[PubMed]

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(13), 2950–2961 (1999).

Prewitt, J.

Rasmussen, J. C.

A. Joshi, W. Bangerth, K. Hwang, J. C. Rasmussen, and E. M. Sevick-Muraca, “Fully adaptive FEM based fluorescence optical tomography from time-dependent measurements with area illumination and detection,” Med. Phys. 33(5), 1299–1310 (2006).
[PubMed]

Ripoll, J.

V. Ntziachristos, J. Ripoll, L. V. Wang, and R. Weissleder, “Looking and listening to light: the evolution of whole-body photonic imaging,” Nat. Biotechnol. 23(3), 313–320 (2005).
[PubMed]

Schotland, J. C.

V. Lukic, V. A. Markel, and J. C. Schotland, “Optical tomography with structured illumination,” Opt. Lett. 34(7), 983–985 (2009).
[PubMed]

S. D. Konecky, G. Y. Panasyuk, K. Lee, V. Markel, A. G. Yodh, and J. C. Schotland, “Imaging complex structures with diffuse light,” Opt. Express 16(7), 5048–5060 (2008).
[PubMed]

J. C. Schotland and V. A. Markel, “Fourier-Laplace structure of the inverse scattering problem for the radiative transport equation,” Inverse Problems and Imaging 1, 147–154 (2007).

G. Y. Panasyuk, J. C. Schotland, and V. A. Markel, “Radiative Transport Equation in Rotated Reference Frames,” J. Phys. A 39(1), 115–137 (2006).

Z. M. Wang, G. Y. Panasyuk, V. A. Markel, and J. C. Schotland, “Experimental demonstration of an analytic method for image reconstruction in optical diffusion tomography with large data sets,” Opt. Lett. 30(24), 3338–3340 (2005).

G. Y. Panasyuk, V. A. Markel, and J. C. Schotland, “Superresolution and corrections to the diffusion approximation in optical tomography,” Appl. Phys. Lett. 87(10), 101111 (2005).

V. A. Markel and J. C. Schotland, “Symmetries, inversion formulas, and image reconstruction for optical tomography,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 70(5), 056616 (2004).
[PubMed]

V. A. Markel, V. Mital, and J. C. Schotland, “Inverse problem in optical diffusion tomography. III. Inversion formulas and singular-value decomposition,” J. Opt. Soc. Am. A 20(5), 890–902 (2003).

V. A. Markel and J. C. Schotland, “Inverse problem in optical diffusion tomography. II. Role of boundary conditions,” J. Opt. Soc. Am. A 19(3), 558–566 (2002).

J. C. Schotland and V. A. Markel, “Inverse scattering with diffusing waves,” J. Opt. Soc. Am. A 18(11), 2767–2777 (2001).

Sevick-Muraca, E. M.

A. Joshi, W. Bangerth, K. Hwang, J. C. Rasmussen, and E. M. Sevick-Muraca, “Fully adaptive FEM based fluorescence optical tomography from time-dependent measurements with area illumination and detection,” Med. Phys. 33(5), 1299–1310 (2006).
[PubMed]

Tromberg, B. J.

D. Abookasis, C. C. Lay, M. S. Mathews, M. E. Linskey, R. D. Frostig, and B. J. Tromberg, “Imaging cortical absorption, scattering, and hemodynamic response during ischemic stroke using spatially modulated near-infrared illumination,” J. Biomed. Opt. 14(2), 024033 (2009).
[PubMed]

D. J. Cuccia, F. Bevilacqua, A. J. Durkin, F. R. Ayers, and B. J. Tromberg, “Quantitation and mapping of tissue optical properties using modulated imaging,” J. Biomed. Opt. 14(2), 024012 (2009).
[PubMed]

B. J. Tromberg, B. W. Pogue, K. D. Paulsen, A. G. Yodh, D. A. Boas, and A. E. Cerussi, “Assessing the future of diffuse optical imaging technologies for breast cancer management,” Med. Phys. 35(6), 2443–2451 (2008).
[PubMed]

D. J. Cuccia, F. Bevilacqua, A. J. Durkin, and B. J. Tromberg, “Modulated imaging: quantitative analysis and tomography of turbid media in the spatial-frequency domain,” Opt. Lett. 30(11), 1354–1356 (2005).
[PubMed]

H. P. Tuan, O. Coquoz, J. B. Fishkin, E. Anderson, and B. J. Tromberg, “Broad bandwidth frequency domain instrument for quantitative tissue optical spectroscopy,” Rev. Sci. Instrum. 71(6), 2500–2513 (2000).

F. Bevilacqua, A. J. Berger, A. E. Cerussi, D. Jakubowski, and B. J. Tromberg, “Broadband absorption spectroscopy in turbid media by combined frequency-domain and steady-state methods,” Appl. Opt. 39(34), 6498–6507 (2000).

J. R. Weber, D. J. Cuccia, A. J. Durkin, and B. J. Tromberg, “Noncontact imaging of absorption and scattering in layered tissue using spatially modulated structured light,” J. Appl. Phys.in press.

Tuan, H. P.

H. P. Tuan, O. Coquoz, J. B. Fishkin, E. Anderson, and B. J. Tromberg, “Broad bandwidth frequency domain instrument for quantitative tissue optical spectroscopy,” Rev. Sci. Instrum. 71(6), 2500–2513 (2000).

Valentini, G.

Wang, L. V.

V. Ntziachristos, J. Ripoll, L. V. Wang, and R. Weissleder, “Looking and listening to light: the evolution of whole-body photonic imaging,” Nat. Biotechnol. 23(3), 313–320 (2005).
[PubMed]

Wang, Z. M.

Weber, J. R.

J. R. Weber, D. J. Cuccia, A. J. Durkin, and B. J. Tromberg, “Noncontact imaging of absorption and scattering in layered tissue using spatially modulated structured light,” J. Appl. Phys.in press.

Weissleder, R.

V. Ntziachristos, J. Ripoll, L. V. Wang, and R. Weissleder, “Looking and listening to light: the evolution of whole-body photonic imaging,” Nat. Biotechnol. 23(3), 313–320 (2005).
[PubMed]

Yodh, A. G.

S. D. Konecky, G. Y. Panasyuk, K. Lee, V. Markel, A. G. Yodh, and J. C. Schotland, “Imaging complex structures with diffuse light,” Opt. Express 16(7), 5048–5060 (2008).
[PubMed]

B. J. Tromberg, B. W. Pogue, K. D. Paulsen, A. G. Yodh, D. A. Boas, and A. E. Cerussi, “Assessing the future of diffuse optical imaging technologies for breast cancer management,” Med. Phys. 35(6), 2443–2451 (2008).
[PubMed]

Appl. Opt. (2)

Appl. Phys. Lett. (1)

G. Y. Panasyuk, V. A. Markel, and J. C. Schotland, “Superresolution and corrections to the diffusion approximation in optical tomography,” Appl. Phys. Lett. 87(10), 101111 (2005).

Inverse Problems and Imaging (1)

J. C. Schotland and V. A. Markel, “Fourier-Laplace structure of the inverse scattering problem for the radiative transport equation,” Inverse Problems and Imaging 1, 147–154 (2007).

J. Appl. Phys. (1)

J. R. Weber, D. J. Cuccia, A. J. Durkin, and B. J. Tromberg, “Noncontact imaging of absorption and scattering in layered tissue using spatially modulated structured light,” J. Appl. Phys.in press.

J. Biomed. Opt. (2)

D. J. Cuccia, F. Bevilacqua, A. J. Durkin, F. R. Ayers, and B. J. Tromberg, “Quantitation and mapping of tissue optical properties using modulated imaging,” J. Biomed. Opt. 14(2), 024012 (2009).
[PubMed]

D. Abookasis, C. C. Lay, M. S. Mathews, M. E. Linskey, R. D. Frostig, and B. J. Tromberg, “Imaging cortical absorption, scattering, and hemodynamic response during ischemic stroke using spatially modulated near-infrared illumination,” J. Biomed. Opt. 14(2), 024033 (2009).
[PubMed]

J. Opt. Soc. Am. A (3)

J. Phys. A (1)

G. Y. Panasyuk, J. C. Schotland, and V. A. Markel, “Radiative Transport Equation in Rotated Reference Frames,” J. Phys. A 39(1), 115–137 (2006).

Med. Phys. (2)

A. Joshi, W. Bangerth, K. Hwang, J. C. Rasmussen, and E. M. Sevick-Muraca, “Fully adaptive FEM based fluorescence optical tomography from time-dependent measurements with area illumination and detection,” Med. Phys. 33(5), 1299–1310 (2006).
[PubMed]

B. J. Tromberg, B. W. Pogue, K. D. Paulsen, A. G. Yodh, D. A. Boas, and A. E. Cerussi, “Assessing the future of diffuse optical imaging technologies for breast cancer management,” Med. Phys. 35(6), 2443–2451 (2008).
[PubMed]

Nat. Biotechnol. (1)

V. Ntziachristos, J. Ripoll, L. V. Wang, and R. Weissleder, “Looking and listening to light: the evolution of whole-body photonic imaging,” Nat. Biotechnol. 23(3), 313–320 (2005).
[PubMed]

Opt. Express (1)

Opt. Lett. (4)

Phys. Med. Biol. (1)

A. P. Gibson, J. C. Hebden, and S. R. Arridge, “Recent advances in diffuse optical imaging,” Phys. Med. Biol. 50(4), R1–R43 (2005).
[PubMed]

Phys. Rev. E Stat. Nonlin. Soft Matter Phys. (1)

V. A. Markel and J. C. Schotland, “Symmetries, inversion formulas, and image reconstruction for optical tomography,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 70(5), 056616 (2004).
[PubMed]

Psychophysiology (1)

J. C. Hebden, “Advances in optical imaging of the newborn infant brain,” Psychophysiology 40(4), 501–510 (2003).
[PubMed]

Rev. Sci. Instrum. (1)

H. P. Tuan, O. Coquoz, J. B. Fishkin, E. Anderson, and B. J. Tromberg, “Broad bandwidth frequency domain instrument for quantitative tissue optical spectroscopy,” Rev. Sci. Instrum. 71(6), 2500–2513 (2000).

Waves Random Media (1)

V. A. Markel, “Modified spherical hamonics method for solving the radiative transport equation,” Waves Random Media 14(1), L13–L19 (2004).

Other (3)

A. C. Kak, and M. Slaney, Principles of Computerized Imaging (IEEE, 1988).

http://www.bli.uci.edu/ntroi/phantoms.php , retrieved April 24th, 2009.

D. J. Cuccia, D. Abookasis, R. D. Frostig, and B. J. Tromberg, “Quantitative in vivo imaging of tissue absorption, scattering, and hemoglobin concentration in rat cortex using spatially-modulated structured light,” in In Vivo Optical Imaging of Brain Function, 2nd ed., R. D. Frostig, ed. (CRC, 2009).

Cited By

OSA participates in CrossRef's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Fig. 1
Fig. 1

Schematic of the experimental setup. Broadband light from a lamp is projected on the sample using a digital micro-mirror device (DMD) and detected at discrete wavelengths using a liquid crystal tunable filter and CCD. Crossed polarizers eliminate specular reflections.

Fig. 2
Fig. 2

Example of a reconstructed image of a single absorbing tube located 3 mm below the surface of a tissue simulating phantom. (a) Schematic of the tissue simulating phantom. (b) Volume rendering of the three-dimensional reconstructed image. (c) Slices through the reconstructed image at depths of d = 1-5 mm.

Fig. 3
Fig. 3

Line profiles of the reconstructed images of a single absorbing tube at depths of (a) 2mm and (b) 3 mm. Solid line corresponds to tomographic reconstruction, and dotted lines correspond to pixel by pixel fitting to a homogeneous model. Curves are normalized to facilitate comparison.

Fig. 4
Fig. 4

(a) Schematic of tissue simulating phantom with four absorbing tubes located at a depth of z0 = 2 mm, and with lateral separations of d1 = 5 mm, d2 = 4 mm, and d3 = 3 mm. (b) Image and line profile of the reconstructed image using the tomographic method. (c) Image and line profile of the image produced by fitting to a homogeneous model.

Fig. 5
Fig. 5

Plot of measured (reconstructed) vs. expected (known) contrast for absorbing tubes whose contrast ranged from 3:1 to 100:1 at depths of (a,b) 2 mm and (c,d) 3 mm. The tomographic method was used in (a) and (c), whereas (b) and (d) were calculated by fitting to a homogeneous model. Error bars denote the standard deviation of voxels within a region of interest (see text). The dotted line represents the known value.

Fig. 6
Fig. 6

Depth profiles of the three-dimensional images. For clarity only contrasts of 5x (dashed line), 10x (dash-dot line), 20x (dotted line), 30x (solid line with dots), and 50x (solid line) are shown. The grey area corresponds to the range in depth were the maximum values of absorption occur.

Equations (18)

Equations on this page are rendered with MathJax. Learn more.

s^I(r,s^)+(μa+μs)I(r,s^)μsd2s'A(s^,s^')I(r,s^')=S(r,s^)  .
I(r,s^)=c4πΦ(r)+34πJs^  .
G(r,s^,r',s^')=c4π(1+s^r)(1s^'r')G(r,r') ,
[D+ cμa] G(r,r') = δ(rr') ,
G(r,r')+n^G(r,r')=0  .
G(r,r')=1(2π)2d2q g(q,z,z')exp[i(q(ρ'ρ))] ,
g(q,z,z')=Dexp[Q(q)|zz'|]Q(q)+1
φ(k,rd,s^d)=d3r d2s I(r,s^) G(r,s^,rd,s^d) δμa(r) ,
I(r,s^)=d3r'd2s' G(r',s^',r,s^)S(r',s^')  .
S(r,s^)=I0exp(ikρ)δ(z)δ(s^+n^)  .
φ(k,ρd)=c I04π(1+)2d3r K(k,ρd;r) c δμa(r) ,
K(k,ρd;r)=1(2π)2d2q κ(k,q;z) exp[i(kq)ρ+iqρd] ,
κ(k,q;z)=(D)2[123(kq+Q(k)Q(q))]exp[(Q(k)+Q(q)) z](Q(k)+1)(Q(q)+1)  .
S(k,ρ,θ)=I02(1+Mcos(kρ+θ)) δ(z) δ(s^+n^)  .
Itot=13(2I2I1I3)+i3(I1I3)  .
A=23{(I1I2)2+(I2I3)2+(I3I1)2}1/2 ,
ϕ=tan1{3(I1I3)/(2I2I1I3)}  .
φRytov=Itot(0){log(A/A(0))+i(ϕϕ(0))}  .

Metrics