CN109118555A - The metal artifacts reduction method and system of computer tomography - Google Patents
The metal artifacts reduction method and system of computer tomography Download PDFInfo
- Publication number
- CN109118555A CN109118555A CN201810925109.4A CN201810925109A CN109118555A CN 109118555 A CN109118555 A CN 109118555A CN 201810925109 A CN201810925109 A CN 201810925109A CN 109118555 A CN109118555 A CN 109118555A
- Authority
- CN
- China
- Prior art keywords
- image
- corrected
- computer tomography
- computer
- metal
- 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/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- 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
Abstract
This application involves a kind of metal artifacts reduction method and systems of computer tomography.The described method includes: obtaining computer tomography image to be corrected, it treats and belongs to the regions of bone tissue intensity value ranges in correcting computer tomographic imaging image and be split, obtain bone segmentation area image, initial prior image is obtained according to computer tomography image to be corrected, final prior image is obtained according to computer tomography image to be corrected, bone segmentation area image and initial prior image, correcting computer tomographic imaging image is obtained according to final prior image and computer tomography image to be corrected.In final prior image include bone tissue information in this method, obtains the higher final prior image of accuracy, the accuracy of metal artifacts reduction can be improved.
Description
Technical field
This application involves Computed tomography fields, pseudo- more particularly to a kind of metal of computer tomography
Shadow bearing calibration, computer tomography metal artifacts reduction system, computer equipment and storage medium.
Background technique
In medical computer tomography (Computed Tomography, CT) imaging, since position implantation is imaged in patient
Metal has very high Decay Rate to X-ray, and there are serious measurement errors for the data that detector receives, and causes computer disconnected
Contain apparent metal artifacts in layer image, seriously affects the quality of computer tomography image.
The metal artifacts reduction method of computer tomography image mainly includes iterative method and projection interpolation method.Iterative method
Due to iteration time-consuming, iterative parameter need to be adjusted, application is restricted.Interpolation method is projected, basic ideas are to data for projection figure
Metal information as in carries out interpolation using the nonmetallic information at its both ends and rebuilds to obtain computer tomography image again, but not
Can effective calibration metal artifact, or even new artifact can be brought into due to interpolation error.
Projection interpolation method is improved, preferable metal artifacts reduction effect is obtained based on the modified method of prior model
Fruit, but still remain limitation.Such as in conventional methods where, using multi-threshold segmentation method, by specifying multiple threshold values manually, will scheme
Other regions as removing metallic region are divided into air, muscle and bone, and every class loading is set to single pixel value；But it is existing not
Foot is: specifying multiple threshold values to need to adjust repeatedly just and can find suitable threshold value and carry out tissue class division manually, due to metal
The CT Distribution value range of artifact may in air, muscle and each region of bone, cause should be soft tissue region, by
Metal artifacts, which seriously affect, is divided into other regions, causes building prior image accuracy low.
Summary of the invention
Based on this, it is necessary to be directed to above-mentioned prior image accuracy low technical problem, provide a kind of computer tomography
Metal artifacts reduction method and system.
A kind of metal artifacts reduction method of computer tomography, comprising the following steps:
Obtain computer tomography image to be corrected；
It treats and belongs to the regions of bone tissue intensity value ranges in correcting computer tomographic imaging image and be split, obtain bone
Tissue segmentation area image；
Initial prior image is obtained according to computer tomography image to be corrected；
It is obtained according to computer tomography image to be corrected, bone segmentation area image and initial prior image final
Prior image；
Correcting computer tomographic imaging image is obtained according to final prior image and computer tomography image to be corrected.
In one embodiment, according to computer tomography image to be corrected, bone segmentation area image and initial
Prior image obtains the step of final prior image, comprising the following steps:
Binary conversion treatment and filtering processing are successively carried out to bone segmentation area image, obtain bone tissue processing image；
It handles image according to bone tissue computer tomography image to be corrected and initial prior image are weighted and are melted
It closes, obtains final prior image.
In one embodiment, binary conversion treatment and filtering processing are successively carried out to bone segmentation area image, obtained
Bone tissue handles the step of image, comprising the following steps:
It sets the pixel value of pixel in the region for belonging to bone tissue intensity value ranges in bone segmentation area image to
1,0 is set by the pixel value of remaining pixel, obtains bone tissue binary image；
According toBone tissue binary image is filtered, is obtained at bone tissue
Image is managed,
Wherein,
W2(i, j) is the pixel that bone tissue handles image the i-th row jth column, and F (i, j) is the judgement of the pixel of the i-th row jth column
Condition, W1(i+a, j+b) is the pixel of bone tissue binary image the i-th+a row jth+b column, and g () is the filtering square of n × n size
Battle array, n is odd number, and the numerical value of filtering matrix is made of 0 and 1,It is the of filtering matrixRow
TheThe element of column, i and j are respectively positive integer, and a is the first increment, and b is the second increment.
In one embodiment, image is handled for computer tomography image to be corrected and initial priori according to bone tissue
The step of image is weighted fusion, obtains final prior image, comprising the following steps:
According to I_prior2 (i, j)=I_original (i, j) × W2(i,j)+I_prior1(i,j)×(1-W2(i,j))
It is weighted fusion, obtains final prior image,
Wherein, I_prior2 (i, j) be final prior image the i-th row jth column pixel, I_original (i, j) be to
The pixel of correcting computer tomographic imaging image the i-th row jth column, W2(i, j) is the picture that bone tissue handles image the i-th row jth column
Element, I_prior1 (i, j) are the pixel of initial prior image the i-th row jth column, and i and j are respectively positive integer.
In one embodiment, the region for belonging to bone tissue intensity value ranges in correcting computer tomographic imaging image is treated
The step of being split, obtaining bone segmentation area image, comprising the following steps:
According to
Bone segmentation area image is obtained,
Wherein, I_bone_artifact (i, j) is the pixel of bone segmentation area image the i-th row jth column, I_
Original (i, j) is the pixel of computer tomography image the i-th row jth to be corrected column, T1For the judgement threshold of metal information
Value, T2For the judgment threshold of bone tissue information, T2<T1, i and j are respectively positive integer.
In one embodiment, the step of initial prior image being obtained according to computer tomography image to be corrected, packet
Include following steps:
It treats metallic region in correcting computer tomographic imaging image to be split, obtains the metal for containing only metal information
Cut zone image；
Correcting computer tomographic imaging image is treated respectively and metal cut zone image carries out forward projection, and respectively
To data for projection to be corrected and metal shadowing's data；
The metallic region that data for projection to be corrected is determined according to metal shadowing's data, according to the metal of data for projection to be corrected
Data outside the both ends of region, the metallic region for treating corrected projection data carry out Interpolate estimation, obtain data for projection to be corrected
Intermediate corrected projection data after removing metal information；
First is obtained according to intermediate corrected projection data and rebuilds computer tomography image, and it is disconnected to rebuild computer to first
Layer image is filtered, and obtains initial prior image.
In one embodiment, correction is obtained according to final prior image and computer tomography image to be corrected to calculate
The step of machine tomographic imaging image, comprising the following steps:
Forward projection is carried out to final prior image, prior image data for projection is obtained, according to prior image data for projection
The metallic region for treating corrected projection data carries out interpolation correction, obtains corrected projection data；
According to corrected projection data obtain second rebuild computer tomography image, by second rebuild computerized tomography at
As image and metal cut zone image co-registration, correcting computer tomographic imaging image is obtained.
A kind of metal artifacts reduction system of computer tomography, comprising:
Image collection module to be corrected, for obtaining computer tomography image to be corrected；
Bone segmentation area image obtains module, belongs to bone tissue for treating in correcting computer tomographic imaging image
The region of intensity value ranges is split, and obtains bone segmentation area image；
Initial prior image obtains module, for obtaining initial priori figure according to computer tomography image to be corrected
Picture；
Final prior image obtains module, for according to computer tomography image to be corrected, bone segmentation region
Image and initial prior image obtain final prior image；
Image collection module is corrected, for obtaining school according to final prior image and computer tomography image to be corrected
Positive computer tomography image.
A kind of computer equipment can be run on a memory and on a processor including memory, processor and storage
Computer program, processor perform the steps of when executing computer program
Obtain computer tomography image to be corrected；
It treats and belongs to the regions of bone tissue intensity value ranges in correcting computer tomographic imaging image and be split, obtain bone
Tissue segmentation area image；
Initial prior image is obtained according to computer tomography image to be corrected；
It is obtained according to computer tomography image to be corrected, bone segmentation area image and initial prior image final
Prior image；
Correcting computer tomographic imaging image is obtained according to final prior image and computer tomography image to be corrected.
A kind of computer readable storage medium is stored thereon with computer program, when computer program is executed by processor
It performs the steps of
Obtain computer tomography image to be corrected；
It treats and belongs to the regions of bone tissue intensity value ranges in correcting computer tomographic imaging image and be split, obtain bone
Tissue segmentation area image；
Initial prior image is obtained according to computer tomography image to be corrected；
It is obtained according to computer tomography image to be corrected, bone segmentation area image and initial prior image final
Prior image；
Correcting computer tomographic imaging image is obtained according to final prior image and computer tomography image to be corrected.
The metal artifacts reduction method and system of above-mentioned computer tomography includes bone tissue in final prior image
Information obtains the higher final prior image of accuracy, the accuracy of metal artifacts reduction can be improved.
Detailed description of the invention
Fig. 1 is the internal structure chart of computer equipment in one embodiment；
Fig. 2 is the flow chart of the metal artifacts reduction method of computer tomography in one embodiment；
Fig. 3 is the flow chart that final prior image obtains in one embodiment；
Fig. 4 is the flow chart of the metal artifacts reduction method of computer tomography in another embodiment；
Fig. 5 is the structural schematic diagram of the metal artifacts reduction system of computer tomography in one embodiment.
Specific embodiment
It is with reference to the accompanying drawings and embodiments, right in order to which the objects, technical solutions and advantages of the application are more clearly understood
The application is further elaborated.It should be appreciated that specific embodiment described herein is only used to explain the application, not
For limiting the application.
The metal artifacts reduction method of computer tomography provided by the present application can be applied to as shown in Figure 1 answer
With in environment, Fig. 1 is the internal structure chart of computer equipment in one embodiment.Wherein, the computer equipment can with but it is unlimited
Then various personal computers, laptop, smart phone, tablet computer and independent server either multiple servers
The server cluster of composition.
In one embodiment, as shown in Fig. 2, Fig. 2 is the metal artifacts school of computer tomography in one embodiment
The flow chart of correction method provides a kind of metal artifacts reduction method of computer tomography in the present embodiment, in this way
Applied to being illustrated for the computer equipment in Fig. 1, comprising the following steps:
Step S210: computer tomography image to be corrected is obtained.
Computer tomography image to be corrected can be include metal information and metal artifacts information CT image.
Step S220: it treats and belongs to the regions of bone tissue intensity value ranges in correcting computer tomographic imaging image and divided
It cuts, obtains bone segmentation area image.
Include bone tissue information in computer tomography image to be corrected, but has the ash of part metals puppet shadow information
Angle value is exactly in bone tissue intensity value ranges, by including bone group in the obtained bone segmentation area image of segmentation
Information and gray value are knitted in the metal artifacts information of bone tissue intensity value ranges.
Step S230: initial prior image is obtained according to computer tomography image to be corrected.
Initial prior image can be used for treating correcting computer tomographic imaging image and be corrected, should but will affect
This is soft tissue area, such as fat and musculature area, other regions are seriously affected and be divided by metal artifacts,
The accuracy of initial prior image is low.
Step S240: according to computer tomography image to be corrected, bone segmentation area image and initial priori figure
As obtaining final prior image.
It should be influence that soft tissue area is subject to reduce, according to bone segmentation area image to initial priori
Image continues further to correct, and obtains final prior image, improves the accuracy of final prior image.
Step S250: correcting computer tomography is obtained according to final prior image and computer tomography image to be corrected
Image.
Correcting computer tomographic imaging image is treated according to the high final prior image of accuracy to be corrected, it is available
More accurate correcting computer tomographic imaging image.
The metal artifacts reduction method of above-mentioned computer tomography includes bone tissue information in final prior image,
The higher final prior image of accuracy is obtained, the accuracy of metal artifacts reduction can be improved.
In one embodiment, as shown in figure 3, Fig. 3 is the flow chart that final prior image obtains in one embodiment, root
Final prior image is obtained according to computer tomography image to be corrected, bone segmentation area image and initial prior image
Step, comprising the following steps:
Step S241: binary conversion treatment and filtering processing are successively carried out to bone segmentation area image, obtain bone tissue
Handle image.
It can be partitioned into bone segmentation region by binary conversion treatment, using filtering processing, optimize obtained bone
Tissue treatment image.
Step S242: according to bone tissue handle image by computer tomography image to be corrected and initial prior image into
Row Weighted Fusion obtains final prior image.
Initial prior image is corrected by way of Weighted Fusion, obtains the higher final priori figure of accuracy
Picture.
The metal artifacts reduction method of above-mentioned computer tomography, by successively successively being carried out to bone segmentation region
Binary conversion treatment, filtering processing and Weighted Fusion, the available higher final prior image of accuracy.
In one embodiment, binary conversion treatment and filtering processing are successively carried out to bone segmentation area image, obtained
Bone tissue handles the step of image, comprising the following steps:
Step S243: by the pixel of pixel in the region for belonging to bone tissue intensity value ranges in bone segmentation area image
Value is set as 1, sets 0 for the pixel value of remaining pixel, obtains bone tissue binary image.
Step S244: according toBone tissue binary image is filtered, is obtained
Image is handled to bone tissue.
Wherein,
W2(i, j) is the pixel that bone tissue handles image the i-th row jth column, and F (i, j) is the judgement of the pixel of the i-th row jth column
Condition, W1(i+a, j+b) is the pixel of bone tissue binary image the i-th+a row jth+b column, and g () is the filtering of n × n size
Matrix, n are odd number, and the numerical value of filtering matrix is made of 0 and 1,It is the of filtering matrix
RowThe element of column, i and j are respectively positive integer, and a is the first increment, and b is the second increment.With
Respectively positive integer.
For example, n value is 3, filtering matrix can be chosen are as follows:
The metal artifacts reduction method of above-mentioned computer tomography obtains bone tissue binary picture in binary conversion treatment
Picture is filtered bone tissue binary image by filtering matrix, obtains bone tissue processing image, filtered bone tissue area
Domain is obvious, in order to improve the accuracy of subsequent final prior image.
In one embodiment, image is handled for computer tomography image to be corrected and initial priori according to bone tissue
The step of image is weighted fusion, obtains final prior image, comprising the following steps:
Step S245:
According to I_prior2 (i, j)=I_original (i, j) × W2(i,j)+I_prior1(i,j)×(1-W2(i,j))
It is weighted fusion, obtains final prior image.
Wherein, I_prior2 (i, j) be final prior image the i-th row jth column pixel, I_original (i, j) be to
The pixel of correcting computer tomographic imaging image the i-th row jth column, W2(i, j) is the picture that bone tissue handles image the i-th row jth column
Element, I_prior1 (i, j) are the pixel of initial prior image the i-th row jth column, and i and j are respectively positive integer.
The metal artifacts reduction method of above-mentioned computer tomography handles image according to bone tissue and obtains Weighted Fusion
Weight, treats correcting computer tomographic imaging image and initial prior image is merged, and obtains including bone tissue information
Final prior image improves the accuracy of final prior image.
In one embodiment, the region for belonging to bone tissue intensity value ranges in correcting computer tomographic imaging image is treated
The step of being split, obtaining bone segmentation area image, comprising the following steps:
Step S221:
According to
Obtain bone segmentation area image.
Wherein, I_bone_artifact (i, j) is the pixel of bone segmentation area image the i-th row jth column, I_
Original (i, j) is the pixel of computer tomography image the i-th row jth to be corrected column, T1For the judgement threshold of metal information
Value, T2For the judgment threshold of bone tissue information, T2<T1, i and j are respectively positive integer.
The metal artifacts reduction method of above-mentioned computer tomography, according to the threshold value of the pixel value range of bone tissue information
Range, it can be determined that area of bone tissue obtains the bone segmentation area image containing bone tissue information, in order to by bone tissue
Packet is contained in final prior image.
In one embodiment, the step of initial prior image being obtained according to computer tomography image to be corrected, packet
Include following steps:
Step S231: it treats metallic region in correcting computer tomographic imaging image and is split, obtain containing only metal
The metal cut zone image of information.
In this step, the metal cut zone image for containing only metal information is obtained, the region of metal is split.
For example, according toObtain metal cut section
Area image, wherein I_metal (i, j) is the pixel of metal cut zone image the i-th row jth column, and I_original (i, j) is
The pixel of computer tomography image the i-th row jth column to be corrected, T1For the judgment threshold of metal information, i and j are positive respectively
Integer.
Step S232: treating correcting computer tomographic imaging image respectively and the progress of metal cut zone image is preceding to throwing
Shadow, and respectively obtain data for projection to be corrected and metal shadowing's data.
It treats correcting computer tomographic imaging image and metal cut zone image does following change process respectively:
I_original_norm=(I_original+Q)/Q,
I_metal_norm=(I_metal+Q)/Q,
Wherein, I_original is computer tomography image to be corrected, and I_original_norm is calculating to be corrected
Machine tomographic imaging limit value image, I_metal are metal cut zone image, and I_metal_norm is metal cut zone limit value figure
Picture, Q are preset limit value coefficient.The value of Q is greater than or equal to 1000, is less than or equal to 5000.
It treats correcting computer tomographic imaging limit value image and carries out forward projection, obtain data for projection and metal to be corrected and throw
Shadow data；Forward projection is carried out to metal cut zone limit value image, obtains metal shadowing's data.When carrying out forward projection, adopt
With identical Ray Tracing Algorithm and parallel beam or fladellum geometric projection mode, data for projection and gold to be corrected are respectively obtained
Belong to data for projection.Data for projection to be corrected and metal shadowing's data are respectively two-dimensional matrix, and size is M × N, wherein M expression pair
Image throw the angle sum of photograph, and N is indicated in each direction for parallel beam used in data for projection or fan beam
Item number, the information after image are received by N number of detection member.
Step S233: determining the metallic region of data for projection to be corrected according to metal shadowing's data, according to projection to be corrected
Data outside the metallic region both ends of data, the metallic region for treating corrected projection data carry out Interpolate estimation, obtain to be corrected
Intermediate corrected projection data after the removal metal information of data for projection.
Nonzero value is metal information in metal shadowing's data, can judge that metal throwing is oriented in operation with simple non-zero
Metallic region in shadow data, so that it is determined that the metallic region in data for projection to be corrected out, remembers under i-th of projection (i.e. projection square
The i-th row of battle array) the beginning and end of metallic region be divided into i_start and i_end, wherein i is greater than 1, and is less than M.
Every a line (the i-th row) in corrected projection data is treated, i-th _ start-1 to the 1st element data are utilized
The value set and i-th _ end+1 data value sets to n-th element, to i-th _ start element to i-th _ end element it
Between data carry out interpolation processing, the intermediate corrected projection data after obtaining interpolation.It is linear that interpolation process method, which can choose,
Interpolation, polynomial interopolation, the methods of B-spline interpolation.
Step S234: first is obtained according to intermediate corrected projection data and rebuilds computer tomography image, to the first weight
It builds computer tomography image to be filtered, obtains initial prior image.
To real projection data reconstruction, the first reconstruction computer tomography image is obtained, which have and is protected
The filtering algorithm of shield picture structure limbic function obtains initial prior image such as based on the Mean Filtering Algorithm of threshold value.Reconstruction side
Method is preferably filtered back-projection method.
The metal artifacts reduction method of above-mentioned computer tomography can rapidly obtain initial prior image, improve
The efficiency of metal artifacts reduction.
In one embodiment, correction is obtained according to final prior image and computer tomography image to be corrected to calculate
The step of machine tomographic imaging image, comprising the following steps:
Step S251: forward projection is carried out to final prior image, prior image data for projection is obtained, according to prior image
The metallic region that data for projection treats corrected projection data carries out interpolation correction, obtains corrected projection data.
When carrying out forward projection to final prior image, using with Ray Tracing Algorithm same in step S232 and parallel
Beam or fladellum geometric projection mode, obtain prior image data for projection.It is thrown using data for projection to be corrected and prior image
The pixel value of the corresponding pixel points of shadow data does division operation, obtains normalized data for projection.In being divided by, if denominator encounter for
0 the case where, enabling denominator value is with lesser positive number.Metal area in normalized data for projection is oriented using metal shadowing's data
Domain, using the data in the data interpolating estimation metallic region of metallic region two sides, process and the Interpolation Process in step S233
It is identical, to obtain the data for projection after interpolation.Utilize the data for projection and prior image data for projection respective pixel after interpolation
The pixel value of point does multiplication operation, and realization goes to normalize, and obtains corrected projection data.
Step S252: second is obtained according to corrected projection data and rebuilds computer tomography image, second is rebuild and is counted
Calculation machine tomographic imaging image and metal cut zone image co-registration, obtain correcting computer tomographic imaging image.
Corrected projection data is rebuild, it is preferred to use filter back-projection reconstruction algorithm obtains I_correct2, and right
I_correct2 carries out limit value and handles, and obtains the second reconstruction computer tomography image.
I_correct=int (I_correct2 × Q-Q)
Wherein, I_correct is the second reconstruction computer tomography image, and int () indicates to carry out data four houses
Five enter to be rounded, and Q is preset limit value coefficient, and Q is identical as the value of Q in step S232.Later, by second rebuild computerized tomography at
As image and metal cut zone image co-registration, i.e., second rebuild the pixel value of the corresponding pixel points of computer tomography image
It is added, the correcting computer tomographic imaging image after obtaining correction of a final proof.
I_final (i, j)=I_correct (i, j)+I_metal (i, j)
Wherein, I_final (i, j) is the pixel of correcting computer tomographic imaging image the i-th row jth column, I_correct
(i, j) is the second pixel for rebuilding computer tomography image the i-th row jth column, and I_metal (i, j) is metal cut zone
The pixel of image the i-th row jth column, i and j are respectively positive integer.
The metal artifacts reduction method of above-mentioned computer tomography is thrown according to the prior image that final prior image obtains
Shadow data carry out interpolation correction, obtain corrected projection data, according to corrected projection data can further with metal cut section
Area image is merged, and the correcting computer tomographic imaging image after obtaining correction of a final proof improves the accurate of metal artifacts reduction
Property.
In another embodiment, as shown in figure 4, Fig. 4 is the metal puppet of computer tomography in another embodiment
The flow chart of shadow bearing calibration, provides a kind of metal artifacts reduction method of computer tomography in the present embodiment, including with
Lower step:
Import a width matrix size be 512 × 512 I_original, I_original be computerized tomography to be corrected at
As image, computer tomography image to be corrected is the computer tomography image containing metal and its artifact.Wherein, it counts
The minimum value of calculation machine tomographic imaging image pixel value is -1000, indicates air；Maximum value reaches+3500, indicates the maximum of metal
Computer tomography pixel value.
T is obtained according to CT image grey level histogram or experience1, T1For the judgment threshold of metal information, T1=+3000.It is logical
It crosses following formula and is partitioned into metal information from I_original, obtain I_metal, I_metal is metal cut zone image:
Wherein, I_metal (i, j) is the pixel of metal cut zone image the i-th row jth column, and I_original (i, j) is
The pixel of computer tomography image the i-th row jth column to be corrected, T1For the judgment threshold of metal information, i and j are positive respectively
Integer.
T is obtained according to CT image grey level histogram or experience2, T2For the judgment threshold of bone tissue information, T2=+800.
Bone tissue information and metal artifacts identical with bone tissue intensity value ranges are partitioned into from I_original by following formula
Information obtains I_bone_artifact, and I_bone_artifact is bone segmentation area image:
Wherein, I_bone_artifact (i, j) is the pixel of bone segmentation area image the i-th row jth column, I_
Original (i, j) is the pixel of computer tomography image the i-th row jth to be corrected column, T1For the judgement threshold of metal information
Value, T2For the judgment threshold of bone tissue information, T2<T1, i and j are respectively positive integer.
It is excessive that there is a situation where data values in a series of subsequent processing in order to prevent, does to I_original and I_metal
Following change process:
I_original_norm=(I_original+1000)/1000
I_metal_norm=(I_metal+1000)/1000
I_original_norm and I_metal_norm is respectively computer tomography limit value image and metal to be corrected
Cut zone limit value image.Then I_original_norm and I_metal_norm minimum value is 0, maximum value 4.5.
When carrying out forward projection to I_original_norm and I_metal_norm, using identical Ray Tracing Algorithm
With parallel beam geometry projection pattern, respectively obtain p_original and p_metal, p_original and p_metal be respectively to
Corrected projection data and metal shadowing's data, above-mentioned data for projection are two-dimensional matrix, and size is 720 × 724, that is, surround image one
The data for projection under 720 directions is uniformly acquired in all ranges, the data for projection number in each direction is 724.
Nonzero value is metal information in p_metal, judges that data for projection p_metal is oriented in operation with simple non-zero
Middle metallic region, so that it is determined that the metallic region in p_original out, remembers (i.e. the i-th row of projection matrix) under i-th of projection
The beginning and end of metallic region is divided into i_start and i_end, and wherein i is greater than 1, and less than 720.
To every a line (the i-th row) in p_original, i-th _ start-1 to the 1st element data value collection are utilized
The data value set with i-th _ end+1 to the 724th element is closed, to i-th _ start element between i-th _ end element
Data carry out interpolation processing, the p_line after obtaining interpolation, p_line be real projection data.Interpolation process method is selected as
Linear interpolation.
Back projection method is filtered to p_line to rebuild, obtains I_correct1, and I_correct1 is the first reconstruction meter
Calculation machine tomographic imaging image carries out I_correct1 to obtain I_ with the filtering algorithm of protection picture structure limbic function
Prior1, I_prior1 are initial prior image, and in the present embodiment, filtering algorithm formula be can choose are as follows:
Wherein, I_prior1 (i, j) is the pixel of initial prior image the i-th row jth column, and I_correct1 (i', j') is
The pixel of first reconstruction the i-th ' row of computer tomography image jth ' column, w is matrix, and v value is that 3, S value is 0.15, sum
(w) for matrix w sum operation, i, j, i' and j ' it is respectively positive integer.
To the I_bone_artifact containing bone tissue information and metal artifacts identical with bone tissue intensity value ranges into
The process of row binary conversion treatment are as follows: enable containing in bone tissue information and metal artifacts identical with bone tissue intensity value ranges region
The pixel value of all pixels point be 1, the pixel value of all pixels point outside region is 0.Image after binaryzation is denoted as W1。
To W1The process being filtered are as follows:
Wherein,
Wherein, W2(i, j) is the pixel that bone tissue handles image the i-th row jth column, and F (i, j) is the pixel of the i-th row jth column
Rule of judgment, W1(i+a, j+b) is the pixel of bone tissue binary image the i-th+a row jth+b column, and g () is n × n size
Filtering matrix, n is odd number, and the numerical value of filtering matrix is made of 0 and 1,It is the of filtering matrixRowThe element of column, a, b, i and j are respectively positive integer.N value is 3, g () selection in the present embodiment
Are as follows:
Utilize W2Instruct the process of the Weighted Fusion of I_original and I_prior1 are as follows:
I_prior2 (i, j)=I_original (i, j) × W2(i,j)+I_prior1(i,j)×(1-W2(i,j))
Wherein, I_prior2 (i, j) be final prior image the i-th row jth column pixel, I_original (i, j) be to
The pixel of correcting computer tomographic imaging image the i-th row jth column, W2(i, j) is the picture that bone tissue handles image the i-th row jth column
Element, I_prior1 (i, j) are the pixel of initial prior image the i-th row jth column, and i and j are respectively positive integer.
When to prior image I_prior2 row forward projection, using with I_original_norm and I_metal_norm into
Same Ray Tracing Algorithm and parallel beam geometry projection pattern when row forward projection, obtain p_prior, p_prior is priori
Image projection data.
Division operation is done using the pixel value of the corresponding pixel points of p_original and p_prior, obtains normalized throwing
Shadow data, are denoted as p_norm1.In being divided by, when denominator encounters the case where being 0, enabling denominator value is with lesser positive number, the present embodiment
In be set as 0.0001.
Metallic region in p_norm1 is oriented using p_metal, estimates metal using the data interpolating of metallic region two sides
Data in region, process is identical as above-mentioned interpolation method process, to obtain p_norm2, p_norm2 is the projection after interpolation
Data.
Multiplication operation is done using the pixel value of p_norm2 and p_prior corresponding pixel points, realization goes to normalize, and obtains p_
Correct2, p_correct2 are corrected projection data.
Backprojection reconstruction is filtered to p_correct2, obtains I_correct2, and I_correct2 is carried out as follows
Processing:
I_correct=int (I_correct2 × 1000-1000)
Wherein, I_correct is the second reconstruction computer tomography image, and int () indicates to carry out data four houses
Five enter to be rounded.
I_correct is merged with I_metal, and the I_final after obtaining correction of a final proof, I_final are disconnected for correcting computer
Layer image.
I_final (i, j)=I_correct (i, j)+I_metal (i, j)
Wherein, I_final (i, j) is the pixel of correcting computer tomographic imaging image the i-th row jth column, I_correct
(i, j) is the second pixel for rebuilding computer tomography image the i-th row jth column, and I_metal (i, j) is metal cut zone
The pixel of image the i-th row jth column, i and j are respectively positive integer.
The metal artifacts reduction method of above-mentioned computer tomography constructs good final prior image, improves correction
As a result；And independent of hardware parameter, applicability and practical.
It should be understood that although each step in the flow chart of Fig. 2 to 4 is successively shown according to the instruction of arrow,
It is these steps is not that the inevitable sequence according to arrow instruction successively executes.Unless expressly stating otherwise herein, these steps
There is no stringent sequences to limit for rapid execution, these steps can execute in other order.Moreover, in Fig. 2 to 4 at least
A part of step may include that perhaps these sub-steps of multiple stages or stage are not necessarily in same a period of time to multiple sub-steps
Quarter executes completion, but can execute at different times, the execution in these sub-steps or stage be sequentially also not necessarily according to
Secondary progress, but in turn or can replace at least part of the sub-step or stage of other steps or other steps
Ground executes.
In one embodiment, as shown in figure 5, Fig. 5 is the metal artifacts school of computer tomography in one embodiment
The structural schematic diagram of positive system provides a kind of metal artifacts reduction system of computer tomography, including image to be corrected
Obtain module 310, bone segmentation area image obtains module 320, initial prior image obtains module 330, final priori figure
As obtaining module 340 and correction image collection module 350, in which:
Image collection module 310 to be corrected, for obtaining computer tomography image to be corrected.
Computer tomography image to be corrected can be include metal information and metal artifacts information CT image.
Bone segmentation area image obtains module 320, belongs to bone for treating in correcting computer tomographic imaging image
The region of tissue intensity value ranges is split, and obtains bone segmentation area image.
Include bone tissue information in computer tomography image to be corrected, but has the ash of part metals puppet shadow information
Angle value is exactly in bone tissue intensity value ranges, by including bone group in the obtained bone segmentation area image of segmentation
Information and gray value are knitted in the metal artifacts information of bone tissue intensity value ranges.Initial prior image obtains module 330, is used for
Initial prior image is obtained according to computer tomography image to be corrected.
Initial prior image can be used for treating correcting computer tomographic imaging image and be corrected, should but will affect
This is soft tissue area, such as bone and musculature area, other regions are seriously affected and be divided by metal artifacts,
The accuracy of initial prior image is low.
Final prior image obtains module 340, for according to computer tomography image to be corrected, bone segmentation area
Area image and initial prior image obtain final prior image.
It should be influence that soft tissue area is subject to reduce, according to bone segmentation area image to initial priori
Image continues further to correct, and obtains final prior image, improves the accuracy of final prior image.
Image collection module 350 is corrected, for obtaining according to final prior image and computer tomography image to be corrected
Take correcting computer tomographic imaging image.
Correcting computer tomographic imaging image is treated according to the high final prior image of accuracy to be corrected, it is available
More accurate correcting computer tomographic imaging image.
The metal artifacts reduction system of above-mentioned computer tomography includes bone tissue information in final prior image,
The higher final prior image of accuracy is obtained, the accuracy of metal artifacts reduction can be improved.
The specific restriction of metal artifacts reduction system about computer tomography may refer to above for calculating
The restriction of the metal artifacts reduction method of machine tomographic imaging, details are not described herein.The metal artifacts of above-mentioned computer tomography
Modules in correction system can be realized fully or partially through software, hardware and combinations thereof.Above-mentioned each module can be hard
Part form is embedded in or independently of in the processor in computer equipment, can also be stored in computer equipment in a software form
Memory in, execute the corresponding operation of above modules in order to which processor calls.
In one embodiment, a kind of computer equipment is provided, which can be server, internal junction
Composition can as shown in FIG. 1, FIG. 1 is the internal structure charts of computer equipment in one embodiment.The computer equipment includes logical
Cross processor, memory and the network interface of system bus connection.Wherein, the processor of the computer equipment is for providing calculating
And control ability.The memory of the computer equipment includes non-volatile memory medium, built-in storage.The non-volatile memories are situated between
Matter is stored with operating system and computer program.The built-in storage is operating system and computer in non-volatile memory medium
The operation of program provides environment.The network interface of the computer equipment is used to communicate with external terminal by network connection.It should
A kind of metal artifacts reduction method of computer tomography is realized when computer program is executed by processor.
It will be understood by those skilled in the art that structure shown in Fig. 1, only part relevant to application scheme is tied
The block diagram of structure does not constitute the restriction for the computer equipment being applied thereon to application scheme, specific computer equipment
It may include perhaps combining certain components or with different component layouts than more or fewer components as shown in the figure.
In one embodiment, a kind of computer equipment is provided, including memory, processor and storage are on a memory
And the computer program that can be run on a processor, processor perform the steps of when executing computer program
Obtain computer tomography image to be corrected；
Initial prior image is obtained according to computer tomography image to be corrected；
It is obtained according to computer tomography image to be corrected, bone segmentation area image and initial prior image final
Prior image；
Correcting computer tomographic imaging image is obtained according to final prior image and computer tomography image to be corrected.
In one embodiment, it is also performed the steps of when processor executes computer program
Binary conversion treatment and filtering processing are successively carried out to bone segmentation area image, obtain bone tissue processing image；
Image is handled according to bone tissue, computer tomography image to be corrected and initial prior image are weighted fusion, obtain most
Whole prior image.
In one embodiment, it is also performed the steps of when processor executes computer program
It sets the pixel value of pixel in the region for belonging to bone tissue intensity value ranges in bone segmentation area image to
1,0 is set by the pixel value of remaining pixel, obtains bone tissue binary image；According to
Bone tissue binary image is filtered, bone tissue processing image is obtained,
Wherein,
W2(i, j) is the pixel that bone tissue handles image the i-th row jth column, and F (i, j) is sentencing for the pixel of the i-th row jth column
Broken strip part, W1(i+a, j+b) is the pixel of bone tissue binary image the i-th+a row jth+b column, and g () is the filter of n × n size
Wave matrix, n are odd number, and the numerical value of filtering matrix is made of 0 and 1,It is the of filtering matrix
RowThe element of column, i and j are respectively positive integer, and a is the first increment, and b is the second increment.
In one embodiment, it is also performed the steps of when processor executes computer program
According to I_prior2 (i, j)=I_original (i, j) × W2(i,j)+I_prior1(i,j)×(1-W2(i,j))
It is weighted fusion, obtains final prior image,
Wherein, I_prior2 (i, j) be final prior image the i-th row jth column pixel, I_original (i, j) be to
The pixel of correcting computer tomographic imaging image the i-th row jth column, W2(i, j) is the picture that bone tissue handles image the i-th row jth column
Element, I_prior1 (i, j) are the pixel of initial prior image the i-th row jth column, and i and j are respectively positive integer.
In one embodiment, it is also performed the steps of when processor executes computer program
According to
Bone segmentation area image is obtained,
Wherein, I_bone_artifact (i, j) is the pixel of bone segmentation area image the i-th row jth column, I_
Original (i, j) is the pixel of computer tomography image the i-th row jth to be corrected column, T1For the judgement threshold of metal information
Value, T2For the judgment threshold of bone tissue information, T2<T1, i and j are respectively positive integer.
In one embodiment, it is also performed the steps of when processor executes computer program
It treats metallic region in correcting computer tomographic imaging image to be split, obtains the metal for containing only metal information
Cut zone image；Correcting computer tomographic imaging image is treated respectively and metal cut zone image carries out forward projection, and
Respectively obtain data for projection to be corrected and metal shadowing's data；The metal of data for projection to be corrected is determined according to metal shadowing's data
Region, according to the data outside the metallic region both ends of data for projection to be corrected, the metallic region for treating corrected projection data is carried out
Interpolate estimation, the intermediate corrected projection data after obtaining the removal metal information of data for projection to be corrected；It is thrown according to centre correction
Shadow data acquisition first rebuilds computer tomography image, is filtered place to the first reconstruction computer tomography image
Reason, obtains initial prior image.
In one embodiment, it is also performed the steps of when processor executes computer program
Forward projection is carried out to final prior image, prior image data for projection is obtained, according to prior image data for projection
The metallic region for treating corrected projection data carries out interpolation correction, obtains corrected projection data；It is obtained according to corrected projection data
Second rebuilds computer tomography image, and the second reconstruction computer tomography image is melted with metal cut zone image
It closes, obtains correcting computer tomographic imaging image.
In one embodiment, a kind of computer readable storage medium is provided, computer program is stored thereon with, is calculated
Machine program performs the steps of when being executed by processor
Obtain computer tomography image to be corrected；
Initial prior image is obtained according to computer tomography image to be corrected；
It is obtained according to computer tomography image to be corrected, bone segmentation area image and initial prior image final
Prior image；
Correcting computer tomographic imaging image is obtained according to final prior image and computer tomography image to be corrected.
In one embodiment, it is also performed the steps of when computer program is executed by processor
Binary conversion treatment and filtering processing are successively carried out to bone segmentation area image, obtain bone tissue processing image；
Image is handled according to bone tissue, computer tomography image to be corrected and initial prior image are weighted fusion, obtain most
Whole prior image.
In one embodiment, it is also performed the steps of when computer program is executed by processor
It sets the pixel value of pixel in the region for belonging to bone tissue intensity value ranges in bone segmentation area image to
1,0 is set by the pixel value of remaining pixel, obtains bone tissue binary image；According to
Bone tissue binary image is filtered, bone tissue processing image is obtained,
Wherein,
W2(i, j) is the pixel that bone tissue handles image the i-th row jth column, and F (i, j) is sentencing for the pixel of the i-th row jth column
Broken strip part, W1(i+a, j+b) is the pixel of bone tissue binary image the i-th+a row jth+b column, and g () is the filter of n × n size
Wave matrix, n are odd number, and the numerical value of filtering matrix is made of 0 and 1,It is the of filtering matrix
RowThe element of column, i and j are respectively positive integer, and a is the first increment, and b is the second increment.
In one embodiment, it is also performed the steps of when computer program is executed by processor
According to I_prior2 (i, j)=I_original (i, j) × W2(i,j)+I_prior1(i,j)×(1-W2(i,j))
It is weighted fusion, obtains final prior image,
Wherein, I_prior2 (i, j) be final prior image the i-th row jth column pixel, I_original (i, j) be to
The pixel of correcting computer tomographic imaging image the i-th row jth column, W2(i, j) is the picture that bone tissue handles image the i-th row jth column
Element, I_prior1 (i, j) are the pixel of initial prior image the i-th row jth column, and i and j are respectively positive integer.
In one embodiment, it is also performed the steps of when computer program is executed by processor
According to
Bone segmentation area image is obtained,
Wherein, I_bone_artifact (i, j) is the pixel of bone segmentation area image the i-th row jth column, I_
Original (i, j) is the pixel of computer tomography image the i-th row jth to be corrected column, T1For the judgement threshold of metal information
Value, T2For the judgment threshold of bone tissue information, T2<T1, i and j are respectively positive integer.
In one embodiment, it is also performed the steps of when computer program is executed by processor
It treats metallic region in correcting computer tomographic imaging image to be split, obtains the metal for containing only metal information
Cut zone image；Correcting computer tomographic imaging image is treated respectively and metal cut zone image carries out forward projection, and
Respectively obtain data for projection to be corrected and metal shadowing's data；The metal of data for projection to be corrected is determined according to metal shadowing's data
Region, according to the data outside the metallic region both ends of data for projection to be corrected, the metallic region for treating corrected projection data is carried out
Interpolate estimation, the intermediate corrected projection data after obtaining the removal metal information of data for projection to be corrected；It is thrown according to centre correction
Shadow data acquisition first rebuilds computer tomography image, is filtered place to the first reconstruction computer tomography image
Reason, obtains initial prior image.
In one embodiment, it is also performed the steps of when computer program is executed by processor
Forward projection is carried out to final prior image, prior image data for projection is obtained, according to prior image data for projection
The metallic region for treating corrected projection data carries out interpolation correction, obtains corrected projection data；It is obtained according to corrected projection data
Second rebuilds computer tomography image, and the second reconstruction computer tomography image is melted with metal cut zone image
It closes, obtains correcting computer tomographic imaging image.
Those of ordinary skill in the art will appreciate that realizing all or part of the process in above-described embodiment method, being can be with
Relevant hardware is instructed to complete by computer program, the computer program can be stored in a non-volatile computer
In read/write memory medium, the computer program is when being executed, it may include such as the process of the embodiment of above-mentioned each method.Wherein,
To any reference of memory, storage, database or other media used in each embodiment provided herein,
Including non-volatile and/or volatile memory.Nonvolatile memory may include read-only memory (ROM), programming ROM
(PROM), electrically programmable ROM (EPROM), electrically erasable ROM (EEPROM) or flash memory.Volatile memory may include
Random access memory (RAM) or external cache.By way of illustration and not limitation, RAM is available in many forms,
Such as static state RAM (SRAM), dynamic ram (DRAM), synchronous dram (SDRAM), double data rate sdram (DDRSDRAM), enhancing
Type SDRAM (ESDRAM), synchronization link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM
(RDRAM), direct memory bus dynamic ram (DRDRAM) and memory bus dynamic ram (RDRAM) etc..
Each technical characteristic of above embodiments can be combined arbitrarily, for simplicity of description, not to above-described embodiment
In each technical characteristic it is all possible combination be all described, as long as however, the combination of these technical characteristics be not present lance
Shield all should be considered as described in this specification.
The several embodiments of the application above described embodiment only expresses, the description thereof is more specific and detailed, but simultaneously
It cannot therefore be construed as limiting the scope of the patent.It should be pointed out that coming for those of ordinary skill in the art
It says, without departing from the concept of this application, various modifications and improvements can be made, these belong to the protection of the application
Range.Therefore, the scope of protection shall be subject to the appended claims for the application patent.
Claims (10)
1. a kind of metal artifacts reduction method of computer tomography, which comprises the following steps:
Obtain computer tomography image to be corrected；
The region for belonging to bone tissue intensity value ranges in the computer tomography image to be corrected is split, bone is obtained
Tissue segmentation area image；
Initial prior image is obtained according to the computer tomography image to be corrected；
According to computer tomography image, the bone segmentation area image and the initial prior image to be corrected
Obtain final prior image；
Correcting computer tomographic imaging is obtained according to the final prior image and the computer tomography image to be corrected
Image.
2. the metal artifacts reduction method of computer tomography according to claim 1, which is characterized in that the basis
The computer tomography image to be corrected, the bone segmentation area image and the initial prior image obtain final
The step of prior image, comprising the following steps:
Binary conversion treatment and filtering processing are successively carried out to the bone segmentation area image, obtain bone tissue processing image；
According to the bone tissue handle image by the computer tomography image to be corrected and the initial prior image into
Row Weighted Fusion obtains final prior image.
3. the metal artifacts reduction method of computer tomography according to claim 2, which is characterized in that described to institute
The step of stating bone segmentation area image and successively carry out binary conversion treatment and filtering processing, obtaining bone tissue processing image, packet
Include following steps:
It sets the pixel value of pixel in the region for belonging to bone tissue intensity value ranges in the bone segmentation area image to
1,0 is set by the pixel value of remaining pixel, obtains bone tissue binary image；
According toThe bone tissue binary image is filtered, the bone group is obtained
Processing image is knitted,
Wherein,
W2(i, j) is the pixel that the bone tissue handles image the i-th row jth column, and F (i, j) is the judgement of the pixel of the i-th row jth column
Condition, W1(i+a, j+b) is the pixel of bone tissue binary image the i-th+a row jth+b column, and g () is n × n size
Filtering matrix, n are odd number, and the numerical value of the filtering matrix is made of 0 and 1,For the filtering matrix
?RowThe element of column, i and j are respectively positive integer, and a is the first increment, and b is the second increment.
4. the metal artifacts reduction method of computer tomography according to claim 2, which is characterized in that the basis
The computer tomography image to be corrected and the initial prior image are weighted and are melted by the bone tissue processing image
The step of closing, obtaining final prior image, comprising the following steps:
According to I_prior2 (i, j)=I_original (i, j) × W2(i,j)+I_prior1(i,j)×(1-W2(i, j)) it carries out
Weighted Fusion obtains the final prior image,
Wherein, I_prior2 (i, j) is the pixel of final prior image the i-th row jth column, and I_original (i, j) is institute
State the pixel of computer tomography image the i-th row jth column to be corrected, W2(i, j) is that the bone tissue handles the i-th row of image the
The pixel of j column, I_prior1 (i, j) are the pixel of initial prior image the i-th row jth column, and i and j are respectively positive integer.
5. the metal artifacts reduction method of computer tomography according to claim 1, which is characterized in that described to institute
It states and belongs to the regions of bone tissue intensity value ranges in computer tomography image to be corrected and be split, obtain bone segmentation
The step of area image, comprising the following steps:
According to
The bone segmentation area image is obtained,
Wherein, I_bone_artifact (i, j) is the pixel of bone segmentation area image the i-th row jth column, I_
Original (i, j) is the pixel of computer tomography image the i-th row jth column to be corrected, T1For sentencing for metal information
Disconnected threshold value, T2For the judgment threshold of bone tissue information, T2<T1, i and j are respectively positive integer.
6. the metal artifacts reduction method of computer tomography according to claim 1, which is characterized in that the basis
The step of computer tomography image to be corrected obtains initial prior image, comprising the following steps:
Metallic region in the computer tomography image to be corrected is split, the metal for containing only metal information is obtained
Cut zone image；
Forward projection is carried out to the computer tomography image to be corrected and the metal cut zone image respectively, and is divided
Data for projection to be corrected and metal shadowing's data are not obtained；
The metallic region that the data for projection to be corrected is determined according to metal shadowing's data, according to the projection number to be corrected
According to metallic region both ends outside data, Interpolate estimation is carried out to the metallic region of the data for projection to be corrected, is obtained described
Intermediate corrected projection data after the removal metal information of data for projection to be corrected；
First is obtained according to the intermediate corrected projection data and rebuilds computer tomography image, is rebuild and is calculated to described first
Machine tomographic imaging image is filtered, and obtains initial prior image.
7. the metal artifacts reduction method of computer tomography according to claim 1, which is characterized in that the basis
The final prior image and the computer tomography image to be corrected obtain the step of correcting computer tomographic imaging image
Suddenly, comprising the following steps:
Forward projection is carried out to the final prior image, prior image data for projection is obtained, is projected according to the prior image
Data carry out interpolation correction to the metallic region of the data for projection to be corrected, obtain corrected projection data；
Second is obtained according to the corrected projection data and rebuilds computer tomography image, and the second reconstruction computer is broken
Layer image and the metal cut zone image co-registration, obtain the correcting computer tomographic imaging image.
8. a kind of metal artifacts reduction system of computer tomography, which is characterized in that the system comprises:
Image collection module to be corrected, for obtaining computer tomography image to be corrected；
Bone segmentation area image obtains module, for belonging to bone tissue in the computer tomography image to be corrected
The region of intensity value ranges is split, and obtains bone segmentation area image；
Initial prior image obtains module, for obtaining initial priori figure according to the computer tomography image to be corrected
Picture；
Final prior image obtains module, for according to the computer tomography image to be corrected, the bone segmentation
Area image and the initial prior image obtain final prior image；
Image collection module is corrected, for obtaining according to the final prior image and the computer tomography image to be corrected
Take correcting computer tomographic imaging image.
9. a kind of computer equipment including memory, processor and stores the meter that can be run on a memory and on a processor
Calculation machine program, which is characterized in that the processor realizes any one of claims 1 to 7 institute when executing the computer program
The step of metal artifacts reduction method for the computer tomography stated.
10. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer program
The metal artifacts reduction method of computer tomography described in any one of claims 1 to 7 is realized when being executed by processor
The step of.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810925109.4A CN109118555B (en) | 2018-08-14 | 2018-08-14 | The metal artifacts reduction method and system of computer tomography |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810925109.4A CN109118555B (en) | 2018-08-14 | 2018-08-14 | The metal artifacts reduction method and system of computer tomography |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109118555A true CN109118555A (en) | 2019-01-01 |
| CN109118555B CN109118555B (en) | 2019-06-14 |
Family
ID=64853262
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810925109.4A Active CN109118555B (en) | 2018-08-14 | 2018-08-14 | The metal artifacts reduction method and system of computer tomography |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109118555B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110335666A (en) * | 2019-05-22 | 2019-10-15 | 平安国际智慧城市科技股份有限公司 | Medical image appraisal procedure, device, computer equipment and storage medium |
| WO2023046193A1 (en) * | 2021-09-27 | 2023-03-30 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for image segmentation |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101044984A (en) * | 2006-03-29 | 2007-10-03 | 西门子公司 | Method for decreasing interstage artifact in heart computer tomography contrast displaying |
| CN102274040A (en) * | 2010-06-10 | 2011-12-14 | 东软飞利浦医疗设备系统有限责任公司 | Post-processing method for removing ring artifacts in CT (computed tomography) image |
| CN103310432A (en) * | 2013-06-25 | 2013-09-18 | 西安电子科技大学 | Computerized Tomography (CT) image uniformization metal artifact correction method based on four-order total-variation shunting |
| US20150030247A1 (en) * | 2013-07-26 | 2015-01-29 | Qualcomm Incorporated | System and method of correcting image artifacts |
| CN104992409A (en) * | 2014-09-30 | 2015-10-21 | 中国科学院苏州生物医学工程技术研究所 | CT image metal artifact correction method |
| CN105528800A (en) * | 2016-01-21 | 2016-04-27 | 上海联影医疗科技有限公司 | Computer tomography imaging pseudo shadow correction method and device |
| CN106296615A (en) * | 2016-08-16 | 2017-01-04 | 广州华端科技有限公司 | CT image corrects the method and system of metal artifacts |
-
2018
- 2018-08-14 CN CN201810925109.4A patent/CN109118555B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101044984A (en) * | 2006-03-29 | 2007-10-03 | 西门子公司 | Method for decreasing interstage artifact in heart computer tomography contrast displaying |
| CN102274040A (en) * | 2010-06-10 | 2011-12-14 | 东软飞利浦医疗设备系统有限责任公司 | Post-processing method for removing ring artifacts in CT (computed tomography) image |
| CN103310432A (en) * | 2013-06-25 | 2013-09-18 | 西安电子科技大学 | Computerized Tomography (CT) image uniformization metal artifact correction method based on four-order total-variation shunting |
| US20150030247A1 (en) * | 2013-07-26 | 2015-01-29 | Qualcomm Incorporated | System and method of correcting image artifacts |
| CN104992409A (en) * | 2014-09-30 | 2015-10-21 | 中国科学院苏州生物医学工程技术研究所 | CT image metal artifact correction method |
| CN105528800A (en) * | 2016-01-21 | 2016-04-27 | 上海联影医疗科技有限公司 | Computer tomography imaging pseudo shadow correction method and device |
| CN106296615A (en) * | 2016-08-16 | 2017-01-04 | 广州华端科技有限公司 | CT image corrects the method and system of metal artifacts |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110335666A (en) * | 2019-05-22 | 2019-10-15 | 平安国际智慧城市科技股份有限公司 | Medical image appraisal procedure, device, computer equipment and storage medium |
| WO2023046193A1 (en) * | 2021-09-27 | 2023-03-30 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for image segmentation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109118555B (en) | 2019-06-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jin et al. | A model-based image reconstruction algorithm with simultaneous beam hardening correction for X-ray CT | |
| Thibault et al. | A three‐dimensional statistical approach to improved image quality for multislice helical CT | |
| La Riviere et al. | Reduction of noise-induced streak artifacts in X-ray computed tomography through spline-based penalized-likelihood sinogram smoothing | |
| De Man et al. | Distance-driven projection and backprojection | |
| US9159122B2 (en) | Image domain de-noising | |
| EP2881039A1 (en) | X-ray computer tomography image pick-up device and image reconstruction method | |
| CN100508891C (en) | Maximum posteriori optimizing image rebuilding method in PET imaging | |
| Zamyatin et al. | Extension of the reconstruction field of view and truncation correction using sinogram decomposition | |
| Peng et al. | An irregular metal trace inpainting network for x‐ray CT metal artifact reduction | |
| Cheng et al. | Learned full-sampling reconstruction from incomplete data | |
| CN109118555B (en) | The metal artifacts reduction method and system of computer tomography | |
| Ma et al. | Generalized Gibbs priors based positron emission tomography reconstruction | |
| CN114387359A (en) | Three-dimensional X-ray low-dose imaging method and device | |
| Michielsen et al. | Patchwork reconstruction with resolution modeling for digital breast tomosynthesis | |
| Kim et al. | Convolutional neural network–based metal and streak artifacts reduction in dental CT images with sparse‐view sampling scheme | |
| Sun et al. | Dual gating myocardial perfusion SPECT denoising using a conditional generative adversarial network | |
| Qiu et al. | New iterative cone beam CT reconstruction software: parameter optimisation and convergence study | |
| Wang et al. | Hybrid pre-log and post-log image reconstruction for computed tomography | |
| CN112204607B (en) | Scattering correction for X-ray imaging | |
| Zeng | Noise‐weighted spatial domain FBP algorithm | |
| Kim et al. | CNN-based CT denoising with an accurate image domain noise insertion technique | |
| Ikuta et al. | A deep recurrent neural network with FISTA optimization for ct metal artifact reduction | |
| US10573029B2 (en) | Fast iterative image reconstruction method for emission tomography | |
| Li et al. | Neural network with maximum entropy constraint for nuclear medicine image restoration | |
| Qiao et al. | Compensating respiratory motion in PET image reconstruction using 4D PET/CT |
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 |