The present invention relates generally to diagnostic imaging and, more particularly, to a method and apparatus for maximizing image quality for an image of a multi-component object while minimizing the absorbed dose by the object on a per-component basis.
Generally, four key properties define the performance of a computed tomography (CT) scan: spatial resolution, temporal resolution, image noise, and radiation dose. Spatial resolution defines the degree of small object detail in an image and is generally affected by a number of factors including detector aperture, number of acquisition views, focal spot size, object magnification, slice thickness, slice sensitivity profile, helical pitch, reconstruction algorithm, pixel matrix, patient motion, and field-of-view. Temporal resolution defines the length of the temporal interval over which the scan data is acquired for a given slice. Generally, it is desirable to increase temporal resolution (i.e., reduce the length of the temporal interval) as it enables improved imaging of anatomy in motion, such as the heart. Image noise is the random error on the reconstructed image pixel values due to quantum noise or electronic noise, and largely depends on scan geometry and protocol, patient-anatomy, and is location dependent. Radiation dose corresponds to the number of x-rays absorbed by the patient during a scan.
There is an increasing desire to reduce radiation dose to a patient during radiographic data acquisition. However, since quantum noise level is inversely proportional to the square root of the number of x-rays, image quality is directly related to radiation exposure. That is, image quality generally improves as higher radiation doses are used for data acquisition. Over the years, radiation profiles have become more and more optimized. Radiation dose is modulated spatially by the use of bowtie filters, resulting in decreased radiation towards the periphery of the field of view to compensate for the reduced path lengths thereat. Radiation dose is modulated temporally by using tube current modulation, resulting in decreased radiation at view angles and z-position where the path lengths are smaller, for example, lower radiation anterio-posterior relative to laterally, or, for example, lower radiation in the head region and higher radiation in the shoulder region. Finally, the energy profile is optimized for a given application by choosing an optimal tube voltage and hardware filtration.
Since some organs are more sensitive than other organs, it is desirable to limit irradiation to sensitive organs as much as possible, for example, minimizing the absorbed dose to the thyroid, the breasts, the eyes, etc. Sensitive anatomical structures generally comprise only a portion of a given region-of-interest of which an image is to be reconstructed. Thus, if the radiation dose is set to the maximum permitted for the sensitive anatomical structures, the entire image will have poor spatial and contrast resolution. In this regard, the radiation experienced by a patient varies during the course of the scan. This variable radiation profile is typically achieved via x-ray tube current modulation, x-ray tube voltage modulation, x-ray pulse width modulation, x-ray filter modulation, x-ray tube focal spot modulation, or a combination thereof.
In conventional CT scans, the variable radiation dose profile is constructed so as to minimize the variance (noise in image) for a given amount of radiation, or vice-versa. In other words, in conventional CT scans, the radiation profile used to define the scan considers the total radiation, but does not consider the effective dose for the patient. That is, conventionally, the optimal radiation profiles for given acceptable noise variances and the manner for achieving those optimal radiation profiles for the several anatomical structures that comprise a given region-of-interest are not considered.
Therefore, it would be desirable to design an apparatus and method for tailoring a radiation dose profile to optimize the radiation dose on a per component structure basis while maintaining image noise below a noise variance level.
The present invention is directed to a dose optimization process that overcomes the aforementioned drawbacks. The present invention includes a methodology to find a spatial and temporal radiation profile that results in a desirable trade-off between image quality and effective patient dose. The effective dose to an object is minimized by determining a segmented component map for the object, parameterizing tube current/energy level/x-ray filtration/x-ray pulse width as a function of time, determining a corresponding absorbed dose map and variance map, and determining an energy level/tube current profile or curve that results in the lowest effective dose to the object for a given constraint on the noise variance, or vice-versa. Therefore, in accordance with an aspect of the invention, an imaging system is disclosed as having a computer that executes a computer program representing a set of instructions that when executed by the computer causes the computer to determine a component map of an object to be imaged. The object has a plurality of identifiable and imageable components. The computer also determines a relationship between coefficients of a radiation profile and resulting effective dose for the object and also determines a relationship between the coefficients of the radiation profile and a measure of the resulting variance in an image of the object. The computer further determines an irradiating profile that results in one of a minimal effective dose for the object without noise in an image of the object exceeding a desired noise variance, a minimal noise variance for an image of the object for a desired effective dose, or a desired effective dose for the object and a desired noise variance for an image of the object without total dose to the object exceeding a prescribed limit and noise in an image of the object not exceeding a noise limit.
In accordance with another aspect, a radiographic imaging system is presented and includes an x-ray source configured to project x-rays towards a detector according to a certain radiation profile, which establishes number of x-rays projected and energy level of the x-rays projected as a function of time and location, and possibly a finite time interval during which x-rays are produced for each view. The detector is configured to output electrical signals in response to a reception of x-rays. The system further has a computer programmed to acquire an organ map for a subject to be imaged and determine a parameterized dose absorption map for the subject to be imaged and determine a parameterized noise variance map for the subject to be imaged. The computer further determines an irradiation profile that minimizes effective dose for each organ of the organ map and maximizes image quality for an image of the subject.
According to another aspect, a method of dose management for a CT scan is disclosed. The method further includes the step profiling anatomical layout of a patient to be scanned wherein the object has a plurality of anatomical structures. The method also includes the steps of determining a relationship between coefficients of a radiation profile and an absorbed dose for each of the plurality of anatomical structures and determining a relationship between the coefficients of the radiation profile and a noise variance for an image of the patient. The method then determines a radiation profile that results in each anatomical structure receiving a minimal radiation dose without exceeding a noise variance for the image of the patient.
Various other features and advantages of the present invention will be made apparent from the following detailed description and the drawings.
The drawings illustrate one preferred embodiment presently contemplated for carrying out the invention.
In the drawings:
The operating environment of the present invention is described with respect to a four-slice computed tomography (CT) system for imaging of a multi-component object, such as a medical patient. However, it will be appreciated by those skilled in the art that the present invention is equally applicable for use with single-slice or other multi-slice configurations. Moreover, the present invention will be described with respect to the detection and conversion of x-rays. However, one skilled in the art will further appreciate that the present invention is equally applicable for the detection and conversion of other types of radiation. The present invention will be described with respect to a “third generation” CT scanner, but is equally applicable with other CT systems. For example, the invention is also applicable with systems having multiple source spots for increased flexibility in determining an optimal radiation profile by individually steering the different sources.
Referring to
Rotation of gantry 12 and the operation of x-ray source 14 are governed by a control mechanism 26 of CT system 10. Control mechanism 26 includes an x-ray controller 28 that provides power and timing signals to an x-ray source 14 and a gantry motor controller 30 that controls the rotational speed and position of gantry 12. A data acquisition system (DAS) 32 in control mechanism 26 samples analog data from detectors 20 and converts the data to digital signals for subsequent processing. An image reconstructor 34 receives sampled and digitized x-ray data from DAS 32 and performs high speed reconstruction. The reconstructed image is applied as an input to a computer 36 which stores the image in a mass storage device 38.
Computer 36 also receives commands and scanning parameters from an operator via console 40 that has a keyboard. An associated cathode ray tube display 42 allows the operator to observe the reconstructed image and other data from computer 36. The operator supplied commands and parameters are used by computer 36 to provide control signals and information to DAS 32, x-ray controller 28 and gantry motor controller 30. In addition, computer 36 operates a table motor controller 44 which controls a motorized table 46 to position patient 22 and gantry 12. Particularly, table 46 moves portions of patient 22 through a gantry opening 48.
The present invention is directed to a process for determining a dose profile that minimizes the effective dose for a certain image quality or optimizes image quality for a given effective dose. For purposes of this application reference will made to mA/kV modulation which establishes the manner in which the x-ray tube is controlled to produce a desired number of x-rays and an energy level for those x-rays as a function of view angle and position. However, it is contemplated that other factors in addition to the energization of the x-ray tube may help define the radiation dose to an object, such as degree and type of x-ray filtration and the length of time a focal spot of a multi-focal spot x-ray tube is energized. Therefore, reference to mA/kV includes the radiation profile that defines the irradiation experienced by a subject as a result of tube current, tube voltage, x-ray filter filtration, focal spot energization, and the like.
Referring now to
Spatial resolution, temporal resolution, image noise, and radiation dose are key parameters for a CT scan. These key parameters can be related to another in the following expression:
σimg˜1/sqrt(D·FWHM3·ST) (Eqn. 1),
where σimg is the standard deviation of the image noise, and D is the radiation dose, FWHM is the full-width-at-half-maximum of the in-plane image point-spread-function, and ST is the slice thickness. While this is a fundamental relationship, the proportionality constant depends strongly on scanner design and efficiency, on the scan protocol, and on the reconstruction technique. Thus, process 50 described above is designed to optimize the number and energy of x-rays generated as a function of time, location, and energy. Thus, for a given scan geometry, an mA value may be established for each view acquisition. For example, for 1000 views, 360 degree acquisition, a radiation value may be established for views 1, 2, 3 . . . 1000. It is recognized that there are some constraints on establishing the radiation values for each view. For example, the radiation settings for each view will be constrained by a maximum value, mAMAX. A parameterized radiation model is then used to compute a dose and a variance map as a function of any possible radiation profile in order to optimize the radiation profile. Therefore, the radiation profile can be modeled as a function of time by the following expression:
mA(τ)=c1·F1(τ)+c2·F2(τ)+ . . . +cN·FN(τ) (Eqn. 2),
where Fi is a basis function for the mA as a function of time time τ and ci is the weight corresponding to this basis function. One skilled in the art will appreciate that by limiting the radiation profile to a fixed number of basis functions Fi, the computational requirements to determine an optimal radiation profile is less demanding because the number of coefficients ci is typically much smaller than the number of views. For example, by using a basis function that constrains tube modulation to operate along a sine curve and a cosine curve, the number of coefficients is limited to two. One skilled in the art that a multitude of coefficients may be used, but the number may be constrained by the physical limitations of the x-ray tube and/or x-ray filter. That is, a fixed number of different tube current modulations may be permitted by the physics of the x-ray tube and/or x-ray filter and, as such, limit the number of coefficients that are considered for the radiation profile. Equation 2 provides a generalized radiation modulation scheme for an exemplary CT system, such as that shown in
Referring again to
Referring now to
Shown in
Referring now to
As set forth with respect to
Effective Dose=i wi·Di (Eqn. 3),
where Di is the average absorbed dose in component i and wi is the weight that is associated with component i. More dose sensitive components are given a higher weight and x-rays to these components will therefore contribute to a larger increase in effective dose. The sum of the weights is assumed to be one. The effective dose is a single value that is desirably minimized and is determined based on the absorbed dose map and the segmented component map.
As shown in
The effective dose formula together with the noise variance map can then be used to optimize dose and image quality on a per-component, per location basis. That is, image noise σ and effective dose D can be calculated as a function of ci or mA(t). Thus, the optimization process can determine D(ci) and σ(x, ci) for a location x. As a result, a constraint can be defined such that σ(x, ci) must be lower than a predefined limit, σlim, in a certain region x ε R and find the ci that minimizes D(ci). On the other hand, the optimization process can similarly require D(ci) to be lower than Dlim and thus minimize the average σ(x, ci) in a certain region x ε R. For example, the result of the optimization process can be a parmetric formula such as D=Σiαi while the noise calculation at the center of the image results in α=Σi βi·exp(ci·γi), where αi, βi, and γi are calculated constants that depend on object composition and scanner geometry, and ci are the coefficients to be chosen in an optimal fashion to minimize D and/or σ.
As a result of the described optimization process, an effective dose profile can be determined for a given noise variance, or vice versa. In the context of medical imaging, the invention advantageously determines a mA/kV/filtration profile that takes into account the anatomical weightings that differentiate sensitive and non-sensitive organs. Thus, sensitive organs can be imaged with the minimum dose required to provide an image with the desired noise variance. As a result, as shown in the schematic of
It is contemplated that the present invention can be used singly or in combination with other dose reduction tools to not only limit radiation exposure to a scan subject but also advantageously prevent detector saturation for those types of detectors that easily saturate in a CT scan, such as photon counting and energy discriminating detectors. Thus, the invention may be used with active filter control techniques that dynamically adjust the degree and shape of filtration during the course of a scan to tailor radiation to the given scan subject so as to reduce dose to the subject as well as prevent detector saturation by non-attenuated or reduced attenuated x-rays.
While the present invention has been described with respect to a “third generation” CT scanner, it is contemplated that the invention is also applicable with other radiographic systems. For example, the invention is equivalently applicable with CT scanners having a rotatable x-ray source and a stationary ring of detectors. Moreover, the invention is applicable with so-called “cine CT” scanners having a stationary ring of detectors and a tungsten ring to generate an imaging electron beam. Further, the invention is applicable with helical CT scanners as well as scanners having multiple detector arrays and/or multiple x-ray sources.
Therefore, in accordance with an embodiment of the invention, an imaging system is disclosed as having a computer that executes a computer program representing a set of instructions that when executed by the computer causes the computer to determine a component map of an object to be imaged. The object has a plurality of identifiable and imageable components. The computer also determines a relationship between coefficients of a radiation profile and resulting effective dose for the object and also determines a relationship between the coefficients of the radiation profile and a measure of the resulting variance in an image of the object. The computer further determines an irradiating profile that results in one of a minimal effective dose for the object without noise in an image of the object exceeding a desired noise variance, a minimal noise variance for an image of the object for a desired effective dose, or a desired effective dose for the object and a desired noise variance for an image of the object without total dose to the object exceeding a prescribed limit and noise in an image of the object not exceeding a noise limit.
In accordance with another embodiment, a radiographic imaging system is presented and includes an x-ray source configured to project x-rays towards a detector according to a certain radiation profile, which establishes number of x-rays projected and energy level of the x-rays projected as a function of time and location, and possibly a finite time interval during which x-rays are produced for each view. The detector is configured to output electrical signals in response to a reception of x-rays. The system further has a computer programmed to acquire an organ map for a subject to be imaged and determine a parameterized dose absorption map for the subject to be imaged and determine a parameterized noise variance map for the subject to be imaged. The computer further determines an irradiation profile that minimizes effective dose for each organ of the organ map and maximizes image quality for an image of the subject.
According to another embodiment, a method of dose management for a CT scan is disclosed. The method further includes the step profiling anatomical layout of a patient to be scanned wherein the object has a plurality of anatomical structures. The method also includes the steps of determining a relationship between coefficients of a radiation profile and an absorbed dose for each of the plurality of anatomical structures and determining a relationship between the coefficients of the radiation profile and a noise variance for an image of the patient. The method then determines a radiation profile that results in each anatomical structure receiving a minimal radiation dose without exceeding a noise variance for the image of the patient. The present invention has been described in terms of the preferred embodiment, and it is recognized that equivalents, alternatives, and modifications, aside from those expressly stated, are possible and within the scope of the appending claims.