1. Field of Invention
The invention relates to a method for computer assisted localization of the distal locking holes of an implanted intramedullary nail and to a device for the computer assisted distal locking of an intramedullary nail.
2. Description of Related Art
It is known that one of the most difficult steps of intramedullary nailing of e.g. femoral shaft fractures is distal locking—the insertion of distal interlocking screws, for which it is necessary to know the positions of the distal locking holes (denoted in the following as DLHs) of the intramedullary nail (IMN). The process of locating and inserting the distal interlocking screws is complicated by the nail deformation during insertion. It is known that deformation occurs in several planes due to medio-lateral (ML) and antero-posterior (AP) flexion of the intramedullary nail during its insertion. The reason for the wide variations of the insertion-related intramedullary nail deformation is due to the fact that the nail has to deform to the shape of the medullary canal upon insertion. The shape of the canal varies widely from person to person. It is not possible to predict how the nail will deform accordingly. In a conventional surgical procedure the surgeon depends on intra-operative X-ray means for providing precise positions of the DLH's. It requires positioning the axis of the fluoroscope coaxial to the locking holes so that these holes appear perfectly circular in the images. This is achieved through a trial-and-error method and requires long time X-ray exposure for both the surgeon and patient. It has been reported that the surgeon's direct exposure to radiation for each conventional surgical procedure was 3-30 min, of which 31%-51% was used for distal locking.
The object to target accurately the DLHs with as little as possible X-ray exposure has led to various attempts to develop image-based methods for recovering the positions of the DLHs.
A computer assisted intramedullary rod surgery system is known from WO 03/043485 KIENZLE. In order to determine the position of the nail axis and of the distal locking holes of the intramedullary nail two fluoroscopic images are to be acquired, one in an approximate antero-posterior direction of the bone and a second in a lateral-medial direction of the bone. Subsequently, image processing algorithms are used to calculate the position of the nail axis as well as of the distal locking holes. This known method shows the drawback that the acquisition of the mentioned two fluoroscopic images requires a repeated and cumbersome displacing and rearranging of other surgical apparatus, devices and instruments in order to rotate e.g. the C-arm of the X-ray device from the antero-posterior direction to the lateral-medial direction.
Further, a method for robot-assisted guide positioning for distal locking of intramedullary nails is disclosed in: Ziv Yaniv, Leo Joskowicz, IEEE TRANSACTIONS ON MEDICAL IMAGING, Vol. 24, No. 5 May 2005. This known method requires only one fluoroscopic image, but in a view exactly parallel to the axes of the distal locking holes, such that the distal locking holes appear as circles (and not as ellipses) in the image. In order to correctly adjust the fluoroscope an X-ray technician must use a try-and-move method several times to achieve the correct adjustment of the C-arm of the X-ray device which typically requires 2-6 images to be taken with the X-ray device. The requirement of a fluoroscopic image in a view exactly parallel to the axes of the distal locking holes therefore shows the drawback of a very time-consuming adjustment procedure of the fluoroscope.
On this point the invention intends to provide remedial measures. The invention is based on the objective of providing a method for solving the above problems using one single calibrated fluoroscopic image taken in a lateral-medial direction of the bone only independent of the orientation of the distal locking holes.
The invention solves the posed problem with a method for computer assisted localization of the distal locking holes of an implanted intramedullary nail and with a device for the computer assisted distal locking of an intramedullary nail.
One of the advantages of the present invention is due to the acquisition of only one single image, which allows to avoid a repeated displacing and rearranging of other surgical apparatus, devices and instruments and which minimizes the X-ray exposure for both the surgeon and the patient.
Another advantage is to be seen in the fact that the positioning of the X-ray device is irrespective to the orientation of the distal locking holes allowing an uncomplicated positioning of the X-ray device.
In a preferred embodiment the method according to the invention comprises the following steps:
In a further embodiment the single image (IM) of the distal end section of said intramedullary nail is taken by means of the X-ray device in a lateral-medial direction of said intramedullary nail implanted in a bone.
In another embodiment a data set defining a geometrical model of said distal end section is entered in the data store of the computer under step i). This allows the advantage that the geometrical model may comprise distal locking holes with a non-circular, e.g. an oval cross section and/or e.g. a conical distal end section instead of a circular cylindrical distal end section.
In still a further embodiment the calculation under step vi) is performed by an iterative closest projection point algorithm.
In another embodiment the steps iv)-vi) are combined to an iterative calculation procedure.
In yet another embodiment the method further comprises the step of attaching a dynamic reference base at an essentially rigid body formed by the intramedullary nail, a proximal and a distal fragment of a bone.
In a further embodiment the method further comprises the step of measuring the position and orientation of the dynamic reference base with a position measurement device allowing to establish a three-dimensional coordinate system fixedly connected with the intramedullary nail or with a bone in which the intramedullary nail is implanted.
In yet a further embodiment the method further comprises the step of displaying a graphic representation of the first and second distal locking hole including trajectories of the axes of the first and second distal locking hole allowing a surgeon to accurately drill the holes in a bone that are necessary to insert the locking screws in the distal locking holes.
In another embodiment the X-ray device is a C-arm X-ray device.
In still another embodiment the X-ray device is a fluoroscope.
One embodiment of the invention will be described in the following by way of example and with reference to the accompanying drawings in which:
In the following an embodiment of the method according to the invention is described with reference to
The proximal fragment of the bone 1, the distal fragment of the bone 1 and the intramedullary nail 2 may be treated as three rigid bodies and registered independently. The rigid transformations between these three rigid bodies are obtained from a navigator or position measurement device 6 such as an optoelectronic tracker, a magnetic tracker, or even a medical robot. According to the present invention it is assumed that the fractured bone 1, e.g. femur has already been reduced and the proximal fragment and distal fragment are kept fixed relative to each other at the time of image acquisition. It is also assumed that the intramedullary nail 2 has been inserted till the distal end of the bone 1, e.g. the bone 1 and has been locked proximally by a bone screw so that the assembly, i.e. the bone 1 and the intramedullary nail 2 can be treated as one rigid body 4. A local coordinate system is established on this rigid body 4 through a common dynamic reference base technique using a dynamic reference base 5. In the following description, this patient coordinate system COS is denoted as a three-dimensional coordinate system 7 (A-COS). All computations effected by the computer 8 are done in this reference A-COS.
To relate a pixel in the two-dimensional (2D) projection image IM acquired by means of an X-ray device 9 to A-COS 7, the acquired image IM has to be calibrated for physical projection properties and be corrected for various types of distortion. A weak-perspective pin-hole camera model may be chosen for modeling the C-arm projection of the X-ray device 9. Using such a camera model, a 2D pixel VI is related to a three-dimensional (3D) point VA by following equations:
where ∥•∥ means to calculate the length of a vector and the vectors fA, rA, cA and pI represent the position of focal point 26, the vector along image row increasing direction, the vector along image column increasing direction, and the 2D position of piercing point, respectively. They are projection parameters used to describe the projection properties of the X-ray device 9 and need to be calibrated preoperatively.
Eq. (1) can be used for both forward and backward projections. For example, in order to calculate the direction sA of the forward projection ray of an image point VI, an additional constraint ∥sA∥=1 can be used together with Eq. (1) to solve it. The forward projection ray of point VI is defined by the focal point 26 and the direction sA.
The position of the imaging plane in the three-dimensional coordinate system 7 (A-COS) and the focal length in our camera model is implicitly determined using the calibrated focal point fA and the vectors rA and cA. Any 2D image point VI corresponds to a 3D spatial point IA in this imaging plane, which is the intersection between its forward projection ray and this plane.
The task of the preprocessing is to determine the projection point 32 of the nail tip 13 and the projection points 33;34 of the centers 21,22 of the distal locking holes 11,14 (
Let's denote the initial position of the three points in the axis, i.e., the nail tip 13, the center 21 of the distal DLH, and the center 22 of the proximal DLH as N0(x0, y0), N1(x1, y1), and N2(x2, y2), respectively. The coordinates of these points can be obtained by assuming that the initial position of this axis is parallel to the X axis and by a rough estimation of the depth of the axis using the ratio between the length of line segment D0AD2A in the imaging plane and the true length from the center of the proximal DLH to the nail tip. Given a realization of (θ, tx, ty), the transformed coordinates are:
The projection of these points onto the X axis are:
The optimal values of (θ, tx, ty) are found by minimizing a joint cost function combining the likelihood energy and the prior energy:
where τ is a parameter that controls the relative weighting between these two terms.
The Iterative Closest Projection Point Algorithm for Pose Recovery of DLHs is schematically illustrated in
Observing the wide variations of the nail deformation, we are apt for recovering the poses of the two DLHs one by one. In the following descriptions, we concentrate on pose recovery of the distal DLH 14 (
Given the estimated nail axis 19, the initial transformation between the local COS of the geometrical model of the DLH and A-COS can be obtained by taking the estimated nail axis 19 as the u axis and the normal of the imaging place as the v axis of the local COS. All points defined in the local COS of the geometrical model of the DLH can then be transformed to A-COS using this initial transformation. The task of recovery of the DLHs is to resolve the rotation a and the translation δ of the DLH along the axis by fitting the geometrical model of the DLH to the image, which is solved iteratively as follows.
Let's denote E be a set of NE detected 2D edge pixels {e1, e2, . . . , eNE} of the DLH projection. Further denote Mt-1 be a set of NM model point {m0t-1, m1t-1, . . . mNmt-1} at iteration step t-1. Now in the iteration step t, we perform following steps:
Simulating X-ray projection: In this step, we simulate the X-ray projection of the geometrical model of the DLH to remove invisible points. Let pt-1 be a set of Np 2D projection points {p1t-1, p2t-1, . . . pNpt-1} obtained by simulating X-ray projection of 3D model into the image. Normally Np<<NM. Thus, for each 2D projection point ptt-1, we know its associated 3D model point mtt-1.
Find closest projection point: In this step, we try to find the closest neighbor edge pixel et of each 2D model projection point Ptt-1.
Establishing 3D-2D correspondence: For each 2D matched pairs (et, ptt-1), calculate the forward projection ray BPt of the 2D edge pixel et. Then for the ray BPt, calculate a 3D point pair PP1t-1=(bett-1, mtt-1), where bett-1 is a point on the line BP1 that is closest to the 3D model point mtt-1 of the model projection point ptt-1.
Estimating pose: For all calculated 3D point pairs PPS(t-1)={PPtt-1}, find an optimal local solution of all pose parameters by minimizing following cost function:
where T(α(t-1),δf(t-1)) is a constrained transformation around the axis of the nail determined by rotation angle α(t-1) and translation δ(t-1).
Updating pose: Update the pose of all model points {m1t-1, m2t-1, . . . , mNmt-1} by T(α(t-1)), δ(t-1)) to {m1t, m2t, . . . , mNmt}.
These steps are repeated until all pose parameters are converged.
In another embodiment the distal end section 25 of the intramedullary nail 2 as illustrated in
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/CH2006/000432 | 8/15/2006 | WO | 00 | 4/15/2008 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2008/019510 | 2/21/2008 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
5584838 | Rona et al. | Dec 1996 | A |
6503249 | Krause | Jan 2003 | B1 |
6856827 | Seeley et al. | Feb 2005 | B2 |
6917827 | Kienzle, III | Jul 2005 | B2 |
20030004518 | Perren et al. | Jan 2003 | A1 |
20040111024 | Zheng et al. | Jun 2004 | A1 |
20040143340 | Tuma et al. | Jul 2004 | A1 |
20040181149 | Langlotz et al. | Sep 2004 | A1 |
20050027304 | Leloup et al. | Feb 2005 | A1 |
20050080427 | Govari et al. | Apr 2005 | A1 |
20050203539 | Grimm et al. | Sep 2005 | A1 |
20050261700 | Tuma et al. | Nov 2005 | A1 |
20060064106 | Fernandez | Mar 2006 | A1 |
20060122495 | Kienzle, III | Jun 2006 | A1 |
20060173293 | Marquart et al. | Aug 2006 | A1 |
Number | Date | Country |
---|---|---|
03043485 | May 2003 | WO |
2005000129 | Jan 2005 | WO |
Entry |
---|
Yaniv, Z. et al.; “Precise robot-assisted guide positioning for distal locking of intramedullary nails;” Database Inspec (online); The Institution of Electrical Engineers; Stevenage, Great Britain; May 2005; XP002433206. |
Number | Date | Country | |
---|---|---|---|
20080281334 A1 | Nov 2008 | US |