Abstract:
For automatically controlling or adjusting the X-ray dose when scanning an examination subject in a CT apparatus, a series of reference measurements are initially made by obtaining CT images of respective phantoms having different attenuation. The tube current for each phantom that produces an image of the phantom having an acceptable noise level, and thus an acceptable image quality, is stored in a table. In a subsequent CT scan of an examination subject, the tube voltage is automatically adjusted when irradiating a region of the subject having attenuation comparable to the attenuation of one of the phantoms, so as to employ the tube current when irradiating the examination subject that produced the image of acceptable quality for the comparable phantom.

Description:
BACKGROUND OF THE INVENTION  
       [0001]     1. Field of the Invention  
         [0002]     The present invention is directed to a method for automatically setting an X-ray dosage for producing a tomographic image, using an X-ray CT device.  
         [0003]     2. Description of the Prior Art  
         [0004]     In the production of images of an examination subject in a computed tomography system, it is unavoidable that a certain level of noise will be present in the image. In general, the level of noise can be reduced by increasing the tube current, and thereby increasing the X-ray dose, however, obviously this cannot be done indiscriminately because overexposure of the examination subject to X-rays must be avoided.  
         [0005]     Moreover, different portions of a human examination subject exhibit different attenuations, due to respectively different compositions and physical dimensions. Particularly in computed tomography examinations wherein a larger portion of the body is scanned, which may encompass such differently dimensioned body regions, it is undesirable, from the point of view of image quality, to select one and only one tube current/X-ray dose for the entire scan, because the dose may be suitable for scanning one portion of the subject, but may not be suitable for scanning another portion. If a single dose is employed for the entire scan, it must be selected at a compromise value.  
         [0006]     Various investigations have been undertaken to analyze the relationship between the level of noise in a CT image, and the respective values of various parameters that are associated with the production of the CT image. One example of such an investigation is the article “Noise Due To Photon Counting Statistics In Computed X-Ray Tomography,” Chesler et al., Journal of Computed Assisted Tomography, Volume 1, No. 1 (1977), pages 64-74.  
       SUMMARY OF THE INVENTION  
       [0007]     It is an object of the present invention to provide a method for adjusting the X-ray dose in the production of a computed tomography image wherein an overall image with constant image quality is produced without overexposing the examination subject to X-rays.  
         [0008]     The above object is achieved in accordance with the principles of the present invention in a method wherein the tube current is limited automatically at higher levels of attenuation of the examination subject, while still insuring appropriate image quality. A tube load computer, preferably mounted on the gantry of the CT apparatus, employs values obtained in reference measurements for different attenuation values and, dependent on the measured attenuation and the particular organ or body region being examined, the tube current is appropriately limited. The reference measurements can be obtained using a water phantom. On the basis of water phantoms having different attenuation, the interrelation of the desired image noise, tube current, maximum attenuation of the subject, and detector voltage are determined. These values are stored in a table for a particular tube voltage and slice thickness. The table is consulted by the load computer and operating parameters are set for examination of a particular organ or body region by selecting, from the table, the parameters for the water phantom having an attenuation most closely approximating the attenuation of the organ or region under examination.  
     
    
     DESCRIPTION OF THE DRAWINGS  
       [0009]      FIG. 1  is a schematic illustration of a CT apparatus for use in explaining the inventive method.  
         [0010]      FIG. 2  is a flowchart illustrating the sources of noise in a CT image.  
         [0011]      FIG. 3  is a graph showing the detector voltage for various modes of operation of a CT apparatus.  
         [0012]      FIG. 4  is a graph showing the ratio of the actual tube current to a reference tube current as a function of attenuation.  
         [0013]      FIG. 5  is a graph showing noise as a function of the subject attenuation.  
         [0014]      FIG. 6  shows the attenuation-dependent tube current for use in the inventive method.  
         [0015]      FIG. 7  schematically illustrates adjustment of the tube current in accordance with the invention for scanning different body regions. 
     
    
     DESCRIPTION OF THE PREFERRED EMBODIMENTS  
       [0016]     First the theoretical interrelation between measurement system, absorption and noise is summarized and is applied to measurements on the Plus 4. Subsequently, procedures for automatic dosage are discussed.  
         [0017]     The number of quanta N0 emitted by the tube during the measurement of a projection and of the slice fade-in h depends on the tube high voltage kV, the tube current IRöhre, the pre-filtering V and the wedge filter W: 
 
 N   0   =N   0 ( kV, I   Röhre   , h, V, W )  (1) 
 
         [0018]     For a system as shown in  FIG. 1 , the pre-filtering can be assumed to be constant. The wedge filter causes a detector-channel-dependent variation of the quanta N 0 . The high voltage, as well as the pre-filtering and wedge filter, determines the spectral energy distribution of the x-ray spectrum.  
         [0019]     For the functional interrelation of N with the high voltage U, the tube voltage I, and slice thickness h, the following holds:  
             N   =       N   0     *       (     U     U   0       )     2.26     *     (     I     I   0       )     *     (     h     h   0       )               (   2   )             
 
 It can be seen that N˜I, N˜h and N˜U 2.26 . If U=const, h=const, then the number of quanta depends in a linear fashion on the tube current I. 
 
         [0021]     After passing through a homogenous medium having thickness d and absorption μ, the following holds for the number of quanta N: 
 
 N=N   0   ·e   −μd   (3) 
 
 The subject attenuation A is given as 
 
 A=N   0   /N   (4) 
 
 N is the number of quanta that impinge on the detector and produce the signal U Det  in a channel of the detector: 
 
 U   Det   =c·N   (5) 
 
         [0024]     The constant c is a property of the detector that describes the absorption of the quanta and their conversion into an electrical signal (Signal/Quant). C depends on the spectral energy distribution of the x-ray spectrum and on the spectral sensitivity of the detector. Without an absorber, N 0  is to be used instead of N in equation (5).  
         [0025]     The detector signal from equation (5) varies due to the quantum noise. In addition, the electronic noise U E  must be taken into account. Thus, the following holds for the detector signal: 
 
 U   Det   =U   Quanten   +U   E   (6) 
 
         [0026]     The quantum noise is based on the Poisson statistic, so that the following holds: 
 
 U   Quanten   =c *( {overscore (N)}±{square root}{square root over (N)} )= c*N   (7) 
 
 In connection with equation (4), the following holds for the quantum noise U Q  after passage through a subject having attenuation A:  
               U   Q     =       c   *         N   0     A         =       U   0         A   *     N   0                     (   8   )             
 
 According to equation (5), U 0  is the detector signal without absorber. From equation (8), it can be seen that there is an attenuation A 0  at which the quantum noise is equal to the electronic noise: 
 
 A   0 =( U   0   /U   E ) 2   /N   0 =( c/U   E ) 2   *N   0   (9) 
 
 Equation (6) can then be transformed to 
 
 U   Det   =c*N   0   /A+c*{square root}{square root over (N     0     /A)}+   U   E   (10) 
 
         [0030]     The detector signal thus depends on system-specific quantities (pre-filtering, wedge filter, detector constant c) that are equal for all modes of operation of the scanner, as well as on quantities (kV, I Röhre , h) that are selected by the user, and of course also on the attenuation A of the subject of measurement. In addition to the noise of the measurement signal, which enters into the reconstruction via the line integral, the noise from the reconstruction itself must be taken into account as shown in  FIG. 2 .  
         [0031]     For linear interpolation, the image noise σ of a homogenous water disk in the center of the measurement field is calculated as follows:  
             σ   =       500       μ   w     ⁢         a   D   3     ⁢   NEQ           ⁢     q   r               (   11   )             
 
 The relative noise factor q r  takes into account the noise of the reconstruction algorithms. According to the definition, q r =1 for the Shepp-Logan core kernel with linear interpolation. a D  is the scanning raster determined by the detector and related to the center of rotation. For the noise equivalent quanta, NEQ, the following holds:  
             NEQ   =         N   R     *     N   0           a   D     *   A               (   12   )             
 
 N R  is the number of readings per scan (number of projections*reading number, e.g. 1.5 sec: 1056 projections·2). From equations (11) and (12), the relation can be seen between the number of quanta (dosage), subject attenuation and image noise.  
               N   0     =       A       a   D   2     *     N   R         *       (       500       μ   w     *   σ       *     q   r       )     2               (   13   )             
 
         [0034]     The following table shows the results of a measurement of the 20 cm water phantom in the Plus4 (UFC detector) mode (Quality/Noise):  
                                                         Parameter   Value   Derived quantity                                Mode   RED   Projections 1056   Reading number N R =  2       KV   140   q r  = 1       mA   146       TI    1.5       SL    10       Kernel   SP90       FOV   250 0/0       Result   σ 140 kV  = 8.76                  
 
 With equation (13), N 0  can thus be calculated (a D =0.34 mm): N 0 =1.68·10 6  
 
         [0036]     For determining the electronic noise:  
                                                                                                           Measurement                Parameter   Value   Evaluation                    Mode   STATIC       Analyze   Load 00001.RCT           KV   140           Row/Proj. Step:   2       MA   146           Format   μV       TI   1.5 sec           STAT/SD       SL   10 mm           WA (19)       SP90               Column/Channel   lb = 350       FOV                   ub = 420            Storage   FPA   Row Stat --&gt;MEAN            Control   X-ray   OFF           platform   Anode Red   OFF       Result           Average electronic noise U E                     in μV                  
 
         [0037]     From the above values, U E =32.93 μV  
         [0038]     For determining the strength of the detector signal:  
                                                                                             Measurement                Parameter   Value   Evaluation                    Absorber   WITHOUT       Analyze   Load 00001.RCT           Mode   RED           Row/Proj. Step:   2       KV   140           Format   μV       mA   146           STAT/Mean       TI   1.5 sec           WA(18)       SL   10 mm           Column/Channel   lb = 350       SP90                   ub = 420            FOV           Row/STAT --&gt;MEAN       Storage   FPA       Control   X-ray   ON       platform   Anode Red   ON       Result           Signal strength U 0  in μV                  
 
         [0039]     From these values, U 0 =4460736 μV  
         [0040]     From equation (8), using the measurements it is thus possible to determine the quantum noise, and also the detector constant c, for each attenuation if a 12000-fold attenuation is assumed, U Q =41 μV results.  
         [0041]     From equation (5), using the measured U 0  (sec. 4.3) and the value N 0  determined in sec. 4.1, the constant c can be determined as c=2.65 μV/Quant  
         [0042]     According to equation (10), U Det =c*N/A holds for the detector signal without noise portions. N can be calculated according to equation (2). The rotation time TI can unproblemmatically be added to the equation. In 4.5, the constant c has been determined for the Plus 4 UFC detector. In  FIG. 3 , the strength of the detector signal is shown for various subject attenuations (water phantom) and measurement protocols. It is clear that for 1 mm slices, beginning with a subject attenuation of 4000 (corresponding to a water phantom 40 cm in diameter), the electronic noise is greater than the detector signal. The attenuation in the shoulder can reach a value of up to 10,000.  FIG. 3  shows the detector voltage for various modes.  
         [0043]     As can be seen from  FIG. 2  noted earlier, the noise is composed of two portions: the noise from the measurement system (quanta and electronic noise) and the noise from the image reconstruction. The latter depends on subject properties (A, μ w ), system parameters (a D ), recording parameters (N R , N 0 ) and parameters of the image reconstruction (interpolation, kernel [q r ]). Subject properties and system parameters cannot be influenced; the parameters of the image reconstruction are normally fixed by the organ mode used, so that N 0 , i.e. the tube current, remains as a variable quantity that has to be optimized online during the measurement. The matching of the tube current to the maximum subject attenuation during rotation is necessary.  
         [0044]     The determination of the subject attenuation is already part of the dosage modulation and can be taken over. For a single line (or row) detector, this involves: 
    1. Offset correction of the raw data     2. Smoothing of the projections (e.g. sliding [or: variable, flexible] averaging)     3. Norming of the data (U Monitor /U k ), U k =Signal of the k-th detector channel     4. Calculation of the maximum of all detector channels.    
 
         [0049]     The optimal dosage is determined by means of reference measurements. Since Equation (11) was calculated for round, homogenous water phantoms, at first only such phantoms need be considered. On the basis of water phantoms having different attenuation, the interrelation of the desired image noise σ, tube current I Ref , maximum attenuation A of the subject and detector voltage U Det  are determined. These values, U Det (A, σ, I), are stored in a table for a particular tube voltage U kV  and slice thickness h.  
         [0050]     If a is fixedly predetermined, i.e., the kernel is not modified, then it is sufficient to consider the detector signal U Det . From equations ( 6 ) and ( 7 ), the following results:  
               U   Det     =       c   *       N   0     A       +     U   E               (   14   )             
 
 For an acceptable image quality, the electronic noise must be negligible. Dependent on the attenuation A, a minimum tube current I is to be predetermined, so that c*N 0 /A&gt;U E  holds. U det  is then simplified to  
               U   Det     =     c   *       N   0     A               (   15   )             
 
         [0052]     U Det   Ref,σ , N o   Ref ,A Ref  are the parameters for the reference measurement of a phantom having attenuation A, so that for a given dosage N 0  an image noise a and a detector signal U Det  is produced. If a subject having a different attenuation A act , i.e., a different diameter, is now measured, then N 0   act  is to be adapted so that the following holds:  
               U   Det     Ref   ,   σ       =     U   Det   act             (   15   )             
 
 Applying equations (2) and (15), the following results:  
                 A   act       A   Ref       =         (       U   act     k   ⁢           ⁢   V         U   Ref     k   ⁢           ⁢   V         )     2.26     *     (       I   act       I   Ref       )     *     (       h   act       h   Ref       )               (   16   )             
 
 In addition, equation (16) takes into account the tube voltage and the slice thickness. If these are not modified, the new tube current I act  can be calculated:  
               I   act     =         A   act       A   Ref       *     I   Ref               (   17   )             
 
         [0055]      FIG. 4 a  was calculated corresponding to equation (17), whereby a round water phantom having a diameter of 20 cm (A=46.5) was used as a reference, and norming took place to I ref =146 mA (see sec. 4). It is clear that for phantoms having a diameter of greater than 25 cm, corresponding to an attenuation of 122, the required tube current can no longer be set! In the shoulder region, the attenuation is greater than 10000! It is therefore necessary to permit a higher degree of noise as the attenuation increases. In  FIG. 5 , the noise is plotted over the attenuation, if the tube current is held constant (blue curve) and the tube current according to  FIG. 4  is increased.  
         [0056]     Because the maximum tube current I max  is predetermined by the system, it is necessary to ensure that I act &lt;I max  continues to hold. However, in clinical application the outlook is not as negative as it appears according to  FIGS. 4 and 5 . It will not be expected to examine a shoulder with the low degree of noise of an abdominal measurement.  
         [0057]     The tube current can be limited by reducing the gradient of the change of the momentary tube current with increasing attenuation. In  FIG. 6 , this is achieved by dividing the attenuation region into three segments having decreasing gradients. Comparison with  FIG. 4  shows that the increase in the tube current is reduced from a factor of 120 to a factor of 5. The division of the attenuation can ensue in a further region having respectively different gradients, so that there results from equation (17):  
               I   act     =       g   ⁡     (     organ   ,     A   organ       )       ⁢       A   act       A   Ref       *     f   user     *     I   Ref               (   18   )             
 
 g (organ, A organ ) is predetermined dependent on the organ mode and the maximum attenuation, if this is known for example from topogram measurements. 
 
         [0059]     F user  is a factor permitting the user to adjust the reference current I Ref  to his preference in order to achieve a lower degree of noise or a higher savings in dosage. The factor should be adjusted once and then left unchanged.  
         [0060]     As in dosage modulation, in automatic dosage the tube current is also modified online during the measurement. In contrast to modulation, which in the present implementation only reduces the tube current, in automatic dosage it is necessary to increase the current.  
         [0061]     One approach is shown in  FIG. 7 , wherein the automatic control in accordance with the invention is designated “auto.” As was also previously the case, the tube load computer takes into account the nominal currents required by the user for the selected organ program, and interactively carries out an optimization process in order to avoid cooling pauses. If the maximum thermal tube load is not achieved using the prescribed scanning protocol, the load computer simulates possible maximum limits. In this way, the nominal currents of all scans can be increased proportionally, or the magnitude of the increase of the currents takes place in an organ-specific manner. These limits are downloaded for control purposes, and predetermine the maximum range of variation for the tube current.  
         [0062]     The following table describes the scans shown in  FIG. 7 :  
                               TABLE 1                       Currents   Remark   Scan 1   Scan 2   Scan 3                   I nom     Selected by user   I nom   1     I nom   2     I nom   3         I TLC max     Calculated by   I TLC max   1     I TLC max   2     I TLC max   3             TLC       I auto     Tube current of   I auto   1     I auto   2     I auto   3             the automatic   (Rotation)   (Rotation)   (Rotation)           dosage                  
 
         [0063]     An online tube load computer can be located on the gantry, which monitors the tube load during the measurement and limits the momentary tube current, so that all planned scans can be measured with a “reasonable” tube current. It is to be avoided that, for example, two of three spirals are measured with an optimal tube current, while the current must be very strongly reduced for the third.  
         [0064]     In examinations using contrast agents, a scan region is measured twice. The control unit notes the maximum attenuation and the modulation function and applies them in the second scan. If necessary, interpolation can take place in the z direction.  
         [0065]     The maximum attenuation is determined ahead of time from two orthogonal topograms. Using the algorithm that is also used in the control unit for the interrelation between attenuation and tube current, a good estimation of the required tube current can be produced ahead of time.  
         [0066]     Although modifications and changes may be suggested by those skilled in the art, it is the intention of the inventors to embody within the patent warranted hereon all changes and modifications as reasonably and properly come within the scope of their contribution to the art.