Abstract:
A method to process radiological images is provided. The method comprises partitioning a radiological image of a region to be treated into a superimposition of layers, the region to be treated comprising at least one first structure and a second structure, wherein one layer solely comprises part of the first structure to be isolated from the remainder of the image, the layer solely comprising that part of the first structure to be isolated from the remainder of the image being determined by means of a parametric model of the first structure. The method further comprises determining an image of the region to be treated from the layering thus obtained, in which the isolated part of the first structure is omitted.

Description:
BACKGROUND OF THE INVENTION 
     1. Field of the Invention 
     Embodiments of the present invention relate to the field of medical imaging, and more particularly to the field of radiology. 
     2. Description of the Prior Art 
     Vascular interventional radiology includes procedures conducted under monitoring by imaging, and notably allows the treatment of arteriosclerosis. 
     Arteriosclerosis causes local narrowing of the diameter of a vessel and also causes obstruction of the vessel, which hinders the circulating blood flow. Angioplasty consists of deploying an inflatable balloon at the point of narrowing for the purpose of widening the diameter of the vessel. An endovascular prosthesis (or stent) can also be implanted in the vessel to keep the vessel open. 
     During angioplasty procedure, the surgeon inserts a guide wire into the vessel to bring a catheter carrying the inflatable balloon or stent as far as the narrowed region. The surgeon must position the stent at a precise point. 
     To monitor conducting of the procedure, the surgeon visualizes a radiological image of the treated region. The radiological image is an image acquired in real time, which enables the surgeon to visualize the positioning of the different instruments used in the region to be treated. 
     In order to be able to control the positioning of the prosthesis in the region to be treated, the surgeon must have an image available in which the stent can be visualized. Yet in the image obtained after deployment of the stent, the guide wire is superimposed over the stent. Therefore, the presence of the guide wire hampers visualization of the stent. However this guide wire must remain in place until the end of the procedure for medical reasons. 
     Because the stent is a thin mesh having very low contrast compared with the remainder of the image, and since the guide wire has greater contrast, the guide wire may mask some parts of the stent which are of interest to the surgeon. Therefore, for the surgeon to have access to a maximum amount of information, the guide wire needs to be deleted so that only the stent is seen in the region to be treated. 
     Techniques are known which allow objects in an image to be removed without impairing the visual aspect of the image. These techniques are known as inpainting. However, these techniques proceed by interpolating the image around the object to reconstruct the removed part. This technique results in a loss of information which is incompatible with diagnosis in the medical sphere. 
     SUMMARY OF THE INVENTION 
     According to one embodiment of the present invention, method to process radiological images is provided. The method comprises partitioning a radiological image of a region to be treated into a superimposition of layers, the region to be treated comprising at least one first structure and a second structure, wherein one layer solely comprises part of the first structure to be isolated from the remainder of the image, the layer solely comprising that part of the first structure to be isolated from the remainder of the image being determined by means of a parametric model of the first structure. The method further comprises determining an image of the region to be treated from the layering thus obtained, in which the isolated part of the first structure is omitted. 
     According to an alternate embodiment of the present invention, a medical imaging system is provided. The medical imaging system comprises: a source configured to emit a beam of rays; a detector positioned facing the source and configured to detect the rays emitted by the source; a support positioned between the source and the detector; a storage unit; an interface unit; and a processing unit. The processing unit it configured to: partition a radiological image of a region to be treated into a superimposition of layers, the region to be treated comprising at least one first structure and a second structure, wherein one layer solely comprises part of the first structure to be isolated from the remainder of the image, the layer solely comprising that part of the first structure to be isolated from the remainder of the image being determined by means of a parametric model of the first structure; and determine an image of the region to be treated from the layering thus obtained, in which the isolated part of the first structure is omitted. 
    
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
       Other characteristics and advantages will become further apparent from the following description which is solely illustrative and non-limiting and is to be read with reference to the appended drawings in which: 
         FIG. 1  illustrates a medical imaging system according to an embodiment of the present invention; 
         FIG. 2  schematically illustrates a method according to an embodiment of the present invention; 
         FIG. 3  schematically illustrates a layer P according to an embodiment of the present invention; 
         FIG. 4  schematically illustrates a superimposition of three layers which model an image of a region to be treated according to an embodiment of the present invention; 
         FIG. 5  schematically depicts a radiological image of a region to be treated according to an embodiment of the present invention; 
         FIG. 6  schematically illustrates a layer estimated in accordance with an embodiment of the present invention; 
         FIG. 7  schematically illustrates the image of the region to be treated, in which the portion of guide wire has been removed, according to an embodiment of the present invention; 
         FIG. 8  is a diagram of the guide wire according to an embodiment of the present invention. 
     
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
     In  FIG. 1 , the radiological imaging equipment comprises a source  1  able to emit a beam  2  of X-rays, a detector  3  arranged facing the source  1  and configured to detect the rays emitted by the source  1 , a support  8  arranged between the source  1  and the detector  3 , a processing unit  4 , a storage unit  5  and an interface unit  6 . 
     The support  8  is intended to receive a patient  7  in whom the surgeon is to conduct a procedure such as an angioplasty to treat atherosclerosis. 
     The processing unit  4  is configured to control the emission of X-rays by the source  1  and movement of the source  1  relative to the detector  3 . 
     Also, the processing unit  4  is configured to command reading of an image by the detector  3  and to receive data acquired by the detector  3 . 
     The processing unit  4  is for example one or more computers, one or more processors, one or more microcontrollers, one or more micro-computers, one or more programmable logic controllers, one or more application-specific integrated circuits, other programmable circuits, or other devices which include a computer such as a work station. 
     The processing unit  4  is coupled with the storage means  5  which may be integrated in or separate from the processing unit  4 . These means can be formed of a hard disk or any other removable storage means (CD-ROM, disk, etc.). These storage means  5  can be used to store an acquired or processed radiological image of the region to be treated. They may be a ROM/RAM memory of the processing unit  4 , a CD-ROM, USB key, memory of a central server. The processing unit  4  may comprise a reading device (not shown) e.g. a disk drive or CD-ROM drive, to read the instructions of a method to process radiological images (described below) from an instruction medium (not shown) such as a floppy disk or CD-ROM. As a variant, the processing unit  4  executes the instructions of the processing method (described below) stored in firmware (not shown). 
     The interface unit  6  comprises a display device  9 . The interface unit  6  provides the surgeon with means to control the procedure. 
     More precisely, during procedure, the surgeon is able to see the vascular structures of the operated region on the display device  9 . 
     The display device  9  is for example a computer screen, a monitor, flat screen, plasma screen or any commercially available display device. The display device  9  enables the surgeon to visualize the vascular structures of the operated region. 
       FIG. 2  schematically illustrates a method according to an embodiment of the present invention. 
     The image processing method uses the principle that a radiological image by nature is transparent, and that it can be modelled as a superimposition of several layers. 
     A “layer” is a plane in which each pixel of the plane corresponds to an absorbance value μ such that: 
             μ   =     log   ⁡     (       I   0     I     )             
in which I 0  is the intensity of the transmitted X-ray beam and in which I is the intensity of the X-ray beam incident to the layer under consideration.
 
       FIG. 3  schematically illustrates a layer P with the incident beam I and transmitted beam I 0 . 
     Starting with the assumption that it is possible to model a radiological image by superimposing several layers, for each pixel, the intensity of the beam transmitted by all the layers is expressed according to Beer-Lambert&#39;s law, given by: 
               I   image     =     I   ·     ⅇ     -       ∑   i     ⁢     μ   i                   
in which I is the intensity of incident X-ray beam, I image  is the intensity of the X-ray beam after passing through each layer i, and in which μ i  is absorbance.
 
       FIG. 4  schematically illustrates the superimposition of three layers which model an image of the region to be treated. 
     If the radiological image is the image of a region to be treated comprising at least one instrument and a vascular prosthesis, the X-rays will pass through the instrument and prosthesis so that each pixel of the radiological image comprises all the data relating to the region to be treated. 
       FIG. 5  schematically illustrates a radiological image of a region to be treated. This image  20  comprises a stent  202 , a guide wire  203  and markers  2030 ,  2031 . These markers  2030 ,  2031  are respectively arranged at the two ends of the stent  202  and allow easy identification of the stent  202  in the image  20 . 
     It is to be noted that the blood vessel in which the stent, the guide wire  203  and the markers  2030 ,  2031  are located, has radiation absorption properties which are similar to those of the surrounding tissues. Thus, it does not generate any contrast in the image and is consequently invisible in the radiological image. 
     The method to process radiological images described herein sets out to omit part of the guide wire  203 . This part of the guide wire  203  is the part located between the markers  2030 ,  2031  i.e. inside the stent  202 . Similarly, it would also be possible to omit any other part, or the entirety of the guide. 
     During an initial step  100 , the radiological image  20  of the region to be treated is acquired, then during a second step  200  the acquired radiological image is broken down into a superimposition of layers, one layer solely comprising that part of the instrument to be isolated from the remainder of the image. 
     More generally, the radiological image processed here can be derived from a sequence of previously acquired and recorded radiological images, or it can be derived from the combining together of several images e.g. thirty images. 
       FIG. 6  illustrates layer  21  with that part of the guide wire  203  arranged between the two markers  2030 ,  2031  isolated from the remainder of the image  20 . 
     Once this layer modelling has been obtained, during a third step  300  the layer  21  comprising the instrument isolated from the remainder of the image is removed from the radiological image  20  of the region to be treated. 
     This leads to obtaining the image  22  of the region to be treated in which the portion of instrument has been removed. 
       FIG. 7  illustrates the image  22  of the region to be treated in which the portion of instrument located between the two markers  2030 ,  2031  is deleted. 
     Therefore visualization of the stent  202  is no longer impaired by the presence of the guide wire  203 . 
     In addition, taking into account that the radiological image of the region to be treated is transparent, with said method no data is lost notably around the region in which part of the guide wire  203  is deleted. 
     The layered modelling of the acquired image is made possible since a layer can be described as per a parametric model. 
     A description is given below of the estimation of the layer solely comprising the portion of guide wire  203  that it is desired to delete. 
     To estimate the layer comprising the part of guide wire  203  to be deleted, a parametric model is used of the layer to be estimated. 
     Said layer (its intensity mapping) is mathematically described by the function ƒ(p) in which p is a parameter vector and ƒ is the function allowing intensity values to be obtained which describe the layer, having knowledge of p. To estimate p from the image, several known estimation techniques can be used such as the Maximum a Posteriori technique—MAP. 
       FIG. 8  illustrates the guide wire of the layer to be estimated. 
     The layer is defined as a parametric image ƒ(p) with p=[c,profile] in which c is the central line of all the pixels defining the guide wire  203  and profile is a function which, at a distance d, associates the value of the guide wire  203  for every pixel at this distance d from the central line c. 
     It is assumed here that the value of the pixels of the guide wire only depends on the distance to the central line c. 
     To determine the curve c describing the guide wire  203 , it is possible to use several known techniques such as the one described in document: Adrian Barbu, Vassilis Athitsos, Bogdan Georgescu, Stefan Böhmn Peter Durlak, Dorin Comaniciu “ Hierarchical Learning of Curves Application to Guidewire Localization in Fluoroscopy”,  2007 IEEE Computer Society Conference on Computer Vision and Pattern recognition (CVPR  2007 ), 18-23 Jun. 2007, Minneapolis, Minn., USA, IEEE Computer Society 2007. 
     Once the central line c of the guide wire  203  has been estimated, it is possible to determine profile(d) as follows. 
     It is to be noted that before performing the estimation properly so-called, it is possible to conduct a step  1000  to pre-process the acquired image so as to eliminate the image component corresponding to its background i.e. to structures of much greater extent than the structures of interest, namely the stent  202  and guide wire  203 . 
     This pre-processing step  1000  consists of estimating the background of the image using low-pass filtering or using a mathematical morphology approach of known type. 
     Once the image background has been estimated, this is subtracted from the acquired image using logarithmic subtraction. This can also be considered as the estimation and subtraction of a layer of the image formed by the background. 
     In the resulting image, the value of the pixels of this image is proportional to the absorbance of the thin structures such as the guide wire  203  and the stent  202 . 
     Consideration is given to the pixels at a distance d from the central line c, and a value is determined which represents all these pixels, such as their mean or median. Therefore, with knowledge of c and profile(d) the values of the pixels p i  are obtained which correspond to the guide wire  203 . 
     The layer containing the portion of guide wire  203  to be deleted is finally subtracted from the image  20  of the region to be treated. It is a logarithmic subtraction. 
     Embodiments of the present invention enable the acquisition of a radiological image of a region to be treated using a prosthesis and an instrument, in which the instrument is omitted from the image, wherein the obtained image is able to be used for diagnosis purposes. Embodiments of the present invention simplify the acquisition using a single image, by describing one of the layers using a parametric model that can be estimated from the acquired image. Embodiments of the present invention enable the practitioner to obtain all necessary information for the procedure, and is easy to implement with known imaging systems.