Reconstruction of in-vivo optical properties for human prostate using interstitial diffuse optical tomography

doi:10.1364/OE.17.011665

« Show journal navigation

Reconstruction of in-vivo optical properties for human prostate using interstitial diffuse optical tomography

Ken Kang-Hsin Wang and Timothy C. Zhu »View Author Affiliations

Optics Express, Vol. 17, Issue 14, pp. 11665-11672 (2009)
http://dx.doi.org/10.1364/OE.17.011665

View Full Text Article

Acrobat PDF (633 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Methods
3. Results and discussion
4. Conclusion
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (14)
Cited By
Figures (6)
Metrics

Abstract

A CW interstitial diffuse optical tomography has been developed to characterize the in-vivo optical properties of prostate gland during photodynamic therapy. The spatial distributions of light fluence rate can be described by the diffusion equation. Optical properties of the prostate are reconstructed by solving the inverse problem with an adjoint method. The 3D reconstructed in-vivo optical properties for a human prostate is illustrated and compared with the results generated by a well-established point-by-point method. Moreover, the calculated fluence rate using the reconstructed optical properties matches the measured data.

1. Introduction

The therapeutic efficacy of photodynamic therapy (PDT) significantly depends on the quantification of light fluence rate distribution, sensitizer concentration, and oxygen transport and consumption during treatment. With the implementation of interstitial diffuse optical tomography (iDOT), the optical properties of tissue can be reconstructed spatially using an adjoint inverse algorithm [1

H. Dehgani, M. E. Eames, P. K. Yalavarthy, S. C. Davis, S. Srinivasan, C. M. Carpenter, B. W. Pogue, and K. D. Paulsen, “Near infrared optical tomogrophy using NIRFAST: Algorithm for numerical model and image reconstruction,” Commun. Numer. Methods Eng. 25, 711–732 (2008). [CrossRef]

]. Accurately characterizing the abnormal optical properties of tumor and thereby rigorously estimate patient-specific light fluence rate distribution is an important subject for PDT light dosimetry.

Prostate cancer is the most common cancer expected to occur in men and is also the second leading cause of cancer death in year 2008 among American male population [2]. In addition, there is a need for fast and economical diagnostic tool. Recently, there is an increasing interest in the field of continuous wave (CW)-DOT. Axelsson et al. [3

J. Axelsson, J. Swartling, and S. Andersson-Engels, “In vivo photosensitizer tomography inside the human prostate,” Opt. Lett. 34, 232–234 (2009). [CrossRef] [PubMed]

] present promising results of reconstructing spatial distributions of photosensitizer in human prostate using iDOT, which is expected to facilitate the use of light-photosensitizer dosimetry model in a PDT clinical environment. A new globally accelerated reconstruction algorithm for CW-DOT is developed by Pantong et al. [4

N. Pantong, J. Su, H. Shan, M. V. Klibanov, and H. Liu, “Globally accelerated reconstruction algorithm for diffusion tomography with continuous-wave source in an arbitrary convex shape domain,” J. Opt. Soc. Am. A. 26, 456–472 (2009). [CrossRef]

]. It has been demonstrated as a highly efficient and stable method for a large data set obtained from an organ with particular shape, such as prostate.

An iDOT system with multiple light diffusers and isotropic detectors has been developed in our group to characterize the optical properties of prostate gland during PDT [5

T. C. Zhu, J. C. Finlay, and S. M. Hahn, “Determination of the distribution of light, optical properties, drug concentration, and tissue oxygenation in-vivo in human prostate during motexafin lutetium-mediated photodynamic therapy,” J. Photochem. Photobiol. B 79, 231–241 (2005). [CrossRef] [PubMed]

, 6

T. C. Zhu, A. Dimofte, J. C. Finlay, D. Stripp, T. Busch, J. Miles, R. Whittington, S. B. Malkowicz, Z. Tochner, E. Glatstein, and S. M. Hahn, “Optical properties of human prostate at 732 nm measured in mediated photodynamic therapy,” Photochem. Photobiol. 81, 96–105 (2005). [CrossRef]

]. Several catheters paralleled with each other are inserted into prostate under ultrasound guidance. Inside the catheters, there are light point sources and detectors. For each source position, the isotropic detectors scan and record the light fluence rate distribution in the direction along with the catheter axis. For our CW-iDOT system, the spatial distributions of light fluence rate can be described by the steady state diffusion equation. Optical properties of the prostate gland are reconstructed by solving the inverse problem with the use of a CW diffusion model.

To exam our methodology, a mathematical phantom including anomalies with known optical properties is examined. Furthermore, the 3D reconstructed in-vivo optical properties for a human prostate is illustrated and compared with the results generated by an established point-by-point method [6

]. Moreover, we show the calculated fluence rates match the measured data using the reconstructed optical properties.

2. Methods

2.1 Description of the patient being studied

A Phase I clinical trial of motexafin lutetium (MLu)-mediated PDT in patients with locally recurrent prostate carcinoma was initiated at the University of Pennsylvania. The protocol was approved by the Cancer Therapy Evaluation Program of the National Cancer Institute, and the detail of the trial was published elsewhere [6

]. Briefly, two weeks prior to the scheduled treatment a transrectal ultrasound (TRUS) was performed for treatment planning. An urologist drew the target prostate volume on each slice of the ultrasound images and these images were spaced 0.5 cm apart. A built-in template with a 0.5-cm grid projected the locations of light sources relative to the prostate. The interstitial-PDT treatment was conducted using cylindrical diffusing fibers (CDF). A treatment plan was then prepared to determine the location and length of the CDFs with active length 1–5 cm to cover the full volume of the prostate. The treatment light sources were spaced one centimeter apart and the light power per unit length was less than or equal to 150 mW/cm. Prior to PDT, point light source (Pioneer optics company, Bloomfield, CT) for optical properties measurement were inserted into patient prostate at the same locations of CDF. For one source position, there were four isotropic detectors (Rare Earth Medical Inc, West Yarmouth, MA) surrounding it (Fig. 1 (A)). These detectors are the same ones used to perform light fluence rate measurements during PDT. For the study patient shown in this work (#16), five source and 12 detector locations were used. A 15-W 732 nm diode laser (model 730, Diomed, Ltd., Cambridge, United Kingdom) was chosen as the light source.

2.2 Reconstruction procedure-DOT method

In the first step, prostate contour is obtained from TRUS image modality. The contours of prostate, urethra, and rectum for a certain z-plane are shown in Fig. 1(A). Circle and cross symbols in Fig. 1(A) represent the locations of detectors and sources, respectively. For the sources at a certain z location, our detectors scan along the z-axis and record light fluence rate. The locations where the isotropic detectors recorded data are set equal to the virtual detector positions assigned in the subsequent reconstruction calculation. Figure 1(B) shows the sources and virtual detectors orientation in a 3D geometry. The separation between each detector along the z-axis is 1 mm. The 2D prostate contour for a particular z plane (Fig. 1(A)) is expanded in z-axis direction to construct a 3D structure (Fig. 1(B)) used in the calculation. This 3D geometry is processed in MATLAB 7.5 (MathWorks, Inc., Natick, MA, United States) and COMSOL 3.5a (COMSOL AB, SE-111 40, Stockholm, Sweden) program to generate mesh for the subsequent finite element method (FEM) calculation. It is worthwhile to mention that for different source locations at z-axis, the above procedure will be repeated to generate the corresponding 3D geometry.

Our goal is to recover the optical properties (absorption and reduced scattering coefficients, µ_a and µ_s ’) at each finite element mesh node using the fluence rate data ϕ_m measured interstitially. Unlike our prior studies [7

X. Zhou and T. C. Zhu, “Interstitial diffuse optical tomography using an adjoint model with linear sources,” Proc. SPIE 6846, 68450C (2008). [CrossRef]

], we concentrate on the use of NIRFAST [1

, 8

K. K.-H. Wang and T. C. Zhu, “Reconstruction of optical properties using a diffusion model for interestitial diffuse optical tomography,” Proc. SPIE 7164, 71640P (2009). [CrossRef]

] for solving the inverse problem. The NIRFAST [1

] has been modified to more accurately calculate light fluence rate per source strength for our CW reconstruction purpose [8

K. K.-H. Wang and T. C. Zhu, “Reconstruction of optical properties using a diffusion model for interestitial diffuse optical tomography,” Proc. SPIE 7164, 71640P (2009). [CrossRef]

]. The fluence rate data in forward calculation (ϕ_c ) is generated using FEM and based on the steady state light diffusion equation,

μ_{a} ϕ - \nabla \cdot (D \nabla ϕ) = S,

(1)

with a set of initial values of µ_a and diffusion coefficient D. In Eq. (1), S is the source term. Next, the Jacobian matrix J in CW scheme is established, which has the form

J = [\begin{array}{c} \frac{\partial \ln ϕ_{1}}{\partial D_{1}} & \dots & \frac{\partial \ln ϕ_{1}}{\partial D_{N}} & ; & \frac{\partial \ln ϕ_{1}}{\partial μ_{a 1}} & \dots & \frac{\partial \ln ϕ_{1}}{\partial μ_{a N}} \\ ⋮ & ⋱ & ⋮ & ; & ⋮ & ⋱ & ⋮ \\ \frac{\partial \ln ϕ_{M}}{\partial D_{1}} & \dots & \frac{\partial \ln ϕ_{M}}{\partial D_{N}} & ; & \frac{\partial \ln ϕ_{M}}{\partial μ_{a l}} & \dots & \frac{\partial \ln ϕ_{M}}{\partial μ_{a N}} \end{array}]

(2)

where M and N are the number of measured data and total mesh nodes, respectively. In our modification, absolute fluence rate per source strength, without renormalization, is used throughout the algorithm to enable separation of µ_a and µ_s ’.

In the next step, the projection error is calculated as the summation of the absolute difference between natural logarithm of ϕ_m and ϕ_c . If the stopping criteria are achieved, we record the reconstructed optical properties and terminate the program. Otherwise, the optical properties will be updated and iteratively used in the above calculation procedures. The stopping criteria are either the change of projection error less than 2 % or negative optical properties value generated during the iterative calculation. The method described in this section is named as DOT method.

Fig. 1. (A) 2D orientation of prostate, urethra, and rectum contours for one cut of the patient prostate. Cross and circle are the locations of point sources and isotropic detectors, respectively. The actual sources and detectors are numbered. (B) Extension of the 2D case (A) in a 3D geometry: The prostate contour in (A) is extended 1 cm in z-axis direction. Circles are the detector locations used in the 3D reconstruction calculation.

2.3 Point-by-point method

The optimization method used to extract optical properties based on the measurement of fluence rate at several distances from point source has been published elsewhere [6

, 9

A. Dimofte, J. C. Finlay, and T. C. Zhu, “A method for determination of the absorption and scattering properties interstitially in turbid media,” Phys. Med. Biol. 50, 2291–2311 (2005). [CrossRef] [PubMed]

]. The accuracy of this method has been established to be 8 % for µ_a and 18 % for µ_s ’ [9

]. Briefly, with the diffusion approximation, the light fluence rate ϕ at a distance r from a point source with source strength, S, can be expressed as

ϕ (r) ⁄ S = \frac{3 {μ'}_{s}}{4 π r} e^{- μ_{eff} r}

(3)

where µ_eff is effective attenuation coefficient and

μ_{eff} = \sqrt{3 μ_{a} {μ'}_{s}}

. For a point source measurement with a given source strength, µ_a and µ_s ’ can be determined separately by the slope of the spatial decay of light fluence rate and the magnitude of fluence rate near the source. If sufficient number of measurements are made in a prostate, a modified differential evolution algorithm developed by our group [9

] is able to fit the measured data using Eq. (3). A map of optical properties can therefore be made by fitting each source detector pair data, individually. This method is named as point-by-point method in this study.

3. Results and discussion

3.1 Reconstruction using 3D Mathematical phantom

A 3D mathematical phantom is used to verify our CW reconstruction algorithm developed for the iDOT system. To simulate the clinical environment, we place 5 point sources and 12 detectors within or close to the prostate (Fig. 1(A)). After the prostate contour has been extended to 1 cm in the z-axis direction as described in Sec. 2.2, the number of detectors used in the calculation is 132, and total simulated data point is 660. Fig. 2(A) shows two anomalies, spheres of radius 0.5 cm, and an outer medium are placed inside and surrounding the prostate. The prostate contour is the same as shown in Fig. 1(A), which is extracted from a patient. The outer boundary is used to avoid an artificial effect, a sudden drop of fluence rate on the prostate boundary if there is no outer medium existed [10

J. Li, T. C. Zhu, and J. C. Finlay, “Study of light fluence rate distribution in photodynamic therapy using finite-element method,” Proc. SPIE 6139, 61390M (2006). [CrossRef]

]. The optical properties of the background are set as µ_s ’=14 (cm^-1) and µ_a =0.3 (cm^-1) for the prostate and the surrounding medium, based on a published result [6

]. For the left anomaly, the µ_s ’ and µ_a are 14 (cm^-1) and 0.6 (cm-1), respectively, and for the right anomaly, µ_s ’ is 28 (cm^-1) and µ_a is 0.3 (cm^-1).

Fig. 2. (A) Schematics of the mathematical phantom for one slide of prostate contour. Two anomalies are inserted into the prostate, and an extra outer medium is applied to enclose the prostate. (B) Corresponding COMSOL-generated mesh used in the forward calculation for the reconstruction procedures.

The simulated data with 1 % noise is calculated on a finer mesh (not shown) using the known optical properties. To reduce the number of unknowns, a coarse mesh is used for the forward calculation in the reconstruction procedures (Fig. 2(B)). To further reduce the reconstructed unknowns during the optical properties updated procedure, a separated second reconstruction basis, pixel basis 12×12×12 provided by NIRFAST is chosen [1

The reconstructed results for the 3D mathematical phantom are shown in Fig. 3. The source and the center of anomalies are located at 0.5 cm at z-axis. Fig. 3 shows the known and reconstructed distribution of µ_a and µ_s ’, where the locations of anomalies are correctly reconstructed. The maximum reconstructed µ_a and µ_s’ are 0.5 and 27 (1/cm) vs. the true optical properties 0.6 and 28 (1/cm) (Fig. 3 (A1) vs. (B1) and Fig. 3 (A2) vs. (B2)). A cross talk between two anomalies is observed in Fig. 3 (B1) and (B2).

Fig. 3. Known distributions for (A1) µ_a and (A2) µ_s ’ at z=0.5 cm and the corresponding reconstructed result for the distribution of (B1) µ_a and (B2) µ_s ’ at the same plane for the 3D mathematical phantom.

3.2 Clinical data reconstruction

For a prostate patient, we have conducted a series of in-vivo measurement for 6 source locations located at z-axis, z=0.5 to 3 cm. The source and detector orientation follow the same scenario shown in Fig. 1(A). However, due to the experimental complexity, only the detectors close to the source recorded data. For example, for source 1 (Fig. 1(A)), only detector 2, 5 and 6 obtained the corresponding fluence rate profile. Therefore, for each plane, total measured data point is ranging from 149 to 160, compared to 660 used in the 3D phantom calculation. For DOT method, approximate 2600 to 4400 nodes and 14500 to 25000 tetrahedral elements are used in the forward model to simulate the fluence rate distribution for different cases while a coarser pixel base mesh (12×12×12) are used in the inverse model to calculate µ_a and µ_s ’. The total computational time for one cut is around 30 mins using a desktop with 2.4 GHz processor and 2 GB of RAM.

Figure 4 shows the comparison between the point-by-point and DOT method for the patient prostate at z=0.5 cm. Significant heterogeneity for µ_a and µ_s ’ are observed in images obtained by both methods (Fig. 4). The hot spot locations of optical properties reconstructed by DOT method is consistent with the results calculated by point-by-point method (Fig. 4(A) vs. (B) and Fig. 4(D) vs. (E)). The central hot-spot may be due to the bleeding within patient prostate. The bleeding is independently verified by examining the catheter after the treatment is completed. Based on our clinical experience of treating 12 patients, it is concluded that if the measured light fluence rate per source strength at 5 mm is less than 0.1 mW/(cm²mW), then bleeding occurs. However, we are not sure if the hot spot for scattering is also related to bleeding because of the extremely low fluence rate signal detected in these regions. The inverse algorithm, therefore, reconstructed high µ_a and µ_s ’ over there regions. We also compared these two methods along x-axis at z=0.5 cm and y=1.75 cm, blue dash lines in Fig. 4, for µ_a and µ_s ’ (Fig. 4(C) and (F)). The spatial-resolved profiles of µ_a are consistent for both methods (Fig. 4 (C)). The comparison for µ_s ’ shows large inconsistency at x=1.75 cm, but the other points are matched acceptably well (Fig. 4 (F)).

Figure 5(A) shows the measured (plus) and calculated (circle) fluence rate for each source and detector pair vs. virtual detector indices assigned for DOT method. The computed fluence rate is based on the reconstructed optical properties from Fig. 4(B) and (E) at z=0.5 cm plane. We notice that for detector indices from 23 to 33 and 78 to 83, there is a large discrepancy. This inconsistency corresponds to the data of source 1 and detector 6 pair, and source 3 and detector 7 pair (Fig. 5(B)) and is attributed to the difference of source and detector separations between reconstruction assumption and the actual measurement condition [9

]. Template with a fixed 0.5-cm grid is used for the position of sources and detectors during the reconstruction. The actual separation between the detector and point source can be different from this assumption. To test this argument, we assign extra detector 13 and 14 to replace the detector 6 and 7 with respect to source 1 and 3 (Fig. 5(B)). Detector 13 and 14 are 0.3 and 0.4 cm closer to the source than the detector 6 and 7. The results of recalculated fluence rate are shown in Fig. 5(C). Interestingly, we notice changing the separation distance not only improve the calculated fluence rate for a certain source and detector pair, but also induce a global optimization for other pairs (Fig. 5(C)).

Fig. 4. Results of in vivo map of µ_a using point-by-point (A) and DOT (B) methods for patient prostate at z=0.5 cm plane. (C) The corresponding comparison of µ_a between the two methods along x-axis (blue dash line in (A) and (B)) at z=0.5 and y=1.75 cm. The solid line and circle represent the results of the DOT and point-by-point method, respectively. (D), (E) and (F) are the corresponding plots for µ_S ’.

Fig. 5. (A) Forward calculated fluence rate using the reconstructed optical properties from Fig. 4 (B) and (E) for z=0.5 cm plane, compared to the measured data at different virtual detector (scanning position) locations. (B) Two extra detectors 13 and 14 (green triangle) replacing the detector 6 and 7 with respect to source 1 and 3. Detector 13 and 14 are 0.3 and 0.4 cm closer to the source than the detector 6 and 7. (C) is the comparison result after correcting the source and detector separation.

We have produced volumetric distributions of µ_a and µ_s ’ for the in-vivo patient prostate (Fig. 6 (A) and (B)), respectively. This hot spot is observed through prostate (Fig. 6 (A) and (B)). For either µ_a or µ_s ’ map, the heterogeneity is observed to penetrate several planes for the same patient. This result emphasizes the importance to explore the patient-specific optical properties for PDT treatment.

Fig. 6. Volumetric distributions of (A) µ_a and (B) µ_s’, respectively, for an in-vivo prostate. The prostate contours (red lines) at different z planes are shown in both cases.

Generally, homogenous optical properties are assumed in PDT treatment planning. The fluence rate will be largely affected by the heterogeneous optical properties (Figs. 4 and 6), as described in detail elsewhere [11

J. Li and T. C. Zhu, “Determination of in vivo light fluence distribution in a heterogeneous prostate during photodynamic therapy,” Phys. Med. Biol. 53, 2103–2114 (2008). [CrossRef] [PubMed]

]. The corresponding optimization algorithm to adjust light source parameters (i.e. strengths, lengths, and locations), based on the optical properties heterogeneity, is therefore necessary to facilitate the current PDT light dosimetry system [12

J. Li, M. D. Altschuler, S. M. Hahn, and T. C. Zhu, “Optimization of light source parameters in the photodynamic therapy of heterogeneous prostate,” Phys. Med. Biol. 53, 4107–4121 (2008). [CrossRef] [PubMed]

]. One good strategy for taking into account the inhomogeneity in optical properties without prior knowledge is to perform in-situ light dosimetry to measure the actual light fluence rate at area of interest to ensure sufficient light delivery time to achieve uniform total light fluence.

4. Conclusion

A 3D reconstruction model has been developed, which is able to recover volumetric µ_a and µ_s ’ of human prostate using the clinical data obtained from a CW-iDOT system. This model maximizes the usage of our scanning data; for each source location, we are able to use the measured data at, below and above the source plan. Moreover, the calculated fluence rate can match measured data using correct source and detector separation. In the other word, correctly predicting light fluence rate distribution within patient prostate based on the reconstructed optical properties can be a significant improvement to the prostate-PDT light dosimetry.

The ongoing direction in our group is to improve the accuracy of the reconstructed optical properties and anomalies location, we propose to implement a regional base reconstruction algorithm [13

Z. Jiang, D. Piao, G. Xu, J. W. Ritchey, G. R. Holyoak, K. E. Bartels, C. F. Bunting, G. Slobodov, and S. Kransinski, “Trans-rectal ultrasound-coupled near-infrared optical tomography of the prostate, Part II: Experimental demonstration,” Opt. Express 16, 17505–17520 (2008). [CrossRef] [PubMed]

, 14

G. Xu, D. Piao, C. H. Musgrove, C. F. Bunting, and H. Dehgani, “Trans-rectal ultrasound-coupled nearinfrared optical tomography of the prostate, Part I: Simulation,” Opt. Express 16, 17484–17504 (2008). [CrossRef] [PubMed]

] with prior information into the current routine. Other image modalities (e.g. MR) will be used to locate tumor within patient prostate, and the regional base routine will enforce a homogeneous optical property within the same zone regime.

Acknowledgments

The authors would like to thank Dr. Hamid Dehghani for helping with NIRFAST programming and both Dr. Brian Pogue and Dehghani for providing the NIRFAST software. We thank Dr. Xiaodong Zhou for helpful discussions. This work is supported by grants from National Institute of Health (NIH) R01 CA 109456 and P01 CA87971.

References and links

1.	H. Dehgani, M. E. Eames, P. K. Yalavarthy, S. C. Davis, S. Srinivasan, C. M. Carpenter, B. W. Pogue, and K. D. Paulsen, “Near infrared optical tomogrophy using NIRFAST: Algorithm for numerical model and image reconstruction,” Commun. Numer. Methods Eng. 25, 711–732 (2008). [CrossRef]
2.	A. Jemal, R. Siegel, E. Ward, Y. Hao, J. Xu, T. Murray, and M. J. Thun, “Cancer statistics, 2008,” Ca-Cancer J. Clin. 58, 71–96 (2008). [CrossRef] [PubMed]
3.	J. Axelsson, J. Swartling, and S. Andersson-Engels, “In vivo photosensitizer tomography inside the human prostate,” Opt. Lett. 34, 232–234 (2009). [CrossRef] [PubMed]
4.	N. Pantong, J. Su, H. Shan, M. V. Klibanov, and H. Liu, “Globally accelerated reconstruction algorithm for diffusion tomography with continuous-wave source in an arbitrary convex shape domain,” J. Opt. Soc. Am. A. 26, 456–472 (2009). [CrossRef]
5.	T. C. Zhu, J. C. Finlay, and S. M. Hahn, “Determination of the distribution of light, optical properties, drug concentration, and tissue oxygenation in-vivo in human prostate during motexafin lutetium-mediated photodynamic therapy,” J. Photochem. Photobiol. B 79, 231–241 (2005). [CrossRef] [PubMed]
6.	T. C. Zhu, A. Dimofte, J. C. Finlay, D. Stripp, T. Busch, J. Miles, R. Whittington, S. B. Malkowicz, Z. Tochner, E. Glatstein, and S. M. Hahn, “Optical properties of human prostate at 732 nm measured in mediated photodynamic therapy,” Photochem. Photobiol. 81, 96–105 (2005). [CrossRef]
7.	X. Zhou and T. C. Zhu, “Interstitial diffuse optical tomography using an adjoint model with linear sources,” Proc. SPIE 6846, 68450C (2008). [CrossRef]
8.	K. K.-H. Wang and T. C. Zhu, “Reconstruction of optical properties using a diffusion model for interestitial diffuse optical tomography,” Proc. SPIE 7164, 71640P (2009). [CrossRef]
9.	A. Dimofte, J. C. Finlay, and T. C. Zhu, “A method for determination of the absorption and scattering properties interstitially in turbid media,” Phys. Med. Biol. 50, 2291–2311 (2005). [CrossRef] [PubMed]
10.	J. Li, T. C. Zhu, and J. C. Finlay, “Study of light fluence rate distribution in photodynamic therapy using finite-element method,” Proc. SPIE 6139, 61390M (2006). [CrossRef]
11.	J. Li and T. C. Zhu, “Determination of in vivo light fluence distribution in a heterogeneous prostate during photodynamic therapy,” Phys. Med. Biol. 53, 2103–2114 (2008). [CrossRef] [PubMed]
12.	J. Li, M. D. Altschuler, S. M. Hahn, and T. C. Zhu, “Optimization of light source parameters in the photodynamic therapy of heterogeneous prostate,” Phys. Med. Biol. 53, 4107–4121 (2008). [CrossRef] [PubMed]
13.	Z. Jiang, D. Piao, G. Xu, J. W. Ritchey, G. R. Holyoak, K. E. Bartels, C. F. Bunting, G. Slobodov, and S. Kransinski, “Trans-rectal ultrasound-coupled near-infrared optical tomography of the prostate, Part II: Experimental demonstration,” Opt. Express 16, 17505–17520 (2008). [CrossRef] [PubMed]
14.	G. Xu, D. Piao, C. H. Musgrove, C. F. Bunting, and H. Dehgani, “Trans-rectal ultrasound-coupled nearinfrared optical tomography of the prostate, Part I: Simulation,” Opt. Express 16, 17484–17504 (2008). [CrossRef] [PubMed]

OCIS Codes
(170.3880) Medical optics and biotechnology : Medical and biological imaging
(170.5180) Medical optics and biotechnology : Photodynamic therapy
(170.6960) Medical optics and biotechnology : Tomography
(170.7230) Medical optics and biotechnology : Urology

ToC Category:
Medical Optics and Biotechnology

History
Original Manuscript: May 4, 2009
Revised Manuscript: June 9, 2009
Manuscript Accepted: June 22, 2009
Published: June 26, 2009

Virtual Issues
Vol. 4, Iss. 9 Virtual Journal for Biomedical Optics

Citation
Ken Kang-Hsin Wang and Timothy C. Zhu, "Reconstruction of in-vivo optical properties for human prostate using interstitial diffuse optical tomography," Opt. Express 17, 11665-11672 (2009)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-14-11665

Sort: Author | Year | Journal | Reset

References

H. Dehgani, M. E. Eames, P. K. Yalavarthy, S. C. Davis, S. Srinivasan, C. M. Carpenter, B. W. Pogue, and K. D. Paulsen, "Near infrared optical tomogrophy using NIRFAST: Algorithm for numerical model and image reconstruction," Commun. Numer. Methods Eng. 25, 711-732 (2008). [CrossRef]
A. Jemal, R. Siegel, E. Ward, Y. Hao, J. Xu, T. Murray, and M. J. Thun, "Cancer statistics, 2008," Ca-Cancer J. Clin. 58, 71-96 (2008). [CrossRef] [PubMed]
J. Axelsson, J. Swartling, and S. Andersson-Engels, "In vivo photosensitizer tomography inside the human prostate," Opt. Lett. 34, 232-234 (2009). [CrossRef] [PubMed]
N. Pantong, J. Su, H. Shan, M. V. Klibanov, and H. Liu, "Globally accelerated reconstruction algorithm for diffusion tomography with continuous-wave source in an arbitrary convex shape domain," J. Opt. Soc. Am. A. 26, 456-472 (2009). [CrossRef]
T. C. Zhu, J. C. Finlay, and S. M. Hahn, "Determination of the distribution of light, optical properties, drug concentration, and tissue oxygenation in-vivo in human prostate during motexafin lutetium-mediated photodynamic therapy," J. Photochem. Photobiol. B 79, 231-241 (2005). [CrossRef] [PubMed]
T. C. Zhu, A. Dimofte, J. C. Finlay, D. Stripp, T. Busch, J. Miles, R. Whittington, S. B. Malkowicz, Z. Tochner, E. Glatstein, and S. M. Hahn, "Optical properties of human prostate at 732 nm measured in mediated photodynamic therapy," Photochem. Photobiol. 81, 96-105 (2005). [CrossRef]
X. Zhou, and T. C. Zhu, "Interstitial diffuse optical tomography using an adjoint model with linear sources," Proc. SPIE 6846, 68450C (2008). [CrossRef]
K. K.-H. Wang, and T. C. Zhu, "Reconstruction of optical properties using a diffusion model for interestitial diffuse optical tomography," Proc. SPIE 7164, 71640P (2009). [CrossRef]
A. Dimofte, J. C. Finlay, and T. C. Zhu, "A method for determination of the absorption and scattering properties interstitially in turbid media," Phys. Med. Biol. 50, 2291-2311 (2005). [CrossRef] [PubMed]
J. Li, T. C. Zhu, and J. C. Finlay, "Study of light fluence rate distribution in photodynamic therapy using finite-element method," Proc. SPIE 6139, 61390M (2006). [CrossRef]
J. Li, and T. C. Zhu, "Determination of in vivo light fluence distribution in a heterogeneous prostate during photodynamic therapy," Phys. Med. Biol. 53, 2103-2114 (2008). [CrossRef] [PubMed]
J. Li, M. D. Altschuler, S. M. Hahn, and T. C. Zhu, "Optimization of light source parameters in the photodynamic therapy of heterogeneous prostate," Phys. Med. Biol. 53, 4107-4121 (2008). [CrossRef] [PubMed]
Z. Jiang, D. Piao, G. Xu, J. W. Ritchey, G. R. Holyoak, K. E. Bartels, C. F. Bunting, G. Slobodov, and S. Kransinski, "Trans-rectal ultrasound-coupled near-infrared optical tomography of the prostate, Part II: Experimental demonstration," Opt. Express 16, 17505-17520 (2008). [CrossRef] [PubMed]
G. Xu, D. Piao, C. H. Musgrove, C. F. Bunting, and H. Dehgani, "Trans-rectal ultrasound-coupled near-infrared optical tomography of the prostate, Part I: Simulation," Opt. Express 16, 17484-17504 (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.

Click here to see a list of articles that cite this paper