Patent ID: 8121380

Claim:
Computerized imaging method for a three-dimensional reconstruction from two-dimensional radiological images of an area of interest, in which two two-dimensional radiological images of a subject are acquired through radiographic means ( 2 , 3 , 3 ′) coupled to a computer equipment, the radiological images being taken successively from different incidence angles relative to the subject, the subject having possibly been deformed between the radiological image acquisitions, characterized in that the following steps are performed: disposing on the surface of the subject, in relation with the area of interest, a set of surface reference marks ( 26 , 26 ′) detectable by a surface sensor ( 4 ) independent of the radiographic means for acquiring three-dimensional surface images of the subject's surface in the computer equipment, the surface reference marks being radio-opaque to be also visible on the acquired radiological images, associating with each acquired radiologic image ( 3 , 3 ′) a corresponding surface image acquired at substantially the same time, implementing an iterative process comprising a phase of incidence angle calculation ( 10 ) and a phase of deformation calculation ( 11 ), the phase of incidence angle calculation being intended to determine ( 15 ) the relative positions of the radiological images and the surface images as well as to perform a spatial positioning ( 16 ) of the radiological images, assuming an absence of deformation of the subject, the phase of deformation calculation being intended to determine ( 19 ) and correct ( 21 ) the interpretation of the radiological images to take into account the deformations of the subject for the relative positions determined in the phase of incidence angle calculation, and in that the criterion for stopping the iteration is based on the deviation ( 23 ) between the 3D coordinates of same reference marks between two successive iterations, that is a previous iteration i−1 and a current iteration i, the current 3D coordinates of the reference marks of a current iteration i being calculated by determining a deformation gradient obtained from the two acquired surface images, applying this gradient to the radiological images of the previous iteration i−1 to obtain current corrected radiological images i, and then calculating the new current 3D coordinates of reference marks from the current corrected radiological images i, the positions of the reference marks of the previous iteration i−1 corresponding to the positions of the reference marks on the initially acquired surface images, and in the case of a first iteration, the radiological images of the previous iteration corresponding to the acquired radiological images.