Registration of head volume images using implantable fiducial markers

MR-guided procedures using contemporaneous imaging frameless stereotaxis in an open-configuration system

Frameless MR-guided procedures have had limited application using conventional closed magnets, due largely to the technical difficulties involved. As a result of in-room MR image-monitoring capabilities, new open-design magnets now allow frameless stereotaxis using contemporaneous imaging to guide more invasive procedures. We evaluate our clinical experience with this new technique. An open-design 0.2 T magnet (Siemens OPEN) combined with an in-room monitor was used for 33 frameless MR-guided procedures (aspiration cytology, biopsy, and/or treatment) in a variety of locations in the head, neck, spine, brain, pelvis, and abdomen. Success of the procedure was based on the ability to accurately position the instrument in the target region to allow biopsy and/or treatment. The open-design magnet allowed the physician to directly access the patient for frameless stereotaxis as the procedure was performed. The in-room monitor provided contemporaneous imaging feedback during the procedure for successful placement of the instrument in the target region. Twenty-eight biopsy and five treatment procedures were performed. In all cases the technique resulted in successful placement of the instrument within the target tissue to complete the procedure. MR-guided procedures using contemporaneous imaging frameless stereotaxis are possible in an open-design magnet with in-room image monitoring and offer exciting possibilities for further development.

Predicting error in rigid-body point-based registration

Guidance systems designed for neurosurgery, hip surgery, and spine surgery, and for approaches to other anatomy that is relatively rigid can use rigid-body transformations to accomplish image registration. These systems often rely on point-based registration to determine the transformation, and many such systems use attached fiducial markers to establish accurate fiducial points for the registration, the points being established by some fiducial localization process. Accuracy is important to these systems, as is knowledge of the level of that accuracy. An advantage of marker-based systems, particularly those in which the markers are bone-implanted, is that registration error depends only on the fiducial localization error (FLE) and is thus to a large extent independent of the particular object being registered. Thus, it should be possible to predict the clinical accuracy of marker-based systems on the basis of experimental measurements made with phantoms or previous patients. This paper presents two new expressions for estimating registration accuracy of such systems and points out a danger in using a traditional measure of registration accuracy. The new expressions represent fundamental theoretical results with regard to the relationship between localization error and registration error in rigid-body, point-based registration. Rigid-body, point-based registration is achieved by finding the rigid transformation that minimizes "fiducial registration error" (FRE), which is the root mean square distance between homologous fiducials after registration. Closed form solutions have been known since 1966. The expected value (FRE2) depends on the number N of fiducials and expected squared value of FLE, (FLE-2, but in 1979 it was shown that (FRE2) is approximately independent of the fiducial configuration C. The importance of this surprising result seems not yet to have been appreciated by the registration community: Poor registrations caused by poor fiducial configurations may appear to be good due to a small FRE value. A more critical and direct measure of registration error is the "target registration error" (TRE), which is the distance between homologous points other than the centroids of fiducials. Efforts to characterize its behavior have been made since 1989. Published numerical simulations have shown that (TRE2) is roughly proportional to (FLE2)/N and, unlike (FRE2), does depend in some way on C. Thus, FRE, which is often used as feedback to the surgeon using a point-based guidance system, is in fact an unreliable indicator of registration-accuracy. In this work we derive approximate expressions for (TRE2), and for the expected squared alignment error of an individual fiducial. We validate both approximations through numerical simulations. The former expression can be used to provide reliable feedback to the surgeon during surgery and to guide the placement of markers before surgery, or at least to warn the surgeon of potentially dangerous fiducial placements; the latter expression leads to a surprising conclusion: Expected registration accuracy (TRE) is worst near the fiducials that are most closely aligned! This revelation should be of particular concern to surgeons who may at present be relying on fiducial alignment as an indicator of the accuracy of their point-based guidance systems.

Surface-based registration of CT images to physical space for image-guided surgery of the spine: a sensitivity study

This paper presents a method designed to register preoperative computed tomography (CT) images to vertebral surface points acquired intraoperatively from ultrasound (US) images or via a tracked probe. It also presents a comparison of the registration accuracy achievable with surface points acquired from the entire posterior surface of the vertebra to the accuracy achievable with points acquired only from the spinous process and central laminar regions. Using a marker-based method as a reference, this work shows that submillimetric registration accuracy can be obtained even when a small portion of the posterior vertebral surface is used for registration. It also shows that when selected surface patches are used, CT slice thickness is not a critical parameter in the registration process. Furthermore, the paper includes qualitative results of registering vertebral surface points in US images to multiple CT slices. The method has been tested with US points and physical points on a plastic spine phantom and with simulated data on a patient CT scan.

Calibration of tracking systems in a surgical environment

The purpose of this paper was to assess to what extent an optical tracking system (OTS) used for position determination in computer-aided surgery (CAS) can be enhanced by combining it with a direct current (dc) driven electromagnetic tracking system (EMTS). The main advantage of the EMTS is the fact that it is not dependent on a free line-of-sight. Unfortunately, the accuracy of the EMTS is highly affected by nearby ferromagnetic materials. We have explored to what extent the influence of the metallic equipment in the operating room (OR) can be compensated by collecting precise information on the nonlinear local error in the EMTS by using the OTS for setting up a calibration look-up table. After calibration of the EMTS and registration of the sensor systems in the OR we have found the average euclidean deviation in position readings between the dc tracker and the OTS reduced from 2.9+/-1.0 mm to 2.1+/-0.8 mm within a half-sphere of 530-mm radius around the magnetic field emitter. Furthermore we have found the calibration to be stable after re-registration of the sensors under varying conditions such as different heights of the OR table and varying positions of the OR equipment over a longer time interval. These results encourage the further development of a hybrid magnetooptical tracker for computer-aided surgery where the electromagnetic tracker acts as an auxiliary source of position information for the optical system. Strategies for enhancing the reliability of the proposed hybrid magnetooptic tracker by detecting artifacts induced by mobile ferromagnetic objects such as surgical tools are discussed.

Registration of head CT images to physical space using a weighted combination of points and surfaces

Most previously reported registration techniques that align three-dimensional image volumes by matching geometrical features such as points or surfaces use a single type of feature. We recently reported a hybrid registration technique that uses a weighted combination of multiple geometrical feature shapes. In this study we use the weighted geometrical feature (WGF) algorithm to register computed tomography (CT) images of the head to physical space using the skin surface only, the bone surface only, and various weighted combinations of these surfaces and one fiducial point (centroid of a bone-implanted marker). We use data acquired from 12 patients that underwent temporal lobe craniotomies for the resection of cerebral lesions. We evaluate and compare the accuracy of the registrations obtained using these various approaches by using as a reference gold standard the registration obtained using three bone-implanted markers. The results demonstrate that a combination of geometrical features can improve the accuracy of CT-to-physical space registration. Point-based registration requires a minimum of three noncolinear points. The position of a bone-implanted marker can be determined much more accurately than that of a skin-affixed marker or an anatomic landmark. A major disadvantage of using bone-implanted markers is that an invasive procedure is required to implant each marker. By combining surface information, the WGF algorithm allows registration to be performed using only one or two such markers. One important finding is that the use of a single very accurate point (a bone-implanted marker) allows very accurate surface-based registration to be achieved using very few surface points. Finally, the WGF algorithm, which not only allows the combination of multiple types of geometrical information but also handles point-based and surface-based registration as degenerate cases, could form the foundation of a "flexible" surgical navigation system that allows the surgeon to use what he considers the method most appropriate for an individual clinical situation.

Measurement of intraoperative brain surface deformation under a craniotomy

OBJECTIVE

Several causes of spatial inaccuracies in image-guided surgery have been carefully studied and documented for several systems. These include error in identifying the external features used for registration, geometrical distortion in the preoperative images, and error in tracking the surgical instruments. Another potentially important source of error is brain deformation between the time of imaging and the time of surgery or during surgery. In this study, we measured the deformation of the dura and brain surfaces between the time of imaging and the start of surgical resection for 21 patients.

METHODS

All patients underwent intraoperative functional mapping, allowing us to measure brain surface motion at two times that were separated by nearly an hour after opening the dura but before performing resection. The positions of the dura and brain surfaces were recorded and transformed to the coordinate space of a preoperative magnetic resonance image, using the Acustar surgical navigation system (manufactured by Johnson & Johnson Professional, Inc., Randolph, MA) (the Acustar trademark and associated intellectual property rights are now owned by Picker International, Highland Heights, OH). This system performs image registration with bone-implanted markers and tracks a surgical probe by optical triangulation.

RESULTS

The mean displacements of the dura and the first and second brain surfaces were 1.2, 4.4, and 5.6 mm, respectively, with corresponding mean volume reductions under the craniotomy of 6, 22, and 29 cc. The maximum displacement was greater than 10 mm in approximately one-third of the patients for the first brain surface measurement and one-half of the patients for the second. In all cases, the direction of brain shift corresponded to a "sinking" of the brain intraoperatively, compared with its preoperative position. Analysis of the measurement error revealed that its magnitude was approximately 1 to 2 mm. We observed two different patterns of the brain surface deformation field, depending on the inclination of the craniotomy with respect to gravity. Separate measurements of brain deformation within the closed cranium caused by changes in patient head orientation with respect to gravity suggested that less than 1 mm of the brain shift recorded intraoperatively could have resulted from the change in patient orientation between the time of imaging and the time of surgery.

CONCLUSION

These results suggest that intraoperative brain deformation is an important source of error that needs to be considered when using surgical navigation systems.

Visual assessment of the accuracy of retrospective registration of MR and CT images of the brain

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.

Effects of acceleration, jerk, and field inhomogeneities on vessel positions in magnetic resonance angiography

Blood flow and magnetic field inhomogeneities lead to distortions in MR angiography (MRA) images that present added risk for stereotactic neurosurgical applications. These effects are demonstrated in an MRA image of a model of cerebrovasculature. Analysis of the effects of velocity, acceleration, jerk, and field inhomogeneities on vessel position is presented; results are used to predict vessel shifts for several cerebral blood vessels. The actual encoded position for flowing spins is shown to be a moment-weighted average position. Maximum shift of 3.11 mm was reduced to 0.05 mm when velocity compensation was added. Velocity compensation applied specifically in the phase-encoding direction reduces flow-dependent shifts to the point that they can be safely ignored even if acceleration and jerk are present. Those prescribing and using MRA images for stereotactic applications must be aware of whether compensation is actually applied along the phase-encoding axis when a flow-compensated sequence is used.

An ultrasonic approach to localization of fiducial markers for interactive, image-guided neurosurgery--Part I: Principles

Fiducial markers are reference points used in the registration of image space(s) with physical (patient) space. As applied to interactive, image-guided surgery, the registration of image space with physical space allows the current location of a surgical tool to be indicated on a computer display of patient-specific preoperative images. This intrasurgical guidance information is particularly valuable in surgery within the brain, where visual feedback is limited. The accuracy of the mapping between physical and image space depends upon the accuracy with which the fiducial markers were located in each coordinate system. To effect accurate space registration for interactive, image-guided neurosurgery, the use of permanent fiducial markers implanted into the surface of the skull is proposed in this paper. These small cylindrical markers are composed of materials that make them visible in the image sets. The challenge lies in locating the subcutaneous markers in physical space. This paper presents an ultrasonic technique for transcutaneously detecting the location of these markers. The technique incorporates an algorithm based on detection of characteristic properties of the reflected A-mode ultrasonic waveform. The results demonstrate that ultrasound is an appropriate technique for accurate transcutaneous marker localization. The companion paper to this article describes an automatic, enhanced implementation of the marker-localization theory described in this article.

An ultrasonic approach to localization of fiducial markers for interactive, image-guided neurosurgery--Part II: Implementation and automation

Registration of image space and physical space lies at the heart of any interactive, image-guided neurosurgery system. This paper, in conjunction with the previous companion paper [1], describes a localization technique that enables bone-implanted fiducial markers to be used for the registration of these spaces. The nature of these subcutaneous markers allows for their long-term use for registration which is desirable for surgical follow-up, monitoring of therapy efficacy, and performing fractionated stereotactic radiosurgery. The major challenge to using implanted markers is determining the location of the markers in physical space after implantation. The A-mode ultrasonic technique described here is capable of determining the three-dimensional (3-D) location of small implanted cylindrical markers. Accuracy tests were conducted on a phantom representing a human head. The accuracy of the system was characterized by comparing the location of a marker analogue as determined with an optically tracked pointer and the location as determined with the ultrasonic localization. Analyzing the phantom in several orientations revealed a mean system accuracy of 0.5 mm with a +/- 0.1-mm 95% confidence interval. These tests indicate that transcutaneous localization of implanted fiducial markers is possible with a high degree of accuracy.

Correcting scaling errors in tomographic images using a nine degree of freedom registration algorithm

PURPOSE

Clinical imaging systems, especially MR scanners, frequently have errors of a few percent in their voxel dimensions. We evaluate a nine degree of freedom registration algorithm that maximizes mutual information for determining scaling errors. We evaluate it by registering MR and CT images for each of five patients (patient scaling) and by registering MR images of a phantom to a computer model of the phantom (phantom scaling).

METHOD

Each scaling method was validated using bone-implanted markers localized in the patient images and also intraoperatively. The root mean square residual in the alignment of the fiducial markers [fiducial registration error (FRE)] was determined without scale correction, with patient scaling, and with phantom scaling.

RESULTS

Each scaling method significantly reduced the average FRE (p < 0.05) for MR to CT registration and for MR to physical space registration, indicating that voxel scaling errors were reduced. The greater reduction in scaling errors was achieved using the phantom scaling method.

CONCLUSION

We have demonstrated that a nine degree of freedom registration algorithm that maximizes mutual information can significantly reduce scaling errors in MR.