A novel approach is presented for obtaining fast robust three-dimensional (3-D) reconstructions of bioluminescent reporters buried deep inside animal subjects from multispectral images of surface bioluminescent photon densities. The proposed method iteratively acts upon the equations relating the multispectral data to the luminescent distribution with high computational efficiency to provide robust 3-D reconstructions. Unlike existing algebraic reconstruction techniques, the proposed method is designed to use adaptive projections that iteratively guide the updates to the solution with improved speed and robustness. Contrary to least-squares reconstruction methods, the proposed technique does not require parameter selection or optimization for optimal performance. Additionally, optimized schemes for thresholding, sampling, and ordering of the bioluminescence tomographic data used by the proposed method are presented. The performance of the proposed approach in reconstructing the shape, volume, flux, and depth of luminescent inclusions is evaluated in a multitude of phantom-based and dual-modality in vivo studies in which calibrated sources are implanted in animal subjects and imaged in a dual-modality optical/computed tomography platform. Statistical analysis of the errors in the depth and flux of the reconstructed inclusions and the convergence time of the proposed method is used to demonstrate its unbiased performance, low error variance, and computational efficiency.
A new method is described for obtaining a 3-D reconstruction of a bioluminescent light source distribution inside a living animal subject, from multispectral images of the surface light emission acquired on charge-coupled device (CCD) camera. The method uses the 3-D surface topography of the animal, which is obtained from a structured light illumination technique. The forward model of photon transport is based on the diffusion approximation in homogeneous tissue with a local planar boundary approximation for each mesh element, allowing rapid calculation of the forward Green's function kernel. Absorption and scattering properties of tissue are measured a priori as input to the algorithm. By using multispectral images, 3-D reconstructions of luminescent sources can be derived from images acquired from only a single view. As a demonstration, the reconstruction technique is applied to determine the location and brightness of a source embedded in a homogeneous phantom subject in the shape of a mouse. The technique is then evaluated with real mouse models in which calibrated sources are implanted at known locations within living tissue. Finally, reconstructions are demonstrated in a PC3M-luc (prostate tumor line) metastatic tumor model in nude mice.
Physiologic systems can be represented by compartmental models which describe the uptake of radio-labeled tracers from blood to tissue and their subsequent washout. Arterial and venous time-activity curves from isolated heart experiments are analyzed using spectral analysis, in which the impulse response function is represented by a sum of decaying exponentials. Resolution and uniqueness tests are conducted by synthesizing isolated heart data with predefined compartmental models, adding noise, and applying the spectral analysis technique. Venous time-activity curves are generated by convolving a typical arterial input function with the predefined spectrum. The coefficients of a set of decaying exponential basis functions are determined using a non- negative least squares algorithm, and results are compared with the predefined spectrum. The uniqueness of spectral method solutions is investigated by computing model covariance matrices, using error propagation and prior knowledge of noise distributions. Coupling between model parameters is illustrated with correlation matrices.