PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Med Phys. Author manuscript; available in PMC 2008 October 7.
Published in final edited form as:
Med Phys. 2008 July; 35(7): 3354–3363.
doi:  10.1118/1.2940603
PMCID: PMC2562618
NIHMSID: NIHMS57088

Three-dimensional fluorescence-enhanced optical tomography using a hand-held probe based imaging system

Abstract

Hand-held based optical imaging systems are a recent development towards diagnostic imaging of breast cancer. To date, all the hand-held based optical imagers are used to perform only surface mapping and target localization, but are not capable of demonstrating tomographic imaging. Herein, a novel hand-held probe based optical imager is developed towards three-dimensional (3-D) optical tomography studies. The unique features of this optical imager, which primarily consists of a hand-held probe and an intensified charge coupled device (ICCD) detector, are its ability to; (i) image large tissue areas (5×10 sq. cm) in a single scan, (ii) perform simultaneous multiple point illumination and collection, thus reducing the overall imaging time; and (iii) adapt to varying tissue curvatures, from a flexible probe head design. Experimental studies are performed in the frequency-domain on large slab phantoms (~ 650 ml) using fluorescence target(s) under perfect uptake (1:0) contrast ratios, and varying target depths (1-2 cm) and X-Y locations The effect of implementing simultaneous over sequential multiple point illumination towards 3-D tomography is experimentally demonstrated. The feasibility of 3-D optical tomography studies has been demonstrated for the first time using a hand-held based optical imager. Preliminary fluorescence-enhanced optical tomography studies are able to reconstruct 0.45 ml target(s) located at different target depths (1-2 cm). However, the depth recovery was limited as the actual target depth increased, since only reflectance measurements were acquired. Extensive tomography studies are currently carried out to determine the resolution and performance limits of the imager on flat and curved phantoms.

I. INTRODUCTION

Diagnostic imaging of early-stage breast cancer is essential for decreasing the cancer death rate in women in the United States. The conventional anatomical based screening techniques are neither comprehensive nor infallible, especially in women with dense breast tissues. Over the past 20 years, near-infrared (NIR) optical imaging approaches have been developed for breast cancer diagnosis, based upon the endogenous absorption contrast owing to the non-specific process of angiogenesis1-10 NIR light between 700-900 nm wavelengths is minimally absorbed and preferentially scattered, allowing its propagation through deep tissues. In recent years, external contrast agents are used to enhance the optical contrast between the diseased and normal tissue regions (i.e. fluorescence-enhanced optical imaging) and thus molecularly target metastatic cancer cells within the breast tissue.11-15

NIR optical imaging (with our without fluorescence) employs different source-detector imaging configurations during actual experimental studies. These imaging configurations can be broadly classified into10,16: (i) compressed tissue-based configuration, (ii) circular configuration, and (iii) sub-surface configuration. To date, most three-dimensional (3-D) optical imaging studies towards breast cancer diagnosis are restricted either to compressed tissue configuration2,17-21 or circular configuration22-36. The compressed tissue configuration is analogous to x-ray mammography, and is disadvantageous due to minimal patient comfort from tissue compression and limited information obtained around the entire breast tissue volume. The circular configuration has minimal patient discomfort, but is limited by the bulky and non-portable instrumentation. Sub-surface configuration is a relatively new method that requires notissue compression, and can be designed to mimic a portable and flexible imaging probe.32,37-50 However, 3-D tomography studies using the sub-surface imaging configuration are limited32,49-50 and challenging due to limited depth information obtained using only reflectance-based measurements during imaging reconstructions.

In recent years, hand-held based optical imaging systems, employing the sub-surface imaging configuration are developed in an attempt to translate the technology to the clinic, with maximum patient comfort and portability (against the bulky optical imagers available). However, the hand-held optical imagers available to date37-48 are not capable of: (i) performing 2-D or 3-D optical tomography, since the source and detector points are not co-registered onto the tissue contours, thus limiting them to only surface imaging and 2-D target localization; and (ii) contouring along tissue curvatures with good surface contact, since the hand-held probes have a flat surface..

Herein, a novel hand-held probe based optical imaging system is developed and 3-D optical tomography using fluorescence contrast agents has been demonstrated under different experimental conditions. A frequency-domain intensified charge-coupled device (ICCD) based detection system was coupled to a unique hand-held optical probe, such that simultaneous multiple point illumination and collection is possible over large areas (5×10 cm2). The effect of simultaneous multiple point illumination measurement geometry over sequential single point illumination geometry, in terms of total area of detected fluorescence amplitude and overall signal strength, has been demonstrated experimentally. Three-dimensional (fluorescence-enhanced) optical tomography has been successfully performed for the first time using a hand-held based optical imager, accounting for the excitation leakage issues as well.

II. MATERIALS AND METHODS

A. Instrumentation: The hand-held probe based optical imager

1. Hand-held probe

A portable and flexible hand-held optical probe was recently developed in our Optical Imaging Laboratory.51 The unique features of this probe are its ability to: (i) perform simultaneous multiple point illumination and collection, thus decreasing the total data acquisition time; and (ii) adapt to varying tissue curvatures, thus obtaining both trans-illumination and reflectance measurements. The major components of the hand-held probe are the 5× 10cm2 probe head constituting the illuminating/collecting optical fibers, a source fiber bundle, and a collection fiber bundle (Romack Inc., Williamsburg, VA) (See Fig. 1). The probe head consists of 6 illumination points and 165 collection points, with all points distributed 0.5 cm apart (see Fig. 1). The appropriate number of illuminating & collecting fibers for the given probe head area, along with their appropriate layout were determined from homogeneous fluorescence simulation studies performed on slab phantoms.51 Herein, the hand held optical probe is coupled to the source end (laser diode) and detector end (intensified charge coupled device, ICCD) via the source fiber bundle and collection fiber bundle, respectively.

FIG. 1
The hand-held probe consisting of a probe head, a source fiber bundle, and a collection fiber bundle. The solid (red) dots on the probe represent the six illumination points.

2. Imaging system

A frequency-domain based imaging system is developed for the current studies (see Fig.2). A modulated (100 MHz) laser light source is generated using a 785nm-wavelength laser diode (80mw, Sanyo DL7140-201S, Thorlabs Inc., Newton, NJ) and a frequency synthesizer (HP 8656B, Hayward, CA). The AC and DC power levels for the laser source are chosen such that the modulation depth (AC/DC) at the source end is maintained at 95%. In our case, the DC current of laser light is set to 217 mA, and the RF signal (i.e. AC signal) is set to +29.7 dBM. The RF signal is provided by a +27 dBM RF amplifier (ENI500L, ENI Inc., Rochester, NY) and the frequency synthesizer (-127 dBM to + 13 dBM range ), with power set to +2.7 dBM. The modulation source light is connected to source fiber bundle of the hand-held probe via a collimator-diffuser package, such that the single modulated laser source output can be split into six equal-intensity simultaneous modulated point sources.51

FIG. 2
(a) Schematic of the instrumentation for the hand-held probe based frequency domain ICCD optical imager. The hand-held probe is coupled to the laser source and ICCD detector for illumination and detection, respectively. (b) Experimental set-up of the ...

A gain modulated ICCD detection system is developed in order to acquire optical signals from various point locations simultaneously. The ICCD detector consists of an image intensifier (FS9910, ITT Night Vision, VA) that is optically coupled to a 16-bit frame transfer CCD camera (PI-SCX 7495-0002, Roper Scientific, Trenton, NJ). The detector is modulated at the image intensifier end (at its photocathode), using a second frequency synthesizer (PTS 310, Programmed Test Sources Inc., Littleton, MA) that is phase locked with the frequency synthesizer at the source end. The DC (-24 V) and AC values (+ 53 dBM) at the image intensifier are chosen such that the modulation depth is ~90% at the detector end. The AC (or RF) signal is provided by a +40 dBM RF amplifier (ENI503L, ENI Inc., Rochester, NY) and the PTS frequency synthesizer (0 dBM to + 13 dBM range), with power set to +13 dBM. The ICCD detector acquires the modulated optical signal from the tissue surface via the collection fiber bundle. The fiber bundle is coupled to the ICCD camera using a focusing lens (50mm, f/1.8D AF model, Nikon Inc., Melville, NY). The weak fluorescence optical signals are filtered from the attenuated, yet strong excitation signal via interference filters.). In the current study, two 830 nm interference filters of blocking OD 4 and 6 (HRF-830.0, blocking OD 6 & F-830.0, blocking OD 4, CVI Laser, Albuquerque, NM) are used together, in order to eliminate the excitation leakage and collect only the fluorescence signals. Preliminary experimental studies on non-fluorescing homogenous phantoms showed that the use of two interference filters (OD6 and OD4) reduced the excitation leakage by 2 orders of magnitude over using a single interference (OD6) filter.

During our instrument assessment studies, it was observed that maximum recovered modulation depth in homogeneous fluorescence phantoms was ~ 90%. This recovered modulation depth was close to the modulation depth at the source and detector ends (95% and 90%, respectively), which indicates that the quality of the recovered frequency-domain signal is high. In other words, the high modulated fluorescence signal is an indicator of the effectiveness of the imaging system to operate as a frequency-domain imaging system and recover the modulated (yet weak) fluorescence signals effectively.

3. Data acquisition

Homodyne based imaging is performed in order to acquire frequency domain measurements, wherein the modulation frequency at the source and detector end are maintained constant. A phase delay is introduced between the two phase locked frequency synthesizers.33,52 Steady-state phase-sensitive images are acquired by the CCD camera, at phase delays (here 32) varying between 0 to 2π. For all experiments, the exposure time of the CCD camera is set to 0.2 sec, and the acquired image of 1024×1024 pixels is hardware binned to 128×128 pixels, in order to improve signal to noise ratio, and also shorten the image read out time. The amplitude and phase (AC & θ) information for each of the 128×128 pixel is extracted using Fast Fourier Transforms (FFT), followed by extracting the AC & θ values for the 165 collection point locations on the hand-held probe head.

A referencing technique is applied to the fluorescence optical measurements (AC & θ), in order to overcome the various instrument effects. The amplitude and phase measurement at each point location (i) are referenced with respect to the point location (ref) with maximum fluorescence amplitude, in each experimental case.

equation M1
(1)

where ACR represents AC ratio (i.e. referenced AC), and RPS represents relative phase shift (i.e. referenced phase shift).

A subtraction-based post-processing technique (or “subtraction technique”) was employed on acquired fluorescence measurements to further eliminate excitation leakage, apart from using two interference filters. Initially, optical measurements (using both interference filters) were obtained from non-fluorescence homogenous phantom, wherein the signals detected are a representation of the excitation leakage (or background noise) . By subtracting this background noise from fluorescence measurements of any given phantom study, the excitation leakage can be removed to a greater extent. .

B. Experimental studies

Experimental studies were performed on heterogeneous phantoms to demonstrate the feasibility of 3-D tomography using the hand-held based imaging system. The heterogeneous phantom studies employed a single fluorescence target under perfect uptake conditions (i.e. no fluorescence in the background). In all these phantom studies, 1% Liposyn solution (Liposyn II, 20%, Henry Schein, Melville, NY) was used to mimic the human breast tissue for its similarities in optical scattering coefficient.53 The 1% Liposyn solution was filled into a transparent acrylic container (10×10×10 cm3) up to a volume of 650 ml. A hollow clear plastic sphere (volume of 0.45cm3) was filled with 1μM Indocyanine green dye (ICG, fluorescing contrast agent) in 1% Liposyn solution, and suspended in the slab phantom to mimic a target (i.e. a tumor) as shown in Fig. 2(b). Sodium polyaspartate (MW 3000-8000) (Sigma-Aldrich Chemical Co., St. Louis, CO) was used as a stabilizer for ICG, since ICG is highly unstable and tends to degrade in aqueous solution.54 The optical properties of 1% Liposyn solution are shown in Table I. The optical properties of the fluorescence agent (ICG) were determined from its concentration and extinction coefficient (also listed in Table I).33 The light transmission loss (~ 10%) due to acrylic sphere was taken into account for the quantum efficiency of ICG (during simulation studies).

TABLE I
Optical properties of fluorescent target and non-fluorescent background for all phantom studies

Two sets of experimental studies were performed. In experimental study I, the effect of implementing the unique simultaneous multiple point illumination geometry for 3-D tomographic imaging was assessed in comparison to the conventionally used sequential multiple point illumination geometry. In this study, the slab phantom had a single target located 1 cm deep from the imaging surface (i.e. at [5 2.5 1] cm). Fluorescence measurements were acquired using both the illumination geometries, under the same experimental conditions. For sequential multiple point illumination, each of the 6 point sources were illuminated sequentially (while the rest of source fibers were covered by light-blocking materials) and the fluorescence measurements from the 165 collection points were acquired simultaneously. During both the simultaneous and sequential illumination geometry based studies, the integration time of the camera was constant (0.2 sec) and only the gain settings at the image intensifier end was adjusted in order to amplify the weak fluorescence signals. The subtraction technique was implemented to the fluorescence measurements obtained from both the illumination geometries and the measurements were referenced (as described in Section II A.3) prior to performing the 3-D image reconstructions. The 3-D tomographic results obtained using both the illumination geometries were compared (as described in Section III A).

In experimental study II, the effect of implementing the unique simultaneous multiple point illumination geometry on target depth recovery was assessed. In this study, the tissue phantom had a single target located at 5 different x-y locations, and under two target depths of 1.5 and 2 cm deep from the imaging surface (see Table II and Fig.3). In these studies as well, the integration time of the camera was constant (0.2 sec) and the gain settings at the image intensifier were varied in order to amplify the weak fluorescence signals. For all phantom studies, the probe head was positioned ~2mm away from the phantom surface without actual contacting it. The fluorescence amplitude (AC) measurements obtained under different experimental conditions are presented as 2-D surface contour plots, post to implementing the subtraction technique. Three-dimensional image reconstructions were performed using the referenced fluorescence measurements (post to implementing the subtraction technique) and the tomographic results are described in Section III B.

FIG. 3
Target’s X-Y location with respect to the slab phantom for different experimental cases (A-E) shown as white hollow circles. The same X-Y locations were used for different target depths of 1.5 and 2.0 cm deep, during the heterogeneous phantom ...
TABLE II
The target location details for different experimental cases (A-E) at two different target depths

C. Three-dimensional image reconstructions

1. Forward problem

A Galerkin finite element based method 55 was used to solve for the forward problem of the coupled diffusion equations. The 3-D phantom of a 10×6.5×10cm3 (since 1% Liposyn was filled only up to 650 ml in the 10×10×10cm3 acrylic container) was generated and discretized into tetrahedral elements using Gambit 2.1.6 software (Fluent Inc., Lebanon, NH). The simulated phantom consisted of 10835 nodes and 53467 volume elements. The illumination/collection surface (at Z=0 cm plane) consisted of 2080 elements and 1160 nodes, among which 171 nodes were corresponding to the 6 illuminating and 165 detecting (prefixed) point locations. Partial-current boundary condition 56 was applied to all the 6 boundary surfaces, of which the top surface (Z=0 cm) is Liposyn/air interface and rest of the surfaces are acrylic/air interfaces. The input parameters for forward simulation studies were: i) the optical properties of the experimental phantom; and ii) the output power of each source fiber, which was measured using an optical power meter (PM 100, Thorlabs Inc., Newton, NJ). Although equal intensity of six sources had been attempted, the actual source strengths were unequal, and hence accounted for in the forward model.

2. Inverse problem

Three-dimensional image reconstructions were performed using a computationally efficient version of the Approximate Extended Kalman Filter (AEKF) algorithm, 17,33 in order to obtain the 3-D optical property map of the fluorescence absorption coefficient at excitation wavelength (μaxf ). The extensive details of the complete AEKF algorithm have been described elsewhere,17,32-35 and hence briefly described here. The AEKF algorithm recursively minimizes the variance of parameter error (ie. error in μaxf for this study), given the estimation of measurement error covariance R, model error covariance Q, and parameter error covariance P.57 The measurement error covariance R was estimated from the variances of the means of 5 repeated experimental measurements from each detector points. From preliminary experimental studies, it was observed that upon increasing the number of repeated measurements from 5 to 10, there was not a marked reduction in the measurement error for the cost of doubling the image acquisition time. Hence, the number of repetitions was limited to 5, in order to minimize the total imaging time. The model error covariance Q was empirically chosen to be equal to ¼ of the measurement error covariance R17,33,57 The covariance of measurement error and model error were implemented in the inversions to weight the updates at each iteration. The parameter error covariance P is defined as the error in the spatially distributed (unknown) parameter values (μaxf), which is used to damp into the inversion for better convergence. The reconstructions were assumed to have converged, when the root mean square output error (RMSE) was less than 1%, or the total number of iterations exceeded 50 times. The initial guess of reconstruction parameter μaxf was arbitrarily set to 0.003 cm-1 for all the reconstruction cases. The same reconstruction algorithm was implemented to measurements obtained from both sequential and simultaneous multiple point illumination geometry, with minor changes to the way the fluorescence measurements were implemented in the Matlab-based inversion code. Three-dimensional iso-surface contour maps of the reconstructed parameter (μaxf) of the entire phantom volume are presented in TecPlot, in order to qualitatively determine the reconstructed target size and location. For quantitative estimation of the reconstructed target with respect to the background, a cut-off value of μaxf was selected based on the first break point of histogram plot of μaxf. The comparison of reconstructed vs. true target details are tabulated in terms of: (i) target centroid location; (ii) target volume, which is the total volume of all the elements whose μaxf is higher than the chosen cut-off value; and (iii) target μaxf , which is the volume weighted μaxf value of all the elements whose μaxf is above the cut-off value as shown in following formula:

equation M2
(4)

where suffix element represents the tetrahedral elements of the finite element mesh.

III. RESULTS AND DISCUSSION

A. Study I: Comparison of sequential vs. simultaneous illumination

Three-dimensional image reconstructions were compared for the two illumination based geometries, i.e. sequential and simultaneous, for the experimental study I employing a single 1-cm deep target under perfect uptake conditions. The qualitative reconstruction results are shown in Fig.4 and the quantitative details are provided in Table III. From Fig. 4, it is evident that the novel hand-held based optical imager employing the simultaneous multiple point illumination geometry has demonstrated the feasibility of 3-D tomographic imaging (for the first time). In addition, the target was comparably reconstructed using both the illumination geometries. This indicates that the simultaneous multiple point illumination geometry tends to maintain the quality of image reconstructions as obtained using the conventionally used sequential multiple point illumination geometry, without introducing any artifacts. This is possibly because the current experiments were carried out under the perfect uptake conditions where there was no fluorescence in the background, thus minimizing the probability of artifacts using either illumination geometry. Extensive fluorescence tomography studies are currently carried out under different experimental conditions of optical contrast ratios (i.e. include fluorescence in the background as well), target size, number and volume in order to validate the effectiveness of simultaneous multiple point illumination geometry over the conventional sequential multiple point illumination geometry.

FIG. 4
Iso-surface contour plots of the 3-D image reconstructions presented in the X-Y and X-Z planes for heterogeneous phantoms under 1:0 fluorescence contrast ratio at a fixed target location (x=5cm, y=2.5cm, z=1cm). The hollow circle in each plot represented ...
TABLE III
Details of 3-D image reconstructions using simultaneous and sequential illumination for the perfect uptake (1:0) phantom study at target 3D location of (x=5cm, y=2.5cm and z=1cm)

Apart from the comparable image quality using both the measurement geometries, the simultaneous illumination based geometry is advantages over the sequential illumination geometry, in terms of the overall imaging time, which was reduced to 1/6th of the actual data acquisition time (without even accounting for the optical switching time and the time lost in altering the gain settings for each source during sequential illumination). During the initial design of this novel hand-held probe51, the source-detector layout on the probe head was appropriately chosen such that a constant setting of the integration time and gain at the detector end was sufficient to recover optical signals from most of the detectors above the noise floor, when implementing simultaneous multiple point illumination geometry.

B. Study II: Phantom studies using simultaneous illumination

1. Two-dimensional surface imaging

The 2-D surface contour plots (in X-Y plane) of fluorescence amplitude (AC) measurements on the illumination/collection surface for 5 different target locations, and two target depths (1.5cm and 2cm deep) are presented in Fig. 5. When the target was positioned at 1.5cm deep (Fig. 5(a)), for any given X-Y location, the target’s presence and 2-D location can be observed close to the true location. When the target was located 2cm deep (Fig. 5(c)), for any given X-Y location, the target’s location was observed, although it was not distinctly differentiated from the background. This is possibly from the weakened fluorescence signals, which exponentially decrease with target depth.

FIG. 5
Two-dimensional surface contour plots of detected fluorescence amplitude on the imaging plane of heterogeneous phantoms under 1:0 fluorescence contrast ratio, at different target depths of 1.5 cm (rows (a)) and 2.0 cm deep (rows (b)), and different experimental ...

With the implementation of a real-time co-registering technique, which is a focus of our ongoing efforts, a co-registered 2-D location of the target on 3-D breast tissues can serve as a real-time imaging tool (similar to ultrasound) in a clinical environment. Besides, the estimated 2-D target location can also be used as a priori information during 3-D tomography studies (part of our future studies) to reduce the computational time and possibly improve target localization, especially in large tissue volumes. However, in order to obtain a more detailed target information such as target size, depth and optical properties, 3-D tomography studies is indispensable.

2. Three-dimensional image reconstructions

Three-dimensional image reconstructions were performed for all the 10 experimental cases from study II, as described in Section II. B. The 3-D iso-surface contour plots are presented at X-Y and X-Z plane in order to illustrate the size and location of reconstructed target with respect to the true target size/location, for the 1.5cm (Fig. 6(a)) and 2cm deep target(s) (Fig. 6 (b)). In all the tomographic reconstructions, a single target was reconstructed without any artifacts, at both the target depths. Details of the reconstructed target location, size and optical property (μaxf) for the 1.5cm and 2cm deep target(s) are given in Table TableIVIV and andV,V, respectively. These quantitative results showed that most of the reconstructed parameters for all experimental cases were recovered, for both the target depth cases.

FIG. 6
Iso-surface contour plots of the 3-D image reconstructions presented in the X-Y and X-Z planes for heterogeneous phantoms under 1:0 contrast ratio, different target depths of (a) 1.5 cm, and (b) 2.0 cm. The five columns in each row represent the five ...
TABLE IV
Details of 3-D image reconstructions for the perfect uptake (1:0) phantom study under a target depth of 1.5 cm
TABLE V
Details of 3-D image reconstructions for the perfect uptake (1:0) phantom study under a target depth of 2.0 cm

However, in all the experimental cases, the reconstructed target depth (Z-axis) was closer to the phantom surface (observed even from qualitative results in Fig. 6). This could possibly be due to the fact that only reflectance measurements are acquired (from flat phantom surfaces) and used towards reconstructions. Simulations studies have been performed prior to experimental studies, using reflectance measurement obtains from similar target volumes (0.45 cc) and depths (1.5 cm and 2.0 cm deep) as in the experimental cases. The simulation tomography results (using simulated measurements without adding any noise) also reconstructed the target(s) closer to the surface (not shown due to brevity). This further confirms that it is probably the underlying physics of reflectance based measurement that affects the depth recovery.

The difficulty in depth recovery from 3-D tomography studies using only reflectance based measurements was also reported by other researchers.32,58-61 Although different researchers used different measurement geometries (sequential multiple point illumination, and area illumination), instrumentation set-ups, and reconstruction algorithms, all the studies commonly acquired only reflectance-based measurements. Consistently it was observed that the reconstructed target was closer to the surface for different depth studies. In our current study, since a different illumination scheme (multiple point simultaneous illumination) has been employed, tomography studies were performed at various target depths (reaching the reported limits of 1.5-2 cm deep) in order to determine if the overall improved signal strength from simultaneous multiple point illumination over sequential illumination34,51 improves the depth recovery. However, the current studies demonstrated that the limitation in depth recovery was independent of the nature of illumination geometry, and was still affected by the physics of reflectance-based measurements.

These tomography studies also imply that trans-illumination based measurements are possibly required apart from reflectance measurements, in order to recover the target depth accurately. Tomographic imaging studies using trans-illumination measurements acquired by the hand-held based optical imager is the focus of our future work on curved phantoms under different experimental conditions. During our instrumentation design and development phase, the probe head was built such that it is capable of contouring to the tissue curvatures (i.e. the curved phantoms) and obtain trans-illumination measurements apart from reflectance measurements.

IV. CONCLUSIONS

A novel hand-held probe based optical imaging system has been developed with the following unique features: (i) a hand-held optical probe with multiple sources and detectors that allows imaging over larger areas (5×10 cm2) in a single scan, (ii) simultaneous illumination and collection based measurement geometry that allows rapid data acquisitions (< 20 sec per location of the probe); (iii) a gain-modulated ICCD detector for rapid data acquisitions of weaker optical (fluorescence) signals; (iv) a flexible hand-held probe that can contour along any tissue curvature (part of future studies to implement the probe on curved surfaces); and (v) a portable set-up to allow bedside imaging in a clinical environment. Three-dimensional fluorescence tomography studies have been successfully demonstrated for the first time using a hand-held based optical imager. From the preliminary tomography studies on large slab phantoms (~ 650 ml), target(s) of 0.45 ml were reconstructed at various target depths (1-2 cm), upon using ICG under perfect uptake conditions. The target depth recovery was limited (as observed by other researchers as well) and appears to be predominantly affected by the physics of reflectance-based measurements.

Excitation leakage significantly contaminated the acquired weak fluorescence measurements, especially when the target was located deep (here, 2cm). The use of appropriate optical (interference) filters along with a subtraction-based post-processing technique remarkably aided in removing the artifacts arising from excitation leakage, both in 2-D surface contour plots and in 3-D image reconstructions. Extensive 3-D tomography studies are currently carried out using different target depths, volumes, number, and target:background optical contrast ratios in order to determine the performance limits and resolution of our optical imager. In addition, future studies will implement the rough estimate of the target location (as obtained from 2-D surface contour plots) as a priori information, in order to improve the quality of reconstructions, and also reduce the computational time. Parallel to these studies, experimental are carried out using the hand-held probe on curved phantoms (since the probe head is flexible to contour along curvatures up to 45°), wherein both reflectance and trans-illumination based measurements can possibly improve the quality of depth reconstructions. Efforts are made towards developing innovative methods to obtain real-time 3-D co-registered images of the tissue surface during multiple scans, using the hand-held optical imager.

In summary, the novel hand-held optical imaging system with its unique features, has potential to significantly impact the clinical translation of the optical imaging technology as a bed-side portable hand-held based imaging device for 3-D target localization without the use of a second imaging modality (e.g. an expensive MRI image) for co-registration. On a long-term, the developing hand-held based optical imaging technology also has applicability towards prognostic breast imaging, image-guided breast cancer therapy, and volumetric analysis of tumors in response to chemotherapeutic trials.

ACKNOWLEDGEMENTS

This work was funded in part by the Florida Department of Health (06BB-08) and National Institutes of Health grant (R15CA119253-01A1)

REFERENCES

1. Franceschini MA, Moesta KT, Fantini S, Gaida G, Gratton E, Jess H, Mantulin WW, Seeber M, Schlag PM, Kaschke M. Frequency-domain techniques enhance optical mammography: Initial clinical results. Proc. Natl. Acad. Sci. USA. 1997;94:6468–6473. [PubMed]
2. Fantini S, Walker SA, Franceschini MA, Kaschke M, Schlag PM, Moesta KT. Assessment of the size, position, and optical properties of breast tumors in vivo by noninvasive optical methods. Appl. Opt. 1998;37:1982–1989. [PubMed]
3. Moesta KT, Fantini S, Jess H, Totkas S, Franceschini MA, Kaschke M, Schlag PM. Contrast features of breast cancer in frequency-domain laser scanning mammography. J. Biomed. Opt. 1998;3:129–136. [PubMed]
4. Grosenick D, Wabnitz H, Rinneberg HH, Moesta KT, Schlag PM. Development of a time-domain optical mammograph and first in vivo applications. Appl. Opt. 1999;38:2927–2943. [PubMed]
5. Colak SB, van der Mark MB, ’t Hooft GW, Hoogenraad JH, van der Linden ES, Kuijpers FA. Clinical Optical Tomography and NIR Spectroscopy for Breast Cancer Detection. IEEE Journal of Selected Topics in Quantum Electronics. 1999;5:1143–1158.
6. Pogue BW, Poplack SP, McBride TO, Wells WA, Osterman KS, Osterberg UL, Paulsen KD. Quantitative hemoglobin tomography with diffuse near-infrared spectroscopy: Pilot results in the breast. Radiology. 2001;218:261–266. [PubMed]
7. Jiang H, Xu Y, Iftimia N, Eggert J, Klove K, Baron L, Fajardo L. Three-dimensional optical tomographic imaging of breast in a human subject. IEEE Trans. Med. Imaging. 2001;20:1334–1340. [PubMed]
8. Li A, Boverman G, Zhang Y, Brooks D, Miller EL, Kilmer ME, Zhang Q, Hillman EM, Boas DA. Optimal linear inverse solution with multiple priors in diffuse optical tomography. Appl. Opt. 2005;44(10):1948–56. [PubMed]
9. Zhang Q, Brukilacchio TJ, Li A, Stott JJ, Chaves T, Hillman E, Wu T, Chorlton M, Rafferty E, Moore RH, Kopans DB, Boas DA. Coregistered tomographic x-ray and optical breast imaging: initial results. J. Biomed. Opt. 2005;10(2):24033. [PubMed]
10. Gibson AP, Hebden JC, Arridge SR. Recent advances in diffuse optical imaging. Phys Med Biol. 2005;50(4):R1–43. [PubMed]
11. Hawrysz DJ, Sevick-Muraca EM. Developments toward diagnostic breast cancer imaging using near-infrared optical measurements and fluorescent contrast agents. Neoplasia. 2000;2:388–417. [PMC free article] [PubMed]
12. Licha K. Contrast agents for optical imaging. Top. Curr. Chem. 2002;222:1–29.
13. Sevick-Muraca EM, Houston JP, Gurfinkel M. Fluorescence-enhanced, near infrared diagnostic imaging with contrast agents. Curr. Opin. Chem. Biol. 2002;6:642–650. [PubMed]
14. Sevick-Muraca EM, Godavarty A, Houston JP, Thompson AB, Roy R. Near-infrared imaging with fluorescent contrast agents. In: Pogue Brian W., Mycek Mary-Ann., editors. Handbook of Biomedical Fluorescence. Marcel Dekker Inc.; New York: 2003.
15. Gurfinkel M, Ke S, Wen X, Li C, Sevick-Muraca EM. Near-infrared fluorescence optical imaging and tomography. Dis. Markers. 2004;19(23):107–121. [PubMed]
16. Pogue BW, McBride TO, Osterberg UL, Paulsen KD. Comparison of imaging geometries for diffuse optical tomography of tissue. Opt. Express. 1999;4(8):270–286. [PubMed]
17. Eppstein MJ, Hawrysz DJ, Godavarty A, Sevick-Muraca EM. Three-dimensional near-infrared fluorescence tomography with Bayesian methodologies for image reconstruction from sparse and noisy data sets. Proc. Natl. Acad. Sci. USA. 2002;99:9619–9624. [PubMed]
18. Milstein B, Stott JJ, Oh S, Boas DA, Millane RP, Bouman CA, Webb KJ. Fluorescence optical diffusion tomography using multiple-frequency data. J. Opt. Soc. Am. A. 2004;21(6):1035–1049. [PubMed]
19. Lee J, Sevick-Muraca EM. 3-D Fluorescence enhanced optical tomography using references frequency-domain photon migration measurements at emission and excitation measurements. J. Opt. Soc. Am. A. 2002;19:759–771. [PubMed]
20. Quaresima V, Matcher SJ, Ferrari M. Identification and quantification of intrinsic optical contrast for near-infrared mammography. Photochem. Photobiol. 1998;67:4–14. [PubMed]
21. Li X, Culver J, Durduran T, Chance B, Yodh AG, Pattanayak DN. Advances in Optical Imaging and Photon Migration (AOIPM) OSA publications; Orlando, FL: 1993. Diffraction tomography with diffuse-photon density waves: clinical studies and background subtraction.
22. 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. [PubMed]
23. Yang Y, Iftimia N, Xu Y, Jiang H. Chance B, Alfano RR, Tromberg BJ, editors. Frequency-domain fluorescent diffusion tomography of turbid media and in vivo tissues. Proc. Soc. Photo-Opt. Instrum. Eng. 2001;4250:537–545. Optical Tomography and Spectroscopy of Tissue IV.
24. Ntziachristos V, Bremer C, Weissleder R. Fluorescence imaging with near-infrared light: new technological advances that enable in vivo molecular imaging. Eur Radiol. 2003;13:195–208. [PubMed]
25. Ntziachristos V, Weissleder R. Charge-coupled-device based scanner for tomography of fluorescent near-infrared probes in turbid media. Med. Phys. 2002;29:803–809. [PubMed]
26. Ntziachristos V, Bremer C, Tung C, Weissleder R. Imaging cathepsin B up-regulation in HT-1080 tumor models using fluorescence-mediated molecular tomography (FMT) Acad. Radiol. 2002;9(suppl 2):S323–S325. [PubMed]
27. Schulz RB, Ripoll J, Ntziachristos V. Noncontact optical tomography of turbid media. Opt. Lett. 2003;28(18):1701–1703. [PubMed]
28. Ntziachristos V, Schellenberger EA, Ripoll J, Yessayan D, Graves E, Bogdanov A, Josephson JL, Weissleder R. Visualization of antitumor treatment by means of fluorescence molecular tomography with an annexin V-Cy5.5 conjugate. Proc. Natl. Acad. Sci. USA. 2004;101(33):12294–12299. [PubMed]
29. Schulz RB, Ripoll J, Ntziachristos V. Experimental fluorescence tomography of tissues with noncontact measurements. IEEE Trans. Med. Imaging. 2004;23(4):492–500. [PubMed]
30. Bremer C, Ntziachristos V, Weitkamp B, Theilmeier G, Heindel W, Weissleder R. Optical imaging of spontaneous breast tumors using protease sensing ‘smart’ optical probes. Invest. Radiol. 2005;40(6):321–7. [PubMed]
31. Colak SB, Papaioannou DG, t’Hooft GW, van der Mark MB, Schomberg H, Paasschens JCJ, Melissen JB, van Asten NJ. Tomographic image reconstruction from optical projections in light-diffusing media. Appl. Opt. 1997;36:180–213. [PubMed]
32. Godavarty A, Thompson AB, Roy R, Eppstein MJ, Zhang C, Gurfinkel M, Sevick-Muraca EM. Diagnostic imaging of breast cancer using fluorescence-enhanced optical tomography: phantom studies. J. Biomed. Opt. Special edition on Biomedical Optics and Women’s Health. 2004;9(3):488–496. [PubMed]
33. Godavarty A, Eppstein MJ, Zhang C, Theru S, Thompson AB, Gurfinkel M, Sevick-Muraca EM. Fluorescence-enhanced optical imaging in large tissue volumes using a gain modulated ICCD camera. Phys. Med. Biol. 2003;48:1701–1720. [PubMed]
34. Godavarty A, Zhang C, Eppstein MJ, Sevick-Muraca EM. Fluorescence-enhanced optical imaging of large phantoms using single and simultaneous dual point illumination geometries. Med. Phys. 2004;31:183–90. [PubMed]
35. Godavarty A, Sevick-Muraca EM, Eppstein MJ. Three-dimensional fluorescence lifetime tomography. Med. Phys. 2005;32:992–1000. [PubMed]
36. Zhu B, Eppstein MJ, Sevick-Muraca EM, Godavarty A. Noise pre-filtering techniques in fluorescence-enhanced optical tomography. Opt. Exp. 2007;15(18):11285–11300. [PubMed]
37. Tromberg BJ, Coquoz O, Fishkin JB, Pham T, Anderson ER, Butler J, Cahn M, Cross JD, Venugopalan V, Pham D. Non-invasive measurements of breast tissue optical properties using frequency-domain photon migration. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 1997;352(1354):661–668. [PMC free article] [PubMed]
38. Pham TH, Coquoz O, Fishkin JB, Anderson E, Tromberg BJ. Broad bandwidth frequency domain instrument for quantitative tissue optical spectroscopy. Rev. Sci. Instrum. 2000;71(6):2500–2513.
39. Shah N, Cerussi A, Eker C, Espinoza J, Butler J, Fishkin J, Hornung R, Tromberg B. Noninvasive functional optical spectroscopy of human breast tissue. Proc. Natl. Acad. Sci. U.S.A. 2001;98(8):4420–4425. [PubMed]
40. Jakubowski DB, Cerussi AE, Bevilacqua F, Shah N, Hsiang D, Butler J, Tromberg BJ. Monitoring neoadjuvant chemotherapy in breast cancer using quantitative diffuse optical spectroscopy: a case study. J. Biomed. Opt. 2004;9(1):230–238. [PubMed]
41. Chen N, Huang M, Xia H, Piao D. Portable near-infrared diffusive light imager for breast cancer detection. J. Biomed. Opt. 2004;9(3):504–510. [PMC free article] [PubMed]
42. Shah N, Cerussi AE, Jakubowski D, Hsiang D, Butler J, Tromberg BJ. Spatial variations in optical and physiological properties of healthy breast tissue. J. Biomed. Opt. 2004;9(3):534–540. [PubMed]
43. Chance B, Nioka S, Zhang J, Conant EF, Hwang E, Briest S, Orel SG, Schnall MD, Czerniecki BJ. Breast cancer detection based on incremental biochemical and physiological properties of breast cancers: A six-year, two-site study. Acad. Radiol. 2005;12(8):925–933. [PubMed]
44. Tromberg BJ. Optical scanning and breast cancer. Acad. Radiol. 2005;12(8):923–924. [PubMed]
45. No KS, Chou PH. Mini-FDPM and Heterodyne Mini-FDPM: Handheld Non-Invasive Breast Cancer Detectors Based on Frequency Domain Photon Migration. IEEE Trans. Circuits Syst. [Circuits and Systems I: Fundamental Theory and Applications] 2005;52(12):2672–2685.
46. Durduran T, Choe R, Yu G, Zhou C, Tchou JC, Czerniecki BJ, Yodh AG. Diffuse optical measurement of blood flow in breast tumors. Opt. Lett. 2005;30(21):2915–2917. [PubMed]
47. Chance B, Zhao Z, Wen S, Chen Y. Simple ac circuit for breast cancer detection and object detection. Rev. Sci. Instr. 2006;77:064301.
48. No KS, Xie Q, Kwong R, Cerussi A, Tromberg BJ, Chou P. HBS: a handheld breast cancer detector based on frequency domain photon migration with full heterodyne. Proc. IEEE BioCAS. 2006
49. Bremer C, Ntziachristos V, Weissleder R. Optical-based molecular imaging: Contrast agents and potential medical applications. Eur. Radiol. 2003;13:231–243. [PubMed]
50. Roy R, Thompson AB, Godavarty A, Sevick-Muraca EM. Tomographic fluorescence imaging in tissue phantoms: A novel reconstruction algorithm and imaging geometry. IEEE Trans. Med. Imaging. 2005;24(2):137–154. [PubMed]
51. Jayachandran B, Ge J, Regalado S, Godavarty A. Design and development of a hand-held optical probe towards fluorescence diagnostic imaging. J. Biomed. Opt. 2007;12(5):054014. [PubMed]
52. Lakowicz JR, Berndt KW. Lifetime-selective fluorescence imaging using an rf phase-sensitive camera. Rev. Sci. Instrum. 1991;62(7):1727–1734.
53. Spinelli L, Torricelli A, Pifferi A, Taroni P, Danesini GM, Cubeddu R. Bulk optical properties and tissue components in the female breast from multiwavelength time-resolved optical mammography. J. Biomed. Opt. 2004;9(6):1137–1142. [PubMed]
54. Rajagopalan R, Uetrecht P, Bugaj JE, Achilefu SA, Dorshow RB. Stabilization of the optical tracer agent indocyanine green using noncovalent interactions. Photochem. Photobiol. 2000;71:347–350. [PubMed]
55. Fedele F, Laible JP, Eppstein MJ. Coupled complex adjoint sensitivities for frequency-domain fluorescence tomography: Theory and vectorized implementation. J. Comput. Phys. 2003;187:597–619.
56. Haskell RC, Scasssand LO, Tsay T-T, Feng T-C, Adams MSM, Tromberg BJ. Boundary conditions for the diffusion equation in radiative transfer. J. Opt. Soc. Am. A. 1994;11(27272741) [PubMed]
57. Eppstein MJ, Dougherty DE, Hawrysz DJ, Sevick-Muraca EM. Three dimensional bayesian optical image reconstruction with domain decomposition. IEEE Trans. Med. Imaging. 2001;20(3):147–163. [PubMed]
58. Kepshire DS, Davis SC, Dehghani H, Paulsen KD, Pogue BW. Subsurface diffuse optical tomography can localize absorber and fluorescent objects but recovered image sensitivity is nonlinear with depth. Appl. Opt. 2007;46(10):1669–1678. [PubMed]
59. Kepshire DS, Davis SC, Dehghani H, Paulsen KD, Pogue BW. Challenges in sub-surface fluorescence diffuse optical imaging. Proc. SPIE. 2007;6434(64340V):1–9.
60. Joshi A, Bangerth W, Hwang K, Rasmussen JC, Sevick-Muraca EM. Fully adaptive FEM based fluorescence optical tomography from time-dependent measurements with area illumination and detection. Med. Phys. 2006;33(5):1299–1310. [PubMed]
61. Joshi A, Bangerth W, Hwang K, Rasmussen JC, Sevick-Muraca EM. Plane-wave fluorescence tomography with adaptive finite elements. Opt. Lett. 2006;31(2):193–195. [PubMed]