Method for reconstructing a 3D image from 2D X-ray images
The present invention relates to a method for reconstructing a 3D image from 2D X-ray images acquired with an X-ray imaging system, said method comprising the steps of: a) receiving a set of 2D X-ray images of a region of a patient with said X-ray imaging system, b) computing an initial 3D image within the coordinate system of the X-ray imaging system by using at least part of said 2D X-ray images with their respective projective geometry data; c) projecting said initial 3D image on at least part of said 2D X-ray images and adjusting the respective projective geometry data of said images, said adjustment comprising registration of said images with the projection of the initial 3D image using an image-to-image registration technique; d) computing an updated 3D image using the complete set of 2D X-ray images with their respective adjusted projective geometry data.
1. A method for reconstructing a 3D image from 2D X-ray images acquired with an X-ray imaging system, said method comprising the steps of:
a) receiving a set of 2D X-ray images of a region of a patient with said X-ray imaging system,
b) computing an initial 3D image within the coordinate system of the X-ray imaging system by using the complete set of 2D X-ray images with their respective projective geometry data;
c) projecting said initial 3D image on each of said 2D X-ray images and adjusting the respective projective geometry data of said 2D X-ray images by registration of each 2D X-ray image with the projection of the initial 3D image in a predefined region of interest using an image-to-image registration technique;
d) computing an updated 3D image using the complete set of 2D X-ray images with their respective adjusted projective geometry data.
2. The method of claim 1 , wherein one or several steps (b) to (d) are iterated.
3. The method of claim 1 , wherein the calculated adjusted projective geometry data of each image are used as the updated nominal projective geometry data for further 3D image reconstruction.
4. The method of claim 1 , wherein the evolution of the projective geometry data is tracked.
5. The method of claim 1 , wherein the image-to-image registration technique implemented in step (c) determines a rigid transform consisting of one 2D translation and one rotation.
6. The method of claim 1 , wherein the image-to-image registration technique implemented in step (c) determines an affine transform or a warping transform that displaces the four corners of the image independently.
7. The method of claim 1 , wherein a contrast criterion is calculated on the reconstructed 3D image and the method is stopped when said criterion is above a given threshold.
8. The method of claim 1 , wherein the definition of the region of interest comprises the following steps:
searching the initial computed 3D image for the voxels that have a high contrast after filtering noise,
identifying the centroid of said high-contrast voxels, and
projecting said centroid on each 2D image with a bounding box around the resulting projected point that has a determined size that is less than the size of the initial 3D image in each dimension.