Abstract:
A method of determining disturbing signal frequencies in a signal representation of radiation image originating from the presence of an anti-scatter grid during image acquisition, the signal representation being generated by an image read out and digitizing process, whereby the digitizing of an image of a simulated anti-scatter grid is simulated and the signal resulting from this simulation is evaluated to determine the disturbing frequencies.

Description:
FIELD OF THE INVENTION  
       [0001]     The present invention relates to method of identifying disturbing frequencies originating from the use of an anti scatter grid during radiation image recording in computed radiography.  
       BACKGROUND OF THE INVENTION  
       [0002]     A commonly used technique to reduce the amount of scattered X-rays in computed radiography, digital radiography as well as classical film-based X-ray systems is the use of anti-scatter grids. These grids are lead foil strips, placed apart at a certain distance in a suitable covering.  
         [0003]     There exist different types of anti-scatter grids. In parallel grids, the lead foil strips are parallel, while in honeycomb grids the strips are placed in a honeycomb pattern. The grids are stationary or moving. The use of these grids effectively reduces the radiation scattering but occasionally introduces artifacts such as grid lines into the image.  
         [0004]     In a moving system, the motion of the grids removes the grid lines in the image. However, in some circumstances e.g. short exposure time or malfunctioning of the system, the artifacts remain in the image. With stationary grids, the grid lines are almost always visible.  
         [0005]     If the image is formed digitally or converted afterwards to a digital image representation, Moiré artifacts may appear when displaying the image at a certain scale. These low frequent Moiré artifacts are mostly disturbing and should be eliminated. Before displaying the image, the grid lines, if present in the image, should be removed.  
         [0006]     Current gridline correction techniques, when having detected a grid, suppress only the fundamental frequency in the image sometimes in combination with frequencies in the neighborhood of the harmonic frequencies. The assumption of the harmonic frequencies is not always true. The number of harmonic frequencies to be suppressed is fixed in the prior art methods.  
         [0007]     It is an object of the present invention to provide a method to determine the disturbing frequencies in a signal representation of a radiation image that overcomes the above-mentioned shortcoming of the prior art.  
       SUMMARY OF THE INVENTION  
       [0008]     The above-mentioned aspects are realised by a method as set out in claim  1 .  
         [0009]     Specific features for preferred embodiments of the invention are set out in the dependent claims.  
         [0010]     Further advantages and embodiments of the present invention will become apparent from the following description and drawings.  
         [0011]     The present invention determines which frequencies to suppress in an image signal representing an image originating from an x ray exposure whereby an anti scatter grid has been used in order to remove the annoying Moiré artifacts.  
         [0012]     Having found the fundamental frequency of the anti scatter grid in a signal representation of an x ray image, the digitization of an image containing only an anti scatter grid is simulated. Next, spectrum analysis, e.g. Fourier analysis of this simulated digitized signal is used to determine which frequencies the image might be disturbed by the presence of the anti scatter grid. If the candidate frequencies are also present in the image, they are selected as frequencies to suppress. 
     
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
       [0013]      FIG. 1  shows the expected signal representation of a grid image at the pixel positions when only the fundamental frequency is taken into account (solid line).  
         [0014]     The dash-dotted line is the signal representation of a simulated grid, drawn at a higher resolution than the scanning resolution (pixel positions) of a read out device used for reading out a stored radiation image.  
         [0015]     The dotted line shows the simulated signal after simulation of the digitization process of the digitizer. It is clear that the simulated signal contains more frequencies than only the fundamental frequency.  
         [0016]      FIG. 2  shows the result of a Fourier analysis of the simulated grid signal. If we want to suppress the frequencies for which the log magnitude is greater than 2, the relative frequency F 1 , at about 70% of the Nyquist frequency, is a candidate frequency for suppression. F f  is the fundamental frequency found in the image of the grid. 
     
    
     DETAILED DESCRIPTION OF THE INVENTION  
       [0017]     A specific embodiment of the method of the present invention will be described hereinafter.  
         [0018]     The input of the method is the fundamental frequency or period of the anti-scatter grid in an x-ray image of the grid.  
         [0019]     In the context of the present invention the x-ray image of the grid is obtained by means of a computed radiography system. The grid is irradiated by means of x-rays. The radiation image of the grid is stored in a radiation detector such as a photostimulable phosphor screen. The image is then read out in a read out device (also called digitizer) by scanning the screen with stimulating radiation. Image-wise modulated light emitted by the exposed screen upon stimulation is detected and converted into a signal representation.  
         [0020]     The signal representation is analysed and the frequency at which the signal of the anti-scatter grid has maximum amplitude, hereinafter called the fundamental frequency F f , with corresponding period P f , is detected. This fundamental frequency does not necessarily coincide with the true grid frequency F g , with corresponding period P g .  
         [0021]     If F s  is the sampling frequency of the digitizer and F n  the corresponding Nyquist frequency  
           F   n     =       F   s     2       ,       
 
 one could try to guess the true grid frequency F g . Having found this true grid frequency, the corresponding h th  harmonic frequency to suppress in the image signal is 
 
| F   h   |=|k·F   s   −h·F   g |,   Equation 1 
 
 where k is choosen to satisfy 
 
| F   h   |≦F   n . 
 
         [0022]     The success of this scheme lies entirely in the exact guess of F g  and the selection of enough harmonic frequencies to suppress in the image signal. Because we do not know how many frequencies to suppress, a different approach is preferably chosen.  
         [0023]     Having detected the fundamental frequency F f  in the signal representation of the radiation image of the grid, a grid is simulated.  
         [0024]     For this purpose a signal representation is generated representing a grid with the distance between two grid lamella&#39;s being equal to P f . This grid signal representation g is constructed at a higher resolution than the sampling frequency F s  of the read out apparatus:  
         g   i     =     {               a   ⁢     :     ⁢   0     ≤     i   ⁢           ⁢   mod   ⁢           ⁢     P   g       &lt;   r                 1   ⁢     :     ⁢   r     ≤     i   ⁢           ⁢   mod   ⁢           ⁢     P   g       &lt;     P   g             ⁢     
     ⁢     {           0   ≤   a   &lt;   1               0   ≤   r   &lt;   1                     
 
         [0025]     If the pixel size P s  is for example 100 μm, the grid is constructed at a resolution where one sample i represents 1 μm,  
         [0026]     a is the attenuation of the radiation at the positions of the grid lamella&#39;s and r is the ratio of the lamella width to the distance between two lamella&#39;s.  
         [0027]     This grid signal representation is then averaged over distance P s  to obtain the simulated grid signal representation G:  
         G   j     =       ∑     j   =     iP   s             (     i   +   1     )     ⁢     P   s       -     1   s         ⁢       g   j     .           
 
         [0028]     This averaging simulates the digitization process of the image read out device.  
         [0029]     More complex simulations may be used, such as averaging the grid after convolution with a Gaussian kernel to simulate the flying spot of the light stimulation in some digitizers.  
         [0030]     For most purposes, a simple averaging of the grid at the high resolution is sufficient.  
         [0031]     Fourier analysis of the simulated grid representation G shows at which frequencies the original image signal is disturbed (see  FIG. 2 ).  
         [0032]     Only the frequencies F h , with an amplitude greater than a fixed threshold, are selected.  
         [0033]     In the example of  FIG. 2 , in addition to the fundamental frequency F f  located at about 0.9 times the Nyquist frequency, we also want to suppress the frequency located at about 0.7 times the Nyquist frequency.  
         [0034]     The frequencies to be suppressed can be stored in a look up table with the fundamental frequency as the indexing variable and retrieved when these frequencies are to be suppressed in order to eliminate the disturbing effect in an image due to the presence of an anti-scatter grid at image recording.  
         [0035]     The entries of the look up table can thus be obtained by the steps described higher, alternatively they can be obtained by image analysis, e.g. by applying equation (1).