The three-year UK NeuroGrid project aims to develop a Grid-based collaborative research environment to support the data and compute needs for a UK Neurosciences community. This paper describes the challenges in developing this architecture and details initial results from the development of its first prototype to support psychosis, dementia and stroke research and the social challenges of such a collaborative research project. The paper discusses approaches being taken to explore the collaborative science process to inform the requirements for follow on prototypes and methods utilized to develop an effective project team.

The main goal of this work is the generation of ground-truth data for the validation of atrophy measurement techniques, commonly used in the study of neurodegenerative diseases such as dementia. Several techniques have been used to measure atrophy in cross-sectional and longitudinal studies, but it is extremely difficult to compare their performance since they have been applied to different patient populations. Furthermore, assessment of performance based on phantom measurements or simple scaled images overestimates these techniques' ability to capture the complexity of neurodegeneration of the human brain. We propose a method for atrophy simulation in structural magnetic resonance (MR) images based on finite-element methods. The method produces cohorts of brain images with known change that is physically and clinically plausible, providing data for objective evaluation of atrophy measurement techniques. Atrophy is simulated in different tissue compartments or in different neuroanatomical structures with a phenomenological model. This model of diffuse global and regional atrophy is based on volumetric measurements such as the brain or the hippocampus, from patients with known disease and guided by clinical knowledge of the relative pathological involvement of regions and tissues. The consequent biomechanical readjustment of structures is modelled using conventional physics-based techniques based on biomechanical tissue properties and simulating plausible tissue deformations with finite-element methods. A thermoelastic model of tissue deformation is employed, controlling the rate of progression of atrophy by means of a set of thermal coefficients, each one corresponding to a different type of tissue. Tissue characterization is performed by means of the meshing of a labelled brain atlas, creating a reference volumetric mesh that will be introduced to a finite-element solver to create the simulated deformations. Preliminary work on the simulation of acquisition artefacts is also presented. Cross-sectional and longitudinal sets of simulated data are shown and a visual classification protocol has been used by experts to rate real and simulated scans according to their degree of atrophy. Results confirm the potential of the proposed methodology.

Measures of overlap of labelled regions of images, such as the Dice and Tanimoto coefficients, have been extensively used to evaluate image registration and segmentation algorithms. Modern studies can include multiple labels defined on multiple images yet most evaluation schemes report one overlap per labelled region, simply averaged over multiple images. In this paper, common overlap measures are generalized to measure the total overlap of ensembles of labels defined on multiple test images and account for fractional labels using fuzzy set theory. This framework allows a single "figure-of-merit" to be reported which summarises the results of a complex experiment by image pair, by label or overall. A complementary measure of error, the overlap distance, is defined which captures the spatial extent of the nonoverlapping part and is related to the Hausdorff distance computed on grey level images. The generalized overlap measures are validated on synthetic images for which the overlap can be computed analytically and used as similarity measures in nonrigid registration of three-dimensional magnetic resonance imaging (MRI) brain images. Finally, a pragmatic segmentation ground truth is constructed by registering a magnetic resonance atlas brain to 20 individual scans, and used with the overlap measures to evaluate publicly available brain segmentation algorithms.

Recent innovations in drug therapies have made it highly desirable to obtain sensitive biomarkers of disease progression that can be used to quantify the performance of candidate disease modifying drugs. In order to measure potential image-based biomarkers of disease progression in an experimental model of rheumatoid arthritis (RA), we present two different methods to automatically quantify changes in a bone in in-vivo serial magnetic resonance (MR) images from the model. Both methods are based on rigid and nonrigid image registration to perform the analysis. The first method uses segmentation propagation to delineate a bone from the serial MR images giving a global measure of temporal changes in bone volume. The second method uses rigid body registration to determine intensity change within a bone, and then maps these into a reference coordinate system using nonrigid registration. This gives a local measure of temporal changes in bone lesion volume. We detected significant temporal changes in local bone lesion volume in five out of eight identified candidate bone lesion regions, and significant difference in local bone lesion volume between male and female subjects in three out of eight candidate bone lesion regions. But the global bone volume was found to be fluctuating over time. Finally, we compare our findings with histology of the subjects and the manual segmentation of bone lesions.

Large-scale clinical studies in neuro-imaging are hampered by several factors including variances in acquisition techniques, quality assurance and access to remote datasets. The Neurogrid project will build on the experience of other UK e-science projects to assemble a grid infrastructure, and apply this to three exemplar areas: stroke, dementia and psychosis, to conduct collaborative neuroscience research.

The hybrid magnetic resonance (MR)/X-ray suite (XMR) is a recently introduced imaging solution that provides new possibilities for guidance of cardiovascular catheterization procedures. We have previously described and validated a technique based on optical tracking to register MR and X-ray images obtained from the sliding table XMR configuration. The aim of our recent work was to extend our technique by providing an improved calibration stage, real-time guidance during cardiovascular catheterization procedures, and further off-line analysis for mapping cardiac electrical data to patient anatomy. Specially designed optical trackers and a dedicated calibration object have resulted in a single calibration step that can be efficiently checked and updated before each procedure. An X-ray distortion model has been implemented that allows for distortion correction for arbitrary c-arm orientations. During procedures, the guidance system provides a real-time combined MR/X-ray image display consisting of live X-ray images with registered recently acquired MR derived anatomy. It is also possible to reconstruct the location of catheters seen during X-ray imaging in the MR derived patient anatomy. We have applied our registration technique to 13 cardiovascular catheterization procedures. Our system has been used for the real-time guidance of ten radiofrequency ablations and one aortic stent implantation. We demonstrate the real-time guidance using two exemplar cases. In a further two cases we show how off-line analysis of registered image data, acquired during electrophysiology study procedures, has been used to map cardiac electrical measurements to patient anatomy for two different types of mapping catheters. The cardiologists that have used the guidance system suggest that real-time XMR guidance could have substantial value in difficult interventional and electrophysiological procedures, potentially reducing procedure time and delivered radiation dose. Also, the ability to map measured electrical data to patient specific anatomy provides improved visualization and a path to investigation of cardiac electromechanical models.

Recent studies have shown that the surface of the brain is deformed by up to 20 mm after the skull is opened during neurosurgery, which could lead to substantial error in commercial image-guided surgery systems. We quantitatively analyze the intraoperative brain deformation of 24 subjects to investigate whether simple rules can describe or predict the deformation. Interventional magnetic resonance images acquired at the start and end of the procedure are registered nonrigidly to obtain deformation values throughout the brain. Deformation patterns are investigated quantitatively with respect to the location and magnitude of deformation, and to the distribution and principal direction of the displacements. We also measure the volume change of the lateral ventricles by manual segmentation. Our study indicates that brain shift occurs predominantly in the hemisphere ipsi-lateral to the craniotomy, and that there is more brain deformation during resection procedures than during biopsy or functional procedures. However, the brain deformation patterns are extremely complex in this group of subjects. This paper quantitatively demonstrates that brain deformation occurs not only at the surface, but also in deeper brain structure, and that the principal direction of displacement does not always correspond with the direction of gravity. Therefore, simple computational algorithms that utilize limited intraoperative information (e.g., brain surface shift) will not always accurately predict brain deformation at the lesion.

TEST 02 - Elsevier's Scopus, the largest abstract and citation database of peer-reviewed literature. Search and access research from the science, technology, medicine, social sciences and arts and humanities fields.

We describe a registration and tracking technique to integrate cardiac X-ray images and cardiac magnetic resonance (MR) images acquired from a combined X-ray and MR interventional suite (XMR). Optical tracking is used to determine the transformation matrices relating MR image coordinates and X-ray image coordinates. Calibration of X-ray projection geometry and tracking of the X-ray C-arm and table enable three-dimensional (3-D) reconstruction of vessel centerlines and catheters from bi-plane X-ray views. We can, therefore, combine single X-ray projection images with registered projection MR images from a volume acquisition, and we can also display 3-D reconstructions of catheters within a 3-D or multi-slice MR volume. Registration errors were assessed using phantom experiments. Errors in the combined projection images (two-dimensional target registration error--TRE) were found to be 2.4 to 4.2 mm, and the errors in the integrated volume representation (3-D TRE) were found to be 4.6 to 5.1 mm. These errors are clinically acceptable for alignment of images of the great vessels and the chambers of the heart. Results are shown for two patients. The first involves overlay of a catheter used for invasive pressure measurements on an MR volume that provides anatomical context. The second involves overlay of invasive electrode catheters (including a basket catheter) on a tagged MR volume in order to relate electrophysiology to myocardial motion in a patient with an arrhythmia. Visual assessment of these results suggests the errors were of a similar magnitude to those obtained in the phantom measurements.

This paper presents a novel method for validation of nonrigid medical image registration. This method is based on the simulation of physically plausible, biomechanical tissue deformations using finite-element methods. Applying a range of displacements to finite-element models of different patient anatomies generates model solutions which simulate gold standard deformations. From these solutions, deformed images are generated with a range of deformations typical of those likely to occur in vivo. The registration accuracy with respect to the finite-element simulations is quantified by co-registering the deformed images with the original images and comparing the recovered voxel displacements with the biomechanically simulated ones. The functionality of the validation method is demonstrated for a previously described nonrigid image registration technique based on free-form deformations using B-splines and normalized mutual information as a voxel similarity measure, with an application to contrast-enhanced magnetic resonance mammography image pairs. The exemplar nonrigid registration technique is shown to be of subvoxel accuracy on average for this particular application. The validation method presented here is an important step toward more generic simulations of biomechanically plausible tissue deformations and quantification of tissue motion recovery using nonrigid image registration. It will provide a basis for improving and comparing different nonrigid registration techniques for a diversity of medical applications, such as intrasubject tissue deformation or motion correction in the brain, liver or heart.

Previous work has suggested the existence of differences between the cerebral cortex of normal individuals, and those of patients with diseases such as epilepsy and schizophrenia. These shape abnormalities may be of developmental origin. Improved shape measures could provide useful tools for neuroscience research and patient diagnosis. We consider the theoretically desirable properties of measures of brain shape. We have implemented seven measures, three from the neuroscience literature, and four new to this field. Three of the measures are zero-order and four are second-order with respect to the surface. We validate the measures using simple geometrical shapes, and a collection of magnetic resonance scans of ten histologically normal ex vivo fetal brains with gestational ages from 19-42 weeks. We then apply the measures to MR scans from two histologically abnormal ex vivo brains. We demonstrate that our implementation of the measures is sensitive to anatomical variability rather than to the discreteness of the image data. All the measures were sensitive to changes in shape during fetal development. Several of the measures could distinguish between the normal and abnormal fetal brains. We propose a multivariate approach to studying the shape of the cerebral cortex, in which both zero-order and second-order measures are used to quantify folding.

Nonrigid registration can automatically quantify small changes in volume of anatomical structures over time by means of segmentation propagation. Here, we use a nonrigid registration algorithm based on optimising normalized mutual information to quantify small changes in brain ventricle volume in magnetic resonance (MR) images of a group of five patients treated with growth hormone replacement therapy and a control group of six volunteers. The lateral ventricles are segmented from each subject image by registering with the brainweb image which has this structure delineated. The mean (standard deviation) volume change measurements are 1.09 (0.73) cm3 for the patient group and 0.08 (0.62) cm3 for the volunteer group; this difference is statistically significant at the 1% level. We validate our volume measurements by determining the precision from three consecutive scans of five volunteers and also comparing the measurements to previously published volume change estimates obtained by visual inspection of difference images. Results demonstrate a precision of sigma < or = 0.52 cm3 (n = 5) and a rank correlation coefficient with assessed difference images of p = 0.7 (n = 11). To determine the level of shape correspondence we manually segmented subject's ventricles and compared them to the propagations using a voxel overlap similarity index, this gave a mean similarity index of 0.81 (n = 7).

This paper describes a quantitative assessment of respiratory motion of the heart and the construction of a model of respiratory motion. Three-dimensional magnetic resonance scans were acquired on eight normal volunteers and ten patients. The volunteers were imaged at multiple positions in the breathing cycle between full exhalation and full inhalation while holding their breath. The exhalation volume was segmented and used as a template to which the other volumes were registered using an intensity-based rigid registration algorithm followed by nonrigid registration. The patients were imaged at inhale and exhale only. The registration results were validated by visual assessment and consistency measurements indicating subvoxel registration accuracy. For all subjects, we assessed the nonrigid motion of the heart at the right coronary artery, right atrium, and left ventricle. We show that the rigid-body motion of the heart is primarily in the craniocaudal direction with smaller displacements in the right-left and anterior-posterior directions; this is in agreement with previous studies. Deformation was greatest for the free wall of the right atrium and the left ventricle; typical deformations were 3-4 mm with deformations of up to 7 mm observed in some subjects. Using the registration results, landmarks on the template surface were mapped to their correct positions through the breathing cycle. Principal component analysis produced a statistical model of the motion and deformation of the heart. We discuss how this model could be used to assist motion correction.

The authors propose a novel method to register two or more optical images to a 3D surface model. The potential applications of such a registration method could be in medicine for example, in image guided interventions, surveillance and identification, industrial inspection, or telemanipulation in remote or hostile environments. Registration is performed by optimizing a similarity measure with respect to the transformation parameters. We propose a novel similarity measure based on "photo-consistency." For each surface point, the similarity measure computes how consistent the corresponding optical image information in each view is with a lighting model. The relative pose of the optical images must be known. We validate the system using data from an optical-based surface reconstruction system and surfaces derived from magnetic resonance (MR) images of the human face. We test the accuracy and robustness of the system with respect to the number of video images, video image noise, errors in surface location and area, and complexity of the matched surfaces. We demonstrate the algorithm working on 10 further optical-based reconstructions of the human head and skin surfaces derived from MR images of the heads of five volunteers. Matching four optical images to a surface model produced a 3D error of between 1.45 and 1.59 mm, at a success rate of 100 percent, where the initial misregistration was up to 16 mm or degrees from the registration position.

The problem of providing surgical navigation using image overlays on the operative scene can be split into four main tasks--calibration of the optical system; registration of preoperative images to the patient; system and patient tracking, and display using a suitable visualization scheme. To achieve a convincing result in the magnified microscope view a very high alignment accuracy is required. We have simulated an entire image overlay system to establish the most significant sources of error and improved each of the stages involved. The microscope calibration process has been automated. We have introduced bone-implanted markers for registration and incorporated a locking acrylic dental stent (LADS) for patient tracking. The LADS can also provide a less-invasive registration device with mean target error of 0.7 mm in volunteer experiments. These improvements have significantly increased the alignment accuracy of our overlays. Phantom accuracy is 0.3-0.5 mm and clinical overlay errors were 0.5-1.0 mm on the bone fiducials and 0.5-4 mm on target structures. We have improved the graphical representation of the stereo overlays. The resulting system provides three-dimensional surgical navigation for microscope-assisted guided interventions (MAGI).

We have evaluated eight different similarity measures used for rigid body registration of serial magnetic resonance (MR) brain scans. To assess their accuracy we used 33 clinical three-dimensional (3-D) serial MR images, with deformable extradural tissue excluded by manual segmentation and simulated 3-D MR images with added intensity distortion. For each measure we determined the consistency of registration transformations for both sets of segmented and unsegmented data. We have shown that of the eight measures tested, the ones based on joint entropy produced the best consistency. In particular, these measures seemed to be least sensitive to the presence of extradural tissue. For these data the difference in accuracy of these joint entropy measures, with or without brain segmentation, was within the threshold of visually detectable change in the difference images.

In this paper we present a new approach for the nonrigid registration of contrast-enhanced breast MRI. A hierarchical transformation model of the motion of the breast has been developed. The global motion of the breast is modeled by an affine transformation while the local breast motion is described by a free-form deformation (FFD) based on B-splines. Normalized mutual information is used as a voxel-based similarity measure which is insensitive to intensity changes as a result of the contrast enhancement. Registration is achieved by minimizing a cost function, which represents a combination of the cost associated with the smoothness of the transformation and the cost associated with the image similarity. The algorithm has been applied to the fully automated registration of three-dimensional (3-D) breast MRI in volunteers and patients. In particular, we have compared the results of the proposed nonrigid registration algorithm to those obtained using rigid and affine registration techniques. The results clearly indicate that the nonrigid registration algorithm is much better able to recover the motion and deformation of the breast than rigid or affine registration algorithms.

In recent years, advances in computer technology and a significant increase in the accuracy of medical imaging have made it possible to develop systems that can assist the clinician in diagnosis, planning, and treatment. This paper deals with an area that is generally referred to as computer-assisted surgery, image-directed surgery, or image-guided surgery. We report the research, development, and clinical validation performed since January 1996 in the European Applications in Surgical Interventions (EASI) project, which is funded by the European Commission in their "4th Framework Telematics Applications for Health" program. The goal of this project is the improvement of the effectiveness and quality of image-guided neurosurgery of the brain and image-guided vascular surgery of abdominal aortic aneurysms, while at the same time reducing patient risks and overall cost. We have developed advanced prototype systems for preoperative surgical planning and intraoperative surgical navigation, and we have extensively clinically validated these systems. The prototype systems and the clinical validation results are described in this paper.

In a previous study we demonstrated that automatic retrospective registration algorithms can frequently register magnetic resonance (MR) and computed tomography (CT) images of the brain with an accuracy of better than 2 mm, but in that same study we found that such algorithms sometimes fail, leading to errors of 6 mm or more. Before these algorithms can be used routinely in the clinic, methods must be provided for distinguishing between registration solutions that are clinically satisfactory and those that are not. One approach is to rely on a human observer to inspect the registration results and reject images that have been registered with insufficient accuracy. In this paper, we present a methodology for evaluating the efficacy of the visual assessment of registration accuracy. Since the clinical requirements for level of registration accuracy are likely to be application dependent, we have evaluated the accuracy of the observer's estimate relative to six thresholds: 1-6 mm. The performance of the observers was evaluated relative to the registration solution obtained using external fiducial markers that are screwed into the patient's skull and that are visible in both MR and CT images. This fiducial marker system provides the gold standard for our study. Its accuracy is shown to be approximately 0.5 mm. Two experienced, blinded observers viewed five pairs of clinical MR and CT brain images, each of which had each been misregistered with respect to the gold standard solution. Fourteen misregistrations were assessed for each image pair with misregistration errors distributed between 0 and 10 mm with approximate uniformity. For each misregistered image pair each observer estimated the registration error (in millimeters) at each of five locations distributed around the head using each of three assessment methods. These estimated errors were compared with the errors as measured by the gold standard to determine agreement relative to each of the six thresholds, where agreement means that the two errors lie on the same side of the threshold. The effect of error in the gold standard itself is taken into account in the analysis of the assessment methods. The results were analyzed by means of the Kappa statistic, the agreement rate, and the area of receiver-operating-characteristic (ROC) curves. No assessment performed well at 1 mm, but all methods performed well at 2 mm and higher. For these five thresholds, two methods agreed with the standard at least 80% of the time and exhibited mean ROC areas greater than 0.84. One of these same methods exhibited Kappa statistics that indicated good agreement relative to chance (Kappa > 0.6) between the pooled observers and the standard for these same five thresholds. Further analysis demonstrates that the results depend strongly on the choice of the distribution of misregistration errors presented to the observers.

All image-guided neurosurgical systems that we are aware of assume that the head and its contents behave as a rigid body. It is important to measure intraoperative brain deformation (brain shift) to provide some indication of the application accuracy of image-guided surgical systems, and also to provide data to develop and validate nonrigid registration algorithms to correct for such deformation. We are collecting data from patients undergoing neurosurgery in a high-field (1.5 T) interventional magnetic resonance (MR) scanner. High-contrast and high-resolution gradient-echo MR image volumes are collected immediately prior to surgery, during surgery, and at the end of surgery, with the patient intubated and lying on the operating table in the operative position. In this paper we report initial results from six patients: one freehand biopsy, one stereotactic functional procedure, and four resections. We investigate intraoperative brain deformation by examining threshold boundary overlays and difference images and by measuring ventricular volume. We also present preliminary results obtained using a nonrigid registration algorithm to quantify deformation. We found that some cases had much greater deformation than others, and also that, regardless of the procedure, there was very little deformation of the midline, the tentorium, the hemisphere contralateral to the procedure, and ipsilateral structures except those that are within 1 cm of the lesion or are gravitationally above the surgical site.

A comparison of six similarity measures for use in intensity-based two-dimensional-three-dimensional (2-D-3-D) image registration is presented. The accuracy of the similarity measures are compared to a "gold-standard" registration which has been accurately calculated using fiducial markers. The similarity measures are used to register a computed tomography (CT) scan of a spine phantom to a fluoroscopy image of the phantom. The registration is carried out within a region-of-interest in the fluoroscopy image which is user defined to contain a single vertebra. Many of the problems involved in this type of registration are caused by features which were not modeled by a phantom image alone. More realistic "gold-standard" data sets were simulated using the phantom image with clinical image features overlaid. Results show that the introduction of soft-tissue structures and interventional instruments into the phantom image can have a large effect on the performance of some similarity measures previously applied to 2-D-3-D image registration. Two measures were able to register accurately and robustly even when soft-tissue structures and interventional instruments were present as differences between the images. These measures were pattern intensity and gradient difference. Their registration accuracy, for all the rigid-body parameters except for the source to film translation, was within a root-mean-square (rms) error of 0.54 mm or degrees to the "gold-standard" values. No failures occurred while registering using these measures.

Registration of intraoperative fluoroscopy images with preoperative three-dimensional (3-D) CT images can be used for several purposes in image-guided surgery. On the one hand, it can be used to display the position of surgical instruments, which are being tracked by a localizer, in the preoperative CT scan. On the other hand, the registration result can be used to project preoperative planning information or important anatomical structures visible in the CT image onto the fluoroscopy image. For this registration task, a novel voxel-based method in combination with a new similarity measure (pattern intensity) has been developed. The basic concept of the method is explained at the example of two-dimensional (2-D)/3-D registration of a vertebra in an X-ray fluoroscopy image with a 3-D CT image. The registration method is described, and the results for a spine phantom are presented and discussed. Registration has been carried out repeatedly with different starting estimates to study the capture range. Information about registration accuracy has been obtained by comparing the registration results with a highly accurate "ground-truth" registration, which has been derived from fiducial markers attached to the phantom prior to imaging. In addition, registration results for different vertebrae have been compared. The results show that the rotation parameters and the shifts parallel to the projection plane can accurately be determined from a single projection. Because of the projection geometry, the accuracy of the height above the projection plane is significantly lower.

We present the use of an entropy focus criterion to enable automatic focusing of motion corrupted magnetic resonance images. We demonstrate the principle using illustrative examples from cooperative volunteers. Our technique can determine unknown patient motion or use knowledge of motion from other measures as a starting estimate. The motion estimate is used to compensate the acquired data and is iteratively refined using the image entropy. Entropy focuses the whole image principally by favoring the removal of motion induced ghosts and blurring from otherwise dark regions of the image. Using only the image data, and no special hardware or pulse sequences, we demonstrate correction for arbitrary rigid-body translational motion in the imaging plane and for a single rotation. Extension to three-dimensional (3-D) and more general motion should be possible. The algorithm is able to determine volunteer motion well. The mean absolute deviation between algorithm and navigator-echo-determined motion is comparable to the displacement step size used in the algorithm. Local deviations from the recorded motion or navigator-determined motion are explained and we indicate how enhanced focus criteria may be derived. In all cases we were able to compensate images for patient motion, reducing blurring and ghosting.

Multiple sources of 3D medical image data can be used to construct detailed patient representations. Typically registration is achieved assuming the validity of rigid body transformation. In many applications, and in particular when updating representations used for guidance during surgery and therapeutic interventions, this assumption is inappropriate. In this paper we describe a general method for 3D deformation, show how registration can incorporate a composite of rigid body and deformation components and illustrate this methodology on 3 example sets of images. The first is a repeated 3D MR scan of the abdomen of a volunteer who purposely changed position between scans; the second is an MR and CT scan of the head and neck, in which the patient was in a different position for the two scans; and the third is a set of MR and CT images of the head taken before and after epilepsy surgery. Non rigid deformation and composite warping showed significant improvement in registration accuracy in each case. 1

This paper presents a fully automated technique for the rigid body registration of clinically acquired MR and CT images of the head. We describe our multiresolution approach to the optimisation of voxel similarity measures and evaluate the performance of a number of measures when presented with clinical images with a small overlapping volume. Results are compared with those derived from manual registration by identification of corresponding point landmarks. We show that by limiting the measures to intra-cranial regions of the images, not containing deformable skin surface features, a greater accuracy may be provided for certain types of truncated image.

In this paper we present our work on using intensity feature spaces to study the relationship between voxel values in registered and unreg-istered medical images. By taking a simple image model we predict structures that we might expect to find in an intensity feature space produced from different modality images of the same scene. We show how this structure will be modified by image noise, misregistration and differing point spread functions of the two modalities. We show examples of such structure in feature spaces created from clinically ac-quired Magnetic Resonance (MR) and Positron Emission Tomography (PET) image data. We show how two simple measures of voxel sim-ilarity based on these feature space observations, a modified variance of intensity ratio and the 3rd order moment of the feature space his-togram, can be used to quantify image misregistration. The 3rd order moment measure is then used with a genetic optimisation algorithm to automatically register pre and post Gadolinium injection MR images. 1

Abstract The combination of multimodality medical images is generating increasing clinical interest, but the absence of robust automatic registration algorithms is an obstacle to these techniques entering routine clinical use. In this paper, we present a new algorithm based on surface fitting that makes use of anatomical knowledge of adjacency of identified anatomical structures to solve the 3D rigid body registration problem. The algorithm has been applied to the registration of MR CT and vascular images of the head, and has achieved accuracies similar to those obtained using interactive location of corresponding landmarks.