WO2021135339A1 - Ct图像中环形伪影的校正方法、装置及计算机程序介质 - Google Patents
Ct图像中环形伪影的校正方法、装置及计算机程序介质 Download PDFInfo
- Publication number
- WO2021135339A1 WO2021135339A1 PCT/CN2020/113215 CN2020113215W WO2021135339A1 WO 2021135339 A1 WO2021135339 A1 WO 2021135339A1 CN 2020113215 W CN2020113215 W CN 2020113215W WO 2021135339 A1 WO2021135339 A1 WO 2021135339A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- ring artifacts
- sinogram
- detection value
- images according
- value
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 121
- 238000012937 correction Methods 0.000 title claims abstract description 43
- 238000004590 computer program Methods 0.000 title claims abstract description 8
- 238000001514 detection method Methods 0.000 claims abstract description 136
- 238000012935 Averaging Methods 0.000 claims abstract description 59
- 238000012545 processing Methods 0.000 claims abstract description 29
- 238000007781 pre-processing Methods 0.000 claims abstract description 16
- 230000008569 process Effects 0.000 claims description 48
- 238000001914 filtration Methods 0.000 claims description 47
- 230000006870 function Effects 0.000 claims description 8
- 239000000654 additive Substances 0.000 claims description 6
- 230000000996 additive effect Effects 0.000 claims description 6
- 230000002238 attenuated effect Effects 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 20
- 230000015654 memory Effects 0.000 description 11
- 230000000694 effects Effects 0.000 description 7
- 238000002591 computed tomography Methods 0.000 description 5
- 108010001267 Protein Subunits Proteins 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000013135 deep learning Methods 0.000 description 2
- 238000003702 image correction Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 102000002067 Protein Subunits Human genes 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000006227 byproduct Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4233—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using matrix detectors
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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
Definitions
- the present invention relates to the field of CT image processing, and more particularly to a method, device and computer program medium for correcting ring artifacts in CT images based on sinograms.
- CT is a commonly used medical imaging equipment, which usually includes a scanning part, a computer system, and an image display and storage system.
- the scanning part includes an X-ray tube, a detector module, and a scanning frame.
- the X-ray tube emits an X-ray beam to scan the object.
- the detector module receives the X-rays transmitted through the scanned object and converts it into visible light, which is further converted into electrical signals through photoelectric conversion;
- the computer system is responsible for storing and computing the electrical signal data collected by the scanning; image display and storage system responsible for displaying the images processed and reconstructed by the computer system on a TV screen, multiple cameras or laser cameras, etc.
- the scanning principle of CT can be simply summarized as detecting the attenuation of X-rays by objects at multiple angles.
- the X-ray tube 1 emits X during CT scanning.
- the ray irradiates the scanned object 3, and the scanned object 3 is placed between the X-ray tube 1 and the X-ray detector 2 module.
- the X-ray detector module 2 includes a plurality of detectors 21 arranged in sequence ( Figure 2) These detectors are arranged regularly in rows or columns, and each row or column of detectors is called a channel.
- the X-ray detector module 2 can detect the intensity of X-rays at different positions that have been attenuated by the scanned object 3, and compare the intensity with the intensity of the X-rays that have not been attenuated to calculate the attenuation coefficient. Subsequently, the X-ray tube 1 and the X-ray detector module 2 rotate 360° around the z-axis at the same time to complete a scan.
- the attenuation coefficients at different angles and positions obtained by scanning can further form a projection map of each angle.
- the projection map can be expressed in the form of a sinogram, as shown in Figure 3.
- the sinogram is a form of CT raw data.
- Each line of the sinogram is the x direction of the detector under a certain scanning angle (as shown in Figure 1).
- the data collected by a row of detector channels.
- the data can be X-ray intensity, attenuation coefficient, or attenuation coefficient after slope filtering.
- Arrange the data of each row in the order of angle to get a two-dimensional matrix Each element in the two-dimensional matrix can be represented by I(s, ⁇ ), and ⁇ represents the scanning angle of the detector.
- the scanning angle can be based on the X-ray irradiation direction during the initial scan or other known directions.
- S represents a certain position in a row of detectors. Since the detector has multiple rows of channels in the x direction, each row of channels can obtain a sinogram after scanning multiple angles, so multiple sinograms can be obtained in one scan.
- the purpose of the present invention is to provide a method, device and computer program medium for correcting ring artifacts in CT images, so as to solve the problem that correction methods in the prior art cannot balance correction quality and image spatial resolution.
- Step S1 preprocessing the original detection data
- Step S2 Mark the bad tracks of the detector according to the data in the preprocessed original sinogram
- Step S3 Use the first averaging process to obtain the substitute detection value corresponding to the bad track
- Step S4 Perform CT image reconstruction using the sinogram data after the first average processing.
- I (s, ⁇ ) represents the actual detection value of the original sinogram at the scanning angle ⁇ and the detector position s
- I r (s, ⁇ ) represents the ideal sine at the scanning angle ⁇ and the detector position s Detection value
- b(s) represents the offset value of the actual detection value corresponding to the detector position s
- ⁇ (s, ⁇ ) represents the original sinogram corresponding to the scanning angle ⁇ and the detector position s Random error.
- the preprocessing includes:
- Step S12 filtering the original sinogram in the s direction of the detector position
- Step S15 Obtain the standard deviation of the noise in the ⁇ direction in the above step S13.
- the number of selected scanning angles ⁇ is not less than 50.
- the filtering process is a weighted filtering process.
- the weighted filtering process is an average filtering process
- the formula used in the average filtering process is as follows:
- the noise is obtained by subtracting the ideal detection value in the ideal sinogram after the filtering from the actual detection value in the original sinogram.
- the second averaging of the noise in the direction of the scanning angle ⁇ may be specifically performed by the following formula:
- the standard deviation of the noise at the detector position s is obtained by the following formula:
- N is the number of scanning angles ⁇ selected in the step S11.
- marking bad tracks is performed by judging whether the following formula holds:
- the correction coefficient d is 1.96.
- the specific step of the first averaging processing is: using the average value of the actual detection values of 2n detectors adjacent to the same straight line on both sides of the bad track as The substitute detection value of the bad track, n is a natural number.
- the specific step of the first averaging processing is: using four detectors adjacent to the bad track The average value of the actual detection value of the detector is used as the substitute detection value of the bad track.
- the specific steps of the first averaging process are: taking the direction perpendicular to the bad track
- the actual detection values of the eight detectors are averaged as the substitute detection value of the bad track.
- the present invention also provides a correction device for ring artifacts in CT images.
- the correction device includes: a preprocessing unit configured to preprocess the original detection data; a marking unit configured to perform preprocessing based on the preprocessed original The data of the sinogram marks the bad sectors; the first averaging processing unit is configured to adopt the first averaging process to obtain the substitute detection value corresponding to the bad sectors; and the image reconstruction unit is configured to adopt the first averaging process The sinogram data is reconstructed by CT image.
- the second averaging subunit is configured to output the offset value of the actual detection value after the second averaging.
- the first averaging process uses the average value of actual detection values of 2n detectors adjacent to the same straight line on both sides of the bad track as the substitute detection value of the bad track, and n is a natural number. .
- the first averaging process uses the average value of the actual detection values of the four detectors adjacent to the bad track as the bad track.
- the substitute detection value of the channel is used.
- the first averaging process uses the average value of the actual detection values of the eight detectors adjacent to the bad track as the bad track.
- the substitute detection value of the channel is used.
- the present invention also provides a computer program medium with program instructions stored on the computer storage medium.
- the program instructions When executed, the following functions are realized: preprocessing the original sinogram; marking from the preprocessed original sinogram A bad track is generated; the first averaging process is used to obtain the corresponding substitute detection value at the bad track; the data after the first averaging process is used for CT image reconstruction.
- the analysis of the additive noise model and the method based on the hypothesis test in the statistical principle are used to detect the bad tracks, which can effectively identify the bad tracks. If there are more obvious bad sectors, the detection value of the bad sectors is replaced with the average of the detection values of the adjacent channels, so as to eliminate the more obvious ring artifacts in the image.
- the ring artifacts that are not particularly obvious in the image are usually caused by the inconsistent detection efficiency of the channel.
- the present invention can reduce the ring artifacts by subtracting the offset value of the channel to determine the detection efficiency.
- the present invention only uses the average value of adjacent channels to replace the bad tracks with a relatively small number, which reduces the loss of image spatial resolution to the greatest extent and ensures the image quality.
- the method, device and computer program medium for correcting ring artifacts in CT images require no training before correction, which can basically eliminate ring artifacts in CT images while ensuring that the spatial resolution of CT images is basically not lost. , It has positive significance for CT image correction, which is beneficial to cost saving.
- Fig. 2 is a schematic diagram of a CT X-ray detector module and bad sectors in the X-ray detector module in the prior art
- Fig. 3 is a schematic diagram of a sinogram obtained by CT scanning in the prior art
- FIG. 4 is a schematic diagram of steps of a method for correcting ring artifacts in CT images according to an embodiment of the present invention
- FIG. 7 is a schematic diagram of noise according to the method for correcting ring artifacts in the CT image of FIG. 6;
- FIG. 9 is a schematic diagram of the offset standard deviation according to the method for correcting ring artifacts in the CT image of FIG. 7;
- FIG. 11 is an average schematic diagram of one of the methods for correcting ring artifacts in CT images according to an embodiment of the present invention.
- FIG. 12 is an average schematic diagram of one of the methods for correcting ring artifacts in CT images according to an embodiment of the present invention.
- 15 is a schematic diagram of the correction effect of a method for correcting ring artifacts in a CT image according to an embodiment of the present invention.
- 16 is a schematic diagram of the correction effect of a method for correcting ring artifacts in a CT image according to an embodiment of the present invention
- Fig. 17 is a schematic diagram of the arrangement of an apparatus for correcting ring artifacts in CT images according to an embodiment of the present invention.
- connection/connection may include electrical and/or mechanical physical connection/connection.
- the term “include/include” refers to the existence or addition of features, steps or components/parts, but does not exclude the existence or addition of one or more other features, steps or components/parts.
- the term “and/or” as used herein includes any and all combinations of one or more of the associated listed items.
- the present invention needs to first establish an additive noise model for CT scan detection:
- I(s, ⁇ ) represents the actual detection value of the original sinogram obtained by the CT scan at the scanning angle ⁇ and the detector position s.
- the scanning angle ⁇ and the detector position s can be customized according to different application scenarios.
- the initial irradiation direction of the X-ray generator can be selected as the reference direction of the scanning angle ⁇ , and the detector position s can be marked according to the array arrangement of the detectors in the X-ray detector module;
- I r (s, ⁇ ) means in The ideal detection value of the ideal sinogram at the scanning angle ⁇ and the detector position s;
- b(s) represents the offset value of the actual detection value of the corresponding detector channel at the detector position s, which has nothing to do with the angle , Which is only related to whether the detector is damaged and the output signal effect of the detector itself, is the root cause of the ring artifact;
- ⁇ (s, ⁇ ) represents the random error corresponding to the original sinogram at the scanning angle ⁇ and the detector position s
- Step S1 preprocessing the original detection data
- Step S2 Mark the bad tracks of the detector according to the data in the preprocessed original sinogram
- Step S3 Use the first averaging process to obtain the substitute detection value corresponding to the bad track
- Step S4 Perform CT image reconstruction using the sinogram data after the first average processing.
- step S1 further includes the following steps:
- Step S12 filtering the original sinogram formed in step S11 in the direction of the detector position s;
- Step S14 Perform a second average on the noise obtained in the above step S13 in the ⁇ direction;
- Step S15 Obtain the standard deviation of the noise at the detector position s in the above step S13.
- the attenuation coefficient is the ratio of the intensity of the attenuated X-ray detected by the detector to the intensity of the initial X-ray emitted by the X-ray.
- the sinogram formed by step S11 is shown in Figure 3, where the horizontal represents the detector position s, and the vertical represents the scanning angle ⁇ .
- the black horizontal line in Figure 3 represents the position of each detector s at a certain scanning angle ⁇ .
- the filtering processing of the original sinogram in the s direction may be average filtering, and the average filtering may be performed by the following formula:
- n is a natural number, that is, the average value of the adjacent 2n+1 actual detection values I(s, ⁇ ) can be used as the ideal detection value I r (s, ⁇ ) of the corresponding position.
- the above formula (2) is:
- the filtering process in step S12 may also be weighted filtering.
- Mean filtering is generally considered to be weighted filtering with equal weights.
- Weighted filtering may also be Gaussian weighted filtering or other types with high center weight and low weight on both sides. Weighted filtering, which is easily understood by those skilled in the art based on the technical enlightenment of the present invention combined with mathematical knowledge, will not be repeated here.
- the sinogram filtered through the above step S12 can be regarded as an ideal sinogram, that is , an ideal sinogram formed by I r (s, ⁇ ), as shown in Figure 6, but at this time, the filtered sinogram cannot be directly used for CT Image reconstruction, so this will lead to a great reduction in the spatial resolution of CT images. It should be understood by those skilled in the art that other filtering methods can also be used for processing in step S12, which will not be listed here. At the same time, the number of n can also be selected according to actual needs or spatial resolution requirements. No longer.
- the noise of the original sinogram is obtained by subtracting the ideal detection value data in the filtered ideal sinogram from the actual detection value data in the original sinogram, that is, by I(s, ⁇ )-I r (s, ⁇ ) Obtain the sum of the offset value of the detector output and the random error b(s) + ⁇ (s, ⁇ ), which is collectively referred to as noise, as shown in Figure 7.
- N is the number of scanning angles ⁇ selected in step S11, Represents the average offset value of all offset values b(s) in the direction of all scanning angles ⁇ .
- the average bias value of all ideal bias values b(s) can be used in this application instead of the ideal bias value b(s), that is, the average bias value Equal to the ideal offset value b(s) corresponding to each detection channel.
- the mean value of the random error ⁇ (s, ⁇ ) should be 0, that is, when the detector position s is fixed, the mean value of the random error ⁇ (s, ⁇ ) in the direction of the scanning angle ⁇ should be Is 0, therefore, the average offset value of the detector output can be obtained after the second average
- Figure 8 shows the average offset value corresponding to the detector position s after the second average
- the abscissa represents the position of the detector s
- the ordinate represents the average offset value
- step S15 the standard deviation of the noise at the detector position s is obtained by the following formula:
- N is the number of scanning angles ⁇ selected in step S11
- the standard deviation std[ ⁇ (s)] of the error ⁇ (s, ⁇ ) which is essentially a random error on the detector position s, is only related to the detector position s and has nothing to do with the scanning angle ⁇ , as shown in Figure 9 , Where the abscissa represents the detector position s, and the ordinate represents the standard deviation std[ ⁇ (s)].
- step S2 marking bad tracks is performed by judging whether the following formula holds:
- Means the average bias value The absolute value of, d represents the correction coefficient, sqrt(N) represents the root mean square of N. Since the bias value b(s) is the root cause of ring artifacts, when the average bias value Too large will cause ring artifacts, so it is necessary to determine the average offset value of which channels is too big. For each detector position s, judge the average offset value according to whether formula (6) is valid Whether it is too large, if formula (6) is established, determine the average offset value It is too large, and the channel corresponding to the detector position s is marked as a bad track. The marked bad track is shown in Figure 2, where the squares with different gray values are all bad tracks.
- the average bias value can be known through the central limit theorem. Satisfy the Gaussian distribution.
- the offset value b(s) should be 0, which is the calculated average offset value
- the average value of is 0, due to the existence of error, the average offset value In fact, it satisfies a Gaussian distribution with a mean of 0 and a standard deviation of std[ ⁇ (s)]ssqrt(N). In this Gaussian distribution, the probability of abs[b(s)] ⁇ d ⁇ std[ ⁇ (s)] ⁇ sqrt(N) is relatively large.
- the correction coefficient is between 1.6-2.6; when the correction coefficient d is between 1.6-1.8, almost all ring artifacts can be removed by the correction method of the present invention, but the spatial resolution will be slightly Decrease; when the correction coefficient d is between 1.8-2.0 (preferably 1.96), the correction method of the present invention can remove almost all ring artifacts, while having almost no impact on the spatial resolution; when the correction coefficient d is between When it is between 2.0 and 2.6, most of the ring artifacts can be removed by the correction method of the present invention, while the spatial resolution is hardly affected.
- the specific method of using the first averaging process to obtain the substitute detection value of the detector corresponding to the bad track is as follows:
- the actual detection values of the detectors corresponding to the 2n channels adjacent to the bad track in the s direction are averaged as the substitute detection value corresponding to the bad track, for example,
- the actual detection values of all eight detectors in the row except the bad track can be averaged as the alternative detection of the bad track.
- n is 4, and the actual detection values of the four detectors are taken on both sides of the bad track; the actual detection values of the four adjacent detectors on the left and right sides of the row can also be averaged as the bad track.
- the substitute detection value of the track, at this time n is 2, and the actual detection value of two detectors is taken on both sides of the bad track.
- the actual detection values of the detectors corresponding to or adjacent to the bad track are averaged as the substitute detection value corresponding to the bad track.
- adjacent in this application refers to the detector adjacent to the bad track, that is, there is no other detector between the bad track detector and the adjacent detector; "near” or “nearby” means the bad track There may be other detectors near the detector, that is, between the detector at the bad track and the adjacent or neighboring detector.
- the actual detection values of eight detectors distributed around the bad track can be averaged as the substitute detection value of the bad track; in Figure 12, the four detectors adjacent to the top, bottom, left, and right of the bad track can be taken.
- the actual detection value of the bad track is averaged as the substitute detection value of the bad track; in Figure 13, the actual detection value of the eight detectors along the perpendicular orthogonal direction to the bad track can be averaged as the substitute detection value of the bad track; in the figure In 14, the actual detection values of twelve detectors adjacent to the bad track can be averaged as the substitute detection value of the bad track.
- the approximate value obtained by subtracting the average offset value from the actual detection value of each detector may be used.
- the first averaging process can use the average value of the approximate values of 2n detectors adjacent to both sides of the bad track as a replacement for the bad track Detection value.
- the weighted actual detection value may also be used as the basis of the first averaging process, which will not be repeated here.
- step S4 the CT image reconstruction using the data after the first averaging process is consistent with the image reconstruction method in the prior art, and will not be repeated here.
- FIG. 15 is a schematic diagram of the correction effect of a method for correcting ring artifacts in CT images according to an embodiment of the present invention.
- the part indicated by the arrow in the left figure is the ring artifact, and the arrow in the right figure is the elimination of the ring artifact.
- the ring artifact is basically completely eliminated by the present invention and the spatial resolution of the CT image is hardly affected.
- FIG. 16 is a schematic diagram of the correction effect of a method for correcting ring artifacts in CT images according to another embodiment of the present invention.
- the part indicated by the arrow in the left figure is the ring artifact, and the arrow in the right figure is After eliminating the ring artifact, the CT image of the corresponding area obviously eliminates the ring artifact completely and the spatial resolution of the CT image is almost unaffected by the present invention.
- An embodiment of the present application also provides a correction device for ring artifacts in CT images.
- the correction device at least includes:
- the preprocessing unit 110 can be configured to preprocess the original detection data and can output the offset value of the actual detection value of each position detector;
- the marking unit 120 may be configured to mark the bad tracks of the detector according to the data in the preprocessed original sinogram
- the first averaging processing unit 130 may be configured to adopt the first averaging process to obtain the substitute detection value corresponding to the bad track, or be configured to subtract the deviation from the actual detection value of each detector before the first averaging process. Set value;
- the image reconstruction unit 140 may be configured to perform CT image reconstruction using the sinogram data after the first average processing.
- the above-mentioned pre-processing unit may further include:
- the data acquisition subunit 111 may be configured to select actual detection values corresponding to multiple scanning angles from the original detection data to form an original sinogram;
- the filtering subunit 112 may be configured to perform filtering processing on the original sinogram in the direction of the detector position s;
- the noise acquisition sub-unit 113 which may be configured to acquire the noise of the original sinogram through filtering processing
- the second averaging subunit 114 may be configured to perform a second average of the noise in the direction of the scanning angle ⁇ , and can output the offset value of the actual detection value of each position detector;
- the standard deviation obtaining sub-unit 115 may be configured to be able to obtain the standard deviation of the noise in the direction of the scanning angle ⁇ .
- the correction device may be a server, an electronic device, etc., or may be any device that can perform data processing based on a sinogram, which is not limited here.
- the functions implemented by the units in the above-mentioned correction device may also be executed by the processor in the computer according to the program instructions stored in the memory.
- the above-mentioned program instructions can be executed by a processor or by other processing devices.
- Volatile memory may include random access memory (RAM) or external cache memory.
- RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous chain Channel (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.
- the devices, units, sub-units, etc. described in the foregoing embodiments may be specifically implemented by computer chips and/or entities, or implemented by products with certain functions.
- the functions are divided into various units and described separately.
- the functions of each unit can be implemented in the same or multiple computer chips.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Pulmonology (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
一种CT图像中环形伪影的校正方法、装置及计算机程序介质,校正方法包括:对原始探测数据进行预处理(S1);根据预处理后的原始正弦图中的数据标记出坏道(S2);采用第一平均处理获取坏道处对应的替代探测值(S3);采用第一平均处理后的正弦图数据进行CT图像重建(S4)。装置包括:被配置为对原始探测数据进行预处理的预处理单元(110),被配置为根据预处理的原始正弦图的数据标记出坏道的标记单元(120),被配置为采用第一平均处理获取坏道处对应的替代探测值的第一平均处理单元(130),以及图像重建单元(140)。该校正方法和装置能够有效地消除CT图像中的环形伪影,同时保证CT图像的空间分辨率基本没有损失。
Description
本公开要求于2019年12月30日提交的中国专利申请201911392508.X的优先权,其全部内容通过援引加入本文。
本发明涉及CT图像处理领域,更具体地涉及一种基于正弦图的CT图像中环形伪影的校正方法、装置及计算机程序介质。
CT是一种常用的医疗成像设备,通常包括扫描部分、计算机系统以及图像显示和存储系统,扫描部分包括X射线管、探测器模块和扫描架等,X射线管发出X射线束以对扫描对象进行扫描,探测器模块接收透过扫描对象的X射线并将其转变为可见光,进一步经过光电转换变为电信号;计算机系统负责对扫描收集到的电信号数据进行贮存运算;图像显示和存储系统负责将经计算机系统处理、重建的图像通过电视屏、多幅照相机或激光照相机等进行显示。具体地,如图1所示,CT的扫描原理可以简单概括为探测多个角度下物体对X射线的衰减情况,比如,对于锥形束CT,CT扫描时,X射线管1发射出的X射线照射被扫描物体3,被扫描物体3放置于X射线管1和X射线探测器2模块之间,通常情况下,X射线探测器模块2中包括多个依序排列的探测器21(图2),这些探测器按照行或者列进行规则排列,每一行或者每一列的探测器称为一个通道。X射线探测器模块2能探测不同位置处经过被扫描物体3衰减后的X射线的强度,该强度与未经过衰减的X射线的强度进行比较可计算出衰减系数。随后X射线管1和X射线探测器模块2同时绕z轴旋转360°完成一 次扫描。扫描所获得的不同角度、不同位置处的衰减系数可以进一步形成每个角度的投影图。该投影图可以通过正弦图的形式表示,如图3所示,正弦图是一种CT原始数据的表现形式,正弦图的每一行是某个扫描角度下探测器在x方向(如图1所示)一行探测器通道采集到的数据,该数据可以是X射线强度,也可以是衰减系数,或者经过斜坡滤波后的衰减系数,按照角度顺序对每一行的数据进行排列可以得到一个二维矩阵,二维矩阵中每一个元素可以用I(s,θ)表示,θ表示探测器的扫描角度,该扫描角度可以初始扫描时X射线的照射方向为基准,也可以以其它已知方向为基准,s表示一行探测器中的某一位置。由于探测器在x方向有多行通道,每一行通道在扫描多个角度后都可以得到一个正弦图,因此一次扫描可得到多个正弦图。
当某个通道无法正常读取信号或者无法读出正确的信号时,即探测器21损坏或存在不均匀的响应时,该通道称为坏道,如图2中不同灰度值的方块22所示,当探测器损坏而无法正常读取信号时,灰度值相对较深;当探测器无法读出正确的信号时,灰度值相对较浅。由于CT扫描时X射线管和X射线探测器模块会一同旋转360度,因此当X射线探测器模块中的部分通道出现坏道或者响应不均匀时(以下统称为坏道),往往会在重建图像中出现环形伪影,这将直接影响阅片师对CT图像进行查阅的质量。然而,实际使用中,若仅仅因为存在少量的坏道而选择更换探测器模块,则成本过于高昂,因此,通常采用图像校正方法进行处理以消除坏道的影响。
目前CT图像的环形伪影的校正方法可以分为基于图像的环形伪影校正、基于正弦图的环形伪影校正以及基于深度学习的环形伪影校正,其中,基于深度学习的环形伪影校正方法效果较好,但是该方法需要大量的医学图像作为训练样本,然而,医学图像样本由于各方面 原因比较难于获取,因此该方法难于推广应用。对于其它两种校正方法,则存在一个难以调和矛盾:若要使图像不损失过多的空间分辨率,则较弱的环形伪影就无法彻底去除;若要使较弱的环形伪影去除干净,则会使得图像的空间分辨率损失较大。
因此,有必要开发一种CT图像中环形伪影的校正方法以克服上述技术问题。
发明内容
本发明的目的是提供一种CT图像中环形伪影的校正方法、装置及计算机程序介质,从而解决现有技术中校正方法无法兼顾校正质量和图像空间分辨率的问题。
本发明提供一种CT图像中环形伪影的校正方法,该校正方法包括:
步骤S1:对原始探测数据进行预处理;
步骤S2:根据预处理后的原始正弦图中的数据标记出探测器的坏道;
步骤S3:采用第一平均处理获取所述坏道处对应的替代探测值;
步骤S4:采用所述第一平均处理后的正弦图数据进行CT图像重建。
根据本发明的一个实施例,所述校正方法还包括在所述步骤S1之前建立加性噪声模型:
I(s,θ)=I
r(s,θ)+b(s)+ε(s,θ),
其中,I(s,θ)表示原始正弦图在扫描角度θ、探测器位置s处的实际探测值,I
r(s,θ)表示理想正弦图在扫描角度θ、探测器位置s处的理想探测值;b(s)表示在探测器位置s处对应的所述实际探测值的偏置值,ε(s,θ)表示所述原始正弦图在扫描角度θ、探测器位置s处对 应的随机误差。
根据本发明的一个实施例,在所述步骤S1中,所述预处理包括:
步骤S11:从所述原始探测数据中选取多个扫描角度θ对应的实际探测值形成原始正弦图;
步骤S12:对所述原始正弦图在探测器位置s方向进行滤波处理;
步骤S13:通过所述滤波处理获取所述原始正弦图的噪声;
步骤S14:对所述噪声在扫描角度θ方向进行第二平均;
步骤S15:获取上述步骤S13中噪声在θ方向的标准差。
根据本发明的一个实施例,在所述步骤S11中,选取的扫描角度θ的数量不低于50个。
根据本发明的一个实施例,在所述步骤S11中,所述实际探测值为探测器测量所得的经过衰减后的X射线的强度、衰减系数或者经过斜坡滤波后的衰减系数。
根据本发明的一个实施例,在所述步骤S12中,所述滤波处理为加权滤波处理。
根据本发明的一个实施例,在所述步骤S12中,所述加权滤波处理为均值滤波处理,均值滤波处理采用的公式如下:
其中,n为自然数。
根据本发明的一个实施例,在所述步骤S13中,所述噪声通过将原始正弦图中的实际探测值减去经过所述滤波后的理想正弦图中的理想探测值获得。
根据本发明的一个实施例,在所述步骤S14中,对所述噪声在扫描角度θ方向进行第二平均具体可以通过以下公式进行:
根据本发明的一个实施例,在所述步骤S15中,所述噪声在探测器位置s处的标准差通过以下公式获得:
其中,N为在所述步骤S11中选取的扫描角度θ的数量。
根据本发明的一个实施例,在所述步骤S2中,标记坏道通过判断以下公式是否成立进行:
根据本发明的一个实施例,所述校正系数d介于1.6-2.6、1.8-2.0或者2.0-2.6之间。
根据本发明的一个实施例,所述校正系数d为1.96。
根据本发明的一个实施例,在所述步骤S3中,所述第一平均处理的具体步骤为:采用所述坏道两侧同一直线方向临近的2n个探测器的实际探测值的平均值作为所述坏道的替代探测值,n为自然数。
根据本发明的一个实施例,在所述步骤S3中,当CT探测器的排列为二维排布时,所述第一平均处理的具体步骤为:采用与所述坏道邻接的四个探测器的实际探测值的平均值作为所述坏道的替代探测值。
根据本发明的一个实施例,在所述步骤S3中,当CT探测器的排列为二维排布时,所述第一平均处理采用与所述坏道邻近的八个探测器的实际探测值的平均值作为所述坏道的替代探测值。
根据本发明的一个实施例,在所述步骤S3中,当CT探测器的排列为二维排布时,所述第一平均处理的具体步骤为:取沿与所述坏道垂直正交方向的八个探测器的实际探测值进行平均作为所述坏道的替代探测值。
根据本发明的一个实施例,在所述步骤S3中,当CT探测器的排列为二维排布时,所述第一平均处理的具体步骤为:取与所述坏道邻近的十二个探测器的实际探测值进行平均作为坏道的替代探测值。
根据本发明的一个实施例,在所述第一平均处理中,各个探测器的实际探测值采用减去对应的平均偏置值的实际探测值替代。
根据本发明的一个实施例,在所述步骤S4的所述正弦图数据中,仅有所述坏道处的数据采用所述替代探测值替代,非坏道处的正弦图数据与所述原始正弦图中保持一致。
根据本发明的一个实施例,在所述步骤S4的所述正弦图数据中,所述坏道处的数据采用所述替代探测值替代,非坏道处的正弦图数据采用所述原始正弦图减去所述平均偏置值后的数据进行替代。
本发明还提供一种CT图像中环形伪影的校正装置,所述校正装置包括:预处理单元,其被配置为对原始探测数据进行预处理;标记单元,其被配置为根据预处理的原始正弦图的数据标记出坏道;第一平均处理单元,其被配置为采用第一平均处理获取坏道处所对应的替代探测值;以及图像重建单元,其被配置为采用第一平均处理后的正弦图数据进行CT图像重建。
根据本发明的一个实施例,所述预处理单元还包括:数据获取子单元,其被配置为从所述原始探测数据中选取多个扫描角度对应的实际探测值形成原始正弦图;滤波子单元,其被配置为对原始正弦图在探测器位置s方向进行滤波处理;噪声获取子单元,其被配置为通过 滤波处理获取所述原始正弦图的噪声;第二平均子单元,其被配置为对所述噪声在所述扫描角度方向进行第二平均;以及标准差获取子单元,其被配置为获取所述噪声在所述扫描角度方向的标准差。
根据本发明的一个实施例,所述第二平均子单元被配置为在所述第二平均后输出实际探测值的偏置值。
根据本发明的一个实施例,所述第一平均处理采用所述坏道两侧同一直线方向临近的2n个探测器的实际探测值的平均值作为所述坏道的替代探测值,n为自然数。
根据本发明的一个实施例,当CT探测器的排列为二维排布时,所述第一平均处理采用与所述坏道邻接的四个探测器的实际探测值的平均值作为所述坏道的替代探测值。
根据本发明的一个实施例,当CT探测器的排列为二维排布时,所述第一平均处理采用与所述坏道邻近的八个探测器的实际探测值的平均值作为所述坏道的替代探测值。
本发明还提供一种计算机程序介质,所述计算机存储介质上存储有程序指令,所述程序指令被执行时实现以下功能:对原始正弦图进行预处理;从预处理后的原始正弦图中标记出坏道;采用第一平均处理获取坏道处对应的替代探测值;采用第一平均处理后的数据进行CT图像重建。
在本发明中,由于图像中较为明显的环形伪影通常是通过坏道产生的,通过对加性噪声模型的分析,用基于统计原理中的假设检验的方法对坏道进行检测,能有效识别出较为明显的坏道,用邻近的通道的探测值的均值对坏道的探测值进行替换,从而消除图像中较为明显的环形伪影。图像中不是特别明显的环形伪影通常是由于通道探测效率不一致所产生的,本发明通过减去该通道的偏置值以判断其探测效 率,起到抑制环形伪影的效果。本发明仅对数量占比较少的坏道运用临近通道的均值取代,在最大程度上减少了图像空间分辨率的损失,保证了图像质量。
总之,本发明提供的CT图像中环形伪影的校正方法、装置及计算机程序介质,校正前无需进行训练,能够基本消除CT图像中的环形伪影,同时保证CT图像的空间分辨率基本没有损失,对于CT图像校正具有积极的意义,有利于节约成本。
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是现有技术中CT的扫描原理示意图;
图2是现有技术中CT的X射线探测器模块以及X射线探测器模块中的坏道示意图;
图3是现有技术中进行CT扫描所获得的正弦图的示意图;
图4是根据本发明一个实施例的CT图像中环形伪影的校正方法的步骤示意图;
图5是根据本发明一个实施例的CT图像中环形伪影的校正方法的步骤S1示意图;
图6根据图5的CT图像中环形伪影的校正方法的校正结果示意图,其中正弦图经过了滤波处理;
图7根据图6的CT图像中环形伪影的校正方法的噪声示意图;
图8是根据图7的CT图像中环形伪影的校正方法的偏置示意 图;
图9是根据图7的CT图像中环形伪影的校正方法的偏置标准差示意图;
图10是根据本发明一个实施例的CT图像中环形伪影的校正方法的平均示意图;
图11是根据本发明一个实施例的CT图像中环形伪影的校正方法的其中一种平均示意图;
图12是根据本发明一个实施例的CT图像中环形伪影的校正方法的其中一种平均示意图;
图13是根据本发明一个实施例的CT图像中环形伪影的校正方法的其中一种平均示意图;
图14是根据本发明一个实施例的CT图像中环形伪影的校正方法的其中一种平均示意图;
图15是根据本发明一个实施例的CT图像中环形伪影的校正方法的校正效果示意图;
图16是根据本发明一个实施例的CT图像中环形伪影的校正方法的校正效果示意图;
图17为根据本发明一个实施例的CT图像中环形伪影的校正装置的布置示意图。
以下结合具体实施例,对本发明做进一步说明。应理解,以下实施例仅用于说明本发明而非用于限制本发明的范围。
需要说明的是,当部件/零件被称为“设置在”另一个部件/零件上,它可以直接设置在另一个部件/零件上或者也可以存在居中的部件/零件。当部件/零件被称为“连接/联接”至另一个部件/零件,它可以是直 接连接/联接至另一个部件/零件或者可能同时存在居中部件/零件。本文所使用的术语“连接/联接”可以包括电气和/或机械物理连接/联接。本文所使用的术语“包括/包含”指特征、步骤或部件/零件的存在,但并不排除一个或更多个其它特征、步骤或部件/零件的存在或添加。本文所使用的术语“和/或”包括一个或多个相关所列项目的任意的和所有的组合。
除非另有定义,本文所使用的所有的技术和科学术语与属于本申请的技术领域的技术人员通常理解的含义相同。本文中所使用的术语只是为了描述具体实施例的目的,而并不是旨在限制本申请。
另外,在本发明的描述中,术语“第一”、“第二”等仅用于描述目的和区别类似的对象,两者之间并不存在先后顺序,也不能理解为指示或暗示相对重要性。此外,在本发明的描述中,除非另有说明,“多个”的含义是两个或两个以上。
在进行CT图像环形伪影的校正之前,本发明需要首先建立CT扫描探测的加性噪声模型:
I(s,θ)=I
r(s,θ)+b(s)+ε(s,θ) (式1)
其中,I(s,θ)表示CT扫描所获取的原始正弦图在扫描角度θ、探测器位置s处的实际探测值,扫描角度θ和探测器位置s具体可以根据不同的应用场景进行自定义,通常可以选择X射线发生器的初始照射方向为扫描角度θ的基准方向,探测器位置s可以根据X射线探测器模块中探测器的行列排布进行标记;I
r(s,θ)表示在理想正弦图在扫描角度θ、探测器位置s处的理想探测值;b(s)表示在探测器位置s处对应的探测器通道的实际探测值的偏置值,该偏置值与角度无关,仅与探测器是否损坏以及探测器自身的输出信号效果有关,是产生环形伪影的根源;ε(s,θ)表示原始正弦图在扫描角度θ、探测器位 置s处对应的随机误差,该随机误差的期望值应当为0。
图4是根据本发明一个实施例的CT图像中环形伪影的校正方法的步骤示意图,由图4可知,本发明提供的CT图像环形伪影的校正方法至少包括以下步骤:
步骤S1:对原始探测数据进行预处理;
步骤S2:根据预处理后的原始正弦图中的数据标记出探测器的坏道;
步骤S3:采用第一平均处理获取坏道处所对应的替代探测值;
步骤S4:采用第一平均处理后的正弦图数据进行CT图像重建。
具体地,如图5所示,上述步骤S1进一步包括以下步骤:
步骤S11:从原始探测数据中选取多个扫描角度θ对应的实际探测值形成原始正弦图;
步骤S12:对步骤S11中形成的原始正弦图在探测器位置s方向进行滤波处理;
步骤S13:通过滤波处理获取原始正弦图的噪声;
步骤S14:对上述步骤S13中获得的噪声在θ方向进行第二平均;
步骤S15:获取上述步骤S13中噪声在探测器位置s处的标准差。
具体地,在上述步骤S11中,原始正弦图通过选取多个扫描角度θ对应的探测器的实际探测值形成,即I(s,θ),通常情况下,扫描角度θ的数量不低于50个。在原始正弦图中,实际探测值可以是经过衰减后的X射线的强度,可以是衰减系数,也可以是经过斜坡滤波后的衰减系数。本领域技术人员应当理解的是,斜坡滤波为本领域常用的一种滤波处理方法,其可以是在频域空间乘以一个斜坡函数,斜坡滤波包括标准斜坡滤波、加窗斜坡滤波,也可以采用卷积法代替频域空 间的斜坡滤波,在此不再赘述。衰减系数为探测器所探测到的经过衰减后的X射线的强度与X射线发出的初始X射线的强度的比值。通过步骤S11所形成的正弦图如图3所示,其中横向表示探测器位置s,纵向表示扫描角度θ,比如,图3中黑色横线表示在某一扫描角度θ下各个探测器位置s处对应的探测器的实际探测值I(s,θ)。
进一步地,在上述步骤S12中,对原始正弦图在s方向进行滤波处理可以是均值滤波,均值滤波可以通过以下公式进行:
其中,n为自然数,即,可以采用临近的2n+1个实际探测值I(s,θ)的均值作为对应位置的理想探测值I
r(s,θ),比如,可以选取某一通道临近的3个(n=1时)实际探测值的均值作为该通道的理想探测值,此时上述公式(2)为:
I
r(s,θ)=[I(s-1,θ)+I(s,θ)+I(s+1,θ)]/3 (式3)
本领域技术人员应当理解的是,步骤S12中的滤波处理还可以是加权滤波,均值滤波通常认为是权重相等的加权滤波,加权滤波还可以是高斯加权滤波或者其它中心权重高两侧权重低的加权滤波,这属于本领域技术人员根据本发明的技术启示结合数学知识容易理解的,在此不再赘述。
通过上述步骤S12进行滤波后的正弦图可以视为理想正弦图,即I
r(s,θ)形成的理想正弦图,如图6所示,但是此时不能直接采用滤波后的正弦图进行CT图像重建,因此这将导致CT图像的空间分辨率极大的降低。本领域技术人员应当理解的是,步骤S12中还可以采用其它的滤波方法进行处理,在此不再一一列举,同时n的数量也可以根据实际需要或者空间分辨率要求等进行选择,在此不再赘述。
在上述步骤S13中,获取原始正弦图的噪声通过将原始正弦图中 的实际探测值数据减去滤波后的理想正弦图中的理想探测值数据获得,即通过I(s,θ)-I
r(s,θ)得到探测器输出的偏置值与随机误差之和b(s)+ε(s,θ),这统称为噪声,如图7所示。
在上述步骤S14中,对噪声在扫描角度θ方向进行第二平均具体可以通过以下公式进行:
其中,N为在步骤S11中选取的扫描角度θ的数量,
表示在所有扫描角度θ方向上所有偏置值b(s)的平均偏置值,根据本发明所建立的加性噪声模型,由于每一个通道的实际探测值I(s,θ)均存在一个理想的偏置值b(s),本申请中可以采用所有理想的偏置值b(s)的平均偏置值
代替理想的偏置值b(s),也就是可以使平均偏置值
等于每一个探测通道对应的理想的偏置值b(s)。
对于某一探测器通道而言,其随机误差ε(s,θ)的均值应当为0,即当探测器位置s固定时,随机误差ε(s,θ)在扫描角度θ方向上的均值应当为0,因此,经过第二平均后可以得到探测器输出的平均偏置值
如图8所示为探测器位置s对应的经过第二平均后的平均偏置值
其中横坐标表示探测器位置s,纵坐标表示平均偏置值
在上述步骤S15中,噪声在探测器位置s处的标准差通过以下公式获得:
其中,N为在步骤S11中选取的扫描角度θ的数量,
表示所有偏置值b(s)的平均偏置值,由于探测器输出的偏置值b(s)在θ方向为常量,因此在探测器位置s处求得的噪声的标准差即为随机误差ε(s,θ)的标准差std[ε(s)],其本质上是探测器位置s上的随机误差,仅与探测器位置s有关,与扫描角度θ无关,如图9所示,其中横坐标 表示探测器位置s,纵坐标表示标准差std[ε(s)]。
在上述步骤S2中,标记坏道通过判断以下公式是否成立进行:
其中,
表示平均偏置值
的绝对值,d表示校正系数,sqrt(N)表示N的均方根。由于偏置值b(s)是引起环形伪影的根源,当平均偏置值
过大时会引起环形伪影,因此需要判定哪些通道的平均偏置值
过大。对于每一个探测器位置s,根据公式(6)是否成立判断平均偏置值
是否过大,若公式(6)成立,则判定平均偏置值
过大,并将该探测器位置s对应的通道标记为坏道,经过标记后的坏道如图2所示,其中不同灰度值的方块均为坏道。
该判定方法的理论依据为假设检验,虽然偏置值b(s)的分布规律是未知的,但是通过中心极限定理可知平均偏置值
满足高斯分布。对于正常的探测器,其偏置值b(s)应当为0,即计算得到的平均偏置值
的均值为0,由于误差的存在,平均偏置值
实际上满足均值为0、标准差为std[ε(s)]ssqrt(N)的高斯分布。在该高斯分布中,abs[b(s)]<d×std[ε(s)]÷sqrt(N)的概率较大,比如当d=1.96时,abs[b(s)]<d×std[ε(s)]÷sqrt(N)的概率为95%,因此,满足公式(6)的概率只有5%,这是一个极小的概率,通常认为小概率事件的发生另有原因,即偏置值b(s)不为0,说明探测器某个位置s处出现了坏道。
进一步地,在上述公式(6)中,校正系数d对应的是一个用于判定是否为小概率事件的概率,比如,当d=1.96时,对应的认为上述公式(6)成立的概率仅有5%,此时认为5%为小概率事件。通常情况下,该校正系数介于1.6-2.6之间;当校正系数d介于1.6-1.8之间时,通过本发明的校正方法能够去除几乎所有的环形伪影,但是空间 分辨率会略有降低;当校正系数d介于1.8-2.0之间(优选地为1.96)时,通过本发明的校正方法能够去除几乎所有的环形伪影,同时对空间分辨率几乎无影响;当校正系数d介于2.0-2.6之间时,通过本发明的校正方法能够去除大部分的环形伪影,同时空间分辨率几乎不受影响。
根据本发明的一个实施例,在上述步骤S3中,当获取坏道的信息以后,采用第一平均处理获取坏道所对应的探测器的替代探测值的具体方法如下:
对于一维通道排列的CT探测器,如图10所示,取s方向与坏道相邻的2n个通道对应的探测器的实际探测值进行平均作为该坏道对应的替代探测值,比如在图10中,当一维排列的探测器为9个且中间一个探测器为坏道时,可以取该行中除了坏道外全部八个探测器的实际探测值进行平均作为该坏道的替代探测值,此时n为4,坏道两侧各取四个探测器的实际探测值;也可以取该行中除了坏道外左右两侧相邻的四个探测器的实际探测值进行平均作为坏道的替代探测值,此时n为2,坏道两侧各取两个探测器的实际探测值。本领域技术人员应当理解的是,对于二维通道,同样适用图10所示的平均方法。
对于二维通道排列的CT探测器,取与坏道邻接或者邻近的若干个通道对应的探测器的实际探测值进行平均作为该坏道对应的替代探测值,本领域技术人员应当理解的是,“邻接”在本申请中表示与坏道处探测器相邻的探测器,即坏道处探测器与邻接的探测器之间不存在其它探测器;“临近”或者“邻近”表示坏道处探测器附近的探测器,即坏道处探测器与临近或者邻近的探测器之间还可以存在其它的探测器。比如,在图11中,可以取沿坏道四周分布的八个探测器 的实际探测值进行平均作为坏道的替代探测值;在图12中,可以取坏道上下左右临近的四个探测器的实际探测值进行平均作为坏道的替代探测值;在图13中,可以取沿与坏道垂直正交方向的八个探测器的实际探测值进行平均作为坏道的替代探测值;在图14中,可以取与坏道邻近的十二个探测器的实际探测值进行平均作为坏道的替代探测值。
根据本发明的另一个实施例,在上述步骤S3中,在进行第一平均处理之前,可以采用每一个探测器的实际探测值减去平均偏置值后的近似值
作为第一平均处理的基础,比如,当CT探测器的排列为一维排布时,第一平均处理可以采用坏道两侧临近的2n个探测器的近似值的平均值作为该坏道的替代探测值。本领域技术人员应当理解的是,在进行第一平均处理时,还可以采用加权的实际探测值作为第一平均处理的基础,在此不再赘述。
在上述步骤S4中,采用第一平均处理后的数据进行CT图像重建与现有技术中进行图像重建的方法一致,在此不再赘述。
图15是根据本发明一个实施例的CT图像环形伪影的校正方法的校正效果示意图,在图15中,左侧图中箭头所指部分即为环形伪影,右侧图中箭头为消除环形伪影后对应区域的CT图像,显然通过本发明基本上完全消除了环形伪影并且CT图像的空间分辨率几乎不受影响。
图16是根据本发明另一个实施例的CT图像中环形伪影的校正方法的校正效果示意图,在图16中,左侧图中箭头所指部分即为环形伪影,右侧图中箭头为消除环形伪影后对应区域的CT图像,显然通过本发明基本上完全消除了环形伪影并且CT图像的空间分辨率几乎不受影响。
在本发明中,由于图像中较为明显的环形伪影通常是通过坏道产生的,通过对加性噪声模型的分析,用基于统计原理中的假设检验的方法对坏道进行检测,能有效识别出较为明显的坏道,用邻近的通道的探测值的均值对坏道的探测值进行替换,从而消除图像中较为明显的环形伪影。图像中不是特别明显的环形伪影通常是由于通道探测效率不一致所产生的,本发明通过减去该通道的偏置值以判断其探测效率,起到抑制环形伪影的效果。本发明仅对数量占比较少的坏道运用临近通道的均值取代,在最大程度上减少了图像空间分辨率的损失,保证了图像质量。
本申请实施例还提供了一种CT图像中环形伪影的校正装置,如图17所示,该校正装置至少包括:
预处理单元110,其可以被配置为对原始探测数据进行预处理,并且能够输出各个位置探测器的实际探测值的偏置值;
标记单元120,其可以被配置为根据预处理的原始正弦图中的数据标记出探测器的坏道;
第一平均处理单元130,其可以被配置为采用第一平均处理获取坏道处所对应的替代探测值,或者被配置为在第一平均处理前先将每一个探测器的实际探测值减去偏置值;以及
图像重建单元140,其可以被配置为采用第一平均处理后的正弦图数据进行CT图像重建。
进一步地,上述预处理单元还可以包括:
数据获取子单元111,其可以被配置为从原始探测数据中选取多个扫描角度对应的实际探测值形成原始正弦图;
滤波子单元112,其可以被配置为对原始正弦图在探测器位置s方向进行滤波处理;
噪声获取子单元113,其可以被配置为通过滤波处理获取原始正弦图的噪声;
第二平均子单元114,其可以被配置为对噪声在扫描角度θ方向进行第二平均,并且能够输出各个位置探测器的实际探测值的偏置值;以及
标准差获取子单元115,其可以被配置为能够获取噪声在扫描角度θ方向的标准差。
关于各个单元、子单元的具体描述,可以参照上述方法实施例中对步骤S1~S4以及步骤S11~S15的描述,在此不再赘述。
该校正装置可以是服务器、电子设备等,也可以是能够基于正弦图进行数据处理的任意装置,在此并不限制。另外,需要说明的是,上述校正装置中各单元所实现的功能也可以由计算机中的处理器根据存储器中所存储的程序指令来执行。
本发明还提供了一种计算机存储介质,该计算机存储介质上存储有程序指令,该程序指令被执行时可以实现以下功能:对原始正弦图进行预处理,并且输出各个位置探测器的实际探测值的偏置值;从预处理后的原始正弦图中标记出坏道;采用第一平均处理获取坏道所对应的替代探测值,或者先将探测器输出的实际探测值减去偏置值,再采用第一平均处理获取坏道所对应的替代探测值;采用第一平均处理后的数据进行CT图像重建。
上述程序指令可以由处理器执行,也可以由其他处理装置执行。
本领域技术人员应当理解的是,实现上述方法实施例中的全部或部分步骤,可以通过计算机程序来指令相关的硬件完成,所述的程序可存储于一非易失性计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,本发明所提供的各实施例 中所使用的对存储器、存储介质、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
上述实施例阐明的装置、单元、子单元等,具体可以由计算机芯片和/或实体实现,或者由具有某种功能的产品来实现。为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个计算机芯片中实现。
以上所述的,仅为本发明的较佳实施例,并非用以限定本发明的范围,本发明的上述实施例还可以做出各种变化。即凡是依据本发明申请的权利要求书及说明书内容所作的简单、等效变化与修饰,皆落入本发明专利的权利要求保护范围。本发明未详尽描述的均为常规技术内容。
Claims (28)
- 一种CT图像中环形伪影的校正方法,其特征在于,所述校正方法包括:步骤S1:对原始探测数据进行预处理;步骤S2:根据预处理后的原始正弦图中的数据标记出探测器的坏道;步骤S3:采用第一平均处理获取所述坏道处对应的替代探测值;步骤S4:采用所述第一平均处理后的正弦图数据进行CT图像重建。
- 根据权利要求1所述的CT图像中环形伪影的校正方法,其特征在于,所述校正方法还包括在所述步骤S1之前建立加性噪声模型:I(s,θ)=I r(s,θ)+b(s)+ε(s,θ),其中,I(s,θ)表示原始正弦图在扫描角度θ、探测器位置s处的实际探测值,I r(s,θ)表示理想正弦图在扫描角度θ、探测器位置s处的理想探测值;b(s)表示在探测器位置s处对应的所述实际探测值的偏置值,ε(s,θ)表示所述原始正弦图在扫描角度θ、探测器位置s处对应的随机误差。
- 根据权利要求2所述的CT图像中环形伪影的校正方法,其特征在于,在所述步骤S1中,所述预处理包括:步骤S11:从所述原始探测数据中选取多个扫描角度θ对应的实际探测值形成原始正弦图;步骤S12:对所述原始正弦图在探测器位置s方向进行滤波处理;步骤S13:通过所述滤波处理获取所述原始正弦图的噪声;步骤S14:对所述噪声在扫描角度θ方向进行第二平均;步骤S15:获取上述步骤S13中噪声在θ方向的标准差。
- 根据权利要求3所述的CT图像中环形伪影的校正方法,其特征在于,在所述步骤S11中,选取的扫描角度θ的数量不低于50个。
- 根据权利要求3所述的CT图像中环形伪影的校正方法,其特征在于,在所述步骤S11中,所述实际探测值为探测器测量所得的经过衰减后的X射线的强度、衰减系数或者经过斜坡滤波后的衰减系数。
- 根据权利要求3所述的CT图像中环形伪影的校正方法,其特征在于,在所述步骤S12中,所述滤波处理为加权滤波处理。
- 根据权利要求3所述的CT图像中环形伪影的校正方法,其特征在于,在所述步骤S13中,所述噪声通过将原始正弦图中的实际探测值减去经过所述滤波后的理想正弦图中的理想探测值获得。
- 根据权利要求11所述的CT图像中环形伪影的校正方法,其特征在于,所述校正系数d介于1.6-2.6、1.8-2.0或者2.0-2.6之间。
- 根据权利要求11所述的CT图像中环形伪影的校正方法,其特征在于,所述校正系数d为1.96。
- 根据权利要求1所述的CT图像中环形伪影的校正方法,其特征在于,在所述步骤S3中,所述第一平均处理的具体步骤为:采用所述坏道两侧同一直线方向临近的2n个探测器的实际探测值的平均值作为所述坏道的替代探测值,n为自然数。
- 根据权利要求1所述的CT图像中环形伪影的校正方法,其特征在于,在所述步骤S3中,当CT探测器的排列为二维排布时,所述第一平均处理的具体步骤为:采用与所述坏道邻接的四个探测器的实际探测值的平均值作为所述坏道的替代探测值。
- 根据权利要求1所述的CT图像中环形伪影的校正方法,其特征在于,在所述步骤S3中,当CT探测器的排列为二维排布时,所述第一平均处理采用与所述坏道邻近的八个探测器的实际探测值的平均值作为所述坏道的替代探测值。
- 根据权利要求1所述的CT图像中环形伪影的校正方法,其特征在于,在所述步骤S3中,当CT探测器的排列为二维排布时,所述第一平均处理的具体步骤为:取沿与所述坏道垂直正交方向的八个探测器的实际探测值进行平均作为所述坏道的替代探测值。
- 根据权利要求1所述的CT图像中环形伪影的校正方法,其特征在于,在所述步骤S3中,当CT探测器的排列为二维排布时,所述第一平均处理的具体步骤为:取与所述坏道邻近的十二个探测器的实际探测值进行平均作为坏道的替代探测值。
- 根据权利要求14-18中任一项所述的CT图像中环形伪影的校正方法,其特征在于,在所述第一平均处理中,各个探测器的实际探测值采用减去对应的平均偏置值的实际探测值替代。
- 根据权利要求1所述的CT图像中环形伪影的校正方法,其特征在于,在所述步骤S4的所述正弦图数据中,仅有所述坏道处的数据采用所述替代探测值替代,非坏道处的正弦图数据与所述原始正弦图中保持一致。
- 根据权利要求9所述的CT图像中环形伪影的校正方法,其特征在于,在所述步骤S4的所述正弦图数据中,所述坏道处的数据采用所述替代探测值替代,非坏道处的正弦图数据采用所述原始正弦图减去所述平均偏置值后的数据进行替代。
- 一种CT图像中环形伪影的校正装置,其特征在于,所述校正装置包括:预处理单元,其被配置为对原始探测数据进行预处理;标记单元,其被配置为根据预处理的原始正弦图的数据标记出坏道;第一平均处理单元,其被配置为采用第一平均处理获取坏道处所 对应的替代探测值;以及图像重建单元,其被配置为采用第一平均处理后的正弦图数据进行CT图像重建。
- 根据权利要求22所述的CT图像中环形伪影的校正装置,其特征在于,所述预处理单元还包括:数据获取子单元,其被配置为从所述原始探测数据中选取多个扫描角度对应的实际探测值形成原始正弦图;滤波子单元,其被配置为对原始正弦图在探测器位置s方向进行滤波处理;噪声获取子单元,其被配置为通过滤波处理获取所述原始正弦图的噪声;第二平均子单元,其被配置为对所述噪声在所述扫描角度方向进行第二平均;以及标准差获取子单元,其被配置为获取所述噪声在所述扫描角度方向的标准差。
- 根据权利要求23所述的CT图像中环形伪影的校正装置,其特征在于,所述第二平均子单元被配置为在所述第二平均后输出实际探测值的偏置值。
- 根据权利要求22所述的CT图像中环形伪影的校正装置,其特征在于,所述第一平均处理采用所述坏道两侧同一直线方向临近的2n个探测器的实际探测值的平均值作为所述坏道的替代探测值,n为自然数。
- 根据权利要求22所述的CT图像中环形伪影的校正装置,其特征在于,当CT探测器的排列为二维排布时,所述第一平均处理采用与所述坏道邻接的四个探测器的实际探测值的平均值作为所述坏 道的替代探测值。
- 根据权利要求22所述的CT图像中环形伪影的校正装置,其特征在于,当CT探测器的排列为二维排布时,所述第一平均处理采用与所述坏道邻近的八个探测器的实际探测值的平均值作为所述坏道的替代探测值。
- 一种计算机程序介质,其特征在于,所述计算机存储介质上存储有程序指令,所述程序指令被执行时实现以下功能:对原始正弦图进行预处理;从预处理后的原始正弦图中标记出坏道;采用第一平均处理获取坏道处对应的替代探测值;采用第一平均处理后的数据进行CT图像重建。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US17/790,168 US20230088126A1 (en) | 2019-12-30 | 2020-09-03 | Method and Device of Correction of Ring Artifact in CT Image and Computer Program Medium |
EP20910772.1A EP4085843A4 (en) | 2019-12-30 | 2020-09-03 | CORRECTION METHOD AND APPARATUS FOR ANNULAR ARTIFACT IN COMPUTED TOMOGRAPHY (CT) IMAGE, AND COMPUTER PROGRAM MEDIUM |
JP2022540633A JP2023514944A (ja) | 2019-12-30 | 2020-09-03 | Ct画像内のリングアーチファクトの補正の方法及びデバイス並びにコンピュータプログラム媒体 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911392508.XA CN111053568B (zh) | 2019-12-30 | 2019-12-30 | Ct图像中环形伪影的校正方法、装置及计算机存储介质 |
CN201911392508.X | 2019-12-30 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2021135339A1 true WO2021135339A1 (zh) | 2021-07-08 |
Family
ID=70304510
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2020/113215 WO2021135339A1 (zh) | 2019-12-30 | 2020-09-03 | Ct图像中环形伪影的校正方法、装置及计算机程序介质 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20230088126A1 (zh) |
EP (1) | EP4085843A4 (zh) |
JP (1) | JP2023514944A (zh) |
CN (1) | CN111053568B (zh) |
WO (1) | WO2021135339A1 (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111053568B (zh) * | 2019-12-30 | 2021-10-08 | 苏州瑞派宁科技有限公司 | Ct图像中环形伪影的校正方法、装置及计算机存储介质 |
CN115797486B (zh) * | 2022-11-04 | 2023-06-30 | 赛诺威盛科技(北京)股份有限公司 | 一种ct图像线状伪影消除方法、装置及存储介质 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7676073B2 (en) * | 2006-08-29 | 2010-03-09 | Siemens Medical Solutions Usa, Inc. | System and method for reducing circular artifacts in tomographic imaging |
JP2010252951A (ja) * | 2009-04-23 | 2010-11-11 | Hitachi Medical Corp | X線ct装置およびこれを用いたデータ取得方法 |
CN102451015A (zh) * | 2010-10-20 | 2012-05-16 | 上海西门子医疗器械有限公司 | Ct设备、以及确定探测器中不稳定通道的方法 |
CN102521801A (zh) * | 2011-11-23 | 2012-06-27 | 中国科学院深圳先进技术研究院 | Ct图像的环形伪影及弧形伪影校正方法 |
US20130101080A1 (en) * | 2011-10-19 | 2013-04-25 | Toshiba Medical Systems Corporation | Method and system for substantially reducing ring artifact based upon ring statistics |
CN108109185A (zh) * | 2017-12-18 | 2018-06-01 | 上海联影医疗科技有限公司 | 一种生成用于消除ct伪影的校正系数的方法,以及一种基于校正系数消除ct伪影的方法 |
CN109801343A (zh) * | 2018-12-16 | 2019-05-24 | 西安电子科技大学 | 基于重建前后图像的环形伪影校正方法、ct控制系统 |
CN110400358A (zh) * | 2019-07-17 | 2019-11-01 | 郑州信大先进技术研究院 | 一种平板探测器环状伪影去除方法 |
CN111053568A (zh) * | 2019-12-30 | 2020-04-24 | 苏州瑞派宁科技有限公司 | Ct图像中环形伪影的校正方法、装置及计算机存储介质 |
Family Cites Families (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4670840A (en) * | 1983-03-09 | 1987-06-02 | Elscint, Inc. | Ring artifact correction for computerized tomography |
JP3253990B2 (ja) * | 1991-10-17 | 2002-02-04 | 株式会社日立メディコ | 多次元画像計測装置 |
JPH11511375A (ja) * | 1996-03-13 | 1999-10-05 | アナロジック コーポレーション | コンピュータ断層撮影用リング抑圧フィルタ |
JP4603086B2 (ja) * | 1999-08-23 | 2010-12-22 | Geヘルスケア・ジャパン株式会社 | プロジェクションデータ補正方法、プロジェクションデータ補正装置、記録媒体および放射線断層撮像装置 |
US6498831B2 (en) * | 2000-12-22 | 2002-12-24 | Ge Medical Systems Global Technology Company, Llc | Panel detector pixel replacement method and apparatus |
DE10142605A1 (de) * | 2001-08-31 | 2003-04-03 | Siemens Ag | Verfahren zur Überprüfung der Qualität der Bilddarstellung von bildgebenden Systemen für die medizinische Diagnostik |
CN100393278C (zh) * | 2004-07-12 | 2008-06-11 | Ge医疗系统环球技术有限公司 | 环状伪影消除方法及x线ct装置 |
CN101909525B (zh) * | 2008-01-11 | 2013-06-05 | 株式会社岛津制作所 | 图像处理方法、装置以及断层摄影装置 |
JP6054036B2 (ja) * | 2012-01-26 | 2016-12-27 | 東芝メディカルシステムズ株式会社 | X線ct装置及びx線ct装置のデータ処理方法 |
CN102800073B (zh) * | 2012-06-28 | 2015-01-28 | 西北工业大学 | 一种锥束ct环形伪影的自动判别与校正方法 |
CN103593826B (zh) * | 2013-10-22 | 2016-08-17 | 南方医科大学 | 图像环形伪影校正方法 |
JP6418091B2 (ja) * | 2015-07-10 | 2018-11-07 | コニカミノルタ株式会社 | 胸部画像表示システム及び画像処理装置 |
CN106551703B (zh) * | 2015-09-30 | 2018-10-30 | 上海联影医疗科技有限公司 | 计算机断层成像方法和计算机断层成像系统 |
JP6651342B2 (ja) * | 2015-12-10 | 2020-02-19 | キヤノンメディカルシステムズ株式会社 | X線診断装置およびx線ct装置 |
CN106910165B (zh) * | 2015-12-23 | 2023-06-30 | 通用电气公司 | 修复原始ct投影数据的方法及装置、ct成像系统 |
CN105787905B (zh) * | 2016-03-24 | 2019-03-26 | 中国人民解放军信息工程大学 | 一种基于动态电流的锥束ct环状伪影校正方法 |
CN106296763B (zh) * | 2016-07-20 | 2019-05-31 | 中国兵器科学研究院宁波分院 | 一种金属材料工业ct图像质量快速校正方法 |
DE102016220347A1 (de) * | 2016-10-18 | 2018-04-19 | Siemens Healthcare Gmbh | Verfahren zum Bereitstellen von korrigierten Röntgenbildern und Röntgensystem |
JP6887030B2 (ja) * | 2017-06-26 | 2021-06-16 | エレクタ、インク.Elekta, Inc. | ディープ畳み込みニューラルネットワークを使用してコーンビームコンピュータ断層撮影画像を改善する方法 |
US10922855B2 (en) * | 2017-11-30 | 2021-02-16 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for determining at least one artifact calibration coefficient |
CN109187591B (zh) * | 2018-06-04 | 2020-10-02 | 东南大学 | 一种x射线超分辨成像方法及其应用 |
CN110264421B (zh) * | 2019-06-13 | 2022-07-12 | 明峰医疗系统股份有限公司 | 一种ct坏通道校正方法 |
-
2019
- 2019-12-30 CN CN201911392508.XA patent/CN111053568B/zh active Active
-
2020
- 2020-09-03 EP EP20910772.1A patent/EP4085843A4/en active Pending
- 2020-09-03 JP JP2022540633A patent/JP2023514944A/ja active Pending
- 2020-09-03 US US17/790,168 patent/US20230088126A1/en active Pending
- 2020-09-03 WO PCT/CN2020/113215 patent/WO2021135339A1/zh unknown
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7676073B2 (en) * | 2006-08-29 | 2010-03-09 | Siemens Medical Solutions Usa, Inc. | System and method for reducing circular artifacts in tomographic imaging |
JP2010252951A (ja) * | 2009-04-23 | 2010-11-11 | Hitachi Medical Corp | X線ct装置およびこれを用いたデータ取得方法 |
CN102451015A (zh) * | 2010-10-20 | 2012-05-16 | 上海西门子医疗器械有限公司 | Ct设备、以及确定探测器中不稳定通道的方法 |
US20130101080A1 (en) * | 2011-10-19 | 2013-04-25 | Toshiba Medical Systems Corporation | Method and system for substantially reducing ring artifact based upon ring statistics |
CN102521801A (zh) * | 2011-11-23 | 2012-06-27 | 中国科学院深圳先进技术研究院 | Ct图像的环形伪影及弧形伪影校正方法 |
CN108109185A (zh) * | 2017-12-18 | 2018-06-01 | 上海联影医疗科技有限公司 | 一种生成用于消除ct伪影的校正系数的方法,以及一种基于校正系数消除ct伪影的方法 |
CN109801343A (zh) * | 2018-12-16 | 2019-05-24 | 西安电子科技大学 | 基于重建前后图像的环形伪影校正方法、ct控制系统 |
CN110400358A (zh) * | 2019-07-17 | 2019-11-01 | 郑州信大先进技术研究院 | 一种平板探测器环状伪影去除方法 |
CN111053568A (zh) * | 2019-12-30 | 2020-04-24 | 苏州瑞派宁科技有限公司 | Ct图像中环形伪影的校正方法、装置及计算机存储介质 |
Non-Patent Citations (1)
Title |
---|
See also references of EP4085843A4 |
Also Published As
Publication number | Publication date |
---|---|
EP4085843A4 (en) | 2024-02-07 |
CN111053568B (zh) | 2021-10-08 |
JP2023514944A (ja) | 2023-04-12 |
CN111053568A (zh) | 2020-04-24 |
EP4085843A1 (en) | 2022-11-09 |
US20230088126A1 (en) | 2023-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11493458B2 (en) | Imaging system and method with scatter correction | |
US8184767B2 (en) | Imaging system and method with scatter correction | |
US8294717B2 (en) | Advanced clustering method for material separation in dual energy CT | |
JP7073380B2 (ja) | 自己較正ct検出器、自己較正を行うためのシステムおよび方法 | |
WO2021135339A1 (zh) | Ct图像中环形伪影的校正方法、装置及计算机程序介质 | |
Meng et al. | Online geometric calibration of cone-beam computed tomography for arbitrary imaging objects | |
US20100232725A1 (en) | Temperature artifact correction | |
JP2014535039A (ja) | 光子数の補正方法および装置{photoncountcorrection} | |
US8660335B2 (en) | Transient pixel defect detection and correction | |
Rashid et al. | An improved method for the removal of ring artifacts in high resolution CT imaging | |
US9773326B2 (en) | Method and device for calibrating computerized tomography (CT) image and CT system including the same | |
CN111553849B (zh) | 基于局部特征匹配的锥束ct几何伪影去除方法及装置 | |
US20100135455A1 (en) | Digital integration with detector correction | |
US20170123079A1 (en) | System and method of acquiring images using an x-ray imaging system | |
US6411671B2 (en) | Method for reducing line artifacts in a CT image and device for implementing the method | |
CN115797487A (zh) | Ct图像环状伪影自适应保结构校正方法、装置及成像设备 | |
CN109754446B (zh) | 一种探测器模块之间的拼接缝宽度估计方法及系统 | |
JP2016138757A (ja) | 放射線測定装置及び方法 | |
CN108663386B (zh) | 基于特征纹理模板的锥束ct系统探测器角度偏差测量方法 | |
Zhang et al. | A new correction method for flat panel detector in Cone-Beam CT | |
KR101938183B1 (ko) | 이동하는 대상체에 대한 방사선 영상의 화질 개선 방법 및 그 프로그램 | |
CN102735169B (zh) | 一种计算剂量曲线带宽度的系统和方法 | |
KR102193662B1 (ko) | 감마 카메라를 이용한 핫스팟 검출 장치 및 이의 핫스팟 검출 방법 | |
WO2020012520A1 (ja) | 医用x線画像処理装置およびx線画像撮影装置 | |
JP2024028042A (ja) | 情報処理装置、情報処理方法、及びプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 20910772 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2022540633 Country of ref document: JP Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
ENP | Entry into the national phase |
Ref document number: 2020910772 Country of ref document: EP Effective date: 20220801 |