US20020012417A1 - Method for reducing line artifacts in a CT image and device for implementing the method - Google Patents
Method for reducing line artifacts in a CT image and device for implementing the method Download PDFInfo
- Publication number
- US20020012417A1 US20020012417A1 US09/797,330 US79733001A US2002012417A1 US 20020012417 A1 US20020012417 A1 US 20020012417A1 US 79733001 A US79733001 A US 79733001A US 2002012417 A1 US2002012417 A1 US 2002012417A1
- Authority
- US
- United States
- Prior art keywords
- image
- filtering
- picture elements
- resulting
- tangents
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/40—Imaging
- G01N2223/419—Imaging computed tomograph
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/612—Specific applications or type of materials biological material
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
Definitions
- the present invention is directed to a method and apparatus for reducing line artifacts in a CT image, wherein the CT image is produced by scanning with an X-ray source which is rotatable around an examination subject, with X-rays from the X-ray source, after being attenuated by an examination subject, being incident on a detector system.
- CT devices which have an X-ray source, e.g. an X-ray tube, which direct a collimated, pyramid-shaped X-ray bundle through the examination subject, e.g. a patient, onto a detector system that is composed of a number of detector channels.
- Each detector channel has at least one detector element and one associated electronic element for reading out and amplifying the signal that is generated in the detector element as a result of the incident radiation.
- a number of detector elements can be allocated to one electronic element.
- the X-ray source and, depending on the construction of the CT device, the detector system as well are attached to a gantry that rotates around the examination subject.
- a support device for the examination subject can be displaced along the system axis relative to the gantry.
- the position along the system axis at which the X-ray bundle penetrates the examination subject, and the angle, under which the X-ray bundle penetrates the examination subject, are continuously modified as a result of the displacement and the rotation of the gantry.
- Each detector element of the detector system struck by the radiation produces a signal representing a measure of the overall transparency of the examination subject for the radiation proceeding from the X-ray source to the detector system.
- the set of output signals of the detector elements of the detector system, which is acquired for a specific position of the X-ray source, is referred to as a projection.
- a scan is composed of a set of projections, which are acquired at different positions of the gantry and/or at different positions of the support device.
- the CT device picks up a number of projections during a scan in order to be able to construct a two-dimensional tomogram of a slice of the examination subject.
- a number of slices can be picked up at the same time by a detector system that is composed of an array having a number of rows and columns of detector elements.
- Such planar-like detector systems frequently contain detector channels which do not supply proper data. It may be that detector system contains faulty detector channels already after the production process, for example as a result of defects in fabrication caused by the high integration density of the electronic elements. Defects of individual detector channels also may arise during the operation of the CT device. Such defects cause circular structures in the acquired CT images, these circular structures being referred to as circle artifacts.
- a disadvantage of such known methods is that they insufficiently eliminate artifacts which arise in a CT device having a defective detector channel.
- An object of the present invention is to provide a method for reducing artifacts in a CT image, so that the obtainable image quality is improved in a CT device having at least one defective detector channel. It is also an object of the invention to provide a CT device for implementing the method.
- the above object is achieved in accordance with the principles of the present invention in a method for reducing line artifacts in a CT image D 1 , as well as in an apparatus for implementing the method, wherein the image has been subjected to interpolated filtering for preventing or correcting faulty values of picture elements represented in a circle K 1 , and wherein the following steps are implemented.
- the picture elements of the CT image D 1 are subjected to a median filtering, orthogonal to the straight line extending through the respective picture element and the center of the circle K 1 , for producing an image M 1 .
- a difference value image F 1 is generated by subtracting the image M 1 from the CT image D 1 .
- Two resulting images G 1 1 and G 2 1 are produced by filtering the picture elements of the difference value image F 1 in the respective directions of tangents t 1 1 and t 2 1 to the circle K 1 extending through the respective picture element. Filtering is conducted along t 1 1 to produce the resulting image G 1 1 and is conducted along t 2 1 for producing the resulting image G 2 1 .
- a correction image D 2 is then obtained by subtracting both of the resulting images G 1 1 and G 2 1 from the image D 1 .
- Defective detector channels of a detector system lead to faulty values for picture elements, which appear as a circle in a CT image acquired by the detector system. Such image errors therefore are referred to as circle artifacts.
- Defective detector channels are not only ones that fail to supply an output signal as a result of the defect, but also are channels with a measuring accuracy that exceeds a specific tolerance value.
- Various methods are known for correcting or preventing circle artifacts, as noted above. These methods are carried out on the measurement data or on the image data and are primarily based on interpolated filtering, and achieve a noticeable weakening of the circle artifacts in the acquired CT images.
- This width should be selected wider than the expected line width of the line artifacts.
- a median filter of the width 5 has proven to be beneficial.
- the sampling distance A 1 must be selected dependent on the convolution kernel used for the reconstruction. Ideally, the image M 1 no longer contains line artifacts.
- a filtering is carried out in each picture element of the difference value image F 1 in the respective directions of the tangents t 1 1 and t 2 1 to the circle K 1 extending through the picture element in question, i.e., the filtering is carried out along t 1 1 for the resulting image G 1 1 and the filtering is carried out along t 2 1 for the resulting image G 2 1 .
- This step is to eliminate, for the most part, the image pixel noise difference value in the image F 1 for the most part.
- the correction image D 2 serves as basis image for correcting the line artifacts caused by other detector channels if the detector system contains further defective detector channels.
- the correction image D 1 serves as a basis image for correcting line artifacts caused by the i-th defective detector channel. Since the image M 1 has been calculated in relation to the center of the circle that is valid for all circles k i , the difference value image F 1 can still be used in the following.
- a filtering is carried out in each picture element of the difference value image F 1 situated outside of the circle K i in the direction of the tangents D 1 i and D 2 i to the circle K i , i.e., the filtering is carried out along D 1 i for the resulting image G 1 i and the filtering is carried out along D 2 i for the resulting image G 2 i .
- the above-described method is relatively time-consuming as a result of the serial processing of the line artifacts caused by the defective detector channels.
- Parallel processing of the image errors caused by the different defective detector channels is preferable for processing the data faster.
- a further embodiment of the inventive method proceeds as described above until the calculation of the difference value image F 1 .
- the resulting images G 1 i and G 2 i are calculated in this manner for all defective detector channels.
- different correction images are not consecutively determined, which respectively serve as an image for calculating the next correction images, but instead all resulting images G 1 i and G 2 1 are subtracted from D 1 .
- the correction image D 2 results.
- the defective detector channels should have a minimum distance (spacing) of ten detector channels as a condition for this parallel processing.
- a median filtering is carried out in each picture element of the CT image D i , orthogonal to the tangents to the circle K i extending through the respective picture element, for producing an image M i ,
- a filtering is carried out in each picture element of the difference value image F 1 situated outside of the circle K 1 in the direction of the tangents t 1 i and t 2 i to the circle K i extending through the respective picture element in order to produce two resulting images G 1 i and G 2 i , with the filtering being carried out along t 1 i for the resulting image G 1 i and the filtering being carried out along t 2 i for the resulting image G 2 i .
- the correction image D i ⁇ 1 D i ⁇ D i ⁇ G 1 i ⁇ G 2 i is determined by subtracting the resulting images G 1 i and G 2 i from D i .
- a running averaging serves as the filtering for producing the resulting images G 1 i and G 2 i .
- a sum operator is implemented as the filtering for producing the resulting images G 1 i and G 2 i .
- a combination of these two versions also can be employed, wherein a running averaging and a sum operator are simultaneously employed for filtering.
- the values of the CT image D 1 are preferably limited to a range 1000 HU (H 2 O) ⁇ . This limitation is expedient for examining soft-tissue parts, in which disturbances caused by defective detector channels are particularly apparent in the image.
- the intensity of the noise signal in the resulting images G 1 i and G 2 i is dependent on the filter width of the running averaging or of the sum operator.
- This noise signal influences the correction images D 1+1 and therefore can lead to undesired noise structures in the resulting image in the case of a number of defective detector channels.
- a high filter width is required as a result.
- the intensity of a line artifact varies, however, so that a limitation of the filter width of the running averaging or of the sum operator is required.
- FIG. 1 schematically illustrates the basic components of an X-ray computed tomography apparatus, constructed and operating in accordance with the invention.
- FIG. 2 is a flow chart for describing an exemplary embodiment of the inventive method for reducing line artifacts.
- FIG. 3 is a schematic representation of a median filtering in the inventive method.
- FIG. 4 is a schematic representation of a running averaging in the inventive method.
- FIG. 5 illustrates the extension of the method according to FIG. 2 given more than one defective detector channel.
- FIG. 6 is a flow chart for the inventive method for reducing line artifacts with an exact calculation of error images.
- FIG. 7 is a flow chart for the inventive method for reducing line artifacts with some of the data processing taking place in parallel.
- FIG. 8 shows a reference CT image of the skull base.
- FIG. 9 shows a CT image having high-contrast line artifacts.
- FIG. 10 shows the CT image of FIG. 9 after median filtering in accordance with the inventive method.
- FIG. 11 shows the difference value image formed from the respective images of FIG. 10 and FIG. 9 in accordance with the inventive method.
- FIG. 12 shows the resulting image G 1 1 formed by a running averaging of the difference value image of FIG. 11 in accordance with the inventive method.
- FIG. 13 shows the resulting image G 2 1 formed by a running averaging of the difference value image of FIG. 11 in accordance with the inventive method.
- FIG. 14 shows the correction image D 2 obtained in accordance with the inventive method.
- FIG. 1 shows a CT device, which is provided for scanning an examination subject 1 , having an X-ray source 2 , such as an X-ray tube, with a focus 3 from which a pyramid-shaped X-ray bundle 4 is emitted, which is gated by a radiation diaphragm (not shown) and which penetrates the examination subject 1 , for example a patient, and strikes a detector system 5 .
- the detector system 5 is an array of parallel rows 6 and parallel columns 7 of detector elements 8 .
- the X-ray source 2 and the detector system 5 form a measuring system, which can be rotated in a direction 6 around a system axis 9 and which can be displaced along the system axis 9 relative to the examination subject 1 , so that the examination subject 1 is irradiated from different projection angles and in different z-positions along the system axis 9 .
- a signal processing unit 10 forms measurement values from the output signals of the detector elements 8 of the detector system 5 . These measurement values are supplied to a computer 11 , which calculates an image of the examination subject 1 , which is reproduced at the monitor 12 .
- the flow chart as shown in FIG. 2 illustrates the inventive method for reducing line artifacts given a defective detector channel.
- the artifacts caused by the defective detector channel are situated on a circle K 1 . They are already corrected in a known manner in the CT image D 1 , so that the circle K 1 is no longer visible therein.
- line artifacts as tangents to the circle K 1 , appear in the CT image D 1 .
- a median filtering orthogonal to the straight line which extends through the respective picture element of the CT image D 1 and the center of the circle K 1 , produces the image M 1 from the CT image D 1 .
- a difference value image F 1 is produced by subtracting the image M 1 from the CT image D 1 .
- Two resulting images G 1 1 and G 2 1 are derived from the difference value image F 1 as a result of a further filter operation.
- the filtering ensues on the basis of each picture element of F 1 situated outside of the circle K 1 in the direction of the tangents to K 1 .
- the filtering is carried out along the tangent t 1 1 for the resulting image G 1 1 and is carried out along the tangent t 2 1 for the resulting image G 2 1 .
- the line artifacts that are present in the error image F 1 thus are emphasized in the resulting images G 1 1 , and G 2 1 .
- the two resulting images G 1 1 and G 2 1 are subsequently subtracted from the original CT image D 1 , and the CT image D 2 is obtained, which no longer contains any visible line artifacts.
- FIG. 3 schematically shows an exemplary median filtering. Shown are the circle K 1 with the center m and radius r, the two picture elements P 1 and P 2 , as well as the straight lines g 1 and g 2 , respectively through P 1 or P 2 and the center m. Furthermore, the positions of adjacent picture elements for P 1 and P 2 is indicated, which are situated on an orthogonal straight line with regard to g 1 and g 2 , respectively. Given the median filtering, a mean value is formed in the example from P 1 and the four illustrated adjacent picture elements, and from P 2 and the four corresponding adjacent picture elements, and is utilized as a new image value of P 1 and P 2 , respectively.
- FIG. 4 schematically shows pixel-dependent averages.
- the directions of the averages in the picture element P 1 are fixed by the tangents t 1 1 and t 2 1 to the circle K 1 having the radius r and the center m, this circle K 1 being defined by the defective detector channel.
- FIG. 5 shows the expansion of the method of FIG. 2 to more than one defective detector channel.
- the two result images G 1 i and G 2 i are produced by filtering the picture elements of F 1 in the direction of the tangents t 1 i , and t 2 1 to the circle K i extending through the respective picture element.
- the filtering along the tangent t 1 i leads to the resulting image G 1 i and the filtering along the tangent t 2 i leads to the resulting image G 2 i .
- the resulting images G 1 1 and G 2 1 are subtracted from the image D i .
- the image D i+1 arises.
- the index i is increased by 1 and the method is repeated until i>M.
- the error image F i is exactly calculated for each defective detector channel in the method shown in FIG. 6.
- the median filtering is carried out for each picture element of D i situated outside of the circle K i in the direction of the tangents to the circle K i extending through the respective picture element.
- the number of error images to be calculated in this method therefore corresponds to the number of defective detector channels.
- the two resulting images G 1 , and G 2 1 therefore are produced anew for each detector channel.
- the remaining method steps corresponds to the method steps described in FIGS. 2 and 5. This method, which is more complicated compared to that of FIG. 2 or FIG. 5, also is repeated until all defective detector channels are processed.
- FIG. 7 shows an exemplary method embodiment, wherein the resulting images G 1 i and G 2 i are not sequentially calculated as in the previously described embodiments for a number of defective detector channels, but instead the result images G 1 i and G 2 i , 1 ⁇ i ⁇ M, are calculated in parallel for the multiple defective detector channels.
- the image M 1 On the basis of the CT image D 1 , the image M 1 also initially arises as a result of a median filtering, and this image M 1 is subtracted from D 1 for calculating the image F 1 .
- the resulting images G 1 i and G 2 i , 1 ⁇ i ⁇ M, which are determined in parallel on the basis of F 1 also are subtracted from the CT image D 1 for producing the CT image D 2 .
- this method has a time advantage for the correction of the line artifacts in the CT images. Since the use of filter for calculating a resulting image has an effect on the adjacent resulting images, the defective detector elements should have a minimum spacing of ten channels from each other in order to be able to reconstruct a CT image that is free of line artifacts.
- FIG. 8 shows image data of the skull basis given a scan carried out in a 4 ⁇ 1 slice modus, axial scan operation.
- the measuring data of a detector row are evaluated, whereby a predefined individual channel does not supply measurement values.
- a linear interpolation eliminates the circle artifacts arising as a result of the “defective” channel.
- the result is the CT image, which serves as an initial CT image D 1 for the inventive method and which has a high-contrast line artifact at the location marked by the arrow, for example.
- the method as shown in FIG. 2 is used for eliminating the line artifacts.
- FIG. 10 shows the median-filtered input image M 1 .
- the difference value image F 1 essentially showing the line artifacts and image noise results as shown in FIG. 11.
- a running averaging in the direction of the tangents of the circle caused by the defective detector channel subsequently results for each pixel of the difference value image F 1 .
- two resulting images G 1 1 and G 2 1 are produced. These are shown in the FIGS. 12 and 13.
- the image noise is mainly suppressed in these images.
- CT devices of the third generation are shown, i.e., the X-ray source and the detector rotate together around the system axis during the image generation.
- the invention also can be used with CT devices of the fourth generation, wherein only the X-ray source rotates and cooperates with a stationary detector ring.
- the described exemplary embodiments relate to the medical application of inventive CT devices.
- the invention can also be applied outside of the medical field baggage inspection or material testing, for example.
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pulmonology (AREA)
- Radiology & Medical Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
Description
- 1. Field of the Invention
- The present invention is directed to a method and apparatus for reducing line artifacts in a CT image, wherein the CT image is produced by scanning with an X-ray source which is rotatable around an examination subject, with X-rays from the X-ray source, after being attenuated by an examination subject, being incident on a detector system.
- 2. Description of the Prior Art
- CT devices are known which have an X-ray source, e.g. an X-ray tube, which direct a collimated, pyramid-shaped X-ray bundle through the examination subject, e.g. a patient, onto a detector system that is composed of a number of detector channels. Each detector channel has at least one detector element and one associated electronic element for reading out and amplifying the signal that is generated in the detector element as a result of the incident radiation. A number of detector elements can be allocated to one electronic element. The X-ray source and, depending on the construction of the CT device, the detector system as well are attached to a gantry that rotates around the examination subject. A support device for the examination subject can be displaced along the system axis relative to the gantry. The position along the system axis at which the X-ray bundle penetrates the examination subject, and the angle, under which the X-ray bundle penetrates the examination subject, are continuously modified as a result of the displacement and the rotation of the gantry. Each detector element of the detector system struck by the radiation produces a signal representing a measure of the overall transparency of the examination subject for the radiation proceeding from the X-ray source to the detector system. The set of output signals of the detector elements of the detector system, which is acquired for a specific position of the X-ray source, is referred to as a projection. A scan is composed of a set of projections, which are acquired at different positions of the gantry and/or at different positions of the support device. The CT device picks up a number of projections during a scan in order to be able to construct a two-dimensional tomogram of a slice of the examination subject. A number of slices can be picked up at the same time by a detector system that is composed of an array having a number of rows and columns of detector elements. Such planar-like detector systems, however, frequently contain detector channels which do not supply proper data. It may be that detector system contains faulty detector channels already after the production process, for example as a result of defects in fabrication caused by the high integration density of the electronic elements. Defects of individual detector channels also may arise during the operation of the CT device. Such defects cause circular structures in the acquired CT images, these circular structures being referred to as circle artifacts. Techniques referred to as “ring-balancing” methods are known from the literature for the purpose of attempting to correct or prevent such artifacts in CT images. Such methods are disclosed in U.S. Pat. No. 4,670,840 and in German OS 198 35 451 (corresponding to U.S. Pat. No. 6,047,039), for example.
- A disadvantage of such known methods is that they insufficiently eliminate artifacts which arise in a CT device having a defective detector channel.
- An object of the present invention is to provide a method for reducing artifacts in a CT image, so that the obtainable image quality is improved in a CT device having at least one defective detector channel. It is also an object of the invention to provide a CT device for implementing the method.
- The above object is achieved in accordance with the principles of the present invention in a method for reducing line artifacts in a CT image D 1, as well as in an apparatus for implementing the method, wherein the image has been subjected to interpolated filtering for preventing or correcting faulty values of picture elements represented in a circle K1, and wherein the following steps are implemented. The picture elements of the CT image D1 are subjected to a median filtering, orthogonal to the straight line extending through the respective picture element and the center of the circle K1, for producing an image M1. A difference value image F1 is generated by subtracting the image M1 from the CT image D1. Two resulting images G1 1 and G2 1 are produced by filtering the picture elements of the difference value image F1 in the respective directions of tangents t1 1 and t2 1 to the circle K1 extending through the respective picture element. Filtering is conducted along t1 1 to produce the resulting image G1 1 and is conducted along t2 1 for producing the resulting image G2 1. A correction image D2 is then obtained by subtracting both of the resulting images G1 1 and G2 1 from the image D1.
- Defective detector channels of a detector system lead to faulty values for picture elements, which appear as a circle in a CT image acquired by the detector system. Such image errors therefore are referred to as circle artifacts. Defective detector channels are not only ones that fail to supply an output signal as a result of the defect, but also are channels with a measuring accuracy that exceeds a specific tolerance value. Various methods are known for correcting or preventing circle artifacts, as noted above. These methods are carried out on the measurement data or on the image data and are primarily based on interpolated filtering, and achieve a noticeable weakening of the circle artifacts in the acquired CT images. Such known methods have the disadvantage, however, that line-like image errors frequently arise in the resulting CT images after such a method has been implemented. These image errors are referred to as line artifacts. They increasingly occur in association with large signal unsteadiness caused by high-contrast areas of an examination subject. The inventive method is particularly advantageous for eliminating such line artifacts, which occur after circle artifacts have been eliminated and which appear as tangents to the circles in the CT image. The size and position data of the circles of the circle artifacts caused by the defective detector channels are assumed to be known. It is sufficient to know the position of the circle center, which is the same for all circles, and the radii of the circles.
- The elimination of the line artifacts is initially described for the case of a single defective detector channel. The faulty values of picture elements caused by the defective detector channel are situated on a circle K 1. These faulty values are corrected by a known ring-balancing method on the measuring data or on the image data. Line artifacts, which appear as tangents to the circle K1, arise in the resulting CT image D1. For producing an image M1, a median filtering is carried out for each picture element of the CT image D1 situated outside of the circle K1, orthogonally to the straight line extending through the respective picture element and the center of the circle K1. The width of the median filtering can be modifiable. This width should be selected wider than the expected line width of the line artifacts. A median filter of the
width 5 has proven to be beneficial. The sampling distance A1 must be selected dependent on the convolution kernel used for the reconstruction. Ideally, the image M1 no longer contains line artifacts. - A difference value image F 1=D1−M1 is generated by subtracting the image M1 from the CT image D1. For producing two resulting images G1 1 and G2 1, a filtering is carried out in each picture element of the difference value image F1 in the respective directions of the tangents t1 1 and t2 1 to the circle K1 extending through the picture element in question, i.e., the filtering is carried out along t1 1 for the resulting image G1 1 and the filtering is carried out along t2 1 for the resulting image G2 1. This step is to eliminate, for the most part, the image pixel noise difference value in the image F1 for the most part. Furthermore, the line artifacts that are present in the error image F1 are emphasized in the result images. It is necessary to calculate two result images G1 1 and G2 1, since there are two possible tangent directions to the circle K1 per picture element. The thus-determined resulting images G1 1 and G2 1 are subtracted from D1 and the correction image D2 is obtained, which ideally no longer contains line artifacts generated by the defective channel.
- The correction image D 2 serves as basis image for correcting the line artifacts caused by other detector channels if the detector system contains further defective detector channels. In general, the correction image D1 serves as a basis image for correcting line artifacts caused by the i-th defective detector channel. Since the image M1 has been calculated in relation to the center of the circle that is valid for all circles ki, the difference value image F1 can still be used in the following. For producing the result images G1 i and G2 i, a filtering is carried out in each picture element of the difference value image F1 situated outside of the circle Ki in the direction of the tangents D1 i and D2 i to the circle Ki, i.e., the filtering is carried out along D1 i for the resulting image G1 i and the filtering is carried out along D2 i for the resulting image G2 i. The correction image Di−1=Di−G1 i−G2 i is determined by subtracting the resulting images G1 i and G2 i from Di. This procedure is repeated until the line artifacts of all defective detector channels are eliminated.
- The above-described method is relatively time-consuming as a result of the serial processing of the line artifacts caused by the defective detector channels. Parallel processing of the image errors caused by the different defective detector channels is preferable for processing the data faster. Accordingly, a further embodiment of the inventive method proceeds as described above until the calculation of the difference value image F 1. The resulting images G1 i and G2 i are calculated in this manner for all defective detector channels. In contrast to the initially described embodiment, different correction images are not consecutively determined, which respectively serve as an image for calculating the next correction images, but instead all resulting images G1 i and G2 1 are subtracted from D1. The correction image D2 results. The defective detector channels should have a minimum distance (spacing) of ten detector channels as a condition for this parallel processing.
- The aforementioned median filtering for producing an image M 1, in each picture element of the initial CT image, takes place orthogonally to the straight line extending through the picture element in question and the center of the circle Ki. This calculation represents an approximation, and it is an advantage of the inventive method that it must be carried out only once for each picture element even given a number of defective detector channels. Better results are obtained, however, by the median filtering orthogonal to the tangents to the circles Ki. This makes it necessary to carry out the pixel-oriented median filtering anew for each defective detector channel, with an associated high computing outlay. Given a single defective detector channel and the faulty values of picture elements on a circle K1 caused as a result thereof, the correction includes the following steps:
- carrying out an interpolated filtering of the artifact-containing CT image, thereby leading to the CT image D 1, in order to avoid or to correct faulty values of picture elements on the circle K1,
- carrying out a median filtering in each picture element of the CT image D 1 situated outside of the circle K1, orthogonal to the tangents to the circle K1 extending through the respective picture element, for producing an image M1,
- generating a difference value image F 1D1−M1 by subtracting the image M1 from the CT image D1,
- carrying out a filtering in each picture element of the difference value image F 1 situated outside of the circle K1 in the direction of the tangents t1 1 and t2 1 to the circle K1 extending through the respective picture element in order to produce two resulting images G1 1 and G2 1, with the filtering being carried out along t1 1 for the resulting image G1 1 and the filtering is carried out along t2 1 for the resulting image G2 1.
- determining the correction image D 2=D1−G1 1−G2 1 by subtracting the resulting images G1 1 and G2 1 are from D1.
- If a number of defective detector channels are present, the following applies for processing the i-th defective detector channel:
- a median filtering is carried out in each picture element of the CT image D i, orthogonal to the tangents to the circle Ki extending through the respective picture element, for producing an image Mi,
- a difference value image F i=Di−Mi is produced by subtracting the image Mi from the correction image Di,
- a filtering is carried out in each picture element of the difference value image F 1 situated outside of the circle K1 in the direction of the tangents t1 i and t2 i to the circle Ki extending through the respective picture element in order to produce two resulting images G1 i and G2 i, with the filtering being carried out along t1 i for the resulting image G1 i and the filtering being carried out along t2 i for the resulting image G2 i.
- the correction image D i−1=Di−Di−G1 i−G2 i is determined by subtracting the resulting images G1 i and G2 i from Di.
- These steps are repeated until all image errors caused by the different defective detector channels are processed.
- In a further version of the invention, a running averaging serves as the filtering for producing the resulting images G 1 i and G2 i. In another version of the invention, a sum operator is implemented as the filtering for producing the resulting images G1 i and G2 i. A combination of these two versions also can be employed, wherein a running averaging and a sum operator are simultaneously employed for filtering.
- At the beginning of each of the aforementioned embodiments, the values of the CT image D 1 are preferably limited to a range 1000 HU (H2O)±Δ. This limitation is expedient for examining soft-tissue parts, in which disturbances caused by defective detector channels are particularly apparent in the image.
- The intensity of the noise signal in the resulting images G 1 i and G2 i is dependent on the filter width of the running averaging or of the sum operator. This noise signal influences the correction images D1+1 and therefore can lead to undesired noise structures in the resulting image in the case of a number of defective detector channels. A high filter width is required as a result. The intensity of a line artifact varies, however, so that a limitation of the filter width of the running averaging or of the sum operator is required.
- The summation of the HU values in G 1 i in the tangent direction D1 i or in G2 i in the tangent direction D2 i represents an additional measure for suppressing noise. This corresponds to the calculation of the radon values in these tangent directions. This operation of the image processing is referred to as HUG transformation. On the basis of a threshold criterion, a noise signal can be principally differentiated from artifact structures and can be eliminated. This threshold must be suitably selected, however, so that low-contrast lines in G1 i and G2 i, that are actually a part of the diagnostically relevant image context, are not erroneously eliminated.
- FIG. 1 schematically illustrates the basic components of an X-ray computed tomography apparatus, constructed and operating in accordance with the invention.
- FIG. 2 is a flow chart for describing an exemplary embodiment of the inventive method for reducing line artifacts.
- FIG. 3 is a schematic representation of a median filtering in the inventive method.
- FIG. 4 is a schematic representation of a running averaging in the inventive method.
- FIG. 5 illustrates the extension of the method according to FIG. 2 given more than one defective detector channel.
- FIG. 6 is a flow chart for the inventive method for reducing line artifacts with an exact calculation of error images.
- FIG. 7 is a flow chart for the inventive method for reducing line artifacts with some of the data processing taking place in parallel.
- FIG. 8 shows a reference CT image of the skull base.
- FIG. 9 shows a CT image having high-contrast line artifacts.
- FIG. 10 shows the CT image of FIG. 9 after median filtering in accordance with the inventive method.
- FIG. 11 shows the difference value image formed from the respective images of FIG. 10 and FIG. 9 in accordance with the inventive method.
- FIG. 12 shows the resulting image G 1 1 formed by a running averaging of the difference value image of FIG. 11 in accordance with the inventive method.
- FIG. 13 shows the resulting image G 2 1 formed by a running averaging of the difference value image of FIG. 11 in accordance with the inventive method.
- FIG. 14 shows the correction image D 2 obtained in accordance with the inventive method.
- FIG. 1 shows a CT device, which is provided for scanning an
examination subject 1, having anX-ray source 2, such as an X-ray tube, with afocus 3 from which a pyramid-shapedX-ray bundle 4 is emitted, which is gated by a radiation diaphragm (not shown) and which penetrates theexamination subject 1, for example a patient, and strikes adetector system 5. Thedetector system 5 is an array ofparallel rows 6 andparallel columns 7 ofdetector elements 8. TheX-ray source 2 and thedetector system 5 form a measuring system, which can be rotated in adirection 6 around asystem axis 9 and which can be displaced along thesystem axis 9 relative to theexamination subject 1, so that theexamination subject 1 is irradiated from different projection angles and in different z-positions along thesystem axis 9. Asignal processing unit 10 forms measurement values from the output signals of thedetector elements 8 of thedetector system 5. These measurement values are supplied to acomputer 11, which calculates an image of theexamination subject 1, which is reproduced at themonitor 12. - The flow chart as shown in FIG. 2 illustrates the inventive method for reducing line artifacts given a defective detector channel. The artifacts caused by the defective detector channel are situated on a circle K 1. They are already corrected in a known manner in the CT image D1, so that the circle K1 is no longer visible therein. As an undesired side effect of this correction, line artifacts, as tangents to the circle K1, appear in the CT image D1. A median filtering orthogonal to the straight line, which extends through the respective picture element of the CT image D1 and the center of the circle K1, produces the image M1 from the CT image D1. A difference value image F1 is produced by subtracting the image M1 from the CT image D1. Two resulting images G1 1 and G2 1 are derived from the difference value image F1 as a result of a further filter operation. The filtering ensues on the basis of each picture element of F1 situated outside of the circle K1 in the direction of the tangents to K1. The filtering is carried out along the tangent t1 1 for the resulting image G1 1 and is carried out along the tangent t2 1 for the resulting image G2 1. The line artifacts that are present in the error image F1 thus are emphasized in the resulting images G1 1, and G2 1. The two resulting images G1 1 and G2 1 are subsequently subtracted from the original CT image D1, and the CT image D2 is obtained, which no longer contains any visible line artifacts.
- FIG. 3 schematically shows an exemplary median filtering. Shown are the circle K 1 with the center m and radius r, the two picture elements P1 and P2, as well as the straight lines g1 and g2, respectively through P1 or P2 and the center m. Furthermore, the positions of adjacent picture elements for P1 and P2 is indicated, which are situated on an orthogonal straight line with regard to g1 and g2, respectively. Given the median filtering, a mean value is formed in the example from P1 and the four illustrated adjacent picture elements, and from P2 and the four corresponding adjacent picture elements, and is utilized as a new image value of P1 and P2, respectively.
- FIG. 4 schematically shows pixel-dependent averages. The directions of the averages in the picture element P 1 are fixed by the tangents t1 1 and t2 1 to the circle K1 having the radius r and the center m, this circle K1 being defined by the defective detector channel.
- FIG. 5 shows the expansion of the method of FIG. 2 to more than one defective detector channel. On the basis of the difference value image F 1, the two result images G1 i and G2 i are produced by filtering the picture elements of F1 in the direction of the tangents t1 i, and t2 1 to the circle Ki extending through the respective picture element. The filtering along the tangent t1 i leads to the resulting image G1 i and the filtering along the tangent t2 i leads to the resulting image G2 i. The resulting images G1 1 and G2 1 are subtracted from the image Di. The image Di+1 arises. The index i is increased by 1 and the method is repeated until i>M. The last-determined image Di with i=M+1 is displayed as the CT image at the monitor.
- In contrast to FIG. 5, the error image F i is exactly calculated for each defective detector channel in the method shown in FIG. 6. For this purpose, the median filtering is carried out for each picture element of Di situated outside of the circle Ki in the direction of the tangents to the circle Ki extending through the respective picture element. The number of error images to be calculated in this method therefore corresponds to the number of defective detector channels. In addition to the error image, the two resulting images G1, and G2 1 therefore are produced anew for each detector channel. The remaining method steps corresponds to the method steps described in FIGS. 2 and 5. This method, which is more complicated compared to that of FIG. 2 or FIG. 5, also is repeated until all defective detector channels are processed.
- FIG. 7 shows an exemplary method embodiment, wherein the resulting images G 1 i and G2 i are not sequentially calculated as in the previously described embodiments for a number of defective detector channels, but instead the result images G1 i and G2 i, 1≦i≦M, are calculated in parallel for the multiple defective detector channels. On the basis of the CT image D1, the image M1 also initially arises as a result of a median filtering, and this image M1 is subtracted from D1 for calculating the image F1. The resulting images G1 i and G2 i, 1≦i≦M, which are determined in parallel on the basis of F1, also are subtracted from the CT image D1 for producing the CT image D2. Assuming the availability of a very powerful computer, this method has a time advantage for the correction of the line artifacts in the CT images. Since the use of filter for calculating a resulting image has an effect on the adjacent resulting images, the defective detector elements should have a minimum spacing of ten channels from each other in order to be able to reconstruct a CT image that is free of line artifacts.
- The efficacy of the method is demonstrated below on the basis of actual image data, acquired with a
SIEMENS SOMATOM PLUS 4 VOLUME ZOOM CT device. - As a reference CT image, FIG. 8 shows image data of the skull basis given a scan carried out in a 4×1 slice modus, axial scan operation. The measuring data of a detector row are evaluated, whereby a predefined individual channel does not supply measurement values. A linear interpolation eliminates the circle artifacts arising as a result of the “defective” channel. As shown in FIG. 9, the result is the CT image, which serves as an initial CT image D 1 for the inventive method and which has a high-contrast line artifact at the location marked by the arrow, for example. The method as shown in FIG. 2 is used for eliminating the line artifacts. FIG. 10 shows the median-filtered input image M1. Given the subtraction of M1 from the original CT image D1, the difference value image F1 essentially showing the line artifacts and image noise results as shown in FIG. 11. A running averaging in the direction of the tangents of the circle caused by the defective detector channel subsequently results for each pixel of the difference value image F1. Since—proceeding from each pixel of the difference value image F1—there are two possible tangents to the circle, two resulting images G1 1 and G2 1 are produced. These are shown in the FIGS. 12 and 13. Ideally, the image noise is mainly suppressed in these images. When the two resulting images G1 1, and G2 1 are subtracted from the original CT image D1 as shown in FIG. 9, the correction image D2 as shown in FIG. 14 results. As can be seen from D2, the line artifacts have been effectively suppressed without increasing the noise contribution.
- For the practical realization of a method according to the invention, only the corrected images are displayed at the monitor. The method process therefore proceeds essentially on the basis of computer-internal “image data”, which are not reconstructed as “images” at the monitor.
- In the described exemplary embodiments, CT devices of the third generation are shown, i.e., the X-ray source and the detector rotate together around the system axis during the image generation. The invention also can be used with CT devices of the fourth generation, wherein only the X-ray source rotates and cooperates with a stationary detector ring.
- The described exemplary embodiments relate to the medical application of inventive CT devices. The invention can also be applied outside of the medical field baggage inspection or material testing, for example.
- Although modifications and changes may be suggested by those skilled in the art, it is the intention of the inventors to embody within the patent warranted hereon all changes and modifications as reasonably and properly come within the scope of their contribution to the art.
Claims (40)
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| DE10009746 | 2000-03-01 | ||
| DE10009746A DE10009746B4 (en) | 2000-03-01 | 2000-03-01 | Method for reducing line artifacts in a CT image |
| DE10009746.4 | 2000-03-01 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| US20020012417A1 true US20020012417A1 (en) | 2002-01-31 |
| US6411671B2 US6411671B2 (en) | 2002-06-25 |
Family
ID=7632956
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| US09/797,330 Expired - Lifetime US6411671B2 (en) | 2000-03-01 | 2001-03-01 | Method for reducing line artifacts in a CT image and device for implementing the method |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US6411671B2 (en) |
| JP (1) | JP2001276055A (en) |
| DE (1) | DE10009746B4 (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030161444A1 (en) * | 2001-08-16 | 2003-08-28 | Alexander Katsevich | Method of reconstructing images for spiral and non-spiral computer tomography |
| US6804321B2 (en) | 2001-08-16 | 2004-10-12 | University Of Central Florida | Filtered back projection (FBP) algorithm for computer tomography |
| US20050190984A1 (en) * | 2004-02-24 | 2005-09-01 | Daniel Fischer | Method for filtering tomographic 3D images after completed reconstruction of volume data |
| US20060029180A1 (en) * | 2001-08-16 | 2006-02-09 | Alexander Katsevich | Exact filtered back projection (FBP) algorithm for spiral computer tomography with variable pitch |
| US20090114826A1 (en) * | 2006-09-29 | 2009-05-07 | Isao Takahashi | Nuclear medical diagnosis apparatus |
| CN102451015A (en) * | 2010-10-20 | 2012-05-16 | 上海西门子医疗器械有限公司 | CT device and method for determining unstable channel in detector |
| CN113610716A (en) * | 2021-07-02 | 2021-11-05 | 中铁二十局集团有限公司 | Image artifact eliminating method, device and equipment |
| CN116165713A (en) * | 2021-11-24 | 2023-05-26 | 中国石油天然气集团有限公司 | Method and device for denoising seismic data of DAS in well |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3847101B2 (en) * | 2001-05-22 | 2006-11-15 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X-ray CT apparatus and method |
| CN1786819B (en) | 2004-12-09 | 2011-08-10 | Ge医疗系统环球技术有限公司 | X ray diaphragm, X ray radiator and X ray imaging apparatus |
| JP6497912B2 (en) * | 2014-12-01 | 2019-04-10 | キヤノン株式会社 | Image processing apparatus, radiation imaging system, control method, and program |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4670840A (en) * | 1983-03-09 | 1987-06-02 | Elscint, Inc. | Ring artifact correction for computerized tomography |
| JPH0323847A (en) * | 1989-06-21 | 1991-01-31 | Toshiba Corp | X-ray ct scanner apparatus |
| IL119714A0 (en) * | 1996-11-28 | 1997-02-18 | Elscint Ltd | CT system with oblique image planes |
| US5960056A (en) * | 1997-07-01 | 1999-09-28 | Analogic Corporation | Method and apparatus for reconstructing volumetric images in a helical scanning computed tomography system with multiple rows of detectors |
| DE19835451B4 (en) * | 1997-08-20 | 2005-03-24 | Siemens Ag | Method for a computer tomograph for post-processing a sectional image and computer tomograph operating according to this method |
| US6115445A (en) * | 1999-01-12 | 2000-09-05 | Analogic Corporation | Progressive correction of ring artifacts in a computed tomography system |
-
2000
- 2000-03-01 DE DE10009746A patent/DE10009746B4/en not_active Expired - Fee Related
-
2001
- 2001-03-01 JP JP2001056557A patent/JP2001276055A/en not_active Withdrawn
- 2001-03-01 US US09/797,330 patent/US6411671B2/en not_active Expired - Lifetime
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7280632B2 (en) | 2001-08-16 | 2007-10-09 | University Of Central Florida Research Foundation, Inc. | Exact filtered back projection (FBP) algorithm for spiral computer tomography with variable pitch |
| US6771733B2 (en) | 2001-08-16 | 2004-08-03 | University Of Central Florida | Method of reconstructing images for spiral and non-spiral computer tomography |
| US6804321B2 (en) | 2001-08-16 | 2004-10-12 | University Of Central Florida | Filtered back projection (FBP) algorithm for computer tomography |
| US6898264B2 (en) | 2001-08-16 | 2005-05-24 | Research Foundation Of The University Of Central Florida | Method of reconstructing images for spiral and non-spiral computer tomography |
| US20030161444A1 (en) * | 2001-08-16 | 2003-08-28 | Alexander Katsevich | Method of reconstructing images for spiral and non-spiral computer tomography |
| US20060029180A1 (en) * | 2001-08-16 | 2006-02-09 | Alexander Katsevich | Exact filtered back projection (FBP) algorithm for spiral computer tomography with variable pitch |
| US7650023B2 (en) * | 2004-02-24 | 2010-01-19 | Siemens Aktiengeśellschaft | Method for filtering tomographic 3D images after completed reconstruction of volume data |
| US20050190984A1 (en) * | 2004-02-24 | 2005-09-01 | Daniel Fischer | Method for filtering tomographic 3D images after completed reconstruction of volume data |
| US20090114826A1 (en) * | 2006-09-29 | 2009-05-07 | Isao Takahashi | Nuclear medical diagnosis apparatus |
| US7795590B2 (en) * | 2006-09-29 | 2010-09-14 | Hitachi, Ltd. | Nuclear medical diagnosis apparatus |
| US20100282974A1 (en) * | 2006-09-29 | 2010-11-11 | Hitachi, Ltd. | Nuclear medical diagnosis apparatus |
| US7964849B2 (en) | 2006-09-29 | 2011-06-21 | Hitachi, Ltd. | Nuclear medical diagnosis apparatus |
| US20110215254A1 (en) * | 2006-09-29 | 2011-09-08 | Hitachi, Ltd. | Nuclear medical diagnosis apparatus |
| US8148695B2 (en) | 2006-09-29 | 2012-04-03 | Hitachi, Ltd. | Nuclear medical diagnosis apparatus |
| CN102451015A (en) * | 2010-10-20 | 2012-05-16 | 上海西门子医疗器械有限公司 | CT device and method for determining unstable channel in detector |
| CN113610716A (en) * | 2021-07-02 | 2021-11-05 | 中铁二十局集团有限公司 | Image artifact eliminating method, device and equipment |
| CN116165713A (en) * | 2021-11-24 | 2023-05-26 | 中国石油天然气集团有限公司 | Method and device for denoising seismic data of DAS in well |
Also Published As
| Publication number | Publication date |
|---|---|
| US6411671B2 (en) | 2002-06-25 |
| DE10009746B4 (en) | 2008-11-20 |
| DE10009746A1 (en) | 2001-09-20 |
| JP2001276055A (en) | 2001-10-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7444010B2 (en) | Method and apparatus for the reduction of artifacts in computed tomography images | |
| US7023951B2 (en) | Method and apparatus for reduction of artifacts in computed tomography images | |
| US4709333A (en) | Method and apparatus for imaging in the presence of multiple high density objects | |
| US7801264B2 (en) | Method for calibrating a dual -spectral computed tomography (CT) system | |
| US5416815A (en) | Adaptive filter for reducing streaking artifacts in x-ray tomographic images | |
| US6115487A (en) | Correction algorithm for bone-induced spectral artifacts in computed tomograph imaging | |
| US8737711B2 (en) | X-ray CT image forming method and X-ray CT apparatus using the same | |
| JP5142664B2 (en) | X-ray computed tomography system | |
| US6493416B1 (en) | Method and apparatus for noise reduction in computed tomographic systems | |
| US5727041A (en) | Methods and apparatus for reducing partial volume image artifacts | |
| US7920672B2 (en) | X-ray detector gain calibration depending on the fraction of scattered radiation | |
| US7929659B2 (en) | System and method for generating computed tomography images | |
| US6411671B2 (en) | Method for reducing line artifacts in a CT image and device for implementing the method | |
| JPWO2005011502A1 (en) | Radiation tomography equipment | |
| US7283605B2 (en) | Methods and apparatus for scatter correction | |
| US7747057B2 (en) | Methods and apparatus for BIS correction | |
| US5761257A (en) | Normalizing projection data in a computed tomography system | |
| JPH08252249A (en) | Method to remove artifact due to deterioration of detector from image data, and system to form laminagram picture of article | |
| US6980681B1 (en) | Methods and apparatus for helical reconstruction for multislice CT scan | |
| JP3484288B2 (en) | X-ray tomography equipment | |
| JP2004531306A (en) | A method for reducing artifacts in target images | |
| US6307908B1 (en) | System and method for data interpolation in a multislice x-ray computed tomography system | |
| JP2006102299A (en) | X-ray dose correcting method and x-ray ct apparatus | |
| CN111053568B (en) | Method and device for correcting ring artifact in CT image and computer storage medium | |
| US6307912B1 (en) | Methods and apparatus for optimizing CT image quality with optimized data acquisition |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| AS | Assignment |
Owner name: SIEMENS AKTIENGESELLSCHAFT, GERMANY Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:BRUDER, DR. HERBERT;FLOHR, DR. THOMAS;STIERSTORFER, DR. KARL;AND OTHERS;REEL/FRAME:012128/0683 Effective date: 20010820 |
|
| 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 |