CN102805634B - Method and device for judging motion of object to be inspected - Google Patents
Method and device for judging motion of object to be inspected Download PDFInfo
- Publication number
- CN102805634B CN102805634B CN201110145206.XA CN201110145206A CN102805634B CN 102805634 B CN102805634 B CN 102805634B CN 201110145206 A CN201110145206 A CN 201110145206A CN 102805634 B CN102805634 B CN 102805634B
- Authority
- CN
- China
- Prior art keywords
- angle
- projection
- probe access
- data group
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The invention relates to the field of medical apparatuses, in particular to a method for judging motion of an object to be inspected. The method includes: acquiring a compensation data set for an original projection data set; calculating a distance matrix for the original projection data set and the compensation data set; and selecting a certain column in the distance matrix to locate a row value corresponding to a minimum in the column, and determining the existence of motion artifact in the original projection data set if an absolute of a difference between the row value and a column value is larger than a preset threshold. The invention further provides a device for judging motion of the object to be inspected. The motion of the object to be inspected in current scan can be judged directly from the original projection data set (0, 180 DEG + 2beta) by the method and the device, the specific position of the motion is obtained, and further dosage of X-ray taken by the object to be inspected is reduced.
Description
Technical field
The present invention relates to field of medical device, particularly a kind ofly judge the method and apparatus that object to be checked moves.
Background technology
Along with the progress of medical level, X ray computer fault imaging (Computed Tomography, CT) equipment is widely used.When carrying out CT scan to object to be checked, rotary frame is around Z axis rotating 360 degrees, and drive the data collection system synchronizing comprising x-ray source and detector around object rotating 360 degrees to be checked simultaneously, with from [0,360 °] projection angle obtain X-ray projection data, wherein Z axis be CT examination bed turnover rotary frame horizontal direction.At each projection angle place, M passage of detector receives X-ray projection data simultaneously, and the image re-construction system being sent to CT is for generating reconstruction image.But inevitably move due to human body, so when carrying out CT scan to human body, rebuilding image can exist motion artifacts because of the motion of human body, thus the quality of image is rebuild in impact.
Publication number is that the application for a patent for invention of 102018524A discloses a kind of artifact detection method and device, it is by calculating [0,360 °] each data for projection judges in reconstructed image with or without artifact with the difference of the data for projection of corresponding same channels in the group corresponding to last projection angle or a rear projection angle respectively in group in scope corresponding to each projection angle, if the arbitrary difference in gained difference is more than or equal to the upper threshold value that pre-sets or is less than or equal to the bottom threshold value pre-set, then judge to there is artifact in reconstructed image.But this application for a patent for invention is applicable to check whether CT scanner has fault, to avoid the error of hand inspection.
Summary of the invention
The object of the present invention is to provide a kind of method and apparatus whether object to be checked moves that judges, it needs data collecting system to carry out [0 around object to be checked, 180 ° of+2 β] project and can judge whether object to be checked moves, to reduce the radiation dose of object reception to be checked.
In view of this, the present invention proposes a kind of method that object to be checked moves that judges, described method comprises: the offset data group obtaining described Raw projection data group; Calculate the distance matrix of described Raw projection data group and described offset data group; Get certain string of described distance matrix, find out the row value in these row corresponding to minima, if the absolute value of the difference of this row value and this train value is greater than a predetermined threshold value, then judge that object to be checked moves.Visible the present invention directly can just judge from Raw projection data group whether object to be checked moves in present scan, accuracy rate and reliability higher.
According to one embodiment of present invention, described distance matrix is Euclidean distance matrix.
According to still a further embodiment, described acquisition offset data group comprises: read the first projection angle corresponding with a data for projection in described Raw projection data group and the first probe access; Calculate one second projection angle and the second probe access according to described first projection angle and described first probe access, the first probe access under the second probe access under wherein said second projection angle and described first projection angle is positioned on same x-ray spread path; The data for projection collected by second probe access under described second projection angle is as the data in described offset data group.
According to still another embodiment of the invention, described calculating second projection angle and the second probe access comprise: the angle calculating described x-ray spread path deviation X-ray beam center according to described first probe access by following formula: θ=β-n1* Δ β, wherein, n1 is described first probe access, Δ β is the angle step pitch of described probe access, θ is the angle at described x-ray spread path and X-ray beam center, and β is a half angle of described X-ray beam; Described second probe access is calculated by following formula: n2=(β+θ)/Δ β according to described angle, wherein, n2 is described second probe access, and described second probe access and described first probe access are specular about the center of described detector; Calculate described second projection angle according to described angle by following formula: α 2=α 1 ± pi-2 θ, wherein, α 2 is described second projection angle, and α 1 is described first projection angle.
According to still another embodiment of the invention, the scope of described first projection angle is [0,180 ° of+2 β].Such data collecting system only need carry out [0 around object to be checked, 180 ° of+2 β] projection can obtain described Raw projection data, the projection pattern comparing 360 ° decreases the radiation dose of object reception to be checked, decrease the data volume of Raw projection data group simultaneously, decrease data space.
According to still another embodiment of the invention, described predetermined threshold value is 20.
Further, described method comprises: if the absolute value of the difference of this row value and this train value is greater than described predetermined threshold value, then obtain according to this row value the particular location that described object to be checked moves.
The present invention also provides a kind of computer program, the method making computer perform described judgement object to be checked to move.
The present invention also provides a kind of computer-readable recording medium, stores described computer program.
The present invention also provides a kind of device that object to be checked moves that judges, described device comprises: an offset data securing component, for obtaining the offset data group of described Raw projection data group; One distance matrix computation module, for calculating the distance matrix of described Raw projection data group and described offset data group; One motion determination assembly, gets certain string of described distance matrix, finds out the row value in these row corresponding to minima, if the absolute value of the difference of this row value and this train value is greater than a predetermined threshold value, then judges that object to be checked moves.Visible the present invention directly can just judge from Raw projection data group whether object to be checked moves in present scan, thus accuracy rate and reliability higher.
According to still a further embodiment, described offset data securing component comprises: a read module, for reading the first projection angle corresponding with a data for projection in described Raw projection data group and the first probe access; One computing module, for calculating one second projection angle and one second probe access according to described first projection angle and described first probe access, the first probe access under the second probe access under wherein said second projection angle and described first projection angle is positioned on same x-ray spread path; One offset data acquisition module, its data for projection collected by the second probe access under described second projection angle is as the data in described offset data group.
According to still another embodiment of the invention, described computing module comprises: an angle calculation unit, for calculating the angle at described x-ray spread path deviation X-ray beam center by following formula according to described first probe access: θ=β-n1* Δ β, wherein, n1 is described first probe access, Δ β is the angle step pitch of described probe access, and θ is the angle at described x-ray spread path and X-ray beam center, and β is a half angle of described X-ray beam; One second channel computing unit, for calculating described second probe access according to described angle by following formula: n2=(β+θ)/Δ β, wherein, n2 is described second probe access, and described second probe access and described first probe access are specular about the center of described detector; One second angle calculation unit, for calculating described second projection angle according to described angle by following formula: α 2=α 1 ± pi-2 θ, wherein, α 2 is described second projection angle, and α 1 is described first projection angle.
Further, described device comprises: a movement position detection components, when the absolute value for the difference in this row value and this train value is greater than described predetermined threshold value, obtains according to this row value the particular location that described object to be checked moves.
As can be seen from such scheme, because the present invention directly can just judge from Raw projection data group whether object to be checked moves in present scan, thus accuracy rate and reliability higher; And described Raw projection data is that data collecting system carries out [0 around object to be checked, 180 ° of+2 β] projection, the projection pattern comparing 360 ° decreases the radiation dose of object reception to be checked, decreases the data volume of Raw projection data group simultaneously, decreases data space.In addition, the present invention can also obtain the particular location that object to be checked moves.
Accompanying drawing explanation
Embodiments of the invention will be described in detail by referring to accompanying drawing below, make clearer above-mentioned and other feature and advantage of the present invention of those skilled in the art, in accompanying drawing:
Fig. 1 is that the present invention judges the method schematic diagram that object to be checked moves.
Fig. 2 is the principle schematic that the present invention obtains offset data.
Fig. 3 a is object to be checked schematic diagram of certain column data in distance matrix of the present invention when not moving.
Fig. 3 b is object to be checked schematic diagram of certain column data in distance matrix of the present invention when moving.
Fig. 4 is that the present invention judges the device schematic diagram that object to be checked moves.
Fig. 5 is the schematic diagram of offset data securing component of the present invention.
Reference numeral
Step 101 ~ 108 judge the step that object to be checked moves
1X radiographic source
3 detectors
10 the present invention judge the device that object to be checked moves
11 offset data securing components
12 distance matrix computation modules
13 motion determination assemblies
14 movement position detection components
111 read modules
112 computing modules
113 offset data acquisition modules
Detailed description of the invention
For making the object, technical solutions and advantages of the present invention clearly, the present invention is described in more detail by the following examples.
Below for CT equipment, illustrate the present invention and judge the method and apparatus that object to be checked moves.For current double CT, the time that usual data collecting system rotates a circle around object synchronization to be checked is less than 1s, as being 0.6s, and human body moves frequency is generally one time 1 second, so in scanning process, not scan the exercise data all containing object to be checked in the Raw projection data group obtained at every turn.Unless object to be checked suffers from parkinsonism, quickly, the high-end CT needing sweep time very short scans it the motion frequency of the patient of this symptom.If all contain the exercise data of object to be checked in Raw projection data group, then rebuilding image can exist motion artifacts because of these exercise datas, thus the quality of image is rebuild in impact, so need to detect the exercise data whether containing object to be checked in original data set.
In view of this, the present invention proposes a kind of method whether object to be checked moves that judges, as shown in Figure 1.Described method comprises:
Step 101, reads the first projection angle α 1 and the first probe access n1 corresponding with a data for projection v1 in Raw projection data group R.
The present invention carries out based on the principle of compensation data as shown in Figure 2.Data collecting system shown in Fig. 2 comprises a kind of x-ray source 1 (F1 or F2) and detector 3, and the port number of detector 3 is M, and its initial probe access is C-1, and end probe access is C-M.Wherein, x-ray source 1 launches an X-ray beam, and X-ray beam subtended angle is 2 β, detector 3 together with x-ray source 1 around the clockwise synchronous rotary of an object to be checked, namely x-ray source is rotated clockwise to F2 place from F1, to gather X-ray projection data from a plurality of viewpoint.In theory, on same x-ray spread road on F1F2, due to X-ray from F1 to F2 direction through object to be checked decay with from F2 to F1 direction through object to be checked decling phase with, so when projection angle is α 1 (rotating to the anglec of rotation of F1 point from X-axis positive axis), the X-ray projection data v1 that a passage n1 of detector 3 receives should equal projection angle when being α 2 (rotating to the anglec of rotation of F2 point from X-axis positive axis), the X-ray projection data v2 of another passage n2 reception of detector 3, i.e. v1=v2.The present invention obtains the data of supplementary data group R_C accordingly.
Step 102, calculates by following formula (1) the angle θ that x-ray spread path F1F2 departs from X-ray beam center according to the first probe access n1:
θ=β-n1*Δβ (1)
Wherein, Δ β is adjacent two interchannel angle step pitches in detector 3, and β is a half angle of described X-ray beam.
Step 103, calculates the second probe access n2 according to angle theta by following formula (2):
n2=(β+θ)/Δβ (2)
Here, formula (1) can be deformed into following formula (3):
n1=(β-θ)/Δβ (3)
Relatively formula (2) and (3) can draw: the first probe access n1 and the second probe access n2 is specular about the center of detector 3.
Step 104, calculates the second projection angle α 2 according to angle theta by following formula (4):
α2=α1±pi-φ (4)
Due to φ=2 θ, so
α2=α1±pi-2θ (5)
α 2=α 1-pi-2 θ in the embodiment of Fig. 2.
Step 105, using the data for projection v2 that collects at the second projection angle α 2 times the second probe access n2 as the data in the offset data group R_C of Raw projection data group R.
Adopting above-mentioned steps 101 ~ 105, rotate a circle (360 °) around object to be checked without the need to comprising data collecting system, and the projection only need carrying out [0,180 ° of+2 β] can obtain the total data in offset data group R_C.So the scope of the first projection angle is [0,180 ° of+2 β] in the present invention, be also [0,180 ° of+2 β] according to the scope of known second projection angle of formula (5).Which reduces the radiation dose of object reception to be checked, decrease the data volume of Raw projection data simultaneously, decrease data space.
Step 106, calculates the distance matrix S of Raw projection data group R and offset data group R_C.
According to one embodiment of present invention, described distance matrix is Euclidean distance matrix.
In CT scan, data collecting system gathers X-ray projection data according to certain angle intervals, i.e. read tablet.As in single pass process, the present invention incites somebody to action [0,180 ° of+2 β] be divided into N number of angle intervals, as N=2000, then data collecting system often rotates ((180 ° of+2 β)/2000) ° around object to be checked, detector 3 just reads primary X-ray data for projection, altogether needs read tablet 2000 times.
Suppose that distance matrix S is the matrix of a N capable N row, wherein the data of the i-th row jth row can be the data of i-th read tablet in Raw projection data group R and Euclid (Euclidean) distance of the data of the secondary read tablet of jth in offset data group R_C.That is, in distance matrix S the data of the i-th row be i-th read tablet in Raw projection data group R data respectively with the Euclidean distance of the data of the 1st ~ N read tablet in offset data group R_C, in distance matrix S the data of jth row be jth time read tablet in offset data group R_C data respectively with the Euclidean distance of the data of the 1st ~ N read tablet in Raw projection data group R.
The data of certain distance matrix S i-th row jth row also can be the data of i-th read tablet in offset data group R_C and Euclid (Euclidean) distance of the data of the secondary read tablet of jth in Raw projection data group R.
Step 107, gets certain string j of distance matrix S, and find out the row value i in these row corresponding to minima, if the absolute value of the difference of this row value i and this train value j is greater than a predetermined threshold value, then described object to be checked moves.
Suppose J=1000, then in distance matrix S the 1000th column data be the 1000th read tablet in offset data group R_C data respectively with the Euclidean distance of the data of the 1st ~ N read tablet in Raw projection data group R.In theory, in distance matrix S, the minimum data of jth row is S
j, j, i.e. i=j.But in fact, because a variety of causes is as machine error (size of such as probe access, x-ray source focus are to the distance of scanning center, x-ray source focus to the error of the distance, X-ray beam subtended angle etc. of detector), accuracy rate that strip artifact detects, the speed etc. that flies focus, rotary frame, make the minimum data of jth row in S might not be S
j, j, but at S
j, jin a neighbouring preset threshold range.According to one embodiment of present invention, described predetermined threshold value is 20.Be the schematic diagram of the 1000th column data in S when object to be checked does not move as shown in Figure 3 a, wherein abscissa is the row value i of S, and vertical coordinate is S
i, jrelative size, near i=1000, have a circle in figure, this circle illustrates a preset threshold range, namely when in S the 1000th row minima in this circle, then think that object to be checked does not move in this scanning process.Otherwise, if the minima of the 1000th row is not in this circle, as shown in Figure 3 b, then think that object to be checked there occurs motion in this scanning process.Visible the present invention directly can just judge from Raw projection data group R whether object to be checked moves in present scan, accuracy rate and reliability higher.
Further, method of the present invention also comprises:
Step 108, if the absolute value of the difference of this row value i and this train value j is greater than described predetermined threshold value, then obtain according to this row value i the particular location that described object to be checked moves.
In S, the 1000th row minima is about i=1250 as shown in Figure 3 b, its circle compared near i=1000 has very large skew, known in present scan process during the about the 1250th read tablet object to be checked there occurs motion, and vertical coordinate has one significantly to jump (sharp break).
The present invention also provides a kind of computer program, the method making computer perform described judgement object to be checked to move.
The present invention also provides a kind of computer-readable recording medium, stores and judges the method instruction that object to be checked moves as described herein for making a machine perform.Particularly, the system or device of being furnished with storage medium can be provided, store the software program code of the function realizing any embodiment in above-described embodiment on the storage medium, and make the computer of this system or device (or CPU or MPU) read and perform the program code be stored in storage medium.In this case, can realize the function of any one embodiment above-described embodiment from the program code of storage medium reading itself, therefore program code and program code stored storage medium constitute a part of the present invention.
Floppy disk, hard disk, magneto-optic disk, CD (as CD-ROM, CD-R, CD-RW, DVD-ROM, DVD-RAM, DVD-RW, DVD+RW), tape, Nonvolatile memory card and ROM is comprised for providing the storage medium embodiment of program code.Selectively, can by communication network from download program code server computer.
In addition, be noted that, not only can by performing the program code read-out by computer, and by making the operating system etc. of calculating hands-operation carry out practical operation partly or completely based on the instruction of program code, thus the function of any one embodiment in above-described embodiment can be realized.
In addition, be understandable that, in the memorizer program code read by storage medium being write in memorizer set in the expansion board inserted in computer or arrange in writing the expanding element that is connected with computer, instruction subsequently based on program code makes the CPU etc. be arranged on expansion board or expanding element perform part and whole practical operation, thus realizes the function of any embodiment in above-described embodiment.
The present invention also provides a kind of device that object to be checked moves that judges, as shown in Figure 4, described device 10 comprises:
One offset data securing component 11, for obtaining the offset data group R_C of Raw projection data group R.
Particularly, offset data acquisition group 11 comprises as shown in Figure 5:
One read module 111, for reading the first projection angle α 1 and the first probe access n1 corresponding with a data for projection v1 in Raw projection data group R.
One computing module 112, for calculating one second projection angle α 2 and the second probe access n2 according to the first projection angle α 1 and the first probe access n1, the second probe access n2 wherein under the second projection angle α 2 is positioned on same x-ray spread path F1F2 with the first probe access n1 under the first projection angle α 1.Comprise particularly:
One angle calculation unit, for calculating by formula (1) the angle θ that x-ray spread path F1F2 departs from X-ray beam center according to the first probe access n1;
One second channel computing unit, for calculating the second probe access n2 according to angle theta by formula (2);
One second angle calculation unit, for calculating the second projection angle α 2 according to angle theta by formula (5).
One offset data acquisition module 113, the data for projection v2 that the second projection angle α that computing module 112 obtains by it 2 times the second probe access n2 collect is as the data in offset data group R_C.
One distance matrix computation module 12, for calculating the distance matrix S of Raw projection data group R and offset data group R_C.
One motion determination assembly 13, gets the jth row in described distance matrix, finds out the row value i in these row corresponding to minima, if the absolute value of the difference of this row value i and this train value j is greater than a predetermined threshold value, then thinks that object to be checked moves in this scanning process.According to one embodiment of present invention, predetermined threshold value is 20.Concrete judgment mode elaborates in the method for the invention, repeats no more here.
Further, described device 10 also comprises:
One movement position detection components 14, when the absolute value for the difference at this row value i and this train value j is greater than described predetermined threshold value, obtains according to this row value the particular location that described object to be checked moves.
As can be seen from such scheme, because the present invention directly can just judge from Raw projection data group whether object to be checked moves in present scan, thus accuracy rate and reliability higher; And described Raw projection data group is that data collecting system carries out [0 around object to be checked, 180 ° of+2 β] projection, the projection pattern comparing 360 ° decreases the radiation dose of object reception to be checked, decreases the data volume of Raw projection data simultaneously, decreases data space.
The present invention relates to field of medical device, particularly a kind ofly judge the method that object to be checked moves, described method comprises: the offset data group obtaining described Raw projection data group; Calculate the distance matrix of described Raw projection data group and described offset data group; Get certain string of described distance matrix, find out the row value in these row corresponding to minima, if the absolute value of the difference of this row value and this train value is greater than a predetermined threshold value, then there is motion artifacts in described Raw projection data group.The present invention also provides a kind of device that object to be checked moves that judges.Employing the present invention can directly from [0,180 ° of+2 β] Raw projection data group in just judge whether object to be checked moves in present scan, and obtain the particular location that moves further, decrease the X-ray radiation dosage of object reception to be checked simultaneously.
The foregoing is only preferred embodiment of the present invention, not in order to limit the present invention, within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (11)
1. judge to it is characterized in that the method that object to be checked moves, described method comprises:
Obtain the offset data group of Raw projection data group;
Calculate the distance matrix of described Raw projection data group and described offset data group;
Get certain string of described distance matrix, find out the row value in these row corresponding to minima, if the absolute value of the difference of the train value of this row value and these row is greater than a predetermined threshold value, then judge that object to be checked moves.
2. method according to claim 1, is characterized in that, described distance matrix is Euclidean distance matrix.
3. method according to claim 1, is characterized in that, described acquisition offset data group comprises:
Read the first projection angle corresponding with a data for projection in described Raw projection data group and the first probe access;
Calculate one second projection angle and the second probe access according to described first projection angle and described first probe access, the first probe access under the second probe access under wherein said second projection angle and described first projection angle is positioned on same x-ray spread path;
The data for projection collected by second probe access under described second projection angle is as the data in described offset data group.
4. method according to claim 3, is characterized in that, described calculating second projection angle and the second probe access comprise:
Calculate the angle at described x-ray spread path deviation X-ray beam center by following formula according to described first probe access:
θ=β-n1*Δβ,
Wherein, n1 is described first probe access, and Δ β is the angle step pitch of described probe access, and θ is the angle at described x-ray spread path and X-ray beam center, and β is a half angle of described X-ray beam;
Described second probe access is calculated by following formula according to described angle:
n2=(β+θ)/Δβ,
Wherein, n2 is described second probe access, and described second probe access and described first probe access are specular about the center of described detector;
Described second projection angle is calculated by following formula according to described angle:
α2=α1±pi-2θ,
Wherein, α 2 is described second projection angle, and α 1 is described first projection angle.
5. method according to claim 4, is characterized in that, the scope of described first projection angle is [0,180 ° of+2 β].
6. method according to claim 1, is characterized in that, described predetermined threshold value is 20.
7. method according to claim 1, is characterized in that, described method comprises further:
If the absolute value of the difference of this row value and this train value is greater than described predetermined threshold value, then obtain according to this row value the particular location that described object to be checked moves.
8. judge to it is characterized in that the device that object to be checked moves, described device (10) comprising:
One offset data securing component (11), for obtaining the offset data group of Raw projection data group;
One distance matrix computation module (12), for calculating the distance matrix of described Raw projection data group and described offset data group;
One motion determination assembly (13), gets certain string of described distance matrix, finds out the row value in these row corresponding to minima, if the absolute value of the difference of the train value of this row value and these row is greater than a predetermined threshold value, then judges that object to be checked moves.
9. device according to claim 8, is characterized in that, described offset data securing component (11) comprising:
One read module (111), for reading the first projection angle corresponding with a data for projection in described Raw projection data group and the first probe access;
One computing module (112), for calculating one second projection angle and one second probe access according to described first projection angle and described first probe access, the first probe access under the second probe access under wherein said second projection angle and described first projection angle is positioned on same x-ray spread path;
One offset data acquisition module (113), its data for projection collected by the second probe access under described second projection angle is as the data in described offset data group.
10. device according to claim 9, is characterized in that, described computing module (112) comprising:
One angle calculation unit, for calculating the angle at described x-ray spread path deviation X-ray beam center by following formula according to described first probe access:
θ=β-n1*Δβ,
Wherein, n1 is described first probe access, and Δ β is the angle step pitch of described probe access, and θ is the angle at described x-ray spread path and X-ray beam center, and β is a half angle of described X-ray beam;
One second channel computing unit, for calculating described second probe access according to described angle by following formula:
n2=(β+θ)/Δβ,
Wherein, n2 is described second probe access, and described second probe access and described first probe access are specular about the center of described detector;
One second angle calculation unit, for calculating described second projection angle according to described angle by following formula:
α2=α1±pi-2θ,
Wherein, α 2 is described second projection angle, and α 1 is described first projection angle.
11. devices according to claim 8, is characterized in that, described device (10) comprises further:
One movement position detection components (14), when the absolute value for the difference in this row value and this train value is greater than described predetermined threshold value, obtains according to this row value the particular location that described object to be checked moves.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110145206.XA CN102805634B (en) | 2011-05-31 | 2011-05-31 | Method and device for judging motion of object to be inspected |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110145206.XA CN102805634B (en) | 2011-05-31 | 2011-05-31 | Method and device for judging motion of object to be inspected |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102805634A CN102805634A (en) | 2012-12-05 |
CN102805634B true CN102805634B (en) | 2015-01-21 |
Family
ID=47229630
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110145206.XA Active CN102805634B (en) | 2011-05-31 | 2011-05-31 | Method and device for judging motion of object to be inspected |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102805634B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102805635B (en) * | 2011-05-31 | 2015-04-08 | 上海西门子医疗器械有限公司 | Method and device for optimizing scanning flow |
CN106821407A (en) * | 2016-12-28 | 2017-06-13 | 上海联影医疗科技有限公司 | For the method for testing motion and device of computed tomography |
CN107481226B (en) * | 2017-07-27 | 2021-06-01 | 东软医疗系统股份有限公司 | Method and device for removing abnormal scanning data and PET system |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1787077A (en) * | 2005-12-13 | 2006-06-14 | 浙江大学 | Method for fast identifying speeking person based on comparing ordinal number of archor model space projection |
CN1895174A (en) * | 2005-06-10 | 2007-01-17 | 通用电气公司 | Step-and-shoot cardiac ct imaging |
CN1918605A (en) * | 2004-02-13 | 2007-02-21 | 皇家飞利浦电子股份有限公司 | Motion artifact compensation |
CN101222703A (en) * | 2007-01-12 | 2008-07-16 | 杭州波导软件有限公司 | Identity verification method for mobile terminal based on voice identification |
CN102755170A (en) * | 2011-04-29 | 2012-10-31 | 上海西门子医疗器械有限公司 | Detector and X-ray projection data acquiring system comprising same |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6904118B2 (en) * | 2002-07-23 | 2005-06-07 | General Electric Company | Method and apparatus for generating a density map using dual-energy CT |
JP3944173B2 (en) * | 2004-02-05 | 2007-07-11 | キヤノン株式会社 | Radiation image processing apparatus and processing method |
US7085342B2 (en) * | 2004-04-22 | 2006-08-01 | Canamet Canadian National Medical Technologies Inc | Method for tracking motion phase of an object for correcting organ motion artifacts in X-ray CT systems |
JP4539186B2 (en) * | 2004-06-17 | 2010-09-08 | コニカミノルタエムジー株式会社 | Medical imaging system, information processing device |
-
2011
- 2011-05-31 CN CN201110145206.XA patent/CN102805634B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1918605A (en) * | 2004-02-13 | 2007-02-21 | 皇家飞利浦电子股份有限公司 | Motion artifact compensation |
CN1895174A (en) * | 2005-06-10 | 2007-01-17 | 通用电气公司 | Step-and-shoot cardiac ct imaging |
CN1787077A (en) * | 2005-12-13 | 2006-06-14 | 浙江大学 | Method for fast identifying speeking person based on comparing ordinal number of archor model space projection |
CN101222703A (en) * | 2007-01-12 | 2008-07-16 | 杭州波导软件有限公司 | Identity verification method for mobile terminal based on voice identification |
CN102755170A (en) * | 2011-04-29 | 2012-10-31 | 上海西门子医疗器械有限公司 | Detector and X-ray projection data acquiring system comprising same |
Non-Patent Citations (2)
Title |
---|
Eberhard Hansis,et al.Motion-Compensated Iterative Sparse Data Reconstruction for Interventional 3-D Coronary Artery Imaging.《2008 IEEE Nuclear Science Symposium Conference Record》.2008,第4280-4284页. * |
Oliver Morin,et al.MEGVOLTAGE CONE-BEAM CT:SYSTEM DESCRIPTION AND CLINICAL APPLICATIONS.《Medical Dosimetry》.2006,第31卷(第1期),第51-61页. * |
Also Published As
Publication number | Publication date |
---|---|
CN102805634A (en) | 2012-12-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8699659B2 (en) | Systems and methods for focal spot motion correction | |
US9427194B2 (en) | Apparatus and methods for computed tomography imaging | |
CN100545745C (en) | X-ray imaging apparatus and control method thereof | |
US9058658B2 (en) | Methods and devices for locating object in CT imaging | |
US20070092127A1 (en) | Method and device for segmenting at least one substance in an x-ray image | |
CN102598059A (en) | Method and apparatus for generating computed tomography images with offset detector geometries | |
CN100502785C (en) | Tomography device and method | |
US20040260176A1 (en) | Systems and methods for correcting a positron emission tomography emission image | |
CN103596502A (en) | Adaptive calibration for tomographic imaging systems | |
CN102123664A (en) | Calibration method for ring artifact correction in non-ideal isocentric 3D rotational x-ray scanner systems using a calibration phantom based rotation center finding algorithm | |
US20070253525A1 (en) | Method for scattered radiation correction of a CT system having at least two focus/detector systems arranged with an angular offset, and a CT system | |
CN101495041B (en) | Gated ct with irregular sampling for slow ct acquisition | |
JP2005058653A (en) | X-ray ct apparatus, information processing method, storage medium and program | |
US8279999B2 (en) | X-ray CT apparatus and a method of controlling the X-ray CT apparatus | |
US7583784B2 (en) | Method for calculating computed tomography pictures from detector data of a CT having at least two radiation sources | |
CN111728632B (en) | Radiation detection device, radiation detection method and CT image reconstruction method | |
US20080174310A1 (en) | Method for operating a hybrid medical imaging unit comprising a first imaging device of high spatial resolution and a second nuclear medicine imaging device of high sensitivity | |
CN102805634B (en) | Method and device for judging motion of object to be inspected | |
US8953739B2 (en) | Image processing apparatus and X-ray CT system | |
US20080240340A1 (en) | Method for scattered radiation correction in x-ray imaging devices | |
US20140211912A1 (en) | Computerized tomography (ct) method and ct system | |
CN102805635B (en) | Method and device for optimizing scanning flow | |
CN102894989B (en) | Method and device for reducing scanning dose of X-ray | |
US10217248B2 (en) | Method for removing streak from detector cell with performance difference | |
CN104797194A (en) | X-ray ct device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |