US20020015476A1 - X-ray computed tomography apparatus - Google Patents

X-ray computed tomography apparatus Download PDF

Info

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
Application number
US09/912,137
Other versions
US6430252B2 (en
Inventor
Mario Reinwand
Karl Stierstorfer
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens AG
Original Assignee
Siemens AG
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Siemens AG filed Critical Siemens AG
Assigned to SIEMENS AKTIENGESELLSCHAFT reassignment SIEMENS AKTIENGESELLSCHAFT ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: REINWAND, MARIO, STIERSTORFER, KARL
Publication of US20020015476A1 publication Critical patent/US20020015476A1/en
Application granted granted Critical
Publication of US6430252B2 publication Critical patent/US6430252B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

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

    BACKGROUND OF THE INVENTION
  • 1. Field of the Invention [0001]
  • The present invention is directed to beam-hardening correction in an X-ray computer tomography apparatus. [0002]
  • 2. Description of the Prior Art [0003]
  • 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. [0004]
  • 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. [0005]
  • 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. [0006]
  • 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. [0007]
  • 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. [0008]
    • SUMMARY OF THE INVENTION
  • An object of the present invention is to provide a computed tomography apparatus with improved beam-hardening correction. [0009]
  • In a first version of the solution, this object is achieved in an X-ray computed tomography apparatus constructed and operating as follows. [0010]
  • 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 [0011]
  • 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; [0012]
  • b) reconstructing an overall image of the body slice proceeding from the overall attenuation values; [0013]
  • 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; [0014]
  • 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; [0015]
  • 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; [0016]
  • 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; [0017]
  • 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 [0018]
  • h) determining an overall attenuation value corrected for beam hardening for each overall attenuation value according to the following equation: [0019]
  • g corr=g−k(t 1 , t 2)   (1).
  • wherein g is the overall attenuation value, g[0020] corr 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 g[0021] ref (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: [0022]
  • 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 [0023]
  • 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; [0024]
  • b) reconstructing an overall image of the body slice proceeding from the overall attenuation values; [0025]
  • 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; [0026]
  • 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; [0027]
  • 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 [0028]
  • f) determining an overall attenuation value corrected for beam hardening for each overall attenuation value according to the following equation: [0029]
  • g corr g+k(t)   (3),
  • wherein g is the overall attenuation value, g[0030] corr 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 g[0031] ref (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 [0032]
  • 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. [0033]
  • 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 K[0034] 2HPO4 and water can be selected as reference material (S.C.E. Cann, Radiology 166, pp. 509-522 (1988)).
    •  DESCRIPTION OF THE DRAWING
  • The single figure is a schematic illustration of a computed tomography apparatus constructed and operating in accordance with the principles of the present invention. [0035]
    • DESCRIPTION OF THE PREFERRED EMBODIMENTS
  • 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 [0036] 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 [0037] wedges 20 and 22, respectively composed of different one of the two reference materials, 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. When the X-radiator 10 and the detector arrangement 12 are displaced along the axis 18, measured values are obtained for a number of different thickness value pairs of the two reference materials. Using these measured values, the evaluation and reconstruction 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 the detector cells 14, it is recommended to also scan the wedge 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 g[0038] ref (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 d[0039] 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.
  • Moreover, a reference attenuation error e[0040] ref (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 s[0041] 1 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 s[0042] 2 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 t[0043] 1 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 g[0044] ref (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. [0045]
  • 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 g[0046] ref(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 and reconstruction 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 and reconstruction 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 [0047] 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 and reconstruction unit 16, it is also possible to directly store the reference attenuation errors eref (s1, s2) in table form in the evaluation and reconstruction 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 [0048]
  • 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 [0049] 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(t[0050] 1, 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 [0051]
  • 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: [0052]
  • 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 e[0053] ref (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 [0054] 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 [0055] 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. [0056]
  • 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. [0057]

Claims (7)

We claim as our invention:
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.
US09/912,137 2000-07-25 2001-07-25 X-ray computed tomography apparatus Expired - Lifetime US6430252B2 (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Cited By (16)

* Cited by examiner, † Cited by third party
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