US20100215233A1 - Method and system for generating a computed tomography image - Google Patents

Method and system for generating a computed tomography image Download PDF

Info

Publication number
US20100215233A1
US20100215233A1 US12/393,816 US39381609A US2010215233A1 US 20100215233 A1 US20100215233 A1 US 20100215233A1 US 39381609 A US39381609 A US 39381609A US 2010215233 A1 US2010215233 A1 US 2010215233A1
Authority
US
United States
Prior art keywords
projections
image
energy level
synthesized
computed tomography
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
US12/393,816
Other versions
US8050479B2 (en
Inventor
Jiang Hsieh
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.)
General Electric Co
Original Assignee
General Electric Co
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 General Electric Co filed Critical General Electric Co
Priority to US12/393,816 priority Critical patent/US8050479B2/en
Assigned to GENERAL ELECTRIC COMPANY reassignment GENERAL ELECTRIC COMPANY ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: HSIEH, JIANG
Priority to NL2004249A priority patent/NL2004249C2/en
Priority to JP2010033080A priority patent/JP5491230B2/en
Priority to CN2010101323581A priority patent/CN101912272B/en
Publication of US20100215233A1 publication Critical patent/US20100215233A1/en
Application granted granted Critical
Publication of US8050479B2 publication Critical patent/US8050479B2/en
Active legal-status Critical Current
Adjusted expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/40Arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/405Source units specially adapted to modify characteristics of the beam during the data acquisition process
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/482Diagnostic techniques involving multiple energy imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
    • A61B6/5235Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
    • A61B6/5241Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT combining overlapping images of the same imaging modality, e.g. by stitching
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/54Control of apparatus or devices for radiation diagnosis
    • A61B6/542Control of apparatus or devices for radiation diagnosis involving control of exposure
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/006Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/408Dual energy

Definitions

  • This disclosure relates generally to the field of medical imaging.
  • the following techniques relate to multiple-energy computed tomography imaging systems and a method of generating a computed tomography image.
  • Computed tomography (CT) imaging systems operate by projecting an x-ray beam from an x-ray source through an attenuating object, such as a patient. The attenuated x-ray beam is then detected by a detector assembly. Some materials, such as bone, are much more x-ray attenuating than other materials, such as soft tissue.
  • Conventional third-generation CT imaging systems acquire attenuation data by projecting a polychromatic x-ray beam from the x-ray source.
  • the polychromatic x-ray beam contains many different frequencies of x-rays and is typically centered around a particular energy level. With sufficient angular coverage around the patient, cross-sectional images can be formed revealing the inner structure of the scanned object.
  • the images are typically displayed on a flat-screen monitor or cathode ray tube.
  • a virtual 3-D image may also be produced based on data acquired during a computed tomography scan.
  • a processor may rapidly switch the output from a generator so that the input voltage to the x-ray source changes from projection to projection.
  • a typical dual-energy CT imaging system may rapidly alternate between acquiring a high-energy projection and a low-energy projection.
  • the time spent transitioning between the high-energy level and the low-energy level must be short relative to the duration of each projection. Otherwise, the effective energy level during the high-energy projections will be less than the desired high-energy setting and the effective energy level during the low-energy projections will be greater than the desired low-energy setting.
  • the generator must transition between the high-energy level and the low-energy level in significantly less than a fraction of a millisecond. Enabling a generator to switch this quickly places an enormous demand on the generator hardware. Predictably, the demands on hardware will further increase as scan speed increases.
  • a method of generating a computed tomography image includes acquiring a first plurality of projections at a first energy level and acquiring a second plurality of projections at a second energy level.
  • the method includes reconstructing an image from the first plurality of projections.
  • the method includes generating a synthesized projection from the image that corresponds to one of the second plurality of projections.
  • the method includes comparing the synthesized projection to the one of the second plurality of projections.
  • the method also includes modifying the image to form an updated image based on the comparing the synthesized projection to the one of the second plurality of projections.
  • a method of generating a computed tomography image includes acquiring a first plurality of projections at a first energy level.
  • the method includes acquiring a second plurality of projections at a second energy level.
  • the method includes reconstruction an image from the first plurality of projections.
  • the method includes generating a plurality of synthesized projections from the image, each of the plurality of synthesized projections corresponding to one of the second plurality of projections.
  • the method includes comparing each of the plurality of synthesized projections to the corresponding one of the second plurality of projections.
  • the method also includes modifying the image to form an updated image based on the comparing each of the plurality of synthesized projections to the corresponding one of the second plurality of projections.
  • a method of generating a computed tomography image includes acquiring a first plurality of projections of an object at a first energy level.
  • the method includes acquiring a second plurality of projections of the object at a second energy level.
  • the method includes reconstructing an image of the object at the first energy level from the first plurality of projections.
  • the method also includes modifying the image to generate an updated image of the object, wherein the updated image represents the object at the second energy level.
  • a computed tomography imaging system includes a gantry, an x-ray source attached to the gantry, a table adapted to translate with respect to the gantry, and a processor in electrical communication with the gantry, the x-ray source, and the table.
  • the processor is configured to acquire a first plurality of projections at a first energy level and to acquire a second plurality of projections at a second energy level.
  • the processor is configured to reconstruct an image from the first plurality of projections.
  • the processor is configured to generate a synthesized projection from the image that corresponds to one of the second plurality of projections.
  • the processor is configured to determine a difference between the synthesized projection and the one of the second plurality of projections.
  • the processor is also configured to generate an updated image by modifying the first image based on the difference between the synthesized projection and the one of the second plurality of projections.
  • FIG. 1 is a schematic representation of a computed tomography imaging system in accordance with an embodiment
  • FIG. 2 is a flow chart illustrating a method of generating a computed tomography image in accordance with an embodiment.
  • the CT imaging system 10 includes a gantry support 12 , a gantry 14 , a table support 15 , a table 16 , an x-ray generator (not shown), an x-ray source 18 , a detector assembly 20 , and a processor 22 .
  • the gantry 14 is configured to rotate within the gantry support 12 .
  • the gantry 14 is adapted to retain the x-ray source 18 and the detector assembly 20 .
  • the x-ray generator is configured to deliver both a high-voltage input and a low-voltage input to the x-ray source 18 .
  • the high-voltage input may comprise a 140 kVp input while the low-voltage input may comprise an 80 kVp input.
  • the x-ray source 18 receives either the high-voltage input or the low-voltage input from the generator and emits an x-ray beam. If the x-ray source 18 receives the high-voltage input, the x-ray beam will comprise a high-energy x-ray beam. If the x-ray source 18 receives the low-voltage input, the x-ray beam will comprise a low-energy x-ray beam. Both the high-energy x-ray beam and the low-energy x-ray beam comprise photons from a wide energy spectrum.
  • high-energy x-ray beam indicates that the highest energy x-ray photons in the high-energy x-ray beam are of high energy and that the term “low-energy x-ray beam” indicates that the highest energy x-ray photons in the low-energy x-ray beam are of low energy.
  • high-energy x-ray beam and low-energy x-ray beam are relative to each other, and they should not be interpreted as indicating an absolute level of energy.
  • the x-ray source 18 is configured to emit either a high-energy x-ray beam or a low-energy x-ray beam through a patient 24 being examined. After passing through the patient 24 , either the high-energy x-ray beam or the low-energy x-ray beam is received at the detector assembly 20 .
  • the detector assembly 20 comprises a plurality of detector elements (not shown). Each of the plurality of detector elements produces an electrical signal that varies based on the intensity of the x-ray beam received during a sampling interval.
  • the table 16 is adapted to translate the patient 24 in a z-direction with respect to the gantry 14 as indicated by a coordinate axis 26 .
  • the processor 22 is configured to control the rotation of the gantry 14 , the position of the table 16 , and the activation of the x-ray source 18 .
  • FIG. 2 is a flow chart representing the method 100 in accordance with an embodiment.
  • the individual blocks 102 - 116 represent steps that may be performed in accordance with the method 100 .
  • the steps 102 - 116 of the method 100 need not be performed in the order shown.
  • the technical effect of the method 100 is the generation of computed tomography image.
  • the processor 22 may be configured to implement the steps illustrated in the method 100 .
  • the processor 22 controls the generator (not shown) in a manner so that the x-ray source 18 emits an x-ray beam at a first energy level.
  • the first energy level may be a relatively high energy level, such as 140 kVp.
  • the processor 22 controls the generator, the x-ray source 18 , the gantry 14 , and the table 16 so that a first plurality of projections are acquired at the first energy level.
  • the first plurality of projections may be acquired as part of a helical acquisition according to an embodiment and as part of an axial acquisition according to another embodiment.
  • a projection may be defined to include data obtained at generally the same time from some or all of the detector elements for a given gantry angle and location along the z-axis.
  • a projection may comprise the output from all of the detector elements at generally the same time.
  • the number of projections that can be acquired depends upon the speed of gantry rotation, the type of hardware used in the detector assembly 20 , and the number of detector elements among other factors.
  • 984 projections may comprise a complete set of projections for one gantry revolution.
  • the processor 22 controls the generator (not shown) in a manner so that the x-ray source 18 emits an x-ray beam at a second energy level.
  • the second energy level may be a relatively low energy level, such as 80 kVp.
  • the processor 22 controls the generator, the x-ray source 18 , the gantry 14 , and the table 16 so that a second plurality of projections are obtained at the second energy level.
  • the second plurality of projections may be acquired as part of a helical acquisition according to an embodiment or as part of an axial acquisition according to another embodiment. Steps 102 and 104 may be performed simultaneously, sequentially, or a combination of simultaneously and sequentially.
  • the first plurality of projections are acquired during generally the same period of time that the second plurality of projections are acquired.
  • the method 100 may alternate between acquiring one or more projections at the first energy level and acquiring one or more projections at the second energy level.
  • the first plurality of projections would comprise all of the projections acquired at the first energy level and the second plurality of projection would comprise all of the projections acquired at the second energy level.
  • an image produced by an x-ray beam at a higher energy level may generally produce less noise than an image produced by an x-ray beam at a lower energy level if all other imaging parameters are held constant. Therefore, according to an embodiment, it may be preferable to collect more projections at the higher energy level than at the lower energy level.
  • the processor 22 may control the generator (not shown) and x-ray source 18 in a manner so that one lower energy projection is acquired for every 24 projections.
  • projections 1 through 23 may be acquired at a first energy level, such as 140 kVp, and projection 24 may be acquired at a second energy level, such as 80 kVp.
  • All 984 projections of the exemplary embodiment may be acquired by alternating between acquiring 23 projections at 140 kVp and acquiring one projection at 80 kVp. Therefore, according to the exemplary embodiment, the result will be 943 projections at 140 kVp and 41 projections at 80 kVp. It should be understood that other embodiments may collect the first plurality of projections and the second plurality of projection in an entirely different manner.
  • the first plurality of projections comprises 943 projection acquired at 140 kVp and the second plurality of projections comprises 41 projections acquired at 80 kVp.
  • the 41 projections acquired at 80 kVp may be separated by a common angle.
  • some or all of the 41 projections acquired at 80 kVp may be separated by an amount that varies.
  • some of the projections may be spaced apart by 2 degrees or less and others of the projections may be spaced apart by more than 10 degrees.
  • the amount of angular separation between projections at a particular energy level may be determined by many factors including hardware limitations of the CT imaging system 10 (shown in FIG. 1 ).
  • a half-scan may be used during the acquisition of the first plurality of projections and the second plurality of projections to improve the temporal resolution of the CT imaging system.
  • Half-scans are well-known by those skilled in the art.
  • a total of 616 projections are acquired at 140 kVp and 27 projections are acquired at 80 kVp.
  • Other embodiments may acquire a different ratio of projections acquired at a first energy level to projections acquired at a second energy level.
  • projections may be acquired at three or more different energy levels.
  • a first plurality of projections may comprise 902 projections at 140 kVp
  • a second plurality of projections may comprise 41 projections at 100 kVp
  • a third plurality of projections may comprise 41 projections at 80 kVp.
  • the first plurality of projections, the second plurality of projections, and the third plurality of projections may be acquired in an interleaved manner.
  • the processor 22 shown in FIG. 1
  • the generator may utilize a slower transition when switching between energy levels.
  • 943 projections are acquired at 140 kVp and 41 projections are acquired at 80 kVp.
  • Each of the 41 projections acquired at 80 kVp may be evenly spaced apart.
  • the projections that are immediately before or after each of the projections acquired at 80 kVp projection may be collected while the generator is in a transition zone.
  • the term “transition zone” includes a time when the generator is between two of the desired voltages.
  • the spectrum of the x-ray beam may be somewhere between the desired levels of 80 kVp and 140 kVp during each of the transition zones. The demand on the generator is significantly reduced by acquiring projections during the transition zones instead of requiring that the generator is at exactly 140 kVp immediately before and after each 80 kVp projection.
  • the x-ray tube current used when acquiring the projections at the first energy level may be substantially different than the x-ray tube current used when acquiring the projections at the second energy level.
  • the processor 22 may utilize a reduced x-ray tube current level while acquiring the projections at the higher energy level.
  • the reduced x-ray tube current level is reduced in comparison to the x-ray tube current level used while acquiring the projections at the lower energy level.
  • the duration of a projection acquired at first energy level may be significantly different than the duration of a view acquired at second energy level.
  • an image is reconstructed from the first plurality of projections.
  • the image is reconstructed from the 943 projections acquired at 140 kVp. Since 943 projections at 140 kVp represents an almost complete set of projections, the first image may be reconstructed by using a conventional reconstruction algorithm, such as a filtered back-projection algorithm. If a filtered back-projection algorithm is used to reconstruct the first image, it may be necessary to estimate any missing projections.
  • the term “image” includes a reconstructed volume from one or more slices of computed tomography data.
  • an image may comprise data from a plurality of slices.
  • 943 projections were acquired 140 kVp and 41 of the projections were acquired at 80 kVp.
  • This embodiment is missing 41 projections from the 140 kVp set of projections.
  • the missing projections may be estimated by performing a linear interpolation or a high-order interpolation based on adjacent projections that were acquired at 140 kVp.
  • the first image may be reconstructed through other techniques according to additional embodiments. For example, an iterative reconstruction technique may be used to generate the first image.
  • a plurality of synthesized projections are generated.
  • the plurality of synthesized projections are based on the image the first time the method 100 reaches step 108 .
  • a plurality of synthesized projections may be generated based on the updated image in subsequent iterations as will be described hereinafter.
  • the term “synthesized projection” is defined to include a forward projection from an image. Normally, a CT imaging system collects a plurality of projections and uses the projections to reconstruct an image.
  • each of the synthesized projections generated at step 108 may correspond to one of the second plurality of projections.
  • a projection and a synthesized projection are considered to correspond if they both represent generally the same location and gantry angle of the same object. Additional normalization steps may be used to adjust the scale of the synthesized projection.
  • one of the plurality of synthesized projections generated at step 108 may represent the same location and gantry angle as one of the first plurality of projections acquired at step 102 . Then, a scaling factor or a plurality of polynomial parameters may be determined based on a comparison of the synthesized projection and the corresponding one of the first plurality of projections. The scaling factor or the plurality of polynomial parameters may be used to adjust each of the synthesized projections generated at step 108 .
  • each of the plurality of synthesized projections generated at step 108 corresponds to one of the second plurality of projections. Therefore, each one of the synthesized projections is generated at the same gantry angle and location in the z-direction as one of the second plurality of projections.
  • the second plurality of projections comprises 41 projections that were acquired at 80 kVp.
  • Each of the plurality of synthesized projections corresponds to one of the 41 projections that were acquired at 80 kVp.
  • the plurality of synthesized projections represent the same locations and angles as the 41 projections that were acquired at 80 kVp during step 104 . It is not necessary for the plurality of synthesized projections and the second plurality of projections to be of the same number in all embodiments.
  • step 108 there may be only four synthesized projections generated during step 108 .
  • the four synthesized projections may be separated by approximately 45 degrees from each other. It should be appreciated that each of the four synthesized projections still corresponds to one of the second plurality of projections acquired during step 104 .
  • step 108 may be replaced with a step where only one synthetic projection is generated for each iteration.
  • each of the plurality of synthesized projections is compared to the corresponding one of the second plurality of projections.
  • a difference and/or a ratio may be calculated between each synthesized projection and the corresponding one of the second plurality of projections. For example, if the synthesized projections were generated from an image acquired at 140 kVp and the second plurality of projections were acquired at 80 kVp, any difference between a synthesized projection and the corresponding one of the second plurality of projections may be assumed to be from the mean energy difference between the x-ray beam produced at 140 kVp and the x-ray beam produced at 80 kVp.
  • the first image is modified based on the comparison performed at step 110 to form an updated image.
  • a difference between each of the synthesized projections generated at step 108 and the corresponding projection acquired as part of the second plurality of projections during step 104 is distributed according to the following equation:
  • ⁇ f ( x, y, z ) ⁇ [ f ( x, y, z )] ⁇ [ ⁇ p ( ⁇ , ⁇ , ⁇ )]
  • f(x, y, z) is a CT number of a voxel in an image
  • ⁇ f(x, y, z) is an update value for this voxel
  • is a function
  • is a function.
  • may be a monotonically increasing function.
  • may be a piecewise continuous function.
  • the function ⁇ is based on CT number due to the observation that almost all CT imaging systems are calibrated based on water. That is, if the scanned object is made of water only, the reconstructed image will have a mean CT number of zero Hounsfield Units (HU) regardless of the energy level.
  • an image of an object made only of water will be the same if acquired at an energy level of 80 kVp or if it was acquired at an energy level of 140 kVp. Therefore, for the portion of the image with CT numbers close to zero, there is very little need for adjustment between the 140 kVp and 80 kVp energy levels. Instead, the adjustment is needed mainly for the material in the scanned object that is significantly different from water.
  • ⁇ p( ⁇ , ⁇ , ⁇ ) is the difference between one of the plurality of synthesized projections and the corresponding one of the second plurality of projections.
  • ⁇ , ⁇ , and a are respectively the fan-angle, projection angle, and cone angle corresponding to the image voxel (x, y, z).
  • the function ⁇ is a mapping function used to adjust the intensity of the updated image based on a path-length at the projection sample ( ⁇ , ⁇ , ⁇ ) For example, for a fixed projection difference, the function ⁇ may result in a smaller amount of adjustment per voxel for a long path-length through the object than for a shorter path-length through the object since the difference needs to be distributed over a larger number of voxels.
  • the following update function may be used to form an updated image at step 112 :
  • ⁇ ⁇ ⁇ f ⁇ ( x , y , z ) [ ⁇ f ⁇ ( x , y , z ) - t , 0 ⁇ m ] ⁇ p ⁇ ( ⁇ , ⁇ , ⁇ ) - p ′ ⁇ ( ⁇ , ⁇ , ⁇ ) ⁇ p ′ ⁇ ( ⁇ , ⁇ , ⁇ ) , s ⁇
  • f(x, y, z) ⁇ t,0 is a function that forces f(x, y, z) to be larger or equal to t,0;
  • f(x, y, z) is a CT number of a voxel in an image;
  • p( ⁇ , ⁇ , ⁇ ) is a projection acquired during step 104 passing through voxel (x, y, z);
  • p′( ⁇ , ⁇ , ⁇ ) is a corresponding synthesized projection;
  • t, m, and s are parameters.
  • One of the goals of the update function may be to generate an image of a patient or an object at the second energy level even though only an incomplete set of projections were acquired at the second energy level. For example, an image may be reconstructed from a first plurality of projections acquired at 140 kVp and it is desired to generate an image at 80 kVp. According to another embodiment, it may be possible to realize a reduction in x-ray dose by initially reconstructing an image at first energy level and then using the update function to generate an image at a second energy level that is higher than the first energy level.
  • modifying the image comprises applying an update function to the image that varies based on a CT number of a voxel in the image.
  • the update function may be weighted so that voxels with relatively high CT numbers, such as bone and iodine, are modified more than voxels with relatively low CT numbers, such as soft tissue.
  • the method 100 By making the update function dependent on the CT number of the voxels in the image, it is possible for the method 100 to converge on the desired image, such as the 80 kVp image, more quickly. Additionally, by using an update function that only makes small changes to portions of the image with CT numbers close to that of water, streaking artifacts caused by dense objects may be significantly reduced or eliminated. While the energy levels of 80 kVp and 140 kVp have been used as exemplary energy levels in this disclosure, it should be appreciated by those skilled in the art that additional embodiments my use energy levels significantly different than 80 kVp and/or 140 kVp.
  • the processor 22 determines if an additional iteration is required. According to an embodiment, the processor 22 determines if the difference between each of the plurality of synthesized projections and the corresponding one of the second plurality of projections is within an acceptable limit. According to another embodiment, the processor 22 determines if the method 100 has looped through steps 108 - 114 a predetermined number of times. Other embodiments may use additional methods of determining if additional iterations are required at step 114 . If an additional iteration is required at step 114 , the method 100 returns to step 108 where a plurality of synthesized projections are generated from the updated image of step 112 . Iteratively modifying the image by looping through steps 108 - 114 from three to five times has produced favorable results according to experimental data.
  • step 116 the updated image is displayed.
  • step 116 may be replaced by a step where the image is stored instead of being displayed.
  • the stored image may then be used in a subsequent step.
  • the updated image represents an image at the second energy level. It should be noted that by using the method 100 , it may be possible to generate an image of acceptable image quality at the second energy level even though too few projections were acquired at step 104 for an image to be reconstructed using a conventional reconstruction algorithm.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Pulmonology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

A method of generating a computed tomography image is disclosed herein. The method includes acquiring a first plurality of projections at a first energy level and acquiring a second plurality of projections at a second energy level. The method includes reconstructing an image from the first plurality of projections and generating a synthesized projection from the image that corresponds to one of the second plurality of projections. The method includes comparing the synthesized projection to the one of the second plurality of projections and modifying the image to form an updated image based on the comparing the synthesized projection to the one of the second plurality of projections. A computed tomography imaging system is also disclosed.

Description

    FIELD OF THE INVENTION
  • This disclosure relates generally to the field of medical imaging. In particular, the following techniques relate to multiple-energy computed tomography imaging systems and a method of generating a computed tomography image.
  • BACKGROUND OF THE INVENTION
  • Computed tomography (CT) imaging systems operate by projecting an x-ray beam from an x-ray source through an attenuating object, such as a patient. The attenuated x-ray beam is then detected by a detector assembly. Some materials, such as bone, are much more x-ray attenuating than other materials, such as soft tissue. Conventional third-generation CT imaging systems acquire attenuation data by projecting a polychromatic x-ray beam from the x-ray source. The polychromatic x-ray beam contains many different frequencies of x-rays and is typically centered around a particular energy level. With sufficient angular coverage around the patient, cross-sectional images can be formed revealing the inner structure of the scanned object. The images are typically displayed on a flat-screen monitor or cathode ray tube. A virtual 3-D image may also be produced based on data acquired during a computed tomography scan.
  • However, some materials share very similar x-ray attenuation characteristics at a particular energy level. For example, bone and iodine contrast agent have similar x-ray attenuation characteristics at some commonly used energy levels. As a result, it can be difficult to differentiate materials with similar x-ray attenuation characteristics at a particular energy level. Recently, multiple-energy CT imaging systems have been developed. By collecting x-ray attenuation data at more than one energy level, it is possible to gain additional insight into the nature of the scanned object.
  • In a conventional third-generation dual-energy CT imaging system, a processor may rapidly switch the output from a generator so that the input voltage to the x-ray source changes from projection to projection. For example, a typical dual-energy CT imaging system may rapidly alternate between acquiring a high-energy projection and a low-energy projection. In order for a conventional dual-energy CT imaging system to work well, the time spent transitioning between the high-energy level and the low-energy level must be short relative to the duration of each projection. Otherwise, the effective energy level during the high-energy projections will be less than the desired high-energy setting and the effective energy level during the low-energy projections will be greater than the desired low-energy setting. Given the fact that the sampling rate for a state-of-the-art third-generation CT imaging system is several kilohertz, the generator must transition between the high-energy level and the low-energy level in significantly less than a fraction of a millisecond. Enabling a generator to switch this quickly places an enormous demand on the generator hardware. Predictably, the demands on hardware will further increase as scan speed increases.
  • For these and other reasons, there is a need for a multiple-energy CT imaging method and system that significantly reduces the demand on generator hardware without negatively affecting scan time or image quality.
  • BRIEF DESCRIPTION OF THE INVENTION
  • The above-mentioned shortcomings, disadvantages and problems are addressed herein which will be understood by reading and understanding the following specification.
  • In an embodiment, a method of generating a computed tomography image includes acquiring a first plurality of projections at a first energy level and acquiring a second plurality of projections at a second energy level. The method includes reconstructing an image from the first plurality of projections. The method includes generating a synthesized projection from the image that corresponds to one of the second plurality of projections. The method includes comparing the synthesized projection to the one of the second plurality of projections. The method also includes modifying the image to form an updated image based on the comparing the synthesized projection to the one of the second plurality of projections.
  • In an embodiment, a method of generating a computed tomography image includes acquiring a first plurality of projections at a first energy level. The method includes acquiring a second plurality of projections at a second energy level. The method includes reconstruction an image from the first plurality of projections. The method includes generating a plurality of synthesized projections from the image, each of the plurality of synthesized projections corresponding to one of the second plurality of projections. The method includes comparing each of the plurality of synthesized projections to the corresponding one of the second plurality of projections. The method also includes modifying the image to form an updated image based on the comparing each of the plurality of synthesized projections to the corresponding one of the second plurality of projections.
  • In an embodiment, a method of generating a computed tomography image includes acquiring a first plurality of projections of an object at a first energy level. The method includes acquiring a second plurality of projections of the object at a second energy level. The method includes reconstructing an image of the object at the first energy level from the first plurality of projections. The method also includes modifying the image to generate an updated image of the object, wherein the updated image represents the object at the second energy level.
  • In an embodiment, a computed tomography imaging system includes a gantry, an x-ray source attached to the gantry, a table adapted to translate with respect to the gantry, and a processor in electrical communication with the gantry, the x-ray source, and the table. The processor is configured to acquire a first plurality of projections at a first energy level and to acquire a second plurality of projections at a second energy level. The processor is configured to reconstruct an image from the first plurality of projections. The processor is configured to generate a synthesized projection from the image that corresponds to one of the second plurality of projections. The processor is configured to determine a difference between the synthesized projection and the one of the second plurality of projections. The processor is also configured to generate an updated image by modifying the first image based on the difference between the synthesized projection and the one of the second plurality of projections.
  • Various other features, objects, and advantages of the invention will be made apparent to those skilled in the art from the accompanying drawings and detailed description thereof.
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • FIG. 1 is a schematic representation of a computed tomography imaging system in accordance with an embodiment; and
  • FIG. 2 is a flow chart illustrating a method of generating a computed tomography image in accordance with an embodiment.
  • DETAILED DESCRIPTION OF THE INVENTION
  • In the following detailed description, reference is made to the accompanying drawings that form a part hereof, and in which is shown by way of illustration specific embodiments that may be practiced. These embodiments are described in sufficient detail to enable those skilled in the art to practice the embodiments, and it is to be understood that other embodiments may be utilized and that logical, mechanical, electrical and other changes may be made without departing from the scope of the embodiments. The following detailed description is, therefore, not to be taken as limiting the scope of the invention.
  • Referring to FIG. 1, a schematic representation of a computed tomography (CT) imaging system 10 according to an embodiment is shown. The CT imaging system 10 includes a gantry support 12, a gantry 14, a table support 15, a table 16, an x-ray generator (not shown), an x-ray source 18, a detector assembly 20, and a processor 22. The gantry 14 is configured to rotate within the gantry support 12. The gantry 14 is adapted to retain the x-ray source 18 and the detector assembly 20. The x-ray generator is configured to deliver both a high-voltage input and a low-voltage input to the x-ray source 18. According to an embodiment, the high-voltage input may comprise a 140 kVp input while the low-voltage input may comprise an 80 kVp input. The x-ray source 18 receives either the high-voltage input or the low-voltage input from the generator and emits an x-ray beam. If the x-ray source 18 receives the high-voltage input, the x-ray beam will comprise a high-energy x-ray beam. If the x-ray source 18 receives the low-voltage input, the x-ray beam will comprise a low-energy x-ray beam. Both the high-energy x-ray beam and the low-energy x-ray beam comprise photons from a wide energy spectrum. One skilled in the art should appreciated that the term “high-energy x-ray beam” indicates that the highest energy x-ray photons in the high-energy x-ray beam are of high energy and that the term “low-energy x-ray beam” indicates that the highest energy x-ray photons in the low-energy x-ray beam are of low energy. Also, for the purposes of this disclosure, the terms “high-energy x-ray beam” and “low-energy x-ray beam” are relative to each other, and they should not be interpreted as indicating an absolute level of energy.
  • The x-ray source 18 is configured to emit either a high-energy x-ray beam or a low-energy x-ray beam through a patient 24 being examined. After passing through the patient 24, either the high-energy x-ray beam or the low-energy x-ray beam is received at the detector assembly 20. The detector assembly 20 comprises a plurality of detector elements (not shown). Each of the plurality of detector elements produces an electrical signal that varies based on the intensity of the x-ray beam received during a sampling interval. The table 16 is adapted to translate the patient 24 in a z-direction with respect to the gantry 14 as indicated by a coordinate axis 26. The processor 22 is configured to control the rotation of the gantry 14, the position of the table 16, and the activation of the x-ray source 18.
  • Having described the structure of the CT imaging system 10, a method 100 will be described hereinafter. FIG. 2 is a flow chart representing the method 100 in accordance with an embodiment. The individual blocks 102-116 represent steps that may be performed in accordance with the method 100. The steps 102-116 of the method 100 need not be performed in the order shown. The technical effect of the method 100 is the generation of computed tomography image.
  • Referring to both FIG. 1 and FIG. 2, the processor 22 may be configured to implement the steps illustrated in the method 100. At step 102, the processor 22 controls the generator (not shown) in a manner so that the x-ray source 18 emits an x-ray beam at a first energy level. According to an embodiment, the first energy level may be a relatively high energy level, such as 140 kVp. The processor 22 controls the generator, the x-ray source 18, the gantry 14, and the table 16 so that a first plurality of projections are acquired at the first energy level. The first plurality of projections may be acquired as part of a helical acquisition according to an embodiment and as part of an axial acquisition according to another embodiment.
  • For third-generation CT systems, a projection may be defined to include data obtained at generally the same time from some or all of the detector elements for a given gantry angle and location along the z-axis. In a third-generation CT imaging system a projection may comprise the output from all of the detector elements at generally the same time. The number of projections that can be acquired depends upon the speed of gantry rotation, the type of hardware used in the detector assembly 20, and the number of detector elements among other factors. According to an embodiment, 984 projections may comprise a complete set of projections for one gantry revolution.
  • At step 104, the processor 22 controls the generator (not shown) in a manner so that the x-ray source 18 emits an x-ray beam at a second energy level. According to an embodiment, the second energy level may be a relatively low energy level, such as 80 kVp. The processor 22 controls the generator, the x-ray source 18, the gantry 14, and the table 16 so that a second plurality of projections are obtained at the second energy level. The second plurality of projections may be acquired as part of a helical acquisition according to an embodiment or as part of an axial acquisition according to another embodiment. Steps 102 and 104 may be performed simultaneously, sequentially, or a combination of simultaneously and sequentially. According to an embodiment where steps 102 and 104 are performed simultaneously, the first plurality of projections are acquired during generally the same period of time that the second plurality of projections are acquired. For example, the method 100 may alternate between acquiring one or more projections at the first energy level and acquiring one or more projections at the second energy level. According to this example, the first plurality of projections would comprise all of the projections acquired at the first energy level and the second plurality of projection would comprise all of the projections acquired at the second energy level.
  • It is well-known by those skilled in the art that an image produced by an x-ray beam at a higher energy level may generally produce less noise than an image produced by an x-ray beam at a lower energy level if all other imaging parameters are held constant. Therefore, according to an embodiment, it may be preferable to collect more projections at the higher energy level than at the lower energy level. According to an embodiment, the processor 22 may control the generator (not shown) and x-ray source 18 in a manner so that one lower energy projection is acquired for every 24 projections. According to an exemplary embodiment, projections 1 through 23 may be acquired at a first energy level, such as 140 kVp, and projection 24 may be acquired at a second energy level, such as 80 kVp. All 984 projections of the exemplary embodiment may be acquired by alternating between acquiring 23 projections at 140 kVp and acquiring one projection at 80 kVp. Therefore, according to the exemplary embodiment, the result will be 943 projections at 140 kVp and 41 projections at 80 kVp. It should be understood that other embodiments may collect the first plurality of projections and the second plurality of projection in an entirely different manner.
  • In accordance with the previously described exemplary embodiment, the first plurality of projections comprises 943 projection acquired at 140 kVp and the second plurality of projections comprises 41 projections acquired at 80 kVp. Several or all of the 41 projections acquired at 80 kVp may be separated by a common angle. For example, according to an embodiment, there may be a common angle of approximately 8.8 degrees separating each of the 41 projections acquired at 80 kVp. According to another embodiment, some or all of the 41 projections acquired at 80 kVp may be separated by an amount that varies. For example, some of the projections may be spaced apart by 2 degrees or less and others of the projections may be spaced apart by more than 10 degrees. The amount of angular separation between projections at a particular energy level may be determined by many factors including hardware limitations of the CT imaging system 10 (shown in FIG. 1).
  • According to an embodiment, a half-scan may be used during the acquisition of the first plurality of projections and the second plurality of projections to improve the temporal resolution of the CT imaging system. Half-scans are well-known by those skilled in the art. In accordance with an embodiment using a half-scan, a total of 616 projections are acquired at 140 kVp and 27 projections are acquired at 80 kVp. Other embodiments may acquire a different ratio of projections acquired at a first energy level to projections acquired at a second energy level.
  • In accordance with other embodiments, projections may be acquired at three or more different energy levels. For example, according to an embodiment, a first plurality of projections may comprise 902 projections at 140 kVp, a second plurality of projections may comprise 41 projections at 100 kVp, and a third plurality of projections may comprise 41 projections at 80 kVp. According to an embodiment, the first plurality of projections, the second plurality of projections, and the third plurality of projections may be acquired in an interleaved manner. In other words, the processor 22 (shown in FIG. 1) may rotate between acquiring one or more projections at each of the three energy levels.
  • In accordance with another embodiment, the generator (not shown) may utilize a slower transition when switching between energy levels. According to one embodiment, 943 projections are acquired at 140 kVp and 41 projections are acquired at 80 kVp. Each of the 41 projections acquired at 80 kVp may be evenly spaced apart. The projections that are immediately before or after each of the projections acquired at 80 kVp projection may be collected while the generator is in a transition zone. For the purposes of this disclosure, the term “transition zone” includes a time when the generator is between two of the desired voltages. In accordance with an exemplary embodiment, the spectrum of the x-ray beam may be somewhere between the desired levels of 80 kVp and 140 kVp during each of the transition zones. The demand on the generator is significantly reduced by acquiring projections during the transition zones instead of requiring that the generator is at exactly 140 kVp immediately before and after each 80 kVp projection.
  • According to another embodiment, the x-ray tube current used when acquiring the projections at the first energy level may be substantially different than the x-ray tube current used when acquiring the projections at the second energy level. For example, according to an embodiment designed to reduce x-ray dose, the processor 22 may utilize a reduced x-ray tube current level while acquiring the projections at the higher energy level. The reduced x-ray tube current level is reduced in comparison to the x-ray tube current level used while acquiring the projections at the lower energy level. According to another embodiment, the duration of a projection acquired at first energy level may be significantly different than the duration of a view acquired at second energy level.
  • Still referring to FIG. 2, at step 106 an image is reconstructed from the first plurality of projections. According to the previously described exemplary embodiment, the image is reconstructed from the 943 projections acquired at 140 kVp. Since 943 projections at 140 kVp represents an almost complete set of projections, the first image may be reconstructed by using a conventional reconstruction algorithm, such as a filtered back-projection algorithm. If a filtered back-projection algorithm is used to reconstruct the first image, it may be necessary to estimate any missing projections. For the purposes of this disclosure, the term “image” includes a reconstructed volume from one or more slices of computed tomography data. For example, for a multi-slice CT imaging system, an image may comprise data from a plurality of slices. According to the previously described exemplary embodiment, 943 projections were acquired 140 kVp and 41 of the projections were acquired at 80 kVp. This embodiment is missing 41 projections from the 140 kVp set of projections. The missing projections may be estimated by performing a linear interpolation or a high-order interpolation based on adjacent projections that were acquired at 140 kVp. It should be appreciated that the first image may be reconstructed through other techniques according to additional embodiments. For example, an iterative reconstruction technique may be used to generate the first image.
  • At step 108, a plurality of synthesized projections are generated. According to an embodiment, the plurality of synthesized projections are based on the image the first time the method 100 reaches step 108. According to an embodiment with multiple iterations, a plurality of synthesized projections may be generated based on the updated image in subsequent iterations as will be described hereinafter. For purposes of this disclosure, the term “synthesized projection” is defined to include a forward projection from an image. Normally, a CT imaging system collects a plurality of projections and uses the projections to reconstruct an image. However, in the case of a synthesized projection, an image is used as the input and the synthesized projection represents what the projection would most likely comprise at a particular location and gantry angle. Synthesized projections are well-known by those skilled in the art. Each of the synthesized projections generated at step 108 may correspond to one of the second plurality of projections. For the purposes of this disclosure, a projection and a synthesized projection are considered to correspond if they both represent generally the same location and gantry angle of the same object. Additional normalization steps may be used to adjust the scale of the synthesized projection. For example, according to another embodiment, one of the plurality of synthesized projections generated at step 108 may represent the same location and gantry angle as one of the first plurality of projections acquired at step 102. Then, a scaling factor or a plurality of polynomial parameters may be determined based on a comparison of the synthesized projection and the corresponding one of the first plurality of projections. The scaling factor or the plurality of polynomial parameters may be used to adjust each of the synthesized projections generated at step 108.
  • According to an embodiment, each of the plurality of synthesized projections generated at step 108 corresponds to one of the second plurality of projections. Therefore, each one of the synthesized projections is generated at the same gantry angle and location in the z-direction as one of the second plurality of projections. According to the previously described exemplary embodiment, the second plurality of projections comprises 41 projections that were acquired at 80 kVp. Each of the plurality of synthesized projections corresponds to one of the 41 projections that were acquired at 80 kVp. In other words, the plurality of synthesized projections represent the same locations and angles as the 41 projections that were acquired at 80 kVp during step 104. It is not necessary for the plurality of synthesized projections and the second plurality of projections to be of the same number in all embodiments.
  • According to another embodiment, there may be only four synthesized projections generated during step 108. The four synthesized projections may be separated by approximately 45 degrees from each other. It should be appreciated that each of the four synthesized projections still corresponds to one of the second plurality of projections acquired during step 104. According to an additional embodiment, step 108 may be replaced with a step where only one synthetic projection is generated for each iteration.
  • At step 110, each of the plurality of synthesized projections is compared to the corresponding one of the second plurality of projections. A difference and/or a ratio may be calculated between each synthesized projection and the corresponding one of the second plurality of projections. For example, if the synthesized projections were generated from an image acquired at 140 kVp and the second plurality of projections were acquired at 80 kVp, any difference between a synthesized projection and the corresponding one of the second plurality of projections may be assumed to be from the mean energy difference between the x-ray beam produced at 140 kVp and the x-ray beam produced at 80 kVp.
  • At step 112, the first image is modified based on the comparison performed at step 110 to form an updated image. According to an embodiment, a difference between each of the synthesized projections generated at step 108 and the corresponding projection acquired as part of the second plurality of projections during step 104 is distributed according to the following equation:

  • Δf(x, y, z)=ξ[f(x, y, z)]×η[Δp(γ, β, α)]
  • where f(x, y, z) is a CT number of a voxel in an image, Δf(x, y, z) is an update value for this voxel, ξ is a function, and η is a function. According to an embodiment, ξ may be a monotonically increasing function. In another embodiment, ξ may be a piecewise continuous function. The function ξ is based on CT number due to the observation that almost all CT imaging systems are calibrated based on water. That is, if the scanned object is made of water only, the reconstructed image will have a mean CT number of zero Hounsfield Units (HU) regardless of the energy level. In other words, an image of an object made only of water will be the same if acquired at an energy level of 80 kVp or if it was acquired at an energy level of 140 kVp. Therefore, for the portion of the image with CT numbers close to zero, there is very little need for adjustment between the 140 kVp and 80 kVp energy levels. Instead, the adjustment is needed mainly for the material in the scanned object that is significantly different from water. Δp(γ, β, α) is the difference between one of the plurality of synthesized projections and the corresponding one of the second plurality of projections. According to an embodiment, γ, β, and a are respectively the fan-angle, projection angle, and cone angle corresponding to the image voxel (x, y, z). The function η is a mapping function used to adjust the intensity of the updated image based on a path-length at the projection sample (γ, β, α) For example, for a fixed projection difference, the function η may result in a smaller amount of adjustment per voxel for a long path-length through the object than for a shorter path-length through the object since the difference needs to be distributed over a larger number of voxels.
  • According to an embodiment, the following update function may be used to form an updated image at step 112:
  • Δ f ( x , y , z ) = [ f ( x , y , z ) - t , 0 m ] p ( γ , β , α ) - p ( γ , β , α ) p ( γ , β , α ) , s
  • where
    Figure US20100215233A1-20100826-P00001
    f(x, y, z)−t,0
    Figure US20100215233A1-20100826-P00002
    is a function that forces f(x, y, z) to be larger or equal to t,0; f(x, y, z) is a CT number of a voxel in an image; p(γ, β, α) is a projection acquired during step 104 passing through voxel (x, y, z); p′(γ, β, α) is a corresponding synthesized projection;
    Figure US20100215233A1-20100826-P00001
    p′(γ, β, α), s
    Figure US20100215233A1-20100826-P00002
    forces p′(γ, β, α) to be equal or larger than s; and t, m, and s are parameters.
  • One of the goals of the update function may be to generate an image of a patient or an object at the second energy level even though only an incomplete set of projections were acquired at the second energy level. For example, an image may be reconstructed from a first plurality of projections acquired at 140 kVp and it is desired to generate an image at 80 kVp. According to another embodiment, it may be possible to realize a reduction in x-ray dose by initially reconstructing an image at first energy level and then using the update function to generate an image at a second energy level that is higher than the first energy level.
  • Still referring to step 112, according to an embodiment, modifying the image comprises applying an update function to the image that varies based on a CT number of a voxel in the image. As was described previously, experimental results have shown that there is a relatively small difference between the soft-tissue or water portion of an image at 80 kVp and the soft-tissue or water portion of an image at 140 kVp. Conversely, it is known that there is a relatively large difference between bone at 80 kVp and bone at 140 kVp. Therefore, if the first image, reconstructed from data acquired at 140 kVp, is used as an initial assumption, it should not be necessary to significantly change the voxel values in regions with CT numbers close to that of water when modifying the 140 kVp image to form an image at 80 kVp. On the other hand, there should be more significant changes to areas with relatively high CT numbers in the image. According to an embodiment, the update function may be weighted so that voxels with relatively high CT numbers, such as bone and iodine, are modified more than voxels with relatively low CT numbers, such as soft tissue. By making the update function dependent on the CT number of the voxels in the image, it is possible for the method 100 to converge on the desired image, such as the 80 kVp image, more quickly. Additionally, by using an update function that only makes small changes to portions of the image with CT numbers close to that of water, streaking artifacts caused by dense objects may be significantly reduced or eliminated. While the energy levels of 80 kVp and 140 kVp have been used as exemplary energy levels in this disclosure, it should be appreciated by those skilled in the art that additional embodiments my use energy levels significantly different than 80 kVp and/or 140 kVp.
  • At step 114, the processor 22 (shown in FIG. 1) determines if an additional iteration is required. According to an embodiment, the processor 22 determines if the difference between each of the plurality of synthesized projections and the corresponding one of the second plurality of projections is within an acceptable limit. According to another embodiment, the processor 22 determines if the method 100 has looped through steps 108-114 a predetermined number of times. Other embodiments may use additional methods of determining if additional iterations are required at step 114. If an additional iteration is required at step 114, the method 100 returns to step 108 where a plurality of synthesized projections are generated from the updated image of step 112. Iteratively modifying the image by looping through steps 108-114 from three to five times has produced favorable results according to experimental data.
  • If no additional iterations are necessary at step 114, the method 100 proceeds to step 116, where the updated image is displayed. According to other embodiments, step 116 may be replaced by a step where the image is stored instead of being displayed. The stored image may then be used in a subsequent step. According to an embodiment, the updated image represents an image at the second energy level. It should be noted that by using the method 100, it may be possible to generate an image of acceptable image quality at the second energy level even though too few projections were acquired at step 104 for an image to be reconstructed using a conventional reconstruction algorithm.
  • This written description uses examples to disclose the invention, including the best mode, and also to enable any person skilled in the art to practice the invention, including making and using any devices or systems and performing any incorporated methods. The patentable scope of the invention is defined by the claims, and may include other examples that occur to those skilled in the art. Such other examples are intended to be within the scope of the claims if they have structural elements that do not differ from the literal language of the claims, or if they include equivalent structural elements with insubstantial differences from the literal language of the claims.

Claims (24)

1. A method of generating a computed tomography image comprising:
acquiring a first plurality of projections at a first energy level;
acquiring a second plurality of projections at a second energy level;
reconstructing an image from the first plurality of projections;
generating a synthesized projection from the image that corresponds to one of the second plurality of projections;
comparing the synthesized projection to said one of the second plurality of projections; and
modifying the image to form an updated image based on said comparing the synthesized projection to said one of the second plurality of projections.
2. The method of claim 1, wherein the first energy level comprises a higher energy level than the second energy level.
3. The method of claim 1, wherein the first energy level comprises a lower energy level than the second energy level.
4. The method of claim 1, wherein said comparing the synthesized projection to said one of the second plurality of projections comprises calculating a difference between the synthesized projection and said one of the second plurality of projections.
5. The method of claim 1, wherein said modifying the image comprises applying an update function that varies based on a CT number.
6. The method of claim 1, further comprising displaying the updated image.
7. The method of claim 1, wherein the updated image represents the second energy level.
8. The method of claim 1, wherein the first plurality of projections contains more than 5 times the number of projections as the second plurality of projections.
9. The method of claim 1, wherein the second plurality of projections comprises two or more projections separated by a common angle.
10. The method of claim 1, wherein the second plurality of projections comprises 4 projections separated by 45 degrees.
11. A method of generating a computed tomography image comprising:
acquiring a first plurality of projections at a first energy level;
acquiring a second plurality of projections at a second energy level;
reconstructing an image from the first plurality of projections;
generating a plurality of synthesized projections from the image, each of the plurality of synthesized projections corresponding to one of the second plurality of projections;
comparing each of the plurality of synthesized projections to the corresponding one of the second plurality of projections; and
modifying the image to form an updated image based on said comparing each of the plurality of synthesized projections to the corresponding one of the second plurality of projections.
12. The method of claim 11, wherein said modifying the image comprises applying an update function that varies based on a CT number.
13. The method of claim 11, wherein said comparing comprises calculating a difference between each of the plurality of synthesized projections and the corresponding one of the second plurality of projections.
14. A method of generating a computed tomography image comprising:
acquiring a first plurality of projections of an object at a first energy level;
acquiring a second plurality of projections of the object at a second energy level;
reconstructing an image of the object at the first energy level from the first plurality of projections;
modifying the image to generate an updated image of the object, wherein the updated image represents the object at the second energy level.
15. The method of claim 14, wherein said modifying the image comprises iteratively modifying the image.
16. The method of claim 14, wherein iteratively modifying the image comprises four or more iterations.
17. The method of claim 14, further comprising acquiring a third plurality of projections at a third energy level.
18. The method of claim 17, further comprising modifying the image to generate a second updated image, wherein the second updated image represents the third energy level.
19. A computed tomography imaging system comprising;
a gantry;
an x-ray source attached to the gantry;
a table adapted to translate with respect to the gantry; and
a processor in electrical communication with the gantry, the x-ray source, and the table, wherein the processor is configured to:
acquire a first plurality of projections at a first energy level;
acquire a second plurality of projections at a second energy level;
reconstruct an image from the first plurality of projections;
generate a synthesized projection from the image that corresponds to one of the second plurality of projections;
determine a difference between the synthesized projection and said one of the second plurality of projections; and
generate an updated image by modifying the first image based on the difference between the synthesized projection and said one of the second plurality of projections.
20. The computed tomography imaging system of claim 19, wherein the processor is further configured to display the updated image.
21. The computed tomography imaging system of claim 19, wherein the first energy level exceeds the second energy level.
22. The computed tomography imaging system of claim 19, wherein the second energy level exceeds the first energy level.
23. The computed tomography imaging system of claim 19, wherein the updated image represents the second energy level.
24. The computed tomography imaging system of claim 19, wherein the processor is further configured to apply an update function that varies based on a CT number of a voxel in the image.
US12/393,816 2009-02-26 2009-02-26 Method and system for generating a computed tomography image Active 2030-06-26 US8050479B2 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
US12/393,816 US8050479B2 (en) 2009-02-26 2009-02-26 Method and system for generating a computed tomography image
NL2004249A NL2004249C2 (en) 2009-02-26 2010-02-16 METHOD AND SYSTEM FOR GENERATING A COMPUTED TOMOGRAPHY IMAGE.
JP2010033080A JP5491230B2 (en) 2009-02-26 2010-02-18 Computerized tomography image forming system.
CN2010101323581A CN101912272B (en) 2009-02-26 2010-02-26 Method and system for generating a computed tomography image

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US12/393,816 US8050479B2 (en) 2009-02-26 2009-02-26 Method and system for generating a computed tomography image

Publications (2)

Publication Number Publication Date
US20100215233A1 true US20100215233A1 (en) 2010-08-26
US8050479B2 US8050479B2 (en) 2011-11-01

Family

ID=42313119

Family Applications (1)

Application Number Title Priority Date Filing Date
US12/393,816 Active 2030-06-26 US8050479B2 (en) 2009-02-26 2009-02-26 Method and system for generating a computed tomography image

Country Status (4)

Country Link
US (1) US8050479B2 (en)
JP (1) JP5491230B2 (en)
CN (1) CN101912272B (en)
NL (1) NL2004249C2 (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140050302A1 (en) * 2012-08-14 2014-02-20 Siemens Aktiengesellschaft Method for the determination of dual energy image data records and x-ray facility
US9247919B2 (en) 2011-02-01 2016-02-02 Koninklijke Philips N.V. Method and system for dual energy CT image reconstruction
US20160117850A1 (en) * 2014-10-24 2016-04-28 General Electric Company Systems and methods for metal artifact reduction
WO2016093917A1 (en) * 2014-12-11 2016-06-16 General Electric Company Methods and systems for spectral ct imaging
US20170365076A1 (en) * 2016-06-21 2017-12-21 Carestream Health, Inc. Virtual projection image method
US20180000430A1 (en) * 2016-06-29 2018-01-04 Siemens Healthcare Gmbh Generating an image sequence
US20180182129A1 (en) * 2015-09-09 2018-06-28 Tsinghua University Spectral ct image reconstructing method and spectral ct imaging system
US11026642B2 (en) * 2019-03-29 2021-06-08 Canon Medical Systems Corporation Apparatuses and a method for artifact reduction in medical images using a neural network

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9517042B2 (en) * 2014-09-12 2016-12-13 General Electric Company Systems and methods for imaging phase selection for computed tomography imaging
US10165996B2 (en) * 2015-09-30 2019-01-01 General Electric Company Systems and methods for dual-energy computed tomography imaging
JP7233911B2 (en) * 2018-12-17 2023-03-07 キヤノンメディカルシステムズ株式会社 X-ray CT system and processing program

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030031289A1 (en) * 2001-07-18 2003-02-13 Jiang Hsieh Methods and apparatus for FOV-dependent aliasing artifact reduction
US20040101086A1 (en) * 2002-11-27 2004-05-27 Sabol John Michael Method and apparatus for quantifying tissue fat content
US6816567B2 (en) * 2002-07-15 2004-11-09 Ge Medical System Global Technology Company, Llc System and method for acquiring x-ray data
US20050105679A1 (en) * 2003-02-12 2005-05-19 Tao Wu Tomosynthesis imaging system and method
US20080310598A1 (en) * 2007-03-14 2008-12-18 Tsinghua University Method for calibrating dual-energy CT system and method of image reconstruction
US20090147919A1 (en) * 2005-05-18 2009-06-11 Taiga Goto Radiographing apparatus and image processing program
US20100128948A1 (en) * 2008-11-26 2010-05-27 General Electric Company Forward projection for the generation of computed tomography images at arbitrary spectra
US20100183117A1 (en) * 2007-07-19 2010-07-22 Hitachi Medical Corporation X-ray generating apparatus and x-ray ct apparatus using the same

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7272429B2 (en) * 2002-11-27 2007-09-18 Ge Medical Systems Global Technology Company, Llc Methods and apparatus for facilitating a reduction in artifacts

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030031289A1 (en) * 2001-07-18 2003-02-13 Jiang Hsieh Methods and apparatus for FOV-dependent aliasing artifact reduction
US6816567B2 (en) * 2002-07-15 2004-11-09 Ge Medical System Global Technology Company, Llc System and method for acquiring x-ray data
US20040101086A1 (en) * 2002-11-27 2004-05-27 Sabol John Michael Method and apparatus for quantifying tissue fat content
US20050105679A1 (en) * 2003-02-12 2005-05-19 Tao Wu Tomosynthesis imaging system and method
US20090147919A1 (en) * 2005-05-18 2009-06-11 Taiga Goto Radiographing apparatus and image processing program
US20080310598A1 (en) * 2007-03-14 2008-12-18 Tsinghua University Method for calibrating dual-energy CT system and method of image reconstruction
US20100183117A1 (en) * 2007-07-19 2010-07-22 Hitachi Medical Corporation X-ray generating apparatus and x-ray ct apparatus using the same
US20100128948A1 (en) * 2008-11-26 2010-05-27 General Electric Company Forward projection for the generation of computed tomography images at arbitrary spectra

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9247919B2 (en) 2011-02-01 2016-02-02 Koninklijke Philips N.V. Method and system for dual energy CT image reconstruction
US20140050302A1 (en) * 2012-08-14 2014-02-20 Siemens Aktiengesellschaft Method for the determination of dual energy image data records and x-ray facility
US9830685B2 (en) * 2012-08-14 2017-11-28 Siemens Aktiengesellschaft Method for the determination of dual energy image data records and X-ray facility
US9655580B2 (en) * 2014-10-24 2017-05-23 General Electric Company Systems and methods for metal artifact reduction
US20160117850A1 (en) * 2014-10-24 2016-04-28 General Electric Company Systems and methods for metal artifact reduction
WO2016093917A1 (en) * 2014-12-11 2016-06-16 General Electric Company Methods and systems for spectral ct imaging
US9585626B2 (en) 2014-12-11 2017-03-07 General Electric Company Methods and systems for spectral CT imaging
US20180182129A1 (en) * 2015-09-09 2018-06-28 Tsinghua University Spectral ct image reconstructing method and spectral ct imaging system
US10489939B2 (en) * 2015-09-09 2019-11-26 Tsinghua University Spectral CT image reconstructing method and spectral CT imaging system
US20170365076A1 (en) * 2016-06-21 2017-12-21 Carestream Health, Inc. Virtual projection image method
US9965875B2 (en) * 2016-06-21 2018-05-08 Carestream Health, Inc. Virtual projection image method
US20180000430A1 (en) * 2016-06-29 2018-01-04 Siemens Healthcare Gmbh Generating an image sequence
US10751007B2 (en) * 2016-06-29 2020-08-25 Siemens Healthcare Gmbh Generating an image sequence
US11026642B2 (en) * 2019-03-29 2021-06-08 Canon Medical Systems Corporation Apparatuses and a method for artifact reduction in medical images using a neural network

Also Published As

Publication number Publication date
NL2004249C2 (en) 2010-11-30
US8050479B2 (en) 2011-11-01
CN101912272A (en) 2010-12-15
CN101912272B (en) 2013-04-24
JP2010194311A (en) 2010-09-09
JP5491230B2 (en) 2014-05-14
NL2004249A (en) 2010-08-30

Similar Documents

Publication Publication Date Title
US8050479B2 (en) Method and system for generating a computed tomography image
US7272429B2 (en) Methods and apparatus for facilitating a reduction in artifacts
US8705822B2 (en) Method for creating images indicating material decomposition in dual energy, dual source helical computed tomography
Buzug Computed tomography
JP4402435B2 (en) Method and apparatus for visualization of soft tissue space
JP5469293B2 (en) Method for image reconstruction using a hybrid computed tomography detector
US6891918B2 (en) Methods and apparatus for acquiring perfusion data
JP5582514B2 (en) X-ray CT system
US9159462B2 (en) Detection apparatus
JP6242631B2 (en) Medical image processing apparatus and X-ray computed tomography apparatus
JP6492005B2 (en) X-ray CT apparatus, reconstruction calculation apparatus, and reconstruction calculation method
US10789738B2 (en) Method and apparatus to reduce artifacts in a computed-tomography (CT) image by iterative reconstruction (IR) using a cost function with a de-emphasis operator
JP2002017716A (en) Method and device for multi-slice ct scanning for region of interest
US11419566B2 (en) Systems and methods for improving image quality with three-dimensional scout
JP7199399B2 (en) Systems and methods for automatic tube potential selection in dual-energy imaging
US20160106386A1 (en) Methods and systems for task-based data generation and weighting for ct spectral imaging
US10383589B2 (en) Direct monochromatic image generation for spectral computed tomography
US20050018889A1 (en) Systems and methods for filtering images
US11361415B2 (en) Material-selective adaptive blending of volumeiric image data
US20120177173A1 (en) Method and apparatus for reducing imaging artifacts
EP3404618B1 (en) Poly-energetic reconstruction method for metal artifacts reduction
JP2023039438A (en) Image generation device, x-ray ct apparatus and image generation method

Legal Events

Date Code Title Description
AS Assignment

Owner name: GENERAL ELECTRIC COMPANY, NEW YORK

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:HSIEH, JIANG;REEL/FRAME:022345/0633

Effective date: 20090220

STCF Information on status: patent grant

Free format text: PATENTED CASE

FPAY Fee payment

Year of fee payment: 4

MAFP Maintenance fee payment

Free format text: PAYMENT OF MAINTENANCE FEE, 8TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1552); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY

Year of fee payment: 8

MAFP Maintenance fee payment

Free format text: PAYMENT OF MAINTENANCE FEE, 12TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1553); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY

Year of fee payment: 12