Method for determining the position and orientation of an object, especially of a catheter, from two-dimensional X-ray images

• US 7,801,342 B2
• Filed: 06/20/2006
• Issued: 09/21/2010
• Est. Priority Date: 06/21/2005
• Status: Expired due to Fees

• Current Assignee: Siemens Healthcare GMBH (Siemens AG)
• Original Assignee: Siemens AG
• Inventors: Rahn, Norbert, Boese, Jan, Sandkamp, Bernhard
• Primary Examiner(s): Le; Vu
• Assistant Examiner(s): Allison; Andrae S
• Application Number:
  US11/471,123
• Publication Number:
  US 20060285638A1
• Time in Patent Office: 1,554 Days
• Field of Search: 382/128, 382/132, 382/154, 600/425, 600/407, 600/424, 378/21, 378/62, 378/64, 378/63
• US Class Current: 382/128
• CPC Class Codes: A61B 2090/364 Correlation of different images or relation of image positions in respect to the body
• CPC Class Codes: A61B 2090/376 using X-rays, e.g. fluoroscopy
• CPC Class Codes: A61B 34/20 Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
• CPC Class Codes: A61B 6/12 Arrangements for detecting or locating foreign bodies
• CPC Class Codes: A61B 90/36 Image-producing devices or illumination devices not otherwise provided for

Abstract

To determine the position and orientation of an object in a x-ray image the x-ray image is initially pre-processed. A three-dimensional template (data record) of the object based on the known constructional features of the object is created. Three parameters for the position and for the orientation respectively are modified iteratively. The three-dimensional template with the parameters for position and orientation is projected in each case onto a two-dimensional plane and the created image is compared on the basis of the generation of a degree of similarity with the pre-processed x-ray image.

Claims

1. A method for determining a three-dimensional position and angular orientation of a catheter from two-dimensional x-ray images, comprising:

a) generating a first x-ray image;

b) pre-processing the first x-ray image to accentuate the imaging structures of the catheter in the x-ray image;

c) generating a three-dimensional template for the catheter based on known dimensional characteristics of the catheter;

d) initializing three parameters for the position and three parameters for the angular orientation of the catheter relative to the pre-processed x-ray image;

e) generating, by a data processing unit, a two-dimensional projection of the template with reference to the six parameters for position and angular orientation of the catheter;

f) comparing, by a data processing unit, the two-dimensional projection with the pre-processed x-ray image;

g) determining, by a data processing unit, a similarity between the two-dimensional projection and the pre-processed x-ray image based on the comparison;

h) aborting, by a data processing unit, the procedure if the similarity exceeds a predetermined threshold and establishing the initialized position and orientation parameters as the actual position and orientation of the catheter in the x-ray image;

i) repeating steps e) through h) in an optimization loop until the degree of similarity exceeds a predetermined threshold,

j) obtaining, by a data processing unit, a first projection line based on the six position and orientation parameters obtained with the first x-ray image via a back-projection of the catheter from the first x-ray image into the three-dimensional space;

k) obtaining, by a data processing unit, a second projection line based on the six position and orientation parameters based on the second x-ray image via a back-projection of the object from the second x-ray image into the three-dimensional space;

l) determining, by a data processing unit, the location with the smallest quadratic distance to the two projection lines; and

m) determining, by a data processing unit, an intersection point that locates the catheter in three-dimensional space.

2. The method in accordance with claim 1, wherein the three parameters for the position of the catheter correspond to the three translational degrees of freedom of the catheter and the three parameters for the angular orientation of the catheter correspond to the three rotational degrees of freedom of the catheter.

3. The method in accordance with claim 1, wherein the pre-processing of the first x-ray image comprises high-pass filtering or conversion of the x-ray image into binary format.

4. The method in accordance with claim 1, further comprising the step of:

generating a second x-ray image from a different recording angle than the first x-ray image,

pre-processing the second x-ray image for accentuating the imaging structures of the catheter object in the second x-ray image, and

further analyzing the second x-ray image in accordance with steps d) through i).

5. The method in accordance with claim 3, wherein the six position and orientation parameters based on the first x-ray image are considered in the determination of the six new position and orientation parameters associated with the second x-ray image.

6. The method in accordance with claim 1, wherein the smallest distance between the two projection lines is used as a measure of the quality of the location determination of the catheter.

7. The method in accordance with claim 1, wherein the three-dimensional template:

identifies a part of the catheter that remains rigid and a part that can be bent,

identifies further parameters that record the bending of the object, and

including the further identified parameters in the optimization loop of steps d) through i).

8. An X-ray imaging system comprising:

an x-ray device for providing x-ray images including a first x-ray image, the first x-ray image is pre-processed to accentuate the imaging structures of the catheter in the x-ray image;

a data processing unit that stores a three-dimensional data record of contours of a catheter to be automatically detected from the x-ray images, wherein the data processing unit;

a) initializes three parameters for the position and three parameters for the angular orientation of the catheter relative to the pre-processed x-ray image;

b) generates a two-dimensional projection of the template with reference to the six parameters for position and angular orientation of the catheter;

c) compares the two-dimensional projection with the pre-processed x-ray image;

d) determines a similarity between the two-dimensional projection and the pre-processed x-ray image based on the comparison;

e) aborts the procedure if the similarity exceeds a predetermined threshold and establishing the initialized position and orientation parameters as the actual position and orientation of the catheter in the x-ray image; and

f) repeats steps b) through e) in an optimization loop until the degree of similarity exceeds a predetermined threshold

g) obtaining a first projection line based on the six position and orientation parameters obtained with the first x-ray image via a back-projection of the catheter from the first x-ray image into the three-dimensional space;

h) obtaining a second projection line based on the six position and orientation parameters based on the second x-ray image via a back-projection of the object from the second x-ray image into the three-dimensional space;

i) determining the location with the smallest quadratic distance to the two projection lines; and

j) determining an intersection point that locates the catheter in three-dimensional space.

9. The imaging system in accordance with claim 8, wherein the three parameters for the position of the catheter correspond to the three translational degrees of freedom of the catheter and the three parameters for the angular orientation of the catheter correspond to the three rotational degrees of freedom of the catheter.

10. The imaging system in accordance with claim 8, wherein the pre-processing of the first x-ray image comprises high-pass filtering or conversion of the x-ray image into binary format.

11. The imaging system in accordance with claim 8, wherein the data processing unit further comprises:

generating a second x-ray image from a different recording angle than the first x-ray image,

pre-processing the second x-ray image for accentuating the imaging structures of the catheter object in the second x-ray image, and

further analyzing the second x-ray image in accordance with steps a) through f).