US20050259784A1 - Methods for spectrally calibrating CT imaging apparatus detectors - Google Patents
Methods for spectrally calibrating CT imaging apparatus detectors Download PDFInfo
- Publication number
- US20050259784A1 US20050259784A1 US10/850,009 US85000904A US2005259784A1 US 20050259784 A1 US20050259784 A1 US 20050259784A1 US 85000904 A US85000904 A US 85000904A US 2005259784 A1 US2005259784 A1 US 2005259784A1
- Authority
- US
- United States
- Prior art keywords
- accordance
- estimated
- dataset
- ideal
- projection
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
- A61B6/582—Calibration
- A61B6/585—Calibration of detector units
Definitions
- This invention relates generally to methods for computed tomography (CT) imaging and more particularly to methods and apparatus for calibrating CT imaging apparatus detectors to produce images having reduced artifacts.
- CT computed tomography
- X-ray CT measures x-ray attenuation through an object with a plurality of detector elements at various angular positions with respect to the object. With appropriate data processing, a cross-sectional image revealing inner structures of the object can be reconstructed. Reconstructed images represent a property of an object in response to x-rays (or other radiation), namely, the x-ray (or other radiation) linear attenuation.
- calibration of an CT imaging system and an appropriate correction must to be applied in the reconstruction process. For example, calibration and correction related detector response to incident x-ray spectrum must be performed.
- X-ray spectra often have wide energy bandwidths, typically from a few tens of keV to more than a hundred keV.
- a spectral calibration can be performed to determine residual spectral response differences among detector elements and to ensure that reconstructed images of water phantoms are uniform and free of artifacts. Inaccurate determination of spectral response differences between detector elements results in rings and bands in images around the iso-center.
- uniform phantoms made of water-like material are used for calibration. By scanning these phantoms, the detector response to an x-ray spectrum attenuated through a given length of water can be determined.
- multiple phantoms of various sizes are often required to cover the attenuation range needed for some types of imaging, including medical imaging, leading to lengthy calibration times.
- Spectral calibration is often performed every few months. Because of the complexity of the calibration procedure, well-trained technicians must usually be employed to perform these calibrations.
- some configurations of the present invention provide a method for calibrating a computed tomographic imaging apparatus having a gantry, a radiation source operable at a plurality of kVp's, and a detector array having a plurality of detector elements.
- the method includes using a system detection function to estimate signals of each detector element that would have been detected through air and through a given thickness of water to determine estimated datasets.
- the estimated datasets are used to determine data pair sets each comprising a normalized water projection value and an ideal projection value for each detector element.
- the method further includes determining and storing a representation of a mapping function of the normalized water projections values to the ideal projection values in a memory of the computed tomographic imaging apparatus as a spectral calibration of the computed tomographic imaging apparatus.
- some configurations of the present invention provide a computed tomographic imaging apparatus having a gantry, a radiation source operable at a plurality of kVp's, and a detector array having a plurality of detector elements.
- the apparatus is configured to use a system detection function to estimate signals of each detector element that would have been detected through air and through a given thickness of water to determine estimated datasets.
- the apparatus is configured to utilize the estimated datasets to determine data pair sets each comprising a normalized water projection value and an ideal projection value for each detector element.
- the apparatus is configured to determine and store a representation of a mapping function of the normalized water projections values to the ideal projection values in a memory as a spectral calibration.
- some configurations of the present invention provide a computed tomographic imaging apparatus having a rotating gantry, a radiation source operable at a plurality of kVp's, and a detector array having a plurality of detector elements.
- the radiation source and detector array are mounted on the rotating gantry and configured to rotate with the gantry.
- the CT imaging apparatus is configured to determine a set of system detection coefficients utilizing air scan measurements at a plurality of kVp's of the radiation source. Utilizing the determined set of system detection coefficients, the CT imaging apparatus is also configured to determine detected signals of each detector element through air and through a given thickness of water. Using the determined detected signals, the CT imaging apparatus is also configured to determine a plurality of data pair sets having a predicted signal and an ideal signal for each detector element and to obtain a spectral correction for the detector elements using the determined data pair sets.
- Configurations of the present invention in principle do not require the scanning of phantoms. Instead, some configurations of the present invention are able to determine x-ray detection efficiency as a function of photon energy for individual detector elements by measuring x-ray projections through air at multiple kVps. With the detection efficiencies thus determined, spectral errors for an x-ray spectrum attenuated through any given size of water-like object can also be determined. Configurations of the present invention thus can be performed in the time required for conventional air calibrations and can be fully automated, making daily spectral calibration feasible. Moreover, some configurations of the present invention enable spectral calibration of arbitrary or dynamically varying beam filters (such as bow tie filters) that can be used to reduce radiation dose to the patient, even through physical calibration is not necessarily performed with a shape of a particular beam filter.
- beam filters such as bow tie filters
- FIG. 1 is a pictorial drawing representative of some configurations of CT imaging apparatus of the present invention.
- FIG. 2 is a functional block diagram representative of the CT imaging apparatus of FIG. 1 .
- FIG. 3 is a flow chart representative of various configurations of calibration procedures of the present invention.
- Example embodiments of methods and apparatus that facilitate calibration of CT imaging apparatus and CT imaging methods and apparatus having reduced image artifacts are described below in detail.
- a technical effect of the methods and apparatus described herein include at least one of calibration of CT imaging apparatus, reduction of image artifacts, and the determination of detector element detection efficiency as a function of photon energy for individual detector elements.
- an x-ray source projects a fan-shaped beam which is collimated to lie within an X-Y plane of a Cartesian coordinate system and generally referred to as an “imaging plane”.
- the x-ray beam passes through an object being imaged, such as a patient.
- the beam after being attenuated by the object, impinges upon an array of radiation detectors.
- the intensity of the attenuated radiation beam received at the detector array is dependent upon the attenuation of an x-ray beam by the object.
- Each detector element of the array produces a separate electrical signal that is a measurement of the beam intensity at the detector location. The intensity measurements from all the detectors are acquired separately to produce a transmission profile.
- the x-ray source and the detector array are rotated with a gantry within the imaging plane and around the object to be imaged such that the angle at which the x-ray beam intersects the object constantly changes.
- a group of x-ray attenuation measurements, i.e., projection data, from the detector array at one gantry angle is referred to as a “view”.
- a “scan” of the object comprises a set of views made at different gantry angles, or view angles, during one revolution of the x-ray source and detector.
- the projection data is processed to construct an image that corresponds to a two-dimensional slice taken through the object.
- One method for reconstructing an image from a set of projection data is referred to in the art as the filtered backprojection technique. This process converts the attenuation measurements from a scan into integers called “CT numbers” or “Hounsfield units” (HU), which are used to control the brightness of a corresponding pixel on a cathode ray tube display.
- CT numbers or “Hounsfield units” (HU)
- a “helical” scan may be performed.
- the patient is moved while the data for the prescribed number of slices is acquired.
- Such a system generates a single helix from a fan beam helical scan.
- the helix mapped out by the fan beam yields projection data from which images in each prescribed slice may be reconstructed.
- Reconstruction algorithms for helical scanning typically use helical weighing algorithms that weight the collected data as a function of view angle and detector channel index. Specifically, prior to a filtered backprojection process, the data is weighted according to a helical weighing factor, which is a function of both the gantry angle and detector angle. The weighted data is then processed to generate CT numbers and to construct an image that corresponds to a two-dimensional slice taken through the object.
- multi-slice CT has been introduced.
- multi-slice CT multiple rows of projection data are acquired simultaneously at any time instant.
- the system When combined with helical scan mode, the system generates a single helix of cone beam projection data. Similar to the single slice helical, weighting scheme, a method can be derived to multiply the weight with the projection data prior to the filtered backprojection algorithm.
- the phrase “reconstructing an image” is not intended to exclude embodiments of the present invention in which data representing an image is generated but a viewable image is not. However, many embodiments generate (or are configured to generate) at least one viewable image.
- a multi-slice scanning imaging system for example, a Computed Tomography (CT) imaging system 10
- CT Computed Tomography
- Gantry 12 has an x-ray tube 14 (also called x-ray source 14 herein) that projects a beam of x-rays 16 toward a detector array 18 on the opposite side of gantry 12 .
- a filter for example, a bowtie beam filter 19 , may be provided. If provided, the bowtie beam filter may be operated either statically or dynamically.
- Detector array 18 is formed by a plurality of detector rows (not shown) including a plurality of detector elements 20 which together sense the projected x-rays that pass through an object, such as a medical patient 22 between array 18 and source 14 .
- a phantom for example, a water-based phantom
- Each detector element 20 produces an electrical signal that represents the intensity of an impinging x-ray beam and hence can be used to estimate the attenuation of the beam as it passes through object or patient 22 .
- gantry 12 and the components mounted therein rotate about a center of rotation 24 .
- multi-slice detector array 18 includes a plurality of parallel detector rows of detector elements 20 such that projection data corresponding to a plurality of quasi-parallel or parallel slices can be acquired simultaneously during a scan.
- Control mechanism 26 includes an x-ray controller 28 that provides power and timing signals to x-ray source 14 and a gantry motor controller 30 that controls the rotational speed and position of components on gantry 12 .
- a data acquisition system (DAS) 32 in control mechanism 26 samples analog data from detector elements 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 image reconstruction. The reconstructed image is applied as an input to a computer 36 , which stores the image in a storage device 38 .
- Image reconstructor 34 can be specialized hardware or computer programs executing on computer 36 .
- 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 (or a suitable display of another type) 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 .
- computer 36 operates a table motor controller 44 , which controls a motorized table 46 to position patient 22 in gantry 12 . Particularly, table 46 moves portions of patient 22 through gantry opening 48 .
- computer 36 includes a device 50 , for example, a floppy disk drive, CD-ROM drive, DVD drive, magnetic optical disk (MOD) device, or any other digital device including a network connecting device such as an Ethernet device for reading instructions and/or data from a computer-readable medium 52 , such as a floppy disk, a CD-ROM, a DVD or another digital source such as a network or the Internet, as well as yet to be developed digital means.
- computer 36 executes instructions stored in firmware (not shown).
- Computer 36 is programmed to perform functions described herein, and as used herein, the term computer is not limited to just those integrated circuits referred to in the art as computers, but broadly refers to computers, processors, microcontrollers, microcomputers, programmable logic controllers, application specific integrated circuits, and other programmable circuits, and these terms are used interchangeably herein.
- x-ray source 14 In many of the example embodiments described herein, it is assumed that an x-ray source 14 is used and that radiation beam 16 is an x-ray beam. However, with the use of suitable detector elements 20 , a different type of radiation source 14 can be substituted for x-ray source 14 . In such configurations, beam 16 can more generally be referred to as a radiation beam 16 .
- Spectral errors that appear in reconstructed CT images, particularly third generation CT images, result from uncontrolled filtration of radiation beam 16 (for example, shadowing of an antiscatter collimator) and differences in signal production efficiency between detector elements as a function of incident photon energy.
- an incident x-ray spectrum attenuated by the same thickness of water will result in unequal signals being produced from different detector elements 20 .
- attenuation of x-rays 16 through water is strongly dependent upon photon energy, and detection efficiency is also a function of photon energy, a difference in detected signal can not be compensated using air measurements.
- rings and bands occur in third generation CT images.
- a detector array 18 in many configurations comprises a plurality of detector modules, i.e., groups of detector elements 20 that comprise a single unit out of a plurality of such units forming a portion of detector array 18 .
- the probability of finding high quality detector modules in all detector array 18 rows decreases and staging time for phantom-based spectral calibration increases.
- Configurations of the present invention thus provide methods and apparatus for rapid, high sensitivity spectral calibration. Detection efficiency is used to determine spectral errors rather than phantom-by-phantom direct spectral error measurements.
- a number of scans are performed at various kVp settings of x-ray tube 14 .
- Detection efficiency for each detector element 20 is determined using changes in detected signal from these scans.
- Spectral errors for various path lengths of water absorber i.e., a water-based phantom used as object 22 ) are determined using knowledge of the initial x-ray spectrum, the materials of a beam filter (not shown in the Figures), and detector element 20 detection efficiencies as a function of photon energy.
- Some configurations of the present invention thus determine detector element 20 detection efficiencies as a function of photon energy.
- Spectral errors originate from several aspects, namely, beam-hardening by water, the uncontrolled beam filter, and usually most importantly, differential efficiency among the plurality of detector elements 20 .
- Beam hardening by water is a physical property that can be calculated precisely. If the incident x-ray spectra are well-controlled, the beam hardening effect can be accurately predicted and removed.
- Spectral errors resulting from the uncontrolled beam filter and differential efficiency among detector elements 20 are produced by imperfections in manufacturing and are thus unpredictable. The combined effect of uncontrolled beam filter and differential efficiency can be combined mathematically into a deviation in detection efficiency.
- Various configurations of the present invention determine the detection efficiency as a function of incident photon energy by comparing relative signal level among several x-ray spectra.
- FIG. 3 is a flowchart 100 illustrating example processes utilized by CT imaging system 10 .
- Many of the steps shown in flowchart 100 can be performed automatically under control of computer 36 under suitable program control (e.g., a firmware program built into computer 36 or provided device 50 ). However, it is not required that the steps be performed in an automatic progression.
- the technical effect of methods and apparatus utilizing processes represented by flowchart 100 is achieved by a user first operating CT imaging system to acquire (i.e., measure) a detected signal dataset ⁇ D(i, kv) ⁇ at 102 , where i is the detector element 20 index, for each kVp setting provided by the CT system, and without any object in the path of x-ray beam 16 except for system beam filter 19 , if present.
- P ideal (i, kv) ⁇ E ⁇ S ( E, kv )* E[ 1 ⁇ exp( ⁇ (scintillator, E )* T L )]*exp( ⁇ b ( E, kv, i ) ⁇ (1)
- the predicted (i.e., estimated) projection values form a dataset ⁇ P(i, kv) ⁇ that takes into account physical properties of detection components of imaging system 10 .
- the dataset of ideal projection coefficients is determined at 106 from the simulation at 104 .
- the “detection components of imaging system 10 ” include detector elements 20 , radiation source 14 , and filter 19 , if present, and other spectral effects introduced by detector array 18 .)
- an air projection dataset ⁇ A(i, kv) ⁇ is estimated at 114 using the values of X n (i).
- Proj(i, kv, L) ⁇ log[ W ( i, kv, L )/ A ( i, kv )].
- the purpose of spectral cal for water is to linearize a measured projection value to water path length.
- a polynomial fit can characterize the mapping.
- the coefficients of the mapping polynomial can be obtained by fitting data pairs (Proj(i, kv, L), Ideal(i, kv, L)) generated at various water path length L.
- Steps 114 , 116 , 118 , and 120 can be iterated with different bowtie filter 19 shapes, if bowtie filter 19 (or another adjustable filter) is present to obtain different mappings for different values of b(E, kv, i).
- the different mappings are especially useful for applying a spectral correction to adaptive or dynamic beam filters that change the beam filter shape before or during patient scanning.
- an ideal spectral effect is modeled by simulation of an x-ray beam spectrum and its interaction with materials such as filters in the beam path and water phantoms. Deviation from the ideal model is determined from the measurements at multiple kVp's. Detector detection efficiency as a function of photon energy and any additional filtration in the beam path is modeled as the absorption coefficients of a polynomial function directly as shown herein, or two distinct materials plus a polynomial, or any other smooth functional form. System detection coefficients X n (i) are determined from measurements of air scans through known beam filters or through air scans without beam filters, at a plurality of kVp's.
- the number of kVp's required in some configurations is equal to or greater than the number of unknown coefficients in the functional form, or, in some other configurations, a least squares method of obtaining coefficients is used.
- the functional form can be further reduced to linear terms of the unknown coefficients, simplifying the process of obtaining the solutions of the unknown coefficients.
- the detected signal of each detector element through the beam filter and through a given thickness of water is determined for a spectrum at a given kVp.
- the spectral errors of any detector element can be determined by comparing the calculated signal, normalized to a calculated air signal, with the negative of the logarithm taken and incorporating the detection functional form, to the ideal signal (or a corrected signal).
- the ideal signal used here is a linear function of the water penetration length.
- a data pair set of predicated signal and ideal signal is built for each detector element.
- This data pair set provides a direct spectral correction for removing beam hardening in water and eliminating detector errors.
- This data pair set can be fitted by a polynomial, providing fast spectral correction.
- the spectral correction can be performed dynamically or statically for any beam filter shape, upon the determination of the system detection coefficients X n (i).
- a spectral correction can be applied to adaptive beam filters that change the beam filter shape before or during patient scanning for radiation dose reduction.
- Configurations of the present invention can be used in combination with other calibration methods.
- spectral calibration using phantoms can be used in combination with configurations of the present invention to achieve improved image quality by accounting for non-spectral related detector efficiency deviation such as x-ray scatter.
- Configurations of the present invention described herein are directed to spectral correction for water or water-like materials.
- configurations of the present invention that determine detection efficiency as a function of photon energy for detector elements can also be used in other applications.
- Such other applications include, but are not limited to, correction of image artifacts originated from spectral changes due to materials other than water or water-like materials, such as bones, contract agents, metal implants, or other materials. Correction factors can be obtained in a manner similar to that described above by including the extra materials in the processing.
- configurations of the present invention provide accurate determination of the spectral errors for individual detector elements across an entire field of view. These spectral errors can be characterized as a high-order polynomial for fast data correction. Besides providing high accuracy, configurations of the present invention do not require calibration phantoms, thereby allowing rapid and simple calibration. The process can be fully automated and performed daily. Configurations of the present invention make possible spectral calibration for adaptive beam filters that reduce radiation dose to the patient. Such an adaptive beam filter can dynamically change filter shape according to patient size, making it difficult and time-consuming to perform spectral calibration on each of the beam filter configurations.
- configurations of the present invention determine detector response to x-rays captured under a given filter, or no beam filter at all, so that the spectral correction for any given beam filter can be dynamically calculated based on the knowledge of the material in the beam filter, and the determination of X n (i).
Abstract
Description
- This invention relates generally to methods for computed tomography (CT) imaging and more particularly to methods and apparatus for calibrating CT imaging apparatus detectors to produce images having reduced artifacts.
- X-ray CT measures x-ray attenuation through an object with a plurality of detector elements at various angular positions with respect to the object. With appropriate data processing, a cross-sectional image revealing inner structures of the object can be reconstructed. Reconstructed images represent a property of an object in response to x-rays (or other radiation), namely, the x-ray (or other radiation) linear attenuation. To reconstruct artifact free images, calibration of an CT imaging system and an appropriate correction must to be applied in the reconstruction process. For example, calibration and correction related detector response to incident x-ray spectrum must be performed. X-ray spectra often have wide energy bandwidths, typically from a few tens of keV to more than a hundred keV. In third generation CT configurations an X-ray tube and detectors rotate together around an object. This configuration makes image quality sensitive to spectral characteristics of detector elements. This sensitivity tends to be greatest for detector elements near iso-center, the axis around which the CT gantry rotates. Thus, effort has to be expended to produce and select detector elements with nearly identical spectral characteristics.
- After detector elements are arranged in a detector array, a spectral calibration can be performed to determine residual spectral response differences among detector elements and to ensure that reconstructed images of water phantoms are uniform and free of artifacts. Inaccurate determination of spectral response differences between detector elements results in rings and bands in images around the iso-center. In at least one known method for spectral calibration, uniform phantoms made of water-like material are used for calibration. By scanning these phantoms, the detector response to an x-ray spectrum attenuated through a given length of water can be determined. However, multiple phantoms of various sizes are often required to cover the attenuation range needed for some types of imaging, including medical imaging, leading to lengthy calibration times. Spectral calibration is often performed every few months. Because of the complexity of the calibration procedure, well-trained technicians must usually be employed to perform these calibrations.
- In one aspect, some configurations of the present invention provide a method for calibrating a computed tomographic imaging apparatus having a gantry, a radiation source operable at a plurality of kVp's, and a detector array having a plurality of detector elements. The method includes using a system detection function to estimate signals of each detector element that would have been detected through air and through a given thickness of water to determine estimated datasets. The estimated datasets are used to determine data pair sets each comprising a normalized water projection value and an ideal projection value for each detector element. The method further includes determining and storing a representation of a mapping function of the normalized water projections values to the ideal projection values in a memory of the computed tomographic imaging apparatus as a spectral calibration of the computed tomographic imaging apparatus.
- In another aspect, some configurations of the present invention provide a computed tomographic imaging apparatus having a gantry, a radiation source operable at a plurality of kVp's, and a detector array having a plurality of detector elements. The apparatus is configured to use a system detection function to estimate signals of each detector element that would have been detected through air and through a given thickness of water to determine estimated datasets. The apparatus is configured to utilize the estimated datasets to determine data pair sets each comprising a normalized water projection value and an ideal projection value for each detector element. Also, the apparatus is configured to determine and store a representation of a mapping function of the normalized water projections values to the ideal projection values in a memory as a spectral calibration.
- In another aspect, some configurations of the present invention provide a computed tomographic imaging apparatus having a rotating gantry, a radiation source operable at a plurality of kVp's, and a detector array having a plurality of detector elements. The radiation source and detector array are mounted on the rotating gantry and configured to rotate with the gantry. The CT imaging apparatus is configured to determine a set of system detection coefficients utilizing air scan measurements at a plurality of kVp's of the radiation source. Utilizing the determined set of system detection coefficients, the CT imaging apparatus is also configured to determine detected signals of each detector element through air and through a given thickness of water. Using the determined detected signals, the CT imaging apparatus is also configured to determine a plurality of data pair sets having a predicted signal and an ideal signal for each detector element and to obtain a spectral correction for the detector elements using the determined data pair sets.
- Configurations of the present invention in principle do not require the scanning of phantoms. Instead, some configurations of the present invention are able to determine x-ray detection efficiency as a function of photon energy for individual detector elements by measuring x-ray projections through air at multiple kVps. With the detection efficiencies thus determined, spectral errors for an x-ray spectrum attenuated through any given size of water-like object can also be determined. Configurations of the present invention thus can be performed in the time required for conventional air calibrations and can be fully automated, making daily spectral calibration feasible. Moreover, some configurations of the present invention enable spectral calibration of arbitrary or dynamically varying beam filters (such as bow tie filters) that can be used to reduce radiation dose to the patient, even through physical calibration is not necessarily performed with a shape of a particular beam filter.
-
FIG. 1 is a pictorial drawing representative of some configurations of CT imaging apparatus of the present invention. -
FIG. 2 is a functional block diagram representative of the CT imaging apparatus ofFIG. 1 . -
FIG. 3 is a flow chart representative of various configurations of calibration procedures of the present invention. - Example embodiments of methods and apparatus that facilitate calibration of CT imaging apparatus and CT imaging methods and apparatus having reduced image artifacts are described below in detail. A technical effect of the methods and apparatus described herein include at least one of calibration of CT imaging apparatus, reduction of image artifacts, and the determination of detector element detection efficiency as a function of photon energy for individual detector elements.
- In some known CT imaging system configurations, an x-ray source projects a fan-shaped beam which is collimated to lie within an X-Y plane of a Cartesian coordinate system and generally referred to as an “imaging plane”. The x-ray beam passes through an object being imaged, such as a patient. The beam, after being attenuated by the object, impinges upon an array of radiation detectors. The intensity of the attenuated radiation beam received at the detector array is dependent upon the attenuation of an x-ray beam by the object. Each detector element of the array produces a separate electrical signal that is a measurement of the beam intensity at the detector location. The intensity measurements from all the detectors are acquired separately to produce a transmission profile.
- In third generation CT systems, the x-ray source and the detector array are rotated with a gantry within the imaging plane and around the object to be imaged such that the angle at which the x-ray beam intersects the object constantly changes. A group of x-ray attenuation measurements, i.e., projection data, from the detector array at one gantry angle is referred to as a “view”. A “scan” of the object comprises a set of views made at different gantry angles, or view angles, during one revolution of the x-ray source and detector.
- In an axial scan, the projection data is processed to construct an image that corresponds to a two-dimensional slice taken through the object. One method for reconstructing an image from a set of projection data is referred to in the art as the filtered backprojection technique. This process converts the attenuation measurements from a scan into integers called “CT numbers” or “Hounsfield units” (HU), which are used to control the brightness of a corresponding pixel on a cathode ray tube display.
- To reduce the total scan time, a “helical” scan may be performed. To perform a “helical” scan, the patient is moved while the data for the prescribed number of slices is acquired. Such a system generates a single helix from a fan beam helical scan. The helix mapped out by the fan beam yields projection data from which images in each prescribed slice may be reconstructed.
- Reconstruction algorithms for helical scanning typically use helical weighing algorithms that weight the collected data as a function of view angle and detector channel index. Specifically, prior to a filtered backprojection process, the data is weighted according to a helical weighing factor, which is a function of both the gantry angle and detector angle. The weighted data is then processed to generate CT numbers and to construct an image that corresponds to a two-dimensional slice taken through the object.
- To further reduce the total acquisition time, multi-slice CT has been introduced. In multi-slice CT, multiple rows of projection data are acquired simultaneously at any time instant. When combined with helical scan mode, the system generates a single helix of cone beam projection data. Similar to the single slice helical, weighting scheme, a method can be derived to multiply the weight with the projection data prior to the filtered backprojection algorithm.
- As used herein, an element or step recited in the singular and proceeded with the word “a” or “an” should be understood as not excluding plural said elements or steps, unless such exclusion is explicitly recited. Furthermore, references to “one embodiment” of the present invention are not intended to be interpreted as excluding the existence of additional embodiments that also incorporate the recited features.
- Also as used herein, the phrase “reconstructing an image” is not intended to exclude embodiments of the present invention in which data representing an image is generated but a viewable image is not. However, many embodiments generate (or are configured to generate) at least one viewable image.
- Referring to
FIGS. 1 and 2 , a multi-slice scanning imaging system, for example, a Computed Tomography (CT)imaging system 10, is shown as including agantry 12 representative of a “third generation” CT imaging system.Gantry 12 has an x-ray tube 14 (also calledx-ray source 14 herein) that projects a beam ofx-rays 16 toward adetector array 18 on the opposite side ofgantry 12. A filter, for example, abowtie beam filter 19, may be provided. If provided, the bowtie beam filter may be operated either statically or dynamically.Detector array 18 is formed by a plurality of detector rows (not shown) including a plurality ofdetector elements 20 which together sense the projected x-rays that pass through an object, such as amedical patient 22 betweenarray 18 andsource 14. For calibration purposes, a phantom (for example, a water-based phantom) is used asobject 22. Eachdetector element 20 produces an electrical signal that represents the intensity of an impinging x-ray beam and hence can be used to estimate the attenuation of the beam as it passes through object orpatient 22. During a scan to acquire x-ray projection data,gantry 12 and the components mounted therein rotate about a center ofrotation 24.FIG. 2 shows only a single row of detector elements 20 (i.e., a detector row). However,multi-slice detector array 18 includes a plurality of parallel detector rows ofdetector elements 20 such that projection data corresponding to a plurality of quasi-parallel or parallel slices can be acquired simultaneously during a scan. - Rotation of components on
gantry 12 and the operation ofx-ray source 14 are governed by acontrol mechanism 26 ofCT system 10.Control mechanism 26 includes anx-ray controller 28 that provides power and timing signals to x-raysource 14 and agantry motor controller 30 that controls the rotational speed and position of components ongantry 12. A data acquisition system (DAS) 32 incontrol mechanism 26 samples analog data fromdetector elements 20 and converts the data to digital signals for subsequent processing. Animage reconstructor 34 receives sampled and digitized x-ray data fromDAS 32 and performs high-speed image reconstruction. The reconstructed image is applied as an input to acomputer 36, which stores the image in astorage device 38.Image reconstructor 34 can be specialized hardware or computer programs executing oncomputer 36. -
Computer 36 also receives commands and scanning parameters from an operator viaconsole 40 that has a keyboard. An associated cathode ray tube display 42 (or a suitable display of another type) allows the operator to observe the reconstructed image and other data fromcomputer 36. The operator supplied commands and parameters are used bycomputer 36 to provide control signals and information toDAS 32,x-ray controller 28, andgantry motor controller 30. In addition,computer 36 operates atable motor controller 44, which controls a motorized table 46 to positionpatient 22 ingantry 12. Particularly, table 46 moves portions ofpatient 22 throughgantry opening 48. - In one embodiment,
computer 36 includes adevice 50, for example, a floppy disk drive, CD-ROM drive, DVD drive, magnetic optical disk (MOD) device, or any other digital device including a network connecting device such as an Ethernet device for reading instructions and/or data from a computer-readable medium 52, such as a floppy disk, a CD-ROM, a DVD or another digital source such as a network or the Internet, as well as yet to be developed digital means. In another embodiment,computer 36 executes instructions stored in firmware (not shown).Computer 36 is programmed to perform functions described herein, and as used herein, the term computer is not limited to just those integrated circuits referred to in the art as computers, but broadly refers to computers, processors, microcontrollers, microcomputers, programmable logic controllers, application specific integrated circuits, and other programmable circuits, and these terms are used interchangeably herein. - Although the specific embodiment mentioned above refers to a third generation CT system, the methods described herein equally apply to fourth generation CT systems (stationary detector—rotating x-ray source) and fifth generation CT systems (stationary detector and x-ray source). Additionally, it is contemplated that the benefits of the invention accrue to imaging modalities other than CT. Thus, a rotating gantry is not present in all configurations of the present invention. Additionally, although the herein described methods and apparatus are described in a medical setting, it is contemplated that the benefits of the invention accrue to non-medical imaging systems such as those systems typically employed in an industrial setting or a transportation setting, such as, for example, but not limited to, a baggage scanning system for an airport or other transportation center.
- In many of the example embodiments described herein, it is assumed that an
x-ray source 14 is used and thatradiation beam 16 is an x-ray beam. However, with the use ofsuitable detector elements 20, a different type ofradiation source 14 can be substituted forx-ray source 14. In such configurations,beam 16 can more generally be referred to as aradiation beam 16. - Spectral errors that appear in reconstructed CT images, particularly third generation CT images, result from uncontrolled filtration of radiation beam 16 (for example, shadowing of an antiscatter collimator) and differences in signal production efficiency between detector elements as a function of incident photon energy. Thus, an incident x-ray spectrum attenuated by the same thickness of water will result in unequal signals being produced from
different detector elements 20. Because attenuation ofx-rays 16 through water is strongly dependent upon photon energy, and detection efficiency is also a function of photon energy, a difference in detected signal can not be compensated using air measurements. Thus, rings and bands occur in third generation CT images. - A
detector array 18 in many configurations comprises a plurality of detector modules, i.e., groups ofdetector elements 20 that comprise a single unit out of a plurality of such units forming a portion ofdetector array 18. As the number ofdetector 18 rows in aCT imaging system 10 increases and the number of operation modes dramatically increases, the probability of finding high quality detector modules in alldetector array 18 rows decreases and staging time for phantom-based spectral calibration increases. Configurations of the present invention thus provide methods and apparatus for rapid, high sensitivity spectral calibration. Detection efficiency is used to determine spectral errors rather than phantom-by-phantom direct spectral error measurements. - In some configurations of the present invention, for a selected beam filter, a number of scans are performed at various kVp settings of
x-ray tube 14. Detection efficiency for eachdetector element 20 is determined using changes in detected signal from these scans. Spectral errors for various path lengths of water absorber (i.e., a water-based phantom used as object 22) are determined using knowledge of the initial x-ray spectrum, the materials of a beam filter (not shown in the Figures), anddetector element 20 detection efficiencies as a function of photon energy. - Some configurations of the present invention thus determine
detector element 20 detection efficiencies as a function of photon energy. Spectral errors originate from several aspects, namely, beam-hardening by water, the uncontrolled beam filter, and usually most importantly, differential efficiency among the plurality ofdetector elements 20. Beam hardening by water is a physical property that can be calculated precisely. If the incident x-ray spectra are well-controlled, the beam hardening effect can be accurately predicted and removed. Spectral errors resulting from the uncontrolled beam filter and differential efficiency amongdetector elements 20 are produced by imperfections in manufacturing and are thus unpredictable. The combined effect of uncontrolled beam filter and differential efficiency can be combined mathematically into a deviation in detection efficiency. Various configurations of the present invention determine the detection efficiency as a function of incident photon energy by comparing relative signal level among several x-ray spectra. -
FIG. 3 is aflowchart 100 illustrating example processes utilized byCT imaging system 10. Many of the steps shown inflowchart 100 can be performed automatically under control ofcomputer 36 under suitable program control (e.g., a firmware program built intocomputer 36 or provided device 50). However, it is not required that the steps be performed in an automatic progression. The technical effect of methods and apparatus utilizing processes represented byflowchart 100 is achieved by a user first operating CT imaging system to acquire (i.e., measure) a detected signal dataset {D(i, kv)} at 102, where i is thedetector element 20 index, for each kVp setting provided by the CT system, and without any object in the path ofx-ray beam 16 except forsystem beam filter 19, if present. In some configurations, there are four kVp settings, so four scans are recorded at 102. In other configurations, a greater or lesser number of scans are recorded. Separately, at 104, ideal projections are simulated. More particularly, in some configurations, the following parameters are defined in the simulation: - S(E, kv) is a spectrum as a function of x-ray kv (kilovolts) and photon energy E;
- b(E, kv, i) is a combined material attenuation of the beam filters for detector i (this function can change with time for a dynamically varying beam filter, and in such cases, is also a function of t);
- μ(A, E) is a linear attenuation coefficients of a material A;
- Xn(i) are a set of unknown coefficients to be solved;
- TL is the scintillator thickness of the detector elements.
- An ideal projection value Pideal(i, kv) is written:
P ideal(i, kv)=ΣE {S(E, kv)*E[1−exp(−μ(scintillator, E)*T L)]*exp(−b(E, kv, i)} (1) - The predicted (i.e., estimated) projection values (i.e., ideal projection coefficients) form a dataset {P(i, kv)} that takes into account physical properties of detection components of
imaging system 10. The dataset of ideal projection coefficients is determined at 106 from the simulation at 104. (The “detection components ofimaging system 10” includedetector elements 20,radiation source 14, and filter 19, if present, and other spectral effects introduced bydetector array 18.) The predicted projection values are written:
P(i, kv)=ΣE {S(E, kv)*E[1−exp(−μ(scintillator, E)*T L)]*exp(−b(E, kv, i)*f(E, i)} (2)
where f(E, i) is the system detection function for photons with energy E. The system detection function f(E, i) can be modeled as a smooth function of E, for example and not by way of limitation, a polynomial with system detection coefficients Xn(i) such as f(E, i)=X0(i)+X1(i)*E+X2(i)*E2+X3(i)*E3. - The acquired dataset {D(i, kv)} (subtracting an electronic offset reading) is normalized at 108, using an expression written Dnorm(i, kv)=D(i, kv)/ΣiD(i, kv), where Dnorm(i, kv) represents the normalized measured data.
- The ideal projection data coefficients are normalized at 110, using an expression written Pnorm(i, kv)=P(i, kv)/ΣiPideal(i, kv).
- In various configurations of the present invention, the various values of system detection coefficients Xn(i) are determined at 112 for each detector element using a relationship written Pnorm(i, kv)=Dnorm(i, kv) by solving a set of linear equations, if system detection function f(E, i) is expressed in a polynomial form as shown herein. Otherwise, if necessary, nonlinear or numerically-derived techniques can be used to determine the system detection function f(E, i).
- In some configurations, an air projection dataset {A(i, kv)} is estimated at 114 using the values of Xn(i). The model A(i, kv) is written:
A(i, kv)=ΣE {S(E, kv)*E[1−exp(−μ(scintillator, E)*T L)]*exp(−b(E, kv, i)*f(E, i)}. - Also, at 116, a projection dataset through water {W(i, kv, L)} for path length L is estimated at 116 using the values of Xn(i):
W(i, kv, L)=ΣE {S(E, kv)*E[1−exp(−μ(scintillator, E)*T L)]*exp(−b(E, kv, i)*f(E, i)*exp(−μ(water, E)*L)}. - Next, the water projection is logarithmically normalized to the air projection at 118 to produce a projection value Proj(i, kv, L) written:
Proj(i, kv, L)=−log[W(i, kv, L)/A(i, kv)]. - An nth order polynomial is fitted to the data pair (normalized water projection, ideal projection)=(Proj(i, kv, L), Ideal(i, kv, L)) to generate the spectral calibration vector at 120. The normalized water projection data Proj(i, kv, L) is mapped to the ideal projection value, which is written as:
Ideal(i, kv, L)=μ*L,
where μ is a constant that is equal to the water linear attenuation at the average x-ray energy for a given kVp, and spectral cal vectors are also generated. The purpose of spectral cal for water is to linearize a measured projection value to water path length. More particularly, a mapping function Proj(i, kv, L)=>Ideal(i, kv, L) is the spectral calibration, which can either be in the form of a look-up table or in some functional form. A polynomial fit can characterize the mapping. The coefficients of the mapping polynomial can be obtained by fitting data pairs (Proj(i, kv, L), Ideal(i, kv, L)) generated at various water path length L. -
Steps different bowtie filter 19 shapes, if bowtie filter 19 (or another adjustable filter) is present to obtain different mappings for different values of b(E, kv, i). The different mappings are especially useful for applying a spectral correction to adaptive or dynamic beam filters that change the beam filter shape before or during patient scanning. - Thus, an ideal spectral effect is modeled by simulation of an x-ray beam spectrum and its interaction with materials such as filters in the beam path and water phantoms. Deviation from the ideal model is determined from the measurements at multiple kVp's. Detector detection efficiency as a function of photon energy and any additional filtration in the beam path is modeled as the absorption coefficients of a polynomial function directly as shown herein, or two distinct materials plus a polynomial, or any other smooth functional form. System detection coefficients Xn(i) are determined from measurements of air scans through known beam filters or through air scans without beam filters, at a plurality of kVp's. The number of kVp's required in some configurations is equal to or greater than the number of unknown coefficients in the functional form, or, in some other configurations, a least squares method of obtaining coefficients is used. The functional form can be further reduced to linear terms of the unknown coefficients, simplifying the process of obtaining the solutions of the unknown coefficients. Upon obtaining the coefficients in the functional form, the detected signal of each detector element through the beam filter and through a given thickness of water is determined for a spectrum at a given kVp. The spectral errors of any detector element can be determined by comparing the calculated signal, normalized to a calculated air signal, with the negative of the logarithm taken and incorporating the detection functional form, to the ideal signal (or a corrected signal). The ideal signal used here is a linear function of the water penetration length. By making such determinations for all penetration lengths of water relevant to imaging (for example, in a medical imaging application), a data pair set of predicated signal and ideal signal is built for each detector element. This data pair set provides a direct spectral correction for removing beam hardening in water and eliminating detector errors. This data pair set can be fitted by a polynomial, providing fast spectral correction. Moreover, the spectral correction can be performed dynamically or statically for any beam filter shape, upon the determination of the system detection coefficients Xn(i). For example and without limitation, a spectral correction can be applied to adaptive beam filters that change the beam filter shape before or during patient scanning for radiation dose reduction.
- Configurations of the present invention can be used in combination with other calibration methods. For example, spectral calibration using phantoms can be used in combination with configurations of the present invention to achieve improved image quality by accounting for non-spectral related detector efficiency deviation such as x-ray scatter.
- Configurations of the present invention described herein are directed to spectral correction for water or water-like materials. However, configurations of the present invention that determine detection efficiency as a function of photon energy for detector elements can also be used in other applications. Such other applications include, but are not limited to, correction of image artifacts originated from spectral changes due to materials other than water or water-like materials, such as bones, contract agents, metal implants, or other materials. Correction factors can be obtained in a manner similar to that described above by including the extra materials in the processing.
- It will thus be appreciated that configurations of the present invention provide accurate determination of the spectral errors for individual detector elements across an entire field of view. These spectral errors can be characterized as a high-order polynomial for fast data correction. Besides providing high accuracy, configurations of the present invention do not require calibration phantoms, thereby allowing rapid and simple calibration. The process can be fully automated and performed daily. Configurations of the present invention make possible spectral calibration for adaptive beam filters that reduce radiation dose to the patient. Such an adaptive beam filter can dynamically change filter shape according to patient size, making it difficult and time-consuming to perform spectral calibration on each of the beam filter configurations. However, configurations of the present invention determine detector response to x-rays captured under a given filter, or no beam filter at all, so that the spectral correction for any given beam filter can be dynamically calculated based on the knowledge of the material in the beam filter, and the determination of Xn(i).
- While the invention has been described in terms of various specific embodiments, those skilled in the art will recognize that the invention can be practiced with modification within the spirit and scope of the claims.
Claims (20)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/850,009 US7086780B2 (en) | 2004-05-20 | 2004-05-20 | Methods for spectrally calibrating CT imaging apparatus detectors |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/850,009 US7086780B2 (en) | 2004-05-20 | 2004-05-20 | Methods for spectrally calibrating CT imaging apparatus detectors |
Publications (2)
Publication Number | Publication Date |
---|---|
US20050259784A1 true US20050259784A1 (en) | 2005-11-24 |
US7086780B2 US7086780B2 (en) | 2006-08-08 |
Family
ID=35375159
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/850,009 Expired - Fee Related US7086780B2 (en) | 2004-05-20 | 2004-05-20 | Methods for spectrally calibrating CT imaging apparatus detectors |
Country Status (1)
Country | Link |
---|---|
US (1) | US7086780B2 (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030199757A1 (en) * | 2002-04-22 | 2003-10-23 | Toth Thomas L. | Method and apparatus of modulating radiation filtering during radiographic imaging |
US20060159223A1 (en) * | 2005-01-14 | 2006-07-20 | General Electric Company | Method and apparatus for correcting for beam hardening in CT images |
US20070172020A1 (en) * | 2006-01-12 | 2007-07-26 | Shuya Nambu | X-ray ct apparatus and method of creating correction data for x-ray ct |
US7551708B2 (en) * | 2007-02-07 | 2009-06-23 | General Electric Company | Method of iterative reconstruction for energy discriminating computed tomography systems |
DE102010042388A1 (en) * | 2010-10-13 | 2012-04-19 | Siemens Aktiengesellschaft | Method for calibrating a CT system with at least one focus-detector combination with a quantum-counting detector |
US8315352B2 (en) * | 2010-09-16 | 2012-11-20 | General Electric Company | System and method of spectral calibration and basis material decomposition for X-ray CT systems |
US9299171B2 (en) | 2011-04-05 | 2016-03-29 | Koninklijke Philips N.V. | Adaptive calibration for tomographic imaging systems |
US20170276620A1 (en) * | 2016-03-23 | 2017-09-28 | Carl Zeiss X-ray Microscopy, Inc. | Method and system for spectral characterization in computed tomography x-ray microscopy system |
US20200330065A1 (en) * | 2019-04-18 | 2020-10-22 | Canon Medical Systems Corporation | Method and apparatus with photon-counting x-ray detectors having for spatially-varying energy bin thresholds |
US20220003881A1 (en) * | 2020-07-06 | 2022-01-06 | Canon Medical Systems Corporation | Correction x-ray detector, x-ray ct apparatus, and detector element determining method |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005019572A1 (en) * | 2005-04-27 | 2006-11-09 | Siemens Ag | Method for evaluating and displaying X-ray projection images and X-ray screening device |
US7801264B2 (en) * | 2007-12-19 | 2010-09-21 | General Electric Company | Method for calibrating a dual -spectral computed tomography (CT) system |
WO2010144094A1 (en) * | 2009-06-12 | 2010-12-16 | Analogic Corporation | Correction for source switching in multi energy scanner |
US8155263B2 (en) * | 2009-07-31 | 2012-04-10 | General Electric Company | Apparatus and method for voltage modulation in X-ray computed tomography |
US7995702B2 (en) * | 2009-08-25 | 2011-08-09 | General Electric Company | System and method of data interpolation in fast kVp switching dual energy CT |
US8611627B2 (en) | 2009-12-23 | 2013-12-17 | General Electric Company | CT spectral calibration |
US9125572B2 (en) | 2012-06-22 | 2015-09-08 | University Of Utah Research Foundation | Grated collimation system for computed tomography |
US10357221B2 (en) * | 2014-08-22 | 2019-07-23 | Samsung Electronics Co., Ltd. | Anatomical imaging system with improved detector calibration process |
US10507005B2 (en) | 2016-09-22 | 2019-12-17 | General Electric Company | Spectral calibration of spectral computed tomography (CT) |
US10433811B2 (en) | 2016-12-23 | 2019-10-08 | General Electric Company | Self-calibrating CT detectors, systems and methods for self-calibration |
WO2020028896A1 (en) | 2018-08-03 | 2020-02-06 | Neurologica Corporation | Ai-based rendered volume auto-correction for fixed and mobile x-ray imaging modalities and other imaging modalities |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5450462A (en) * | 1993-11-19 | 1995-09-12 | General Electric Company | Modulation of x-ray tube current during CT scanning with modulation limit |
US5579359A (en) * | 1995-12-21 | 1996-11-26 | General Electric Company | Methods and apparatus for calibrating detector cell output signals |
US6115487A (en) * | 1998-01-08 | 2000-09-05 | General Electric Company | Correction algorithm for bone-induced spectral artifacts in computed tomograph imaging |
US6266434B1 (en) * | 1998-07-17 | 2001-07-24 | General Electric Company | Methods and apparatus for reducing spectral artifacts in a computed tomograph system |
US6325539B1 (en) * | 1998-12-31 | 2001-12-04 | General Electric Company | Calibration simplification for a computed tomograph system |
US6327329B1 (en) * | 1998-08-25 | 2001-12-04 | General Electric Company | Methods and apparatus for monitoring detector image quality |
US20040022364A1 (en) * | 2002-05-31 | 2004-02-05 | Karl Stierstorfer | Method for determining correction coefficients for detector channels for a computed tomograph |
US20040179646A1 (en) * | 2003-03-14 | 2004-09-16 | Jianying Li | Imaging systems and methods |
US20040218728A1 (en) * | 2003-02-07 | 2004-11-04 | Bjoern Heismann | Calibration of the transformation of spectral X-ray attenuation values in density and atomic number information |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6198791B1 (en) | 1998-08-25 | 2001-03-06 | General Electric Company | Scalable multislice imaging system |
US7636413B2 (en) | 2002-04-16 | 2009-12-22 | General Electric Company | Method and apparatus of multi-energy imaging |
US6775352B2 (en) | 2002-08-16 | 2004-08-10 | Ge Medical Systems Global Technology Company, Llc | Method and system for implementing variable x-ray intensity modulation schemes for imaging systems |
-
2004
- 2004-05-20 US US10/850,009 patent/US7086780B2/en not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5450462A (en) * | 1993-11-19 | 1995-09-12 | General Electric Company | Modulation of x-ray tube current during CT scanning with modulation limit |
US5579359A (en) * | 1995-12-21 | 1996-11-26 | General Electric Company | Methods and apparatus for calibrating detector cell output signals |
US6115487A (en) * | 1998-01-08 | 2000-09-05 | General Electric Company | Correction algorithm for bone-induced spectral artifacts in computed tomograph imaging |
US6266434B1 (en) * | 1998-07-17 | 2001-07-24 | General Electric Company | Methods and apparatus for reducing spectral artifacts in a computed tomograph system |
US6327329B1 (en) * | 1998-08-25 | 2001-12-04 | General Electric Company | Methods and apparatus for monitoring detector image quality |
US6325539B1 (en) * | 1998-12-31 | 2001-12-04 | General Electric Company | Calibration simplification for a computed tomograph system |
US20040022364A1 (en) * | 2002-05-31 | 2004-02-05 | Karl Stierstorfer | Method for determining correction coefficients for detector channels for a computed tomograph |
US20040218728A1 (en) * | 2003-02-07 | 2004-11-04 | Bjoern Heismann | Calibration of the transformation of spectral X-ray attenuation values in density and atomic number information |
US20040179646A1 (en) * | 2003-03-14 | 2004-09-16 | Jianying Li | Imaging systems and methods |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030199757A1 (en) * | 2002-04-22 | 2003-10-23 | Toth Thomas L. | Method and apparatus of modulating radiation filtering during radiographic imaging |
US20060159223A1 (en) * | 2005-01-14 | 2006-07-20 | General Electric Company | Method and apparatus for correcting for beam hardening in CT images |
US7391844B2 (en) * | 2005-01-14 | 2008-06-24 | General Electric Company | Method and apparatus for correcting for beam hardening in CT images |
US20070172020A1 (en) * | 2006-01-12 | 2007-07-26 | Shuya Nambu | X-ray ct apparatus and method of creating correction data for x-ray ct |
US8057098B2 (en) | 2006-01-12 | 2011-11-15 | Kabushiki Kaisha Toshiba | X-ray CT apparatus and method of creating correction data for X-ray CT |
US7551708B2 (en) * | 2007-02-07 | 2009-06-23 | General Electric Company | Method of iterative reconstruction for energy discriminating computed tomography systems |
US8315352B2 (en) * | 2010-09-16 | 2012-11-20 | General Electric Company | System and method of spectral calibration and basis material decomposition for X-ray CT systems |
US20120093282A1 (en) * | 2010-10-13 | 2012-04-19 | Siemens Aktiengesellschaft | Method For Calibrating A CT System With At Least One Focus-Detector Combination With A Quanta-Counting Detector |
DE102010042388A1 (en) * | 2010-10-13 | 2012-04-19 | Siemens Aktiengesellschaft | Method for calibrating a CT system with at least one focus-detector combination with a quantum-counting detector |
US8628241B2 (en) * | 2010-10-13 | 2014-01-14 | Siemens Aktiengesellschaft | Method for calibrating a CT system with at least one focus-detector combination with a quanta-counting detector |
US9299171B2 (en) | 2011-04-05 | 2016-03-29 | Koninklijke Philips N.V. | Adaptive calibration for tomographic imaging systems |
US20170276620A1 (en) * | 2016-03-23 | 2017-09-28 | Carl Zeiss X-ray Microscopy, Inc. | Method and system for spectral characterization in computed tomography x-ray microscopy system |
US10859515B2 (en) * | 2016-03-23 | 2020-12-08 | Carl Zeiss X-ray Microscopy, Inc. | Method and system for spectral characterization in computed tomography x-ray microscopy system |
US20200330065A1 (en) * | 2019-04-18 | 2020-10-22 | Canon Medical Systems Corporation | Method and apparatus with photon-counting x-ray detectors having for spatially-varying energy bin thresholds |
JP2020175169A (en) * | 2019-04-18 | 2020-10-29 | キヤノンメディカルシステムズ株式会社 | Imaging apparatus, method, and program |
US10952698B2 (en) * | 2019-04-18 | 2021-03-23 | Canon Medical Systems Corporation | Method and apparatus with photon-counting x-ray detectors having for spatially-varying energy bin thresholds |
US20220003881A1 (en) * | 2020-07-06 | 2022-01-06 | Canon Medical Systems Corporation | Correction x-ray detector, x-ray ct apparatus, and detector element determining method |
US11874409B2 (en) * | 2020-07-06 | 2024-01-16 | Canon Medical Systems Corporation | Correction X-ray detector, X-ray CT apparatus, and detector element determining method |
Also Published As
Publication number | Publication date |
---|---|
US7086780B2 (en) | 2006-08-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7086780B2 (en) | Methods for spectrally calibrating CT imaging apparatus detectors | |
US7391844B2 (en) | Method and apparatus for correcting for beam hardening in CT images | |
US6987833B2 (en) | Methods and apparatus for identification and imaging of specific materials | |
US7409043B2 (en) | Method and apparatus to control radiation tube focal spot size | |
US7747057B2 (en) | Methods and apparatus for BIS correction | |
US6848827B2 (en) | Method and apparatus for calibrating detector spectral response | |
EP2149284B1 (en) | Methods and systems to facilitate correcting gain fluctuations in image | |
US7801264B2 (en) | Method for calibrating a dual -spectral computed tomography (CT) system | |
US20130202079A1 (en) | System and Method for Controlling Radiation Dose for Radiological Applications | |
US20050180541A1 (en) | Method and apparatus for soft-tissue volume visualization | |
US8611627B2 (en) | CT spectral calibration | |
US7283605B2 (en) | Methods and apparatus for scatter correction | |
JP7042806B2 (en) | Spectral Computed Tomography (CT) Spectral Calibration | |
US10593070B2 (en) | Model-based scatter correction for computed tomography | |
US20180174335A1 (en) | X- Ray CT Device and Sequential Correction Parameter Determination Method | |
US7379527B2 (en) | Methods and apparatus for CT calibration | |
US10098603B2 (en) | Method for estimation and correction of grid pattern due to scatter | |
JP4079488B2 (en) | Method and system for detecting blockage of a reference channel in a computed tomography system | |
JP2022113115A (en) | Beam hardening calibration method, x-ray ct apparatus and beam hardening calibration program | |
US6904120B2 (en) | Method and apparatus for correcting bone induced spectral artifacts | |
CN107019518B (en) | Signal processing method for scatter correction in computed tomography and imaging system | |
JP2004531306A (en) | A method for reducing artifacts in target images | |
JP2022145494A (en) | Image processing device, correction method, and program | |
US20230132514A1 (en) | System and method for controlling errors in computed tomography number | |
US20050201605A1 (en) | Methods and apparatus for CT smoothing to reduce artifacts |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FEPP | Fee payment procedure |
Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
AS | Assignment |
Owner name: GENERAL ELECTRIC COMPANY, NEW YORK Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:WU, XIAOYE;BROMBERG, NEIL BARRY;KOLKER, ALEXANDER;AND OTHERS;REEL/FRAME:017790/0636;SIGNING DATES FROM 20040430 TO 20040512 |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
FPAY | Fee payment |
Year of fee payment: 8 |
|
FEPP | Fee payment procedure |
Free format text: MAINTENANCE FEE REMINDER MAILED (ORIGINAL EVENT CODE: REM.) |
|
LAPS | Lapse for failure to pay maintenance fees |
Free format text: PATENT EXPIRED FOR FAILURE TO PAY MAINTENANCE FEES (ORIGINAL EVENT CODE: EXP.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
STCH | Information on status: patent discontinuation |
Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362 |
|
FP | Lapsed due to failure to pay maintenance fee |
Effective date: 20180808 |