US7668360B2 - Image processing method and X-ray CT system - Google Patents
Image processing method and X-ray CT system Download PDFInfo
- Publication number
- US7668360B2 US7668360B2 US11/531,899 US53189906A US7668360B2 US 7668360 B2 US7668360 B2 US 7668360B2 US 53189906 A US53189906 A US 53189906A US 7668360 B2 US7668360 B2 US 7668360B2
- Authority
- US
- United States
- Prior art keywords
- vectors
- wavelet
- vector
- threshold
- classified
- 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.)
- Active, expires
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 12
- 239000013598 vector Substances 0.000 claims abstract description 175
- 239000011159 matrix material Substances 0.000 claims description 14
- 238000009499 grossing Methods 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 4
- 230000001131 transforming effect Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 4
- 238000002591 computed tomography Methods 0.000 description 17
- 238000001514 detection method Methods 0.000 description 13
- 230000009467 reduction Effects 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 3
- 238000000034 method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000007688 edging Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/10—Image enhancement or restoration by non-spatial domain filtering
-
- G06T5/70—
-
- G06T5/73—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20064—Wavelet transform [DWT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Definitions
- the present invention relates to an image processing method and an X-ray computed tomography (CT) system. More particularly, the present invention is concerned with a method of wavelet-transforming an image, reducing noises in the wavelet domain, and then inverse-wavelet-transforming the wavelet-domain image so as to restore the image, and an X-ray CT system that performs the image processing.
- CT computed tomography
- X-irradiation/detection equipment that rotates in a gantry acquires transmitted X-ray signals, which represent a plurality of views, from a subject, and tomographic images are reconstructed based on the transmitted X-ray signals.
- the rotation of the X-irradiation/detection equipment may be referred to as a scan.
- the acquired transmitted X-ray signal may be referred to as scan data.
- a reconstructed image has noises therein reduced.
- One of noise reduction methods is such that an image is wavelet-transformed, the resultant wavelet-domain signal has noises thereof reduced, and the resultant signal is then inverse-wavelet-transformed in order to restore the image.
- the noise reduction in the wavelet domain is achieved by replacing signal components, of which levels fall below a predetermined threshold, with 0s (refer to, for example, Japanese Unexamined Patent Publication No. 2002-133399 (p.3, FIG. 1 to FIG. 7).
- the foregoing noise reduction method of simply replacing signal components, of which levels fall below a predetermined threshold, with 0s causes a restored image to suffer a shift in a CT number, deformation of edges, an abrupt change, or any other side effect.
- An object of the present invention is to provide an image processing method for reducing noises without causing a side effect and to provide an X-ray CT system that performs the image processing.
- an image processing method in which: an input image is wavelet-transformed using a gradient form; resultant wavelet vectors are classified into categories of strongly edged objects, weakly edged objects, radial lines, and noises; wavelet vectors classified into the category of strongly edged objects are left intact; wavelet vectors classified into the category of weakly edged objects have the directions thereof filtered; and wavelet vectors classified into the categories of radial lines and noises are smoothed; the resultant wavelet vectors are inverse-wavelet-transformed in order to produce an output image.
- an X-ray CT system comprising a data acquisition means for acquiring data by scanning a subject with X-rays, an image reconstruction means for reconstructing an image according to the acquired data, and an image processing means for processing the reconstructed image.
- the image processing means comprises: a transforming means for wavelet-transforming an input image using a gradient form; a classifying means for classifying resultant wavelet vectors into categories of strongly edged objects, weakly edged objects, radial lines, and noise; a processing means for leaving wavelet vectors, which are classified into the category of strongly edged objects, intact, filtering the directions of wavelet vectors classified into the category of weakly edged objects, and smoothing wavelet vectors that are classified into the categories of radial lines and noise; and an image production means for inverse-wavelet-transforming the resultant wavelet vectors so as to produce an output image.
- the classification employs a first threshold concerning the magnitude of a vector, a second threshold smaller than the first threshold, and a third threshold concerning the regulation in the direction of a vector.
- Wavelet vectors whose magnitudes are larger than the first threshold are classified into the category of strongly edged objects.
- Wavelet vectors whose magnitudes are equal to or smaller than the first threshold and larger than the second threshold and whose regulations in directions thereof are larger than the third threshold are classified into the category of weakly edged objects.
- Wavelet vectors whose magnitudes are equal to or smaller than the first threshold and whose regulations in directions thereof are equal to or smaller than the third threshold are classified into the category of radial lines.
- Wavelet vectors whose magnitudes are equal to or smaller than the second threshold are classified into the category of noises.
- the regulation in the direction of a vector is expressed as a ratio of a mean vector, which is a mean of a center vector located in the center of a window defined with a matrix of 3 by 3 and two ambient vectors whose directions from the center vector is the smallest, to a means vector that is a mean of all vectors within the window.
- a mean vector which is a mean of a center vector located in the center of a window defined with a matrix of 3 by 3 and two ambient vectors whose directions from the center vector is the smallest
- the filtering of directions is achieved by sampling a combination of a center vector, which is located in the center of a window defined with a matrix of 3 by 3, and two ambient vectors the sum of whose differences from the center vector is the smallest, and replacing the center vector with a mean vector that is a mean of the two ambient vectors and the center vector.
- the smoothing of vectors is achieved by replacing a center vector, which is located in the center of a window defined with a matrix of 3 by 3, with a mean vector that is a mean of all vectors in the window.
- an input image is wavelet-transformed using a gradient form.
- Resultant wavelet vectors are classified into categories of strongly edged objects, weakly edged objects, radial lines, and noises.
- the wavelet vectors classified into the category of strongly edged objects are left intact.
- the wavelet vectors classified into the category of weakly edged objects have the directions thereof filtered.
- the wavelet vectors classified into the categories of radial lines and noises are smoothed.
- the resultant wavelet vectors are inverse-wavelet-transformed in order to produce an output image. Consequently, there are provided an image processing method for reducing noises without causing a side effect and an X-ray CT system that performs the image processing.
- FIG. 1 shows the configuration of an X-ray CT system that is an example of the best mode for implementing the present invention.
- FIG. 2 shows the configuration of an X-ray CT system that is an example of the best mode for implementing the present invention.
- FIG. 3 shows the configuration of X-irradiation/detection equipment.
- FIG. 4 shows the structure of an X-ray incidence surface of an X-ray detector.
- FIG. 5 describes a flow of image processing that is an example of the best mode for implementing the present invention.
- FIG. 6 shows the classification of wavelet vectors.
- FIG. 7 shows directions to be filtered.
- FIGS. 8 a and 8 b show the concept of smoothing of vectors.
- FIGS. 9 a , 9 b , 9 c , and 9 d show the concept of image processing that is an example of the best mode for implementing the present invention.
- FIG. 10 shows comparison between an original image and a processed image using partly half-tone photographs.
- FIG. 11 shows comparison between an original image and a processed image using partly half-tone photographs.
- FIG. 12 shows comparison between an original image and a processed image using partly half-tone photographs.
- FIGS. 13 a - 13 c show comparison between an original image and a processed image using partly half-tone photographs.
- FIG. 1 illustratively shows the configuration of an X-ray CT system.
- the X-ray CT system is an example of the best mode for implementing the present invention.
- the configuration of the X-ray CT system presents an example of the best mode for implementing the present invention in an X-ray CT system.
- actions to be performed in the X-ray CT system present an example of the best mode for implementing the present invention in an image processing method.
- An X-ray CT system comprises a gantry 100 , a table 200 , and an operator console 300 .
- the gantry 100 uses X-irradiation/detection equipment 110 to scan a subject 10 who is carried while lying down on the table 200 , acquires transmitted X-ray signals (scan data items) representing a plurality of views, and transfers the views to the operator console 300 .
- the gantry 100 is an example of a data acquisition means included in the present invention.
- the operator console 300 reconstructs an image according to the scan data received from the gantry 100 , reduces noises contained in the reconstructed image, and displays a processed image on a display 302 .
- the operator console 300 is an example of an image reconstruction means included in the present invention.
- the operator console 300 is an example of an image processing means included in the present invention.
- the operator console 300 controls the movements of the gantry 100 and the table 200 alike. Under the control of the operator console 300 , the gantry 100 scans a subject under predetermined scan conditions, and the table 200 positions the subject 10 so that a predetermined region will be scanned. The positioning is such that an incorporated alignment mechanism adjusts the height of a tabletop 202 and the distance of a horizontal movement made by a cradle 204 placed on the tabletop.
- axial scan is achieved.
- helical scan is achieved. Every time the cradle 204 is brought to a halt while intermittently moved, if a subject is scanned, cluster scan is achieved.
- a support column 206 is swung with the root of the support column fixed to a base 208 as a center. Along with the swing of the support column 206 , the tabletop 202 is displaced in vertical and horizontal directions. The cradle 204 is moved horizontally on the tabletop 202 , whereby the horizontal displacement of the tabletop 202 is canceled out.
- a subject is scanned with the gantry 100 tilted. The gantry 100 is tilted by an incorporated tilt mechanism.
- the table 200 may be, as shown in FIG. 2 , of a type having the tabletop 202 thereof lifted or lowered vertically with respect to the base 208 .
- the tabletop 202 is lifted or lowered by an incorporated lifting/lowering mechanism.
- the tabletop 202 of the table 200 will not be horizontally moved along with the lifting or lowering.
- FIG. 3 illustratively shows the configuration of the X-irradiation/detection equipment 110 .
- the X-irradiation/detection equipment 110 uses an X-ray detector 150 to detect X-rays 134 radiated from a focal spot 132 in an X-ray tube 130 .
- the X-rays 134 are recomposed by a collimator that is not shown into cone-beam or fan-beam X-rays that are laterally symmetrical.
- the X-ray detector 150 has an X-ray incidence surface 152 that spreads two-dimensionally in line with the spread of X-rays.
- the X-ray incidence surface 152 is curved like part of a cylinder whose center axis passes through the focal spot 132 .
- the X-irradiation/detection equipment 110 rotates about a center axis passing through an isocenter that is a center of radiography.
- the center axis is parallel to the center axis of the cylinder part of which is formed by the X-ray detector 150 .
- the direction of the center axis of the rotation shall be regarded as a z direction
- the direction of a line linking the isocenter O and the focal spot 132 shall be regarded as a y direction
- the direction perpendicular to the z and y directions shall be regarded as an x direction.
- the axes extending in the x, y, and z directions respectively serve as three axes of a rotating coordinate system having the z axis as a center axis.
- FIG. 4 is a plan view illustratively showing the X-ray incidence surface 152 of the X-ray detector 150 .
- the X-ray incidence surface 152 has detection cells 154 two-dimensionally arranged in the x and z directions. In other words, the X-ray incidence surface 152 is realized with a two-dimensional array of detection cells 154 . When fan-beam X-rays are employed, the X-ray incidence surface 152 may be realized with a one-dimensional array of detection cells 154 .
- the individual detection cells 154 provide the X-ray detector 150 with detection channels.
- the X-ray detector 150 is therefore a multi-channel X-ray detector.
- the detection cells 154 are each composed of, for example, a scintillator and a photodiode.
- FIG. 5 describes a flow of noise reduction. The described noise reduction is achieved by a computer incorporated in the operator console 300 .
- an input image is wavelet-transformed in stage 501 .
- the computer that wavelet-transforms the input image in stage 501 is an example of a transforming means included in the present invention.
- Wavelet transform is achieved according to the formula below.
- the above formula expresses wavelet transform that employs a two-dimensional gradient form.
- the result of the transform is a two-dimensional distribution of wavelet vectors expressing edges of objects contained in an input image.
- the two-dimensional distribution of wavelet vectors is an expression of an input image in the wavelet domain.
- the magnitude of a wavelet vector (representing an amplitude) indicates the strength of edging, and the direction of the wavelet vector indicates the direction of the edging.
- the wavelet vector may simply be called a vector.
- stage 503 vectors are classified.
- the computer for classifying vectors in stage 503 is an example of a classifying means included in the present invention.
- Two kinds of thresholds are employed in classification of vectors.
- One of the thresholds is concerned with the magnitude of a vector, and the other is concerned with the regulation in the direction of a vector.
- Two different values are determined as the values of the constant K, whereby two different thresholds Tg and Ts are determined for the magnitude of a vector.
- the threshold Tg is an example of a first threshold employed in the present invention, while the threshold Ts is an example of a second threshold employed in the present invention.
- the threshold Ts is smaller than the threshold Tg.
- a threshold Tr is determined for the regulation in the direction of a vector. Assuming that a ratio of one vector to another vector is adopted as the measure of the regulation in the direction of a vector, the threshold Tr is determined for the ratio.
- the threshold Tr is an example of a third threshold employed in the present invention.
- the ratio of vectors is the ratio of a mean vector, which is a mean of a center vector, which is located in the center of a window defined with a matrix of 3 by 3, and two ambient vectors whose differences from the center vector are the smallest to a mean vector that is a mean of all vectors in the window.
- vectors are classified into four categories of strongly edged objects, weakly edged objects, radial lines, and noises, as shown in FIG. 6 .
- the radial line is an artifact.
- the thus classified vectors are processed in different manners associated with the categories in stage 505 .
- the computer for processing the vectors in stage 505 is an example of a processing means included in the present invention.
- stage 505 vectors classified into the category of strongly edged objects are kept, that is, left intact.
- the vectors classified into the category of strongly edged objects have the magnitudes and directions thereof left intact.
- Vectors classified into the category of weakly edged objects have the directions thereof filtered.
- the filtering of directions is achieved by sampling a combination of a center vector, which is located in the center of a window defined with a matrix of 3 by 3, and two ambient vectors the sum of whose differences from the center vector is the smallest, and replacing the center vector with a mean vector that is a mean of the two ambient vectors and the center vector.
- a center vector exists at the center position 0
- ambient vectors exist at eight ambient positions 1 to 8 .
- FIG. 7 ( 1 ) to FIG. 7 ( 28 ) twenty-eight combinations of the center vector and two ambient vectors are conceivable.
- a combination of the center vector and two ambient vectors the sum of whose differences from the center vector is the smallest is sampled, and the center vector is replaced with a mean vector that is a mean of the two ambient vectors and the center vector.
- a mean vector with which a center vector is replaced may be multiplied by a variable coefficient.
- the coefficient is made larger than 1, the clearness in the edges of the objects can be intensified.
- the coefficient is smaller than 1, the clearness in the edges of the objects can be lessened.
- Vectors classified into the categories of radial lines and noises are smoothed.
- the smoothing of vectors is achieved by replacing a center vector, which is located in the center of a window defined with a matrix of 3 by 3, with a mean vector that is a mean of all vectors in the window.
- FIG. 8 is a conceptual diagram of the smoothing of vectors. Assuming that vectors within a window defined with a matrix of 3 by 3 are defined as shown in FIG. 8( a ), a mean vector of all the vectors in the window is expressed as a blank vector shown in FIG. 8( b ). A center vector is replaced with the mean vector. The smoothing of vectors lessens radial lines and noises.
- the four kinds of vectors having been processed as mentioned above are inverse-wavelet-transformed in stage 507 .
- the computer for performing inverse wavelet transform in stage 507 is an example of an image production means included in the present invention.
- the inverse wavelet transform is inverse wavelet transform using a two-dimensional gradient form.
- the inverse wavelet transform restores a real-domain image.
- FIG. 9 conceptually shows the foregoing signal processing.
- an image signal is regarded as a one-dimensional signal.
- FIG. 9( a ) indicates the strength of an input image signal.
- the strengths of the components of the signal representing objects are relatively large, but the strengths of the components thereof representing noises and radial lines are relatively small.
- Wavelet transform provides a wavelet-domain image like the one shown in FIG. 9( b ).
- an image portion defined in one scale j is shown.
- the components of a signal representing the image portion whose strengths are larger than a threshold Tg represent strongly edged objects.
- the component of the signal whose strength is equal to or smaller than the threshold Tg and larger than a threshold Ts and whose regulation is larger than a threshold Tr represents a weakly edged object.
- the component of the signal whose strength is equal to or smaller than the threshold Tg and whose regulation is equal to or smaller than the threshold Tr represents a radial line.
- the components of the signal whose strengths are equal to or smaller than the threshold Ts represent noises.
- the signal is recomposed into a signal shown in FIG. 9( c ) through classification and processing.
- the strengths of the components of the signal representing the strongly edged object remain unchanged.
- the strength of the component thereof representing the weakly edged object is improved so that the weakly edged object will become sharp.
- the components thereof representing the radial line and noises are lessened.
- the signal is inverse-wavelet-transformed to finally represent a real-domain output image in the manner as shown in FIG. 9( d ).
- the strengths of the components of the resultant signal representing the output image which represent objects remain unchanged compared with those of the corresponding components of a signal representing an input image, but the strengths of the components thereof which represent a radial line and noises are lower than those of the corresponding components of the signal representing the input image. Namely, the radial line and noises in the image can be reduced without occurrence of a shift in a CT number, deformation of edges, or an abrupt change.
- a component of a wavelet-domain signal expressing a weakly edged object has the directions thereof filtered, the fine structure of the object becomes clear.
- FIG. 10 to FIG. 13 show examples of improvements in image quality derived from the foregoing image processing.
- FIG. 10 shows an original image produced by projecting a uniform phantom, a processed image, and a comparison of standard deviations sd of pixels at five positions in the original image with those of corresponding pixels in the processed image.
- a mean of the standard deviations in the processed image is 28% smaller. Thus, noises are reduced markedly.
- FIG. 11 shows an original image produced by projecting a low-contrast density (LCD) phantom, a processed image, and a comparison of contrast resolutions (LCDs) detected at three points in the original image with those detected in the processed image.
- the LCDs in the processed images are 1.4% larger than those in the original image.
- the contrast resolution is markedly improved due to noise reduction.
- FIG. 12 shows an original image produced by projecting a spatial-resolution phantom, a processed image, and a comparison of a spatial resolution (or a degree of line partition (LP)) in the original image with that in the processed image.
- the spatial resolution (LP) in the processed image is identical to that in the original image. This means that the spatial resolution is not degraded.
- FIG. 13 shows an original image (a) produced by projecting a Modulation Transfer Function (MTF) phantom, a processed image (b), and a comparison (c) of Modulation Transfer Functions describing the relationships between values detected in the original image and values detected in the processed image.
- the Modulation Transfer Functions are retained at specified values. This means that the Modulation Transfer Functions are not fluctuated.
Abstract
Description
T=Kσ2−(j−1)/2 [Formula 2]
where K denotes a constant, σ denotes a standard deviation of modules for 1st scale wavelet vectors in a area which is uniform image, and j denotes a scale.
Claims (10)
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200510109922.7 | 2005-09-15 | ||
CNB2005101099227A CN100548222C (en) | 2005-09-15 | 2005-09-15 | Image processing method and X ray CT system |
CN200510109922 | 2005-09-15 |
Publications (2)
Publication Number | Publication Date |
---|---|
US20070058883A1 US20070058883A1 (en) | 2007-03-15 |
US7668360B2 true US7668360B2 (en) | 2010-02-23 |
Family
ID=37855172
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/531,899 Active 2028-08-14 US7668360B2 (en) | 2005-09-15 | 2006-09-14 | Image processing method and X-ray CT system |
Country Status (3)
Country | Link |
---|---|
US (1) | US7668360B2 (en) |
JP (1) | JP4950607B2 (en) |
CN (1) | CN100548222C (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009087777A1 (en) * | 2008-01-11 | 2009-07-16 | Shimadzu Corporation | Image processing method, its device and laminagraph device |
CN101915739B (en) * | 2010-06-30 | 2012-11-14 | 李华 | Multi-dimensional chromatograph vector sorting method in water quality detection based on machine vision |
CN106447684B (en) * | 2016-08-10 | 2019-04-16 | 中国兵器科学研究院宁波分院 | Weak marginal dimension measurement method in INDUSTRIAL CT IMAGE |
CN108664980A (en) * | 2018-05-14 | 2018-10-16 | 昆明理工大学 | A kind of sun crown ring structure recognition methods based on guiding filtering and wavelet transformation |
US11842531B2 (en) * | 2020-02-28 | 2023-12-12 | United States Postal Service | System and method for image compression |
WO2023043036A1 (en) * | 2021-09-14 | 2023-03-23 | 엘지전자 주식회사 | Image processing method and apparatus |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5657085A (en) | 1993-02-18 | 1997-08-12 | Nec Corporation | Wavelet transform coding method |
US6141452A (en) * | 1996-05-13 | 2000-10-31 | Fujitsu Limited | Apparatus for compressing and restoring image data using wavelet transform |
US6233357B1 (en) * | 1997-07-14 | 2001-05-15 | Sharp Laboratories Of America, Inc. | Arbitrary shape wavelet transform with phase alignment |
US6361501B1 (en) | 1997-08-26 | 2002-03-26 | Seiko Epson Corporation | Pulse wave diagnosing device |
JP2002133399A (en) | 2000-10-23 | 2002-05-10 | Hitachi Medical Corp | Image processor and x-ray ct apparatus using it |
US6567081B1 (en) | 2000-01-21 | 2003-05-20 | Microsoft Corporation | Methods and arrangements for compressing image-based rendering (IBR) data using alignment and 3D wavelet transform techniques |
US6643406B1 (en) * | 1999-07-28 | 2003-11-04 | Polaroid Corporation | Method and apparatus for performing linear filtering in wavelet based domain |
US6658158B2 (en) * | 1998-12-29 | 2003-12-02 | Sony Corporation | Wavelet encoding method and apparatus and wavelet decoding method and apparatus |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3813999B2 (en) * | 1995-08-25 | 2006-08-23 | ジーイー横河メディカルシステム株式会社 | Image processing method and image processing apparatus |
JP2001111853A (en) * | 1999-10-13 | 2001-04-20 | Seiko Epson Corp | Color correction device |
JP4679710B2 (en) * | 2000-10-25 | 2011-04-27 | 富士フイルム株式会社 | Noise suppression processing apparatus and recording medium |
US6956975B2 (en) * | 2001-04-02 | 2005-10-18 | Eastman Kodak Company | Method for improving breast cancer diagnosis using mountain-view and contrast-enhancement presentation of mammography |
-
2005
- 2005-09-15 CN CNB2005101099227A patent/CN100548222C/en not_active Expired - Fee Related
-
2006
- 2006-09-12 JP JP2006246324A patent/JP4950607B2/en active Active
- 2006-09-14 US US11/531,899 patent/US7668360B2/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5657085A (en) | 1993-02-18 | 1997-08-12 | Nec Corporation | Wavelet transform coding method |
US6141452A (en) * | 1996-05-13 | 2000-10-31 | Fujitsu Limited | Apparatus for compressing and restoring image data using wavelet transform |
US6233357B1 (en) * | 1997-07-14 | 2001-05-15 | Sharp Laboratories Of America, Inc. | Arbitrary shape wavelet transform with phase alignment |
US6361501B1 (en) | 1997-08-26 | 2002-03-26 | Seiko Epson Corporation | Pulse wave diagnosing device |
US6658158B2 (en) * | 1998-12-29 | 2003-12-02 | Sony Corporation | Wavelet encoding method and apparatus and wavelet decoding method and apparatus |
US6643406B1 (en) * | 1999-07-28 | 2003-11-04 | Polaroid Corporation | Method and apparatus for performing linear filtering in wavelet based domain |
US6567081B1 (en) | 2000-01-21 | 2003-05-20 | Microsoft Corporation | Methods and arrangements for compressing image-based rendering (IBR) data using alignment and 3D wavelet transform techniques |
JP2002133399A (en) | 2000-10-23 | 2002-05-10 | Hitachi Medical Corp | Image processor and x-ray ct apparatus using it |
Non-Patent Citations (1)
Title |
---|
Stephane Mallat and Sifen Zhong; Characterization of Signals from Multiscale Edges; IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 14, No. 7, Jul. 1992; 710-732. |
Also Published As
Publication number | Publication date |
---|---|
US20070058883A1 (en) | 2007-03-15 |
JP2007080268A (en) | 2007-03-29 |
CN100548222C (en) | 2009-10-14 |
CN1931096A (en) | 2007-03-21 |
JP4950607B2 (en) | 2012-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7865006B2 (en) | Methods and apparatus for artifact reduction | |
US8009890B2 (en) | Image display apparatus and X-ray CT apparatus | |
US20050254721A1 (en) | Image processing method, image processing system, and X-ray CT system | |
US7444010B2 (en) | Method and apparatus for the reduction of artifacts in computed tomography images | |
JP4152649B2 (en) | Method and apparatus for CT scout image processing | |
US7978895B2 (en) | X-ray CT system | |
US6493416B1 (en) | Method and apparatus for noise reduction in computed tomographic systems | |
US7668360B2 (en) | Image processing method and X-ray CT system | |
US20070019851A1 (en) | Image processing apparatus and X-ray CT apparatus | |
EP1716809A1 (en) | Tomogram reconstruction method and tomograph | |
Chen et al. | Three-dimensional point spread function measurement of cone-beam computed tomography system by iterative edge-blurring algorithm | |
US20080019607A1 (en) | System and method for correcting for ring artifacts in an image | |
Abu Anas et al. | Comparison of ring artifact removal methods using flat panel detector based CT images | |
US8731267B2 (en) | X-ray computed tomography apparatus and image processing apparatus | |
US8798350B2 (en) | Method and system for reconstruction algorithm in cone beam CT with differentiation in one direction on detector | |
US6292526B1 (en) | Methods and apparatus for preprocessing volumetric computed tomography data | |
US5963614A (en) | Filters for single slice helical image reconstruction in a computed tomography system | |
JP6987352B2 (en) | Medical image processing equipment and medical image processing method | |
EP1468650B1 (en) | X-ray CT system and method with beam-hardening correction | |
US20070036263A1 (en) | Radiation CT method and X-ray CT apparatus | |
US8385620B2 (en) | Method and system for multi-detector-row computed tomography utilizing projection data up-sampling with shift | |
US6647084B1 (en) | Method and apparatus for filtering projection data of a helical scan | |
US6091840A (en) | Methods and apparatus for single slice helical image reconstruction in a computed tomography system | |
US6408088B1 (en) | Methods and apparatus for single slice helical image reconstruction in a computed tomography system | |
JP4607476B2 (en) | Radiographic image diagnostic apparatus and data processing method of radiographic image diagnostic apparatus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: GE HANGWEI MEDICAL SYSTEMS CO. LTD.,CHINA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:XING, ZHANFENG;REEL/FRAME:018253/0876 Effective date: 20050812 Owner name: GE HANGWEI MEDICAL SYSTEMS CO. LTD., CHINA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:XING, ZHANFENG;REEL/FRAME:018253/0876 Effective date: 20050812 Owner name: GE MEDICAL SYSTEMS GLOBAL TECHNOLOGY COMPANY, LLC, Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:GE HANGWEI MEDICAL SYSTEMS CO. LTD.;REEL/FRAME:018253/0900 Effective date: 20050819 |
|
FEPP | Fee payment procedure |
Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
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) 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 |