# Improving the accuracy of a solid spherical source radius and depth estimation using the diffusion equation in fluorescence reflectance mode

- Marjaneh Hejazi†
^{1, 2}, - Florian Stuker†
^{3}, - Divya Vats
^{3}and - Markus Rudin
^{3, 4}Email author

**9**:28

https://doi.org/10.1186/1475-925X-9-28

© Hejazi et al; licensee BioMed Central Ltd. 2010

**Received: **19 October 2009

**Accepted: **19 June 2010

**Published: **19 June 2010

## Abstract

### Background

Non-invasive planar fluorescence reflectance imaging (FRI) is used for accessing physiological and molecular processes in biological tissue. This method is efficiently used to detect superficial fluorescent inclusions. FRI is based on recording the spatial radiance distribution (SRD) at the surface of a sample. SRD provides information for measuring structural parameters of a fluorescent source (such as radius and depth). The aim of this article is to estimate the depth and radius of the source distribution from SRD, measured at the sample surface. For this reason, a theoretical expression for the SRD at the surface of a turbid sample arising from a spherical light source embedded in the sample, was derived using a steady-state solution of the diffusion equation with an appropriate boundary condition.

### Methods

The SRD was approximated by solving the diffusion equation in an infinite homogeneous medium with solid spherical sources in cylindrical geometry. Theoretical predications were verified by experiments with fluorescent sources of radius 2-6 mm embedded at depths of 2-4 mm in a tissue-like phantom.

### Results

The experimental data were compared with the theoretical values which shows that the root mean square (RMS) error in depth measurement for nominal depth values d = 2, 2.5, 3, 3.5, 4 mm amounted to 17%, 5%, 2%, 1% and 5% respectively. Therefore, the average error in depth estimation was ≤ 4% for depths larger than the photon mean free path.

### Conclusions

An algorithm is proposed that allows estimation of the location and radius of a spherical source in a homogeneous tissue-like phantom by accounting for anisotropic light scattering effect using FRI modality. Surface SRD measurement enabled accurate estimates of fluorescent depth and radius in FRI modality, and can be used as an element of a more general tomography reconstruction algorithm.

## Background

Optical imaging methods such as fluorescence reflectance imaging (FRI), bioluminescent imaging and fluorescent tomography are used both in the areas of clinical research, such as dermatological and intra-operative imaging and in pre-clinical research for small animal imaging [1]. A conventional FRI system records the spatial distribution of the scattered fluorescent light at the surface of a sample. This light distribution contains sufficient information for determining structural parameters such as radius and depth of fluorescence probe within biological tissue assuming a distribution of spherical sources. Information on depth and radius is required for deriving quantitative information on the source intensity, which may be translated into information of molecular and cellular processes within the tissue.

Recently, methods of the depth reconstruction for fluorescence molecular tomography have been widely studied [2–4]. In transmission mode, D'Andrea et al. [5] reported an analytical algorithm to recover the coordinates of a fluorescent inclusion embedded in a turbid medium with accuracy better than 2 mm. Kuo et al.[6] developed a three dimensional reconstruction method for localizing and determining the photon flux of the sources only from a single view. The experimental evaluation revealed that the average error in depth reconstruction for a point source in a sample tissue was 3% for depth larger than 7 mm and exceeded 50% for the superficial depths (≤ 2.5 mm) [7]. The lack of reconstruction accuracy was attributed to limitations in the assumption that photons were scattered isotropically in biological tissue. Eisdath et al. [8] proposed a theoretical model of photon migration in reflectance mode for localizing the fluorescent inclusion deeply embedded in a turbid media. The source coordinates were obtained with an average error of 4% at 3 mm and less than 10% for depths larger than 6 mm. Up to our knowledge, the radius estimation for superficial fluorescent inclusions in tissue like phantoms has not been reported in FRI modality to date.

The objective of this article is to improve the accuracy of the estimation of source depth and radius in reflection mode by accounting for anisotropic photon scattering. The proposed algorithm can also be regarded as important part of tomographic reconstruction algorithms.

## Methods

### Theoretical Background

The intensity distribution at the surface of a turbid medium is described by a spatial radiance distribution (SRD). This SRD was estimated by solving the diffusion equation in an infinite homogeneous medium with solid spherical sources of radius *R* with its center located at depth *z*.

*μ*

_{ a }and a reduced scattering coefficient , both with values corresponding to biological tissue.

where Φ(**r**) is the photon fluence (photons/cm^{2}) at location **r** with the radius given by
for *z* ≤ *d*, Ω describes the region of interest, *S*(**r**) an isotropic source term, while the diffusion coefficient *D*(**r**) = 3(*μ*
_{
a
}(**r**) + *μ*
_{
s
}'(**r**)), accounts for diffusive photon propagation. This equation is derived from the adiative transfer equation [11] by applying the diffusion approximation and is valid for highly scattering media of a thickness greater than the photon mean-free path length.

*R*embedded in an infinite homogenous medium has been formulated as [12]:

where *P* is the source power
denotes the effective attenuation coefficient, while *α*(*R*) stands for the exponential terms given in brackets.

**r**,

*R*) at the interface between two media with different refractive indices can be approximated by an index-mismatched Robin-type boundary condition. This condition for Φ(

**r**,

*R*) on ∂Ω is described in the following formula [13, 14]:

**n**is the unit vector normal to a surface element ∂Ω and

*A*is a parameter governing the internal reflection at the boundary. The value of

*A*= (1+R)/(1-R) depends on the relative refractive index mismatch between tissue and air and can be derived from Fresnel's law, with R as [14]:

*n*=

*n*

_{ a }/

*n*

_{ t },

*n*

_{ a }and

*n*

_{ t }being the refractive indices for air and tissue, respectively. This equation yields

*A*= 2.51 for the air-phantom interface with typical value of

*n*= 1.33 [15]. Based on the boundary condition the measured fluence rate on ∂Ω is [13, 16].

*L*(

**r**,

*R*, ) is introduced which accounts for anisotropic optical properties of the tissue at

**r**in direction of , i.e the expression for the radiance (Eq. (6)) is expanded according to [17]:

**F**(

**r**) = -

*D*(

**r**)▽Φ(

**r**,

*R*) at a distance r from the sphere in direction [16]. Combining the various equations one obtains:

which for **r** ∈ *S* corresponds to the radiance distribution at the surface *S* of the sample that is generated by a spherical light source of radius *R*, i.e to the SRD.

### Experimental Evaluation

*λ*= 500-950 nm, 30 W), a standard excitation and emission filter set and in addition a liquid crystal tunable filter (LCTF). The LCTF transmits a narrow bandwidth (10 nm) of the emitted fluorescence light to a 1.3 megapixel CCD camera. A 615-665 nm excitation filter and a 700 nm long-pass emission filter were used for the evaluation.

The phantom was made from 6 g agarose (BioGene, Kimbolton, UK), 24 ml Intralipid 20% (Fresenius SE, Bad Homburg, Germany) and 18 μl Indian ink (Pelikan Holding, Schindellegi, Switzerland) [18] dissolved in 600 ml water. The optical properties of the phantom (with radius 10 cm) were measured with a frequency-domain tissue oximeter (oxiplexTS, ISS, Champaign, USA). The reduced scattering and absorption coefficients, at 692 nm, were then found to be
= 8 cm^{-1} and μ_{a} = 0.11 cm^{-1} respectively. The same mixture was poured into four plexiglass molds of 6 × 6 cm and thickness ranging from 2 mm to 4 mm. The width of the mould was more than 10 penetration depths (0.6 cm) which fulfilled the requirement for a semi-infinite boundary medium [19]. After solidification, the agarose layers were removed from the molds and positioned by a fixed holder on the top of solid spherical source. The sources with radius 2-6 mm contained 5 mg/l near-infrared quantum dot (Qdot 705 ITK, Invitrogen, Basel, Switzerland) with emission wavelength at 705 nm.

The optical images were recorded for varying thickness *d* and denoised by a bandpass filter [20]. The bandpass filter consists of a Gaussian low-pass filter and a boxcar kernel which was used as a high-pass filter. First, a low pass filtered image was generated by convolving the diffuse image with the Gaussian filter, a high-pass filtered image was then obtained by convolving the diffuse image with the boxcar function. The difference between low-pass filtered image and high-pass filtered image allowed extraction of the final image for quantitative analysis.

## Results

**r**, R)(Eq. (5)) as function of d and radius (Fig. 4). Source radii range from 2 mm to 6 mm, while d was varied between 2 mm and 4 mm. In general, good agreement between experimental and theoretical FWHM values has been found. Error estimation revealed that the RMS of depth measurement errors for depth values d = 2, 2.5, 3, 3.5, 4 mm were 17%, 5%, 2%, 1% and 5% respectively, for the source radii 2-6 mm covered in our experiment (Fig. 5).

## Discussion and Conclusions

In this work, we have shown that the structural parameters of a spherical source (depth and radius) can be estimated by measuring the spatial radiance distribution (SRD) at the sample surface. The SRD at the surface of a turbid medium was derived by calculating a steady-state solution of the diffusion equation using the Robin boundary condition.

The accuracy of the theoretical model was validated by fluorescence reflectance imaging of a hemispace homogeneous phantom containing fluorescent source. The comparison showed the RMS of depth measurement errors for depth values d = 2, 2.5, 3, 3.5, 4 mm were 17%, 5%, 2%, 1% and 5% respectively. As discussed in previous studies, the relatively large (17%) overestimation at a depth of 2 mm was attributed to the inaccuracy of the diffusion approximation for distances comparable to the mean free path (1/(
+ μ_{a}) of the photon propagation in a turbid medium [15]. For depth values exceeding 4 mm the decrease in fluorescence intensity at the sample surface became a limiting factor; as a consequence of decreased signal-to-noise ratios extraction of accurate SRDs using the denoising method described became difficult and significant errors in FWHM estimation arose. Comparing our method with the procedure described earlier [6] we obtained improved depth estimates close to the sample surface (5% versus 45% average error). This might indicate the importance of accounting for the anisotropic distribution of the scattered light in highly turbid media when deriving the SRD, in particular when analyzing source locations within a few millimeter from the sample surface. Eisdath et al. [8] have developed a model for reflectance mode imaging for localizing a fluorescent source in three dimensions. Their results showed that the average error in source localization was less than 10% for depths larger than 6 mm. However, when using a fluorescence reflectance system SNR becomes a limiting factor for deriving accurate information on source position and dimension within a turbid medium for depth values large than 5 mm.

The algorithm in its current form is not suited for the accurate determination of source location and radius in irregular biological tissues, in particular when considering reconstruction from a single view. The major problem when translating the model to an irregular and heterogeneous sample is that the solution of the inverse problem for the fluorescent source structural parameter does not yield a unique SRD from measurements at a single wavelength [4, 6]. Next steps in the model development are therefore extension to multiple wavelengths, which will improve the accuracy of the reconstructed source distributions [4, 21]. A multi-wavelength algorithm can recover the source depth and radius due to the wavelength dependence of SRD (Eqs. (3) and (5)).

In conclusion, an algorithm is proposed that allows estimation of the location and radius of a spherical source in a homogeneous tissue-like phantom by accounting for anisotropic light scattering effect using FRI modality. Experimental validation revealed the accuracy and limitations of the current model. It is obvious that heterogeneity of optical parameters as found in biological tissues will affect the accuracy of the source reconstruction algorithm. These effects will be analyzed in future work.

## Notes

## Declarations

### Acknowledgements

The authors would like to thank the Swiss National Science Foundation (SNF 310000-112835) for generous funding.

## Authors’ Affiliations

## References

- Weissleder R, Ntziachristos V:
**Shedding light onto live molecular targets.***Nat Med*2003,**9:**123–128. 10.1038/nm0103-123View ArticleGoogle Scholar - Ntziachristos V, Tung C, Bremer C, Weissleder R:
**Fluorescence molecular tomography resolves protease activity in vivo.***Nat Med*2002,**8:**757–760. 10.1038/nm729View ArticleGoogle Scholar - Ntziachristos V, Weissleder R:
**Experimental three-dimensional fluorescence reconstruction of diffuse media by use of a normalized Born approximation.***Opt Lett*2001,**26:**893–895. 10.1364/OL.26.000893View ArticleGoogle Scholar - Dehghani H, Davis SC, Jiang SD, Pogue BW, Paulsen KD, Patterson MS:
**Spectrally resolved bioluminescence optical tomography.***Opt Lett*2006,**31:**365–367. 10.1364/OL.31.000365View ArticleGoogle Scholar - D'Andrea C, Spinelli L, Comelli D, Valentini G, Cubeddu R:
**Localization and quantification of fluorescent inclusions embedded in a turbid medium.***Phys Med Biol*2005,**50:**2313–2327. 10.1088/0031-9155/50/10/009View ArticleGoogle Scholar - Kuo C, Coquoz O, Troy TL, Xu H, Rice BW:
**Three-dimensional reconstruction of in vivo bioluminescent sources based on multispectral imaging.***J Biomed Opt*2007,**12:**1–12. 024007 10.1117/1.2717898View ArticleGoogle Scholar - Virostko J, Powers AC, Jansen DE:
**Validation of luminescent source reconstruction using single-view spectrally resolved bioluminescence images.***Appl Opt*2007,**46:**2540–2547. 10.1364/AO.46.002540View ArticleGoogle Scholar - Eidsath A, Chernomordik V, Gandjbakhche A, Smith P, Russo A:
**Three-dimensional localization of fluorescent masses deeply embedded in tissue.***Phys Med Biol*2002,**47:**4079–4092. 10.1088/0031-9155/47/22/311View ArticleGoogle Scholar - Schultz R, Ripoll J, Ntziachristos V:
**Experimental fluorescence tomography of tissues with noncontact measurements.***IEEE Trans Med Imag*2004,**23:**492–500. 10.1109/TMI.2004.825633View ArticleGoogle Scholar - Arridge SR, Schweiger M, Hiraoka M, Delpy DT:
**A finite element approach for modeling photon transport in tissue.***Med Phys*1993,**20:**299–309. 10.1118/1.597069View ArticleGoogle Scholar - Wang LV:
*Biomedical optics: Principles and imaging.*John Wiley and Sons; 2007.Google Scholar - Shen H, Cong W, Qian X, Durairaj K, Wang G:
**Numerical study on the validity of the diffusion approximation for computational optical biopsy.***J Opt Soc Am A*2007,**24:**423–429. 10.1364/JOSAA.24.000423View ArticleGoogle Scholar - Cong WX, Wang G, Kumar D, Liu Y, Jiang M, Wang LH, Hoffman EA, McLennan G, McCray PB, Zabner J, Cong A:
**A practical reconstruction method for bioluminescence tomography.***Opt Express*2005,**13:**6756–6771. 10.1364/OPEX.13.006756View ArticleGoogle Scholar - Schweiger M, Arridge SR, Hiraoka M, Delpy DT:
**The finite element method for the propagation of light in scattering media: Boundary and source conditions.***Med Phys*1995,**22:**1779–1792. 10.1118/1.597634View ArticleGoogle Scholar - Hawrysz J, Roy R, Sevick-Muraca EM:
**Influence of the refractive index-mismatch at the boundaries measured in fluorescenceenhanced frequency-domain photon migration imaging Anuradha Godavarty.***Opt Express*2002,**10:**653–662.View ArticleGoogle Scholar - Comsa DC, Farrell TJ, Patterson MS:
**Quantification of bioluminescence images of point source objects using diffusion theory models.***Phys Med Biol*2006,**51:**3733–3746. 10.1088/0031-9155/51/15/009View ArticleGoogle Scholar - Welch AJ, Gemert MJCV:
*Optical thermal response of laser-irradiated tissue.*Plenum Press, New York; 1995.View ArticleGoogle Scholar - Cubeddu R, Pifferi A, Taroni P, Torricelli A, Valentini G:
**A solid tissue phantom for photon migration studies.***Phys Med Biol*1997,**42:**1971–1979. 10.1088/0031-9155/42/10/011View ArticleGoogle Scholar - Lin SP, Wang L, Jacques SL, Tittel FK:
**Measurement of tissue optical properties by the use of oblique-incidence optical fiber reflectometry.***Appl Opt*1997,**36:**136–143. 10.1364/AO.36.000136View ArticleGoogle Scholar - Crocker JC, Grier DG:
**:Methods of digital video microscopy for colloidal studies.***J Colloid Interface Sci*1996,**179:**298–310. 10.1006/jcis.1996.0217View ArticleGoogle Scholar - Wang G, Li Y, Jiang M:
**Uniqueness theorems in bioluminescence tomography.***Med Phys*2004,**31:**2289–299. 10.1118/1.1766420View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.