CN105469366A - Analytic method for eliminating metal artifact of CT image - Google Patents

Analytic method for eliminating metal artifact of CT image Download PDF

Info

Publication number
CN105469366A
CN105469366A CN201510808114.3A CN201510808114A CN105469366A CN 105469366 A CN105469366 A CN 105469366A CN 201510808114 A CN201510808114 A CN 201510808114A CN 105469366 A CN105469366 A CN 105469366A
Authority
CN
China
Prior art keywords
image
metal
artifact
function
artifacts
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
Application number
CN201510808114.3A
Other languages
Chinese (zh)
Other versions
CN105469366B (en
Inventor
陈明
李刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University of Science and Technology
Original Assignee
Shandong University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shandong University of Science and Technology filed Critical Shandong University of Science and Technology
Priority to CN201510808114.3A priority Critical patent/CN105469366B/en
Publication of CN105469366A publication Critical patent/CN105469366A/en
Application granted granted Critical
Publication of CN105469366B publication Critical patent/CN105469366B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

The invention provides an analytic method for eliminating a metal artifact of a CT image. With the method, the iterative operation can be avoided. The method comprises: a metal image is segmented from an image reconstructed based on projection data on the multi-energy-spectrum-ray condition; a projection image is solved; an approximate value of an energy spectrum function is searched; an expression of an artifact function is obtained and an artifact image is obtained; the artifact image is subtracted from an original image to obtain a corrected image as shown in an attached drawing. The method is characterized in that the band shaped artifact can be eliminated and the structural feature of original information can be kept; and no priori knowledge of the energy spectrum and original projection data are needed according to the correction method and the reconstruction speed is fast.

Description

A kind of analytic method cutting down CT image metal artifact
Technical field
The present invention relates generally to medical X-ray Computed tomography (ComputedTomography, be called for short CT) field, specifically, relating to medicine CT is detected as in picture, when detected human body a part (such as oral cavity, thoracic cavity etc.) implanted metal thing, there is metal artifacts situation in the CT image of reconstruction.
Background technology
CT detects the acquisition object faultage image that imaging technique can can't harm, irreplaceable vital role is played in medical science auxiliary diagnosis. in the ideal case, the data for projection that the scanning system of CT equipment records is the line integral of material linear attenuation coefficient, the reconstructing system of CT equipment carries out image reconstruction by these line integral data, obtain the CT image that may be used for clinical diagnosis and use. but in actual CT scan, the data for projection that system obtains comprises various error usually, the error that such as metal object causes, this can cause occurring artifact in reconstruction image, reduce picture quality, having a strong impact on diagnostic result. the metal implant in human body has intraoral tooth filling usually, surgical clips and joint prosthesis etc. in body, these metal objects have high Decay Rate, when CT image reconstruction direct by data for projection, rebuild in image and there will be radial or banding artifact etc., generally be called metal artifacts.
The situation comprising metal object is often run in actual CT scan, such as oral cavity CT scan. because the complex genesis of metal artifacts, and the material of metal and geomery can directly cause the power of metal artifacts in image different with distribution, although therefore have some data about metal artifacts school removing method and patent at present, even if but still there is no a kind of method of more common elimination or abatement metal artifacts. for the metal artifacts reduction software for orthopaedics that Philip is released, also gratifying correction result cannot be provided for multiple situation. therefore, the method researching and developing some blanket metal artifacts abatements has very important significance for Medical CT system.
Metal artifacts reduction (Metalartifactreduction, MAR) algorithm is mainly divided into method of interpolation, projection correction, iterative approximation, the method such as bearing calibration and composed correction based on prior image. also there is the patent of a lot of metal artifacts reduction aspect in recent years, such as:
A kind of method for reducing metal artifacts in CT image data set is disclosed by the Chinese patent CN103190928A of Siemens application. the method by rebuilding the CT image data set through metal artifacts reduction and the 2nd CT image data set without metal artifacts reduction, and make use of data weighting method etc.
The Chinese patent CN103186889A applied for by GE company of the U.S. discloses a kind of method and apparatus for reducing the metal artifacts in medical image. and the method considers the artifact that same tomography comprises large metal object and the little metal object generation with different size simultaneously, processes the artifact that in same faultage image, large metal object and little metal object produce in a different manner.
Chinese patent CN103679642A discloses the bearing calibration of a kind of CT image metal artifact, device and CT equipment, metal shadowing's scope that the method original image determination objects interfered corresponding according to Raw projection data causes, then removes the diagnosis object data for projection after objects interfered and Modling model image based on the metal shadowing's data acquisition within the scope of metal shadowing successively; According to the data for projection of model image, secondary correction, reconstruction are carried out to described Raw projection data again, obtain the target image not having metal artifacts.
Chinese patent CN103745440A discloses a kind of CT system metal artifact correction method, create the matrix of coefficients of described CT system according to initiation parameter in the method, based on the object attenuation coefficient f in the process of described system matrix computation CT image reconstruction, optimize described object attenuation coefficient f again, realize the metal artifacts reduction to CT system.
Chinese patent CN103440636A discloses a kind of metal artifacts removing method of CT system, the method obtains new CT scan data through a series of data processing operation, and then rebuild image. high pressure and low pressure CT scan are converted simultaneously, thus realize the minimizing of calculated amount.
Above-mentioned technology has related to various metals artifact correction method, but all do not have method and algorithm to relate to the Analytical Expression of metal artifacts. The present invention gives a kind of analytic method of effective removal CT image metal artifact, avoid iterative computation, save reconstruction time.
Summary of the invention
The present invention, mainly for metal artifacts problem common in Medical CT, proposes a kind of analytic method of effective removal metal artifacts.
The present invention is achieved through the following technical solutions:
The analytic method of abatement CT image metal artifact, mainly comprise: first CT scan system obtains the data for projection of detected object f (x), then image reconstruction is carried out by CT reconstructing system, metal object is partitioned into from reconstruction image, again reprojection is carried out to metal objects, this data for projection is processed, carry out reconstruction again and obtain artifacts, finally from original image, deduct artifacts, obtain target image, namely the image after correcting. it is characterized in that, while eliminating banding artifact, remain the architectural feature of original information; The priori of X-ray energy spectrum and original data for projection is not needed in bearing calibration.
The data for projection P of the detected object described in acquisition, these data obtain in multi-power spectrum ray situation;
CT image P described in acquisition, from reconstructs projection data CT image f e, E is the energy parameter of ray;
Metal object image described in acquisition, goes out the image of metal object by Threshold segmentation (general adopt threshold value be 3000HU) or Level Set Method from the CT Iamge Segmentation rebuild;
The data for projection of the image of the metal object described in acquisition, utilizes projection matrix to obtain its data for projection, and namely its Radon converts R Ψ metal;
Metal artifacts image described in acquisition, the approximate function Asin (ω (E+)) of design power spectrum, the artifact expression formula K=R of proposition -1(P-Rf e0)=-R -1(ln (π 2cosh (ρ R Ψ metal)/(π 2+ 4 ρ 2(R Ψ metal) 2))) to process data R Ψ metal, then rebuild, obtain the approximate image K of artifact, f e0that hypothesis is in certain ENERGY E 0under objective function, ρ is parameter, Ψ metalit is the fundamental function of metallic region;
Image f after namely target image described in acquisition corrects e0, original weight builds image f ededuct artifacts K (x) and obtain f e0.
The present invention relates to a kind of method of abatement CT image metal artifact newly, the method directly obtains the image of image artifacts, very effective. compared with metal artifacts reduction method in the past, beneficial effect of the present invention has: (1) remains the architectural feature of original information while eliminating banding artifact; (2) without the need to the priori of ray energy spectrum; (3) in bearing calibration without the need to original data for projection; (4) adopt analytic reconstruction method, speed of rebuilding is fast.
Accompanying drawing explanation
The two-dimentional lower jaw Voxel Phantom of improvement of Fig. 1 for using in embodiment;
Fig. 2 is the original image rebuild;
Fig. 3 is the artifacts rebuild;
Fig. 4 is the CT image after correcting.
Embodiment
Invention broadly provides a kind of method of cutting down CT image metal artifact, comprising and obtain data for projection, image reconstruction, segmentation metal objects, reprojection, data processing, artifacts reconstruction and obtain the image after correcting. target is the architectural feature that elimination banding artifact retains original information, obtain high-resolution CT image, for clinical diagnosis provides information more accurately. embodiments of the present invention: first obtain the data for projection under multi-power spectrum and carry out image reconstruction; Then carry out the image that Iamge Segmentation obtains metal objects, and ask its Radon to convert; By the analytic method proposed, this Radon transform data is processed, and the data after reconstruction process, obtain artifacts; Deduct artifacts by original image and obtain the CT image after correcting; Finally providing embodiment. concrete steps are as follows:
(1) original CT image is obtained
Carry out detected object scanning the data for projection obtained under multi-power spectrum X ray, by FBP algorithm or BPF algorithm reconstructed projection data, obtain original CT image, containing metal artifacts in this image.
(2) data for projection of metal image is obtained
Usually utilize threshold method to carry out Iamge Segmentation to original image, obtain metal image, the general threshold value adopted is 3000HU herein; And then obtain the data for projection of metal image.
(3) artifacts is obtained
First the approximate of energy spectrum function is found, provide the approximate expression of artifact function, obtain the artifact function determined. because in actual CT scan system, energy spectrogram be unknown, therefore approximate function can only be found according to its common characteristic, the function of this approximate power spectrum should meet: in energy range, integration is 1. such as: approximate function can select rectangular function, sine function, cosine function or other approximate function. according to can the feature of spectrogram, sine function is close to energy spectrum function. through being designed with related parameter, the approximate expression K=R of artifact function can be obtained -1(P-Rf e0). by this expression formula, the data for projection in step 2 is processed; Then ask Radon inverse transformation, utilize FBP or BPF algorithm realization herein, thus obtain artifacts. wherein relevant parameters can utilize some regularization method to ask its optimum solution.
(4) the CT image after correcting is obtained
Deduct artifacts in step 3 by the original image in step 1 and obtain the CT image after correcting.
(5) embodiment
It is the two-dimentional lower jaw Voxel Phantom improved that the present embodiment is chosen, comprise 512 × 512 pixels, the size of each pixel is 0.5mm × 0.5mm, as shown in Figure 1. the sclerotin composition in this body mould and soft tissue etc. are all arranged according to ICRU44 Report Parameters, their X-ray attenuation coefficient is obtained by XCOM software, and implants 3 mercury alloy objects.
This example adopts this body mould to produce artificial projections data. we used multi-power spectrum X ray in emulation. and regarding system parameter is as follows: x-ray source and rotation center distance are 1100mm, x-ray source and detector distance are 1600mm, the length of detector cells is that 0.87. uses the equidistant fan beam CT scan pattern of Circular test to produce the data for projection of this die body, wherein detector cells number is 512, and 360 degree of full scan projection angle numbers are 660.
Directly rebuild image according to step 1 in embodiment by data for projection, there is the artifacts such as obvious band shape in figure, have a strong impact on clinical diagnosis, as shown in Figure 2.
According to step 2 in embodiment above, selected threshold is 3000HU, to the Image Segmentation Using in Fig. 2, and its data for projection of recycling Radon conversion acquisition.
According to step 3 in embodiment above, choosing voltage range in this example is 0.01-1.2Kev, simulation process uses normalization energy spectrogram. power spectrum approximate function π sin (π (E-E0+h)/(2h))/(4h) chosen, wherein E0 is some energy, and h is an adjustable parameter; In this example, E0=0.65Kev, h=0.55Kev. are according to artifact expression formula-R -1(ln (π 2cosh (ρ R Ψ metal)/(π 2+ 4 ρ 2(R Ψ metal) 2))) obtain artifacts as shown in Figure 3, can see that this figure is the presenting of banding artifact in Fig. 2, the parameter choose ρ=2.31. wherein optimized
According to step 4 in above-mentioned embodiment, original reconstruction figure image subtraction artifacts obtains the image after correcting, as shown in Figure 4.

Claims (7)

1. cut down an analytic method for CT image metal artifact, mainly comprise: from reconstruction image, be partitioned into metal object, reprojection is carried out to metal object; The approximate function of design X-ray energy spectrum; By processing the data for projection of metal object, carrying out reconstruction and obtaining artifacts; Finally from original image, deducting artifacts, obtain target image. the feature of described bearing calibration is: (1) employs the analytical expression of artifacts; (2) while eliminating banding artifact, remain the architectural feature of original information; (3) without the need to the priori of X-ray energy spectrum; (4) original data for projection is not used.
2. the analytic method of abatement CT image metal artifact according to claim 1, it is characterized in that: the data for projection obtained from multi-power spectrum ray, rebuild image by FBP algorithm or BPF algorithm, this image contains metal artifacts.
3. the analytic method of abatement CT image metal artifact according to claim 1, is characterized in that: utilize threshold method to be partitioned into metal image from reconstruction image, and utilize projection matrix to obtain the data for projection of metal image.
4. the analytic method of abatement CT image metal artifact according to claim 1, it is characterized in that: design a function be similar to X-ray energy spectrum, it is 1. can select rectangular function or sine function that this function meets integration in energy range. according to can the feature of spectrogram, use sine function Asin (ω (E+)) herein, it is closer to X-ray energy spectrum function.
5. the analytic method of abatement CT image metal artifact according to claim 1, is characterized in that: the approximate function utilizing design, processes the data for projection of metal image, obtain expression formula K (x)=-R of artifact function -1(ln (π 2cosh (ρ R Ψ metal)/(π 2+ 4 ρ 2(R Ψ metal) 2))), i.e. artifacts, wherein ρ is a parameter, is R -1radon inverse transformation,
Ψ metalrepresent the fundamental function (be 1 in region, other are 0) of metallic region.
6. the analytic method of abatement CT image metal artifact according to claim 1, is characterized in that: parameter ρ can utilize regularization method to ask its optimum solution.
7. the analytic method of abatement CT image metal artifact according to claim 1, is characterized in that: build figure from original weight and deduct artifacts, obtains the target image after correcting.
CN201510808114.3A 2015-11-23 2015-11-23 A kind of analytic method of abatement CT image metal artifacts Active CN105469366B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510808114.3A CN105469366B (en) 2015-11-23 2015-11-23 A kind of analytic method of abatement CT image metal artifacts

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510808114.3A CN105469366B (en) 2015-11-23 2015-11-23 A kind of analytic method of abatement CT image metal artifacts

Publications (2)

Publication Number Publication Date
CN105469366A true CN105469366A (en) 2016-04-06
CN105469366B CN105469366B (en) 2018-09-04

Family

ID=55607027

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510808114.3A Active CN105469366B (en) 2015-11-23 2015-11-23 A kind of analytic method of abatement CT image metal artifacts

Country Status (1)

Country Link
CN (1) CN105469366B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106485680A (en) * 2016-10-13 2017-03-08 上海联影医疗科技有限公司 Method for correcting image and device
JP2017221339A (en) * 2016-06-14 2017-12-21 国立大学法人信州大学 X-ray CT image reconstruction method and computer program
CN108937995A (en) * 2017-05-03 2018-12-07 西门子医疗有限公司 For generating the adaptive approach for reducing the CT image data of artifact
CN109146994A (en) * 2018-09-17 2019-01-04 南京航空航天大学 A kind of metal artifacts reduction method towards the imaging of multi-power spectrum X ray CT
CN110599559A (en) * 2018-06-13 2019-12-20 西门子医疗有限公司 Multi-energy metal artifact reduction
CN111110260A (en) * 2019-12-24 2020-05-08 沈阳先进医疗设备技术孵化中心有限公司 Image reconstruction method and device and terminal equipment
CN111415315A (en) * 2020-03-17 2020-07-14 无锡鸣石峻致医疗科技有限公司 Radial collection diffusion weighted imaging motion artifact correction method
WO2023178527A1 (en) * 2022-03-22 2023-09-28 中国科学院深圳先进技术研究院 Generation method and generation apparatus for tumor radiotherapy region

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060062443A1 (en) * 2004-09-23 2006-03-23 General Electric Company System and method for reconstruction of cone beam tomographic projections with missing data
US20070140407A1 (en) * 2005-10-12 2007-06-21 Siemens Corporate Research Inc Reduction of Streak Artifacts In Low Dose CT Imaging through Multi Image Compounding
CN101777177A (en) * 2009-12-29 2010-07-14 上海维宏电子科技有限公司 Attenuation filter-based metal artifact removing mixed reconstruction method for CT images
CN102567958A (en) * 2010-12-31 2012-07-11 东软飞利浦医疗设备***有限责任公司 Image postprocessing method for removing metal artifact from computed tomography (CT) image
US20140363069A1 (en) * 2012-08-23 2014-12-11 General Electric Company System and method for correcting for metal artifacts using multi-energy computed tomography

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060062443A1 (en) * 2004-09-23 2006-03-23 General Electric Company System and method for reconstruction of cone beam tomographic projections with missing data
US20070140407A1 (en) * 2005-10-12 2007-06-21 Siemens Corporate Research Inc Reduction of Streak Artifacts In Low Dose CT Imaging through Multi Image Compounding
CN101777177A (en) * 2009-12-29 2010-07-14 上海维宏电子科技有限公司 Attenuation filter-based metal artifact removing mixed reconstruction method for CT images
CN102567958A (en) * 2010-12-31 2012-07-11 东软飞利浦医疗设备***有限责任公司 Image postprocessing method for removing metal artifact from computed tomography (CT) image
US20140363069A1 (en) * 2012-08-23 2014-12-11 General Electric Company System and method for correcting for metal artifacts using multi-energy computed tomography

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017221339A (en) * 2016-06-14 2017-12-21 国立大学法人信州大学 X-ray CT image reconstruction method and computer program
CN106485680B (en) * 2016-10-13 2017-11-28 上海联影医疗科技有限公司 Method for correcting image and device
CN106485680A (en) * 2016-10-13 2017-03-08 上海联影医疗科技有限公司 Method for correcting image and device
CN108937995B (en) * 2017-05-03 2022-12-16 西门子医疗有限公司 Adaptive method for generating artifact-reduced CT image data
CN108937995A (en) * 2017-05-03 2018-12-07 西门子医疗有限公司 For generating the adaptive approach for reducing the CT image data of artifact
CN110599559B (en) * 2018-06-13 2023-07-14 西门子医疗有限公司 Multi-energy metal artifact reduction
CN110599559A (en) * 2018-06-13 2019-12-20 西门子医疗有限公司 Multi-energy metal artifact reduction
CN109146994B (en) * 2018-09-17 2023-06-20 南京航空航天大学 Metal artifact correction method for multi-energy spectrum X-ray CT imaging
CN109146994A (en) * 2018-09-17 2019-01-04 南京航空航天大学 A kind of metal artifacts reduction method towards the imaging of multi-power spectrum X ray CT
CN111110260A (en) * 2019-12-24 2020-05-08 沈阳先进医疗设备技术孵化中心有限公司 Image reconstruction method and device and terminal equipment
CN111110260B (en) * 2019-12-24 2023-09-26 沈阳先进医疗设备技术孵化中心有限公司 Image reconstruction method and device and terminal equipment
CN111415315A (en) * 2020-03-17 2020-07-14 无锡鸣石峻致医疗科技有限公司 Radial collection diffusion weighted imaging motion artifact correction method
CN111415315B (en) * 2020-03-17 2023-09-26 无锡鸣石峻致医疗科技有限公司 Radial acquisition diffusion weighted imaging motion artifact correction method
WO2023178527A1 (en) * 2022-03-22 2023-09-28 中国科学院深圳先进技术研究院 Generation method and generation apparatus for tumor radiotherapy region

Also Published As

Publication number Publication date
CN105469366B (en) 2018-09-04

Similar Documents

Publication Publication Date Title
CN105469366A (en) Analytic method for eliminating metal artifact of CT image
Gjesteby et al. Metal artifact reduction in CT: where are we after four decades?
Park et al. Metal artifact reduction for polychromatic X-ray CT based on a beam-hardening corrector
Zhang et al. A hybrid metal artifact reduction algorithm for x‐ray CT
Van Slambrouck et al. Metal artifact reduction in computed tomography using local models in an image block‐iterative scheme
Wang et al. Metal artifact reduction in CT using fusion based prior image
Meilinger et al. Metal artifact reduction in cone beam computed tomography using forward projected reconstruction information
US9202296B2 (en) Metal artifacts reduction for cone beam CT
US20060159223A1 (en) Method and apparatus for correcting for beam hardening in CT images
CN102024251B (en) System and method for multi-image based virtual non-contrast image enhancement for dual source CT
Joemai et al. Metal artifact reduction for CT: Development, implementation, and clinical comparison of a generic and a scanner‐specific technique
US10083543B2 (en) Metal artifacts reduction for cone beam CT using image stacking
CN105225208A (en) A kind of computer tomography metal artifacts reduction method and device
Klotz et al. Algorithms for the reduction of CT artifacts caused by metallic implants
Byl et al. Photon‐counting normalized metal artifact reduction (NMAR) in diagnostic CT
CN112017131B (en) CT image metal artifact removing method and device and computer readable storage medium
Wu et al. Iterative CT shading correction with no prior information
Park et al. Machine-learning-based nonlinear decomposition of CT images for metal artifact reduction
Bayaraa et al. A two-stage approach for beam hardening artifact reduction in low-dose dental CBCT
Chen et al. Novel method for reducing high-attenuation object artifacts in CT reconstructions
WO2008065394A1 (en) Method and apparatus for reducing distortion in a computed tomography image
Lee et al. A direct sinogram correction method to reduce metal-related beam-hardening in computed tomography
JP5317612B2 (en) Tomographic image processing apparatus, X-ray CT apparatus, and program
US20060251313A1 (en) Method of producing a cross-sectional image
CN109685871B (en) Image reconstruction method, computer device, and computer-readable storage medium

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