CN105528771B - The bearing calibration of cupping artifact in a kind of Cone-Beam CT using energy function method - Google Patents
The bearing calibration of cupping artifact in a kind of Cone-Beam CT using energy function method Download PDFInfo
- Publication number
- CN105528771B CN105528771B CN201610035576.0A CN201610035576A CN105528771B CN 105528771 B CN105528771 B CN 105528771B CN 201610035576 A CN201610035576 A CN 201610035576A CN 105528771 B CN105528771 B CN 105528771B
- Authority
- CN
- China
- Prior art keywords
- image
- cone
- artifact
- correction
- reconstruction image
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000012937 correction Methods 0.000 claims abstract description 49
- 230000008859 change Effects 0.000 claims abstract description 6
- 239000000126 substance Substances 0.000 abstract description 4
- 238000012800 visualization Methods 0.000 abstract description 3
- 238000011161 development Methods 0.000 abstract description 2
- 238000007794 visualization technique Methods 0.000 abstract description 2
- 238000002591 computed tomography Methods 0.000 description 23
- 210000000988 bone and bone Anatomy 0.000 description 6
- 238000000342 Monte Carlo simulation Methods 0.000 description 5
- 210000003625 skull Anatomy 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000007408 cone-beam computed tomography Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000004445 quantitative analysis Methods 0.000 description 3
- 239000004576 sand Substances 0.000 description 3
- 238000000333 X-ray scattering Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003702 image correction Methods 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000010603 microCT Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30008—Bone
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
It is used through the corrected method of cupping artifact in optimized energy function pair Cone-Beam CT the invention discloses a kind of, this method is for the sectioning image after rebuilding, it is directly facing user, it does not make any change to original cone-beam CT equipment, it can complete correction work, the CT value uniformities of the substance of the same race of reconstruction image can also be improved while the cupping artifact correction that can be carried out efficiently Cone-Beam CT, to contribute in reconstruction image, the development of perfect volume visualization and the visualization technique based on threshold value.This method is applied to cone-beam CT reconstruction image (i.e. image area) alignment technique field.
Description
Technical field
The present invention relates to technical field of medical image processing more particularly to the correction of pyramidal CT image cupping artifact, gray scales not
Uniformity correction technical field of medical image processing.
Background technology
Cone-Beam CT is usually used in image guided therapy, upper abdomen inspection as the medical treatment and industrial detection instrument developed in recent years
It looks into, examination of mouth, industrial detection etc..Cone-Beam CT (CBCT) based on flat panel detector is compared with traditional two-dimensional ct, tool
Have the advantages that prominent, be mainly manifested in circular scanning period of Cone-Beam CT, can obtain completing hundreds of or even thousands of it is disconnected
The projection of tomographic image has higher sweep speed and radiation utilization rate, and effectively reduces the load output of X-ray tube, drop
Low scanning cost.The many because being known as of cone-beam CT reconstruction picture quality are influenced, such as x-ray scattering, noise, geometric error, power spectrum, spy
Survey unit non_uniform response etc..But since cone-beam tablet CT uses large-scale X-ray flat panel detector, this makes image quality
Be more vulnerable to compared with traditional CT X-ray scattering and beam hardening influence.Because of the artifact for scattering with being formed due to beam hardening
(main includes cupping artifact and streak artifact) seriously affects analysis and judgement to reconstruction image.In the cone-beam of medical grade
In CT reconstruction images, cupping artifact occupies very big proportion, in terms of these artifacts are for the visualization display based on threshold value and is based on
It is influenced in terms of the pyramidal CT image segmentation of threshold value very serious.And the correction of cupping artifact can provide instead for other artifact corrections
Feedback reference, being other, cone beam computed tomography (CT) scattering and beam hardening correction without priori provide verification information.Therefore the present invention is about cone-beam
The correction of cupping artifact in CT is very significant.
In order to reduce cupping artifact (i.e.:CT value inhomogeneities artifact) influence, currently, the prior art or documents and materials
Research is concentrated mainly on to the scatter correction in projected image.The artifact correction of early stage is mainly reflected in hardware based correction,
Such as X-ray filter line device, collimator or metal grate, air-gap method, scanning slit technology and leads or stereotype technology.Recently
Artifact correction research in several years is mainly reflected in based on monte carlo method, scattering analysis method of estimation and is based on partial dispersion ray
The scatter correction method of measurement.Monte Carlo simulation is effectively method in CBCT scatter corrections, but calculation amount is huge.
In recent years, some improved Monte Carlo simulation algorithms are also suggested, and such as use GPU acceleration techniques, the body based on model extensive
Compound method etc..These thoughts based on Monte Carlo simulation all attempt to establish one preferably on simulation precision and calculating cost
Equalization point.But it is too high to be confined to calculating cost always, without easy to use.
There are BSA scatter removal methods in method based on the Cone-Beam CT artifact correction that part ray blocks, early stage, are logical
It crosses measurement and blocks the scattered quantum below array in ray, carry out interpolation and go out the whole scatter distributions reached on detector.Then again into
Row blocks the normal scan of array without ray, and scatter distributions image is subtracted from the projected image of scanning, is just corrected
Image afterwards.This method will carry out twice sweep, and x-ray irradiation dose is also increased while increasing sweep time.Afterwards
Fast method is blocked to occur moving, can solve the problems, such as twice sweep.
Also there is the research to the image progress artifact correction after cone-beam CT reconstruction.The research is mainly by the solution in CT images
Cut open structure.Such methods are completely dependent on deformable registration precision, and need CT image datas.
The shortcomings that prior art includes mainly:
(1) currently available technology is primarily directed to projection image correction, not directly against the school of sectioning image after reconstruction
Just, such as patent CN104408753A.
(2) most of prior art concentrates in the method for artifact correction caused by due to scattering, and needs to add mostly
Add hardware device, such as:Patent 200710019084 and 201310039298, the two patents are required for setting in expensive Cone-Beam CT
Standby upper addition hardware device, increases the complexity of operation and causes potential security risk to equipment.Especially patent
200710019084 need twice sweep testee, so the undoubtedly amount of radiation of increased measured object.
In conclusion in the method for the prior art or documents and materials, Monte-carlo Simulation Method expends the time very much, just
Correction result is limited to the structure of modulation panel itself in grade ray modulation method, is based on partial dispersion radionetric survey method, needs
Increase exposure dose, existing method is not high to the accuracy of estimation of scatter distributions.And the present invention can well solve above
Problem.
Invention content
Present invention aims at solve above-mentioned the deficiencies in the prior art, it is proposed that a kind of use passes through optimized energy letter
It is several to the corrected method of cupping artifact in Cone-Beam CT.This method is directly facing user, no for the sectioning image after rebuilding
It makes any change to original cone-beam CT equipment, so that it may which, to complete correction work, the cup-shaped that can be carried out efficiently Cone-Beam CT is pseudo-
Shadow can also improve the CT value uniformities of the substance of the same race of reconstruction image while correction, complete to contribute in reconstruction image
The development of kind volume visualization and the visualization technique based on threshold value.This method is applied to cone-beam CT reconstruction image (i.e. image area)
Alignment technique field.
The technical scheme adopted by the invention to solve the technical problem is that:Weight in Cone-Beam CT is used the present invention provides a kind of
The optimized energy functional based method that image carries out cupping artifact correction is built, this method has very strong robustness, need not repeat
Testee is scanned, does not increase the complexity of cone-beam CT system.
Method flow:
Step 1:Obtain the sectioning image after rebuilding.
Step 2:Cupping artifact indicates and builds energy function
Wherein G (x)=(g1(x),…gM(x))TFor smooth basis functions, ciFor constant, meets and work as x ∈ ΩiWhen, ui(x)=
1;WhenWhen, ui(x)=0, w=(w1,…wM)
Step 4:Newer w and u is fixed and uses, withFor F (u, c, w) minimum value of variable
Xie Wei:
Step 5:Newer w and c is fixed and uses, with u=(u1,…uN)TFor variable F (u, c, w) minimum value solution when,Meet following condition:
Wherein, imin(x)=argmin { f (x)-ci-wTG(x)}。
Step 6:If w is stable or iterations are more than 10 times, 7 are thened follow the steps, step 3 is otherwise returned to.
Further, the present invention is directly facing Cone-Beam CT slice of data, is not carried out to original Cone-Beam CT existing equipment any
Change, does not need the prior information of user and measured target.
Further, present invention demonstrates that cupping artifact is to decomposite to come from reconstruction image, i.e.,:
Wherein fpIndicate true artifact-free sectioning image, fsIndicate that scattering and beam hardening cause the sectioning image of artifact,
fnIndicate the sectioning image that noise n is formed.
Advantageous effect:
1, the present invention is the cupping artifact correction directly against the sectioning image after reconstruction, and this method calculation amount is relatively
It is small, while capable of efficiently carrying out the correction of Cone-Beam CT sectioning image cupping artifact, improve the CT of substance reconstruction image of the same race
It is worth uniformity.
2, the present invention is directly facing CT slice demand users, does not make any change to original Cone-Beam CT existing equipment, no
The prior information for needing user and measured target, completes correction work well.
3, the present invention can increase picture contrast, and image more accurately shows illuminated object letter originally after enabling correction
Breath.
4, the present invention improves and realizes CT image viewings, segmentation and lesion detection based on threshold value well.
Description of the drawings
Fig. 1 is flow chart of the method for the present invention.
Fig. 2 is human skull's sample axial view on (left side) (right side) afterwards before the present invention is corrected by reconstruction image.
Fig. 3 is the horizontal profile value for human skull's die body that the present invention measures:The longitudinal axis is image profile.
Fig. 4 is the region of interest area image schematic diagram of the present inventor's skull die body selection.
Fig. 5 (a), Fig. 5 (b) cut for two differences of CTP486 reconstruction images in CatPhan500 before and after scatter correction of the present invention
The sample axial view of piece.
Fig. 6 is the horizontal sectional drawing of die body measured by Fig. 5 in the present invention.
Fig. 7 is the selection zoning image schematic diagram of Mouse Bone.
Fig. 8 (a), Fig. 8 (b) are respectively the sample axial view of Mouse Bone reconstruction image before and after scatter correction.
Specific implementation mode
The invention is described in further detail with reference to the accompanying drawings of the specification.
As shown in Figure 1, the present invention provides a kind of correction sides of cupping artifact in Cone-Beam CT using energy function method
Method, this method comprises the following steps:
Step 1:According to FDK algorithms obtain rebuild after sectioning image, and can from theoretical proof cupping artifact can from weight
Rear picture breakdown is built to come out.
Step 2:Cupping artifact indicates and builds energy function
Wherein G (x)=(g1(x),…gM(x))TFor smooth basis functions, ciFor constant, meets and work as x ∈ ΩiWhen, ui(x)=
1;WhenWhen, ui(x)=0, w=(w1,…wM)
Step 4:Newer w and u is fixed and uses, withFor F (u, c, w) minimum value of variable
Xie Wei:
Step 5:Newer w and c is fixed and uses, with u=(u1,…uN)TFor variable F (u, c, w) minimum value solution when,Meet following condition:
Wherein, imin(x)=argmin { f (x)-ci-wTG(x)}。
Step 6:If w is stable or iterations are more than 10 times, 7 are thened follow the steps, step 3 is otherwise returned to.
Cupping artifact of the present invention can be decomposited from reconstruction image to be carried out detailed process and includes:
For reconstruction image based on FDK algorithms, reconstruction image collection f can be written as following form in Cone-Beam CT:
Wherein dsoIndicate radiographic source to the distance of rotary shaft, I3D(t, z (r), φ) indicates the sequence of projected image.Projection
Image I3DIt can be analyzed to following form:
I3D=P3D+S3D+ n, formula 2
Wherein P3DRepresentation theory real projection image, S3DIt is the artifact ingredient as caused by scattering and beam hardening, n is equal
The additive noise that value is zero.Formula 1 is expressed as form:
Wherein fpIndicate true artifact-free sectioning image, fsIndicate that scattering and beam hardening cause the sectioning image of artifact,
fnIndicate the sectioning image that noise n is formed.
By formula 3 it is found that reconstruction image can be expressed as three independent elements additions.Wherein fsIt is by S3DIt obtains, and S3D
It is a smooth low frequency projection signal, all fsMain region (shows as the Physical Zone of the same race of large area, cupping artifact
It is its main component) it should also be a smooth low frequency sectioning image.In order to effectively use fsAnd fpProperty, the present invention will
fsIt is expressed as a known smooth basis functions collection g1,…gMLinear combination, this adapt to cupping artifact smooth change property.
By finding linear combinationIn optimum coefficient w=(w1,…wM) come estimate cupping artifact (present invention
M=10 is taken in the experimental data enumerated).The present invention is fs(x) it is expressed as fs(x)=wTThe vector form of G (x), wherein G (x)=
(g1(x),…gM(x))T。
Assuming that in image area ΩiIn there are N kind tissues, then true sectioning image fp(x) it is in actually being organized at i-th
A constant c about xi.Each ΩiIt can be with its membership function come uiIt indicates.Under ideal conditions, uiIt is one two
System membership function meets and works as x ∈ ΩiWhen, ui(x)=1;WhenWhen, ui(x)=0.As membership function uiWith constant ci
When known, fpIt may be expressed as following form:
Energy functional expression formula includes:
In this model, the present invention considers how in reconstruction image f, obtains when following energy function F being made to minimize
fsAnd fp:
If obvious variable fsAnd fpThere is no any constraints, the minimum of F is an ill-conditioning problem.In fact, working as fs
And fpTo meet condition fp=f-fsArbitrary value when, energy F (fs,fp) it can all obtain minimum value.Use true picture and cup-shaped
The expression formula of artifacts, energy F (fs,fp) following form can be represented as:
Optimization includes:
All variable u=(u of energy F (u, c, w)1,…uN), c=(c1,…cN) and w be all convex, this property is true
Protected F (u, c, w) has unique optimal solution for any variable.Pass through F (u, c, w) minimum value under interleaved computation difference variable
Solution, you can reach solution purpose.
First, fixed c and u, the present invention can be by solving equationTo obtain the minimum value of F (u, c, w).It crosses
Journey is as follows:
Above-mentioned equation can be rewritten as following form:
Aw=v, formula 8
It is easy to prove, matrix A is nonsingular matrix.Therefore, vectorialIt can be represented as:
Newer fs, can be calculated and be obtained by following formula:
Newer w and u is fixed and uses, withIt is for F (u, c, w) minimum value solution of variable:
Newer w and c is fixed and uses, with u=(u1,…uN)TFor variable F (u, c, w) minimum value solution when,Meet following condition:
Wherein, imin(x)=argmin { f (x)-ci-wTG(x)}。
Experiment and result
Quantitative analysis index definition includes:
Define cupping artifact τcup=100 (uM,edge-uM,center)/uM,edge, wherein uM,centerAnd uM,edgeIt is die body center
With the CT values (HU) at edge..
The Cone-Beam CT slice cupping artifact correction of human skull, which is tested, includes:
Experimental image size is 211 × 211.As shown in Figures 2 and 3, the CT value uniformities of image significantly improve after correction.
Analysis result is as shown in table 1.The correction course of a slice takes 1.89 seconds (CPU in this method:i5-2450,RAM:6GB,
GPU:NVIDA GeForce 610M)。
Fig. 3 illustrates the front and back image of horizontal profile correction.It can be seen that observing the cupping artifact of image after correction
It is greatly reduced.
The slice cupping artifact correction of CTP486 die body Cone-Beam CTs, which is tested, includes:
It is tested using CTP486 die bodys.It is as shown in Figure 5 to correct image.CTP486 modules are containing 2% (0- by CT numbers
20H) homogeneous material of water is cast.Image size is 229 × 229.The correction course of one slice takes 1.93 seconds (CPU:
i5-2450,RAM:6GB,GPU:NVIDA GeForce610M)。
The front and back image of horizontal profile correction is as shown in Figure 6.As can be seen from Figure 6 substance slice map CT of the same race before correcting
It is worth uneven, i.e., (shows as cupping artifact), and it is subtle to human eye to use after the method for the present invention cupping artifact to eliminate
Degree.
Table 1:The quantitative analysis of skull die body.Reconstruction image (RI_BC) before cupping artifact correction, after cupping artifact correction
Reconstruction image (RI_AC)
Mouse Bone Cone-Beam CT is sliced cupping artifact correction experiment, specifically includes:
Mouse Bone scattering data is obtained by Hiscan M1000 (Micro-CT).The acquisition of reconstruction image includes 80kVp,
360 projections of 200uA, 30ms.Fig. 8 (a) and Fig. 8 (b) illustrates the sample of the Mouse Bone reconstruction image before and after scatter correction
Axial view.Image size is 339 × 339.The correction course of one slice takes 2.7 seconds (CPU:i5-2450,RAM:6GB,
GPU:NVIDA GeForce610M).As shown in fig. 7, cupping artifact is reduced to 9.8% by this method by 23.8%, analysis result
As shown in table 2.
Table 2:The quantitative analysis of Mouse Bone.Reconstruction image (RI_BC) before cupping artifact correction, after cupping artifact correction
Reconstruction image (RI_AC).
Claims (4)
1. a kind of optimized energy functional based method being carried out cupping artifact correction using reconstruction image in Cone-Beam CT, feature are existed
In described method includes following steps:
Step 1:Obtain the sectioning image after rebuilding;
Step 2:Cupping artifact indicates and builds energy function:
Wherein fpIndicate true artifact-free sectioning image, fsIt indicates to scatter the sectioning image for causing artifact with beam hardening, f is
Reconstruction image without correction, G (x)=(g1(x),…gM(x))TFor smooth basis functions, ciFor constant, meets and work as x ∈ ΩiWhen, ui
(x)=1;WhenWhen, ui(x)=0, variable u=(u1,…uN), c=(c1,…cN) and w=(w1,…wM) all it is convex;
Step 3:Fixed c and u, by solving equationThe minimum value of F (u, c, w) is obtained,Wherein
Step 4:Newer w and u is fixed and uses, withIt is for F (u, c, w) minimum value solution of variable:
Step 5:
Newer w and c is fixed and uses, with u=(u1,…uN)TFor variable F (u, c, w) minimum value solution when,Meet following condition, i.e.,:
Wherein, imin(x)=arg min { f (x)-ci-wTG(x)};
Step 6:If w is stable or iterations are more than 10 times, 7 are thened follow the steps, step 3 is otherwise returned to;
Step 7:Image is after correction
2. a kind of optimized energy carrying out cupping artifact correction using reconstruction image in Cone-Beam CT according to claim 1
Functional based method, which is characterized in that the method is directly facing Cone-Beam CT slice of data, is not carried out to original Cone-Beam CT existing equipment
Any change does not need the prior information of user and measured target.
3. a kind of optimized energy carrying out cupping artifact correction using reconstruction image in Cone-Beam CT according to claim 1
Functional based method, which is characterized in that the cupping artifact of the method is to decomposite to come from reconstruction image f, i.e.,:
Wherein, dsoIt is distance of the radiographic source to rotary shaft, fpIndicate true artifact-free sectioning image, fsIndicate scattering and beam
Hardening causes the sectioning image of artifact, fnIndicate the sectioning image that noise n is formed.
4. a kind of optimized energy carrying out cupping artifact correction using reconstruction image in Cone-Beam CT according to claim 1
Functional based method, which is characterized in that the method is applied to cone-beam CT reconstruction image rectification.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610035576.0A CN105528771B (en) | 2016-01-19 | 2016-01-19 | The bearing calibration of cupping artifact in a kind of Cone-Beam CT using energy function method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610035576.0A CN105528771B (en) | 2016-01-19 | 2016-01-19 | The bearing calibration of cupping artifact in a kind of Cone-Beam CT using energy function method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105528771A CN105528771A (en) | 2016-04-27 |
CN105528771B true CN105528771B (en) | 2018-09-18 |
Family
ID=55770976
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610035576.0A Active CN105528771B (en) | 2016-01-19 | 2016-01-19 | The bearing calibration of cupping artifact in a kind of Cone-Beam CT using energy function method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105528771B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018126434A1 (en) * | 2017-01-06 | 2018-07-12 | 深圳先进技术研究院 | Ct image shadow correction method and apparatus, and electronic device |
EP3704667B1 (en) * | 2017-10-31 | 2021-05-26 | Koninklijke Philips N.V. | Deep learning based motion artifact prediction during magnetic resonance image data acquisition |
CN107997780B (en) * | 2018-01-19 | 2020-11-06 | 重庆大学 | Cone beam CT instantaneous scanning device and reconstruction method |
CN109919868B (en) * | 2019-02-27 | 2022-10-04 | 西北工业大学 | Method for detecting and projection weighting correction of cone beam CT beam hardening curve |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102768759A (en) * | 2012-07-04 | 2012-11-07 | 深圳安科高技术股份有限公司 | Intraoperative CT (Computed Tomography) image beam hardening artifact correction method and device |
CN103020928A (en) * | 2012-11-21 | 2013-04-03 | 深圳先进技术研究院 | Metal artifact correcting method of cone-beam CT (computed tomography) system |
CN103310432A (en) * | 2013-06-25 | 2013-09-18 | 西安电子科技大学 | Computerized Tomography (CT) image uniformization metal artifact correction method based on four-order total-variation shunting |
CN103745440A (en) * | 2014-01-08 | 2014-04-23 | 中国科学院苏州生物医学工程技术研究所 | Metal artifact correction method for CT (computerized tomography) systems |
CN104408753A (en) * | 2014-10-27 | 2015-03-11 | 浙江大学 | Self-adaptive iteration scattering correction method of cone beam CT |
CN104778667A (en) * | 2015-04-14 | 2015-07-15 | 南京邮电大学 | Level-set-based correction method for cupping artifact in cone-beam CT |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102012206714A1 (en) * | 2011-08-10 | 2013-02-14 | Friedrich-Alexander-Universität Erlangen-Nürnberg | Method, arithmetic unit, CT system and C-arm system for the reduction of metal artifacts in CT image datasets |
-
2016
- 2016-01-19 CN CN201610035576.0A patent/CN105528771B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102768759A (en) * | 2012-07-04 | 2012-11-07 | 深圳安科高技术股份有限公司 | Intraoperative CT (Computed Tomography) image beam hardening artifact correction method and device |
CN103020928A (en) * | 2012-11-21 | 2013-04-03 | 深圳先进技术研究院 | Metal artifact correcting method of cone-beam CT (computed tomography) system |
CN103310432A (en) * | 2013-06-25 | 2013-09-18 | 西安电子科技大学 | Computerized Tomography (CT) image uniformization metal artifact correction method based on four-order total-variation shunting |
CN103745440A (en) * | 2014-01-08 | 2014-04-23 | 中国科学院苏州生物医学工程技术研究所 | Metal artifact correction method for CT (computerized tomography) systems |
CN104408753A (en) * | 2014-10-27 | 2015-03-11 | 浙江大学 | Self-adaptive iteration scattering correction method of cone beam CT |
CN104778667A (en) * | 2015-04-14 | 2015-07-15 | 南京邮电大学 | Level-set-based correction method for cupping artifact in cone-beam CT |
Non-Patent Citations (5)
Title |
---|
低能X射线工业CT图像杯状伪影校正;李岭 等;《强激光与粒子束》;20140531;第26卷(第5期);第1-7页 * |
基于SART算法的CL硬化伪影校正方法研究;曹大泉 等;《原子能科学技术》;20140731;第48卷(第7期);第1314-1320页 * |
基于射束衰减网格的锥束CT散射校正方法;张定华 等;《中国机械工程》;20090331;第20卷(第6期);第639-643页 * |
基于自适应点扩散函数的锥束CT散射校正;谢世朋 等;《中国医学影像技术》;20151231;第31卷(第11期);第1763-1767页 * |
基于重建图像全角度前投影的硬化校正方法;徐礼胜 等;《东北大学学报(自然科学版)》;20120831;第33卷(第8期);第1111-1114、1124页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105528771A (en) | 2016-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10111638B2 (en) | Apparatus and method for registration and reprojection-based material decomposition for spectrally resolved computed tomography | |
US8818065B2 (en) | Methods and apparatus for scatter correction for CBCT system and cone-beam image reconstruction | |
CN106491151B (en) | PET image acquisition method and system | |
US20130202080A1 (en) | System and Method for Denoising Medical Images Adaptive to Local Noise | |
US10593070B2 (en) | Model-based scatter correction for computed tomography | |
US9947101B2 (en) | Radiographic image analysis device and method, and recording medium having program recorded therein | |
Zhao et al. | Patient-specific scatter correction for flat-panel detector-based cone-beam CT imaging | |
CN104408758A (en) | Low-dose processing method of energy spectrum CT image | |
CN105528771B (en) | The bearing calibration of cupping artifact in a kind of Cone-Beam CT using energy function method | |
US20160213345A1 (en) | Methods and systems for estimating scatter | |
Ouyang et al. | Effects of the penalty on the penalized weighted least-squares image reconstruction for low-dose CBCT | |
CN104166962A (en) | Cone beam CT scattering correction method by use of scattering nucleus method | |
Chang et al. | Optimization of dose and image quality in adult and pediatric computed tomography scans | |
Xue et al. | Accurate multi-material decomposition in dual-energy CT: A phantom study | |
Shi et al. | X‐ray scatter correction for dedicated cone beam breast CT using a forward‐projection model | |
Lee et al. | A single scatter model for x-ray CT energy spectrum estimation and polychromatic reconstruction | |
Morin et al. | Physical performance and image optimization of megavoltage cone‐beam CT | |
Xie et al. | Scatter correction for cone-beam computed tomography using self-adaptive scatter kernel superposition | |
Pawałowski et al. | Quality evaluation of monoenergetic images generated by dual-energy computed tomography for radiotherapy: A phantom study | |
US11986337B2 (en) | Dose reduction for cardiac computed tomography | |
Marimón et al. | A semi-empirical model for scatter field reduction in digital mammography | |
Thongvigitmanee et al. | Cone-beam CT for breast specimens in surgery: the phantom study | |
Jacobson et al. | Abbreviated on-treatment CBCT using roughness penalized mono-energization of kV-MV data and a multi-layer MV imager | |
CN104778667B (en) | The bearing calibration of cupping artifact in a kind of Cone-Beam CT based on level set | |
WO2022218441A1 (en) | Systems and methods for imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20160427 Assignee: Nanjing causal Artificial Intelligence Research Institute Co., Ltd Assignor: Nanjing Post & Telecommunication Univ. Contract record no.: X2019320000168 Denomination of invention: Cone beam CT cupping artifact correction method by utilizing energy function method Granted publication date: 20180918 License type: Common License Record date: 20191028 |
|
EE01 | Entry into force of recordation of patent licensing contract |