CN108460740A - CT spiral reconstruction image artifacts minimizing technologies - Google Patents
CT spiral reconstruction image artifacts minimizing technologies Download PDFInfo
- Publication number
- CN108460740A CN108460740A CN201810182895.3A CN201810182895A CN108460740A CN 108460740 A CN108460740 A CN 108460740A CN 201810182895 A CN201810182895 A CN 201810182895A CN 108460740 A CN108460740 A CN 108460740A
- Authority
- CN
- China
- Prior art keywords
- row
- data
- detector
- projection
- rows
- 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
- 238000005516 engineering process Methods 0.000 title claims abstract description 12
- 238000000034 method Methods 0.000 claims abstract description 8
- 238000002591 computed tomography Methods 0.000 description 15
- 238000010586 diagram Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 210000004209 hair Anatomy 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001360 synchronised effect 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/73—Deblurring; Sharpening
-
- 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]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The present invention provides a kind of CT spiral reconstructions image artifacts minimizing technology, i.e.,：To in multi-detector, the data for projection of thin row's detector cells acquisition is without processing；By in multi-detector, the data for projection of thickness row's detector cells acquisition is expanded.Specific method is：The data for projection of every row's thickness row's detector cells acquisition is extended for N number of rows evidences；Wherein N is：The thickness of thickness row's detector cells divided by the thickness of thin row's detector cells；The data for projection of thick row's detector cells acquisition is expanded into N number of rows evidences, per number of rows according to associated with immediate two rows of former data for projection of detector cells acquisition, immediate data weighting accounting is big in the data for projection of two rows of original detector cells acquisitions, secondary close data weighting accounting is small, weight and be 1.Advantages of the present invention：Processing speed is fast, the windmill artifact that can be effectively removed in image, while the noise level of image does not have significant change.
Description
Technical field
The present invention relates to a kind of CT spiral reconstructions image artifacts minimizing technologies.
Background technology
In third generation CT (Computed Tomography) machine, spiral is generally carried out using the detector of multiple rows of structure
Track scanning acquires the data for projection of patient, and then, reconstruction obtains the faultage image of patient.In data acquisition, CT machines
Light source (i.e. bulb) is synchronous with detector to be rotated around CT bed boards, meanwhile, light source sends out X-ray, the constraint of collimator, only
Ray in the range of needing passes through collimator to scan patient, reaches the receiving area of detector.
In actual clinical scanning, as shown in Figure 1A, when using small opening collimator, if the only centre of detector
Dry row's detector cells can receive X-ray signal.As shown in Figure 1B, when using big opening collimator when, detector it is complete
Portion's row's detector cells can receive the signal of X-ray.
Currently, third generation CT machines generally use multi-detector, number of rows to have 16 rows, 24 rows, 32 rows, 64 rows, 128 rows
Deng.In general, the multi-detector that third generation CT machines use respectively arranges the third that the thickness of detector cells is consistent, but also has
For CT machines multi-detector, row and the thickness of both sides row's detector cells are different among it, such as 24 row shown in Fig. 2 visits
Survey device.24 row detector shown in Fig. 2, each small lattice represent a physical detector unit in figure, each detector cells according to
The secondary transversely arranged different channels for constituting detector, longitudinal arrangement constitute the different rows of detector, in total 24 row.Dark grey
Partly width of each detector cells of (16 intermediate rows) and both sides bright gray parts (each 4 row up and down) in channel orientation
It is identical, still, the thickness of detector cells is different on longitudinally row direction.Respectively row visits intermediate row (dark gray section)
It is identical to survey device unit thickness, and it is intermediate row (Dark grey portion that both sides, which arrange (bright gray parts) and often arrange the thickness of detector cells,
Point) often arrange 2 times of detector cells thickness.
For small opening collimator as shown in Figure 1A, due to the X-ray that light source is sent out, through small opening collimator
Constraint passes through patient, is only received by that identical a few row's detector cells of multi-detector interior thickness, multi-detector acquisition
After the data for projection of patient, the patient's faultage image rebuild is than more visible.But it is accurate for big opening as shown in Figure 1B
For straight device, the X-ray that light source is sent out, the constraint through big opening collimator passes through patient, is visited by whole rows of multi-detector
Device unit is surveyed to receive, since the thickness that detector cells are arranged on multi-detector both sides is 2 times of intermediate row detector cells thickness,
Therefore according to common practice, in the image reconstruction under big opening collimator as shown in Figure 1B, the image smallest tier of system reconstructing
Thickness is 2 times of reconstruction image thickness under small opening collimator shown in Figure 1A；Meanwhile in data acquisition, centre is incorporated
The data of thin layer detector cells acquisition make its every two layers of adjoining platelet detector cells generate a thickness layer data.This
In the case of, the axial data collection rate of system declines, and the image reconstructed is caused to will appear sternly at the position of tissue change distance
The windmill artifact of weight, the position as shown in white arrow in Fig. 3.
Invention content
In view of the foregoing, the object of the present invention is to provide a kind of CT spiral reconstructions image artifacts minimizing technologies.This method
The puppet generated when can be effectively removed using the intermediate row multi-detector gathered data reconstruction image different from both sides row's thickness
Shadow.
To achieve the above object, the present invention uses following technical scheme：A kind of CT spiral reconstructions image artifacts minimizing technology,
The artifact minimizing technology is：To in multi-detector, the data for projection of thin row's detector cells acquisition is without processing；It will be multiple rows of
In detector, the data for projection of thickness row's detector cells acquisition is expanded.
Preferably, the data for projection by the row's detector cells acquisition of every row's thickness is extended for N number of rows evidences；Wherein N is：Thickness row visits
Survey the thickness of device unit divided by the thickness of thin row's detector cells；The data for projection of thick row's detector cells acquisition is expanded into N
Number of rows evidence, per number of rows according to, two rows of original detector lists associated with immediate two rows of former data for projection of detector cells acquisition
Member acquisition data for projection in immediate data weighting accounting it is big, secondary close data weighting accounting is small, weight and be 1.
Preferably, the N is 2.
Description of the drawings
Figure 1A is bulb detector light path schematic diagram under small opening collimator；
Figure 1B is bulb detector light path schematic diagram under big opening collimator；
Fig. 2 is the intermediate row 24 row panel detector structure schematic diagrames different from both sides row's thickness；
Fig. 3 is the CT faultage images with windmill artifact；
Fig. 4 is the schematic diagram of present invention removal image artifacts；
Fig. 5 is the image using the pseudo- movie queen of present invention removal.
Specific implementation mode
The structure and feature of the present invention are described in detail with reference to the accompanying drawings and examples.It should be noted that can
To make various modifications to disclosed embodiments, therefore, embodiment disclosed in specification should not be considered as to the present invention
Limitation, and only as the example of embodiment, the purpose is to keep the feature of the present invention apparent.
Artifact is generated when the different multi-detector gathered data reconstruction image of analysis and utilization intermediate row and both sides row's thickness
The reason of, the inventors discovered that：When carrying out data scanning image reconstruction under big opening collimator, in order to ensure the projection of acquisition
The consistency of thickness that data are often arranged, it is common practice that merge thin intermediate data, this can cause axial sample rate to decline, to filter
Apparent windmill artifact is generated in the image of wave backprojection reconstruction.
For the artifact in removal image, the method for removal artifact provided by the invention：
S1, in multi-detector, the data for projection of thin row's detector cells acquisition is without processing；
S2, by multi-detector, the data for projection of thickness row's detector cells acquisition is expanded, and the purpose of expansion is to increase
Add axial data sampling rate；
Specific method is：
S21, the data for projection of every row's thickness row's detector cells acquisition is extended for N number of rows evidences；Wherein N is：Thickness row's detection
The thickness of the thickness of device unit divided by thin row's detector cells；
That is, if the thickness of multi-detector thickness row's detector cells is 2 times of thin row's detector cells thickness
(i.e. N=2), such as 24 row's multi-detectors, the 1st the-the 4 row of row of both sides and the 21st-the 24 row's detector cells of row are (light grey
Part) thickness be 2 times of intermediate row the 5th row-the 20 row's detector cells (dark gray section) thickness；Then, for centre the 5th
The data for projection of the thin row's detector cells acquisition of the-the 20 row (totally 16 row) is arranged without processing；To the 1st the-the 4 row of row of both sides and
The data for projection of 21-the 24 row's thickness of row row's detector cells acquisition expands, by the acquisition of the 1st-the 4 row's detector cells of row
Data for projection is expanded into original 2 times, is expanded into 8 row's data for projection, and similarly, the 21st-the 24 row's detector cells of row are acquired
Data for projection be also expanded into original 2 times, be expanded into 8 row's data for projection.8 row detector data for projection original so just expands
It is original 2 times at 16 row's detector data for projection, originally 24 row's data for projection are expanded into 8+16+8=32 row's data for projection,
The 32 row's data for projection obtained after expansion are reconstructed into patient's CT faultage images according to usual way.
S22, the data for projection of thick row's detector cells acquisition is expanded into N rows (i.e. N is 2) data, per number of rows evidence with most
The data for projection of the former detector cells acquisition of close two rows is associated, in the data for projection of two rows of original detector cells acquisitions most
Close data weighting accounting is big, and secondary close data weighting accounting is small, weight and be 1, i.e.,：
B4+l=Dl4 ＜ l ＜ 21
Wherein,
w1+w2=1, and w1＞ w2
DlIndicate the data of original l row's detector cells；
B2l-1、B2lIndicate the detector cells data after original l row's detector cells expand.
It illustrates below in conjunction with the accompanying drawings, how by multi-detector, the data of thickness row's detector cells acquisition carry out
Expand.The first row indicates the initial data of 24 row's detectors acquisition in Fig. 4, uses D1、D2、D3、D4、D5、D6、……、D19、D20、
D21、D22、D23、D24Represent 24 row's detectors acquisition the 1st row, the 2nd row, the 3rd row, the 4th row, the 5th row, the 6th row ..., the 19th
Row, the 20th row, the 21st row, the 22nd row, the 23rd row, the 24th number of rows evidence.
Since the thickness of 24 row's detector both sides thickness row's detector cells (the 1st the-the 4 row of row and the 21st the-the 24 row of row) is
Intermediate thin arranges 2 times of detector cells (the 5th the-the 20 row of row) thickness, therefore, as shown in the second row in Fig. 4, intermediate thin is arranged and is detected
The data of the 5th row of row-the 20th of device unit are not handled, and both sides thickness is arranged the 1st the-the 4 row of row of detector cells and the 21st row-the 24
1 times of the data extending of row, uses B1、B2、B3、B4、B5、B6、B7、B8、B9、B10、B11、B12、……、B21、B22、B23、B24、B25、B26、
B27、B28、B29、B30、B31、B32Represent the 1st row --- the 32nd number of rows evidence after expanding.
Specifically, by former 1st number of rows according to D1It is extended for B1、B2；Former 2nd number of rows is according to D2It is extended for B3、B4；Former 3rd number of rows
According to D3It is extended for B5、B6；Former 4th number of rows is according to D4It is extended for B7、B8；The data of former 5th row of row-the 20th are not handled, D5=B9、D6=
B10、D7=B11、……、D19=B23、D20=B24；Former 21st number of rows is according to D21It is extended for B25、B26；Former 22nd number of rows is according to D22Expand
For B27、B28；Former 23rd number of rows is according to D23It is extended for B29、B30；Former 24th number of rows is according to D24It is extended for B31、B32。
According to formula：
B4+l=Dl4 ＜ l ＜ 21
In the specific embodiment of the invention, W1=0.75, W2=0.25, it is of course also possible to other numerical value be taken, as long as w1+w2
=1, and w1＞ w2.
As shown in figure 4, the present invention is by D1It is extended for B1And B2, B1=D1, B2=0.75D1+0.25D2；By D2It is extended for B3With
B4, B3=0.75D2+0.25D1, B4=0.75D2+0.25D3；
By D3It is extended for B5And B6, B5=0.75D3+0.25D2, B6=0.75D3+0.25D4；
By D4It is extended for B7And B8, B7=0.75D4+0.25D3, B8=0.75D4+0.25D5；
B9=D5；B10=D6；B11=D7；……；B22=D18；B23=D19；B24=D20；
By D21It is extended for B25And B26, B25=0.75D21+0.25D20, B26=0.75D21+0.25D22；
By D22It is extended for B27And B28, B27=0.75D22+0.25D21, B28=0.75D22+0.25D23；
By D23It is extended for B29And B30, B29=0.75D23+0.25D22, B30=0.75D23+0.25D24；
By D24It is extended for B31And B32, B31=0.75D24+0.25D23, B32=D24。
For intermediate row and both sides row's thickness, 16 different rows, 32 rows, 64 rows, 128 row's detectors can be according to this hairs
Artifact in the method removal reconstruction image of bright offer, i.e.,：For thin row's detector cells acquisition data for projection without place
Reason；The data for projection of thick row's detector cells acquisition is expanded.Due to 16 row currently used in the market, 32 rows, 64 rows,
128 row's detectors, the thickness of thickness row's detector cells are 2 times of thin row's detector cells thickness, therefore, it usually will often arrange thick row
The data extending of detector cells acquisition makes the thickness of every number of rows evidence after expansion be adopted with thin row's detector cells at 2 number of rows evidences
The data thickness of collection is identical.
The present invention need not carry out any modification to existing CT machines structure, also need not carry out any change to CT algorithm for reconstructing
Become, it is only necessary to which the data collected are expanded.
Advantages of the present invention：
1, as shown in figure 5, the windmill artifact that can be effectively removed in image, while the noise level of image does not obviously become
Change.
2, method provided by the invention is only handled projection numeric field data, need not be handled image domain data,
Therefore processing speed is fast, does not influence the reconstruction speed of image, is suitable for scan data, the image reconstruction of all different parts.
3, simple linear interpolation method is used during expanding projection numeric field data, ensure that the Shandong of algorithm
Stick does not lose resolution ratio substantially to reconstruction image.
It should be noted that the present invention is suitable for all non-uniform thickness panel detector structures, it is not limited to 24 rows (intermediate 16 rows
Each 4 row thick-layer in thin layer both sides) detector.
Finally it should be noted that：Above-described embodiments are merely to illustrate the technical scheme, rather than to it
Limitation；Although the present invention is described in detail referring to the foregoing embodiments, it will be understood by those of ordinary skill in the art that：
It can still modify to the technical solution recorded in previous embodiment, or to which part or all technical features into
Row equivalent replacement；And these modifications or substitutions, it does not separate the essence of the corresponding technical solution various embodiments of the present invention technical side
The range of case.
Claims (3)
1. a kind of CT spiral reconstructions image artifacts minimizing technology, it is characterised in that：The artifact minimizing technology is：To multi-detector
In, the data for projection of thin row's detector cells acquisition is without processing；By in multi-detector, thickness row's detector cells acquire
Data for projection is expanded.
2. CT spiral reconstructions image artifacts minimizing technology according to claim 1, it is characterised in that：The extending method
For：
S21, the data for projection of every row's thickness row's detector cells acquisition is extended for N number of rows evidences；Wherein N is：Thickness row's detector list
The thickness of member divided by the thickness of thin row's detector cells；
S22, the thick data for projection for arranging detector cells acquisition is expanded into N number of rows evidences, evidence and immediate two rows of originals per number of rows
The data for projection of detector cells acquisition is associated, immediate data power in the data for projection of two rows of original detector cells acquisitions
Weight accounting it is big, secondary close data weighting accounting is small, weight and be 1.
3. CT spiral reconstructions image artifacts minimizing technology according to claim 2, it is characterised in that：The N is 2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810182895.3A CN108460740B (en) | 2018-03-06 | 2018-03-06 | CT spiral reconstruction image artifact removing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810182895.3A CN108460740B (en) | 2018-03-06 | 2018-03-06 | CT spiral reconstruction image artifact removing method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108460740A true CN108460740A (en) | 2018-08-28 |
| CN108460740B CN108460740B (en) | 2021-12-14 |
Family
ID=63217239
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810182895.3A Active CN108460740B (en) | 2018-03-06 | 2018-03-06 | CT spiral reconstruction image artifact removing method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108460740B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111000578A (en) * | 2019-12-25 | 2020-04-14 | 东软医疗系统股份有限公司 | Image reconstruction method and device, CT (computed tomography) equipment and CT system |
| CN112884855A (en) * | 2021-01-13 | 2021-06-01 | 中广核贝谷科技有限公司 | Processing method and device for security check CT reconstructed image |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1669528A (en) * | 2004-03-19 | 2005-09-21 | 深圳安科高技术股份有限公司 | Image reconstruction method in double-line or multi-line helical CT |
| CN100481130C (en) * | 2004-01-29 | 2009-04-22 | 皇家飞利浦电子股份有限公司 | Method and device for reducing windmill artifact in multi-slice CT reconstruction |
| US20090110257A1 (en) * | 2007-10-29 | 2009-04-30 | Kabushiki Kaisha Toshiba | Interpolation interlacing based data upsampling algorithm for cone-beam x-ray ct flying focal spot projection data |
| CN104605881A (en) * | 2014-12-31 | 2015-05-13 | 沈阳东软医疗系统有限公司 | Parameter optimizing method and medical equipment |
| CN104978717A (en) * | 2015-06-11 | 2015-10-14 | 沈阳东软医疗系统有限公司 | CT reconstruction image processing method, apparatus and device |
-
2018
- 2018-03-06 CN CN201810182895.3A patent/CN108460740B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100481130C (en) * | 2004-01-29 | 2009-04-22 | 皇家飞利浦电子股份有限公司 | Method and device for reducing windmill artifact in multi-slice CT reconstruction |
| CN1669528A (en) * | 2004-03-19 | 2005-09-21 | 深圳安科高技术股份有限公司 | Image reconstruction method in double-line or multi-line helical CT |
| US20090110257A1 (en) * | 2007-10-29 | 2009-04-30 | Kabushiki Kaisha Toshiba | Interpolation interlacing based data upsampling algorithm for cone-beam x-ray ct flying focal spot projection data |
| CN104605881A (en) * | 2014-12-31 | 2015-05-13 | 沈阳东软医疗系统有限公司 | Parameter optimizing method and medical equipment |
| CN104978717A (en) * | 2015-06-11 | 2015-10-14 | 沈阳东软医疗系统有限公司 | CT reconstruction image processing method, apparatus and device |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111000578A (en) * | 2019-12-25 | 2020-04-14 | 东软医疗系统股份有限公司 | Image reconstruction method and device, CT (computed tomography) equipment and CT system |
| CN111000578B (en) * | 2019-12-25 | 2023-05-02 | 东软医疗系统股份有限公司 | Image reconstruction method, device, CT equipment and CT system |
| CN112884855A (en) * | 2021-01-13 | 2021-06-01 | 中广核贝谷科技有限公司 | Processing method and device for security check CT reconstructed image |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108460740B (en) | 2021-12-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8705822B2 (en) | Method for creating images indicating material decomposition in dual energy, dual source helical computed tomography | |
| US7599464B2 (en) | Method and computed tomography system for producing tomograms of an object | |
| US20080019474A1 (en) | X-ray ct scanner and data processing method of x-ray ct scanner | |
| CN102631210B (en) | Method, image data record processing facility, x-ray system and for correcting image data of an examination object | |
| CN102144928B (en) | CT measurement with multiple X-ray sources | |
| US9271681B2 (en) | X-ray computed tomography apparatus | |
| US7020234B2 (en) | Method for producing tomograms of a periodically moving object with the aid of a focus detector combination | |
| US20090232269A1 (en) | Methods and apparatus for noise estimation for multi-resolution anisotropic diffusion filtering | |
| CN102576468B (en) | Method for artifact reduction in cone-beam CT images | |
| US7298813B2 (en) | X-ray computed tomographic apparatus, image processing apparatus, and image processing method | |
| CN108460740A (en) | CT spiral reconstruction image artifacts minimizing technologies | |
| US10426417B2 (en) | Computed tomography (CT) hybrid data acquisition | |
| CN105761226B (en) | A kind of compensated reconstruction method of ultraphotic open country CT scan image | |
| US20160307340A1 (en) | Multispectral ct imaging | |
| EP2490180B1 (en) | Medical image processing apparatus and medical image imaging method | |
| US20050175141A1 (en) | Method for producing tomograms of a periodically moving object with the aid of a focus/detector combination | |
| JP4280018B2 (en) | X-ray computed tomography system | |
| JP2010142478A (en) | X-ray ct apparatus | |
| RU2708816C2 (en) | Detector and visualization system for x-ray phase-contrast imaging of tomosynthesis | |
| Kalender et al. | Spiral CT: medical use and potential industrial applications | |
| US7822467B2 (en) | Method for producing CT images of a cyclically moving object to be examined | |
| JP6143533B2 (en) | Nuclear medicine image reconstruction device, nuclear medicine image reconstruction method, and program | |
| CN102670230B (en) | Method for reducing motion artifact in dect | |
| US7924970B2 (en) | Method and device for generating a CT image with a high time resolution | |
| JP4550179B2 (en) | X-ray computed tomography system |
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: 100176 1st floor, building 8, 11 Kangding street, Beijing Economic and Technological Development Zone, Daxing District, Beijing Applicant after: Sinovision Technology (Beijing) Co.,Ltd. Address before: 100176 3 / F, building 1, Yuehong building, No.13, Yongchang North Road, Beijing Economic and Technological Development Zone, Daxing District, Beijing Applicant before: SAINUO WEISHENG TECHNOLOGY (BEIJING) Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |