Registration of challenging pre-clinical brain images

The size and complexity of brain imaging studies in pre-clinical populations are increasing, and automated image analysis pipelines are urgently required. Pre-clinical populations can be subjected to controlled interventions (e.g., targeted lesions), which significantly change the appearance of the brain obtained by imaging. Existing systems for registration (the systematic alignment of scans into a consistent anatomical coordinate system), which assume image similarity to a reference scan, may fail when applied to these images. However, affine registration is a particularly vital pre-processing step for subsequent image analysis which is assumed to be an effective procedure in recent literature describing sophisticated techniques such as manifold learning. Therefore, in this paper, we present an affine registration solution that uses a graphical model of a population to decompose difficult pairwise registrations into a composition of steps using other members of the population. We developed this methodology in the context of a pre-clinical model of stroke in which large, variable hyper-intense lesions significantly impact registration performance. We tested this technique systematically in a simulated human population of brain tumour images before applying it to pre-clinical models of Parkinson's disease and stroke.

CT-PET weighted image fusion for separately scanned whole body rat

PURPOSE

The limited resolution and lack of spatial information in positron emission tomography (PET) images require the complementary anatomic information from the computed tomography (CT) and/or magnetic resonance imaging (MRI). Therefore, multimodality image fusion techniques such as PET/CT are critical in mapping the functional images to structural images and thus facilitate the interpretation of PET studies. In our experimental situation, the CT and PET images are acquired in separate scanners at different times and the inherent differences in the imaging protocols produce significant nonrigid changes between the two acquisitions in addition to dissimilar image characteristics. The registration conditions are also poor because CT images have artifacts due to the limitation of current scanning settings, while PET images are very blurry (in transmission-PET) and have vague anatomical structure boundaries (in emission-PET).

METHODS

The authors present a new method for whole body small animal multimodal registration. In particular, the authors register whole body rat CT image and PET images using a weighted demons algorithm. The authors use both the transmission-PET and the emission-PET images in the registration process emphasizing particular regions of the moving transmission-PET image using the emission-PET image. After a rigid transformation and a histogram matching between the CT and the transmission-PET images, the authors deformably register the transmission-PET image to the CT image with weights based on the intensity-normalized emission-PET image. For the deformable registration process, the authors develop a weighted demons registration method that can give preferences to particular regions of the input image using a weight image.

RESULTS

The authors validate the results with nine rat image sets using the M-Hausdorff distance (M-HD) similarity measure with different outlier-suppression parameters (OSP). In comparison with standard methods such as the regular demons and the normalized mutual information (NMI)-based nonrigid free-form deformation (FFD) registration, the proposed weighted demons registration method shows average M-HD errors: 3.99 ± 1.37 (OSP = 10), 5.04 ± 1.59 (OSP = 20) and 5.92 ± 1.61 (OSP = ∞) with statistical significance (p < 0.0003) respectively, while NMI-based nonrigid FFD has average M-HD errors: 5.74 ± 1.73 (OSP = 10), 7.40 ± 7.84 (OSP = 20) and 9.83 ± 4.13 (OSP = ∞), and the regular demons has average M-HD errors: 6.79 ± 0.83 (OSP = 10), 9.19 ± 2.39 (OSP = 20) and 11.63 ± 3.99 (OSP = ∞), respectively. In addition to M-HD comparisons, the visual comparisons on the faint-edged region between the CT and the aligned PET images also show the encouraging improvements over the other methods.

CONCLUSIONS

In the whole body multimodal registration between CT and PET images, the utilization of both the transmission-PET and the emission-PET images in the registration process by emphasizing particular regions of the transmission-PET image using an emission-PET image is effective. This method holds promise for other image fusion applications where multiple (more than two) input images should be registered into a single informative image.

Hybrid small animal imaging system combining magnetic resonance imaging with fluorescence tomography using single photon avalanche diode detectors

The high sensitivity of fluorescence imaging enables the detection of molecular processes in living organisms. However, diffuse light propagation in tissue prevents accurate recovery of tomographic information on fluorophore distribution for structures embedded deeper than 0.5 mm. Combining optical with magnetic resonance imaging (MRI) provides an accurate anatomical reference for fluorescence imaging data and thereby enables the correlation of molecular with high quality structural/functional information. We describe an integrated system for small animal imaging incorporating a noncontact fluorescence molecular tomography (FMT) system into an MRI detector. By adopting a free laser beam design geometrical constraints imposed by the use of optical fibers could be avoided allowing for flexible fluorescence excitation schemes. Photon detection based on a single-photon avalanche diode array enabled simultaneous FMT/MRI measurements without interference between modalities. In vitro characterization revealed good spatial accuracy of FMT data and accurate quantification of dye concentrations. Feasibility of FMT/MRI was demonstrated in vivo by simultaneous assessment of protease activity and tumor morphology in murine colon cancer xenografts.

Quantitative fluorescence tomography using a trimodality system: in vivo validation

A fully integrated trimodality fluorescence, diffuse optical, and x-ray computed tomography (FT/DOT/XCT) system for small animal imaging is reported in this work. The main purpose of this system is to obtain quantitatively accurate fluorescence concentration images using a multimodality approach. XCT offers anatomical information, while DOT provides the necessary background optical property map to improve FT image accuracy. The quantitative accuracy of this trimodality system is demonstrated in vivo. In particular, we show that a 2-mm-diam fluorescence inclusion located 8 mm deep in a nude mouse can only be localized when functional a priori information from DOT is available. However, the error in the recovered fluorophore concentration is nearly 87%. On the other hand, the fluorophore concentration can be accurately recovered within 2% error when both DOT functional and XCT structural a priori information are utilized together to guide and constrain the FT reconstruction algorithm.

Registration of planar bioluminescence to magnetic resonance and x-ray computed tomography images as a platform for the development of bioluminescence tomography reconstruction algorithms

The procedures we propose make possible the mapping of two-dimensional (2-D) bioluminescence image (BLI) data onto a skin surface derived from a three-dimensional (3-D) anatomical modality [magnetic resonance (MR) or computed tomography (CT)] dataset. This mapping allows anatomical information to be incorporated into bioluminescence tomography (BLT) reconstruction procedures and, when applied using sources visible to both optical and anatomical modalities, can be used to evaluate the accuracy of those reconstructions. Our procedures, based on immobilization of the animal and a priori determined fixed projective transforms, should be more robust and accurate than previously described efforts, which rely on a poorly constrained retrospectively determined warping of the 3-D anatomical information. Experiments conducted to measure the accuracy of the proposed registration procedure found it to have a mean error of 0.36+/-0.23 mm. Additional experiments highlight some of the confounds that are often overlooked in the BLT reconstruction process, and for two of these confounds, simple corrections are proposed.

Simultaneous in vivo dynamic magnetic resonance-diffuse optical tomography for small animal imaging

We present simultaneous measurement of enhancement kinetics of an optical and a magnetic resonance (MR) contrast agent in a small animal breast tumor model (R3230 ac) using a combined MR-diffuse optical tomographic (MR-DOT) imaging system. A mixture of a small molecular-weight MR contrast agent gadolinium-diethylene-triamine-pentaacetic acid (Gd-DTPA) and a large molecular-weight optical contrast agent indocyanine green (ICG) was administered intravenously for multimodal dynamic imaging. Coregistration of optical and MR images was accomplished using agar-water-based markers. Using T(2) and dynamic T(1) weighted MR images, we divided the entire tumor into two regions of interest (ROI): a viable and a nonviable region. The absorption enhancements in the ROIs were calculated. An enhancement of the ICG was observed in the viable region. On the contrary, there was a lower enhancement in the nonviable region.

Automatic nonrigid registration of whole body CT mice images

Three-dimensional intra- and intersubject registration of image volumes is important for tasks that include quantification of temporal/longitudinal changes, atlas-based segmentation, computing population averages, or voxel and tensor-based morphometry. While a number of methods have been proposed to address this problem, few have focused on the problem of registering whole body image volumes acquired either from humans or small animals. These image volumes typically contain a large number of articulated structures, which makes registration more difficult than the registration of head images, to which the majority of registration algorithms have been applied. This article presents a new method for the automatic registration of whole body computed tomography (CT) volumes, which consists of two main steps. Skeletons are first brought into approximate correspondence with a robust point-based method. Transformations so obtained are refined with an intensity-based nonrigid registration algorithm that includes spatial adaptation of the transformation's stiffness. The approach has been applied to whole body CT images of mice, to CT images of the human upper torso, and to human head and neck CT images. To validate the authors method on soft tissue structures, which are difficult to see in CT images, the authors use coregistered magnetic resonance images. They demonstrate that the approach they propose can successfully register image volumes even when these volumes are very different in size and shape or if they have been acquired with the subjects in different positions.

Rodent brain imaging with SPECT/CT

We evaluated methods of imaging rat models of stroke in vivo using a single photon emission computed tomography (SPECT) system dedicated to small animal imaging (X-SPECT, Gamma Medica-Ideas, Northridge, CA). An animal model of ischemic stroke was developed for in vivo SPECT/CT imaging using the middle cerebral artery occlusion (MCAO) technique. The presence of cerebral ischemia was verified in ex vivo studies using triphenyltetrazolium chloride (TTC) staining. In vivo radionuclide imaging of cerebral blood flow was performed in rats following MCAO using dynamic planar imaging of 99mTc-exametazime with parallel hole collimation. This was followed immediately by in vivo radionuclide imaging of cerebral blood flow with 99mTc-exametazime in the same animals using 1-mm pinhole SPECT. Correlated computed tomography imaging was performed to localize radiopharmaceutical uptake. The animals were allowed to recover and ex vivo autoradiography was performed with separate administration of 99mTc-exametazime. Time activity curve of 99mTc-exametazime showed that the radiopharmaceutical uptake could be maintained for over 9 min. The activity would be expected to be relatively stable for a much longer period, although the data were only obtained for 9 min. TTC staining revealed sizable infarcts by visual observation of inexistence of TTC stain in infracted tissues of MCAO rat brains. In vivo SPECT imaging showed cerebral blood flow deficit in the MCAO model, and the in vivo imaging result was confirmed with ex vivo autoradiography. We have demonstrated a capability of imaging regions of cerebral blood flow deficit in MCAO rat brains in vivo using a pinhole SPECT dedicated to small animal imaging.

Broad-spectrum multi-modality image registration: from PET, CT, and MRI to autoradiography, microscopy, and beyond

Image registration and fusion are increasingly important components of both clinical and small-animal imaging and have lead to the development of a variety of pertinent hardware and software tools, including multi-modality, e.g. PET-CT, devices. At the same time, advances in microscopic imaging, including phosphor-plate digital autoradiography and immunohistochemistry, now allow ultra-high (sub-100 microm)-resolution molecular characterization of tissue sections. To date, however, in vivo imaging of intact subjects and ex vivo imaging of harvested tissues sections have remained separate and distinct, making it difficult to reliably inter-compare the former and the latter. The Department of Medical Physics and the Radiation Biophysics Laboratory at Memorial Sloan-Kettering Cancer Center, under the direction of Dr. Clifton Ling, has now designed, fabricated, and tested a stereotactic imaging system for so-called "broad-spectrum" image registration, from coarser-resolution in vivo imaging modalities such as PET, CT, and MRI to ultra-high-resolution ex vivo imaging techniques such as histology, autoradiography, and immunohistochemistry.

Introduction to clinical and laboratory (small-animal) image registration and fusion

Imaging has long been a vital component of clinical medicine and, increasingly, of biomedical research in small-animals. Clinical and laboratory imaging modalities can be divided into two general categories, structural (or anatomical) and functional (or physiological). The latter, in particular, has spawned what has come to be known as "molecular imaging". Image registration and fusion have rapidly emerged as invaluable components of both clinical and small-animal imaging and has lead to the development and marketing of a variety of multi-modality, e.g. PET-CT, devices which provide registered and fused three-dimensional image sets. This paper briefly reviews the basics of image registration and fusion and available clinical and small-animal multi-modality instrumentation.

Differentiation between chondroma and squamous cell carcinoma: evaluation of subtraction image of MRI

The image processing procedure has become widely applied as a visual aid in imaging diagnosis. The subtraction image of MRI obtained by digitally subtracting an unenhanced image from a contrast-enhanced image, depicted a discrete distribution of the contrast agent. Subtraction images in the clinical cases were compared. The subtraction image in a case of chondroma demonstrated a relatively uniform distribution of a contrast agent, with a well delineated neoplastic lesion. The subtraction image in a case of squamous cell carcinoma, the contrast agent was distributed heterogeneously in the carcinoma and well permeated into the muscle. Different patterns of subtraction image were obtained in the clinical cases of chondroma and squamous cell carcinoma in this study. The findings suggest the potential usefulness of the subtraction image for diagnosing the degree of malignancy.

Regional perfusion of the equine carpus for antibiotic delivery

Regional perfusion of carpal tissues by forced intramedullary administration of fluids was evaluated in 10 horses. Results of subtraction radiography after perfusion with a contrast medium demonstrated that perfusate was delivered to the carpal tissues by the venous system. Perfused India ink was distributed uniformly in the antebrachiocarpal and middle carpal synovial membranes. Histologically, the ink was within the venules of the synovial villi. Immediately after perfusion with gentamicin sulfate (1 g), the gentamicin concentrations in the synovial fluid and synovial membrane of the antebrachiocarpal joint were 349 +/- 240 micrograms/mL and 358 +/- 264 micrograms/g, respectively. When gentamicin concentrations in the synovial fluid of the antebrachiocarpal joint and serum were measured 0, 0.5, 1, 4, 8, 12, and 24 hours after carpal perfusion, the mean peak gentamicin concentration in the synovial fluid was 589 +/- 429 micrograms/mL. At hour 24, the mean gentamicin concentration in the synovial fluid was 4.8 +/- 2.0 micrograms/mL. The resulting peak gentamicin concentration in the serum was 23.7 +/- 14.5 micrograms/mL immediately after the perfusion; it decreased below the desired trough level of 1 micrograms/mL between hours 4 and 8.