US20020015476A1 - X-ray computed tomography apparatus - Google Patents
X-ray computed tomography apparatus Download PDFInfo
- Publication number
- US20020015476A1 US20020015476A1 US09/912,137 US91213701A US2002015476A1 US 20020015476 A1 US20020015476 A1 US 20020015476A1 US 91213701 A US91213701 A US 91213701A US 2002015476 A1 US2002015476 A1 US 2002015476A1
- Authority
- US
- United States
- Prior art keywords
- attenuation
- value
- overall
- partial
- image
- 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
- 238000002591 computed tomography Methods 0.000 title claims abstract description 18
- 238000012937 correction Methods 0.000 claims abstract description 57
- 239000012925 reference material Substances 0.000 claims abstract description 33
- 239000000126 substance Substances 0.000 claims abstract description 26
- 239000000463 material Substances 0.000 claims abstract description 24
- 210000000988 bone and bone Anatomy 0.000 claims abstract description 15
- 238000011156 evaluation Methods 0.000 claims description 30
- 230000001419 dependent effect Effects 0.000 claims description 17
- 210000004872 soft tissue Anatomy 0.000 claims description 10
- 239000000284 extract Substances 0.000 claims 3
- 230000035515 penetration Effects 0.000 claims 2
- 230000006870 function Effects 0.000 description 12
- 238000000034 method Methods 0.000 description 7
- 230000003595 spectral effect Effects 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 238000010521 absorption reaction Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 210000000746 body region Anatomy 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000013558 reference substance Substances 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- ZPWVASYFFYYZEW-UHFFFAOYSA-L dipotassium hydrogen phosphate Chemical compound [K+].[K+].OP([O-])([O-])=O ZPWVASYFFYYZEW-UHFFFAOYSA-L 0.000 description 1
- 229910000396 dipotassium phosphate Inorganic materials 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- 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 or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S378/00—X-ray or gamma ray systems or devices
- Y10S378/901—Computer tomography program or processor
Definitions
- the present invention is directed to beam-hardening correction in an X-ray computer tomography apparatus.
- Standard algorithms for instance, polynomial correction
- Polynomial correction are known for the correction of such image artifacts caused by beam hardening. These produce satisfactory results as long as the spectral absorption or attenuation behavior of the transirradiated body substances does not significantly differ from the spectral attenuation behavior of a reference substance for which the correction algorithm was developed.
- Water is used as the reference substance in the standard case since water exhibits a spectral attenuation behavior comparable to soft tissue in the human body and the human body is largely composed of soft tissue. Beam-hardening errors then can be eliminated to a significant extent in body regions where essentially only soft tissue is encountered.
- the algorithm is no longer accurate since bone tissue exhibits a spectral attenuation behavior that deviates substantially from water.
- the same is also true, for example, of vessels filled with contrast agent. Since the extent to which soft tissue and bone tissue were responsible for the beam attenuation is initially unknown for the measured values acquired in the course of the examination of a patient, a satisfactory beam hardening correction is not possible based solely on knowledge of the measured values.
- correction values that are added to the originally measured overall attenuation values are determined for the partial attenuation values of each sub-image.
- the correction values are taken from correction characteristics that were separately determined in advance for the respective body substances on the basis of reference materials with comparable attenuation.
- An overall image—which is now corrected for beam hardening—of the transirradiated body slice is reconstructed a second time from the corrected, overall attenuation values.
- An object of the present invention is to provide a computed tomography apparatus with improved beam-hardening correction.
- this object is achieved in an X-ray computed tomography apparatus constructed and operating as follows.
- a radiator/detector arrangement supplies a set of measured intensity values for each X-ray projection of a body slice of a patient under examination, each measured value thereof being representative of the intensity of the X-rays that have passed through the body slice in a respective partial projection region of the overall projection region.
- An electronic evaluation and reconstruction unit is connected to the radiator/detector arrangement and is configured for
- g corr g ⁇ k ( t 1 , t 2 ) (1).
- g is the overall attenuation value
- g corr is the overall attenuation value corrected for beam hardening
- t 1 is the first attenuation partial value
- t 2 is the second attenuation partial value
- k(t 1 , t 2 ) is the correction value dependent on t 1 and t 2 .
- a set of reference overall attenuation values g ref (s 1 , s 2 ) is determined for a material combination of a first reference material and a second reference material different therefrom.
- This set of reference overall attenuation values g ref (s 1 , s 2 ) is representative of the actual overall attenuation of the X-rays produced by this material combination at various respective thicknesses of the first material and the second reference material.
- s 1 references a first individual attenuation value that is representative of the theoretical linear attenuation of the X-rays by the first reference material for the respective thickness of the first reference material
- s 2 references a second individual attenuation value that is representative of the theoretical linear attenuation of the X-rays by the second reference material for the respective thickness of the second reference material.
- the evaluation and reconstruction unit determines (and uses) the aforementioned correction according to the following equation:
- the inventive X-ray computed tomography apparatus is constructed and operates as follows:
- a radiator/detector arrangement supplies a set of measured intensity values for each X-ray projection of a body slice of a patient under examination, each measured value thereof being representative of the intensity of the X-rays that have passed through the body slice in a respective partial projection region of the overall projection region.
- An electronic evaluation and reconstruction unit is connected to the radiator/detector arrangement and is configured for
- g is the overall attenuation value
- g corr is the overall attenuation value corrected for beam hardening
- k(t) is the correction value dependent on t.
- a set of reference overall attenuation values g ref (s 1 , s 2 ) is determined for a material combination of a first reference material and a second reference material different therefrom.
- This set of reference overall attenuation values g ref (s 1 , s 2 ) is representative of the actual overall attenuation of the X-rays produced by this material combination at various respective thicknesses of the first material and the second reference material.
- s is a first individual attenuation value that is representative of the theoretical linear attenuation of the X-rays by the first reference material for the respective thickness of the first reference material
- s 2 references a second individual attenuation value that is representative of the theoretical linear attenuation of the X-rays by the second reference material for the respective thickness of the second reference material.
- the evaluation and reconstruction unit is configured for determining (and using) the overall attenuation value dependent on the reference overall attenuation values according to the following equation applies:
- the two versions have in common the use of a correction value that takes the attenuation by a combination of two different materials into consideration. It has been shown in the human body that the beam hardening by one substance (for instance, bone tissue) is not independent of whether other substances (for instance, soft tissue) are additionally present in the beam path. However, the known retrospective correction methods are based precisely on the premise that this precondition of independency exits, by taking only the attenuation by a single (generalized) substance into consideration. By employing a correction value dependent on the attenuation of two materials, it is possible to come very close to the actual conditions in the human body. Images that are very low in disturbing image artifacts thus can be generated, particularly given exposures of body regions having a comparatively high proportion of bone.
- references materials Materials whose spectral attenuation behavior is similar to the body substances that are to be taken into consideration in the partial images are expediently selected as the reference materials.
- water For a partial image that should essentially show only soft tissue, it is expedient to select water as reference material.
- a mixture of K 2 HPO 4 and water can be selected as reference material (S.C.E. Cann, Radiology 166, pp. 509-522 (1988)).
- the single figure is a schematic illustration of a computed tomography apparatus constructed and operating in accordance with the principles of the present invention.
- the reference overall attenuation values for use in the inventive computed tomography apparatus can be determined either by computer simulation or by measurement. For the measurement, for example, one can proceed such as shown in the figure.
- This shows a CT scanner with an X-radiator 10 and a detector arrangement 12 that is composed of a number of detector cells 14 arranged side-by-side along a circular arc.
- the X-radiator 10 emits a fan-shaped X-ray beam in a plane.
- the detector cells 14 cover an overall projection region extending through the angular width of the beam fan, with each individual detector cell 14 detecting the intensity of the incident X-rays in the partial projection region it respectively covers.
- Each cell 14 supplies a corresponding measured intensity signal to an electronic evaluation and reconstruction unit 16 .
- the X-radiator 10 and the detector arrangement 12 can be moved in a direction on normal to the fan plane without rotation along an axis 18 .
- Two wedges 20 and 22 are situated in the beam path.
- the wedges 20 , 22 are arranged such that the thickness of the wedge 20 increases along the axis 18 and the thickness of the wedge 22 increases along the direction of the fan angle.
- measured values are obtained for a number of different thickness value pairs of the two reference materials.
- the evaluation and reconstruction unit 16 uses these measured values to calculate the reference overall attenuation values. These indicate the actual attenuation affected with beam hardening that the X-rays experiences for the respective thicknesses of the two reference materials.
- one of the reference overall attenuation values g ref (s 1 , s 2 ) can be unambiguously allocated to each pair of individual attenuation values s 1 and s 2 (potentially with the assistance of interpolations).
- the individual attenuation values s 1 , s 2 thereby respectively reference the theoretical linear attenuation that the X-rays would experience in case of energy-independent absorption in the first and second reference material. They are linked with the thickness of the material via the following relationship:
- d 1,2 is the thickness of the first material or the second reference material and ⁇ 1,2 is an absorption coefficient of the first or the second reference material that is effective for linear attenuation.
- a reference attenuation error e ref (s 1 , s 2 ) can then also be unambiguously allocated to each pair of individual attenuation values s 1 and s 2 , this deriving according to the following equation
- This reference attenuation error e ref (s 1 , s 2 ) represents the beam hardening error by which the reference overall attenuation value g ref (s 1 , s 2 ) is lower than the sum of the individual attenuation values s 1 and s 2 as a consequence of beam hardening.
- the other individual attenuation value s 2 or s can be unambiguously determined for each pair of a reference overall attenuation value g ref (s 1 , s 2 ) and one of the individual attenuation values s 1 and s 2 .
- the appertaining reference attenuation error e ref (s 1 , s 2 ) then can also be unambiguously determined.
- One or more gray scale value thresholds according to which the overall image is resolved into its various gray scale value regions can, for example, be defined in order to extract the partial images from the overall image. It has proven expedient when one of the two partial images in the first version, or the single partial image in the second version, essentially shows only bone substance that is present in the transirradiated region.
- the correction value is then determined dependent on the re-projected bone attenuation and a further attenuation partial value acquired by re-projection, preferably of the soft tissue attenuation; whereas, in the second solution, the correction value is determined dependent on the re-projected bone attenuation and in the overall attenuation value.
- the re-projected attenuation partial values are preferably attenuation values that indicate the theoretical linear attenuation in the respective body substance. Details regarding how the attenuation partial values can be determined from the partial images by re-projection can be derived, for example, from the previously cited literature.
- the identified reference overall attenuation values g ref (S 1 , s 2 ) can be stored in the form of a look-up table in a memory of the evaluation and reconstruction unit 16 dependent on the individual attenuation values s 1 and s 2 .
- the evaluation and reconstruction unit 16 would also execute the arithmetic operation according to Equation (2) in order to obtain the correction value.
- the reference attenuation errors e ref (s 1 , s 2 ) can be directly stored in table form in the evaluation and reconstruction unit 16 dependent on the individual attenuation values s 1 and s 2 .
- a look-up table can be stored in the evaluation and reconstruction unit 16 that indicates the second individual attenuation value s 2 dependent on the first individual attenuation value s 1 and on the reference overall attenuation value g ref (s 1 , s 2 ). So that the arithmetic operation according to Equation (4) need not be constantly carried out by the evaluation and reconstruction unit 16 , it is also possible to directly store the reference attenuation errors e ref (s 1 , s 2 ) in table form in the evaluation and reconstruction unit 16 dependent on the reference overall attenuation value g ref (s 1 , s 2 ) and on the first individual attenuation value s 1 .
- a function u(x) that is dependent on a variable x can be determined for the determination of the beam hardening correction information, this function allocating a function value u(x) to every value of x, with
- this function value u(x) at least approximately corresponds to the difference between the sum of the respective, two individual attenuation values and the respective reference overall attenuation value, whereby A, B and C are constants.
- the evaluation and reconstruction unit 16 is configured for determining the correction value according to the following equation:
- Equation (10) makes it possible to reduce the dependency of the correction value k(t 1 , t 2 ) on two parameters to the dependency on one parameter.
- a function v(y) dependent on a variable y can be determined for the determination of the beam hardening correction information, this function allocating a function value v(y) to each value of y, with
- the function value v(y) at least approximately corresponds to the difference between the sum of the respective, two individual attenuation values and the respective reference overall attenuation value, whereby D, E and F are constants.
- the evaluation and reconstruction unit is configured for determining the correction value according to the following equation:
- the two functions u(x), v(y) either can be implemented as a look-up table in the evaluation and reconstruction unit 16 or in the form of a mathematical algorithm, insofar as a suitable approximation equation for the function u(x) or v(y) can be found.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
In an X-ray computed tomography apparatus with retrospective beam hardening correction, an overall image of a body slice under examination is determined from overall attenuation values that are obtained from the body slice. At least one partial image that shows essentially only one body substance, such as bone substance, is extracted from this overall image. Attenuation partial values are employed for determining a correction value. The attenuation values are determined for each overall attenuation value from the at least one partial image by re-projection. A correction value is derived from the beam hardening error that is determined for a material combination of two different reference materials.
Description
- 1. Field of the Invention
- The present invention is directed to beam-hardening correction in an X-ray computer tomography apparatus.
- 2. Description of the Prior Art
- In X-ray computed tomography, a shift of the average energy of the X-rays toward higher values occurs as a consequence of the polychromatic spectrum of the radiation emitted by the X-ray source and as a consequence of the energy-dependent absorption of the X-rays in the body of the patient under examination. This effect is called beam hardening. This effect becomes more pronounced as the transirradiated path in the body becomes longer. In the reconstructed image of the transirradiated body slice, this beam-hardening effect leads to unwanted image artifacts that negatively affect the precise medical interpretation of the image.
- Standard algorithms (for instance, polynomial correction) are known for the correction of such image artifacts caused by beam hardening. These produce satisfactory results as long as the spectral absorption or attenuation behavior of the transirradiated body substances does not significantly differ from the spectral attenuation behavior of a reference substance for which the correction algorithm was developed. Water is used as the reference substance in the standard case since water exhibits a spectral attenuation behavior comparable to soft tissue in the human body and the human body is largely composed of soft tissue. Beam-hardening errors then can be eliminated to a significant extent in body regions where essentially only soft tissue is encountered. When, however, the X-rays also passes through bone tissue, the algorithm is no longer accurate since bone tissue exhibits a spectral attenuation behavior that deviates substantially from water. The same is also true, for example, of vessels filled with contrast agent. Since the extent to which soft tissue and bone tissue were responsible for the beam attenuation is initially unknown for the measured values acquired in the course of the examination of a patient, a satisfactory beam hardening correction is not possible based solely on knowledge of the measured values.
- Methods referred to as retrospective correction methods were therefore developed wherein an overall image of the transirradiated body slice is first reconstructed from the measured, overall attenuation values, and this overall image is subsequently analyzed and resolved into sub-images. Each of the sub-images shows only a part of the various body substances. In the standard case, a bone image and a soft tissue image are generated. Partial attenuation values that indicate the beam attenuation by the appertaining part of the body substances, i.e., for example, bone tissue or soft tissue, are then calculated from the individual sub-images by re-projection. Subsequently, correction values that are added to the originally measured overall attenuation values are determined for the partial attenuation values of each sub-image. For example, the correction values are taken from correction characteristics that were separately determined in advance for the respective body substances on the basis of reference materials with comparable attenuation. An overall image—which is now corrected for beam hardening—of the transirradiated body slice is reconstructed a second time from the corrected, overall attenuation values.
- More detailed information about retrospective (post construction) correction methods may be found, for example, in “A Comparative Study of two Postreconstruction Beam Hardening Correction Methods” by G. T. Herman, S. S. Trivedi, IEEE Transactions on Medical Imaging, MI-2, 1983, pp. 128 ff., and in “A Method for Correcting Bone Induced Artifacts in Computer Tomography Scanners” by P. M. Joseph, R. D. Spital, Journal of Computer Assisted Tomography, No. 2, 1978, pp. 100 ff.
- It has been shown in practice that the known retrospective correction methods can in fact clearly reduce image artifacts caused by beam hardening compared to traditional, standard algorithms, however, image artifacts continue to be observed and elimination or at least reduction thereof is desirable.
- An object of the present invention is to provide a computed tomography apparatus with improved beam-hardening correction.
- In a first version of the solution, this object is achieved in an X-ray computed tomography apparatus constructed and operating as follows.
- A radiator/detector arrangement supplies a set of measured intensity values for each X-ray projection of a body slice of a patient under examination, each measured value thereof being representative of the intensity of the X-rays that have passed through the body slice in a respective partial projection region of the overall projection region. An electronic evaluation and reconstruction unit is connected to the radiator/detector arrangement and is configured for
- a) determining an overall attenuation value for each measured intensity value, the overall attenuation value being representative of the actual overall attenuation of the X-rays produced in the body slice in the appertaining partial projection region;
- b) reconstructing an overall image of the body slice proceeding from the overall attenuation values;
- c) extracting a first partial image from this overall image wherein essentially only those image parts of the overall image are contained that correspond to a first part of the various substances occurring in the body slice;
- d) determining respective first attenuation partial values allocated respectively to each overall attenuation value on the basis of this first partial image, the first attenuation partial value being a criterion for the attenuation that the X-rays experiences in the respective projection partial region due to the first part of the substances;
- e) extracting a second partial image from the overall image of the body slice wherein essentially only those image parts of the overall image are contained that correspond to a second part of the substances in the body slice differing from the first part;
- f) determining respective, second attenuation partial values on the basis of this second partial image allocated respectively to each overall attenuation value, the second attenuation sub-value being a criterion for the attenuation that the X-radiation experiences in the respective projection partial region due to the second part of the substances;
- g) determining a correction value for every overall attenuation value on the basis of previously determined beam-hardening correction information stored in the evaluation and reconstruction unit and dependent on the two attenuation partial values; and
- h) determining an overall attenuation value corrected for beam hardening for each overall attenuation value according to the following equation:
- g corr=g−k(t 1 , t 2) (1).
- wherein g is the overall attenuation value, gcorr is the overall attenuation value corrected for beam hardening, t1 is the first attenuation partial value, t2 is the second attenuation partial value and k(t1, t2) is the correction value dependent on t1 and t2.
- For determining the beam hardening correction information, in accordance with the invention a set of reference overall attenuation values gref (s1, s2) is determined for a material combination of a first reference material and a second reference material different therefrom. This set of reference overall attenuation values gref (s1, s2) is representative of the actual overall attenuation of the X-rays produced by this material combination at various respective thicknesses of the first material and the second reference material. For this determination, s1 references a first individual attenuation value that is representative of the theoretical linear attenuation of the X-rays by the first reference material for the respective thickness of the first reference material, and s2 references a second individual attenuation value that is representative of the theoretical linear attenuation of the X-rays by the second reference material for the respective thickness of the second reference material. The evaluation and reconstruction unit determines (and uses) the aforementioned correction according to the following equation:
- k(t 1 , t 2)=t 1 +t 2 −g ref(s 1=t1 , s 2 =t 2) (2).
- In an alternative, second version, the inventive X-ray computed tomography apparatus is constructed and operates as follows:
- A radiator/detector arrangement supplies a set of measured intensity values for each X-ray projection of a body slice of a patient under examination, each measured value thereof being representative of the intensity of the X-rays that have passed through the body slice in a respective partial projection region of the overall projection region. An electronic evaluation and reconstruction unit is connected to the radiator/detector arrangement and is configured for
- a) determining an overall attenuation value for each measured intensity value, this overall attenuation value being representative of the actual overall attenuation of the X-rays produced in the body slice in the appertaining partial projection region;
- b) reconstructing an overall image of the body slice proceeding from the overall attenuation values;
- c) extracting a partial image from this overall image wherein essentially only those image parts of the overall image are contained that correspond to a first part of the various substances occurring in the body slice;
- d) determining respective attenuation partial values allocated respectively to each overall attenuation value on the basis of this partial image, the attenuation partial values being a criterion for the attenuation that the X-rays experiences in the respective projection partial region due to the first part of the substances;
- e) determining a correction value for every overall attenuation value on the basis of previously determined beam-hardening correction information stored in the evaluation and reconstruction unit and dependent on the respective attenuation sub-value; and
- f) determining an overall attenuation value corrected for beam hardening for each overall attenuation value according to the following equation:
- g corr g+k(t) (3),
- wherein g is the overall attenuation value, gcorr is the overall attenuation value corrected for beam hardening, and k(t) is the correction value dependent on t.
- For determining the beam hardening correction information in the second version of the invention, a set of reference overall attenuation values gref (s1, s2) is determined for a material combination of a first reference material and a second reference material different therefrom. This set of reference overall attenuation values gref (s1, s2) is representative of the actual overall attenuation of the X-rays produced by this material combination at various respective thicknesses of the first material and the second reference material. For this determination s, is a first individual attenuation value that is representative of the theoretical linear attenuation of the X-rays by the first reference material for the respective thickness of the first reference material, and s2 references a second individual attenuation value that is representative of the theoretical linear attenuation of the X-rays by the second reference material for the respective thickness of the second reference material. The evaluation and reconstruction unit is configured for determining (and using) the overall attenuation value dependent on the reference overall attenuation values according to the following equation applies:
- k(g, t)=t+s 2 −g ref(s 1 =t, s 2) (4),
- wherein
- g ref(s 1 =t, s 2)=g (5)
- applies for g ref(s 1 =t, s 2).
- The two versions have in common the use of a correction value that takes the attenuation by a combination of two different materials into consideration. It has been shown in the human body that the beam hardening by one substance (for instance, bone tissue) is not independent of whether other substances (for instance, soft tissue) are additionally present in the beam path. However, the known retrospective correction methods are based precisely on the premise that this precondition of independency exits, by taking only the attenuation by a single (generalized) substance into consideration. By employing a correction value dependent on the attenuation of two materials, it is possible to come very close to the actual conditions in the human body. Images that are very low in disturbing image artifacts thus can be generated, particularly given exposures of body regions having a comparatively high proportion of bone.
- Materials whose spectral attenuation behavior is similar to the body substances that are to be taken into consideration in the partial images are expediently selected as the reference materials. For a partial image that should essentially show only soft tissue, it is expedient to select water as reference material. For a partial image that should essentially show only bone substance, for example, a mixture of K2HPO4 and water can be selected as reference material (S.C.E. Cann, Radiology 166, pp. 509-522 (1988)).
- The single figure is a schematic illustration of a computed tomography apparatus constructed and operating in accordance with the principles of the present invention.
- The reference overall attenuation values for use in the inventive computed tomography apparatus can be determined either by computer simulation or by measurement. For the measurement, for example, one can proceed such as shown in the figure. This shows a CT scanner with an X-radiator10 and a
detector arrangement 12 that is composed of a number ofdetector cells 14 arranged side-by-side along a circular arc. The X-radiator 10 emits a fan-shaped X-ray beam in a plane. Thedetector cells 14 cover an overall projection region extending through the angular width of the beam fan, with eachindividual detector cell 14 detecting the intensity of the incident X-rays in the partial projection region it respectively covers. Eachcell 14 supplies a corresponding measured intensity signal to an electronic evaluation andreconstruction unit 16. The X-radiator 10 and thedetector arrangement 12 can be moved in a direction on normal to the fan plane without rotation along anaxis 18. - Two
wedges wedges wedge 20 increases along theaxis 18 and the thickness of thewedge 22 increases along the direction of the fan angle. When the X-radiator 10 and thedetector arrangement 12 are displaced along theaxis 18, measured values are obtained for a number of different thickness value pairs of the two reference materials. Using these measured values, the evaluation andreconstruction unit 16 then calculates the reference overall attenuation values. These indicate the actual attenuation affected with beam hardening that the X-rays experiences for the respective thicknesses of the two reference materials. In order to compensate individual errors of thedetector cells 14, it is recommended to also scan thewedge 22 once with a thickness gradient inverted along the fan angle direction, this being indicated in broken lines at 22′. - After the reference overall attenuation values are determined with the measuring structure according to the figure or by simulation, one of the reference overall attenuation values gref (s1, s2) can be unambiguously allocated to each pair of individual attenuation values s1 and s2 (potentially with the assistance of interpolations). The individual attenuation values s1, s2 thereby respectively reference the theoretical linear attenuation that the X-rays would experience in case of energy-independent absorption in the first and second reference material. They are linked with the thickness of the material via the following relationship:
- s 1,2 =d 1,2·μ1,2 (6),
- wherein d1,2 is the thickness of the first material or the second reference material and μ1,2 is an absorption coefficient of the first or the second reference material that is effective for linear attenuation.
- Moreover, a reference attenuation error eref (s1, s2) can then also be unambiguously allocated to each pair of individual attenuation values s1 and s2, this deriving according to the following equation
- e ref(s 1 , s 2)=s 1 +s 2 −g ref(s 1 , s 2) (7)
- from the difference between the sum of the individual attenuation values s1 and s2 and the appertaining reference overall attenuation value gref (s1, s2). This reference attenuation error eref (s1, s2) represents the beam hardening error by which the reference overall attenuation value gref (s1, s2) is lower than the sum of the individual attenuation values s1 and s2 as a consequence of beam hardening.
- In an analogous way, the other individual attenuation value s2 or s, can be unambiguously determined for each pair of a reference overall attenuation value gref (s1, s2) and one of the individual attenuation values s1 and s2. The appertaining reference attenuation error eref (s1, s2) then can also be unambiguously determined.
- The above considerations are utilized in the invention in order to find the respectively correct correction value during operation of the CT scanner given examination of a patient. In the first version of the invention, the two respectively identified attenuation values t1 and t2 are employed as parameters therefor; in the second version, one identified attenuation value t and the overall attenuation value g are employed. When s1=t1 and s2=t2 are set in the first version, then an appertaining reference overall attenuation value gref (s1=t1, s2=t2) and thus an appertaining reference attenuation error eref (s1=t1, s2=t2), thus can be immediately unambiguously determined. The value of this reference attenuation error is then employed as correction value k(t1, t2), i.e.
- k(t 1, t2)=e ref(s 1 =t 1 , s 2 =t 2)
- =t 1 +t 2 −g ref(s 1 =t 1 , s 2 =t 2) (8)
- The comparable case applies given the second version of the invention. When s1=t and gref (s1=t, s2) are set therein, an appertaining individual attenuation value s2, and thus an appertaining reference attenuation error eref (gref=g, s1=t), can be immediately unambiguously determined. The value of this reference attenuation error is then employed as correction value k(g, t), i.e
- k(g, t)=e ref(g ref =g, s 1 =t)
- =t=s 2 −g ref(s 1 =t, s 2) (9).
- One or more gray scale value thresholds according to which the overall image is resolved into its various gray scale value regions can, for example, be defined in order to extract the partial images from the overall image. It has proven expedient when one of the two partial images in the first version, or the single partial image in the second version, essentially shows only bone substance that is present in the transirradiated region. In the first version, the correction value is then determined dependent on the re-projected bone attenuation and a further attenuation partial value acquired by re-projection, preferably of the soft tissue attenuation; whereas, in the second solution, the correction value is determined dependent on the re-projected bone attenuation and in the overall attenuation value. The re-projected attenuation partial values are preferably attenuation values that indicate the theoretical linear attenuation in the respective body substance. Details regarding how the attenuation partial values can be determined from the partial images by re-projection can be derived, for example, from the previously cited literature.
- There are various possibilities regarding the concrete implementation of the beam hardening correction information in the inventive computed tomography apparatus. In the first version, the identified reference overall attenuation values gref(S1, s2) can be stored in the form of a look-up table in a memory of the evaluation and
reconstruction unit 16 dependent on the individual attenuation values s1 and s2. In this case, the evaluation andreconstruction unit 16 would also execute the arithmetic operation according to Equation (2) in order to obtain the correction value. Instead of the reference overall attenuation values, alternatively, the reference attenuation errors eref(s1, s2) can be directly stored in table form in the evaluation andreconstruction unit 16 dependent on the individual attenuation values s1 and s2. - In the second version, a look-up table can be stored in the evaluation and
reconstruction unit 16 that indicates the second individual attenuation value s2 dependent on the first individual attenuation value s1 and on the reference overall attenuation value gref (s1, s2). So that the arithmetic operation according to Equation (4) need not be constantly carried out by the evaluation andreconstruction unit 16, it is also possible to directly store the reference attenuation errors eref (s1, s2) in table form in the evaluation andreconstruction unit 16 dependent on the reference overall attenuation value gref (s1, s2) and on the first individual attenuation value s1. - In an optional embodiment of the first version, a function u(x) that is dependent on a variable x can be determined for the determination of the beam hardening correction information, this function allocating a function value u(x) to every value of x, with
- x=A(Bs 1 +Cs 2) (10)
- this function value u(x) at least approximately corresponds to the difference between the sum of the respective, two individual attenuation values and the respective reference overall attenuation value, whereby A, B and C are constants. The evaluation and
reconstruction unit 16 is configured for determining the correction value according to the following equation: - k(t 1 , t 2)=u(x=A(Bt 1 +Ct 2)) (11)
- The linear combination according to Equation (10) makes it possible to reduce the dependency of the correction value k(t1, t2) on two parameters to the dependency on one parameter. The constants A, B and C are defined such that the error between the reference attenuation error eref (s1=t1, s2=t2) and the function value u(x=A(Bt1+Ct2)) becomes optimally small.
- In an optional embodiment of the second version, analogously, a function v(y) dependent on a variable y can be determined for the determination of the beam hardening correction information, this function allocating a function value v(y) to each value of y, with
- y=D(Eg ref +Fs 1) (12)
- The function value v(y) at least approximately corresponds to the difference between the sum of the respective, two individual attenuation values and the respective reference overall attenuation value, whereby D, E and F are constants. The evaluation and reconstruction unit is configured for determining the correction value according to the following equation:
- k(g, t)=v(y=D(Eg+Ft)) (13)
- The constants D, E and F are defined in this case such that the error between the reference attenuation error eref (gref=g,t1, s1=t) and the function value v(y=D(Eg+ft)) becomes optimally small.
- The two functions u(x), v(y) either can be implemented as a look-up table in the evaluation and
reconstruction unit 16 or in the form of a mathematical algorithm, insofar as a suitable approximation equation for the function u(x) or v(y) can be found. - It is self-evident. moreover, that not only the calibration measurements for the determination of the reference overall attenuation values can be implemented at the CT scanner shown in the figure, but also that the patient examination can ensue in the scanner when, additionally, a rotation of the
radiator 10 and the detector arrangement 12 (insofar as this is not fashioned as ring detector) is provided in the direction of the fan angle. - It should be noted that it is possible without further difficulty to correct the overall attenuation values acquired in the examination of a patient a priori with the assistance of a standard algorithm, and to reconstruct the overall image from the overall attenuation values corrected in this way. The retrospectively determined correction values then will either be added to the original, non-corrected overall attenuation values, or the correction values will be reduced by the amount of the standard correction.
- 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.
Claims (7)
1. An X-ray computed tomography apparatus comprising:
a radiator/detector arrangement that supplies respective sets of measured intensity values for X-ray projections of a body slice of a patient, each measured intensity value of each set representing an intensity of the X-rays after penetration of the body slice in a partial projection region of an overall projection region; and
an electronic evaluation and reconstruction unit connected to the radiator/detector arrangement for
a) determining an overall attenuation value for each of the measured intensity values, said overall attenuation value representing an actual overall attenuation of the X-rays in the body slice in the appertaining partial projection region;
b) reconstructing an overall image of the body slice proceeding from the overall attenuation values;
c) extracting a first partially image of a first substance of the body slice from said overall image;
d) determining a respective first attenuation partial value from said first partial image for each partial region of the respective X-ray projection, each first attenuation partial value being a criterion for attenuation of the X-rays in the appertaining projection partial region in the first substance of the body slice, and allocating said first attenuation partial value to the overall attenuation value of the corresponding partial region;
e) extracting a second partial image from said overall image of a second substance in the body slice;
f) determining a respective second attenuation partial value from said second partial image for each partial region of the respective X-ray projection, each second attenuation partial value being a criterion for attenuation of the X-rays in the appertaining projection partial region in the second substance of the body slice, and allocating said second attenuation partial value to the overall attenuation value of the corresponding partial region;
g) determining a correction value for each overall attenuation value from the first and second attenuation partial values from predetermined beam-hardening correction information stored in the evaluation and reconstruction unit;
h) determining an overall attenuation value corrected for beam hardening for each overall attenuation value according to the equation:
g corr =g+k(t 1 , t 2)
wherein gcorr is the overall attenuation value corrected for beam hardening, g is the overall attenuation value, t1 is the first attenuation partial value, t2 is the second attenuation partial value and k(t1, t2) is the correction value determined from t1 and t2, and for determining the correction value k (t1, t2), said evaluation and reconstruction unit
i) determining a set of reference overall attenuation values gref (s1, s2) for a material sequence composed of a first reference material and a second reference material different therefrom for different thickness combinations, each reference overall attenuation value gref of the set representing the overall attenuation of the X-rays for a specific thickness combination of the material sequence;
j) calculating a theoretical linear attenuation of the X-rays for each thickness of the first reference material from the thickness combinations of the material sequence represented by a first individual attenuation value s1;
k) calculating a theoretical linear attenuation of the X-rays for each thickness of the second reference material from the thickness combinations of the material sequence represented by a second individual attenuation value s2; and
l) determining the correction value for each thickness combination of the material sequence according to the equation
k(t 1 , t 2)=t 1 +t 2 −g ref(s 1 =t 1 , s 2 =t 2),
wherein gref (s1=s2=t2) is the reference overall attenuation value for a thickness combination of the material sequence for which s1=t1 and s2=t2.
2. A computed tomography apparatus as claimed in claim 1 wherein the evaluation and reconstruction unit extracts substantially only image parts of the overall image that correspond to a bone structure in the body slice in one of the two partial images.
3. A computed tomography apparatus as claimed in claim 2 wherein the evaluation and reconstruction unit extracts substantially only image parts of the overall image that correspond to a soft tissue substance in the body slice in the other of the two partial images.
4. A computed tomography apparatus as claimed in claim 1 wherein, for determining the beam hardening correction information, the evaluation and reconstruction unit stores a function u(x) dependent on a variable x, which allocates a function value u(x) to each value of x, with
x=A(Bs 1 +Cs 2)1
said function value substantially corresponding to a difference between a sum of the respective, two individual attenuation values and the respective reference overall attenuation value, wherein A, B and C are constants; and wherein the evaluation and reconstruction unit determines the correction value according to the equation:
k(t 1 , t 2)=u(x=A(Bt 1 +Ct 2)).
5. An X-ray computed tomography apparatus comprising:
a radiator/detector arrangement that supplies respective sets of measured intensity values for X-ray projections of a body slice of a patient under examination, each measured intensity value of the set representing an intensity of X-rays after penetration of the body slice in a partial projection region of an overall projection region; and
an electronic evaluation and reconstruction unit connected to the radiator/detector arrangement for
a) determining an overall attenuation value for each of the measured intensity values, said overall attenuation value representing an actual overall attenuation of the X-rays in the body slice in the appertaining partial projection region;
b) reconstructing an overall image of the body slice proceeding from the overall attenuation values;
c) extracting a partial image of a selected substance of the body slice from said overall image;
d) determining an attenuation partial value from said partial image for each partial region of the respective X-ray projection, said attenuation partial value being a criterion for the attenuation of the X-rays in the appertaining projection partial region in the selected substance of the body slice, and allocating said attenuation partial value to the overall attenuation value of the corresponding partial region;
e) determining a correction value for each overall attenuation value from the attenuation partial value from predetermined beam-hardening correction information stored in the evaluation and reconstruction unit;
f) determining an overall attenuation value corrected for beam hardening for each overall attenuation value according to the equation:
g corr =g+k(t)
wherein gcorr is the overall attenuation value corrected for beam hardening, g is the overall attenuation value, t is the attenuation partial value, and k(t) is the correction value determined from t1 and, for determining the correction value k (t), said evaluation and reconstruction unit
g) determining a set of reference overall attenuation values for a material sequence composed of a first reference material and a second reference material different therefrom for different thickness combinations, each reference overall attenuation value gref of the set representing the overall attenuation of the X-rays for a specific thickness combination of the material sequence;
h) calculating a theoretical linear attenuation of the X-rays for each thickness of the first reference material from the thickness combinations of the material sequence represented by a first individual attenuation value s1;
i) calculating a theoretical linear attenuation of the X-rays for each thickness of the second reference material from the thickness combinations of the material sequence represented by a second individual attenuation value s2; and
j) determining the correction value for each thickness combination of the material sequence according to the equation
k(g, t)=t+s 2 +g ref(s 1 =t, s 2),
wherein gref (s1=t, s2) is the reference overall attenuation value for a thickness combination of the material sequence for which s1=t given s2 and, for which
g ref(s 1 =t, s 2).
6. A computed tomography apparatus as claimed in claim 5 wherein, for determining the beam hardening correction information, the evaluation and reconstruction unit stores a function v(y) dependent on a variable y which allocates a function value v(y) to each value of y, with
y=D(Eg ref 30 Fs1),
said function value substantially corresponding to a difference between a sum of the respective, two individual attenuation values and the respective reference overall attenuation value, wherein D, E and F are constants; and wherein the evaluation and reconstruction unit determines the correction value according to the equation
k(g, t)=v(y=D(Eg+Ft)).
7. A computed tomography apparatus as claimed in claim 5 wherein the evaluation and reconstruction unit extracts substantially only image parts of the overall image that correspond to a bone structure in the body slice in said partial image.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE10036142 | 2000-07-25 | ||
DE10036142.0 | 2000-07-25 | ||
DE10036142A DE10036142B4 (en) | 2000-07-25 | 2000-07-25 | X-ray computed tomography device |
Publications (2)
Publication Number | Publication Date |
---|---|
US20020015476A1 true US20020015476A1 (en) | 2002-02-07 |
US6430252B2 US6430252B2 (en) | 2002-08-06 |
Family
ID=7650115
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US09/912,137 Expired - Lifetime US6430252B2 (en) | 2000-07-25 | 2001-07-25 | X-ray computed tomography apparatus |
Country Status (4)
Country | Link |
---|---|
US (1) | US6430252B2 (en) |
JP (1) | JP2002095657A (en) |
CN (1) | CN1198538C (en) |
DE (1) | DE10036142B4 (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1389444A1 (en) * | 2002-08-16 | 2004-02-18 | GE Medical Systems Global Technology Company LLC | Method and system for implementing variable x-ray intensity modulation schemes for imaging systems |
WO2005076222A1 (en) * | 2004-02-09 | 2005-08-18 | Philips Intellectual Property & Standards Gmbh | Correction of artifacts caused by high attenuating objects |
WO2007081662A2 (en) * | 2006-01-05 | 2007-07-19 | Koninklijke Philips Electronics, N.V. | Adjustable phantom |
US20070293200A1 (en) * | 2004-02-20 | 2007-12-20 | Brian Roundtree | Call intercept methods, such as for customer self-support on a mobile device |
US20080273656A1 (en) * | 2005-12-15 | 2008-11-06 | Koninklijke Philips Electronics N. V. | Advanced Convergence for Multiple Iterative Algorithm |
US20090257549A1 (en) * | 2008-04-10 | 2009-10-15 | Bjorn Heismann | Method for preparing reconstructed CT image data records and CT system |
BE1020452A5 (en) * | 2011-05-04 | 2013-10-01 | Materialise Nv | CALIBRATION DEVICE FOR IMAGE. |
CN103971387A (en) * | 2013-01-29 | 2014-08-06 | 上海联影医疗科技有限公司 | Method for reconstructing CT image |
CN110636796A (en) * | 2017-04-20 | 2019-12-31 | 皇家飞利浦有限公司 | Beam hardening correction in X-ray dark field imaging |
US11094094B2 (en) | 2017-08-28 | 2021-08-17 | Shanghai United Imaging Healthcare Co., Ltd | System and method for removing hard tissue in CT image |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6516045B2 (en) | 2001-05-04 | 2003-02-04 | The Regents Of The University Of California | Device and method for determining proportions of body materials |
DE10202732A1 (en) * | 2002-01-24 | 2003-08-07 | Fraunhofer Ges Forschung | Device and method for creating a correction characteristic for reducing artifacts in a tomography |
JP4142482B2 (en) * | 2003-04-04 | 2008-09-03 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X-ray CT system |
JP3950811B2 (en) * | 2003-04-17 | 2007-08-01 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X-ray CT apparatus and beam hardening post-processing method |
WO2005059592A1 (en) * | 2003-12-16 | 2005-06-30 | Philips Intellectual Property & Standards Gmbh | Correction of artifacts caused by the heel effect |
US7391844B2 (en) * | 2005-01-14 | 2008-06-24 | General Electric Company | Method and apparatus for correcting for beam hardening in CT images |
US8301123B2 (en) * | 2005-06-24 | 2012-10-30 | Nuance Communications, Inc. | Local intercept methods, such as applications for providing customer assistance for training, information calls and diagnostics |
US8682298B2 (en) * | 2005-10-12 | 2014-03-25 | Nuance Communications, Inc. | Message intercept methods, such as for customer self-support on a mobile device |
EP1969886B1 (en) | 2005-12-13 | 2017-09-20 | Nuance Communications, Inc. | Method for performing interactive services on a mobile device, such as time or location initiated interactive services |
CN100365664C (en) * | 2006-03-02 | 2008-01-30 | 西安交通大学 | CT projection data beam hardening effect correcting method based on IIL consistency condition |
US8744414B2 (en) * | 2007-01-05 | 2014-06-03 | Nuance Communications, Inc. | Methods of interacting between mobile devices and voice response systems |
WO2008128119A1 (en) | 2007-04-12 | 2008-10-23 | Snapin Software Inc. | System and method for detecting mutually supported capabilities between mobile devices |
DE102007046359B4 (en) * | 2007-09-27 | 2016-02-04 | Siemens Aktiengesellschaft | Method and device for creating material-selective volume images |
US8589149B2 (en) | 2008-08-05 | 2013-11-19 | Nuance Communications, Inc. | Probability-based approach to recognition of user-entered data |
DE102009010287B4 (en) | 2009-02-24 | 2016-06-30 | Siemens Healthcare Gmbh | Method and image reconstruction device for the reconstruction of image data |
DE102009014726A1 (en) | 2009-03-25 | 2010-10-07 | Siemens Aktiengesellschaft | Method and image reconstruction device for the reconstruction of image data |
US8594759B2 (en) * | 2009-07-30 | 2013-11-26 | Nellcor Puritan Bennett Ireland | Systems and methods for resolving the continuous wavelet transform of a signal |
US9374565B2 (en) | 2014-11-06 | 2016-06-21 | General Electric Company | System and method for multi-material correction of image data |
US10497153B2 (en) | 2018-02-08 | 2019-12-03 | FMI Medical Systems Co., Ltd. | Heel effect correction in computed tomography |
CN113796879B (en) * | 2021-09-27 | 2024-03-08 | 北京万东医疗科技股份有限公司 | Bulb tube emergent energy spectrum verification method and device, electronic equipment and storage medium |
CN115196250B (en) * | 2022-09-19 | 2023-01-24 | 山西潞安环保能源开发股份有限公司五阳煤矿 | Foreign matter identification method, device and system and storage medium |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5953444A (en) * | 1997-10-22 | 1999-09-14 | University Of Pennsylvania | Method for improved correction of spectrum hardening artifacts in computed tomography images |
DE10035984C1 (en) * | 2000-07-24 | 2002-01-31 | Siemens Ag | X-ray computer tomography apparatus, determines the attenuation correction function appropriate to patient, from measurements on filter and reference materials |
-
2000
- 2000-07-25 DE DE10036142A patent/DE10036142B4/en not_active Expired - Fee Related
-
2001
- 2001-07-19 JP JP2001219090A patent/JP2002095657A/en not_active Withdrawn
- 2001-07-24 CN CNB011246499A patent/CN1198538C/en not_active Expired - Fee Related
- 2001-07-25 US US09/912,137 patent/US6430252B2/en not_active Expired - Lifetime
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1389444A1 (en) * | 2002-08-16 | 2004-02-18 | GE Medical Systems Global Technology Company LLC | Method and system for implementing variable x-ray intensity modulation schemes for imaging systems |
WO2005076222A1 (en) * | 2004-02-09 | 2005-08-18 | Philips Intellectual Property & Standards Gmbh | Correction of artifacts caused by high attenuating objects |
US7676221B2 (en) | 2004-02-20 | 2010-03-09 | Snapin Software Inc. | Call intercept methods, such as for customer self-support on a mobile device |
US20070293200A1 (en) * | 2004-02-20 | 2007-12-20 | Brian Roundtree | Call intercept methods, such as for customer self-support on a mobile device |
US7961839B2 (en) * | 2005-12-15 | 2011-06-14 | Koninklijke Philips Electronics N.V. | Advanced convergence for multiple iterative algorithm |
US20080273656A1 (en) * | 2005-12-15 | 2008-11-06 | Koninklijke Philips Electronics N. V. | Advanced Convergence for Multiple Iterative Algorithm |
WO2007081662A3 (en) * | 2006-01-05 | 2008-02-07 | Koninkl Philips Electronics Nv | Adjustable phantom |
US20090052755A1 (en) * | 2006-01-05 | 2009-02-26 | Koninklijke Philips Electronics N. V. | Adjustable phantom |
US7738624B2 (en) | 2006-01-05 | 2010-06-15 | Koninklijke Philips Electronics N.V. | Adjustable phantom |
WO2007081662A2 (en) * | 2006-01-05 | 2007-07-19 | Koninklijke Philips Electronics, N.V. | Adjustable phantom |
US20090257549A1 (en) * | 2008-04-10 | 2009-10-15 | Bjorn Heismann | Method for preparing reconstructed CT image data records and CT system |
US7889834B2 (en) * | 2008-04-10 | 2011-02-15 | Siemens Aktiengesellschaft | Method for preparing reconstructed CT image data records and CT system |
BE1020452A5 (en) * | 2011-05-04 | 2013-10-01 | Materialise Nv | CALIBRATION DEVICE FOR IMAGE. |
CN103971387A (en) * | 2013-01-29 | 2014-08-06 | 上海联影医疗科技有限公司 | Method for reconstructing CT image |
CN110636796A (en) * | 2017-04-20 | 2019-12-31 | 皇家飞利浦有限公司 | Beam hardening correction in X-ray dark field imaging |
US11094094B2 (en) | 2017-08-28 | 2021-08-17 | Shanghai United Imaging Healthcare Co., Ltd | System and method for removing hard tissue in CT image |
Also Published As
Publication number | Publication date |
---|---|
DE10036142B4 (en) | 2004-04-29 |
DE10036142A1 (en) | 2002-02-21 |
CN1198538C (en) | 2005-04-27 |
JP2002095657A (en) | 2002-04-02 |
US6430252B2 (en) | 2002-08-06 |
CN1334064A (en) | 2002-02-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6430252B2 (en) | X-ray computed tomography apparatus | |
US9036879B2 (en) | Multi-material decomposition using dual energy computed tomography | |
CN102048552B (en) | Beam hardening correction for ct perfusion measurements | |
US7760848B2 (en) | Method and system for generating a multi-spectral image of an object | |
EP1928318B1 (en) | Quantitative material decomposition for spectral ct | |
US7260169B2 (en) | Device and method for computer tomography | |
JP5237812B2 (en) | Data processing and analysis in computed tomography with multiple energy windows | |
US20070092056A1 (en) | Method for generating CT displays in x ray computed tomography | |
US9211104B2 (en) | Selection method for two contrast media for use in a dual-energy CT examination, contrast media combination and generation of CT images using a contrast media combination and different energy spectra | |
US8155422B2 (en) | Dynamic optimization of the signal-to-noise ratio of dual-energy attenuation data for reconstructing images | |
Cann | Low-dose CT scanning for quantitative spinal mineral analysis. | |
Elleaume et al. | Performance of computed tomography for contrast agent concentration measurements with monochromatic x-ray beams: comparison of K-edge versus temporal subtraction | |
AU2016222642B2 (en) | Methods for physiological state determination in body scans | |
US20090257549A1 (en) | Method for preparing reconstructed CT image data records and CT system | |
US6438197B2 (en) | X-ray computed tomography apparatus with correction for beam hardening | |
US20170311917A1 (en) | Spectrum-neutral determination of calcareous deposits in blood vessels | |
US9943279B2 (en) | Methods and systems for task-based data generation and weighting for CT spectral imaging | |
WO2010015953A2 (en) | Spectral imaging | |
Fuchs et al. | Coronary artery calcium quantification from contrast enhanced CT using gemstone spectral imaging and material decomposition | |
Guy et al. | DETECT-dual energy transmission estimation CT-for improved attenuation correction in SPECT and PET | |
US11517278B2 (en) | System and method for basis material decomposition with general physical constraint for multi-energy computed tomography | |
Kojima et al. | A novel fast kilovoltage switching dual-energy CT with deep learning: accuracy of CT number on virtual monochromatic imaging and iodine quantification | |
US6324240B1 (en) | Method for beam hardening correction in quantitative computed X-ray tomography | |
Zhang et al. | Accurate measurement of bone mineral density using clinical CT imaging with single energy beam spectral intensity correction | |
Ritchings et al. | A technique for simultaneous dual energy S |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: SIEMENS AKTIENGESELLSCHAFT, GERMANY Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:REINWAND, MARIO;STIERSTORFER, KARL;REEL/FRAME:012020/0487;SIGNING DATES FROM 20010720 TO 20010723 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
FEPP | Fee payment procedure |
Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
FPAY | Fee payment |
Year of fee payment: 8 |
|
FPAY | Fee payment |
Year of fee payment: 12 |