Search

The role of the lymphatics in the clearance of cerebrospinal fluid (CSF) from the brain has been implicated in multiple neurodegenerative conditions. In premature infants, intraventricular hemorrhage causes increased CSF production and, if clearance is impeded, hydrocephalus and severe developmental disabilities can result. In this work, we developed and deployed near-infrared fluorescence (NIRF) tomography and imaging to assess CSF ventricular dynamics and extracranial outflow in similarly sized, intact non-human primates (NHP) following microdose of indocyanine green (ICG) administered to the right lateral ventricle. Fluorescence optical tomography measurements were made by delivering ~10 mW of 785 nm light to the scalp by sequential illumination of 8 fiber optics and imaging the 830 nm emission light collected from 22 fibers using a gallium arsenide intensified, charge coupled device. Acquisition times were 16 seconds. Image reconstruction used the diffusion approximation and hard-priors obtained from MRI to enable dynamic mapping of ICG-laden CSF ventricular dynamics and drainage into the subarachnoid space (SAS) of NHPs. Subsequent, planar NIRF imaging of the scalp confirmed extracranial efflux into SAS and abdominal imaging showed ICG clearance through the hepatobiliary system. Necropsy confirmed imaging results and showed that deep cervical lymph nodes were the routes of extracranial CSF egress. The results confirm the ability to use trace doses of ICG to monitor ventricular CSF dynamics and extracranial outflow in NHP. The techniques may also be feasible for similarly-sized infants and children who may suffer impairment of CSF outflow due to intraventricular hemorrhage.

Although Blood Oxygenation Level Dependent (BOLD) functional MRI (fMRI) is widely used to examine brain function in adults, the need for general anesthesia limits its practical utility in infants and small children. Functional Near-Infrared Spectroscopy - Diffuse Optical Tomography (fNIRS-DOT) imaging promises to be an alternative brain network imaging technique. Yet current versions of continuous-wave fNIRS-DOT systems are restricted to the cortical surface measurements and do not probe deep structures that are frequently injured especially in premature infants. Herein we report a transcranial near infrared optical imaging system, called Cap-based Transcranial Optical Tomography (CTOT) able to image whole brain hemodynamic activity with 3 seconds of data acquisition time. We show the system is capable of whole brain oxygenation mapping in an awake child, and that tomographically reconstructed static CTOT-derived oxy- and deoxygenated blood volumes are spatially correlated with the time-averaged BOLD fMRI volumes. By removing time bottlenecks in the current system, dynamic CTOT mapping should be possible, which would then enable evaluation of functional connectivity in awake infants.

Recently, fluorescence imaging using shortwave infrared light (SWIR, 1,000-2,000 nm) has been proposed as having advantage over conventional near-infrared fluorescence (NIRF) imaging due to the reduced tissue scattering, negligible autofluorescence, comparable tissue absorption, and the discovery that indocyanine green (ICG), used clinically as a NIRF contrast agent, also has fluorescence emission in SWIR regime. Images of ICG in small animals acquired by commercial Si-based and InGaAs-based imaging cameras have been qualitatively compared, however the lack of working standards to quantify performance of these imaging systems limits quantitative comparison. Without quantification using a traceable in vitro test, clinical adoption of rapidly evolving advances in both NIRF and SWIR imaging devices will become limited. In this work, we developed an ICG based fluorescent solid working standard calibrated with SI units (mW [Formula: see text]cm [Formula: see text]sr -1 ) for quantification of measurement sensitivity of Si, GaAs-intensified Si, and InGaAs based camera systems, their signal-to-noise ratio (SNR), and contrast in non-clinical tests. In addition, we present small animal and large animal imaging with ICG for qualitative comparison of the same SWIR fluorescence and NIRF imaging systems. Results suggest that SWIR fluorescence imaging of ICG may have superior resolution in small animal imaging compared to NIRF imaging, but lack of measurement sensitivity, SNR, contrast, as well as water absorption limits deep penetration in large animals.

To date, no emerging preclinical or clinical near-infrared fluorescence (NIRF) imaging devices for noninvasive and/or surgical guidance have their performances validated on working standards with SI units of radiance that enable comparison or quantitative quality assurance. In this work, we developed and deployed a methodology to calibrate a stable, solid phantom for emission radiance with International System of Units (SI) units of mW ·sr(-1) ·cm(-2) for use in characterizing the measurement sensitivity of ICCD and IsCMOS detection, signal-to-noise ratio, and contrast. In addition, at calibrated radiances, we assess transverse and lateral resolution of ICCD and IsCMOS camera systems. The methodology allowed demonstration of superior SNR of the ICCD over the IsCMOS technology and superior resolution of the IsCMOS over the ICCD approach. Contrast depended upon the camera settings (binning and integration time) and gain of intensifier. Finally, because the architecture of CMOS and CCD camera systems results in vastly different performance, we comment on the utility of these technologies for small animal imaging as well as clinical applications for noninvasive and surgical guidance.

The performance evaluation of a variety of small animal tomography measurement approaches and algorithms for recovery of fluorescent absorption cross section has not been conducted. Herein, we employed an intensified CCD system installed in a commercial small animal CT (Computed Tomography) scanner to compare image reconstructions from time-independent, continuous wave (CW) measurements and from time-dependent, frequency domain (FD) measurements in a series of physical phantoms specifically designed for evaluation. Comparisons were performed as a function of (1) number of projections, (2) the level of preprocessing filters used to improve the signal-to-noise ratio (SNR), (3) endogenous heterogeneity of optical properties, as well as in the cases of (4) two fluorescent targets and (5) a mouse-shaped phantom. Assessment of quantitative recovery of fluorescence absorption cross section was performed using a fully parallel, regularization-free, linear reconstruction algorithm with diffusion approximation (DA) and high order simplified spherical harmonics ( SPN) approximation to the radiative transport equation (RTE). The results show that while FD measurements may result in superior image reconstructions over CW measurements, data acquisition times are significantly longer, necessitating further development of multiple detector/source configurations, improved data read-out rates, and detector technology. FD measurements with SP3 reconstructions enabled better quantitative recovery of fluorescent target strength, but required increased computational expense. Despite the developed parallel reconstruction framework being able to achieve more than 60 times speed increase over sequential implementation, further development in faster parallel acceleration strategies for near-real time and real-time image recovery and more precise forward solution is necessary.

A novel image reconstruction algorithm has been developed and demonstrated for fluorescence-enhanced frequency-domain photon migration (FDPM) tomography from measurements of area illumination with modulated excitation light and area collection of emitted fluorescence light using a gain modulated image-intensified charge-coupled device (ICCD) camera. The image reconstruction problem was formulated as a nonlinear least-squares-type simple bounds constrained optimization problem based upon the penalty/modified barrier function (PMBF) method and the coupled diffusion equations. The simple bounds constraints are included in the objective function of the PMBF method and the gradient-based truncated Newton method with trust region is used to minimize the function for the large-scale problem (39919 unknowns, 2973 measurements). Three-dimensional (3-D) images of fluorescence absorption coefficients were reconstructed using the algorithm from experimental reflectance measurements under conditions of perfect and imperfect distribution of fluorophore within a single target. To our knowledge, this is the first time that targets have been reconstructed in three-dimensions from reflectance measurements with a clinically relevant phantom.

Many approaches to fluorescence tomography utilize some form of regularized nonlinear least-squares algorithm for data inversion, thus requiring repeated computation of the Jacobian sensitivity matrix relating changes in observable quantities, such as emission fluence, to changes in underlying optical parameters, such as fluorescence absorption. An exact adjoint formulation of these sensitivities comprises three terms, reflecting the individual contributions of 1) sensitivities of diffusion and decay coefficients at the emission wavelength, 2) sensitivities of diffusion and decay coefficients at the excitation wavelength, and 3) sensitivity of the emission source term. Simplifying linearity assumptions are computationally attractive in that they cause the first and second terms to drop out of the formulation. The relative importance of the three terms is thus explored in order to determine the extent to which these approximations introduce error. Computational experiments show that, while the third term of the sensitivity matrix has the largest magnitude, the second term becomes increasingly significant as target fluorophore concentration or volume increases. Image reconstructions from experimental data confirm that neglecting the second term results in overestimation of sensitivities and consequently overestimation of the value and volume of the fluorescent target, whereas contributions of the first term are so low that they are probably not worth the additional computational costs.

A three-dimensional image reconstruction for fluorescence-enhanced frequency-domain photon migration (FDPM) measurements in turbid media is developed and investigated for three different simulated measurement types: 1) absolute emission measurement, or emission measurements of phase and amplitude attenuation made for a given incident point source of excitation light; 2) referenced emission measurements made relative to an excitation measurement conducted at a single reference point away from the incident source; and 3) referenced emission measurements made relative to the excitation measurement conducted at identical points of detection. The image reconstruction algorithm employs a gradient-based constrained truncated Newton (CONTN) method which implements a bounding parameter, which can be used to govern the level of contrast used to discriminate tissue volumes from heterogeneous background tissues. Reverse differentiation technique is used to calculate the gradients. Using simulated data with superimposed noise to achieve a signal-to-noise ratio of 55 and 35 dB to mimic experimental excitation and emission FDPM measurements, respectively, we show the robustness of emission measurements referenced to excitation light. We investigate the performance of algorithm CONTN using these measurement techniques and show that the absorption coefficients due to fluorophore are reconstructed by CONTN accurately and efficiently. Furthermore, we demonstrate the performance of the bounding parameter for rejection of background artifacts owing to background tissue heterogeneity.

Most current efforts in near-infrared optical tomography are effectively limited to two-dimensional reconstructions due to the computationally intensive nature of full three-dimensional (3-D) data inversion. Previously, we described a new computationally efficient and statistically powerful inversion method APPRIZE (automatic progressive parameter-reducing inverse zonation and estimation). The APPRIZE method computes minimum-variance estimates of parameter values (here, spatially variant absorption due to a fluorescent contrast agent) and covariance, while simultaneously estimating the number of parameters needed as well as the size, shape, and location of the spatial regions that correspond to those parameters. Estimates of measurement and model error are explicitly incorporated into the procedure and implicitly regularize the inversion in a physically based manner. The optimal estimation of parameters is bounds-constrained, precluding infeasible values. In this paper, the APPRIZE method for optical imaging is extended for application to arbitrarily large 3-D domains through the use of domain decomposition. The effect of subdomain size on the performance of the method is examined by assessing the sensitivity for identifying 112 randomly located single-voxel heterogeneities in 58 3-D domains. Also investigated are the effects of unmodeled heterogeneity in background optical properties. The method is tested on simulated frequency-domain photon migration measurements at 100 MHz in order to recover absorption maps owing to fluorescent contrast agent. This study provides a new approach for computationally tractable 3-D optical tomography.