Methods and apparatus for registering CT-scan data to multiple fluoroscopic images

• US 5,951,475 A
• Filed: 09/25/1997
• Issued: 09/14/1999
• Est. Priority Date: 09/25/1997
• Status: Expired due to Term

• Current Assignee: Integrated Surgical Systems, Inc., International Business Machines Corporation
• Original Assignee: International Business Machines Corporation
• Inventors: Kazanzides, Peter, Taylor, Russell H., Gueziec, Andre Pierre
• Primary Examiner(s): Kamm, William E.
• Assistant Examiner(s): SHAW, SHAWNA JEANNINE
• Application Number:
  US08/936,935
• Time in Patent Office: 719 Days
• Field of Search: 600/425, 600/407, 600/473, 600/414, 600/417, 378/205, 378/206, 378/207, 378/42, 378/190, 606/130, 382/128, 382/132, 382/285, 382/294, 901/14, 356/395-397, 356/375
• US Class Current: 600/425
• CPC Class Codes: A61B 90/36 Image-producing devices or illumination devices not otherwise provided for
• CPC Class Codes: G06T 3/14 Transformations for image registration, e.g. adjusting or mapping for alignment of images
• CPC Class Codes: Y10S 128/922 including image analysis

Abstract

A method and system is disclosed for registering two dimensional fluoroscopic images with a three dimensional model of a surgical tissue of interest. The method includes steps of: (a) generating, from CT or MRI data, a three dimensional model of a surgical tissue of interest; (b) obtaining at least two, two dimensional, preferably fluoroscopic, x-ray images representing at least two views of the surgical tissue of interest, the images containing radio-opaque markers for associating an image coordinate system with a surgical (robot) coordinate system; (c) detecting the presence of contours of the surgical tissue of interest in each of the at least two views; (d) deriving bundles of three dimensional lines that pass through the detected contours; and (e) interactively matching three dimensional points on three dimensional silhouette curves obtained from the three dimensional model with the bundles of three dimensional lines until the three dimensional model is registered within the surgical coordinate system to a predetermined level of accuracy. The step of iteratively matching includes steps of: defining a distance between surfaces of the model and a beam of three dimensional lines that approach the surfaces; and finding a pose of the surfaces that minimizes a distance to the lines using, preferably, a statistically robust method, thereby providing a desired registration between a surgical robot and a preoperative treatment plan.

Claims

1. A method for registering two dimensional fluoroscopic images with a three dimensional model of a surgical tissue of interest, comprising steps of:

(a) generating a three dimensional model of a surgical tissue of interest using a plurality of images of the surgical tissue of interest;

(b) obtaining at least two, two dimensional fluoroscopic images representing at least two views of the surgical tissue of interest, the images containing markers for associating an image coordinate system with a surgical coordinate system;

(c) detecting the presence of first contours of the surgical tissue of interest in each of the at least two views;

(d) deriving bundles of lines in three dimensional space that pass through the detected first contours; and

(e) iteratively matching points in three dimensional space on second contours obtained from the three dimensional model with the bundles of lines in three dimensional space until the three dimensional model is registered within the surgical coordinate system to a predetermined level of accuracy.

2. A method as in claim 1, wherein the step of generating a three dimensional model includes preliminary steps of performing a plurality of scans through the surgical tissue of interest to generate a plurality of image slices, and processing the image slices to generate a three dimensional model of the surface of the surgical tissue of interest.

3. A method as in claim 1, wherein the step of obtaining at least two, two dimensional fluoroscopic images includes steps of positioning a calibration device within a field of view of a fluoroscopic imaging system, the calibration device containing radio-opaque markers having a predetermined shape, and processing the two dimensional images to locate edges of images of the markers.

4. A method as in claim 3, wherein the step of positioning includes a step of translating the calibration device through the field of view using a robot manipulator, the calibration device being translated along a first ruled surface and along a second ruled surface that is generally parallel to the first ruled surface, the first and second ruled surfaces each being generally perpendicular to an optical axis of the fluoroscopic imaging system.

5. A method as in claim 4, wherein the calibration device has a generally circular cylindrical rod-like shape having the radio-opaque markers spaced apart along a longitudinal axis of the calibration device.

6. A method as in claim 3, wherein the step of positioning includes a step of mounting the calibration device within a field of view of a camera that forms a portion of the fluoroscopic imaging system.

7. A method as in claim 4, wherein the calibration device has a first set of radio-opaque markers located in a first plane, and a second set of larger radio-opaque markers located in a second plane, and wherein the first plane is located further from a source of imaging radiation than the second plane and is generally parallel to the second plane, the first and second planes each being generally perpendicular to an optical axis of the fluoroscopic imaging system.

8. A method as in claim 1, wherein the step of iteratively matching includes steps of:

defining a distance between points on the surfaces of the model and a beam of lines in three dimensional space that approach the surfaces; and

finding a pose of the surfaces that minimizes a distance to the lines from the points.

9. A method as in claim 1, wherein the step of iteratively matching includes a step of minimizing a sum of positive measures of a Euclidean distance between each point and its corresponding line.

10. A method as in claim 9, wherein the step of minimizing employs one of a non-linear Levenberg Marquardt minimizer, or a Linear minimizer followed by constrained minimization, using a sum of squared distances between points and lines.

11. A method as in claim 9, wherein the step of minimizing employs a statistically Robust minimizer that weights each squared distance between a point and a line in the sum with a function that favors distances that are similar to a median of the distances.

12. A method as in claim 1, wherein the first contours are two dimensional contours and the second contours are apparent contours.

13. A method as in claim 1, wherein the surgical coordinate system is comprised of a robot coordinate system.

14. A robotically-assisted surgical system, comprising:

a robot having an effector that is controllably positioned within a robot coordinate system;

a first imaging device for obtaining two dimensional radiographic images a tissue of interest, the images containing markers obtained from a calibration device that is located within a field of view of the first imaging device and that is separate from the tissue of interest the markers being used for associating an image coordinate system with the robot coordinate system;

a second imaging device for obtaining images of slices through the tissue of interest;

a data processor having an input coupled to an output of the first and second imaging devices and an output coupled to the robot, said data processor generating a three dimensional model of the tissue of interest from the output of the second imaging device, said data processor further detecting a presence of two dimensional contours of the tissue of interest in each of at least two views of the tissue and deriving bundles of lines in three dimensional space that pass through the detected two dimensional contours;

said data processor operating to iteratively match points in three dimensional space, that are located on three dimensional apparent contours that are associated with a surface of the three dimensional model, with the bundles of lines in three dimensional space until the three dimensional model is registered with the robot coordinate system to a predetermined level of accuracy, wherein said surface apparent contours are boundaries between visible and invisible surface triangles.

15. A robotically-assisted surgical system as in claim 14, wherein said calibration device is adapted for being robotically translated, external to a patient who comprises the tissue of interest, and within a field of view of said first imaging device, said calibration device having a rod-like shape and being comprised of a radiolucent material containing a plurality of radio-opaque markers that are spaced apart along a longitudinal axis of said calibration device.

16. A robotically-assisted surgical system as in claim 15, wherein said calibration device is robotically translated along a first plane and along a second plane that is generally parallel to the first plane, the first and second planes each being generally perpendicular to an optical axis of the first imaging device, wherein said radio-opaque markers are each sized so as to be capable of producing an image when said calibration device is translated within the field of view of the first imaging device.

17. A robotically-assisted surgical system as in claim 14, wherein said calibration device is adapted for being positioned within a field of view of said first imaging device and external to a patient who comprises the tissue of interest, said calibration device being comprised of a radiolucent body portion containing a first set of radio-opaque markers located in a first plane and a second set of larger radio-opaque markers located in a second plane, and wherein the first plane is located further from a source of imaging radiation of the first imaging device than the second plane and is generally parallel to the second plane, the first and second planes each being generally perpendicular to an optical axis of the first imaging device, said calibration device further comprising a plurality a robot reference locations on a surface thereof.

18. A system for registering two dimensional fluoroscopic images with a three dimensional model of a tissue of interest, comprising:

a model processor for generating a three dimensional model of a tissue of interest using as input a plurality of images taken through the tissue of interest;

an imaging system for obtaining, in cooperation with a calibration device, at least two, two dimensional images representing at least two views of the tissue of interest, the images containing calibration markers for associating an image coordinate system with a second coordinate system; and

an image processor having inputs coupled to outputs of said model processor and said imaging system for detecting the presence of first contours of the tissue of interest in each of the at least two views, said image processor operating to derive bundles of lines in three dimensional space that pass through the detected first contours, and further operating to match points on second contours obtained from the three dimensional model with the bundles of lines in three dimensional space until the three dimensional model is registered within the second coordinate system to a predetermined degree of accuracy.

19. A system as in claim 18, wherein said model processor uses a plurality of image slices obtained from a plurality of scans through the tissue of interest, and processes the image slices to generate a three dimensional model of the surface of the tissue of interest.

20. A system as in claim 18, wherein said system further comprises a robot manipulator for translating said calibration device within a field of view of said imaging system, said calibration device comprising said calibration markers, and wherein said image processor processes the two dimensional images to locate images of said calibration markers within said two dimensional images.

21. A system as in claim 18, wherein said system further comprises a robot manipulator for translating said calibration device within a field of view of said imaging system, said calibration device comprising said calibration markers, wherein said robot manipulator is controlled for translating said calibration device within a first ruled surface and within a second ruled surface that is spaced apart from the first ruled surface.

22. A system as in claim 20, wherein a composite image of said calibration markers is comprised of a plurality of images of said calibration device located at a plurality of different locations within the field of view of said imaging system.

23. A system as in claim 18, wherein said calibration device has a generally circular cylindrical rod-like shape, and wherein said calibration markers are comprised of radio-opaque material and are spaced apart from one another along a longitudinal axis of said calibration device.

24. A system as in claim 18, wherein said calibration markers are comprised of radio-opaque material, wherein said calibration device comprises a first set of said calibration markers that are disposed in a first plane and a second set of larger calibration markers that are located in a second plane, wherein said first plane is located closer to a source of imaging radiation of said imaging system than said second plane and is generally parallel to said second plane, and wherein said first plane and said second plane are each generally perpendicular to an optical axis of said imaging system.

25. A system as in claim 18, wherein the second coordinate system is comprised of a robot coordinate system.

26. A system as in claim 18, wherein the tissue of interest is comprised of bone, and wherein said imaging system comprises means for obtaining images of bone through overlying tissue.

27. A method for generating a three dimensional model of a tissue of interest and for aligning the three dimensional model with two dimensional images of the tissue of interest, comprising steps of:

performing a plurality of imaging scans through the tissue of interest to generate a plurality of two dimensional image slices;

processing the image slices to generate a three dimensional model of a surface of the tissue of interest;

obtaining a plurality of two dimensional images representing at least two views of the tissue of interest, the at least two views also comprising images of calibration markers for associating an image coordinate system with a second coordinate system;

detecting first contours of the tissue of interest in each of the at least two views;

deriving bundles of lines in three dimensional space that pass through the detected first contours; and

matching points on second contours obtained from the surface of the three dimensional model with the bundles of lines until the three dimensional model is aligned within the second coordinate system to a predetermined level of accuracy.

28. A method as in claim 27, wherein the step of processing the image slices to generate a three dimensional model is comprised of the steps of:

extracting outer contours of the tissue of interest for each image slice to obtain at least one and typically a plurality of polygonal curves representing contours of the tissue of interest.

replacing the polygonal curves with approximating polygonal curves containing fewer vertices, such that the approximating polygonal curves do not deviate from the original polygonal curves by more than a pre-specified threshold, wherein each polygonal curve is represented using an ordered array of vertices, and each vertex is indexed with its position in the array;

building a surface model comprised of triangles that contains every vertex and every edge of the approximating polygonal curves, the step of building a surface model examining in turn pairs of consecutive image slices and constructing a surface slab that contains the approximating polygonal curves extracted from both image slices; and

combining the surface slabs to form a resulting surface by removing duplicate references to vertices belonging to polygonal curves shared by any two slabs.

29. A method as in claim 28, and further comprising a step of:

approximating the resulting surface with a surface containing fewer triangles by the use of a surface simplification technique.

30. A method as in claim 28, wherein the step of extracting comprises steps of:

selecting a plurality of points within an image slice in the vicinity of a structure of interest;

constructing a polygonal curve linking the selected points; and

modifying the polygonal curve to minimize an expression that combines a measure of curve length, a measure of average curvature of the curve, a measure of an image potential average curvature of the curve, and a measure of an image potential.

31. A method as in claim 30, wherein the image potential equals the squared norm of the image gradient, where if I(x,y) is the image intensity of a pixel of location (x,y), the image gradient grad (I) is a vector whose first coordinate is the derivative of I with respect to x, and whose second coordinate is the derivative of I with respect to y.

32. A method as in claim 28, wherein the step of replacing comprises steps of:

for each polygonal curve, iterating the following steps (a)-(c) until a maximum deviation threshold is respected,

(a) computing the maximum deviation between any vertex and a line segment obtained by joining the last vertex with the first vertex;

(b) if the maximum deviation computed is larger than the threshold, splitting the ordered array of vertices into two arrays; and

(c) considering two polygonal curves defined with the resulting split arrays.

33. A method as in claim 28, wherein the second contours are surface apparent contours defined as boundaries between visible and invisible surface triangles, wherein a visible triangle is one that a ray cast from a focal point of the lines to the triangle centroid makes an obtuse angle with the triangle normal direction, wherein surface apparent contours can be constructed as a set of closed, oriented, non-planar polygonal curves.

34. A method as in claim 27, wherein multiple images of calibration markers are superimposed to form a calibration marker grid, the images being obtained for a plurality of different locations of a calibration device within a field of view of a tissue imaging system.

35. A method as in claim 27, wherein images of calibration markers are processed using Thin Plate Spline functions to determine lines in three dimensional space representing a path of imaging radiation.