CN107730569A - A kind of medical image artifact bearing calibration and device - Google Patents
A kind of medical image artifact bearing calibration and device Download PDFInfo
- Publication number
- CN107730569A CN107730569A CN201711244741.4A CN201711244741A CN107730569A CN 107730569 A CN107730569 A CN 107730569A CN 201711244741 A CN201711244741 A CN 201711244741A CN 107730569 A CN107730569 A CN 107730569A
- Authority
- CN
- China
- Prior art keywords
- die body
- projection
- data
- artifact
- correction
- 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
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/41—Medical
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The invention discloses a kind of medical image artifact bearing calibration, including：Die body scan data is received, die body scan data is pre-processed, obtains die body measurement projection value；Rebuild die body measurement projection and be worth to the image with artifact；The SECTION EQUATION of the die body is obtained according to the image with artifact；Die body preferred view value is determined according to die body SECTION EQUATION and flying spot equation；Projection value and preferred view value are measured according to die body, using fitting of a polynomial, it is determined that correction coefficient corresponding to each detector；According to correction coefficient, checking list is determined.
Description
Technical field
The present invention relates to field of medical image processing, the bearing calibration of more particularly to a kind of image artifacts and device.
Background technology
In computed tomography (CT) imaging process, artifact be always influence CT reconstructed image qualities it is important because
Element.Artifact refers to exist in CT reconstruction images not existing miscellaneous shadow and interference in actual object, i.e., scanned with reality
The unrelated image of body.Wherein, ring artifact is one of most common artifact.The formation of ring artifact is primarily due to detector list
Member is to the inconsistent of transmitted intensity response and to caused by the non_uniform response of different-energy photon.Ring artifact substantially reduces
The quality of CT images, it is embodied in ring artifact and has had a strong impact on that (CT values represent X ray and pass through tissue for the CT values of CT images
Pad value after being absorbed) accuracy, the uniformities of CT images is influenceed, so as to influence clinical diagnosis.Therefore, remove or most
Mitigating ring artifact to big degree must solve the problems, such as one.
Existing ring artifact bearing calibration, first by the die body with circular cross-section be placed on X-ray tube and X detectors it
Between scanning field in, scan the manikin from multiple directions, obtain multiple views, be then the first projection using scanning result
Information obtains preferred view value i.e. the second projection information by Mathematical treatment (such as fitting or filtering), is then thrown according to first
Shadow information and the second projection information calculate correction coefficient be used for correct from person under inspection's projection information for obtaining there.This method obtains
To preferred view value process dependent on processing to the first projection information, a clear and definite physical model does not produce,
Due to the difference of the various errors and mathematical processing methods of the first projection information of measurement, easily cause caused preferred view value
Inaccuracy, so as to influence to correct result.
The content of the invention
In view of the above-mentioned problems, it is an object of the invention to calculate school using phantom projection data and die body preferred view data
Positive coefficient, and by the correction coefficient for correcting person under inspection's data for projection, correction coefficient can not be accurate in the prior art for solution
Ground corrects the problem of person under inspection's projection information.
To achieve the above object of the invention, one aspect of the present invention provides a kind of medical imaging artifact correction method, including：Connect
Die body scan data is received, die body scan data is pre-processed, obtains die body measurement projection value；Rebuild die body measurement projection value
Obtain the image with artifact；The SECTION EQUATION of the die body is obtained according to the image with artifact；According to die body SECTION EQUATION
Die body preferred view value is determined with flying spot equation；Projection value and preferred view value are measured according to die body, intended using multinomial
Close, it is determined that correction coefficient corresponding to each detector；According to correction coefficient, checking list is determined.
In the present invention, the medical imaging artifact correction method further comprises die body being placed in scanning range
Diverse location, obtain the die body scan data of diverse location.
In the present invention, the die body is cylindroid die body, and the determination die body preferred view value includes：
According to cylindroid die body SECTION EQUATION and flying spot solving simultaneous equation, including following three kinds of situations：
A. ray does not have intersection point with oval, and for the equation group without solution, preferred view value is 0.
B. ray and ellipse are tangent, and the equation group has unique solution, and preferred view value is 0；
C. ray intersects with ellipse, and the equation group has two solutions, is then determined according to the distance of two intersection points when intersecting
Preferred view value.
In the present invention, the medical image artifact bearing calibration, further comprise correcting detected disease by checking list
People's data for projection, the correction include：
Receive and be detected patient's data for projection；Artifact correction table is obtained, the artifact correction table includes one or more schools
Positive coefficient；According to checking list, detected patient's data for projection is corrected, the data for projection after being corrected；According to
Data for projection reconstruction image after correction.
In the present invention, the die body further comprises following shape：Cylindroid, cylinder, circular cone and oval wimble structure.
Present invention also offers a kind of medical image artifact means for correcting, including：Data reception module, for receiving ellipse
Post phantom projection data；Pretreatment module, for being pre-processed to die body scan data；Module is rebuild, for rebuilding medical science
Image；Correction module, for generating artifact correction table, the generation artifact correction table includes：
Cylindroid die body scan data is pre-processed by pretreatment module, obtains cylindroid die body measurement projection value；
Image with artifact is worth to by the measurement projection of reconstruction remodelling cylindroid die body；Institute is obtained according to the image with artifact
State the SECTION EQUATION of die body；Cylindroid die body preferred view is determined according to cylindroid die body SECTION EQUATION and flying spot equation
Value；Projection value and preferred view value are measured according to cylindroid die body, using fitting of a polynomial, it is determined that school corresponding to each detector
Positive coefficient；According to correction coefficient, checking list is determined.
In the present invention, the data reception module can be further used for obtaining is placed in scanning model by cylindroid die body
Enclose the die body scan data of interior diverse location.
In the present invention, the correction module is determined for die body preferred view value, and the determination die body is preferable to be thrown
Shadow value includes：
According to cylindroid die body SECTION EQUATION and flying spot solving simultaneous equation, including following three kinds of situations：
A. ray does not have intersection point with oval, and for the equation group without solution, preferred view value is 0.
B. ray and ellipse are tangent, and the equation group has unique solution, and preferred view value is 0；
C. ray intersects with ellipse, and the equation group has two solutions, is then determined according to the distance of two intersection points when intersecting
Preferred view value.
In the present invention, the correction module can be used for being detected patient's data for projection using checking list correction, described
Correction includes：Receive and be detected patient's data for projection；Artifact correction table is obtained, the artifact correction table includes one or more schools
Positive coefficient；According to checking list, detected patient's data for projection is corrected, the data for projection after being corrected；According to
Data for projection reconstruction image after correction.
In the present invention, the medical image artifact means for correcting further comprises：By the pretreatment module to described
It is detected patient's data for projection to be pre-processed, obtains initial projection data.
Compared with prior art, beneficial effects of the present invention performance is as follows：
First, preferred view value is produced by the modeling to oval die body, makes the calculating of preferred view value have clear and definite foundation,
And the Mathematical treatment to observing projection value is not depended solely on, correction coefficient is further calculated using this preferred view value,
Higher precision can be reached, improve the accuracy of correction.
2nd, for the shape of oval die body closer to body shape, the physical effect (such as scattering) in correction is closer
In the actual conditions of clinical scanning, so as to reach more preferable calibration result.
Brief description of the drawings
Fig. 1 is a kind of schematic diagram of imaging system provided by the invention；
Fig. 2 is a kind of schematic diagram of ring artifact means for correcting provided by the invention；
Fig. 3 is a kind of exemplary process diagram for establishing checking list provided by the invention；
Fig. 4 is a kind of imaging system scanning schematic diagram provided by the invention；
Fig. 5 is a kind of exemplary process diagram of ring artifact correction provided by the invention.
Fig. 1 is marked：100 be imaging system, and 110 be cavity, and 120 be bedstead, and 130 be high pressure generator, and 140 be operation control
Computer equipment processed, 150 be image composer, and 160 be control display device, and 170 be detector, and 180 be radiation generator；
Fig. 2 is marked：200 be ring artifact means for correcting, and 210 be data reception module, and 220 be pretreatment module, and 230 are
Correction module, 240, which attach most importance to, models block, 250 memory modules；
Fig. 4 is marked：401 be X-ray tube, and 402 be X-ray detector, and 403 be die body.
Embodiment
The present invention is described further below by specific embodiment and with reference to accompanying drawing.
Fig. 1 is a kind of schematic diagram of imaging system provided by the invention.In certain embodiments, imaging system 100 can be with
Target is scanned, scan data is obtained and generates associated image.In certain embodiments, imaging system 100 can
To be a medical image system, for example, PET (Positron Emission Tomography) equipment, a CT
(Computed Tomography) equipment, MRI (Magnetic resonance imaging) equipment etc..
In certain embodiments, imaging system 100 can include a high pressure of the bedstead 120, one of cavity 110, one
The raw image composer 150 of operational control computer equipment 140, one of device 130, one and a control display device 160.Chamber
The component for producing and detecting radioactive ray can be housed inside body 110.In certain embodiments, cavity 110 can be received
Hold a radiation generator 180 and a detector 170.Radiation generator 180 can launch radioactive ray.Radioactive ray
It can be transmitted into and be placed at the object in cavity 110, and be received through object by detector 170.As an example, the radiation hair
Raw device 180 can be an X-ray tube.X-ray tube can launch X ray, the thing that the Radiolucent is placed in inside cavity 110
Body, and received by detector 170.In certain embodiments, detector 170 can be circular detector, square detector or arc
Shape detector etc..
Bedstead 120 can support object to be detected (for example, patient to be detected, die body etc.).In certain embodiments,
Bedstead 120 can move in detection process inside cavity 110.As shown in figure 1, in detection process, bedstead 120 can edge
Z-direction moves.According to the needs of detection, patient can lie on the back, prostrate, head in preceding or foot preceding.In some embodiments
In, bedstead 120 can be moved with constant speed inside cavity 110.Bedstead 120 move speed can with sweep time,
The factors such as scanning area are related.In certain embodiments, the speed that bedstead 120 moves can be system default value, can also be by
User sets.
High pressure generator 130 can produce high pressure or heavy current.In certain embodiments, caused high pressure or forceful electric power
Stream can be transmitted to radiation generator 180.
Operational control computer equipment 140 can be with cavity 110, radiation generator 180, detector 170, high pressure generator
130th, bedstead 120, image composer 150 and/or control display device 160 are associated.Can be by direct between the said equipment
Or indirect mode is connected.In certain embodiments, operational control computer equipment 140 can control cavity 110 to rotate
To a certain position.The position can be system default value, can also be set by user (such as doctor, nurse etc.).In some realities
Apply in example, operational control computer equipment 140 can control high pressure generator 130.For example, operational control computer equipment 140
The intensity of voltage or electric current caused by high-voltage generator 130 can be controlled.In certain embodiments, operational control computer equipment
140 can control display device 160.For example, operational control computer equipment 140 can control the parameter related to showing.Institute
Display size, displaying ratio, display order, display number etc. can be included by stating parameter.
Image composer 150 can generate image.In certain embodiments, it is pre- can to carry out image for image composer 150
The operation such as processing, image reconstruction, and/or artifact correction.Image composer 150 can be with detector 170, operational control computer
Equipment 140, display device 160 and/or external data source (not embodied in figure) are associated.In certain embodiments, image generates
Device 150 can receive data from detector 170 or external data source, and based on received data generation image.At some
In embodiment, the image of generation can be sent to display device 160 and be shown by image composer 150.
Display device 160 can show received data or image.Display device 160 can be with operational control computer
Equipment 140 is connected with image composer 150.In certain embodiments, display device 160 can be shown by image composer
150 images generated.In certain embodiments, display device 160 can be to image composer 150 and/or operational control meter
Calculate machine equipment 140 and send instruction.For example, user can set imaging parameters by display device 160, the imaging parameters can be sent out
Deliver to operational control computer equipment 140.The imaging parameters can include data acquisition parameters and image reconstruction parameter etc..
It should be noted that the description for imaging system 100 above, only for convenience of description, can not limit the present invention
System is within the scope of illustrated embodiment.It is appreciated that for those skilled in the art, understanding the principle of the system
Afterwards, modules may be combined in the case of without departing substantially from this principle, or forms subsystem and other moulds
Block connects, to implementing the various modifications and variations in the above method and systematic difference field form and details.
Fig. 2 is a kind of schematic diagram of ring artifact means for correcting provided by the invention.In certain embodiments, ring artifact
Means for correcting 200 may be embodied in image composer 150.
Ring artifact means for correcting 200 can include a pretreatment module 220, one of data reception module 210, one
Correction module 230, one rebuilds module 240 and a memory module 250.Data reception module 210 can receive and measured object
Related data.The data related to measured object can include scan data, essential information (such as name, age, sex, body
High, body weight, medical history etc.), sweep parameter etc..In certain embodiments, the scan data can be collected and passed by detector 170
Deliver to data reception module 210.In certain embodiments, after scan data is collected by detector 170, it can first be sent to and deposit
Module 250 is stored up, then data reception module 210 is sent to by memory module 250.In certain embodiments, data reception module 210
Sweep parameter data can be received from operational control computer equipment 140.In certain embodiments, data reception module 210 can
To receive data (such as essential information of patient) from external data source (not embodied in figure).
Pretreatment module 220 can analyze and process to the data received.Pretreatment module 220 can connect from data
Module 210, memory module 250 and/or external data source is received to receive data and analyzed and processed.In certain embodiments, in advance
Processing module 220 can carry out pretreatment operation to received data.As an example, pretreatment module 220 can carry out sky
Gas corrects, beam hardening correction, the pretreatment such as defocus correction.
Correction module 230 can carry out artifact correction to the data after processing.Correction module 230 can be from pretreatment mould
Block 220, memory module 250 and/or external data source receive data and analyzed and processed.In certain embodiments, straightening die
Block 230 can be used for calculating die body preferred view value.The die body is the model similar with organization of human body.Launched using radioactive source
Radiation exposure die body, and the figure that the different degrees of decay by being arranged on ray on the imaging panel after die body is formed
As being for die body in medical domain to determine the dose of radiation being applied in the Case treatment of reality in patient particular organization
It is commonly used.In certain embodiments, correction module 230 can obtain according to die body scan data and die body preferred view data
To correction coefficient.In certain embodiments, correction module 230 can obtain correction coefficient, and according to correction coefficient to handling it
The scan data of detected patient afterwards is corrected.
Image can be generated and/or handle by rebuilding module 240.In certain embodiments, image processing module can receive
Image is generated as the scan data handled by correction module 230, and according to the scan data of processing.In certain embodiments, weight
Modeling block 240 can carry out image reconstruction.
Memory module 250 can be with data storage, image and/or relevant parameter etc..The data of storage can be various forms
Data.For example, the one or more in the relevant information of numerical value, signal, image, set objective, order, algorithm, program etc..Make
For example, the data for projection after die body scan data, die body preferred view data, artifact correction can be stored in memory module
In 250.In certain embodiments, memory module 250 can include fixed storage system (for example, disk), portable deposit
Storage system etc..In certain embodiments, memory module 250 can store the die body image with artifact, by artifact correction it
Image afterwards etc..Further, memory module 250 can be the interim storage of data, data processing unloading as next time
Data；It can also be the long-term storage of data, that is, store final data processed result.
Fig. 3 is a kind of exemplary process diagram for establishing checking list provided by the invention.
In step 302, data reception module 210 can receive die body scan data.The scan data is alternatively referred to as
Data for projection.The die body can have the structure similar with human body, be cut as an example, have selected herein with ellipse
The cylindroid die body in face, because the section of most people's body region is approximately ellipse, and when the interface of die body is oval, can
, again can be with so that each detector receives the X-ray intensity by die body different-thickness during frame rotation is turned around
Mitigate die body weight by increasing oval major axis/minor axis ratio.In certain embodiments, die body can be by a kind of material institute
The uniform die body of composition, for example, the material can be polypropylene, polyethylene, polytetrafluoroethylene (PTFE) etc..In certain embodiments,
The module can be the structure of big cylindrical sleeve small column, wherein be a kind of material inside small column, small column and big cylinder it
Between be another material.In certain embodiments, die body scan mode can be that die body is placed on into imaging system scanning range
Interior different position, to simulate human body different parts position different in imaging system scanning range.For example, die body is put
It is placed on the center of imaging system scanning range, partial center 50mm, partial center 100mm etc..Then it is different that die body is obtained respectively
Putting position at least once tomographic data and for the later stage image rectification, for example, the die body according to different putting positions
Scan data obtains correction coefficient, and the scan data for being detected patient is corrected according to the correction coefficient in the later stage.
In certain embodiments, Fig. 4 is a kind of imaging system scanning schematic diagram provided by the invention, and 401 be X-ray tube, and 402 penetrate for X
Line detector, 403 be die body.Shown in figure be by die body 403 be placed in imaging system scanning range center position (even if
Oval center is in frame pivot), X-ray tube 401 is to the radiation X ray of X-ray detector 402, X-ray detector
402 detections obtain the scan data that die body is placed in center through the X ray of die body 403.Die body putting position need to make mould
Body is exposed among the X ray propagated with sector 404 from X-ray tube 401 to X-ray detector 402 completely, the X detectors
402 can include multiple detection channels.During scanning, X-ray tube 401 and X-ray detector 402 can be surrounded by picture
The frame pivot of system 100 is rotated, and after X-ray tube 401 and X-ray detector 402 rotate a circle, obtains die body
It is placed in the scan data of center.Then die body is placed in position different in imaging system scanning range, obtained not
With the die body scan data of putting position.For example, partial center 50mm die body scan data, partial center 100mm die body scanning
Data.
In step 304, pretreatment module 220 can handle the die body scan data, obtain the measurement projection of die body
Value.It is pre- that the processing can be that the die body scan data that will receive carries out Air correction, beam hardening correction, defocus correction etc.
Correction, influence of the physical factor to data known to elimination, data after treatment are designated as measuring projection value ProjMi,j,k, its
In,
I=1,2 ..., nChannel；J=1,2 ..., nView；K=1,2 ..., nScan；NChannel is detector list
First number；NView is the projection angle of each tomoscan；NScan is tomoscan number.
Within step 306, correction module 230 can obtain die body preferred view value.Obtain die body preferred view value it
Before, module 240 is rebuild first by the measurement projection value ProjM of each tomoscani,j,kImage reconstruction is carried out, obtains carrying artifact
Image Imagek.Then image Image of the extraction with artifactkIn die body and air interface, and with common ellipse
Equation is fitted the pixel of interface, obtains the elliptic equation in die body section, as shown in formula (1)：
Wherein, (x0,y0) it is the center that die body is put, t is die body major cross-sectional axis and X-axis in Fig. 1 and Y-axis composition
Inclined angle between coordinate system.
For sending the X ray being irradiated to after die body is decayed on detector, the straight line of each ray from X-ray tube
Equation can pass through the focal coordinates (X of X-ray tubeS,YS) and detector I coordinate (XdI,YdI) uniquely determine.The straight line of ray
Shown in equation such as formula (2)：
Ax+By+C=0 (2)
Wherein, A, B, C are to include the parameter of X-ray rotating information, A=YdI-Ys；
B=Xs-XdI；C=XdI*Ys-Xs*YdI。
Simultaneous equations (1) and (2) solve, and have three kinds of situations：
1st, equation group illustrates that now straight line (X ray) does not have intersection point with oval, corresponding die body preferred view value is 0 without solution
2nd, equation group has unique solution, illustrates that now straight line (X ray) and ellipse are tangent, and corresponding die body preferred view value is also
0
3rd, equation group has two groups of solutions, is designated as (X respectively1, Y1) and (X2, Y2), illustrate now straight line (X ray) and oval phase
Hand over.Die body preferred view value can be determined by the distance of two intersection points, as shown in formula (3)：
Wherein μ0It is E to energy for die body0The linear attenuation coefficient of kev ray, E0For the energy size of ray, kev
(kilo electron volt) be energy unit, ProjII,j,kFor die body preferred view value.
Here, by the focal coordinates and detector I coordinates of X-ray tube establish corresponding to linear equation, and with ellipse side
Journey simultaneous solution, obtain die body preferred view value corresponding to detector I.
In step 308, correction module 230 can determine correction coefficient according to preferred view value and measurement projection value.
The measurement projection value detected in some embodiments according to each detector and its corresponding preferred view value, it is determined that each
Correction coefficient corresponding to detector.Specifically, for each detector I, by detector I all measurement projection values
ProjMI,j,kAs independent variable, preferred view value ProjII,j,kAs dependent variable, it is fitted using N rank multinomials, obtains correction system
Number αI,n, as shown in formula (4)：
J=1,2 ..., nView；K=1,2 ..., nScan；NView is the projection angle of each tomoscan；NScan is
Tomoscan number.
In the step 310, correction module 230 can generate checking list according to correction coefficient.Finally, by the checking list
It is stored in any one heretofore described storage device, is further used for artifact correction.In some implementations, can incite somebody to action
Artifact correction table is stored in memory module 250 or external data source in the form of a file.
The description for determining correction coefficient is only specific example above, is not considered as the implementation of unique feasible
Scheme., may be in the feelings without departing substantially from this principle after general principle is understood for one of skill in the art
Under condition, the various modifications and variations in form and details are carried out to embodiment and step, but these modifications and variations
Still within the scope of above description.As an example, the die body structure can also be cylinder, circular cone, oval wimble structure etc..
Fig. 5 is a kind of exemplary process diagram of ring artifact correction provided by the invention.
In step 502, data reception module 210 can receive the scan data of detected patient.The detected material
Can be human body different parts, for example, it may be head, chest, belly, heart, liver, upper limbs, lower limb, vertebra, bone etc.,
Or any combination at above-mentioned position.In certain embodiments, detected material scan data can come from detector 170, storage
Module 250 or external data source.
In step 504, pretreatment module 220 can handle the scan data of detected patient, obtain initial projections value.
The processing can be that the detected material scan data that will be received carries out Air correction, beam hardening correction, defocus correction etc.
Precorrection, influence of the physical factor to data known to elimination, data after treatment are designated as initial projections value ProjOrig.
In step 506, correction module 230 can obtain checking list.The checking list includes one or more corrections
Number.In certain embodiments, the correction coefficient can be determined in step 308 in figure 3.In certain embodiments, institute
Memory module 250 or external data source can be derived from by stating correction coefficient.
In step 508, the projection value after correction module 230 can be corrected according to checking list.Projection value after correction
ProjCorrIIt can be obtained by formula (5)：
Wherein, ProjOrig is initial projections value, αI,nFor correction coefficient corresponding to detector I.
In step 510, module 240 is rebuild according to projection value reconstruction image after correction.In certain embodiments, for scheming
As the algorithm rebuild can include filtered back projection (Filtered Back Projection, FBP), ordered subset expectation maximization
The one or more being worth in (Ordered Subsets Expectation Maximization, OSEM), FDK algorithms etc..
During reconstruction image, the one or more in above-mentioned method for reconstructing can be used, finally obtain artifact-free image.
Above example only illustrates is applied to computer tomography by ring artifact bearing calibration provided by the present invention
Equipment, it will be understood by those skilled in the art that, such as the C-arm system using X ray, combined type medical image system is (for example, group
Box-like positron emission tomography-computer tomography), or use the Laminographic device of other type rays etc.,
Applicable bearing calibration of the present invention and device, the present invention is not done to the type of computer tomography equipment with structure to be had
Body limits.
Being preferable to carry out for the present invention is the foregoing is only, is not intended to limit the invention, for the technology of this area
For personnel, the present invention can have various modifications and variations.Within the spirit and principles of the invention, that is made any repaiies
Change, equivalent substitution, improvement etc., should be included in the scope of the protection.
Claims (10)
1. a kind of medical image artifact bearing calibration, including：
Die body scan data is received,
Die body scan data is pre-processed, obtains die body measurement projection value；
Rebuild die body measurement projection and be worth to the image with artifact；
The SECTION EQUATION of the die body is obtained according to the image with artifact；
Die body preferred view value is determined according to die body SECTION EQUATION and flying spot equation；
Projection value and preferred view value are measured according to die body, using fitting of a polynomial, it is determined that correction system corresponding to each detector
Number；And
According to correction coefficient, checking list is determined.
2. according to the method for claim 1, further comprising the diverse location being placed in die body in scanning range, obtain
The die body scan data of diverse location.
3. according to the method for claim 1, the die body is cylindroid die body, the determination die body preferred view value bag
Include：
According to cylindroid die body SECTION EQUATION and flying spot solving simultaneous equation, including following three kinds of situations：
A. ray does not have intersection point with oval, and for the equation group without solution, preferred view value is 0.
B. ray and ellipse are tangent, and the equation group has unique solution, and preferred view value is 0；
C. ray intersects with ellipse, and the equation group has two solutions, is then determined according to the distance of two intersection points when intersecting preferable
Projection value.
4. the method described in claim 1, further comprise correcting detected patient's data for projection, the correction by checking list
Including：
Receive and be detected patient's data for projection；
Artifact correction table is obtained, the artifact correction table includes one or more correction coefficient；
According to checking list, detected patient's data for projection is corrected, the data for projection after being corrected；And
According to the data for projection reconstruction image after correction.
5. method according to claim 1, it is characterised in that the die body further comprises following shape：Cylindroid, circle
Post, circular cone and oval wimble structure.
6. a kind of medical image artifact means for correcting, including：
One data reception module, for receiving cylindroid phantom projection data；
One pretreatment module, for being pre-processed to die body scan data；
One reconstruction module, for rebuilding medical image；
One correction module, for generating artifact correction table, the generation artifact correction table includes：
Cylindroid die body scan data is pre-processed by the pretreatment module, obtains cylindroid die body measurement projection value；
Image with artifact is worth to by the reconstruction remodelling cylindroid die body measurement projection；
The SECTION EQUATION of the die body is obtained according to the image with artifact；
Cylindroid die body preferred view value is determined according to cylindroid die body SECTION EQUATION and flying spot equation；
Projection value and preferred view value are measured according to cylindroid die body, using fitting of a polynomial, it is determined that corresponding to each detector
Correction coefficient；And
According to correction coefficient, checking list is determined.
7. device according to claim 6, it is characterised in that the data reception module is further used for obtaining ellipse
Post die body is placed in the die body scan data of diverse location in scanning range.
8. device according to claim 7, it is characterised in that the determination die body preferred view value includes：
According to cylindroid die body SECTION EQUATION and flying spot solving simultaneous equation, including following three kinds of situations：
A. ray does not have intersection point with oval, and for the equation group without solution, preferred view value is 0.
B. ray and ellipse are tangent, and the equation group has unique solution, and preferred view value is 0；
C. ray intersects with ellipse, and the equation group has two solutions, is then determined according to the distance of two intersection points when intersecting preferable
Projection value.
9. device according to claim 8, it is characterised in that the correction module is further used for correcting using checking list
Patient's data for projection is detected, the correction includes：
Receive and be detected patient's data for projection；
Artifact correction table is obtained, the artifact correction table includes one or more correction coefficient；
According to checking list, detected patient's data for projection is corrected, the data for projection after being corrected；And
According to the data for projection reconstruction image after correction.
10. the device described in claim 9, it is characterised in that number is projected to the detected patient by the pretreatment module
According to being pre-processed, initial projection data is obtained.
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711244741.4A CN107730569B (en) | 2017-11-30 | 2017-11-30 | Medical image artifact correction method and device |
| US16/205,152 US10922855B2 (en) | 2017-11-30 | 2018-11-29 | Systems and methods for determining at least one artifact calibration coefficient |
| US17/172,113 US11557071B2 (en) | 2017-11-30 | 2021-02-10 | Systems and methods for determining at least one artifact calibration coefficient |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711244741.4A CN107730569B (en) | 2017-11-30 | 2017-11-30 | Medical image artifact correction method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107730569A true CN107730569A (en) | 2018-02-23 |
| CN107730569B CN107730569B (en) | 2021-01-01 |
Family
ID=61220598
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711244741.4A Active CN107730569B (en) | 2017-11-30 | 2017-11-30 | Medical image artifact correction method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107730569B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108670282A (en) * | 2018-03-28 | 2018-10-19 | 上海联影医疗科技有限公司 | A kind of osteosclerosis artifact correction method |
| CN111798534A (en) * | 2020-07-17 | 2020-10-20 | 沈阳先进医疗设备技术孵化中心有限公司 | Image reconstruction method and device, console equipment and CT system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103310432A (en) * | 2013-06-25 | 2013-09-18 | 西安电子科技大学 | Computerized Tomography (CT) image uniformization metal artifact correction method based on four-order total-variation shunting |
| CN103315763A (en) * | 2013-07-04 | 2013-09-25 | 沈阳东软医疗系统有限公司 | Method and device for normal correction of scanning data in imaging equipment |
| CN105118030A (en) * | 2015-08-11 | 2015-12-02 | 上海联影医疗科技有限公司 | Medical image metal artifact correction method and device |
| CN105608721A (en) * | 2016-01-30 | 2016-05-25 | 上海联影医疗科技有限公司 | Computer tomography pseudo shadow correction method and apparatus |
| US20160306056A1 (en) * | 2015-04-17 | 2016-10-20 | General Electric Company | Temperature correction system and method for x-ray detector |
-
2017
- 2017-11-30 CN CN201711244741.4A patent/CN107730569B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103310432A (en) * | 2013-06-25 | 2013-09-18 | 西安电子科技大学 | Computerized Tomography (CT) image uniformization metal artifact correction method based on four-order total-variation shunting |
| CN103315763A (en) * | 2013-07-04 | 2013-09-25 | 沈阳东软医疗系统有限公司 | Method and device for normal correction of scanning data in imaging equipment |
| US20160306056A1 (en) * | 2015-04-17 | 2016-10-20 | General Electric Company | Temperature correction system and method for x-ray detector |
| CN105118030A (en) * | 2015-08-11 | 2015-12-02 | 上海联影医疗科技有限公司 | Medical image metal artifact correction method and device |
| CN105608721A (en) * | 2016-01-30 | 2016-05-25 | 上海联影医疗科技有限公司 | Computer tomography pseudo shadow correction method and apparatus |
Non-Patent Citations (1)
| Title |
|---|
| 张学松 等: "基于Radon变换的CT图像杯状伪影校正", 《CT 理论与应用研究》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108670282A (en) * | 2018-03-28 | 2018-10-19 | 上海联影医疗科技有限公司 | A kind of osteosclerosis artifact correction method |
| CN108670282B (en) * | 2018-03-28 | 2021-12-17 | 上海联影医疗科技股份有限公司 | Method for correcting bone hardening artifact |
| CN111798534A (en) * | 2020-07-17 | 2020-10-20 | 沈阳先进医疗设备技术孵化中心有限公司 | Image reconstruction method and device, console equipment and CT system |
| CN111798534B (en) * | 2020-07-17 | 2024-03-08 | 东软医疗系统股份有限公司 | Image reconstruction method, device, console device and CT system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107730569B (en) | 2021-01-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Srinivasan et al. | Applications of linac-mounted kilovoltage Cone-beam Computed Tomography in modern radiation therapy: A review | |
| CN103961125B (en) | A kind of CT value correcting method for Cone-Beam CT | |
| CN102697514B (en) | Optimal viewing angle is selected to dissect structures visualization and patient skin dose to optimize | |
| US8818058B2 (en) | Method for determining a correction function for correcting computed tomographic numbers of a small target object in a CT image | |
| US9968318B2 (en) | Estimating apparatus, X-ray diagnosis apparatus, and estimating method | |
| US20120278055A1 (en) | Motion correction in radiation therapy | |
| JP6329490B2 (en) | X-ray CT apparatus and image reconstruction method | |
| CN109671128B (en) | Data processing in image reconstruction process, image reconstruction method and device | |
| US8358738B2 (en) | Respiration-correlated radiotherapy | |
| JP2002191591A (en) | Method for correcting radiation hardening to output image obtained by ct equipment | |
| JP2003232855A (en) | Device for preparing tomogram, method for preparing tomogram, and device for radiographic examination | |
| CN104093360A (en) | Nuclear imaging system | |
| JP2004195121A (en) | X-ray ct apparatus and exposure dose calculation method | |
| JPWO2005011502A1 (en) | Radiation tomography equipment | |
| US10098603B2 (en) | Method for estimation and correction of grid pattern due to scatter | |
| Men et al. | Dual-energy imaging method to improve the image quality and the accuracy of dose calculation for cone-beam computed tomography | |
| US20070297575A1 (en) | Systems and methods for determining object position | |
| JP6971537B2 (en) | Treatment planning device and treatment planning method | |
| CN107730569A (en) | A kind of medical image artifact bearing calibration and device | |
| CN104700377B (en) | Obtain the method and apparatus that the beam hardening correction coefficient of beam hardening correction is carried out to computed tomography data | |
| US20230320688A1 (en) | Systems and methods for image artifact mitigation with targeted modular calibration | |
| Loose et al. | On few‐view tomographic reconstruction with megavoltage photon beams | |
| Xiu et al. | An Innovative Beam Hardening Correction Method for Computed Tomography Systems. | |
| US8526757B2 (en) | Imaging system and imaging method for imaging a region of interest | |
| CN110121300B (en) | Zoom radiographic reconstruction |
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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |