PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of transaThe Royal Society PublishingPhilosophical Transactions AAboutBrowse by SubjectAlertsFree Trial
 
Philos Trans A Math Phys Eng Sci. 2011 November 28; 369(1955): 4531–4557.
PMCID: PMC3263784

Implicit and explicit prior information in near-infrared spectral imaging: accuracy, quantification and diagnostic value

Abstract

Near-infrared spectroscopy (NIRS) of tissue provides quantification of absorbers, scattering and luminescent agents in bulk tissue through the use of measurement data and assumptions. Prior knowledge can be critical about things such as (i) the tissue shape and/or structure, (ii) spectral constituents, (iii) limits on parameters, (iv) demographic or biomarker data, and (v) biophysical models of the temporal signal shapes. A general framework of NIRS imaging with prior information is presented, showing that prior information datasets could be incorporated at any step in the NIRS process, with the general workflow being: (i) data acquisition, (ii) pre-processing, (iii) forward model, (iv) inversion/reconstruction, (v) post-processing, and (vi) interpretation/diagnosis. Most of the development in NIRS has used ad hoc or empirical implementations of prior information such as pre-measured absorber or fluorophore spectra, or tissue shapes as estimated by additional imaging tools. A comprehensive analysis would examine what prior information maximizes the accuracy in recovery and value for medical diagnosis, when implemented at separate stages of the NIRS sequence. Individual applications of prior information can show increases in accuracy or improved ability to estimate biochemical features of tissue, while other approaches may not. Most beneficial inclusion of prior information has been in the inversion/reconstruction process, because it solves the mathematical intractability. However, it is not clear that this is always the most beneficial stage.

Keywords: spectroscopy, molecular, tomography, reconstruction, near-infrared, optical

1. Introduction

Imaging using near-infrared spectroscopy (NIRS) has diversified and expanded into several fields of medical research and science, and pre-packaged biomedical instruments are now commonplace. While the rate at which new instruments have been introduced outpaces their utility, research in this field continues to expand. This expansion is driven by a rapidly increasing pace of discoveries in molecular diagnostic imaging at the pre-clinical stage, with potential to profoundly impact basic biology and clinical medicine. Most implementations of NIRS use a model-based estimation strategy for quantifying tissue parameters based upon spectral features and thus incorporate some form of prior information. For example, NIRS applications include diffuse optical tomography as well as a wide variety of other tissue interrogation approaches where biochemical information is derived using approaches in which spatial tissue sampling is more limited. One aspect that is shared by all NIRS methodologies is that they rely on the acquisition of diffuse near-infrared (NIR) light signals followed by the solution of an inverse problem in the form of either a tomography problem or a simpler type of least-squares problem, such as curve fitting. In recent years, many sources of prior information have been examined, taking several different forms and being analysed quantitatively in terms of their potential to make NIRS amenable to a level of accuracy that is closer to, if not better than, that achievable with concurrent imaging modalities that are providing functional and/or molecular tissue information. The conceptual framework examined in this study is to identify explicit prior information as a multi-parametric factor that can be inserted into the NIRS process at several stages. This is illustrated in the schematic of figure 1. The use of prior information in measuring, processing and interpreting NIRS data is examined with the goal of identifying key examples in which explicit prior information benefits the ability to quantify and/or diagnose.

Figure 1.
The inclusion of prior information of several different types into the NIRS imaging process has been implicit in most applications. The reason for using prior information is to improve the final interpretation of the results, through more accurate quantification, ...

The general framework illustrated in figure 1 is conceptual and can be used to systematically examine at what stages prior information improves NIRS. The NIRS process consists of the following steps: (i) data acquisition, (ii) calibration and pre-processing of the data, (iii) forward light transport modelling, (iv) inverse problem, (v) post-processing, and (vi) interpretation of the parameters. Not all of these steps are necessarily used, but this generically describes the possible steps, and into each of these steps it is feasible to involve prior information from any number of data streams. These possible streams may consist of (i) spatial structure information obtained from other imaging modalities, (ii) spectral information of the molecules being assayed, (iii) physical constraints on parameters in the model and data, (iv) subject or tissue biomarkers and/or demographic information, and (v) biophysical models to modify or interpret the NIRS process data. These are not fundamentally separate prior information features but rather general categories.

Historically, a majority of NIRS systems rely on spectral prior information as a basis for interpreting the signals in the process of using this for diagnosis. NIRS spectral features measured in most cases are broad, provide a direct transformation to physiological parameters, and the prior measurement of them is possible without any change in system hardware. Because of these characteristics, spectral priors arguably have the most effective impact on NIRS quantification for medical applications. Implementation of spatial prior information as either region definitions, soft priors or shape-based priors has been a steady area of investigation, albeit with niche conclusions for specific examples that provide synergy between NIRS and another imaging system. In the sections below, prior information streams for NIRS are identified, and the details of their implementation and potential to impact NIRS quantification and possible diagnosis are discussed. An attempt is also made to quantify the proven and potential benefits presented by different prior–NIRS convolutions.

2. Prior information: structure

(a) External shape and internal fluence

Without question, the most dominant factors in NIRS measurement accuracy are the tissue shape and source–detector coupling to the tissue [1,2]. The external shape affects the internal light transport distribution and thereby affects the amplitude of the detected light signal at points on the surface. See figure 2a for an illustration of how the shape of the tissue may influence the light distribution. A large number of papers have been published describing logistical approaches to modelling light in tissues of different shapes and source–detector arrangements [3], and it is generally thought that using as accurate a light propagation model as possible is required in order to interpret the data or accurately estimate optical properties through an inversion of the measurements. In this case, the prior information is implicit in the design of the system and the forward modelling. Prior information can be used in the pre-processing of data in order to make the data more accurately match the model. For example, the first committed step for light transport modelling is determining, for each optical measurement, the spatial locations of where light enters tissue and where it is detected. This is achieved by first reconstructing the outer surface of the specimen using, for example, an optical profilometer or by co-registering the coordinate system of the NIRS system with that of an imaging modality providing structural images from which the outer surface of the tissue can be extracted. Then, a light transport modelling approach must be devised and strategies designed for finding numerical or analytical solutions that are subsequently used in the scope of the inverse problem allowing interpretation of the data. Improving the fidelity of the match between model and data through the use of prior information for source and detector positioning generally relies on algorithm designs that combine geometrical hardware considerations for the NIRS system design with co-registration tools allowing registration with the coordinate system of the instrument from which the prior information is obtained. As discussed in more detail in §2b, another important factor in the modelling chain that can significantly affect the match between model and data is the light propagation model itself.

Figure 2.
Major limitations of NIRS imaging. There are biophysical factors and internal tissue factors that limit the capability of NIRS imaging. (a,b) With respect to spatial issues, it is well known that the shape of the tissue changes the internal fluence ( ...

The modelling stage of the NIRS process is critical because it is reasonably well established that model–data mismatch errors can dominate the inversion process and lead to large bias errors [2,46]. However, measurement approaches can be used that minimize the need for accurate modelling. Such approaches can be categorized into several different data types, as outlined in table 1. These include: (i) ratiometric or derivative data at two or more different wavelengths, (ii) multiple-distance ratio or derivative data, (iii) small spatial volumes that either limit the effect of physical boundaries through scatter and/or absorption, or allow simpler empirical modelling, and (iv) temporal signals that are less sensitive to boundaries and/or more robustly insensitive to shape changes. The use of prior information about the tissue to be sampled is still essential in the design process with these systems, but can be implemented in the very first step of the NIRS process, namely data acquisition and calibration. This type of prior information is often used to eliminate the need for modelling completely, allowing simpler and perhaps more robust analysis of just the raw data or some processed version of the data. Examples of and references to these are shown in table 1.

Table 1.
The types of data acquisition implemented that implicitly reduce the need for model accuracy or can utilize simpler empirical models.

Approaches to analytically model NIRS light transport and allow direct inversion or fast inversion have dominated the successful commercial applications and led to several practical biological discoveries [3942]. Many advanced solutions to the diffusion equation have been developed for regular geometries [43] and are widely used in practical systems, although contact variation or tissue shape, which cannot be accurately accounted for through modelling, can be a major problem for quantitative parameter estimation. The ideal systems use data types that are robust and less sensitive to boundary shape changes and fibre contact issues [35,44], as described in table 1.

Systems that encode shape, location and external boundaries have been developed in many research and commercial systems. Figure 3 illustrates a handful of possible combinations, showing how inclusion of spatial information is done with additional hardware, and spectral prior inclusion is done with careful choice of the wavelength and source–detector bands.

Figure 3.
The use of spatial and spectral prior information has been widely studied. Spatial information can come from imaging systems such as X-ray, magnetic resonance imaging (MRI), ultrasound or surgical imaging. The integration of NIRS into these systems has ...

One of the earliest hybrid systems was demonstrated by Ntziachristos et al. [45], where concurrent magnetic resonance imaging (MRI) and NIRS of breast tumours were used to improve the estimation of indocyanine green uptake concentration in tumours. The examples shown in this study stimulated substantial interest in this approach.

(b) Internal shapes and physical limitations of light transport

One of the largest limitations of NIRS in tissue is the physical spread of light during transmission as a result of multiple scattering effects. This multiple scattering over large distances leads to diffuse propagation and the inability to resolve structures near the size of the transport scattering length (near 1 mm). While recovery of heterogeneous regions on the surface can approach the resolution of microscopy systems, the resolution degrades by orders of magnitude with depth into the tissue, and so measurements taken over several centimetres of tissue effectively become averages over larger volumes of tissue. See figure 2b for an illustration of this effect. This phenomenon places size limits on objects that can be resolved with NIRS [46,47]. Moon et al. [46] developed a useful metric that estimates that the photon bundle path width, as illustrated in figure 2b, is approximately 20 per cent of the distance between the source and detector, when sources and detectors are arranged in transmission geometry. This suggests that the sensitivity profile (sometimes called the photon measurement density function, or photon path width) is 2 mm for a 1 cm separation, and 2 cm for a 10 cm separation. These numbers provide a rough guide to the type of resolution that could be achieved, although it is generally accepted that regularized inversion solvers can improve this by nearly four times, such that a fundamental resolution for optical tomography in a 10 cm domain might be closer to 5 mm. This improvement in resolution comes in part from multiple overlapping source–detector projections through the domain.

The forward light transport model must have the flexibility to accurately model any spatial heterogeneity. There are several approaches used to solve this problem, largely separated into deterministic models, such as the Boltzmann equation for neutral particle transport, or stochastic approaches, such as a Monte Carlo model. Table 2 summarizes the major modelling approaches used, with a relative ranking of their ability to incorporate heterogeneity. Several solutions for Monte Carlo models exist, and the ability to incorporate arbitrary geometries and arbitrary internal heterogeneities with ease has improved dramatically in recent years, as developed by Fang [48]. The inversion of Monte Carlo data is still somewhat crippled by the need to run repeated forward solutions, which takes significant amounts of time unless very regular geometries are used [49,50]. However, the recent developments in graphics card processing units have made this a more realistic potential even for complex geometries [48,51]. A more common approach is the deterministic modelling of the diffusion equation, a second-order differential equation that is an approximation to the Boltzmann equation (a more complicated integro-differential equation). The transport equation itself can be directly solved with discrete ordinates solvers but requires substantial specialized computer resources and programming. As a consequence, the diffusion equation is simpler to solve than the transport equation in terms of both the mathematical complexity and the computational burden involved. Moreover, the diffusion equation depends on two rather than three phenomenological tissue parameters, namely, the absorption coefficient and the reduced scattering parameter. In fact, in the mathematical derivation leading to the diffusion equation, the scattering coefficient and the so-called phase function are absorbed into one parameter called the reduced scattering parameter. A critical point here, though, is that it is even more difficult, if not impossible, to devise generic approaches that would allow disentangling phase and scattering. However, accomplishing this would be the minimal step required for using the full benefits of the transport equation for forward modelling in NIRS. Moreover, the ability to incorporate accurate boundary conditions and internal heterogeneities also comes at a much greater difficulty level for approaches based on the Boltzmann equation.

Table 2.
The different models that have been used for optical radiation transport modelling are listed (left column) with relative rankings (using + marks) of their applicability and utility for accurate use with different levels of prior information: ...

An aspect that is sometimes regarded as an advantage for using the diffusion equation is the fact that many useful analytical solutions can be derived from it, while the number of solutions associated with the transport equation is much more limited. The transport solutions also typically involve several mathematical simplifying assumptions that limit the scope of their application. However, even in the case of the diffusion-based analytical equations, solving the internal heterogeneity must be accomplished using perturbation theory approaches, where the perturbed optical backgrounds are usually defined from analytical solutions, making such approaches reliable only for simple imaging geometries (e.g. tissue slabs in transmission, semi-infinite turbid media for epi-illumination imaging). These linear perturbation approaches are associated with the other drawback that the actual optical signal contrast rarely corresponds to optical property variations that satisfy the linear approximation from which the model is derived. Nevertheless, this has formed the basis for many reconstruction algorithms, allowing fast inversion with large datasets. The inclusion of arbitrary boundaries and arbitrary internal heterogeneity is arguably best achieved with a numerical model such as the finite element approach, where the domain is solved as a discrete defined set of points, linked by a mesh of finite basis functions. This approach was introduced for NIRS about 15 years ago, and has formed the basis of several available solvers. Some of these solvers have been developed to facilitate the use of a variety of priors, such as structural priors from MRI or X-ray tomosynthesis [52,53], spectral priors and biochemical priors.

(c) Explicit prior shape inclusion applications and analysis

The incorporation of internal structures in the inversion problem has been an area of increased academic interest in recent years. The prior information approach that is the less involved from both computational and mathematical complexity standpoints consists in pre-defining the imaging domain into a small number of specific regions. In this case, reduction of the dimensionality of the inversion is achieved by assuming the regions to be optically homogeneous. This approach is often called ‘hard prior’ information and was demonstrated in breast imaging by Ntziachristos et al. [45]. The composite approach to this consists in relaxing the homogeneity assumption in order to retain the original large dimensionality of the inversion domain. However, increased flexibility is conferred to the inversion problem by applying appropriately penalized inversion methods that are encoding the prior information but still allow optical property gradients to exist within each composite region of the domain. A common implementation of this latter approach used by several groups is to use a Laplacian-type regularization matrix that smoothes the parameter values within the pre-defined regions, but allows separate regularization between different tissue types [54,55]. This has been used in breast imaging with MRI-guided NIRS by Brooksby et al. [56] to more accurately estimate the concentrations of haemoglobin, oxygen saturation, water and scattering values in fibroglandular and adipose tissues in the normal breast. This ‘soft priors’ approach was validated in multiple tissue-simulating phantom studies as shown in figure 4 and demonstrated for tumour imaging by Carpenter et al. [57]. The fluorescence tomography version of this has been adopted by several groups [58], and has been used in several pre-clinical murine tumour imaging studies of contrast agent uptake [59,60].

Figure 4.
(a) Spatial prior inclusion can improve quantitative recovery, as shown in this MRI-guided phantom study, where images show the increasing level of spatial prior information used in the reconstruction algorithm. (b) Spectral prior inclusion can improve ...

The posterior incorporation of imaging shapes can be done in the analysis of NIRS data, and this is perhaps the most conventional way to approach combined modality data. A commercial prototype NIRS breast imager was shown to have significant discriminatory power between benign and malignant tumours [61]. The system was designed to be used in concert with mammogram analysis, thereby allowing posterior added information, for the improvement of the diagnostic reading, similar to the way ultrasound and mammography are interpreted together [62]. This added specificity has been confirmed in the summaries of two academic studies that show promise for interpretation of NIRS imaging in posterior combination with mammography [63,64]. Poplack et al. [64] showed that the combination of NIRS and mammography could increase the positive predictive power significantly. The explicit use of X-ray structure in NIRS inversion was demonstrated by Fang et al. [65] with a hybrid X-ray tomosynthesis system. This approach could be used to increase the sensitivity and/or specificity of the tomosynthesis exam depending upon the approach to use the prior information. Similarly, NIRS breast tumour imaging to quantify response from neoadjuvant chemotherapy with a posterior comparison to MRI by Choe et al. [66] demonstrated that contrast reduction from NIRS added to the diagnostic value of MRI. Additionally, Jiang et al. [67] examined the role of defining the region of interest based upon the size of the lesion in the initial pre-treatment image, and concluded that the quantitative estimate of tumour haemoglobin is significantly affected by the posterior estimate of lesion size. Thus, while posterior analysis is perhaps the least complex combination of NIRS and imaging data, it may become the most utilized as NIRS systems enter the clinic.

(d) Summary of benefits and limits

Several approaches were presented for which the inclusion of prior spatial information can increase the biological information content that can be derived from data acquired with NIRS instruments. Significant gains can be realized by using methodologies where the impact of tissue shape, source and detector coupling and light transport are minimized by an appropriate choice of data type, while still ensuring that sufficient tissue sampling is achieved to provide useful spatial resolution and quantification. Moreover, three approaches for inclusion of spatial information have been extensively examined to include internal tissue structures into the inverse problem, referred to here as (i) hard prior, (ii) soft prior and (iii) posterior use of the structural information. Of these, the one that is less technically involved is the third one, since it does not require that the prior information be fed directly to the NIRS reconstruction algorithm, thereby eliminating any potential reconstruction bias or the propagation of registration errors in the NIRS images. However, a significant limitation of this approach is that the information it can convey is limited by the intrinsic spatial resolution of diffuse optical imaging, namely from a few millimetres to a few centimetres depending on the interrogated volume size. In other words, although a tomography dataset might encode the information associated with a small area of optical contrast (i.e. smaller than the spatial resolution), NIRS images might not reveal it, although contrast can be resolved on the structural image. Alternatively, ‘hard prior’ and ‘soft prior’ methods can potentially reveal the biological information contained in the NIRS dataset precisely because the prior spatial information is directly encoded in the inverse problem. Of those two approaches, the method based on soft priors is preferable because of the reduced likelihood that spatial biases will be introduced during the inversion process. The implementation of soft prior constraints requires that it be tested carefully and the approach be well calibrated, otherwise the use of hard priors is perhaps a safer implementation. Yet, even with hard priors, it is critical to know that the regions identified are accurate, otherwise it is easy to recover false values in the regions.

3. Prior information: molecular and biochemical limits

(a) Molecular and biochemical limits

The inclusion of spectral fitting into data analysis has been ubiquitous throughout most of NIRS's development owing to its ability to recover concentrations of biochemical species rather than relative values or absorption coefficients. Systems with multiple wavelengths have the ability to directly measure the spectral constituent contributions to the signal independently, and it often does not require additional hardware. The value of spectral fitting has been established in many studies, and in bulk tissue spectroscopy without model-based inversion, this is often done either at the pre-processing stage or at the post-processing stage, without a significant difference in quantification between the two approaches. This is because of the linear relationship between the concentration of chromophores and optical absorption coefficients measured by NIRS in a bulk tissue measurement paradigm. However, in the NIRS imaging problem based on matrix inversion, the explicit inclusion of the pre-determined spectra in the inversion algorithm, rather than in a post-processing step, has been shown to dramatically improve the quantification and final interpretation. This is generally considered to be important because the spectral constraint inversion limits the solution space of the ill-posed problem [6870] and thereby improves accuracy spatially. This can be seen in figure 4b where a phantom containing higher blood concentration was imaged with high noise, resulting in blotchy image features when spectral information was applied after the absorption coefficient images had been recovered. The inclusion of the spectral prior constraint directly in the imaging algorithm improved the recovery significantly for this noisy dataset [71]. When compared with a spatial prior implementation, it was determined that spectral prior information was superior in terms of quantitative accuracy, but that the inclusion of both provided the most accurate quantification, in the setting of breast imaging [56,69]. Some optimally fitted values for breast tissues are shown in figure 4c.

The type and quality of prior information should be a key consideration in the design phase of NIRS systems. In addition to the obvious design considerations necessary for building hybrid systems, the use of an optimal wavelength set for particular chromophores should also factor strongly in the system design [72,73]. While this is a subject of ongoing analysis, it is probable that systems that can adaptively change wavelengths or map to key wavelengths will allow more accurate quantification of chromophores in tissue.

Spectral constraints in fluorescence imaging have been widely adopted in both research and commercial luminescence NIRS systems. Prior knowledge of the luminescent spectra can be used to separate the signals from the analyte of interest from contaminating background molecules or to image multiple analytes simultaneously. Major commercial systems incorporate algorithms to fit images at multiple wavelength bands to recover fluorophore concentrations, or use the approach to suppress or remove non-specific background signals [74]. These approaches are essential tools for low-concentration imaging, or multiple analyte imaging in vivo. Additionally, a large body of work has been produced to either correct for or exploit the distortion of luminescent spectra by the absorbing spectra of the tissue. This can correct quantification of signals that would otherwise be incorrect [7577] or be used to make an intractable problem solvable, such as in the case of bioluminescence tomography where spectral constraints are used to reduce the ill-posedness of the problem [7782]. However, an inherent limitation of such approaches is that quantification improvements heavily rely on prior knowledge not only of the absorption spectra of the tissue chromophores but also of the magnitude of their relative contributions. In cases where this information is not available, modelling errors can propagate into the inversion process, potentially negating any gains attributable to the spectral constraints. Research attempting to quantify this effect and minimize its impact in applications such as diffuse fluorescence tomography is ongoing.

The inclusion of physiological limits in the fitting or estimation process is often essential in the presence of noise, either from stochastic processes (e.g. photon-detection shot noise) or from modelling errors (e.g. neglecting to account for the presence of a chromophore in spectral fitting applications). Such constraints are often included in fitting algorithms without mention, yet it can be the key factor in obtaining accurate estimates, especially when multiple overlapping species are being estimated, such as in endogenous NIRS imaging. Most systems explicitly incorporate this in the inversion process, requiring a constrained inversion algorithm. The importance of this problem can be traced to the fact that, for any inversion problem, there generally exist several equivalent solutions that match the data for a given set of convergence criteria determined by a residual expression (e.g. the two-norm of the difference between experimental data and model) or a generalized norm that includes explicit prior information. In tomography, this problem is often referred to as ill-posedness and, because of noise, leads to the paradoxical situation that solutions that fit the measurements perfectly are always non-physical. Often, these non-physical solutions will, for example, include negative values for parameters that can only be positive (e.g. concentration of molecules) and lead to blood concentrations beyond or below what is physiologically possible. In such cases, inversion algorithms must be implemented where prior information in the form of bound constraints (e.g. non-negativity, lower and upper limits for recovered parameters) are introduced.

Endogenous NIRS in vivo breast imaging requires sparse sampling to cover the large volume of tissue, creating the ill-posedness problem discussed previously. From a practical standpoint, prior information has become commonplace in diffuse optical breast imaging, and is quintessential for the recovery of prognostically valuable images. Figure 5 gives an illustrative example of a reconstruction of data from a hybrid MRI–optical imaging system. In this example, both the hardware design and the reconstruction algorithm use spatial and spectral prior information. This system incorporates MRI images by collecting optical and MRI data simultaneously, providing a high-resolution spatial image from which to generate a boundary mesh. Additionally, adipose and fibroglandular tissue can be separated and assigned to specific nodes in the optical domain. By generating the optical solution space from the MRI images, the datasets are easily co-registered and a ‘hard priors’ reconstruction can be applied. Figure 5a shows an axial MRI image, which would be used to separate tissue types. The imaging system also uses spectral prior information in two ways. The six imaging wavelengths are optimized in the NIR window to maximize sensitivity to spectral features. As a result, each of the five chromophores seen in figure 5bg can be more robustly decoupled during image reconstruction. The reconstruction algorithm in this case also constrains the solution to be physiologically valid, as discussed earlier. In this example, the use of spectral and spatial priors reduces the number of unknowns significantly and allows the optical solution to more accurately quantify the tissue properties.

Figure 5.
Images of a single tumour imaged with NIR spectroscopy within the MRI: (a) axial MR image used for segmentation of fat, fibroglandular and tumour regions; (b) total haemoglobin overlay; (c) three-dimensional image of (b); (d) water image overlay on the ...

Posterior or post-processing analysis of spectra has been suggested in the interpretation of pathlength based upon assumed concentrations of absorbers, such as water, in tissue [83]. While this approach is not widely used, it is possible that such post-processing based upon spectra could provide easier or more accurate quantification.

(b) Summary of benefits and limits

For the situation of single source–detector measurements it is unclear whether prior knowledge in pre- versus post-processing is more beneficial, and it seems probable that the two are substantially similar, owing to the linearity of the data during processing. For imaging applications where matrix inversion is involved, which is ill-posed, there seems to be a significant benefit in the inversion step, owing to the magnitude of the regularization parameter relative to the matrix diagonal values. The evidence to date would support the fact that spectral prior information is perhaps even more important than spatial structure, but this is probably case-specific. More development and analysis in this area is warranted.

4. Prior information: patient demographics and biomarkers

(a) Patient demographics and biomarkers

Few studies have investigated incorporating demographic information about patients in NIRS; however, there are some niche applications where it may prove to be beneficial. NIRS data in breast imaging has shown strong correlations with demographic factors such as age, body mass index (BMI), radiographic breast density and menopausal/hormonal status. This correlation implies that there could be some benefit in using this prior information, perhaps to make systems simpler or more reliable. Examples of this approach include the use of radiographic density or BMI to estimate lipid concentration in the breast. Radiographic density is also highly correlated to optical scattering, indicating that this dominant factor affecting optical signals may be derived from a clinical radiograph. In the event that this was shown to be feasible, this would reduce or eliminate the need for accurate frequency- or time-domain NIRS measurement systems required for recovering optical scattering properties. Simpler, less expensive continuous-wave systems would be adequate to quantify absorbing and/or luminescent NIR signals.

The use of explicit biomarker prior information was recently demonstrated with a hybrid MRI–NIRS system [84]. This study demonstrated that incorporating fat/water values estimated with MRI fat–water imaging into the NIRS image recovery algorithm improved the accuracy of haemoglobin and oxygen saturation fitting. This hybrid information combination has been largely untested but will probably have increasing relevance as systems are developed to have mutual, complementary information sets. Additional data that would be synergistic with NIRS include blood haematocrit level and concentration of tracers in the blood.

Risk prediction of breast cancer is a recent example of how the intersection of patient demographics with NIRS has the potential for clinical impact. Lilge and colleagues (e.g. [85]) have systematically examined the potential for NIRS to quantify the parenchymal density of the breast. The density is a known correlate to the risk of incidence of breast cancer, and so NIRS is being examined as a research tool to help understand the nature of the correlation, and to determine if NIRS measurement systems could provide a low-cost screening platform for population risk analysis. In this paradigm, NIRS may be used to test and analyse the demographic community with a known genetic predisposition to high incidence of cancer, particularly younger women in whom higher numbers of mammograms might have adverse effects [86]. In addition to parenchymal density, biochemical and structural changes related to NIRS data have been shown to be associated with the risk of cancer [87], further indicating that NIRS data combined with patient demographic information might provide diagnostic information about breast cancer risk. This work is still ongoing, but serves as a key example of the potential synergy that might exist with the right application.

(b) Summary of benefits and limits

Strong correlations exist between NIRS data and demographics. Datasets from different diagnostic modalities may be exchanged or used to reduce uncertainty in quantification, when applied a priori rather than analysed posterior. The prior implementation of demographic or biomarker data is largely undeveloped in the NIRS prior implementation.

5. Prior information: biophysical models

(a) Biophysical models

Biophysical models to interpret the NIRS data have been used widely for a range of applications. While these could also be incorporated into the NIRS process at earlier steps, it is most common at the post-processing or interpretation stage. Examples of biophysical models are listed in table 3. Tissue oxygen extraction models have perhaps been most widely incorporated into NIRS endogenous imaging of oxygen saturation because of the potential to estimate tissue consumption of oxygen in key organs such as the brain or in diseases such as cancerous tumours. Generally, biophysical model information that could be included is in three different temporal domains: (i) static information (i.e. fat, water, blood volumes), (ii) harmonic information (i.e. blood pulsation or breathing or hormonal cycling), or (iii) dynamic inflow/outflow data (i.e. contrast agent injection and clearance). The latter two are illustrated in figure 6.

Table 3.
Some of the major biophysical models used in conjunction with NIRS, with applications and key references.
Figure 6.
The type of data collected for transmitted or emitted signals are illustrated as either (a) harmonic or (b) dynamic transient. If these patterns are known and can be modelled with an explicit numerical function, then it can be incorporated into the prefiltering, ...

Post-processing of fluorescence data has been widely analysed to study pharmacokinetics in vivo as a way to improve data interpretation and extract additional information from the measured signals. Temporal fitting can be used with compartmental models to estimate the pharmacokinetics of injected dyes [101,123,124], drug production or bleaching rates [125,126], or drug-binding rates in vivo [118]. Temporal analysis was used to illustrate how fluorescence signals can identify different organs in vivo by Hillman et al. [102]. This is a post-processing approach to include prior information about the molecules. Incorporation of the drug kinetics in the forward or inverse photon propagation models has been shown in just one or two studies [95], and there is a promise that including an accurate kinetic model in the inversion problem may improve accuracy in quantification. More complex kinetic or tissue modelling procedures have been studied by some groups [103]. Additionally, several studies have indicated that better knowledge of the drug concentration a priori would improve the ability to define the best acquisition parameters [127], indicating that inclusion in the data acquisition or system design level could also be beneficial.

Incorporating viscoelastic tissue modelling in NIRS is an approach that has seen some recent activity. Dynamic measurements induced by changes in pressure applied to the tissue can provide a rich dataset for robust, often ratiometric, interpretation of oxygen saturation and total haemoglobin values [106111]. These mechanical models are almost uniformly applied in post-processing or forward modelling [112], owing to the complexity of solving the viscoelastic problem in irregular domains. While they may introduce additional uncertainty in the inversion process, if implemented in the right way, they have the potential to reduce the ill-posedness of the problem and thus offer more accurate diagnostic performance.

Model-based interpretation of injected tracers is one area where there may be substantial gains, which would allow NIRS to provide fundamentally new information that cannot be obtained by other modalities. The key factor in this is that NIRS imaging allows quantification of more than one analyte in vivo so multiple species can be used to reference against one another. This has been shown for vascular permeability studies [119] and for drug-binding studies [118]. In this type of analysis, the binding rate constant of the targeted probe is directly related to the drug affinity and the concentration of binding sites within the tissue. Thus, by combining NIRS measurements or images with compartment-binding modelling, it is feasible to quantify molecular receptor expression and affinity in vivo. The processing of this information has been in post-processing to date, but explicit incorporation of these linear models into the inversion algorithm with a temporal dataset is readily done. An example of explicit incorporation of a forward model is in the way spectral constraints are added into a multiple-wavelength inversion problem, to reduce the ill-posedness of the inversion and provide more accurate recovery properties. A similar approach could be tried with forward physical models, in the inversion problem.

(b) Summary of benefits and limits

Proven benefits have been widespread in post-processing of NIRS data. Inclusion of biophysical modelling into forward NIRS modelling is probably somewhat beneficial. However, similar to spectral constraints, there is probably a major benefit in inclusion into the inverse algorithms.

6. Prior information: impact on diagnosis

Ultimately, any convolution of prior information data streams with NIRS must be assessed for diagnostic performance. This is a critical component of the innovation process often missing from NIRS studies. The prior–NIRS convolution must be assessed against other convolution schemes as well as against clinical and pre-clinical diagnostic standards. An important consideration in this analysis is the trade-off between performance and technical complexity. The value of modest gains in diagnostic power realized through extensive, difficult and time-consuming imaging strategies must be considered carefully.

Figure 7 illustrates how a receiver operating characteristic (ROC) curve can be used to compare the diagnostic performance of different prior implementations that can be used for MRI-guided NIRS imaging. In this example, a coronal MRI image of a human breast lightly compressed between two plates (figure 7a) was segmented and used to produce a finite element mesh containing two distinct tissue regions, an adipose and a fibroglandular region. Each compression plate is lined with eight optical source–detector fibres for transmitting and receiving light through the tissue, as illustrated by grey rectangles in figure 7b. The ROC curve analysis began by generating simulated test domains from the finite element mesh. A suspicious circular legion was numerically added to the mesh in one of two locations, either near the centre of the tissue, as shown in figure 7b, or a specified distance from one edge, as illustrated in figure 7c. The lesion simulates a region of gadolinium (Gd) enhancement in the MRI image, but may not necessarily be malignant. The radius of this lesion was varied over the range illustrated in figure 7b,c, and the optical absorption contrast was varied from 0 to 5 : 1 for each lesion size. Lesions with zero contrast in optical absorption represented an MRI false positive. For each lesion size/contrast, forward data were generated numerically and noise added to the measurements. A total of 600 datasets were considered in this example.

Figure 7.
(a) A T1-weighted coronal MRI image of a breast compressed between two plates. Tissue types were segmented and used to generate a two-tissue-region finite element mesh for NIRS modelling as shown in (b,c). Sixteen fibre probes contact the tissue surface ...

The synthetic optical data were then used to recover and interpret absorption images using one of two prior-imaging convolutions. The first method applied the structural prior in the interpretation step, while the second approach considered the hard-prior method, which incorporated the MRI information directly in the image reconstruction step. In the former approach, the optical data were used to perform a ‘no-priors’ image formation. Once the image was recovered, the average value in the region defined by the lesion in the corresponding MRI image was extracted as the diagnostic parameter. The hard-prior approach, on the other hand, encodes the internal structure of the tissue in optically homogeneous regions. Thus, the value recovered in this lesion was extracted directly from the image.

The value of optical absorption in the lesion region was used as the diagnostic parameter in the ROC curve analysis for both methods shown in figure 7e. This result demonstrates an improvement in ROC curve performance when the prior was applied in the reconstruction step as opposed to the interpretation step. Area-under-the-curve metrics were 0.85 and 0.73 for the hard-prior and interpretation-only approaches, respectively. However, this analysis assumed perfect image segmentation of the fibroglandular region. Errors in segmenting this region will affect the performance of the hard-prior approach but should not affect the approach that uses the MRI information only in the interpretation step. To explore to what degree the performance is degraded, the segmentation of the fibroglandular region used in the hard-prior image reconstruction was dilated, as shown in figure 7d, and the analysis repeated. ROC curves for the two approaches under imperfect segmentation conditions are shown in figure 7f, and demonstrate a significant decrease in diagnostic performance for the hard-prior approach due to the segmentation error. Analyses such as this can help determine what imaging strategies are most beneficial under different experimental conditions.

The strategy described above is an effective way to evaluate and compare different implementations of prior information. ROC curve analysis is the preferred method for assessing diagnostic performance in most cases (though specific applications, such as glucose monitoring, have their own standard diagnostic scales). Other types of analyses are commonly reported for NIRS imaging, such as spatial resolution, contrast resolution and contrast detail analysis. While useful for assessing system performance, these analyses should be considered an intermediate step on the path towards a full diagnostic assessment using an ROC curve, which also accounts for biological variability. Establishing truth can be challenging and must be accomplished using accurate gold standards of conventional medicine. Also, the metrics assessed should match the objectives of the diagnostic test. In some cases, the aim of the NIRS technique is to quantify rather than detect diseased tissue, and thus an assessment of detectability would be inappropriate. This is especially true for paradigms in which structural imaging modalities are used to identify suspicious lesions and a prior–NIRS convolution strategy is used to characterize and diagnose the abnormalities. In most cases, a comprehensive ROC curve study should be the ultimate goal of the NIRS researcher to determine whether adoption of the technique in the clinic or research setting is warranted.

7. Conclusions

The observations of improved quantification or improved diagnostic values are a complex range of niche examples. A few applications show clear and significant benefits, others demonstrate only modest improvements in diagnosis, and many have not been tested using the conventions of diagnostic assessment. Table 4 provides a listing of the prior information types from figure 1, along with the NIRS process steps. The table includes a qualitative ranking of which prior information streams have benefits in each of the NIRS steps, using − to indicate no clear benefit, + to indicate possible benefit, ++ to indicate probable benefit with some demonstrated evidence, and +++ to indicate a major demonstrated benefit in quantification or diagnosis.

Table 4.
Listing of the prior information options and NIRS process steps outlined in figure 1, with relative scoring of the proven benefit, where notation is as follows: −, no benefit; +, possible benefit; ++, demonstrated some benefit; and +++, ...

While these scores are admittedly relative and somewhat subjective, this chart is an initial attempt to think about this framework for determining how the benefit from prior information might be maximized in the NIRS process, and how systems can be designed to maximize quantitative accuracy and ultimately diagnostic value in terms of the sensitivity and specificity trade-off.

Acknowledgements

This work was funded by the National Cancer Institute research grants RO1CA069544, RO1CA109558, K25CA138578 and R01CA120368.

References

1. Stott J. J., Culver J. P., Arridge S. R., Boas D. A. 2003. Optode positional calibration in diffuse optical tomography Appl. Opt. 423154–3162.3162 (doi:10.1364/AO.42.003154) doi: 10.1364/AO.42.003154. [PubMed] [Cross Ref]
2. Schweiger M., Nissila I., Boas D. A., Arridge S. R. 2007. Image reconstruction in optical tomography in the presence of coupling errors Appl. Opt. 462743–2756.2756 (doi:10.1364/AO.46.002743) doi: 10.1364/AO.46.002743. [PubMed] [Cross Ref]
3. Pogue B., McBride T., Osterberg U., Paulsen K. 1999. Comparison of imaging geometries for diffuse optical tomography of tissue Opt. Express 4270–286.286 (doi:10.1364/OE.4.000270) doi: 10.1364/OE.4.000270. [PubMed] [Cross Ref]
4. Pogue B. W., Song X., Tosteson T. D., McBride T. O., Jiang S., Paulsen K. D. 2002. Statistical analysis of nonlinearly reconstructed near-infrared tomographic images. I. Theory and simulations IEEE Trans. Med. Imaging 21755–763.763 (doi:10.1109/TMI.2002.801155) doi: 10.1109/TMI.2002.801155. [PubMed] [Cross Ref]
5. Song X., Pogue B. W., Tosteson T. D., McBride T. O., Jiang S., Paulsen K. D. 2002. Statistical analysis of nonlinearly reconstructed near-infrared tomographic images. II. Experimental interpretation IEEE Trans. Med. Imaging 21764–772.772 (doi:10.1109/TMI.2002.801158) doi: 10.1109/TMI.2002.801158. [PubMed] [Cross Ref]
6. Noponen T. E. J., Kotilahti K., Nissila I., Kajava T., Merilainen P. T. 2010. Effects of improper source coupling in frequency-domain near-infrared spectroscopy Phys. Med. Biol. 552941–2960.2960 (doi:10.1088/0031-9155/55/10/010) doi: 10.1088/0031-9155/55/10/010. [PubMed] [Cross Ref]
7. Aoyagi T. 2003. Pulse oximetry: its invention, theory, and future J. Anesth. 17259–266.266 (doi:10.1007/s00540-003-0192-6) doi: 10.1007/s00540-003-0192-6. [PubMed] [Cross Ref]
8. Jobsis F. F. 1977. Noninvasive, infrared monitoring of cerebral and myocardial oxygen sufficiency and circulatory parameters Science 1981264–1267.1267 (doi:10.1126/science.929199) doi: 10.1126/science.929199. [PubMed] [Cross Ref]
9. Huppert T. J., Diamond S. G., Franceschini M. A., Boas D. A. 2009. HomER: a review of time-series analysis methods for near-infrared spectroscopy of the brain Appl. Opt. 48D280–D298.D298 (doi:10.1364/AO.48.00D280) doi: 10.1364/AO.48.00D280. [PMC free article] [PubMed] [Cross Ref]
10. Quaresima V., Springett R., Cope M., Wyatt J. T., Delpy D. T., Ferrari M., Cooper C. E. 1998. Oxidation and reduction of cytochrome oxidase in the neonatal brain observed by in vivo near-infrared spectroscopy Biochim. Biophys. Acta 1366291–300.300 (doi:10.1016/S0005-2728(98)00129-7) doi: 10.1016/S0005-2728(98)00129-7. [PubMed] [Cross Ref]
11. Arakaki L. S., Burns D. H., Kushmerick M. J. 2007. Accurate myoglobin oxygen saturation by optical spectroscopy measured in blood-perfused rat muscle Appl. Spectrosc. 61978–985.985 (doi:10.1366/000370207781745928) doi: 10.1366/000370207781745928. [PubMed] [Cross Ref]
12. Dehghani H., Leblond F., Pogue B. W., Chauchard F. 2010. Application of spectral derivative data in visible and near-infrared spectroscopy Phys. Med. Biol. 553381–3399.3399 (doi:10.1088/0031-9155/55/12/008) doi: 10.1088/0031-9155/55/12/008. [PMC free article] [PubMed] [Cross Ref]
13. Ntziachristos V., Turner G., Dunham J., Windsor S., Soubret A., Ripoll J., Shih H. A. 2005. Planar fluorescence imaging using normalized data J. Biomed. Opt. 10064007. (doi:10.1117/1.2136148) doi: 10.1117/1.2136148. [PubMed] [Cross Ref]
14. Vinegoni C., Razansky D., Figueiredo J. L., Nahrendorf M., Ntziachristos V., Weissleder R. 2009. Normalized Born ratio for fluorescence optical projection tomography Opt. Lett. 34319–321.321 (doi:10.1364/OL.34.000319) doi: 10.1364/OL.34.000319. [PMC free article] [PubMed] [Cross Ref]
15. Leblond F., et al. Analytic expression of fluorescence-ratio detection provides an estimate of depth in multi-spectral sub-surface imaging. Phys. Med. Biol. In press [PubMed]
16. Fantini S., Franceschini M. A., Gratton E. 1997. Effective source term in the diffusion equation for photon transport in turbid media Appl. Opt. 36156–163.163 (doi:10.1364/AO.36.000156) doi: 10.1364/AO.36.000156. [PubMed] [Cross Ref]
17. Franceschini M. A., Fantini S., Paunescu L. A., Maier J. S., Gratton E. 1998. Influence of a superficial layer in the quantitative spectroscopic study of strongly scattering media Appl. Opt. 377447–7458.7458 (doi:10.1364/AO.37.007447) doi: 10.1364/AO.37.007447. [PubMed] [Cross Ref]
18. Pogue B. W., Paulsen K. D., Abele C., Kaufman H. 2000. Calibration of near-infrared frequency-domain tissue spectroscopy for absolute absorption coefficient quantitation in neonatal head-simulating phantoms J. Biomed. Opt. 5185–193.193 (doi:10.1117/1.429985) doi: 10.1117/1.429985. [PubMed] [Cross Ref]
19. Amelink A., Sterenborg H. J. 2004. Measurement of the local optical properties of turbid media by differential path-length spectroscopy Appl. Opt. 433048–3054.3054 (doi:10.1364/AO.43.003048) doi: 10.1364/AO.43.003048. [PubMed] [Cross Ref]
20. Kanick S. C., Sterenborg H. J., Amelink A. 2009. Empirical model of the photon path length for a single fiber reflectance spectroscopy device Opt. Express 17860–871.871 (doi:10.1364/OE.17.000860) doi: 10.1364/OE.17.000860. [PubMed] [Cross Ref]
21. Freeberg J. A., et al. 2007. Fluorescence and reflectance device variability throughout the progression of a phase II clinical trial to detect and screen for cervical neoplasia using a fiber optic probe J. Biomed. Opt. 12034015. (doi:10.1117/1.2750332) doi: 10.1117/1.2750332. [PubMed] [Cross Ref]
22. Myakov A., Nieman L., Wicky L., Utzinger U., Richards-Kortum R., Sokolov K. 2002. Fiber optic probe for polarized reflectance spectroscopy in vivo: design and performance J. Biomed. Opt. 7388–397.397 (doi:10.1117/1.1483314) doi: 10.1117/1.1483314. [PubMed] [Cross Ref]
23. Mourant J. R., Bigio I. J., Jack D. A., Johnson T. M., Miller H. D. 1997. Measuring absorption coefficients in small volumes of highly scattering media: source–detector separations for which path lengths do not depend on scattering properties Appl. Opt. 365655–5661.5661 (doi:10.1364/AO.36.005655) doi: 10.1364/AO.36.005655. [PubMed] [Cross Ref]
24. Bigio I. J., Mourant J. R., Los G. Noninvasive, in-situ measurement of drug concentrations in tissue using optical spectroscopy. J. Gravit. Physiol. 1999;6:P173–P175. See http://www.ncbi.nlm.nih.gov/pubmed/11543008 . [PubMed]
25. Pogue B. W., Burke G. 1998. Fiber-optic bundle design for quantitative fluorescence measurement from tissue Appl. Opt. 377429–7436.7436 (doi:10.1364/AO.37.007429) doi: 10.1364/AO.37.007429. [PubMed] [Cross Ref]
26. Pogue B. W., Chen B., Zhou X., Hoopes P. J. 2005. Analysis of sampling volume and tissue heterogeneity on the in vivo detection of fluorescence J. Biomed. Opt. 10041206. (doi:10.1117/1.2002978) doi: 10.1117/1.2002978. [PubMed] [Cross Ref]
27. Edwards A. D., Wyatt J. S., Richardson C., Delpy D. T., Cope M., Reynolds E. O. 1988. Cotside measurement of cerebral blood flow in ill newborn infants by near infrared spectroscopy Lancet 2770–771.771 (doi:10.1016/S0140-6736(88)92418-X) doi: 10.1016/S0140-6736(88)92418-X. [PubMed] [Cross Ref]
28. Elwell C. E., Owen-Reece H., Wyatt J. S., Cope M., Reynolds E. O., Delpy D. T. 1996. Influence of respiration and changes in expiratory pressure on cerebral haemoglobin concentration measured by near infrared spectroscopy J. Cereb. Blood Flow Metab. 16353–357.357 (doi:10.1097/00004647-199603000-00022) doi: 10.1097/00004647-199603000-00022. [PubMed] [Cross Ref]
29. Redmond R. W., Kochevar I. E. 2006. Spatially resolved cellular responses to singlet oxygen Photochem. Photobiol. 821178–1186.1186 (doi:10.1562/2006-04-14-IR-874) doi: 10.1562/2006-04-14-IR-874. [PubMed] [Cross Ref]
30. Aveline B. M., Redmond R. W. 1999. Can cellular phototoxicity be accurately predicted on the basis of sensitizer photophysics? Photochem. Photobiol. 69306–316.316 (doi:10.1111/j.1751-1097.1999.tb03291.x) doi: 10.1111/j.1751-1097.1999.tb03291.x. [PubMed] [Cross Ref]
31. Andersson-Engels S., Berg R., Svanberg S., Jarlman O. 1990. Time-resolved transillumination for medical diagnostics Opt. Lett. 151179–1181.1181 (doi:10.1364/OL.15.001179) doi: 10.1364/OL.15.001179. [PubMed] [Cross Ref]
32. Turner G. M., Soubret A., Ntziachristos V. 2007. Inversion with early photons Med. Phys. 341405–1411.1411 (doi:10.1118/1.2437103) doi: 10.1118/1.2437103. [PubMed] [Cross Ref]
33. Pifferi A., et al. 2008. Time-resolved diffuse reflectance using small source–detector separation and fast single-photon gating Phys. Rev. Lett. 100138101. (doi:10.1103/PhysRevLett.100.138101) doi: 10.1103/PhysRevLett.100.138101. [PubMed] [Cross Ref]
34. Leblond F., Dehghani H., Kepshire D., Pogue B. W. 2009. Early-photon fluorescence tomography: spatial resolution improvements and noise stability considerations J. Opt. Soc. Am. A 261444–1457.1457 (doi:10.1364/JOSAA.26.001444) doi: 10.1364/JOSAA.26.001444. [PubMed] [Cross Ref]
35. Schweiger M., Arridge S. R. 1999. Application of temporal filters to time resolved data in optical tomography Phys. Med. Biol. 441699–1717.1717 (doi:10.1088/0031-9155/44/7/310) doi: 10.1088/0031-9155/44/7/310. [PubMed] [Cross Ref]
36. Madsen S. J., Anderson E. R., Haskell R. C., Tromberg B. J. 1994. Portable, high-bandwidth frequency-domain photon migration instrument for tissue spectroscopy Opt. Lett. 191934–1936.1936 (doi:10.1364/OL.19.001934) doi: 10.1364/OL.19.001934. [PubMed] [Cross Ref]
37. Bevilacqua F., Berger A. J., Cerussi A. E., Jakubowski D., Tromberg B. J. 2000. Broadband absorption spectroscopy in turbid media by combined frequency-domain and steady-state methods Appl. Opt. 396498–6507.6507 (doi:10.1364/AO.39.006498) doi: 10.1364/AO.39.006498. [PubMed] [Cross Ref]
38. Akers W., Lesage F., Holten D., Achilefu S. In vivo resolution of multiexponential decays of multiple near-infrared molecular probes by fluorescence lifetime-gated whole-body time-resolved diffuse optical imaging. Mol. Imaging. 2007;6:237–246. See http://www.ncbi.nlm.nih.gov/pubmed/17711779 . [PubMed]
39. Gratton E., Fantini S., Franceschini M. A., Gratton G., Fabiani M. 1997. Measurements of scattering and absorption changes in muscle and brain Phil. Trans. R. Soc. Lond. B 352727–735.735 (doi:10.1098/rstb.1997.0055) doi: 10.1098/rstb.1997.0055. [PMC free article] [PubMed] [Cross Ref]
40. Casavola C., Paunescu L. A., Fantini S., Franceschini M. A., Lugara P. M., Gratton E. Application of near-infrared tissue oxymetry to the diagnosis of peripheral vascular disease. Clin. Hemorheol. Microcirc. 1999;21:389–393. See http://www.ncbi.nlm.nih.gov/pubmed/10711775 . [PubMed]
41. Casavola C., Paunescu L. A., Fantini S., Gratton E. 2000. Blood flow and oxygen consumption with near-infrared spectroscopy and venous occlusion: spatial maps and the effect of time and pressure of inflation J. Biomed. Opt. 5269–276.276 (doi:10.1117/1.429995) doi: 10.1117/1.429995. [PubMed] [Cross Ref]
42. Franceschini M. A., Toronov V., Filiaci M., Gratton E., Fantini S. 2000. On-line optical imaging of the human brain with 160-ms temporal resolution Opt. Express 649–57.57 (doi:10.1364/OE.6.000049) doi: 10.1364/OE.6.000049. [PubMed] [Cross Ref]
43. Arridge S. R., Cope M., Delpy D. T. 1992. The theoretical basis for the determination of optical pathlengths in tissue: temporal and frequency analysis Phys. Med. Biol. 371531–1560.1560 (doi:10.1088/0031-9155/37/7/005) doi: 10.1088/0031-9155/37/7/005. [PubMed] [Cross Ref]
44. Pineda A. R., Schweiger M., Arridge S. R., Barrett H. H. 2006. Information content of data types in time-domain optical tomography J. Opt. Soc. Am. A 232989–2996.2996 (doi:10.1364/JOSAA.23.002989) doi: 10.1364/JOSAA.23.002989. [PMC free article] [PubMed] [Cross Ref]
45. Ntziachristos V., Yodh A. G., Schnall M., Chance B. 2000. Concurrent MRI and diffuse optical tomography of breast after indocyanine green enhancement Proc. Natl Acad. Sci. USA 972767–2772.2772 (doi:10.1073/pnas.040570597) doi: 10.1073/pnas.040570597. [PubMed] [Cross Ref]
46. Moon J. A., Mahon R., Duncan M. D., Reintjes J. 1993. Resolution limits for imaging through turbid media with diffuse light Opt. Lett. 181591–1593.1593 (doi:10.1364/OL.18.001591) doi: 10.1364/OL.18.001591. [PubMed] [Cross Ref]
47. Boas D. A., Chen K., Grebert D., Franceschini M. A. 2004. Improving the diffuse optical imaging spatial resolution of the cerebral hemodynamic response to brain activation in humans Opt. Lett. 291506–1508.1508 (doi:10.1364/OL.29.001506) doi: 10.1364/OL.29.001506. [PubMed] [Cross Ref]
48. Fang Q., Boas D. A. 2009. Monte Carlo simulation of photon migration in 3D turbid media accelerated by graphics processing units Opt. Express 1720178–20190.20190 (doi:10.1364/OE.17.020178) doi: 10.1364/OE.17.020178. [PMC free article] [PubMed] [Cross Ref]
49. Alerstam E., Andersson-Engels S., Svensson T. 2008. White Monte Carlo for time-resolved photon migration J. Biomed. Opt. 13041304. (doi:10.1117/1.2950319) doi: 10.1117/1.2950319. [PubMed] [Cross Ref]
50. Kienle A., Patterson M. S. 1996. Determination of the optical properties of turbid media from a single Monte Carlo simulation Phys. Med. Biol. 412221–2227.2227 (doi:10.1088/0031-9155/41/10/026) doi: 10.1088/0031-9155/41/10/026. [PubMed] [Cross Ref]
51. Alerstam E., Svensson T., Andersson-Engels S. 2008. Parallel computing with graphics processing units for high-speed Monte Carlo simulation of photon migration J. Biomed. Opt. 13060504. (doi:10.1117/1.3041496) doi: 10.1117/1.3041496. [PubMed] [Cross Ref]
52. Schweiger M., Zhukov L., Arridge S., Johnson C. 1999. Optical tomography using the SCIRun problem solving environment: preliminary results for three-dimensional geometries and parallel processing Opt. Express 4263–269.269 (doi:10.1364/OE.4.000263) doi: 10.1364/OE.4.000263. [PubMed] [Cross Ref]
53. Dehghani H., Eames M. E., Yalavarthy P. K., Davis S. C., Srinivasan S., Carpenter C. M., Pogue B. W., Paulsen K. D. 2008. Near infrared optical tomography using NIRFAST: algorithm for numerical model and image reconstruction Commun. Numer. Meth. Eng. 25711–732.732 (doi:10.1002/cnm.1162) doi: 10.1002/cnm.1162. [PMC free article] [PubMed] [Cross Ref]
54. Brooksby B., Jiang S. D., Dehghani H., Pogue B. W., Paulsen K. D., Weaver J., Kogel C., Poplack S. P. 2005. Combining near-infrared tomography resonance imaging to study in vivo and magnetic breast tissue: implementation of a Laplacian-type regularization to incorporate magnetic resonance structure J. Biomed. Opt. 10051504. (doi:10.1117/1.2098627) doi: 10.1117/1.2098627. [PubMed] [Cross Ref]
55. Yalavarthy P. K., Pogue B. W., Dehghani H., Carpenter C. M., Jiang S., Paulsen K. D. 2007. Structural information within regularization matrices improves near infrared diffuse optical tomography Opt. Express 158043–8058.8058 (doi:10.1364/OE.15.008043) doi: 10.1364/OE.15.008043. [PubMed] [Cross Ref]
56. Brooksby B., et al. 2006. Imaging breast adipose and fibroglandular tissue molecular signatures by using hybrid MRI-guided near-infrared spectral tomography Proc. Natl Acad. Sci. USA 1038828–8833.8833 (doi:10.1073/pnas.0509636103) doi: 10.1073/pnas.0509636103. [PubMed] [Cross Ref]
57. Carpenter C. M., et al. 2007. Image-guided optical spectroscopy provides molecular-specific information in vivo: MRI-guided spectroscopy of breast cancer hemoglobin, water, and scatterer size Opt. Lett. 32933–935.935 (doi:10.1364/OL.32.000933) doi: 10.1364/OL.32.000933. [PubMed] [Cross Ref]
58. Axelsson J., Svensson J., Andersson-Engels S. 2007. Spatially varying regularization based on spectrally resolved fluorescence emission in fluorescence molecular tomography Opt. Express 1513574–13584.13584 (doi:10.1364/OE.15.013574) doi: 10.1364/OE.15.013574. [PubMed] [Cross Ref]
59. Lin Y., Barber W. C., Iwanczyk J. S., Roeck W., Nalcioglu O., Gulsen G. 2010. Quantitative fluorescence tomography using a combined tri-modality FT/DOT/XCT system Opt. Express 187835–7850.7850 (doi:10.1364/OE.18.007835) doi: 10.1364/OE.18.007835. [PMC free article] [PubMed] [Cross Ref]
60. Ale A., Schulz R. B., Sarantopoulos A., Ntziachristos V. 2010. Imaging performance of a hybrid X-ray computed tomography–fluorescence molecular tomography system using priors Med. Phys. 371976–1986.1986 (doi:10.1118/1.3368603) doi: 10.1118/1.3368603. [PubMed] [Cross Ref]
61. Intes X. 2005. Time-domain optical mammography SoftScan: initial results Acad. Radiol. 12934–947.947 (doi:10.1016/j.acra.2005.05.006) doi: 10.1016/j.acra.2005.05.006. [PubMed] [Cross Ref]
62. Azar F. S., Lee K., Khamene A., Choe R., Corlu A., Konecky S. D., Sauer F., Yodh A. G. 2007. Standardized platform for coregistration of nonconcurrent diffuse optical and magnetic resonance breast images obtained in different geometries J. Biomed. Opt. 12051902. (doi:10.1117/1.2798630) doi: 10.1117/1.2798630. [PubMed] [Cross Ref]
63. Chance B., Nioka S., Zhang J., Conant E. F., Hwang E., Briest S., Orel S. G., Schnall M. D., Czerniecki B. J. 2005. Breast cancer detection based on incremental biochemical and physiological properties of breast cancers: a six-year, two-site study Acad. Radiol. 12925–933.933 (doi:10.1016/j.acra.2005.04.016) doi: 10.1016/j.acra.2005.04.016. [PubMed] [Cross Ref]
64. Poplack S. P., et al. 2007. Electromagnetic breast imaging: results of a pilot study in women with abnormal mammograms Radiology 243350–359.359 (doi:10.1148/radiol.2432060286) doi: 10.1148/radiol.2432060286. [PubMed] [Cross Ref]
65. Fang Q., Selb J., Carp S. A., Boverman G., Kopans D. B., Moore R. H., Boas D. A. 2007. Proc. SPIE 643164310H. (doi:10.1117/12.701623) doi: 10.1117/12.701623. [Cross Ref]
66. Choe R., et al. 2005. Diffuse optical tomography of breast cancer during neoadjuvant chemotherapy: a case study with comparison to MRI Med. Phys. 321128–1139.1139 (doi:10.1118/1.1869612) doi: 10.1118/1.1869612. [PubMed] [Cross Ref]
67. Jiang S., et al. 2009. Evaluation of breast tumor response to neoadjuvant chemotherapy with tomographic diffuse optical spectroscopy: case studies of tumor region-of-interest changes Radiology 252551–560.560 (doi:10.1148/radiol.2522081202) doi: 10.1148/radiol.2522081202. [PubMed] [Cross Ref]
68. Corlu A., Durduran T., Choe R., Schweiger M., Hillman E. M., Arridge S. R., Yodh A. G. 2003. Uniqueness and wavelength optimization in continuous-wave multispectral diffuse optical tomography Opt. Lett. 282339–2341.2341 (doi:10.1364/OL.28.002339) doi: 10.1364/OL.28.002339. [PubMed] [Cross Ref]
69. Srinivasan S., Pogue B. W., Brooksby B., Jiang S. D., Dehghani H., Kogel C., Wells W. A., Poplack S. P., Paulsen K. D. Near-infrared characterization of breast tumors in vivo using spectrally-constrained reconstruction. Technol. Cancer Res. Treat. 2005;4:513–526. See http://www.ncbi.nlm.nih.gov/pubmed/16173822 . [PubMed]
70. D'Andrea C., Spinelli L., Bassi A., Giusto A., Contini D., Swartling J., Torricelli A., Cubeddu R. 2006. Time-resolved spectrally constrained method for the quantification of chromophore concentrations and scattering parameters in diffusing media Opt. Express 141888–1898.1898 (doi:10.1364/OE.14.001888) doi: 10.1364/OE.14.001888. [PubMed] [Cross Ref]
71. Brooksby B., Srinivasan S., Jiang S., Dehghani H., Pogue B. W., Paulsen K. D., Weaver J., Kogel C., Poplack S. P. 2005. Spectral priors improve near-infrared diffuse tomography more than spatial priors Opt. Lett. 301968–1970.1970 (doi:10.1364/OL.30.001968) doi: 10.1364/OL.30.001968. [PubMed] [Cross Ref]
72. Corlu A., Choe R., Durduran T., Lee K., Schweiger M., Arridge S. R., Hillman E. M., Yodh A. G. 2005. Diffuse optical tomography with spectral constraints and wavelength optimization Appl. Opt. 442082–2093.2093 (doi:10.1364/AO.44.002082) doi: 10.1364/AO.44.002082. [PubMed] [Cross Ref]
73. Eames M. E., Wang J., Pogue B. W., Dehghani H. 2008. Wavelength band optimization in spectral near-infrared optical tomography improves accuracy while reducing data acquisition and computational burden J. Biomed. Opt. 13054037. (doi:10.1117/1.2976425) doi: 10.1117/1.2976425. [PubMed] [Cross Ref]
74. Leblond F., Davis S. C., Valdes P. A., Pogue B. W. 2010. Pre-clinical whole-body fluorescence imaging: review of instruments, methods and applications J. Photochem. Photobiol. B 9877–94.94 (doi:10.1016/j.jphotobiol.2009.11.007) doi: 10.1016/j.jphotobiol.2009.11.007. [PubMed] [Cross Ref]
75. Pogue B. W., Patterson M. S. 1994. Frequency-domain optical absorption spectroscopy of finite tissue volumes using diffusion theory Phys. Med. Biol. 391157–1180.1180 (doi:10.1088/0031-9155/39/7/008) doi: 10.1088/0031-9155/39/7/008. [PubMed] [Cross Ref]
76. Durkin A. J., Richards-Kortum R. 1996. Comparison of methods to determine chromophore concentrations from fluorescence spectra of turbid samples Lasers Surg. Med. 1975–89.89 (doi:10.1002/(SICI)1096-9101(1996)19:1<75::AID-LSM9>3.0.CO;2-N) doi: 10.1002/(SICI)1096-9101(1996)19:1%3C75::AID-LSM9%3E3.0.CO;2-N. [PubMed] [Cross Ref]
77. Davis S. C., Pogue B. W., Tuttle S. B., Dehghani H., Paulsen K. D. 2009. Spectral distortion in diffuse molecular luminescence tomography in turbid media J. Appl. Phys. 105102024. (doi:10.1063/1.3116130) doi: 10.1063/1.3116130. [PubMed] [Cross Ref]
78. Dehghani H., Davis S. C., Pogue B. W. 2008. Spectrally resolved bioluminescence tomography using the reciprocity approach Med. Phys. 354863–4871.4871 (doi:10.1118/1.2982138) doi: 10.1118/1.2982138. [PubMed] [Cross Ref]
79. Dehghani H., Davis S. C., Jiang S., Pogue B. W., Paulsen K. D., Patterson M. S. 2006. Spectrally resolved bioluminescence optical tomography Opt. Lett. 31365–367.367 (doi:10.1364/OL.31.000365) doi: 10.1364/OL.31.000365. [PubMed] [Cross Ref]
80. Sterenborg H. J., Thomsen S., Jacques S. L., Motamedi M. 1995. In vivo autofluorescence of an unpigmented melanoma in mice: correlation of spectroscopic properties to microscopic structure Melanoma Res. 5211–216.216 (doi:10.1097/00008390-199508000-00002) doi: 10.1097/00008390-199508000-00002. [PubMed] [Cross Ref]
81. Hyde D. E., Farrell T. J., Patterson M. S., Wilson B. C. 2001. A diffusion theory model of spatially resolved fluorescence from depth-dependent fluorophore concentrations Phys. Med. Biol. 46369–383.383 (doi:10.1088/0031-9155/46/2/307) doi: 10.1088/0031-9155/46/2/307. [PubMed] [Cross Ref]
82. Swartling J., Pifferi A., Enejder A. M., Andersson-Engels S. 2003. Accelerated Monte Carlo models to simulate fluorescence spectra from layered tissues J. Opt. Soc. Am. A 20714–727.727 (doi:10.1364/JOSAA.20.000714) doi: 10.1364/JOSAA.20.000714. [PubMed] [Cross Ref]
83. Matcher S. J., Cope M., Delpy D. T. 1994. Use of the water absorption spectrum to quantify tissue chromophore concentration changes in near-infrared spectroscopy Phys. Med. Biol. 39177–196.196 (doi:10.1088/0031-9155/39/1/011) doi: 10.1088/0031-9155/39/1/011. [PubMed] [Cross Ref]
84. Carpenter C. M., Pogue B. W., Jiang S., Wang J., Hargreaves B. A., Rakow-Penner R., Daniel B. L., Paulsen K. D. 2011. MR water quantitative priors improves the accuracy of optical breast imaging IEEE Trans. Med. Imaging 30159–168.168 (doi:10.1109/TMI.2010.2071394) doi: 10.1109/TMI.2010.2071394. [PubMed] [Cross Ref]
85. Simick M. K., Lilge L. 2005. Optical transillumination spectroscopy to quantify parenchymal tissue density: an indicator for breast cancer risk Br. J. Radiol. 781009–1017.1017 (doi:10.1259/bjr/14696165) doi: 10.1259/bjr/14696165. [PubMed] [Cross Ref]
86. Blackmore K. M., Knight J. A., Jong R., Lilge L. 2007. Assessing breast tissue density by transillumination breast spectroscopy (TIBS): an intermediate indicator of cancer risk Br. J. Radiol. 80545–556.556 (doi:10.1259/bjr/26858614) doi: 10.1259/bjr/26858614. [PubMed] [Cross Ref]
87. Knight J. A., Blackmore K. M., Wong J., Tharmalingam S., Lilge L. 2010. Optical spectroscopy of the breast in premenopausal women reveals tissue variation with changes in age and parity Med. Phys. 37419–426.426 (doi:10.1118/1.3276737) doi: 10.1118/1.3276737. [PubMed] [Cross Ref]
88. Huang P., Chance B., Wang X., Kime R., Nioka S., Chance E. M. Modeling of oxygen diffusion and metabolism from capillary to muscle. Adv. Exp. Med. Biol. 2003;540:325–330. See http://www.springer.com/978-0-306-48035-5 . [PubMed]
89. Hamaoka T., et al. Muscle oxygen consumption at onset of exercise by near infrared spectroscopy in humans. Adv. Exp. Med. Biol. 2003;530:475–483. See http://www.springer.com/978-0-306-47774-4 . [PubMed]
90. Boas D. A., Strangman G., Culver J. P., Hoge R. D., Jasdzewski G., Poldrack R. A., Rosen B. R., Mandeville J. B. 2003. Can the cerebral metabolic rate of oxygen be estimated with near-infrared spectroscopy? Phys. Med. Biol. 482405–2418.2418 (doi:10.1088/0031-9155/48/15/311) doi: 10.1088/0031-9155/48/15/311. [PubMed] [Cross Ref]
91. Culver J. P., Durduran T., Furuya D., Cheung C., Greenberg J. H., Yodh A. G. 2003. Diffuse optical tomography of cerebral blood flow, oxygenation, and metabolism in rat during focal ischemia J. Cereb. Blood Flow Metab. 23911–924.924 (doi:10.1097/01.WCB.0000076703.71231.BB) doi: 10.1097/01.WCB.0000076703.71231.BB. [PubMed] [Cross Ref]
92. Dysart J. S., Patterson M. S. 2006. Photobleaching kinetics, photoproduct formation, and dose estimation during ALA induced PpIX PDT of MLL cells under well oxygenated and hypoxic conditions Photochem. Photobiol. Sci. 573–81.81 (doi:10.1039/b511807g) doi: 10.1039/b511807g. [PubMed] [Cross Ref]
93. Wang K. K., Mitra S., Foster T. H. 2007. A comprehensive mathematical model of microscopic dose deposition in photodynamic therapy Med. Phys. 34282–293.293 (doi:10.1118/1.2401041) doi: 10.1118/1.2401041. [PubMed] [Cross Ref]
94. Foster T. H., Murant R. S., Bryant R. G., Knox R. S., Gibson S. L., Hilf R. 1991. Oxygen consumption and diffusion effects in photodynamic therapy Radiat. Res. 126296–303.303 (doi:10.2307/3577919) doi: 10.2307/3577919. [PubMed] [Cross Ref]
95. Alacam B., Yazici B., Serdaroglu A., Intes X., Chance B., Nioka S. 2006. Reconstruction of spatially resolved pharmacokinetic rate images of fluorescence agents in FDOT Conf. Proc. IEEE Eng. Med. Biol. Soc. 15627–5630.5630 (doi:10.1109/IEMBS.2006.259382) doi: 10.1109/IEMBS.2006.259382. [PubMed] [Cross Ref]
96. Huppert T. J., Allen M. S., Benav H., Jones P. B., Boas D. A. 2007. A multicompartment vascular model for inferring baseline and functional changes in cerebral oxygen metabolism and arterial dilation J. Cereb. Blood Flow Metab. 271262–1279.1279 (doi:10.1038/sj.jcbfm.9600435) doi: 10.1038/sj.jcbfm.9600435. [PMC free article] [PubMed] [Cross Ref]
97. Elwell C. E., Cope M., Edwards A. D., Wyatt J. S., Delpy D. T., Reynolds E. O. Quantification of adult cerebral hemodynamics by near-infrared spectroscopy. J. Appl. Physiol. 1994;77:2753–2760. See http://www.ncbi.nlm.nih.gov/pubmed/7896617 . [PubMed]
98. Brown D. W., Picot P. A., Naeini J. G., Springett R., Delpy D. T., Lee T. Y. 2002. Quantitative near infrared spectroscopy measurement of cerebral hemodynamics in newborn piglets Pediatr. Res. 51564–570.570 (doi:10.1203/00006450-200205000-00004) doi: 10.1203/00006450-200205000-00004. [PubMed] [Cross Ref]
99. Diamond S. G., Huppert T. J., Kolehmainen V., Franceschini M. A., Kaipio J. P., Arridge S. R., Boas D. A. 2006. Dynamic physiological modeling for functional diffuse optical tomography Neuroimage 3088–101.101 (doi:10.1016/j.neuroimage.2005.09.016) doi: 10.1016/j.neuroimage.2005.09.016. [PMC free article] [PubMed] [Cross Ref]
100. Fang Q., Sakadzic S., Ruvinskaya L., Devor A., Dale A. M., Boas D. A. 2008. Oxygen advection and diffusion in a three-dimensional vascular anatomical network Opt. Express 1617530–17541.17541 (doi:10.1364/OE.16.017530) doi: 10.1364/OE.16.017530. [PMC free article] [PubMed] [Cross Ref]
101. Alacam B., Yazici B., Intes X., Nioka S., Chance B. 2008. Pharmacokinetic-rate images of indocyanine green for breast tumors using near-infrared optical methods Phys. Med. Biol. 53837–859.859 (doi:10.1088/0031-9155/53/4/002) doi: 10.1088/0031-9155/53/4/002. [PubMed] [Cross Ref]
102. Hillman E. M., Moore A. 2007. All-optical anatomical co-registration for molecular imaging of small animals using dynamic contrast Nat. Photonics 1526–530.530 (doi:10.1038/nphoton.2007.146) doi: 10.1038/nphoton.2007.146. [PMC free article] [PubMed] [Cross Ref]
103. Alacam B., Yazici B., Intes X., Chance B. 2006. Extended Kalman filtering for the modeling and analysis of ICG pharmacokinetics in cancerous tumors using NIR optical methods IEEE Trans. Biomed. Eng. 531861–1871.1871 (doi:10.1109/TBME.2006.881796) doi: 10.1109/TBME.2006.881796. [PubMed] [Cross Ref]
104. Pogue B. W. 2006. Near-infrared characterization of disease via vascular permeability probes Acad. Radiol. 131–3.3 (doi:10.1016/j.acra.2005.11.002) doi: 10.1016/j.acra.2005.11.002. [PubMed] [Cross Ref]
105. Unlu M. B., Birgul O., Gulsen G. 2008. A simulation study of the variability of indocyanine green kinetics and using structural a priori information in dynamic contrast enhanced diffuse optical tomography (DCE-DOT) Phys. Med. Biol. 533189–3200.3200 (doi:10.1088/0031-9155/53/12/008) doi: 10.1088/0031-9155/53/12/008. [PMC free article] [PubMed] [Cross Ref]
106. Nioka S., Wen S., Zhang J., Du J., Intes X., Zhao Z., Chance B. 2005. Simulation study of breast tissue hemodynamics during pressure perturbation Adv. Exp. Med. Biol. 56617–22.22 (doi:10.1007/0-387-26206-7_3) doi: 10.1007/0-387-26206-7_3. [PubMed] [Cross Ref]
107. Darling A. L., Yalavarthy P. K., Doyley M. M., Dehghani H., Pogue B. W. 2007. Interstitial fluid pressure in soft tissue as a result of an externally applied contact pressure Phys. Med. Biol. 524121–4136.4136 (doi:10.1088/0031-9155/52/14/007) doi: 10.1088/0031-9155/52/14/007. [PubMed] [Cross Ref]
108. Jiang S., Pogue B. W., Laughney A. M., Kogel C. A., Paulsen K. D. 2009. Measurement of pressure–displacement kinetics of hemoglobin in normal breast tissue with near-infrared spectral imaging Appl. Opt. 48D130–D136.D136 (doi:10.1364/AO.48.00D130) doi: 10.1364/AO.48.00D130. [PMC free article] [PubMed] [Cross Ref]
109. Carp S. A., Kauffman T., Fang Q., Rafferty E., Moore R., Kopans D., Boas D. 2006. Compression-induced changes in the physiological state of the breast as observed through frequency domain photon migration measurements J. Biomed. Opt. 11064016. (doi:10.1117/1.2397572) doi: 10.1117/1.2397572. [PubMed] [Cross Ref]
110. Boverman G., Fang Q., Carp S. A., Miller E. L., Brooks D. H., Selb J., Moore R. H., Kopans D. B., Boas D. A. 2007. Spatio-temporal imaging of the hemoglobin in the compressed breast with diffuse optical tomography Phys. Med. Biol. 523619–3641.3641 (doi:10.1088/0031-9155/52/12/018) doi: 10.1088/0031-9155/52/12/018. [PubMed] [Cross Ref]
111. Athanasiou A., Vanel D., Fournier L., Balleyguier C. 2007. Optical mammography: a new technique for visualizing breast lesions in women presenting non palpable BIRADS 4–5 imaging findings: preliminary results with radiologic–pathologic correlation Cancer Imaging 734–40.40 (doi:10.1102/1470-7330.2007.0006) doi: 10.1102/1470-7330.2007.0006. [PMC free article] [PubMed] [Cross Ref]
112. Dehghani H., Doyley M. M., Pogue B. W., Jiang S., Geng J., Paulsen K. D. 2004. Breast deformation modelling for image reconstruction in near infrared optical tomography Phys. Med. Biol. 491131–1145.1145 (doi:10.1088/0031-9155/49/7/004) doi: 10.1088/0031-9155/49/7/004. [PubMed] [Cross Ref]
113. Vo T. V., Hammer P. E., Hoimes M. L., Nadgir S., Fantini S. 2007. Mathematical model for the hemodynamic response to venous occlusion measured with near-infrared spectroscopy in the human forearm IEEE Trans. Biomed. Eng. 54573–584.584 (doi:10.1109/TBME.2006.890123) doi: 10.1109/TBME.2006.890123. [PubMed] [Cross Ref]
114. Leung T. S., Tisdall M. M., Tachtsidis I., Smith M., Delpy D. T., Elwell C. E. 2008. Cerebral tissue oxygen saturation calculated using low frequency haemoglobin oscillations measured by near infrared spectroscopy in adult ventilated patients Adv. Exp. Med. Biol. 614235–244.244 (doi:10.1007/978-0-387-74911-2_27) doi: 10.1007/978-0-387-74911-2_27. [PubMed] [Cross Ref]
115. Edgington L. E., Berger A. B., Blum G., Albrow V. E., Paulick M. G., Lineberry N., Bogyo M. 2009. Noninvasive optical imaging of apoptosis by caspase-targeted activity-based probes Nat. Med. 15967–973.973 (doi:10.1038/nm.1938) doi: 10.1038/nm.1938. [PMC free article] [PubMed] [Cross Ref]
116. Li W., Li F., Huang Q., Frederick B., Bao S., Li C. Y. 2008. Noninvasive imaging and quantification of epidermal growth factor receptor kinase activation in vivo Cancer Res. 684990–4997.4997 (doi:10.1158/0008-5472.CAN-07-5984) doi: 10.1158/0008-5472.CAN-07-5984. [PubMed] [Cross Ref]
117. Bremer C., Bredow S., Mahmood U., Weissleder R., Tung C. H. 2001. Optical imaging of matrix metalloproteinase-2 activity in tumors: feasibility study in a mouse model Radiology 221523–529.529 (doi:10.1148/radiol.2212010368) doi: 10.1148/radiol.2212010368. [PubMed] [Cross Ref]
118. Pogue B. W., Samkoe K. S., Hextrum S., O'Hara J. A., Jermyn M., Srinivasan S., Hasan T. 2010. Imaging targeted-agent binding in vivo with two probes J. Biomed. Opt. 15030513. (doi:10.1117/1.3449109) doi: 10.1117/1.3449109. [PubMed] [Cross Ref]
119. Fischer T., et al. 2006. Assessment of unspecific near-infrared dyes in laser-induced fluorescence imaging of experimental arthritis Acad. Radiol. 134–13.13 (doi:10.1016/j.acra.2005.07.010) doi: 10.1016/j.acra.2005.07.010. [PubMed] [Cross Ref]
120. Kiessling F., et al. 2009. RGD-labeled USPIO inhibits adhesion and endocytotic activity of αv β3-integrin-expressing glioma cells and only accumulates in the vascular tumor compartment Radiology 253462–469.469 (doi:10.1148/radiol.2532081815) doi: 10.1148/radiol.2532081815. [PubMed] [Cross Ref]
121. Thurber G. M., Wittrup K. D. 2008. Quantitative spatiotemporal analysis of antibody fragment diffusion and endocytic consumption in tumor spheroids Cancer Res. 683334–3341.3341 (doi:10.1158/0008-5472.CAN-07-3018) doi: 10.1158/0008-5472.CAN-07-3018. [PMC free article] [PubMed] [Cross Ref]
122. Thurber G. M., Schmidt M. M., Wittrup K. D. 2008. Factors determining antibody distribution in tumors Trends Pharmacol. Sci. 2957–61.61 (doi:10.1016/j.tips.2007.11.004) doi: 10.1016/j.tips.2007.11.004. [PMC free article] [PubMed] [Cross Ref]
123. Alacam B., Yazici B., Intes X., Chance B. Analysis of ICG pharmacokinetics in cancerous tumors using NIR optical methods. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2005;1:62–65. 10.1109/IEMBS.2005.1616342. [PubMed]
124. Qian H., Gu Y., Wang M., Achilefu S. 2009. Optimization of the near-infrared fluorescence labeling for in vivo monitoring of a protein drug distribution in animal model J. Fluoresc. 19277–284.284 (doi:10.1007/s10895-008-0413-3) doi: 10.1007/s10895-008-0413-3. [PubMed] [Cross Ref]
125. Vaidyanathan V. V., Rastegar S., Fossum T. W., Flores P., van der Breggen E. W., Egger N. G., Jacques S. L., Motamedi M. 2000. A study of aminolevulinic acid-induced protoporphyrin IX fluorescence kinetics in the canine oral cavity Lasers Surg. Med. 26405–414.414 (doi:10.1002/(SICI)1096-9101(2000)26:4<405::AID-LSM9>3.0.CO;2-E) doi: 10.1002/(SICI)1096-9101(2000)26:4%3C405::AID-LSM9%3E3.0.CO;2-E. [PubMed] [Cross Ref]
126. Star W. M., Aalders M. C., Sac A., Sterenborg H. J. 2002. Quantitative model calculation of the time-dependent protoporphyrin IX concentration in normal human epidermis after delivery of ALA by passive topical application or lontophoresis Photochem. Photobiol. 75424–432.432 (doi:10.1562/0031-8655(2002)0750424QMCOTT2.0.CO2) doi: 10.1562/0031-8655(2002)0750424QMCOTT2.0.CO2. [PubMed] [Cross Ref]
127. Davis S. C., Pogue B. W., Dehghani H., Paulsen K. D. 2009. Tissue drug concentration determines whether fluorescence or absorption measurements are more sensitive in diffuse optical tomography of exogenous contrast agents Appl. Opt. 48D262–D272.D272 (doi:10.1364/AO.48.00D262) doi: 10.1364/AO.48.00D262. [PMC free article] [PubMed] [Cross Ref]

Articles from Philosophical transactions. Series A, Mathematical, physical, and engineering sciences are provided here courtesy of The Royal Society