|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Technical Developments |
1 From the Mallinckrodt Institute of Radiology (B.R.W., K.T.B.) and Department of OtolaryngologyHead and Neck Surgery (M.W.S.), Washington University School of Medicine, 510 S Kingshighway Blvd, St Louis, MO 63110. Received January 10, 2001; revision requested March 5; revision received April 3; accepted May 1. Supported in part by the Mallinckrodt Institute of Radiology and Electronic Radiology Lab, National Institute on Deafness and Other Communication Disorders (grant R01-DC00581), and the Whitaker Foundation (Biomedical Engineering Program). Address correspondence to B.R.W. (e-mail: whitingb@mir.wustl.edu).
| ABSTRACT |
|---|
|
|
|---|
Index terms: Computed tomography (CT), image processing, 21.12115, 21.12117, 21.1215 Computed tomography (CT), technology, 21.12115, 21.12117, 21.1215 Computed tomography (CT), three-dimensional, 21.12117 Ear, CT, 21.12115, 21.12117, 21.1215 Ear, labyrinth, 21.12115, 21.12117, 21.1215 Phantoms Radiography, digital, 21.1215
| INTRODUCTION |
|---|
|
|
|---|
Either computed tomography (CT) or computed radiography (CR) can be used to determine the position of cochlear implants in an individuals ear, but neither is sufficient to locate implants to submillimeter accuracy in three-dimensional (3D) space. In fact, the ideal imaging technique would have the CT attributes of 3D location information and sensitivity to low-contrast objects along with the high spatial resolution of CR. The purpose of our study was to evaluate a method that coregisters CT and CR images with use of commonly available diagnostic equipment.
| Materials and Methods |
|---|
|
|
|---|
For CT data, we found that a threshold value of 2,500 HU reliably selected voxels associated with metal objects. An automatic search for all corresponding voxels was made to form a list of the x-y coordinates and section number for metal points. This list was manually edited to remove points associated with the interior of the receiver-stimulator package and to identify an initial estimate of the position of the electrode tip and the entry point of the lead wires into this package. The accuracy requirements for these estimates were very modest; varying the inputs by several millimeters (tens of voxels) resulted in the same final parameter values. Typical numbers of metal voxel points used were 50,000 points for 0.1-mm voxels. Because the CT gantry was tilted for most patients by some angle, the position indexes (i,j,k) corresponded to a nonorthogonal coordinate system and had to be converted into x-y-z coordinates in real space with an appropriate trigonometric transformation as follows:
For the CR data (Fig 1), metal electrode shadows corresponded to substantial negative modulation over a 5 x 5-pixel (0.5-mm2) area on a very nonuniform background (due to anatomic structures), which precluded simple thresholding. Therefore, the image was processed by applying a high pass filter of variable kernel sizes and manually identifying the resulting matched-filter points associated with electrodes (Fig 2). Again, the transformation calculation was not sensitive to the selection of individual points. Typically, a list of 5,000 pixels was generated.
|
|
|
Patient Data Acquisition
Data sets of CT and CR images were collected in 12 patients (six men and six women; age range, 3076 years; mean age, 57.8 years) following cochlear implantation surgery (8). As part of the clinical care, an electrode array approved by the U.S. Food and Drug Administration for clinical use (Nucleus 22) was implanted in each patient. The institutional human studies committee approved the participation of the patients in the present research study, and informed consent was obtained. Both data sets were acquired within 2 hours on the same day, with the CR examination followed by the CT examination. The equipment we used is typical of that available in general hospitals.
Imaging Protocols
CT scans were obtained with helical CT scanners (Somatom Plus S or Somatom Plus 4; Siemens Medical Systems, Iselin, NJ) with use of a high-spatial-resolution head protocol (120 kVp, 215 mA, 1-mm collimation, 1 mm/sec table speed, 1 second per rotation) with submillimeter reconstruction (0.1-mm isotropic voxels, extended Hounsfield scale). The scanner gantry was tilted between -7° and +15° to align image sections with the modiolar axis of the cochlea in each patient. Image sections were reconstructed; images and raw projection data were saved on an optical disk and were transmitted over the network to research workstations.
Digital radiographic images were collected by using a protocol defined by Cohen et al (9). A head position similar to that in a modified Stenver view was used; that is, the head was positioned such that the central ray traveled parallel to the modiolar axis of the cochlea. An x-ray neurologic unit (Neurostand; Siemens Medical Systems) was used, with a 915-mm focal spotto-detector distance, 0.5-mm focal spot, Bucky tray with a moving 9:1 grid, and table-to-cassette distance of 1.75 inches (4.4 cm). Set-up images were obtained with either screen-film film (Lanex Regular, T-MAT G; Eastman Kodak, Rochester, NY) or computed radiographic film (Ektascan storage phosphor reader model 400 with an 18 x 24-cm storage phosphor high-resolution cassette; Eastman Kodak). When the patient was positioned correctly, the electrode array appeared as a nominal circular pattern, and each electrode was distinctly separate from the other electrodes. The exposure parameters for data collection with the storage phosphor plate were 81 kVp and 80 mAs. CR images were used as final radiographs for analysis.
| Results |
|---|
|
|
|---|
For the 12 patients, the mean of all individual average distances from projected CT points to the nearest CR point was 0.184 mm, with an SD of 0.044 mm. Anatomic landmarks (eg, modiolar axis, apical turn) on the CT image sets were projected onto the CR image. The images were then analyzed for the implant insertion length, number of electrodes in the cochlea, and presence of kinks or compressions of the array. This information was valuable for programming individual patients speech processors to optimize speech recognition. For example, the assignment of frequency boundaries of incoming sound to electrodes needs to be an appropriate compromise between matching the estimated characteristic frequency of neurons near each electrode and the number of electrodes that provide different pitch percepts to deliver acoustic cues in each frequency range (eg, first and second formant speech energy) (5,10,11). For less than full insertion of the electrodes into the inner ear, it is important to know postoperatively how many electrodes are inside the inner ear so that only those will be stimulated. Because there is a modest correlation between word recognition and the 3D functional insertion depth of the electrode array as a percentage of the total cochlear length (5), knowledge of this depth will provide insight into a patients speech recognition ability as a basis for auditory training. Electrodes that are kinked cause an unpleasant sound when stimulated, and they need to be removed from the speech processor program. This information is valuable for programming in children.
| Discussion |
|---|
|
|
|---|
With spatial resolution of four line pairs per millimeter, CR can readily depict these high-contrast metal objects, but overlying shadows often obscure low-contrast anatomy and structures, with the result that 3D relationships of objects cannot be detected on radiographs. On the other hand, CT provides excellent information about the 3D location of low-contrast objects but, unfortunately, with limited spatial detail (Fig 4). The spatial resolution at CT is governed by the acquisition apertures: minimum longitudinal collimation is usually set to 1 mm, while in-plane detector arrays have an effective aperture of typically greater than 0.5 mm (eg, 768 detectors subtending a field of view of 500 mm gives an effective aperture of 0.68 mm at the isocenter of our scanner). With standard filtered back-projection reconstruction, the cochlear implant electrode signals result in images blurred 0.5 mm or more, with significant blooming or spread of the metal objects into surrounding imagery. While the centroids of symmetric objects can be determined with accuracies that approach 0.3 mm, attenuation numbers (Hounsfield units) are severely distorted in the vicinity of electrode metal for more than 1 mm. (In fact, individual electrodes are not depicted at CT.)
|
|
Initially, we unsuccessfully attempted to determine the perspective transformation linking x-ray and CT data coordinate systems by using external fiducial markers attached to the patients skin. The markers consisted of six 0.5-mm-diameter platinum spheres, mounted several millimeters apart on a thin plastic base, that were positioned on the skin over a patients zygomatic arch on the side of the ear with the cochlear implant. The markers remained attached during the two postoperative imaging examinations. Patients were at different positions for the examinations, however, being upright for the x-ray examination and supine for the CT data acquisition. Because the skin is somewhat mobile and does not provide a solid anchor point for physical attachment, the fiducial markers moved slightly. Originally, we planned to uniquely identify six reference points on each image and determine the parameters of the perspective transformation, by using the 12 equations for the measured radiograph coordinates to solve for nine unknowns. Because of the poor conditioning of the inversion problem in the perspective geometry, small imprecision in coordinate values yielded large errors in the transformation estimate, which precluded robust transformation determination. For the method to work, it was clear that reference fiducials must be attached in an immovable fashion, extend over several centimeters to provide a good "leverage arm" for the location determination equations, and consist of multiple points to achieve better statistical averaging. Our key insight was that the implant array itself was well suited as a fiducial marker system.
The electrode array and its receiver-stimulator package possess enough asymmetry and spatial extension to allow unambiguous transformations. The electrodes and receiver hardware are surgically implanted in a patients skull and usually remain immobile thereafter. The titanium package that surrounds the receiver-stimulator electronics is implanted above and behind the middle and inner ear at a distance of several centimeters. Lead wires extend from this package through the middle ear and travel inside the electrodes in the array, which extends as much as 25 mm inside the cochlea. In the implant array, for example, the centers of the 22 electrodes and 10 supporting rings are manufactured to be 0.75 mm apart, and the platinum-iridium wires are encircled by the ring electrodes. Given the 0.30-mm electrode length, the distance between electrode edges is 0.45 mm, and the total length of the array is roughly 25 mm (15,16). Thus, a fixed set of extended reference objects is present on the images for each modality.
With use of implants as internal markers, our method combines the high spatial resolution of radiographic data with 3D structures from CT data to help locate electrode arrays in the cochlea. Anatomic features can be selected in CT data and projected onto radiographic images that are familiar to radiologists. Our method potentially can be used to provide positional information to model the electrical characteristics, to determine the 3D orientation of an individuals inner ear with external landmarks, and to facilitate a correlative evaluation between patient performance and electrode position. In addition to the cochlear implant application, our technique might be extended to other medical applications that require the coregistration of 3D and 2D imaging data, such as CT angiography, oral maxillary surgery, conformal radiation therapy, and imaging-guided surgical navigation.
|
|
|
| APPENDIX |
|---|
|
|
|---|
Perspective Transformations
To coregister CT and CR images, a mathematic function is needed that will relate the 3D coordinates of a point in the CT volume to the coordinates of the same point on a CR image. Inverse mapping is also possible (ie, construction of a perspective ray, or line) in the 3D CT space that connects a point on the CR image to the x-ray focal spot. A radiographic image is formed from a 3D object by means of a perspective projection that is characterized by nine parameters: three coordinates {x0CT, y0CT, z0CT} for the center of the CT coordinate system relative to the x-ray coordinate system, three coordinates {xFS, yFS, zFS} of the x-ray focal spot in the x-ray coordinate systems, and three Euler angles {
,
,
} (21) that rotate the CT coordinate axes into alignment with the x-ray coordinate axes (Fig 5). Operationally, CT coordinates must first be transformed into the physical 3D coordinate system of the x-ray source or detector through a series of rotations and translations, expressed as the following:
Algorithm for Determining Transformation Parameters
If the coordinate values of four nonplanar noncolinear corresponding points on CT and CR images were known exactly, Equations (A1)(A5) could be solved for the nine parameters needed. Unfortunately, there are limited unique corresponding points that can be identified between the images. The closest candidates are the tip of the electrode and the junction of the lead wires with the receiver-stimulator package (Fig 1). To overcome this, an algorithm was devised that selected projection parameters based on minimization of the ensemble distance of all projected CT electrode points to the nearest electrode point on the CR image, without requiring a specific one-to-one correspondence between any individual points.
Specifically, the cost function or error metric that we seek to minimize is the sum over all CT points of the minimum distance between the projection of each CT data point on the CR image to the nearest metal point on the CR image (Fig 6). This minimum distance was determined by using a computer spiral search pattern (Fig 7), checking neighboring CR pixels for the presence of a metal point, and terminating the search when one point was located. This approach has the characteristic that as the alignment becomes better and search distances become smaller, the convergence rate speeds up. A MATLAB function fminsearch, which incorporates a Nelder-Mead unconstrained nonlinear gradient-search algorithm to minimize the cost function (22), was used to determine the best-fit condition. It was found that convergence of the search routine depended on initial values of parameters; in fact, the algorithm could converge to an incorrect local minimum under some conditions. Therefore, one of the authors (B.R.W.) selected initial parameters to visually place the electrode tip and lead-to-receiver connection point in an approximate alignment between the CR and projected CT data. Convergence typically took less than 5 minutes. After convergence to a transformation set, typical distances to the nearest CR metal point were less than 2 pixels (0.2 mm) per voxel.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Author contributions: Guarantors of integrity of entire study, B.R.W., K.T.B., M.W.S.; study concepts and design, B.R.W., K.T.B., M.W.S.; literature research, B.R.W., K.T.B., M.W.S.; clinical studies, B.R.W., K.T.B., M.W.S.; experimental studies, B.R.W., K.T.B., M.W.S.; data acquisition, data analysis/interpretation, B.R.W., K.T.B., M.W.S.; statistical analysis, B.R.W., K.T.B., M.W.S.; manuscript preparation, manuscript definition of intellectual content, manuscript editing, manuscript revision/review, manuscript final version approval, B.R.W., K.T.B., M.W.S..
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
B.M. Verbist, R.M.S. Joemai, W.M. Teeuwisse, W.J.H. Veldkamp, J. Geleijns, and J.H.M. Frijns Evaluation of 4 Multisection CT Systems in Postoperative Imaging of a Cochlear Implant: A Human Cadaver and Phantom Study AJNR Am. J. Neuroradiol., August 1, 2008; 29(7): 1382 - 1388. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. M. Verbist, J. H. M. Frijns, J. Geleijns, and M. A. van Buchem Multisection CT as a Valuable Tool in the Postoperative Assessment of Cochlear Implant Patients AJNR Am. J. Neuroradiol., February 1, 2005; 26(2): 424 - 429. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. A. Mecca, W. Wagle, A. Lupinetti, and S. Parnes Complication of Cochlear Implantation Surgery AJNR Am. J. Neuroradiol., November 1, 2003; 24(10): 2089 - 2091. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| RADIOLOGY | RADIOGRAPHICS | RSNA JOURNALS ONLINE |