CN108876730A - The method, device and equipment and storage medium of correction of movement artifact - Google Patents
The method, device and equipment and storage medium of correction of movement artifact Download PDFInfo
- Publication number
- CN108876730A CN108876730A CN201810509389.0A CN201810509389A CN108876730A CN 108876730 A CN108876730 A CN 108876730A CN 201810509389 A CN201810509389 A CN 201810509389A CN 108876730 A CN108876730 A CN 108876730A
- Authority
- CN
- China
- Prior art keywords
- data
- projection
- checked object
- mass center
- sweep time
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000009825 accumulation Methods 0.000 claims description 12
- 238000013507 mapping Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000005055 memory storage Effects 0.000 claims description 2
- 238000003384 imaging method Methods 0.000 description 24
- 230000008859 change Effects 0.000 description 8
- 238000001514 detection method Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 239000011159 matrix material Substances 0.000 description 6
- 238000002591 computed tomography Methods 0.000 description 5
- 238000013170 computed tomography imaging Methods 0.000 description 5
- 238000010408 sweeping Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 210000004072 lung Anatomy 0.000 description 3
- 241000406668 Loxodonta cyclotis Species 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012636 positron electron tomography Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000010349 pulsation Effects 0.000 description 2
- 230000002792 vascular Effects 0.000 description 2
- 241000167880 Hirundinidae Species 0.000 description 1
- 238000012879 PET imaging Methods 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 210000005252 bulbus oculi Anatomy 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000004424 eye movement Effects 0.000 description 1
- 230000002496 gastric effect Effects 0.000 description 1
- 230000005176 gastrointestinal motility Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008855 peristalsis Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000001568 sexual effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000009885 systemic effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/77—Retouching; Inpainting; Scratch removal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10104—Positron emission tomography [PET]
-
- 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/30196—Human being; Person
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
This application discloses the method, device and equipment and storage medium of a kind of correction of movement artifact, the method includes：Obtain the data for projection of checked object；By tracking the mass center of acquired data for projection, the motion feature that the mass center of the checked object changes with sweep time is mapped out；Based on the motion feature, the data for projection in acquired data for projection with motion feature is corrected, wherein the data for projection with motion feature includes：The checked object is scanned resulting data in its center of mass motion；Image reconstruction is carried out according to the data for projection after correction, obtains the reconstruction image of the checked object.Implement the application, there is the data for projection of motion feature by correcting, influence of the gained scan data to reconstruction image when checked object movement can be effectively reduced, so that the motion artifacts in reconstruction image are improved.
Description
Technical field
This application involves the method, device and equipment of technical field of image processing more particularly to correction of movement artifact and deposit
Storage media.
Background technique
CT is to be scanned using X-ray beam to certain position of subject, and pass through a series of places to scanning result
Reason, obtains a kind of technology of CT image being made of pixel.Wherein, CT image is a kind of reconstruction image, also known as CT reconstruction image.
By taking subject is human body as an example, when certain positions to human body carry out CT scan, since human body autokinetic movement or physiological are transported
Dynamic reason is likely to occur motion artifacts (structural fuzzy or ghost image of scanned position) in reconstruction image.
And the appearance of motion artifacts can reduce the clarity of CT reconstruction image, and then influence to be obtained according to the CT reconstruction image
The accuracy of the diagnostic result obtained.
Summary of the invention
The application provides the method, device and equipment of correction of movement artifact, poor to solve existing reconstruction image clarity
The problem of.
According to the embodiment of the present application in a first aspect, provide a kind of method of correction of movement artifact, the method includes：
Obtain the data for projection of checked object；
By tracking the mass center of acquired data for projection, the mass center for mapping out the checked object changes with sweep time
Motion feature；
Based on the motion feature, the data for projection in acquired data for projection with motion feature is corrected, wherein tool
The data for projection for having motion feature includes：The checked object is scanned resulting data in its center of mass motion；
Image reconstruction is carried out according to the data for projection after correction, obtains the reconstruction image of the checked object.
In one embodiment, acquired data for projection includes multiple groups data for projection, and one group of data for projection is tested pair
As scanned primary resulting data, the mass center by tracking acquired data for projection maps out the checked object
The motion feature that changes with sweep time of mass center, including：
Determine the corresponding location parameter of each group data for projection and sweep time；
By every group of data for projection and its corresponding location parameter, this group of data for projection is converted in its corresponding sweep time
Mass center；
By every group of data for projection corresponding sweep time and its mass center in the sweep time, it is mapped as described tested pair
As the mass center in corresponding sweep time；
Mass center to the checked object in each sweep time carries out data fitting, obtains the mass center phase of the checked object
To the variation relation of sweep time.
In one embodiment, the location parameter is the location parameter in the channel of acquired projections data.
In one embodiment, the mass center by tracking acquired data for projection, maps out the checked object
The motion feature that changes with sweep time of mass center, further include：
In one embodiment, the mass center by tracking acquired data for projection, maps out the checked object
The motion feature that changes with sweep time of mass center, further include：
If during the checked object is scanned, the equipment for carrying the checked object is kept in motion, then counts
The time difference of maximum scan time determined by calculating and minimum sweep time；
The checked object is obtained in the accumulation result of the mass center of each sweep time；
It calculates accumulation result and the ratio of the time difference is equipment error；
The equipment error is eliminated from the variation relation, obtains revised variation relation.
In one embodiment, the mass center by tracking acquired data for projection, maps out the checked object
The motion feature that changes with sweep time of mass center, further include：
Based on revised variation relation, the main movement shaft of the mass center of the checked object is determined, with current coordinate system
The angle value of axis of abscissas；
By revised variation relation, coordinate rotation is carried out according to the angle value, obtains postrotational variation relation.
In one embodiment, described to be based on the motion feature, correcting in acquired data for projection has movement special
The data for projection of sign, including：
Based on the corresponding relationship and the motion feature between sweep time and scan position, the checked object is determined
Motion intense degree of the mass center in each scan position；
According to the mass center of the checked object in the motion intense degree of each scan position, determine the checked object each
The reconstruction weight of the corresponding data for projection of scan position, wherein the more violent scan position of center of mass motion, corresponding data for projection
Reconstruction weight it is smaller；
Obtain the corresponding reconstruction weight product of the data for projection of each scan position, the data for projection after being corrected.
In one embodiment, the corresponding relationship based between sweep time and scan position and the movement are special
Sign, determines motion intense degree of the mass center in each scan position of the checked object, including：
According to the corresponding relationship between sweep time and scan position, the motion feature is converted into the checked object
The movement relation of relatively each scan position of mass center；
Based on the movement relation, motion intense degree of the mass center in each scan position of the checked object is calculated.
In one embodiment, the corresponding relationship based between sweep time and scan position and the movement are special
Sign, determines motion intense degree of the mass center in each scan position of the checked object, including：
Based on the motion feature, motion intense degree of the mass center in each sweep time of the checked object is calculated；
Based on the corresponding relationship between sweep time and scan position, by the mass center of the checked object in each sweep time
Motion intense degree is converted to motion intense degree of the mass center in each scan position of the checked object.
According to the second aspect of the embodiment of the present application, a kind of device of correction of movement artifact is provided, including：
Data acquisition module, for obtaining the data for projection of checked object；
Characteristic extracting module maps out the checked object for the mass center by the acquired data for projection of tracking
The motion feature that mass center changes with sweep time；
Data correction module, for being based on the motion feature, correcting has motion feature in acquired data for projection
Data for projection, wherein include with the data for projection of motion feature：The checked object is scanned institute in its center of mass motion
The data obtained；
Image reconstruction module obtains the checked object for carrying out image reconstruction according to the data for projection after correction
Reconstruction image.
In one embodiment, acquired data for projection includes multiple groups data for projection, and one group of data for projection is tested pair
As scanned primary resulting data, the characteristic extracting module includes：
Parameter determination module, for determining the corresponding location parameter of each group data for projection and sweep time；
Centroid calculation module, for being converted to this group of data for projection for every group of data for projection and its corresponding location parameter
In the mass center of its corresponding sweep time；
Mass center mapping block, for by every group of data for projection corresponding sweep time and its sweep time matter
The heart is mapped as the checked object in the mass center of corresponding sweep time；
First fitting module carries out data fitting for the mass center to the checked object in each sweep time, obtains institute
State the variation relation of the mass center relative scanning time of checked object.
In one embodiment, the location parameter is the location parameter in the channel of acquired projections data.
In one embodiment, the characteristic extracting module further includes：
Characteristic determination module, for during the checked object is scanned, the equipment for carrying the checked object to be in
When stationary state, with the variation relation of the mass center relative scanning time of the checked object, the mass center of the checked object is indicated
In the motion feature of each sweep time；
Time difference determining module, at the checked object scanned period, the equipment for carrying the checked object
When motion state, the time difference of identified maximum scan time and minimum sweep time are calculated；
Mass center accumulator module, for obtaining the checked object in the accumulation result of the mass center of each sweep time；
Error calculating module, the ratio for calculating accumulation result and the time difference are equipment error；
Error correction module obtains revised variation and closes for eliminating the equipment error from the variation relation
System.
In one embodiment, the characteristic extracting module further includes：
Angle determining module determines the main motion of the mass center of the checked object for being based on revised variation relation
Axis, the angle value with the axis of abscissas of current coordinate system；
Coordinate rotary module, for carrying out revised variation relation coordinate rotation according to the angle value, being revolved
Variation relation after turning.
In one embodiment, the Data correction module includes：
Movement degree determining module, for based between sweep time and scan position corresponding relationship and the movement it is special
Sign, determines motion intense degree of the mass center in each scan position of the checked object；
Rebuild weight determination module, for the mass center according to the checked object each scan position motion intense journey
Degree, determines the checked object in the reconstruction weight of the corresponding data for projection of each scan position, wherein center of mass motion is more violent
The reconstruction weight of scan position, corresponding data for projection is smaller；
Data correction submodule, the corresponding reconstruction weight product of data for projection for obtaining each scan position, obtains
Data for projection after to correction.
In one embodiment, the movement degree determining module is configured as：
According to the corresponding relationship between sweep time and scan position, the motion feature is converted into the checked object
The movement relation of relatively each scan position of mass center；
Based on the movement relation, motion intense degree of the mass center in each scan position of the checked object is calculated.
In one embodiment, the movement degree determining module is configured as：
Based on the motion feature, motion intense degree of the mass center in each sweep time of the checked object is calculated；
Based on the corresponding relationship between sweep time and scan position, by the mass center of the checked object in each sweep time
Motion intense degree is converted to motion intense degree of the mass center in each scan position of the checked object.
According to the third aspect of the embodiment of the present application, a kind of computer equipment is provided, including：
Processor；
The memory of storage processor executable instruction；
Wherein, the processor is coupled in the memory, for reading the program instruction of the memory storage, and makees
For response, the operation in method as described above is executed.
According to the fourth aspect of the embodiment of the present application, one or more machine readable storage mediums are provided, are stored thereon with
Instruction, when executed by one or more processors, so that computer equipment executes the operation in method as described above.
Using the embodiment of the present application, the mass center of the data for projection by tracking checked object can map out checked object
The motion feature that changes with sweep time of mass center, the motion feature is then based on, to having in acquired data for projection
The data of motion feature are corrected, then carry out image reconstruction based on the data for projection after correction.What it is due to correction is mass center fortune
Resulting data for projection when dynamic, therefore, shadow of the gained data for projection to reconstruction image when checked object movement can be effectively reduced
It rings, and then based on the data for projection reconstruction image after correction, the motion artifacts in reconstruction image can be made to be improved, reach and mention
The purpose of the quality of high reconstruction image, and accurate foundation is provided for the subsequent diagnosis carried out based on reconstruction image.
It should be understood that above general description and following detailed description be only it is exemplary and explanatory, not
The application can be limited.
Detailed description of the invention
The drawings herein are incorporated into the specification and forms part of this specification, and shows the implementation for meeting the application
Example, and together with specification it is used to explain the principle of the application.
Figure 1A is the reconstruction image containing motion artifacts shown in one exemplary embodiment of the application；
Figure 1B is the flow chart of the method for the correction of movement artifact shown in one exemplary embodiment of the application；
Fig. 1 C is the architecture diagram of the scanning device shown in one exemplary embodiment of the application；
Fig. 1 D is the process schematic of the adjacent twice sweep shown in one exemplary embodiment of the application；
Fig. 1 E is the schematic diagram of the parallel beam shown in one exemplary embodiment of the application；
Fig. 2 is the flow chart of the method for the correction of movement artifact shown in the application another exemplary embodiment；
Fig. 3 is the flow chart of the method for the correction of movement artifact shown in the application another exemplary embodiment；
Fig. 4 A is the movement of the mass center relative scanning time change of the checked object shown in one exemplary embodiment of the application
Feature schematic diagram；
Fig. 4 B is the movement relation of the mass center versus scan location of the checked object shown in one exemplary embodiment of the application
Schematic diagram；
Fig. 4 C is to determine the schematic diagram for rebuilding weight shown in one exemplary embodiment of the application；
Fig. 5 is the hardware structural diagram of the computer equipment shown in one exemplary embodiment of the application；
Fig. 6 is the block diagram of the device of the correction of movement artifact shown in one exemplary embodiment of the application.
Specific embodiment
It is only to be not intended to be limiting the application merely for for the purpose of describing particular embodiments in term used in this application.
It is also intended in the application and the "an" of singular used in the attached claims, " described " and "the" including majority
Form, unless the context clearly indicates other meaning.It is also understood that term "and/or" used herein refers to and wraps
It may be combined containing one or more associated any or all of project listed.
The invention relates to motion artifacts, refer to imaging system scanning checked object during, checked object
Macroscopic motion caused by motion artifacts, autokinetic movement artifact and physiological motion artifacts can be divided into, wherein autokinetic movement is pseudo-
Shadow, which refers to, not to be had periodically, and checked object is capable of artifact caused by the movement of autonomous control, such as：It swallows, Rotation of eyeball, limb
Artifact caused by body movement etc.；Physiological motion artifacts refer to artifact caused by the non-autonomous movement not controlled by checked object,
Such as：The aperiodic motion of artifact caused by the cycle movement of heartbeat, vascular pulsation etc or gastrointestinal motility etc
Caused by artifact.
In practical application, pass through CT (Computed Tomography, i.e. CT scan), PET
(Positron Emission computed Tomography, PET-Positron emission computed tomography scanning) and PET-CT
Etc imaging system when scanning checked object, due to the macroscopic motion of checked object, in the reconstruction image that imaging system is constituted,
There may be above-mentioned motion artifacts.Such as：Shown in Figure 1A, be abdomen CT reconstruction image, the structural fuzzy of the CT reconstruction image,
Clarity is poor, will affect the accuracy of diagnostic result.
Above-mentioned imaging system can carry out image based on the combined result of multiple groups data for projection during reconstruction image
It rebuilds.Therefore, the data for projection with motion feature is more in the data for projection that reconstruction image uses, and reconstruction image moves
A possibility that artifact, is higher, and the severity of motion artifacts is higher.In view of this, the embodiment of the present application proposes a kind of correction of movement
The method of artifact, by imaging system, in the center of mass motion of checked object, (checked object is in movement to this method when center of mass motion
State), the resulting data of checked object are scanned, referred to as with the data for projection of motion feature.
In practical application, the mass center of the data for projection by tracking checked object can map out the mass center of checked object
With the motion feature that sweep time changes, it is then based on the motion feature, it is special to having movement in acquired data for projection
The data of sign are corrected, then carry out image reconstruction based on the data for projection after correction.By when that is center of mass motion institute of correction
The data for projection obtained, therefore, influence of the gained data for projection to reconstruction image when checked object movement can be effectively reduced, in turn
Based on the data for projection reconstruction image after correction, the motion artifacts in reconstruction image can be made to be improved, reach raising and rebuild
The purpose of the quality of image, and accurate foundation is provided for the subsequent diagnosis carried out based on reconstruction image.It is described in detail below in conjunction with attached drawing
The method of correction of movement artifact provided by the present application.
Figure 1B is please referred to, Figure 1B is the process of the method for the correction of movement artifact shown in one exemplary embodiment of the application
Figure, the embodiment may comprise steps of S101-S104：
Step S101, the data for projection of checked object is obtained.
Step S102, by tracking the mass center of acquired data for projection, the mass center of the checked object is mapped out with sweeping
Retouch the motion feature of time change.
Step S103, it is based on the motion feature, corrects the projection number in acquired data for projection with motion feature
According to, wherein include with the data for projection of motion feature：The checked object is scanned resulting number in its center of mass motion
According to.
Step S104, image reconstruction is carried out according to the data for projection after correction, obtains the reconstruction image of the checked object.
The present embodiments relate to data for projection, be referred to as raw data, imaging system signal source to tested
After object emits corresponding signal, the signal of subject release is received by the detector of imaging system and received signal is carried out
It is generated after corresponding conversion.Imaging system mentioned herein can be CT imaging system, PET imaging system, PET-CT imaging system,
Or the imaging system of other combined result reconstruction images based on multiple groups data for projection.Below by taking CT system as an example, in conjunction with figure
1C specifically introduces the generating process of lower data for projection：
CT imaging system may include scanning device as shown in Figure 1 C, which may include scanning support 11, peace
The equipment 15 of X-ray emitter and X-ray detector 16 and carrying checked object on scanning support 11.Wherein, X is penetrated
Line transmitter may include x-ray source 12 and collimator set 13, equipment 15 can be it is scanning bed, X-ray detector 16 can be with
Channel including multiple detection X-rays.
In practical application, equipment 15 can advance along the direction axis R for being parallel to the direction z, and checked object (is not shown in Fig. 1 C
Take corresponding detection zone to out), scanning support 11 can be rotated around axis R, drive X-ray emitter to multiple positions, to detection
The checked object in area takes multiple scan.
Every time when scanning, from X-ray emitter to the checked object emitting x-ray 14 of detection zone, in emitting x-ray
When 14, the X-ray that x-ray source 12 generates is adjusted to parallel beam, fan-ray beam or taper ray by collimator set 13
Shu Hou, then emit to checked object.
The X-ray beam 14 of transmitting can penetrate checked object, be transmitted on X-ray detector 16, by X-ray detector 16
It receives X-ray and becomes electric signal through photoelectric conversion, then (Analog/Digital Converter, analog/digital turn through ADC
Parallel operation) be converted to digital signal.These digital signals are properly termed as raw data, or alternatively referred to as data for projection.
After X-ray detector 16 exports data for projection, the embodiment of the present application is resulting in order to avoid obscuring each scanning
Data for projection can be grouped by data for projection, and one group of data for projection is scanning checked object once resulting data, correspond to this
At least one of in sweep time of scanning, scan position and scanning angle.Wherein, scanning angle can be by Current Scan number
And the rotation angle of adjacent twice sweep determines.Fig. 1 D shows adjacent twice sweep, and scan position twice is respectively
W1 and W2, rotation angle are θ, and scanning can project as shown in the figure parallel to checked object by X-ray emitter 21 every time
X-ray beam receives X-ray beam by X-ray detector 22 and generates data for projection.
Further, it is also possible to the data for projection with group is distinguished with the different channels in X-ray detector, it will be each in same group
Data for projection corresponds to each channel storage, for example, the location parameter in corresponding each channel, stores the data for projection of each channel output.
When due to actual scanning detected object, scanning support can drive X-ray emitter and X-ray detector along Fig. 1 C
Shown in the rotation of R axis direction advance, therefore, the location parameter in each channel can with the variation of sweep time and scanning times and
Variation.Every time when scanning, it is referred to three-dimensional system of coordinate shown in Fig. 1 C, three-dimensional coordinate is configured for each channel, as each
The location parameter in channel.
In other examples, it is contemplated that when the X-ray detector rotation of the R axis direction shown in Fig. 1 C is advanced, X-ray detection
The change in location in each channel of device can be decomposed into straight line in the z-direction and advance and the angle rotation in x/y plane.Therefore, every time
When scanning, it is referred to x-y coordinate system, is that center channel configures x coordinate and y-coordinate is its location parameter, the position in other channels
Set parameter by：The location parameter of the relative position and central passage of the channel and central passage determines.For example, X-ray detection
When device contains 700 channels, the location parameter of central passage is set to 349.5, the location parameter in other channels is set to channel rope
Draw the location parameter for subtracting central passage.
As referring to figure 1E, X-ray emitter to checked object projection be parallel beam when, central passage can be with
Refer to：The channel of the most intermediate ray 31 of beam is received in X-ray detector.
The embodiment of the present application can receive the projection number of the detector output of above-mentioned imaging system in correction of movement artifact
According to, data for projection can also be transferred from the storage region of above-mentioned imaging system, other modes can also be taken to obtain data for projection,
The embodiment of the present application is without limitation.Specific which data for projection for obtaining checked object, can be according to required correction essence
At least one in the period of motion of degree, the reconstructed positions of reconstruction image and checked object determines.Such as：If checked object packet
Include at least partly heart, the data for projection in the available N number of heart motion cycle corresponding period, in turn, if correction essence
Spend more demanding, N can be configured to biggish numerical value, such as 10 or more integer, if correction accuracy require it is lower, can will
N is configured to lesser numerical value, such as less than 10 positive integer.
Due to, with the data of motion feature, easily leading in the data for projection of acquisition and occurring motion artifacts in reconstruction image,
Therefore, the embodiment of the present application is before the data for projection reconstruction image based on acquisition, for correction of movement artifact, needs to tested pair
Resulting data for projection is corrected when as movement, to reduce influence of these data to reconstruction image.
Furthermore, it is contemplated that in actual imaging scene, no matter checked object swallow, eye movement, limb motion etc
Autokinetic movement, or occur heartbeat, vascular pulsation, gastrointestinal peristalsis etc physiological movement when, be difficult to by straight
The change in location for connecing monitoring checked object obtains its motion feature.Therefore, the embodiment of the present application obtains the projection number of checked object
According to rear, by tracking the mass center of acquired data for projection, the movement that the mass center of checked object changes with sweep time is mapped out
Feature.
In certain examples, after checked object is scanned by imaging system every time, if getting the projection number of checked object
According to, can first track the mass center of data for projection, then according between the mass center of checked object and the mass center of data for projection mapping close
System maps resulting mass center is tracked, obtains the motion feature that the mass center of the checked object changes with flyback time.This Shen
Please embodiment, checked object can reflect the variation of the mass center relative scanning time of checked object in the mass center of each sweep time
Relationship, therefore, the motion feature can be indicated by checked object in the mass center of each sweep time.Wherein it is possible to mass center
Location parameter identifies mass center, such as center-of-mass coordinate.
In other embodiments, when can also be to each scanning of multiple sweep times, checked object in multiple sweep times
Between mass center, carry out data fitting, the variation relation of the mass center relative scanning time of checked object is obtained, to indicate checked object
The motion feature that changes over time of mass center.When carrying out data fitting here, the quantity of the sweep time of fitting can be according to institute
Need correction accuracy, reconstruction image the range of the corresponding scan position of reconstructed positions and in the period of motion of checked object extremely
One item missing determines.Such as：Required precision is higher, then the quantity for the sweep time being fitted is more.
In this example, the data of acquisition may include the multiple groups data for projection of checked object, and one group of data for projection is tested
The scanned primary resulting data of object, when tracking the mass center of the data for projection of acquisition, can track the matter of each group data for projection
The heart obtains every group of data for projection in the mass center of its corresponding sweep time.It, can be by each group data for projection in its correspondence when mapping
Sweep time mass center, be each mapped to checked object in the mass center of respective scanned time.Specific implementation can be refering to 2, Fig. 2
Shown method may comprise steps of S201-S204：
Step S201, the corresponding location parameter of each group data for projection and sweep time are determined.
Step S202, by every group of data for projection and its corresponding location parameter, this group of data for projection is converted in its correspondence
Sweep time mass center.
Step S203, by every group of data for projection corresponding sweep time and its in the mass center of the sweep time, be mapped as
Mass center of the checked object in each sweep time.
Step S204, the mass center to the checked object in each sweep time carries out data fitting, obtains described tested pair
The variation relation of the mass center relative scanning time of elephant.
In step s 201, location parameter and sweep time can be determined simultaneously, can also be separately determined.Location parameter and
Sweep time specifically may refer to the related content in the embodiment that Figure 1B to Fig. 1 E is related to.
In addition, location parameter can also be that other can indicate the parameter of the relative positional relationship between data for projection, this Shen
Please embodiment it is without limitation.
In step S202, can according to the predetermined corresponding relationship between mass center and data for projection and its corresponding position parameter,
By every group of data for projection and its corresponding location parameter, this group of data for projection is converted in the mass center of its corresponding sweep time.
In other examples, it is public every group of data for projection and its corresponding location parameter can also to be substituted into scheduled centroid calculation
Formula obtains this group of data for projection in the mass center of its corresponding sweep time.
When imaging system is CT imaging system, and the X-ray beam of transmitting is parallel beam, the centroid calculation formula can
With such as formula (1)：
Wherein, nChannelNum indicates the number in channel contained by detector, p (θi, tj) indicate and scanning angle θiWith
Location parameter tjCorresponding data for projection, pcom(θi) indicate that corresponding scanning angle is θiData for projection mass center, i ∈ [0,
nHalfViewPerRot）, j ∈ [0, nChannnelNum), nHalfViewPerRot indicate scanning support half rotation when
In, scanning times of the scanning device to checked object.
In step S203, according to the mapping relations of the mass center of the checked object and its projected centroids, every group is projected
Data corresponding sweep time and its mass center in the sweep time, are mapped as the checked object in corresponding sweep time
Mass center.
In other examples, every group of data for projection and its corresponding location parameter can also be substituted into scheduled mapping equation,
This group of data for projection is obtained in the mass center of its corresponding sweep time.
When imaging system is CT imaging system, and the X-ray beam of transmitting is parallel beam, the mapping equation can be as
Formula (2)：
pcom(θi)=xcom(t)*cosθi+ycom(t)*sinθi(2)；
Wherein, (xcom(t), ycom(t)) indicate checked object in the center-of-mass coordinate of t time.
By each group data for projection and its corresponding location parameter, after substituting into the mapping equation, obtain as shown in formula (3)
Result：
Wherein, t- Δ t < ti< t+ Δ t, i ∈ [0, n], Δ t > 0 and specific value can be by this programme designer's root
It is set according to actual scene situation, if imaging system is CT system, checked object includes heart, needs to keep away as far as possible when Δ t is arranged
The equalization for exempting from heart movement rule is set to 7s to the numerical value between 10s.
In the embodiment of the present application, formula (3) can also reflect the variation relation of the mass center relative scanning time of checked object, because
This, can indicate the motion feature of the mass center relative scanning time change of the checked object with formula (3).It can also be to above-mentioned public affairs
Formula (3) carries out data fitting, obtains the variation relation of the mass center relative scanning time of the checked object, to indicate described tested
The motion feature of the mass center relative scanning time change of object.And data fitting mode can there are many, apply embodiment pair
This is with no restrictions.
In certain examples, variation relation shown in following formula can be obtained by least square method：
X (t)=(A (t)TA(t))-1A(t)TY(t) (4)；
Wherein, A (t), Y (t), X (t) distinguish shown in following formula：
In other embodiments, in the mass center by the acquired data for projection of tracking, the mass center of checked object is mapped with sweeping
When retouching the motion feature of time change, it can also will directly be obtained in the data for projection after getting checked object and being scanned every time
The data for projection taken substitutes into formula (4) as above, obtains the variation relation of the mass center relative scanning time of the checked object.
For the variation relation that any of the above-described embodiment obtains, if held during the checked object is scanned
The equipment for carrying the checked object remains static, and can be closed with the variation of the mass center relative scanning time of the checked object
System, indicates motion feature of the mass center in each sweep time of the checked object.The carrying checked object mentioned herein
Equipment, can be different with the difference of imaging system, as imaging system be CT system when, which can be scanning bed.
If during the checked object is scanned, the equipment for carrying the checked object is kept in motion, needs
Correct movement deviation caused by centroid calculation result of the equipment.It, can be by checked object in each scanning in one example
Between mass center, subtract mass center average value of the sweep time nearby within the scope of certain time, correct error, specific amendment can lead to
Following operation is crossed to realize：
Determine mass center accumulated value of the checked object within the scope of the predetermined time before and after each sweep time.
The half of the ratio of the mass center accumulated value and the predetermined time range is calculated, is equipment error.
By the checked object in the mass center of the sweep time and the difference of the equipment error, the checked object is obtained
In the amendment mass center of the sweep time.Wherein, the predetermined time range can refer to that above-mentioned Δ t, the amendment mass center are to disappear
Except the mass center after equipment error.
This example can indicate the matter of the checked object with the checked object in the amendment mass center of each sweep time
Motion feature of the heart in each sweep time.It can also be quasi- in the amendment mass center progress data of each sweep time to the checked object
It closes, obtains the variation relation of the amendment mass center relative scanning time of the checked object, indicate that the mass center of the checked object exists
The motion feature of each sweep time.
In another example, the variation relation of the mass center relative scanning time for the checked object that data are fitted, such as
Formula (4) can correct movement deviation caused by centroid calculation result of the equipment by operating as follows：
The time difference of maximum scan time determined by calculating and minimum sweep time, maximum scan time here and most
Small sweep time can refer to the maximum value and minimum value of each data for projection corresponding sweep time.
The checked object is obtained in the accumulation result of the mass center of each sweep time.
It calculates accumulation result and the ratio of the time difference is equipment error.
The equipment error is eliminated from the variation relation, obtains revised variation relation.
In other examples, when removing movement deviation caused by centroid calculation result of the equipment, can also directly by by
The mass center that object is examined in each sweep time substitutes into following formula, obtains revised variation relation：
Wherein, X (ti) it is checked object in tiThe mass center at moment.
In practical application, it is contemplated that be more convenient for analyzing the motion feature of mass center as reference axis using main movement shaft, the application is real
Example is applied, variation relation resulting for any of the above-described embodiment can be changing coordinates in the main movement shaft of the checked object
Any reference axis of system indicates described tested pair with the variation relation of the amendment mass center relative scanning time of the checked object
Motion feature of the mass center of elephant in each sweep time, changing coordinates system mentioned herein can be by designer's roots of this programme
It is set according to practical application scene, the X-Y coordinate in three-dimensional system of coordinate as shown in Figure 1 C.
It, can be with base in an example when the main movement shaft of the checked object is not any reference axis of current coordinate system
In revised variation relation, the main movement shaft of the mass center of the checked object is determined, with the axis of abscissas of current coordinate system
Angle value；By revised variation relation, coordinate rotation is carried out according to the angle value, obtains postrotational variation relation.
In another example, revised variation relation can be substituted into following formula, obtain inertial matrix：
Wherein, Ixx、Ixy、IyyIt is as follows respectively：
Wherein, xucom(ti) lateral coordinates of the mass center of checked object, y in representation formula (8)ucom(ti) representation formula (8)
The longitudinal coordinate of the mass center of middle checked object.
The orthogonal available diagonal matrix being shown below of variation is carried out to formula (9)：
Wherein, I 'xx≥I′xy。
Carrying out the corresponding rotation angle of orthogonal transformation above is, by the transverse direction of the mass center of checked object in formula (8)
Coordinate and longitudinal coordinate are rotated according to the rotation angle, available main movement shaft institute in a coordinate system, the mass center of checked object
The variation relation of related scans time, as shown by the following formula：
As can be seen from the above embodiments, by the operation of any of the above embodiment, checked object is likely to be obtained in each scanning
Between mass center, checked object each sweep time amendment mass center, checked object the mass center relative scanning time variation relation,
Revised variation relation or postrotational variation relation can determine any expression according to practical application scene
The motion feature of the mass center relative scanning time of checked object.
After motion feature when determining the mass center relative scanning of checked object, it is contemplated that when scanning checked object, be detected
The motion intense degree of object is higher, and the motion artifacts in reconstruction image are more serious, and in an example, the application can be based on described
Motion feature determines motion intense degree of the mass center in each sweep time of the checked object, motion intense degree is higher than
Gained data for projection, is rejected from acquired data for projection in the time of predetermined value, and remaining data for projection is after correcting
Data for projection.Wherein, motion intense degree can be indicated by velocity magnitude.This predetermined value can be according to actual correction need
Setting is asked, the predetermined value is smaller if correction accuracy is higher.
In motion intense degree of the mass center in each sweep time of the checked object, if the motion feature is by quilt
The mass center or amendment mass center that object is examined in each sweep time indicate, can determine motion intense degree by following operation：
Calculate the difference of adjacent sweep time.
It calculates the mass center of adjacent twice sweep time or corrects the difference of mass center.
Then the ratio for calculating above-mentioned two difference obtains mass center or corrects mass center in the average speed of adjacent sweep time
Degree, indicates motion intense degree of the mass center of checked object between adjacent sweep time.
If the motion feature is by being fitted resulting variation relation, revised variation relation or postrotational variation
Relationship indicates, can directly calculate the derivative of the variation relation, obtain the systemic velocity relative scanning time of checked object
Changing rule, to indicate motion intense degree of the mass center in each sweep time of checked object.
In other embodiments, can also be in such a way that weight be rebuild in each data for projection configuration to obtain, correction tool
There is the data for projection for moving a feature, specifically may refer to Fig. 3, method shown in Fig. 3 may include steps of S301-
S303：
Step S301, based on the corresponding relationship and the motion feature between sweep time and scan position, determine described in
Motion intense degree of the mass center of checked object in each scan position.Scan position mentioned herein can be with scheduled coordinate
The difference of system and it is different, if imaging system be CT system, scan position can the Z coordinate as shown in Fig. 1 C indicate.
Step S302, the quilt is determined in the motion intense degree of each scan position according to the mass center of the checked object
Object is examined in the reconstruction weight of the corresponding data for projection of each scan position, wherein the more violent scan position of center of mass motion, it is corresponding
Data for projection reconstruction weight it is smaller.The specific corresponding relationship of the two can be set according to actual correction demand, such as at line
Sexual intercourse, Cos power, trapezoidal relationship etc..
Step S303, the corresponding reconstruction weight product of the data for projection of each scan position is obtained, after being corrected
Data for projection.
In one example, can by it is following operation based between sweep time and scan position corresponding relationship and the fortune
Dynamic feature, determines motion intense degree of the mass center in each scan position of the checked object, including：
According to the corresponding relationship between sweep time and scan position, the motion feature is converted into the checked object
The movement relation of relatively each scan position of mass center.
Based on the movement relation, motion intense degree of the mass center in each scan position of the checked object is calculated.
In this example, if the motion feature is by being fitted resulting variation relation, revised variation relation or rotation
Variation relation after turning indicates, the variation relation can be converted to relatively each scan position of mass center of the checked object
Then movement relation carries out derivation to the movement relation, obtains movement of the mass center in each scan position of the checked object
Speed indicates mass center in the motion intense degree of each scan position.
If mass center or amendment mass center of the motion feature by checked object in each sweep time indicate, can be based on sweeping
Retouch the corresponding relationship between time and scan position, by checked object each sweep time mass center or amendment mass center, be converted to by
Object is examined in the mass center of each scan position or amendment mass center, then calculates the difference of adjacent scan position；Calculate adjacent scan position
Mass center or correct mass center difference；Then the ratio for calculating above-mentioned two difference, obtains mass center or amendment mass center is swept adjacent
The average speed of position is retouched, indicates mass center in the motion intense degree of each scan position.
It, can be by following operation based on corresponding relationship between sweep time and scan position and described in another example
Motion feature determines motion intense degree of the mass center in each scan position of the checked object, including：
Based on the motion feature, motion intense degree of the mass center in each sweep time of the checked object is calculated.
Based on the corresponding relationship between sweep time and scan position, by the mass center of the checked object in each sweep time
Motion intense degree is converted to motion intense degree of the mass center in each scan position of the checked object.
In this example, if the motion feature is by being fitted resulting variation relation, revised variation relation or rotation
Variation relation after turning indicates, can obtain the checked object in the movement of each sweep time to the variation relation derivation
Speed；It is converted to according to the corresponding relationship of sweep time and scan position by checked object in the movement velocity of each sweep time
The mass center of checked object indicates motion intense of the mass center in each scan position of checked object in the movement velocity of each scan position
Degree.
If mass center or amendment mass center of the motion feature by checked object in each sweep time indicate, phase can be calculated
The difference of adjacent sweep time；It calculates the mass center of adjacent twice sweep time or corrects the difference of mass center；Then it calculates above-mentioned two
The ratio of difference obtains mass center or corrects average speed of the mass center in adjacent sweep time, then based on sweep time and scanning position
Corresponding relationship between setting, the average speed by mass center or amendment mass center in adjacent sweep time, is converted to the mass center of checked object
Or average speed of the amendment mass center between adjacent scan position, indicate that the mass center of checked object is acute in the movement of each scan position
Strong degree.
As can be seen from the above embodiments, based on the original for reducing the influence for having the data for projection of motion feature to reconstruction image
Reason after correcting the data for projection with motion feature, can rebuild the image of checked object, weight based on the data for projection after correction
The process of building can be different with imaging system, the difference of method for reconstructing, specifically may refer to the image weight of corresponding imaging system
Technology is built, below only by taking CT system as an example, sketches lower image reconstruction process：
The x-ray attenuation coefficient or absorption coefficient of each voxel are obtained based on the data for projection after correction.By X-ray attenuation
Coefficient is arranged in character matrix (digital matrix), and each number in the character matrix is represented by each pixel
(pixel) the image value such as gray scale.CT image is constituted by all pixels that the character matrix generates.
It is illustrative to introduce lower the application correction of movement artifact below in conjunction with attached drawing 4A to 4C and specific application scenarios
Process：
The imaging system that this application scene uses scans the lung of human body for CT system, since lung is close to heart, institute
It may include lung, heart and other tissues with checked object.Pass through tested pair of the scanning of scanning device shown in Fig. 1 C
As.
When scanning support 11 surrounds the axis R rotation traveling for being parallel to the direction z, X-ray emitter is repeatedly thrown to checked object
Parallel beam is penetrated, X-ray detector 16 receives one group of X-ray every time, and this group of X-ray is changed into telecommunications through photoelectric conversion
Number, then one group of data for projection is converted to through ADC (Analog/Digital Converter, analog/digital converter), input CT
The image reconstruction device of system.
If image reconstruction device directly base each group data for projection is carried out image reconstruction and is likely to be obtained due to heartbeat
The reconstruction image containing motion artifacts for being such as Figure 1A.
For correction of movement artifact, image reconstruction device can project the mass center of number by tracking each group, map out described
The motion feature that the mass center of checked object changes with sweep time specifically can obtain the movement by formula (1) to (14)
Feature.
Under this scene, if human body does not move when scanning, obtained motion feature can reflect the motion feature of heart, such as
Shown in Fig. 4 A, abscissa is sweep time, and ordinate is coordinate of the mass center in main movement shaft, and curve indicates the matter of checked object
The motion feature of heart relative scanning time is marked at the time that the derivative of ascent direction is 0 in figure by ECG signal R blob detection
For Rn, this process is equivalent to the R blob detection process of ECG signal.
In order to reduce the motion artifacts of reconstruction image, need to correct acquired data for projection based on the motion feature
In with motion feature data for projection；Image reconstruction is carried out further according to the data for projection after correction, obtains the checked object
Reconstruction image.
Wherein, when corrected projection data, can according to the corresponding relationship of motion feature and sweep time and scan position,
The movement relation of the mass center versus scan location of checked object is obtained, as shown in Figure 4 B, abscissa is scan position, and ordinate is
Coordinate of the mass center in main movement shaft, curve indicate the motion feature of the mass center versus scan location of checked object.Wherein, z0-z
The corresponding scan position range of data for projection needed for needing to rebuild the image at z0 to z0+z expression.
The movement relation of mass center versus scan location based on checked object can be respectively sweeping within the scope of z0-z to z0+z
Retouch the corresponding data for projection in position, weight is rebuild in configuration, when configuration, can be configured according to corresponding relationship shown in Fig. 4 C, Fig. 4 C
In, phaseDIndicate that the relative phase positions at the peak R, abscissa are scan position, curve 1 is the mass center relative scanning of checked object
The movement relation of position, curve 2 are the reconstruction weight of data for projection and the corresponding relationship of scan position.
It is configured to the corresponding data for projection of each scan position within the scope of z0-z to z0+z according to relationship shown in Fig. 4 C, is matched
It sets reconstruction weight, and data for projection is rebuild into the product of weight as the data for projection after correction with it, after carrying out image reconstruction,
Motion artifacts in reconstruction image can achieve preferable improvement.
Corresponding with the example of method of the application correction of movement artifact, present invention also provides the dresses of correction of movement artifact
The example set.The device of correction of movement artifact can be applied to various computer equipments, such as the imaging of CT, PET, PET-CT etc
In system, carry out the image reconstruction device of image procossing and/or image reconstruction, be also possible to imaging system it is outer other with data
The equipment of processing capacity.As shown in figure 5, the hardware configuration of the computer equipment for the device application of the application correction of movement artifact
Schematic diagram, the computer equipment may include processor 510, memory 520, nonvolatile memory 530.Wherein, 520 He of memory
Nonvolatile memory 530 is machine readable storage medium, and processor 510 and machine readable storage medium 520,530 can be by
Internal bus 540 is connected with each other.In other possible implementations, the computer equipment is also possible that network interface
550, can be communicated with other equipment or component.In addition to processor 510 shown in fig. 5, memory 520, network interface
550 and nonvolatile memory 530 except, which needs to include other hardware according to actual functional capability, in Fig. 5
No longer show one by one.
In different examples, the machine readable storage medium 520,530 can be ROM (Read-Only Memory,
Read-only memory), volatile memory, nonvolatile memory, flash memory, memory driver (such as hard disk drive), solid state hard disk,
Any kind of storage dish (such as CD, DVD) perhaps similar storage medium or their combination.
Further, machine readable storage medium, can be specially the device being stored on memory 520 with correction of movement artifact
600 corresponding machine-executable instructions.It functionally divides, as shown in fig. 6, the device 600 of correction of movement artifact may include number
According to acquisition module 610, characteristic extracting module 620, Data correction module 630 and image reconstruction module 640.
Wherein, data acquisition module 610, for obtaining the data for projection of checked object.
Characteristic extracting module 620 maps out the checked object for the mass center by the acquired data for projection of tracking
The motion feature that changes with sweep time of mass center.
Data correction module 630, for being based on the motion feature, correcting in acquired data for projection has movement special
The data for projection of sign, wherein include with the data for projection of motion feature：The checked object is scanned in its center of mass motion
Resulting data.
Image reconstruction module 640 obtains the checked object for carrying out image reconstruction according to the data for projection after correction
Reconstruction image.
In some examples, acquired data for projection includes multiple groups data for projection, and one group of data for projection is checked object quilt
The resulting data of run-down, characteristic extracting module 620 may include：
Parameter determination module, for determining the corresponding location parameter of each group data for projection and sweep time.
Centroid calculation module, for being converted to this group of data for projection for every group of data for projection and its corresponding location parameter
In the mass center of its corresponding sweep time.
Mass center mapping block, for by every group of data for projection corresponding sweep time and its sweep time matter
The heart is mapped as the checked object in the mass center of corresponding sweep time.
First fitting module carries out data fitting for the mass center to the checked object in each sweep time, obtains institute
State the variation relation of the mass center relative scanning time of checked object.
As an example, the location parameter is the location parameter in the channel of acquired projections data.
As an example, characteristic extracting module 620 can also include：
Characteristic determination module, for during the checked object is scanned, the equipment for carrying the checked object to be in
When stationary state, with the variation relation of the mass center relative scanning time of the checked object, the mass center of the checked object is indicated
In the motion feature of each sweep time.
Time difference determining module, at the checked object scanned period, the equipment for carrying the checked object
When motion state, the time difference of identified maximum scan time and minimum sweep time are calculated.
Mass center accumulator module, for obtaining the checked object in the accumulation result of the mass center of each sweep time.
Error calculating module, the ratio for calculating accumulation result and the time difference are equipment error.
Error correction module obtains revised variation and closes for eliminating the equipment error from the variation relation
System.
As an example, characteristic extracting module 620 can also include：
Angle determining module determines the main motion of the mass center of the checked object for being based on revised variation relation
Axis, the angle value with the axis of abscissas of current coordinate system.
Coordinate rotary module, for carrying out revised variation relation coordinate rotation according to the angle value, being revolved
Variation relation after turning.
In other examples, Data correction module 630 may include：
Movement degree determining module, for based between sweep time and scan position corresponding relationship and the movement it is special
Sign, determines motion intense degree of the mass center in each scan position of the checked object.
Rebuild weight determination module, for the mass center according to the checked object each scan position motion intense journey
Degree, determines the checked object in the reconstruction weight of the corresponding data for projection of each scan position, wherein center of mass motion is more violent
The reconstruction weight of scan position, corresponding data for projection is smaller.
Data correction submodule, the corresponding reconstruction weight product of data for projection for obtaining each scan position, obtains
Data for projection after to correction.
As an example, the movement degree determining module is configured as：
According to the corresponding relationship between sweep time and scan position, the motion feature is converted into the checked object
The movement relation of relatively each scan position of mass center.
Based on the movement relation, motion intense degree of the mass center in each scan position of the checked object is calculated.
As an example, the movement degree determining module is configured as：
Based on the motion feature, motion intense degree of the mass center in each sweep time of the checked object is calculated.
Based on the corresponding relationship between sweep time and scan position, by the mass center of the checked object in each sweep time
Motion intense degree is converted to motion intense degree of the mass center in each scan position of the checked object.
Apparatus above embodiment is corresponding to embodiment of the method, and details are not described herein.
Those skilled in the art will readily occur to its of the application after considering specification and practicing this disclosure
Its embodiment.This application is intended to cover any variations, uses, or adaptations of the application, these modifications, purposes or
Person's adaptive change follows the general principle of the application and including the undocumented common knowledge in the art of the application
Or conventional techniques.The description and examples are only to be considered as illustrative, and the true scope and spirit of the application are by following
Claim is pointed out.
It should be understood that the application is not limited to the precise structure that has been described above and shown in the drawings, and
And various modifications and changes may be made without departing from the scope thereof.Scope of the present application is only limited by the accompanying claims.
Claims (18)
1. a kind of method of correction of movement artifact, which is characterized in that the method includes：
Obtain the data for projection of checked object；
By tracking the mass center of acquired data for projection, the fortune that the mass center of the checked object changes with sweep time is mapped out
Dynamic feature；
Based on the motion feature, the data for projection in acquired data for projection with motion feature is corrected, wherein there is fortune
The data for projection of dynamic feature includes：The checked object is scanned resulting data in its center of mass motion；
Image reconstruction is carried out according to the data for projection after correction, obtains the reconstruction image of the checked object.
2. the method according to claim 1, wherein acquired data for projection includes multiple groups data for projection, one
Group data for projection is the scanned primary resulting data of checked object, the mass center by tracking acquired data for projection,
The motion feature that the mass center of the checked object changes with sweep time is mapped out, including：
Determine the corresponding location parameter of each group data for projection and sweep time；
By every group of data for projection and its corresponding location parameter, this group of data for projection is converted in the matter of its corresponding sweep time
The heart；
By every group of data for projection corresponding sweep time and its mass center in the sweep time, it is mapped as the checked object and exists
The mass center of corresponding sweep time；
Mass center to the checked object in each sweep time carries out data fitting, and the mass center for obtaining the checked object is swept relatively
Retouch the variation relation of time.
3. according to the method described in claim 2, it is characterized in that, the location parameter is the position in the channel of acquired projections data
Set parameter.
4. according to the method described in claim 2, it is characterized in that, the mass center by tracking acquired data for projection,
The motion feature that the mass center of the checked object changes with sweep time is mapped out, further includes：
If the equipment for carrying the checked object is kept in motion during the checked object is scanned, then institute is calculated
The time difference of determining maximum scan time and minimum sweep time；
The checked object is obtained in the accumulation result of the mass center of each sweep time；
It calculates accumulation result and the ratio of the time difference is equipment error；
The equipment error is eliminated from the variation relation, obtains revised variation relation.
5. according to the method described in claim 4, it is characterized in that, the mass center by tracking acquired data for projection,
The motion feature that the mass center of the checked object changes with sweep time is mapped out, further includes：
Based on revised variation relation, the main movement shaft of the mass center of the checked object is determined, the horizontal seat with current coordinate system
The angle value of parameter；
By revised variation relation, coordinate rotation is carried out according to the angle value, obtains postrotational variation relation.
6. the method according to any one of claims 1 to 5, which is characterized in that described to be based on the motion feature, correction
With the data for projection of motion feature in acquired data for projection, including：
Based on the corresponding relationship and the motion feature between sweep time and scan position, the mass center of the checked object is determined
In the motion intense degree of each scan position；
According to the mass center of the checked object in the motion intense degree of each scan position, determine the checked object in each scanning
The reconstruction weight of the corresponding data for projection in position, wherein the more violent scan position of center of mass motion, the weight of corresponding data for projection
It is smaller to build weight；
Obtain the corresponding reconstruction weight product of the data for projection of each scan position, the data for projection after being corrected.
7. according to the method described in claim 6, it is characterized in that, the corresponding pass based between sweep time and scan position
System and the motion feature, determine motion intense degree of the mass center in each scan position of the checked object, including：
According to the corresponding relationship between sweep time and scan position, the motion feature is converted to the mass center of the checked object
The movement relation of relatively each scan position；
Based on the movement relation, motion intense degree of the mass center in each scan position of the checked object is calculated.
8. according to the method described in claim 6, it is characterized in that, the corresponding pass based between sweep time and scan position
System and the motion feature, determine motion intense degree of the mass center in each scan position of the checked object, including：
Based on the motion feature, motion intense degree of the mass center in each sweep time of the checked object is calculated；
Based on the corresponding relationship between sweep time and scan position, by the mass center of the checked object each sweep time movement
Severe degree is converted to motion intense degree of the mass center in each scan position of the checked object.
9. a kind of device of correction of movement artifact, which is characterized in that including：
Data acquisition module, for obtaining the data for projection of checked object；
Characteristic extracting module maps out the mass center of the checked object for the mass center by the acquired data for projection of tracking
The motion feature changed with sweep time；
Data correction module corrects the throwing in acquired data for projection with motion feature for being based on the motion feature
Shadow data, wherein include with the data for projection of motion feature：The checked object is scanned resulting in its center of mass motion
Data；
Image reconstruction module obtains the reconstruction of the checked object for carrying out image reconstruction according to the data for projection after correction
Image.
10. device according to claim 9, which is characterized in that acquired data for projection includes multiple groups data for projection, and one
Group data for projection is the scanned primary resulting data of checked object, and the characteristic extracting module includes：
Parameter determination module, for determining the corresponding location parameter of each group data for projection and sweep time；
Centroid calculation module, for being converted to this group of data for projection at it for every group of data for projection and its corresponding location parameter
The mass center of corresponding sweep time；
Mass center mapping block, for by every group of data for projection corresponding sweep time and its in the mass center of the sweep time, reflect
Penetrate the mass center for the checked object in corresponding sweep time；
First fitting module carries out data fitting for the mass center to the checked object in each sweep time, obtains the quilt
Examine the variation relation of the mass center relative scanning time of object.
11. device according to claim 10, which is characterized in that the location parameter is the channel of acquired projections data
Location parameter.
12. device according to claim 10, which is characterized in that the characteristic extracting module further includes：
Time difference determining module, the equipment for during the checked object is scanned, carrying the checked object are in fortune
When dynamic state, the time difference of identified maximum scan time and minimum sweep time are calculated；
Mass center accumulator module, for obtaining the checked object in the accumulation result of the mass center of each sweep time；
Error calculating module, the ratio for calculating accumulation result and the time difference are equipment error；
Error correction module obtains revised variation relation for eliminating the equipment error from the variation relation.
13. device according to claim 12, which is characterized in that the characteristic extracting module further includes：
Angle determining module, for determining the main movement shaft of the mass center of the checked object based on revised variation relation, with
The angle value of the axis of abscissas of current coordinate system；
Coordinate rotary module, for coordinate rotation being carried out according to the angle value, after obtaining rotation for revised variation relation
Variation relation.
14. the device according to any one of claim 9 to 13, which is characterized in that the Data correction module includes：
Movement degree determining module, for based on the corresponding relationship and the motion feature between sweep time and scan position,
Determine motion intense degree of the mass center in each scan position of the checked object；
Rebuild weight determination module, for the mass center according to the checked object each scan position motion intense degree, really
Reconstruction weight of the fixed checked object in the corresponding data for projection of each scan position, wherein the more violent scanning of center of mass motion
The reconstruction weight of position, corresponding data for projection is smaller；
Data correction submodule, the corresponding reconstruction weight product of data for projection for obtaining each scan position, obtains school
Data for projection after just.
15. device according to claim 14, which is characterized in that the movement degree determining module is configured as：
According to the corresponding relationship between sweep time and scan position, the motion feature is converted to the mass center of the checked object
The movement relation of relatively each scan position；
Based on the movement relation, motion intense degree of the mass center in each scan position of the checked object is calculated.
16. device according to claim 14, which is characterized in that the movement degree determining module is configured as：
Based on the motion feature, motion intense degree of the mass center in each sweep time of the checked object is calculated；
Based on the corresponding relationship between sweep time and scan position, by the mass center of the checked object each sweep time movement
Severe degree is converted to motion intense degree of the mass center in each scan position of the checked object.
17. a kind of computer equipment, which is characterized in that including：
Processor；
The memory of storage processor executable instruction；
Wherein, the processor is coupled in the memory, for reading the program instruction of the memory storage, and as sound
It answers, executes such as the operation in any one of claim 1-8 the method.
18. one or more machine readable storage mediums, which is characterized in that instruction is stored thereon with, when by one or more
When managing device execution, so that computer equipment is executed such as the operation in any one of claim 1-8 the method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810509389.0A CN108876730B (en) | 2018-05-24 | 2018-05-24 | Method, device and equipment for correcting motion artifact and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810509389.0A CN108876730B (en) | 2018-05-24 | 2018-05-24 | Method, device and equipment for correcting motion artifact and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108876730A true CN108876730A (en) | 2018-11-23 |
| CN108876730B CN108876730B (en) | 2022-03-04 |
Family
ID=64333818
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810509389.0A Active CN108876730B (en) | 2018-05-24 | 2018-05-24 | Method, device and equipment for correcting motion artifact and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108876730B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110443866A (en) * | 2019-07-31 | 2019-11-12 | 东软医疗系统股份有限公司 | Image rebuilding method, device, terminal device and PET system |
| CN110473269A (en) * | 2019-08-08 | 2019-11-19 | 上海联影医疗科技有限公司 | A kind of image rebuilding method, system, equipment and storage medium |
| CN111462168A (en) * | 2020-04-22 | 2020-07-28 | 上海联影医疗科技有限公司 | Motion parameter estimation method and motion artifact correction method |
| CN111528825A (en) * | 2020-05-14 | 2020-08-14 | 浙江大学 | Photoelectric volume pulse wave signal optimization method |
| CN111544020A (en) * | 2020-04-17 | 2020-08-18 | 北京东软医疗设备有限公司 | Geometric correction method and device for X-ray imaging equipment |
| CN112001979A (en) * | 2020-07-31 | 2020-11-27 | 上海联影医疗科技股份有限公司 | Motion artifact processing method, system, readable storage medium and device |
| WO2021170147A1 (en) * | 2020-02-28 | 2021-09-02 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for correcting motion artifacts in images |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101454801A (en) * | 2006-02-28 | 2009-06-10 | 皇家飞利浦电子股份有限公司 | Local motion compensation based on list mode data |
| CN101936720A (en) * | 2010-07-30 | 2011-01-05 | 北京航空航天大学 | Method for calibrating detector torsion angle applied to cone-beam XCT system |
| CN202049120U (en) * | 2011-03-04 | 2011-11-23 | 首都师范大学 | System for eliminating geometric artifacts in CT (computed tomography) image |
| CN102317975A (en) * | 2009-02-17 | 2012-01-11 | 皇家飞利浦电子股份有限公司 | Functional imaging |
| CN102426696A (en) * | 2011-10-24 | 2012-04-25 | 西安电子科技大学 | Optical projection tomography motion artifact correction method |
| CN102652674A (en) * | 2011-03-04 | 2012-09-05 | 首都师范大学 | Method and system for eliminating geometrical artifacts in CT (Computerized Tomography) image |
| US20140126796A1 (en) * | 2012-11-08 | 2014-05-08 | Siemens Corporation | Mri reconstruction with motion-dependent regularization |
| CN203776924U (en) * | 2014-03-06 | 2014-08-20 | 北京锐视康科技发展有限公司 | Calibration device for geometric position of cone-beam CT (computed tomography) system |
| CN106251380A (en) * | 2016-07-29 | 2016-12-21 | 上海联影医疗科技有限公司 | Image rebuilding method |
| US20170091962A1 (en) * | 2015-09-24 | 2017-03-30 | General Electric Company | Radiation Tomographic Imaging Method, Apparatus, and Program |
| CN107303184A (en) * | 2016-04-22 | 2017-10-31 | 上海联影医疗科技有限公司 | CT scan x-ray source tube current modulator approach and computed tomographic apparatus |
-
2018
- 2018-05-24 CN CN201810509389.0A patent/CN108876730B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101454801A (en) * | 2006-02-28 | 2009-06-10 | 皇家飞利浦电子股份有限公司 | Local motion compensation based on list mode data |
| CN102317975A (en) * | 2009-02-17 | 2012-01-11 | 皇家飞利浦电子股份有限公司 | Functional imaging |
| CN101936720A (en) * | 2010-07-30 | 2011-01-05 | 北京航空航天大学 | Method for calibrating detector torsion angle applied to cone-beam XCT system |
| CN202049120U (en) * | 2011-03-04 | 2011-11-23 | 首都师范大学 | System for eliminating geometric artifacts in CT (computed tomography) image |
| CN102652674A (en) * | 2011-03-04 | 2012-09-05 | 首都师范大学 | Method and system for eliminating geometrical artifacts in CT (Computerized Tomography) image |
| CN102426696A (en) * | 2011-10-24 | 2012-04-25 | 西安电子科技大学 | Optical projection tomography motion artifact correction method |
| US20140126796A1 (en) * | 2012-11-08 | 2014-05-08 | Siemens Corporation | Mri reconstruction with motion-dependent regularization |
| CN203776924U (en) * | 2014-03-06 | 2014-08-20 | 北京锐视康科技发展有限公司 | Calibration device for geometric position of cone-beam CT (computed tomography) system |
| US20170091962A1 (en) * | 2015-09-24 | 2017-03-30 | General Electric Company | Radiation Tomographic Imaging Method, Apparatus, and Program |
| CN107303184A (en) * | 2016-04-22 | 2017-10-31 | 上海联影医疗科技有限公司 | CT scan x-ray source tube current modulator approach and computed tomographic apparatus |
| CN106251380A (en) * | 2016-07-29 | 2016-12-21 | 上海联影医疗科技有限公司 | Image rebuilding method |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110443866A (en) * | 2019-07-31 | 2019-11-12 | 东软医疗系统股份有限公司 | Image rebuilding method, device, terminal device and PET system |
| CN110473269A (en) * | 2019-08-08 | 2019-11-19 | 上海联影医疗科技有限公司 | A kind of image rebuilding method, system, equipment and storage medium |
| WO2021170147A1 (en) * | 2020-02-28 | 2021-09-02 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for correcting motion artifacts in images |
| CN111544020A (en) * | 2020-04-17 | 2020-08-18 | 北京东软医疗设备有限公司 | Geometric correction method and device for X-ray imaging equipment |
| CN111544020B (en) * | 2020-04-17 | 2023-08-01 | 北京东软医疗设备有限公司 | Geometric correction method and device for X-ray imaging equipment |
| CN111462168A (en) * | 2020-04-22 | 2020-07-28 | 上海联影医疗科技有限公司 | Motion parameter estimation method and motion artifact correction method |
| CN111462168B (en) * | 2020-04-22 | 2023-09-19 | 上海联影医疗科技股份有限公司 | Motion parameter estimation method and motion artifact correction method |
| CN111528825A (en) * | 2020-05-14 | 2020-08-14 | 浙江大学 | Photoelectric volume pulse wave signal optimization method |
| CN112001979A (en) * | 2020-07-31 | 2020-11-27 | 上海联影医疗科技股份有限公司 | Motion artifact processing method, system, readable storage medium and device |
| CN112001979B (en) * | 2020-07-31 | 2024-04-26 | 上海联影医疗科技股份有限公司 | Motion artifact processing method, system, readable storage medium and apparatus |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108876730B (en) | 2022-03-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108876730A (en) | The method, device and equipment and storage medium of correction of movement artifact | |
| CN110533597B (en) | Artifact processing method, artifact processing device, rotation center determining device, artifact processing equipment and storage medium | |
| US7558439B2 (en) | Motion artifact correction of tomographical images | |
| US7260252B2 (en) | X-ray computed tomographic apparatus, image processing apparatus, and image processing method | |
| CN110751702B (en) | Image reconstruction method, system, device and storage medium | |
| CN102125439B (en) | Image processing apparatus, X-ray computed tomography apparatus, and image processing method | |
| US20060285632A1 (en) | Method and device for reconstructing a 3D image data set of a moving object | |
| US20170079608A1 (en) | Spatial Registration Of Positron Emission Tomography and Computed Tomography Acquired During Respiration | |
| JP2007089674A (en) | Shape of appearance measuring apparatus and x-ray ct apparatus | |
| KR20060135560A (en) | X-ray ct apparatus | |
| US6381298B2 (en) | Method of combining reconstruction images | |
| CN107928694A (en) | CT dose modulations, device, CT scan method and CT system | |
| CN110866959B (en) | Image reconstruction method, system, device and storage medium | |
| CN105615912A (en) | CT (computed tomography) scanning method and system | |
| CN104605882B (en) | Image acquiring method, device and radiotherapy system in radiotherapy system | |
| CN107945850A (en) | Method and apparatus for handling medical image | |
| CN110473269A (en) | A kind of image rebuilding method, system, equipment and storage medium | |
| JP2020103902A (en) | Medical image diagnosis system and generation method for learned model | |
| CN106780651B (en) | Cardiac spiral CT image reconstruction method and device | |
| CN114511497A (en) | Imaging method and device applied to cone beam CT sparse sampling | |
| KR20050042727A (en) | Ct image producing method and x-ray ct apparatus | |
| US20210407159A1 (en) | Device and method for editing a panoramic radiography image | |
| JP5317612B2 (en) | Tomographic image processing apparatus, X-ray CT apparatus, and program | |
| CN102013116B (en) | Cerebral rotational angiography-based three-dimensional reconstruction method | |
| US20220414832A1 (en) | X-ray imaging restoration using deep learning algorithms |
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 |
Address after: 110167 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Applicant after: DongSoft Medical System Co., Ltd. Address before: 110167 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Applicant before: Dongruan Medical Systems Co., Ltd., Shenyang |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |