# Patent

## Title
Method and apparatus for correcting the contrast density of a radiography image

## Abstract

A method and apparatus for correcting the density of an image in radiography comprising an X-ray source, an image recording placed facing the source, and a support on which an object and surroundings to be X-rayed are intended to be positioned. The method comprises acquiring an image as a mask image; acquiring at least one image as an injected image after injecting a contrast agent into the object to be imaged; and determining corrected images of the mask image and of the injected image from acquisition parameters of the apparatus, from a composition of the imaged object and surroundings and from a composition of the injected contrast agent.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This application claims the benefit of a priority under 35 USC 119 to French Patent Application No. 01 10924 filed Aug. 20, 2001, the entire contents of which are hereby incorporated by reference.
BACKGROUND OF THE INVENTION
The invention relates to the processing of the contrast density of an image in a radiography device. More particularly, the invention is directed to interpretation or reducing the appearance of artifacts for image reconstruction.
In numerous X-ray imaging applications where the object of which images are being taken has a very low contrast, a contrast agent is injected for the purpose of improving the contrast of the object. For example, a contrast agent of this sort is used for cardiology examinations or for mammographs or for DSA (Digital Subtraction Angiography) examinations. In many cases, the visibility of the object is further improved by acquiring an image called a mask image before injecting the contrast agent. This mask image is then subtracted from the image made after the contrast agent is injected in order to obtain a final image that will be used. In this case, a logarithmic subtraction is carried out. A logarithm of the mask image and the logarithm of the image made with the contrast medium. The logarithm of the mask image is subtracted from the logarithm of the image made with the contrast medium. The resulting image gives a good idea of the distribution of the concentration of the contrast agent but has drawbacks. Because of the beam hardening of the X-ray radiation, the contrast of an image made with a given density of contrast agents depends on the amount of tissue that is superimposed. In the case of DSA examinations, which require high concentrations of contrast agents, the beam hardening of the X-ray radiation may also be due to the contrast agent itself. This drawback is problematic when trying to interpret the images resulting from the subtraction and generally leads to artifacts when using the images for image reconstruction, in particular, three-dimensional construction.
BRIEF DESCRIPTION OF THE INVENTION
An embodiment of the invention provides a method and apparatus for correcting the contrast density of images enabling a good interpretation or reducing the appearance of artifacts. In particular, an embodiment of the present invention can be applied during subsequent calculations for three-dimensional image reconstruction.
An embodiment of the invention is directed to a method 20 (see FIG. 2) for correcting the density of an image in a radiography apparatus comprising means for providing a radiation source, such as an X-ray source; means for taking or recording or acquiring an image placed facing the source; and means for support on which an object and surroundings to be imaged are intended to be positioned. In an embodiment of the invention, the method comprises the steps of:
acquiring 22 (see FIG. 2) an image which image, will be further indicated as a mask image;
acquiring 24 (see FIG. 2) at least one image, which will further indicated as an injected image, after injecting a contrast agent into the object to be imaged; and
determining 26 (see FIG. 2) corrected images of the mask image and of the injected image from acquisition parameters of the apparatus, from the composition of the surroundings of the imaged object and from the composition of the injected contrast agent.
In an embodiment of the invention, the method for correcting the density of an image has at least one of the following characteristics:
the composition of the surroundings of the object is estimated by means of a typical composition;
during the step of determining corrected images, a radiation spectrum is determined for each of the images having to be corrected;
the radiation spectrum is determined for each of the pixels making up the image;
when determining corrected images, gray level values are determined for each of the images to be corrected;
when determining corrected images, a thickness of tissue dtissue traversed by the radiation and a thickness of contrast medium dcontrast traversed by the radiation are determined;
dtissue and dcontrast are approximated by polynomial functions;
all the possible values of dtissue and of dcontrast are calculated and the value of dcontrast-dtissue is stored in a two-dimensional table for each of the possible pairs (mask image, injected image), the corrected images of the mask image and of the injected image being determined using the corresponding value of dconstrast-dtissue read from the two-dimensional table when acquiring an image, the acquisition parameters of the apparatus used for this image are combined with the image—during the acquisition of an image, the acquisition parameters of the device used to acquire this image are combined with the image.
In an embodiment of the invention a radiography apparatus comprising means for providing an radiation source, such as an X-ray source; means for recording an image placed facing the source; and means for support on which an object and surroundings to be X-rayed are intended to be positioned is provided, and means for implementing the method having one at least of the aforementioned characteristics.



BRIEF DESCRIPTION OF THE DRAWING
The invention will become apparent from the following description with reference to the appended drawings in which:
FIG. 1 depicts a schematic view of a gantry of an apparatus for acquiring X-ray images; and
FIG. 2 depicts in flowchart form an exemplary method for correcting contrast density of an image in accordance with embodiments of the invention.



DETAILED DESCRIPTION OF THE INVENTION
With reference to the figure, a radiography apparatus 1, such as may be used with X-rays, comprises means 2 for recording or taking or acquiring radiographic images and means 3 for emission, such as an X-ray source. In this case, the means 2 for taking images may be a digital camera. The source 3 and the digital camera 2 are attached at each end of a carrying arm acting as a gantry, in this case in the form of a semicircle. The semicircular arm 7 is connected by sliding to a second arm 8. The second arm 8 itself is connected by sliding and by rotation to the stand 9 of the apparatus. The arm 8 is mainly capable of rotational movements 6 around its own axis. The semicircular arm 7 itself is capable of sliding with respect to the arm 8, so that the semicircular arm 7 makes a rotational movement 5 with respect to the center of the semicircle forming the arm 7.
In use, an object 4 to be imaged, such as a body of a patient is positioned between the source 3 and the digital camera 2, such that the object to be imaged is in a field 10 of the apparatus.
In order to be able to provide images of optimum quality of imaged objects having a low contrast during radiography, it is common practice to inject a contrast medium into the object to be imaged. This contrast agent spreads within the object to be imaged thereby making it possible to increase this contrast on the images taken by the apparatus. In order to increase the quality of the result, the operator takes an image or images before having injected the contrast agent. This image will hereinafter be indicated as a mask image. The image or images taken after injecting the contrast agent will hereinafter be indicated as an injected images.
In order to optimise the quality of the images which will be subsequently used, for example, for diagnosis by viewing or for a three-dimensional image reconstruction of the imaged object or to improve the subtraction from the corrected injected images of the corrected mask image, which makes it possible to completely isolate the imaged object from its surroundings, the apparatus will determine the radiation spectrum, i.e., of the X-rays, at each point of the image taken. To do this, the apparatus will use the acquisition parameters that it has used to take the images. Furthermore, the apparatus is acquainted with constitution or the composition of the object and with its surroundings that has been imaged. Preferably, the composition of the imaged object and/or its surroundings is generally estimated and approximated by a simple typical composition. For example, in the case of an examination targeting the abdominal region of a human body, the typical composition corresponding to water can be selected as an estimate of the surroundings of the object to be imaged in an abdominal region. Within the scope of a neurological examination targeting the brain, a typical composition of the surroundings of the object to be imaged will be bony tissue. This is because, within the scope of this examination, the majority of the attenuation comes from the bone constituting the skull. These various typical compositions can be pre-recorded directly within the apparatus allowing the rapid selection by the operator of a typical composition depending on the region of the body of a patient of which images are being taken.
Moreover, within the context of calculating the radiation spectrum for X-rays, the injected images, the apparatus uses the composition of the contrast agent that has been injected within the object to be imaged for the calculation. Generally the contrast agent used in such examinations is iodine, but the apparatus can have available a multitude of contrast agents which can be used within the context of radiography. The choice of contrast agent, and therefore its composition, is made available as above for a chosen typical composition, that is to say according to the typical composition used by the operator to carry out the injection.
Once determining the radiation spectrum of each of the images made or of each of the points or pixels of the images made (this depends on the degree of accuracy which is required during later use of the corrected images), the apparatus calculates the effective attenuation coefficient of the contrast agent or of the imaged object. The attenuation coefficient may vary throughout the image because, on the one hand, of X-ray beam hardening and because, on the other hand, of the amount of tissue which surrounds the imaged object which the X-rays have to traverse to reach the means for taking the image.
These assumptions permit calculation of the X-ray spectrum of each of the images produced, so as to write the equations relating the X-ray spectrum to the gray level values for each the images produced:




   
  
    {
    
      
        
          
            
            ⁢
            
              
                I
                mask
              
              =
              
                
                  ∫
                  0
                  
                    E
                    max
                  
                
                ⁢
                
                  
                    η
                    ⁡
                    
                      (
                      E
                      )
                    
                  
                  ·
                  
                    N
                    ⁡
                    
                      (
                      E
                      )
                    
                  
                  ·
                  
                    e
                    
                      
                        -
                        
                          
                            μ
                            tissue
                          
                          ⁡
                          
                            (
                            E
                            )
                          
                        
                      
                      ·
                      
                        d
                        tissue
                      
                    
                  
                  ·
                  
                    
                  
                  ⁢
                  
                    d
                    E
                  
                
              
            
          
        
      
      
        
          
            
              I
              injected
            
            =
            
              
                ∫
                0
                
                  E
                  max
                
              
              ⁢
              
                
                  η
                  ⁡
                  
                    (
                    E
                    )
                  
                
                ·
                
                  N
                  ⁡
                  
                    (
                    E
                    )
                  
                
                ·
                
                  e
                  
                    
                      -
                      
                        
                          μ
                          tissue
                        
                        ⁡
                        
                          (
                          E
                          )
                        
                      
                    
                    ·
                    
                      d
                      tissue
                    
                  
                
                ·
                
                  e
                  
                    
                      -
                      
                        
                          μ
                          contrast
                        
                        ⁡
                        
                          (
                          E
                          )
                        
                      
                    
                    ·
                    
                      d
                      contrast
                    
                  
                
                ·
                
                  
                
                ⁢
                
                  d
                  E
                
              
            
          
        
      
    
  




In which:

Imask represents the gray level values in the mask image;
Iinjected represents the gray level values of the injected image in question;
η(E) represents the conversion factor of the digital camera sensor;
(E) represents the number of photons having energy E (this actually represents the X-ray spectrum);
dtissue represents the thickness of the tissue traversed by the incident X-ray beam;
dconstrast represents the thickness of the contrast medium traversed by the incident X-ray beam;
μtissue represents the estimated typical composition of the surrounding tissues and/or of the X-rayed object; and
μcontrast represents the composition of the contrast agent used to produce the injected images.
On the basis of these two equations, the two unknowns dtissue and dcontrast, can easily be calculated by solving the system formed by the two aforementioned equations. Knowledge of dtissue and dconstrast provides the corrected images of the mask image and of the injected image.
Knowledge of dtissue and dcontrast provides the following advantages:
the same amount of contrast agent creates the same contrast on an image, independently of the structure and of the amount of these structures surrounding the object to be imaged;
when the corrected mask image is subtracted from the corrected injected image, the contrast of the vessels is not degraded when it is behind the bone, for example;
when carrying out the three-dimensional image reconstruction from the corrected images, this three-dimensional image reconstruction has considerably fewer artifacts due to the X-ray radiation beam hardening.
Solving the aforementioned equation system provides a first function to obtain dtissue and a second function to obtain dcontrast. These two functions are implemented within the apparatus to calculate the corrected images. To do this, it is possible to use a polynomial approximation which is simple to program.
Another implementation or embodiment of the method comprises calculating a two-dimensional table containing all the possible combinations of Imask and Iinjected using the X-ray spectrum calculated on the basis of the acquisition parameters of a given image. For a given spectrum, the two-dimensional table can be constructed as follows:
for a given value of Imask (from all the possible values which are finite in number), dtissue is calculated as described above, by using the first equation of the above system of equations;
for a given value of Iinjected (among all the possible values which are finite in number), dconstrast is calculated as described above, by using the second equation of the above system of equations, and this for each of the values of dtissue calculated at the previous step; and
the value dconstrast-dtissue corresponding to the images Imask and Iinjected in question are stored in the two-dimensional table, at the location (Imask, Iinjected).
In this case, the apparatus that has just taken an image, in order to determine the corrected image thereof, would be placed at the intersection (Imask, Iinjected) of a two-dimensional table corresponding to the Imask of the mask effectively taken and to the Iinjected of the injected image effectively taken.
By using this method for correcting images in an radiography apparatus, the corrected images obtained are independent of the composition and of the amount of tissue surrounding the imaged object. Moreover, the corrected images obtained are quantitative and can be compared firstly between different examinations and secondly between different viewing angles are also independent of the acquisition parameters used between each of the views or of the examinations. During subtraction of the mask image and of the injected image or of the three-dimensional reconstruction on the basis of the corrected images, the result has considerably fewer artifacts introducing errors that can be damaging to this type of calculation.
Various modifications in the structure and/or the steps and/or the function or equivalents thereof of the disclosed embodiments may be made by one skilled in the art without departing from the scope and extent of the invention as recited in the claims.

## Claims

1. A method for correcting contrast density of an image comprising:
acquiring an image as a mask image;
acquiring at least one image as an injected image, after injecting a contrast agent into an object to be imaged; and
determining corrected images of the mask image and of the injected image from acquisition parameters of an apparatus that produces the image, from a composition of the imaged object and surroundings and from a composition of the injected contrast agent.

2. The method according to claim 1 wherein the composition of the surroundings and of the object is estimated and approximated by means of a typical composition.

3. The method according to claim 1 wherein when determining corrected images, a radiation spectrum is determined for each of the acquired images.

4. The method according to claim 2 wherein when determining corrected images, a radiation spectrum is determined for each of the acquired images.

5. The method according to claim 1 wherein a radiation spectrum is determined for each of the pixels of each of the acquired images.

6. The method according to claim 3 wherein the radiation spectrum is determined for each of the pixels of each of the acquired images.

7. The method according to claim 4 wherein the radiation spectrum is determined for each of the pixels of each of the acquired images.

8. The method according to claim 1 wherein when determining corrected images, gray level values are determined for each of the acquired images.

9. The method according to claim 2 wherein when determining corrected images, gray level values are determined for each of the acquired images.

10. The method according to claim 3 wherein when determining corrected images, gray level values are determined for each of the acquired images.

11. The method according to claim 5 wherein when determining corrected images, gray level values are determined for each of the acquired images.

12. The method according claim 1 wherein when determining corrected images, a thickness of tissue dtissue traversed by radiation and a thickness of contrast medium dcontrast traversed by radiation are determined.

13. The method according claim 2 wherein when determining corrected images, a thickness of tissue dtissue traversed by radiation and a thickness of contrast medium dcontrast traversed by radiation are determined.

14. The method according claim 3 wherein when determining corrected images, a thickness of tissue dtissue traversed by radiation and a thickness of contrast medium dcontrast traversed by radiation are determined.

15. The method according claim 5 wherein when determining corrected images, a thickness of tissue dtissue traversed by radiation and a thickness of contrast medium dcontrast traversed by radiation are determined.

16. The method according claim 8 wherein when determining corrected images, a thickness of tissue dtissue traversed by radiation and a thickness of contrast medium dcontrast traversed by radiation are determined.

17. The method according to claim 12 wherein dtissue and dcontrast are approximated by polynomial functions.

18. The method according to claim 13 wherein dtissue and dcontrast are approximated by polynomial functions.

19. The method according to claim 14 wherein dtissue and dcontrast are approximated by polynomial functions.

20. The method according to claim 15 wherein in that dtissue and dcontrast are approximated by polynomial functions.

21. The method according to claim 16 wherein dtissue and dcontrast are approximated by polynomial functions.

22. The method according to claim 12 wherein all the possible values of dtissue and of dcontrast are calculated and the value of dtissue-dcontrast is stored in a two-dimensional table for each of the possible pairs (mask image, injected image), the corrected images of the mask image and of the injected image being determined using the corresponding value of dcontrast-dtissue read from the two-dimensional table.

23. The method according to claim 13 wherein all the possible values of dtissue and of dcontrast are calculated and the value of dcontrast-dtissue is stored in a two-dimensional table for each of the possible pairs (mask image, injected image), the corrected images of the mask image and of the injected image being determined using the corresponding value of dcontrast-dtissue read from the two-dimensional table.

24. The method according to claim 14 wherein all the possible values of dtissue and of dcontrast are calculated and the value of dcontrast-dtissue is stored in a two-dimensional table for each of the possible pairs (mask image, injected image), the corrected images of the mask image and of the injected image being determined using the corresponding value of dcontrast-dtissue read from the two-dimensional table.

25. The method according to claim 15 wherein all the possible values of dtissue and of dcontrast are calculated and the value of dcontrast-dtissue is stored in a two-dimensional table for each of the possible pairs (mask image, injected image), the corrected images of the mask image and of the injected image being determined using the corresponding value of dcontrast-dtissue read from the two-dimensional table.

26. The method according to claim 16 wherein all the possible values of dtissue and of dcontrast are calculated and the value of dcontrast-dtissue is stored in a two-dimensional table for each of the possible pairs (mask image, injected image), the corrected images of the mask image and of the injected image being determined using the corresponding value of dcontrast-dtissue read from the two-dimensional table.

27. The method according to claim 1 wherein when acquiring an image, the acquisition parameters of the apparatus used to acquire the image are associated with this image.

28. The method according to claim 2 wherein when acquiring an image, the acquisition parameters of the apparatus used to acquire the image are associated with this image.

29. The method according to claim 3 wherein when acquiring an image, the acquisition parameters of the apparatus used to acquire the image are associated with this image.

30. The method according to claim 5 wherein when acquiring an image, the acquisition parameters of the apparatus used to acquire the image are associated with this image.

31. The method according to claim 8 wherein when acquiring an image, the acquisition parameters of the apparatus used to acquire the image are associated with this image.

32. The method according to claim 17 wherein when acquiring an image, the acquisition parameters of the apparatus used to acquire the image are associated with this image.

33. The method according to claim 22 wherein when acquiring an image, the acquisition parameters of the apparatus used to acquire the image are associated with this image.

34. The method according to claim 1 wherein when acquiring an image, the acquisition parameters of the apparatus used to acquire the image are combined with this image.

35. The method according to claim 2 wherein the typical composition corresponds to water or of bony tissue.

36. The method according to claim 2 wherein the typical composition is recorded in the apparatus.

37. The method according to claim 1 wherein the contrast agent is multitude of contrast agents.

38. The method according to claim 1 wherein the method solves the equation




   
  
    {
    
      
        
          
            
            ⁢
            
              
                I
                mask
              
              =
              
                
                  ∫
                  0
                  
                    E
                    max
                  
                
                ⁢
                
                  
                    η
                    ⁡
                    
                      (
                      E
                      )
                    
                  
                  ·
                  
                    N
                    ⁡
                    
                      (
                      E
                      )
                    
                  
                  ·
                  
                    e
                    
                      
                        -
                        
                          
                            μ
                            tissue
                          
                          ⁡
                          
                            (
                            E
                            )
                          
                        
                      
                      ·
                      
                        d
                        tissue
                      
                    
                  
                  ·
                  
                    
                  
                  ⁢
                  
                    d
                    E
                  
                
              
            
          
        
      
      
        
          
            
              I
              injected
            
            =
            
              
                ∫
                0
                
                  E
                  max
                
              
              ⁢
              
                
                  η
                  ⁡
                  
                    (
                    E
                    )
                  
                
                ·
                
                  N
                  ⁡
                  
                    (
                    E
                    )
                  
                
                ·
                
                  e
                  
                    
                      -
                      
                        
                          μ
                          tissue
                        
                        ⁡
                        
                          (
                          E
                          )
                        
                      
                    
                    ·
                    
                      d
                      tissue
                    
                  
                
                ·
                
                  e
                  
                    
                      -
                      
                        
                          μ
                          contrast
                        
                        ⁡
                        
                          (
                          E
                          )
                        
                      
                    
                    ·
                    
                      d
                      contrast
                    
                  
                
                ·
                
                  
                
                ⁢
                
                  d
                  E
                
              
            
          
        
      
    
  




in which:
Imask represents gray level values in the mask image;
Iinjected represents gray level values of the injected image in question;
η(E) represents a conversion factor of a means for taking or recording or acquiring the image;
N(E) represents a number of photons having energy E;
dtissue represents thickness of the tissue traversed by the incident radiation beam;
dcontrast represents thickness of the contrast medium traversed by the incident radiation beam;
μtissue represents estimated typical composition of the surrounding tissues and/or of the imaged object; and
μcontrast represents composition of the contrast agent used to produce the injected images.

39. A radiography apparatus comprising:
means for providing a radiation source;
means for taking or recording or acquiring an image placed facing the means for providing a radiation source;
means for providing a support on which an object and surroundings to be imaged are intended to be positioned;
means for implementing a method for correcting contrast density of an image comprising:
acquiring an image as a mask image;
acquiring at least one image as an injected image, after injecting a contrast agent into an object to be imaged; and
determining corrected images of the mask image and of the injected image from acquisition parameters of an apparatus that produces the image, from a composition of the imaged object and surroundings and from a composition of the injected contrast agent.

40. The apparatus according to claim 39 wherein the means for implementing solves the equation:




   
  
    {
    
      
        
          
            
            ⁢
            
              
                I
                mask
              
              =
              
                
                  ∫
                  0
                  
                    E
                    max
                  
                
                ⁢
                
                  
                    η
                    ⁡
                    
                      (
                      E
                      )
                    
                  
                  ·
                  
                    N
                    ⁡
                    
                      (
                      E
                      )
                    
                  
                  ·
                  
                    e
                    
                      
                        -
                        
                          
                            μ
                            tissue
                          
                          ⁡
                          
                            (
                            E
                            )
                          
                        
                      
                      ·
                      
                        d
                        tissue
                      
                    
                  
                  ·
                  
                    
                  
                  ⁢
                  
                    d
                    E
                  
                
              
            
          
        
      
      
        
          
            
              I
              injected
            
            =
            
              
                ∫
                0
                
                  E
                  max
                
              
              ⁢
              
                
                  η
                  ⁡
                  
                    (
                    E
                    )
                  
                
                ·
                
                  N
                  ⁡
                  
                    (
                    E
                    )
                  
                
                ·
                
                  e
                  
                    
                      -
                      
                        
                          μ
                          tissue
                        
                        ⁡
                        
                          (
                          E
                          )
                        
                      
                    
                    ·
                    
                      d
                      tissue
                    
                  
                
                ·
                
                  e
                  
                    
                      -
                      
                        
                          μ
                          contrast
                        
                        ⁡
                        
                          (
                          E
                          )
                        
                      
                    
                    ·
                    
                      d
                      contrast
                    
                  
                
                ·
                
                  
                
                ⁢
                
                  d
                  E
                
              
            
          
        
      
    
  




in which:
Imask represents gray level values in the mask image;
Iinjected represents gray level values of the injected image in question;
η(E) represents a conversion factor of a means for taking or recording or acquiring the image;
N(E) represents a number of photons having energy E;
dtissue represents thickness of the tissue traversed by the incident radiation beam;
dcontrast represents thickness of the contrast medium traversed by the incident radiation beam;
μtissue represents estimated typical composition of the surrounding tissues and/or of the imaged object; and
μcontrast represents composition of the contrast agent used to produce the injected images.