CN109345608A - A kind of pyramidal CT image method for reconstructing based on asymmetric scatter removal - Google Patents
A kind of pyramidal CT image method for reconstructing based on asymmetric scatter removal Download PDFInfo
- Publication number
- CN109345608A CN109345608A CN201811230321.5A CN201811230321A CN109345608A CN 109345608 A CN109345608 A CN 109345608A CN 201811230321 A CN201811230321 A CN 201811230321A CN 109345608 A CN109345608 A CN 109345608A
- Authority
- CN
- China
- Prior art keywords
- data
- projection
- pyramidal
- asymmetric
- scatter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- 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
-
- 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/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/421—Filtered back projection [FBP]
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
The invention belongs to pyramidal CT image reconstruction technique fields, in particular to a kind of pyramidal CT image method for reconstructing based on asymmetric scatter removal, this method is first by the correcting plate occlusion area data zero setting after the Cone-Beam CT data for projection scatter correction acquired based on asymmetric scatter removal;Then original global filtering operator is replaced using part filter operator and is weighted filtering processing, effectively prevent gibbs artifact influence;It is last rationally to carry out Windowed back projection reconstruction using the redundancy consistency of data for projection, improve reconstructed image quality, and only need the data for projection of single circular scanning that can rebuild to obtain the second best in quality CT image, effectively reduce scanning times, scan efficiency is greatly improved, dose of radiation is reduced.
Description
Technical field
The invention belongs to pyramidal CT image reconstruction technique field, in particular to a kind of cone based on asymmetric scatter removal
Beam CT image rebuilding method.
Background technique
In cone-beam CT imaging system, due to being influenced by factors such as ray cone angle, the irradiation visual field, signal cross-talks, plate
Often contain more scattered ray in the received incident ray of detector, so that the CT value of reconstruction image deviates ideal true value,
Cup-shaped, bar shaped or banding artifact, shade are formed, picture contrast is greatly reduced, seriously affects the quality of reconstruction image.Cause
This, it usually needs scatter artefacts correction is carried out before three-dimensional reconstruction to the CT data of acquisition, to promote reconstructed image quality, is improved
Reconstruction image identification.
Passed through based on the direct measurement bearing calibration of asymmetric scatter removal (structure is as shown in Figure 2) larger using length
Highly attenuating effect of the leads to X-ray, can directly measure the high-precision scatter distributions of projection, avoid and consider X-ray light
The complex process of son and object interaction can obtain preferable calibration result compared to other algorithms to complicated object.
But the primary signal as caused by the blocking of lead is lost, and directly often has gibbs artifact using FDK algorithm progress three-dimensional reconstruction
It generates, seriously affects reconstructed image quality.
Summary of the invention
Aiming at the problems existing in the prior art, the present invention provides a kind of Cone-Beam CT figure based on asymmetric scatter removal
As method for reconstructing, redundancy of this method based on part filter operator and data for projection, it is only necessary to which carrying out single sweep operation can press down
Gibbs artifact processed influences to complete to rebuild the scatter correction of data for projection, obtains the preferable three-dimensional reconstruction image of quality.
To achieve the goals above, the following technical solution is employed by the present invention:
A kind of pyramidal CT image method for reconstructing based on asymmetric scatter removal provided by the invention includes following step
It is rapid:
Step 1, zero-padding pre-processing is carried out to the data for projection after scatter correction;
Step 2, the pretreated data for projection of zero setting is weighted and is corrected, be fan-beam number by Cone beam projection data Approximate revision
According to;
Step 3, to revised data for projection part filter line by line;
Step 4, the redundancy otherness of data for projection is eliminated to filtered backprojection data weighting using window function;
Step 5, backprojection reconstruction obtains cone-beam CT three-dimensional image.
Further, before the step 1 further include:
Asymmetric scatter removal is placed between object and detector, and is acquired using cone-beam CT system single sweep operation
Data for projection;
Correction is scattered to the data for projection of acquisition.
Further, the realization process that the described pair of data for projection acquired is scattered correction is as follows:
Step 201, effective measured value is extracted
Taking the data in projection in the one third width of leads occlusion area center is effective measured value, and by each leads
The longitudinal direction of effective coverage takes mean value as the measured value at each longitudinal center position of the leads;
Step 202, smooth noise reduction
Mean filter is carried out to effective measured value of extraction;
Step 203, cubic spline interpolation
Based on the discrete measured values after smoothing processing, longitudinal data is carried out first using cubic spline function slotting
Then value carries out interpolation to lateral data, obtains the whole scatter distributions projected under the angle;
Step 204, scatter distributions are deducted
The whole scatter distributions value that interpolation obtains directly is deducted from data for projection;
Step 205, BM3D flat volatility
It is smoothly fluctuated as directly deducting data for projection caused by scatter distributions using BM3D algorithm.
Further, zero-padding pre-processing is data to asymmetric scatter removal leads occlusion area in the step 1
Zero setting defines shown in zero setting template such as formula (1):
Wherein, D is the regional ensemble that blocks of leads in projection, and (i, j) be the cell coordinate in projection, by zero setting template and
Data for projection convolution is multiplied.
Further, the step 2 weights modified process such as formula (2) institute to the pretreated data for projection of zero setting
Show:
Wherein, d is distance of the light source to rotation center, and (u (j), v (i)) is detector plane coordinate,Indicate the
Data after the projection zero-padding pre-processing of θ angle, P 'θ(i, j) indicates to weight modified result to data for projection.
Further, part filter is filtered line by line to revised whole picture data for projection in the step 3, is gone forward side by side
Row truncated error inhibition processing, shown in filtering such as formula (3):
Wherein, P "θ(i, j) is the projection after part filter;It is handled by truncated error inhibition, by filtered backprojection data
Value zero setting in leads occlusion area.
Further, window function is designed based on asymmetric scatter removal in the step 4, defines window function such as formula
(4) shown in:
Further, the redundancy that the step 4 eliminates data for projection to filtered backprojection data weighting using window function is poor
Shown in anisotropic process such as formula (5):
Wherein,For the projection for eliminating redundancy otherness.
Further, step 5 backprojection reconstruction obtains shown in the process such as formula (6) of cone-beam CT three-dimensional image:
Wherein, f (x, y, z) indicates to be reconstructed point, and u (x, y, θ), v (x, y, z, θ) respectively indicate point (x, y, z) and project to
Coordinate on detector.
Compared with prior art, the invention has the following advantages that
The invention discloses a kind of pyramidal CT image method for reconstructing based on asymmetric scatter removal, this method first will
Correcting plate occlusion area data zero setting after Cone-Beam CT data for projection scatter correction based on the acquisition of asymmetric scatter removal;So
Original global filtering operator is replaced using part filter operator afterwards and is weighted filtering processing, effectively prevents gibbs artifact shadow
It rings;It is last rationally to carry out Windowed back projection reconstruction using the redundancy consistency of data for projection, reconstructed image quality is improved, and
It only needs the data for projection of single circular scanning that can rebuild to obtain the second best in quality CT image, effectively reduces scanning times, pole
Scan efficiency is improved greatly, reduces dose of radiation.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
It obtains other drawings based on these drawings.
Fig. 1 is a kind of flow chart of pyramidal CT image method for reconstructing based on asymmetric scatter removal of the invention;
Fig. 2 is the structural schematic diagram of asymmetric scatter removal;
Fig. 3 is that collective modes projection, the two dimension slicing figure rebuild are adopted using the method for the present invention;
Fig. 4 is that asymmetric scatter removal full scan is not added to adopt collective modes projection, two with Fig. 3 identical layer rebuild
Tie up slice map;
Fig. 5 is the hatching line figure of Fig. 3 transversal centerline;
Fig. 6 is the hatching line figure of Fig. 4 transversal centerline.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
Core of the invention is to provide a kind of pyramidal CT image method for reconstructing based on asymmetric scatter removal, effectively keeps away
Exempt from gibbs artifact influence, obtains the preferable three-dimensional reconstruction image of quality.
Present invention will be described in further detail below with reference to the accompanying drawings and specific embodiments.
Embodiment one
As shown in Figure 1, a kind of pyramidal CT image method for reconstructing based on asymmetric scatter removal of the present embodiment, includes
Following steps:
Asymmetric scatter removal is placed between object and detector by step S101, and utilizes cone-beam CT system list
Secondary circular scanning acquired projections data.
Step S102, the data for projection acquired to step 1 are scattered correction, the specific implementation process is as follows:
Step S1021 extracts effective measured value
To reduce measurement error, taking the data in projection in the w/3 width of leads occlusion area center is effective measured value, and
Take mean value as the measured value at each longitudinal center position of the leads longitudinal direction of each leads effective coverage;
Step S1022, smooth noise reduction
Mean filter is carried out to reduce the noise jamming in measured value to effective measured value of extraction;
Step S1023, cubic spline interpolation
Based on the discrete measured values after smoothing processing, longitudinal data is carried out first using cubic spline function slotting
Then value carries out interpolation to lateral data, obtains the whole scatter distributions projected under the angle;
Step S1024 deducts scatter distributions
The whole scatter distributions value that interpolation obtains directly is deducted from data for projection;
Step S1025, BM3D flat volatility
After directly deducting scatter distributions, the fluctuation of data for projection is larger, influences subsequent reconstruction as a result, therefore further sharp
With BM3D algorithm come the fluctuation of smooth data for projection.
Step S103 carries out zero-padding pre-processing to the data for projection after scatter correction;
Due to the shielding action of correcting plate leads, X-ray cannot all obtain just in the leads blocked area of every data for projection
Partial data participation reconstruction process, can be obtained the reconstructed results of mistake, needed the part by the data of true irradiation object
Data zero setting rebuilds to obtain correct result merely with the projection information of its redundancy angle.Therefore right first before filtered back projection
Data for projection after scatter correction carries out zero-padding pre-processing, defines shown in zero setting template such as formula (1):
Wherein, D is the regional ensemble that blocks of leads in projection, and (i, j) be the cell coordinate in projection, by zero setting template and
Data for projection convolution is multiplied.
Step S104 weights the pretreated data for projection of zero setting and corrects, and is fan by Cone beam projection data Approximate revision
Beam data, shown in weighting makeover process such as formula (2):
Wherein, d is distance of the light source to rotation center, (u (j), v (i)))) it is detector plane coordinate,Table
Show data after the projection zero-padding pre-processing of θ angle, Pθ' (i, j) indicates to weight modified result to data for projection.
Step S105, to revised data for projection, part filter, inhibition block caused transversal cross-section by leads line by line
Artifact, shown in filtering such as formula (3):
Wherein, P "θ(i, j) is the projection after part filter;It is handled by truncated error inhibition, by filtered backprojection data
Value zero setting in leads occlusion area D reduces transmitting of the filtering error to reconstruction image.
Step S106 eliminates the redundancy otherness of data for projection using window function to filtered backprojection data weighting.
It is influenced by leads shielding, the redundancy of each point is inconsistent in data for projection, and directly carrying out back projection can be due to
Redundancy difference generates artifact and needs to be weighted really the data of different location before back projection to eliminate this artifact
Protect the anti-consistency for throwing redundancy of data for projection;It defines shown in window function such as formula (4):
Then data for projection is weighted and is eliminated shown in the process such as formula (5) of redundancy otherness:
Wherein,For the projection for eliminating redundancy otherness.
Step S107, backprojection reconstruction obtains cone-beam CT three-dimensional image, shown in backprojection reconstruction process such as formula (6):
Wherein, f (x, y, z) indicates to be reconstructed point, and u (x, y, θ), v (x, y, z, θ) respectively indicate point (x, y, z) and project to
Coordinate on detector.
Effect to further illustrate the method for the present invention is handled specific data for projection below using the method for the present invention.
It adopts collective modes using the method for the present invention to project, the two dimension slicing rebuild is as shown in figure 3, be not added asymmetric scattering
Correcting plate full scan adopt collective modes projection, it is rebuilding with Fig. 3 identical layer two dimension slicing as shown in figure 4, Fig. 5 be Fig. 3 cross
Hatching line figure to the midline, Fig. 6 are the hatching line figures of Fig. 4 transversal centerline.From slice map 3 is rebuild, it is apparent that the method for the present invention
The influence of truncated data is effectively inhibited, reconstruction has obtained the second best in quality CT image;Comparison diagram 3 and Fig. 4 are can be found that simultaneously
The contrast of the method for the present invention reconstruction image is obviously improved compared to the contrast of full scan reconstruction image, and the method for the present invention is obvious
Improve the lower part of contrast (position of black arrow meaning in such as Fig. 3) and edge mould in image caused by scattering shade
Paste problem distinguish the object of two kinds of materials of ectonexine can intuitively.In addition from hatching line figure comparison it can be found that
The CT value of the method for the present invention remains basically stable in material of the same race, and cupping artifact has also obtained effective inhibition.
Further to evaluate the method for the present invention reconstructed results, the quantized data of contrast and cupping artifact is calculated
Quantitative analysis is carried out, specific quantized data is as shown in table 1, wherein contrast=| μ1-μ2|, μ1It is in interior of articles rectangle frame
Mean CT-number, μ2It is the mean CT-number near interior of articles rectangle frame in rectangle frame.μaIt is on the outside of object
The mean CT-number of border circular areas, μbIt is the mean CT-number of object central circular, rectangle frame and border circular areas such as Fig. 3 and Fig. 4
Shown in.It can clearly be seen that the method for the present invention effectively inhibits the cupping artifact in reconstruction image from quantization table, improve
Picture contrast avoids the influence of gibbs artifact.
The quantitative analysis table of the correction of table 1 front and back body mould
A kind of pyramidal CT image method for reconstructing based on asymmetric scatter removal of the invention, to the data for projection of acquisition
Windowed back projection reconstruction is carried out using part filter operator, gibbs artifact influence is effectively prevented, takes full advantage of data for projection
Redundancy, improve reconstructed image quality, it is only necessary to single pass can realize the good reconstruction of scatter correction CT image,
Scan efficiency is improved, dose of radiation is reduced.
Finally, it should be noted that embodiment described above, only a specific embodiment of the invention, to illustrate the present invention
Technical solution, rather than its limitations, scope of protection of the present invention is not limited thereto, although with reference to the foregoing embodiments to this hair
It is bright to be described in detail, those skilled in the art should understand that: anyone skilled in the art
In the technical scope disclosed by the present invention, it can still modify to technical solution documented by previous embodiment or can be light
It is readily conceivable that variation or equivalent replacement of some of the technical features;And these modifications, variation or replacement, do not make
The essence of corresponding technical solution is detached from the spirit and scope of technical solution of the embodiment of the present invention, should all cover in protection of the invention
Within the scope of.Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (9)
1. a kind of pyramidal CT image method for reconstructing based on asymmetric scatter removal, which is characterized in that comprise the steps of:
Step 1, zero-padding pre-processing is carried out to the data for projection after scatter correction;
Step 2, the pretreated data for projection of zero setting is weighted and is corrected, be fan beam data by Cone beam projection data Approximate revision;
Step 3, to revised data for projection part filter line by line;
Step 4, the redundancy otherness of data for projection is eliminated to filtered backprojection data weighting using window function;
Step 5, backprojection reconstruction obtains cone-beam CT three-dimensional image.
2. the pyramidal CT image method for reconstructing according to claim 1 based on asymmetric scatter removal, which is characterized in that
Before the step 1 further include:
Asymmetric scatter removal is placed between object and detector, and utilizes cone-beam CT system single sweep operation acquired projections
Data;
Correction is scattered to the data for projection of acquisition.
3. the pyramidal CT image method for reconstructing according to claim 2 based on asymmetric scatter removal, which is characterized in that
The realization process that the data for projection of described pair of acquisition is scattered correction is as follows:
Step 201, effective measured value is extracted
Taking the data in projection in the one third width of leads occlusion area center is effective measured value, and each leads is effective
The longitudinal direction in region takes mean value as the measured value at each longitudinal center position of the leads;
Step 202, smooth noise reduction
Mean filter is carried out to effective measured value of extraction;
Step 203, cubic spline interpolation
Based on the discrete measured values after smoothing processing, interpolation is carried out to longitudinal data first using cubic spline function, so
Interpolation is carried out to lateral data afterwards, obtains the whole scatter distributions projected under the angle;
Step 204, scatter distributions are deducted
The whole scatter distributions value that interpolation obtains directly is deducted from data for projection;
Step 205, BM3D flat volatility
It is smoothly fluctuated as directly deducting data for projection caused by scatter distributions using BM3D algorithm.
4. the pyramidal CT image method for reconstructing according to claim 3 based on asymmetric scatter removal, which is characterized in that
Zero-padding pre-processing is the data zero setting to asymmetric scatter removal leads occlusion area in the step 1, defines zero setting template
As shown in formula (1):
Wherein, D is the regional ensemble that leads blocks in projection, and (i, j) is the cell coordinate in projection, by zero setting template and projection
Data convolution is multiplied.
5. the pyramidal CT image method for reconstructing according to claim 4 based on asymmetric scatter removal, which is characterized in that
The step 2 weights shown in modified process such as formula (2) the pretreated data for projection of zero setting:
Wherein, d is distance of the light source to rotation center, and (u (j), v (i)) is detector plane coordinate,Indicate the angle θ
Data after the projection zero-padding pre-processing of degree, Pθ' (i, j) indicates to weight modified result to data for projection.
6. the pyramidal CT image method for reconstructing according to claim 5 based on asymmetric scatter removal, which is characterized in that
Part filter is to be filtered line by line to revised whole picture data for projection, and carry out at truncated error inhibition in the step 3
Reason, shown in filtering such as formula (3):
Wherein, P "θ(i, j) is the projection after part filter;It is handled by truncated error inhibition, by filtered backprojection data in lead
Value zero setting in occlusion area.
7. the pyramidal CT image method for reconstructing according to claim 6 based on asymmetric scatter removal, which is characterized in that
Window function is designed based on asymmetric scatter removal in the step 4, is defined shown in window function such as formula (4):
8. the pyramidal CT image method for reconstructing according to claim 7 based on asymmetric scatter removal, which is characterized in that
The step 4 eliminates the process such as formula of the redundancy otherness of data for projection using window function to filtered backprojection data weighting
(5) shown in:
Wherein,For the projection for eliminating redundancy otherness.
9. the pyramidal CT image method for reconstructing according to claim 8 based on asymmetric scatter removal, which is characterized in that
Step 5 backprojection reconstruction obtains shown in the process such as formula (6) of cone-beam CT three-dimensional image:
Wherein, f (x, y, z) indicates to be reconstructed point, and u (x, y, θ), v (x, y, z, θ) respectively indicate point (x, y, z) and project to detection
Coordinate on device.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811230321.5A CN109345608B (en) | 2018-10-22 | 2018-10-22 | Asymmetric scattering correction plate-based cone beam CT image reconstruction method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811230321.5A CN109345608B (en) | 2018-10-22 | 2018-10-22 | Asymmetric scattering correction plate-based cone beam CT image reconstruction method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109345608A true CN109345608A (en) | 2019-02-15 |
CN109345608B CN109345608B (en) | 2022-10-28 |
Family
ID=65311410
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811230321.5A Active CN109345608B (en) | 2018-10-22 | 2018-10-22 | Asymmetric scattering correction plate-based cone beam CT image reconstruction method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109345608B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109949233A (en) * | 2019-02-18 | 2019-06-28 | 深圳蓝韵医学影像有限公司 | Filter out the method, system, equipment and storage medium of scattered rays in radioscopic image |
CN111553959A (en) * | 2020-03-27 | 2020-08-18 | 中国人民解放军战略支援部队信息工程大学 | Cone beam CT truncation artifact suppression method based on projection hyperbolic extrapolation |
CN112053409A (en) * | 2020-07-24 | 2020-12-08 | 重庆真测科技股份有限公司 | Asymmetric data reconstruction method based on double-rotating-table CT scanning system |
CN112446931A (en) * | 2019-09-02 | 2021-03-05 | 上海联影医疗科技股份有限公司 | Reconstruction data processing method and device, medical imaging system and storage medium |
CN116531015A (en) * | 2023-07-04 | 2023-08-04 | 中国科学院深圳先进技术研究院 | Image reconstruction method, device, electronic equipment and storage medium |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101158653A (en) * | 2007-11-16 | 2008-04-09 | 西北工业大学 | Diffuse transmission measuring and correcting method of cone-beam CT system |
US20120207370A1 (en) * | 2010-12-20 | 2012-08-16 | Benjamin Pooya Fahimian | Systems and Methods for Simultaneous Acquisition of Scatter and Image Projection Data in Computed Tomography |
CN103578082A (en) * | 2012-08-09 | 2014-02-12 | 江苏超惟科技发展有限公司 | Cone beam CT scatter correction method and system |
-
2018
- 2018-10-22 CN CN201811230321.5A patent/CN109345608B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101158653A (en) * | 2007-11-16 | 2008-04-09 | 西北工业大学 | Diffuse transmission measuring and correcting method of cone-beam CT system |
US20120207370A1 (en) * | 2010-12-20 | 2012-08-16 | Benjamin Pooya Fahimian | Systems and Methods for Simultaneous Acquisition of Scatter and Image Projection Data in Computed Tomography |
CN103578082A (en) * | 2012-08-09 | 2014-02-12 | 江苏超惟科技发展有限公司 | Cone beam CT scatter correction method and system |
Non-Patent Citations (2)
Title |
---|
谢世朋等: "基于自适应点扩散函数的锥束CT散射校正", 《中国医学影像技术》 * |
郝佳等: "探测器偏置CBCT系统加权校正重建方法", 《清华大学学报(自然科学版)》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109949233A (en) * | 2019-02-18 | 2019-06-28 | 深圳蓝韵医学影像有限公司 | Filter out the method, system, equipment and storage medium of scattered rays in radioscopic image |
CN112446931A (en) * | 2019-09-02 | 2021-03-05 | 上海联影医疗科技股份有限公司 | Reconstruction data processing method and device, medical imaging system and storage medium |
CN111553959A (en) * | 2020-03-27 | 2020-08-18 | 中国人民解放军战略支援部队信息工程大学 | Cone beam CT truncation artifact suppression method based on projection hyperbolic extrapolation |
CN111553959B (en) * | 2020-03-27 | 2023-03-03 | 中国人民解放军战略支援部队信息工程大学 | Cone beam CT truncation artifact suppression method based on projection hyperbolic extrapolation |
CN112053409A (en) * | 2020-07-24 | 2020-12-08 | 重庆真测科技股份有限公司 | Asymmetric data reconstruction method based on double-rotating-table CT scanning system |
CN112053409B (en) * | 2020-07-24 | 2024-05-28 | 重庆真测科技股份有限公司 | Asymmetric data reconstruction method based on double-rotating-table CT scanning system |
CN116531015A (en) * | 2023-07-04 | 2023-08-04 | 中国科学院深圳先进技术研究院 | Image reconstruction method, device, electronic equipment and storage medium |
CN116531015B (en) * | 2023-07-04 | 2023-10-03 | 中国科学院深圳先进技术研究院 | Image reconstruction method, device, electronic equipment and storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN109345608B (en) | 2022-10-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109345608A (en) | A kind of pyramidal CT image method for reconstructing based on asymmetric scatter removal | |
US8989469B2 (en) | Systems and methods for simultaneous acquisition of scatter and image projection data in computed tomography | |
US7382853B2 (en) | Method and system of CT data correction | |
US7283605B2 (en) | Methods and apparatus for scatter correction | |
US8280135B2 (en) | System and method for highly attenuating material artifact reduction in x-ray computed tomography | |
US10013780B2 (en) | Systems and methods for artifact removal for computed tomography imaging | |
US10395397B2 (en) | Metal artifacts reduction for cone beam CT | |
US8553959B2 (en) | Method and apparatus for correcting multi-modality imaging data | |
JP2018110866A (en) | Medical image generation device and medical image generation method | |
US20140328450A1 (en) | System and method for reducing high density artifacts in computed tomography imaging | |
US9483851B2 (en) | Systems and methods for filtering for image generation | |
US10595808B2 (en) | System and method of computed tomography signal restoration via noise reduction | |
KR101255224B1 (en) | X-ray computed tomography system and scatter correction method using the same | |
US9001961B2 (en) | Methods of scatter correction of X-ray projection data 1 | |
JP5214110B2 (en) | X-ray CT system | |
EP1433137B1 (en) | A method of reducing artifacts in object images | |
US9427205B1 (en) | Systems and methods for artifact removal for computed tomography imaging | |
US9901314B2 (en) | Adjustable bow-tie filter for achieving optimal SNR in helical computed tomography | |
US20120177173A1 (en) | Method and apparatus for reducing imaging artifacts | |
KR101914522B1 (en) | Cone Beam Computed Tomography System for Correcting Scatter Using Binary Moving Blocker | |
TWI613998B (en) | Reduction method for boundary artifact on the tomosynthesis | |
Niu et al. | Single-scan scatter correction for cone-beam CT using a stationary beam blocker: a preliminary study | |
Zhang et al. | An improved scatter correction method for cone beam computed tomography | |
Mizuta et al. | Implementation of histogram based soft-tissue segmentation for single spiral transmission scanning in whole body PET | |
Tang et al. | X-ray flat-panel imager (FPI)-based cone-beam volume CT (CBVCT) under a circle-plus-two-arc data acquisition orbit |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |