CN106910164A - A kind of method and device being filtered to CT data for projection - Google Patents
A kind of method and device being filtered to CT data for projection Download PDFInfo
- Publication number
- CN106910164A CN106910164A CN201510979443.4A CN201510979443A CN106910164A CN 106910164 A CN106910164 A CN 106910164A CN 201510979443 A CN201510979443 A CN 201510979443A CN 106910164 A CN106910164 A CN 106910164A
- Authority
- CN
- China
- Prior art keywords
- data
- projection
- high frequency
- filtering
- filtered
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000001914 filtration Methods 0.000 claims description 97
- 238000002591 computed tomography Methods 0.000 description 57
- 238000010586 diagram Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 239000012634 fragment Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 210000000746 body region Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G06T5/70—
-
- 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]
Abstract
The present invention relates to a kind of method and device being filtered to CT data for projection, including：Obtain the high frequency trajectory in CT data for projection；And the trend based on high frequency trajectory is filtered to CT data for projection.
Description
Technical field
The present invention relates to CT technical fields, more particularly to a kind of side being filtered to CT data for projection
Method and device.
Background technology
In computed tomography (CT) equipment, detector (detector) can collect those and wear
The X-ray of scanned object is crossed, and converts thereof into electric signal.After these electric signals are recorded
Just constitute CT initial data, also referred to as CT data for projection.CT is thrown by corresponding algorithm for reconstructing
Shadow data are rebuild, it is possible to obtain CT reconstruction images.
Generally, before image reconstruction is carried out using data for projection, treatment can be filtered to data for projection,
To suppress noise therein and to strengthen useful signal.
The existing method being filtered to CT data for projection, is along the pay-out angle in data for projection
Passage (channel) direction of degree (view) direction and/or detector is filtered.
But, the effect of existing method is barely satisfactory, for the puppet of the types such as smudge, center dot
The treatment effect of shadow is also bad.
It is therefore desirable to a kind of new method and device being filtered to CT data for projection is provided, so as to
The noise in CT data for projection can preferably be filtered and useful signal is kept.
The content of the invention
An embodiment provides a kind of method being filtered to CT data for projection, including：
Obtain the high frequency trajectory in CT data for projection；And the trend based on high frequency trajectory is to CT data for projection
It is filtered.
Another embodiment of the present invention provides a kind of device being filtered to CT data for projection, including：
High frequency trajectory acquisition module, for obtaining the high frequency trajectory in CT data for projection；And high frequency trajectory filter
Ripple module, is filtered for the trend based on high frequency trajectory to CT data for projection.
Brief description of the drawings
It is described for embodiments of the invention by with reference to accompanying drawing, the present invention may be better understood,
In the accompanying drawings：
Fig. 1 show the stream of one embodiment of the method being filtered to CT data for projection of the invention
Journey schematic diagram；
Fig. 2 show it is of the invention during being filtered to CT data for projection based on high frequency trajectory
The schematic flow sheet of one embodiment that is filtered to CT data for projection of trend；
Fig. 3 show it is of the invention during being filtered to CT data for projection based on high frequency trajectory
The schematic flow sheet of another embodiment that is filtered to CT data for projection of trend；
Fig. 4 show showing for one embodiment of the device being filtered to CT data for projection of the invention
Meaning property block diagram.
Specific embodiment
Specific embodiment of the invention explained below, it should be pointed out that in these implementation methods
Specific descriptions during, in order to carry out brief and concise description, this specification can not possibly be to actual
All features of implementation method make description in detail.It is to be understood that in any one reality
In applying the actual implementation process of mode, as the process in any one engineering project or design object
In, in order to realize the objectives of developer, the limit related or commercially related in order to meet system
System, can usually make various specific decision-makings, and this also can be from a kind of implementation method to another kind
Changed between implementation method.Although moreover, it is to be understood that institute in this development process
The effort made is probably complicated and tediously long, but related for present disclosure
For one of ordinary skill in the art, carried out on the basis of the technology contents that the disclosure is disclosed
A little designs, the change such as manufacture or production is conventional technology, is not construed as the disclosure
Content it is insufficient.
Unless otherwise defined, the technical term or scientific terminology for being used in claims and specification
It should be the ordinary meaning that the personage with general technical ability is understood in the technical field of the invention.This
" first ", " second " and the class used in application for a patent for invention specification and claims
As word be not offered as any order, quantity or importance, and be used only to distinguish different groups
Into part.The similar word such as " one " or " one " is not offered as quantity limitation, but expression is deposited
At least one.The similar word such as " including " or "comprising" mean to appear in " including " or
Element or object before person's "comprising" cover and appear in " including " or "comprising" behind arrange
The element or object and its equivalent element of act, it is not excluded that other elements or object." connection "
Or the similar word such as " connected " is not limited to physics or machinery connection, is also not necessarily limited to
It is direct or indirectly connects.
To make the object, technical solutions and advantages of the present invention clearer, below in conjunction with of the invention specific
Embodiment and corresponding accompanying drawing are clearly and completely described to technical solution of the present invention.Obviously, retouched
The embodiment stated only is a part of embodiment of the invention, rather than whole embodiments.Based in the present invention
Embodiment, it is all that those of ordinary skill in the art are obtained under the premise of creative work is not made
Other embodiment, belongs to the scope of protection of the invention.
A kind of embodiments in accordance with the present invention, there is provided method being filtered to CT data for projection.
With reference to Fig. 1, Fig. 1 show the method 100 being filtered to CT data for projection of the invention
The schematic flow sheet of one embodiment.Method 100 can comprise the following steps 101 to 102.
As shown in figure 1, in a step 101, obtaining the high frequency trajectory in CT data for projection.
In one embodiment of the invention, height can be extracted by way of CT reconstruction images are guided
Frequency track.Specifically, CT data for projection first can be obtained into reconstruction image, Ran Hou by reconstruction
The region that density higher building body region or object density are changed greatly is identified in reconstruction image, is entered
And orthographic projection is carried out to these regions and can be obtained by high frequency trajectory.
In another embodiment of the present invention, high frequency trajectory directly can be extracted on CT sinograms.
Specifically, can directly identify that brightness is larger on sinogram according to threshold value set in advance
With less curve, these curves can be considered to be high frequency trajectory.
In a step 102, the trend based on high frequency trajectory is filtered to CT data for projection.
In one embodiment of the invention, CT data for projection can be entered along the trend of high frequency trajectory
Row filtering.That is, when filtering box (box) is applied into CT data for projection being filtered,
Filtering box can be filtered along the direction of moving towards of high frequency trajectory to whole high frequency trajectory, or can edge
The trend for high frequency trajectory needs the fragment for being filtered treatment to be filtered some in high frequency trajectory,
And being not as prior art is filtered along View or Channel directions.Needs are filtered
The fragment of ripple treatment can be can be with the weaker piece of predetermined useful signal according to CT application scenarios
Section.The data for projection that the embodiment can be used for being filtered is probably useful signal and needs to increase it
Strong scene.
In another embodiment of the present invention, can also be along the side perpendicular with the tangent line of high frequency trajectory
It is filtered to CT data for projection.That is, work as being filtered to the certain point in high frequency trajectory
When, the trend of filtering box can be perpendicular to the tangent line of the point.By that analogy, to being needed in high frequency trajectory
When each to be filtered point is filtered, the trend of filtering box can cutting perpendicular to the point
Line.The data for projection that the embodiment can be used for being filtered is probably noise and needs to suppress it
Scene.
As for the part in CT data for projection in addition to high frequency trajectory, low frequency portion is generally considered
Divide, also referred to as background data.For these data, non-directional filtering can be done, i.e.,：In each filtering
Equal weight is imposed on direction.
In one embodiment of the invention, with reference to Fig. 2, step 102 can also include following sub-step
201 to 203.
In sub-step 201, selected on 360 degree of directions on the sinogram being made up of CT data for projection
Take multiple filtering directions.
For example, on 0 to 360 direction, 0 degree, 45 degree, the multiple directions such as 90 degree can be chosen
As filtering direction.
In sub-step 202, in multiple filtering directions, pair immediate filter is moved towards with high frequency trajectory
Ripple direction sets larger weighted value and other filtering directions in multiple filtering direction is set compared with little Quan
Weight values.
In the multiple filtering directions selected by sub-step 201, at least there is filtering direction and a high frequency
Track to move towards direction closest or even overlap, then larger filtering power can be set to the filtering direction
Weight, and for other filtering directions in multiple filtering direction, some less weighted values can be set.
I.e.：The weighted value moved towards on the closest filtering direction in direction with high frequency trajectory can be more than other
Weighted value on one filtering direction.
In sub-step 203, CT data for projection is filtered on multiple filtering direction according to weighted value
Ripple.
In this sub-step, the weighted value that sub-step 202 is set can be applied to corresponding filtering
It is filtered on direction.
In another embodiment of the present invention, with reference to Fig. 3, step 102 can also include following sub-step
Rapid 301 to 303.
In sub-step 301, selected on 360 degree of directions on the sinogram being made up of CT data for projection
Take multiple filtering directions.
For example, on 0 to 360 direction, 0 degree, 45 degree, the multiple directions such as 90 degree can be chosen
As filtering direction.
In sub-step 302, in multiple filtering directions, the vertical line side pair with the tangent line of high frequency trajectory
Greater weight value is set and to other filtering directions in multiple filtering direction to immediate filtering direction
Smaller weighted value is set.
In the multiple filtering directions selected by sub-step 301, at least there is a filtering direction, the filter
Ripple direction is closest with the vertical line direction of the tangent line of high frequency trajectory or even overlaps, then can be to the filtering
Direction sets larger filtering weighting, and for other filtering directions in multiple filtering direction, Ke Yishe
Put some less weighted values.I.e.：Power on the filtering direction closest with above-mentioned vertical line direction
Weight values can be more than the weighted value on other any filtering directions.
In sub-step 303, CT data for projection is filtered on multiple filtering direction according to weighted value
Ripple.
In this sub-step, the weighted value that sub-step 302 is set can be applied to corresponding filtering
It is filtered on direction.
So far the method being filtered to CT data for projection according to embodiments of the present invention is described.Due to
The inventive method is to be based on the trend of high frequency trajectory therein and enter when being filtered to CT data for projection
Capable, therefore, the inventive method can preferably keep the high fdrequency component in CT data for projection and balance
Noise, can reduce the granularity of CT images and reduce the artifacts such as smudge, center dot.
It is similar with the method, present invention also offers corresponding device.
Fig. 4 show showing for one embodiment of the device being filtered to CT data for projection of the invention
Meaning property block diagram.
As shown in figure 4, device 400 can include：High frequency trajectory acquisition module 401, for obtaining CT
High frequency trajectory in data for projection；And high frequency trajectory filtration module, for the trend based on high frequency trajectory
CT data for projection is filtered.
In one embodiment of the invention, high frequency trajectory filtration module 402 may further include：Edge
Filtration module is moved towards in track, and CT data for projection is filtered for the trend along high frequency trajectory.
In one embodiment of the invention, high frequency trajectory filtration module 402 may further include：Edge
Vertical line trend pass filtering module, for projecting number to CT along the direction perpendicular with the tangent line of high frequency trajectory
According to being filtered.
In one embodiment of the invention, high frequency trajectory filtration module 402 may further include：It is many
Module is chosen in filtering direction, for 360 degree of directions on the sinogram being made up of CT data for projection
Choose multiple filtering directions；First weight setting module, for multiple filtering directions in, pair with it is vertical
Maximum weighted value is set in the immediate filtering direction in the direction of the tangent line of high frequency trajectory and multiple is filtered
Other filtering directions in ripple direction set smaller weighted value；And many filtering trend pass filtering modules, it is used for
CT data for projection is filtered on multiple filtering direction according to weighted value.
In one embodiment of the invention, high frequency trajectory filtration module 402 may further include：It is many
Module is chosen in filtering direction, for 360 degree of directions on the sinogram being made up of CT data for projection
Choose multiple filtering directions；Second weight setting module, for multiple filtering directions in, pair and high frequency
The immediate filtering direction in vertical line direction of the tangent line of track sets greater weight value and multiple is filtered
Other filtering directions in direction set smaller weighted value；And many filtering trend pass filtering modules, for root
CT data for projection is filtered on multiple filtering direction according to weighted value.
So far the device being filtered to CT data for projection according to embodiments of the present invention is described.With it is upper
State method to be similar to, be based on height therein when being filtered to CT data for projection due to apparatus of the present invention
The trend of frequency track and carry out, therefore, apparatus of the present invention can be preferably kept in CT data for projection
High fdrequency component and balances noise, the granularity of CT images can be reduced and smudge, center dot is reduced
Etc. artifact.
Embodiments of the invention are the foregoing is only, is not intended to limit the invention, for this area
Technical staff for, the present invention can have various modifications and variations.It is all in the spirit and principles in the present invention
Within, any modification, equivalent substitution and improvements made etc. should be included in claim of the invention
Within the scope of.
Claims (10)
1. a kind of method being filtered to CT data for projection, it is characterized in that, including：
Obtain the high frequency trajectory in the CT data for projection；And
Trend based on the high frequency trajectory is filtered to the CT data for projection.
2. method according to claim 1, it is characterized in that, it is described based on the high frequency trajectory
The step of trend is filtered to the CT data for projection further includes：
The CT data for projection is filtered along the trend of the high frequency trajectory.
3. method according to claim 1, it is characterized in that, it is described based on the high frequency trajectory
The step of trend is filtered to the CT data for projection further includes：
The CT data for projection is carried out along the direction perpendicular with the tangent line of the high frequency trajectory
Filtering.
4. method according to claim 1, it is characterized in that, it is described based on the high frequency trajectory
The step of trend is filtered to the CT data for projection further includes：
Multiple filters are chosen on 360 degree of directions on the sinogram being made up of the CT data for projection
Ripple direction；
In the multiple filtering direction, pair immediate filtering side is moved towards with the high frequency trajectory
Set to the larger weighted value of setting and to other filtering directions in the multiple filtering direction
Put smaller weighted value；And
The CT data for projection is filtered on the multiple filtering direction according to the weighted value
Ripple.
5. method according to claim 1, it is characterized in that, it is described based on the high frequency trajectory
The step of trend is filtered to the CT data for projection further includes：
Multiple filters are chosen on 360 degree of directions on the sinogram being made up of the CT data for projection
Ripple direction；
In the multiple filtering direction, pair with the vertical line direction of the tangent line of the high frequency trajectory most
Close filtering direction sets greater weight value and to other in the multiple filtering direction
Filtering direction sets smaller weighted value；And
The CT data for projection is filtered on the multiple filtering direction according to the weighted value
Ripple.
6. a kind of device being filtered to CT data for projection, it is characterized in that, including：
High frequency trajectory acquisition module, for obtaining the high frequency trajectory in the CT data for projection；With
And
High frequency trajectory filtration module, projects for the trend based on the high frequency trajectory to the CT
Data are filtered.
7. device according to claim 6, it is characterized in that, the high frequency trajectory filtration module enters
One step includes：
Filtration module is moved towards along track, the CT is thrown for the trend along the high frequency trajectory
Shadow data are filtered.
8. device according to claim 6, it is characterized in that, the high frequency trajectory filtration module enters
One step includes：
Along vertical line trend pass filtering module, for along the side perpendicular with the tangent line of the high frequency trajectory
It is filtered to the CT data for projection.
9. device according to claim 6, it is characterized in that, the high frequency trajectory filtration module enters
One step includes：
Module is chosen in many filtering directions, in the sinogram being made up of the CT data for projection
360 degree of directions on choose multiple filtering directions；
First weight setting module, in the multiple filtering direction, pair with perpendicular to described
The weighted value of the immediate filtering direction setting maximum in direction of the tangent line of high frequency trajectory is simultaneously right
Other filtering directions in the multiple filtering direction set smaller weighted value；And
Many filtering trend pass filtering modules, for according to the weighted value in the multiple filtering direction
The CT data for projection is filtered.
10. device according to claim 6, it is characterized in that, the high frequency trajectory filtration module enters
One step includes：
Module is chosen in many filtering directions, in the sinogram being made up of the CT data for projection
360 degree of directions on choose multiple filtering directions；
Second weight setting module, in the multiple filtering direction, pair with the high frequency rail
The immediate filtering direction in vertical line direction of the tangent line of mark sets greater weight value and to institute
Other filtering directions stated in multiple filtering directions set smaller weighted value；And
Many filtering trend pass filtering modules, for according to the weighted value in the multiple filtering direction
The CT data for projection is filtered.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510979443.4A CN106910164A (en) | 2015-12-23 | 2015-12-23 | A kind of method and device being filtered to CT data for projection |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510979443.4A CN106910164A (en) | 2015-12-23 | 2015-12-23 | A kind of method and device being filtered to CT data for projection |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106910164A true CN106910164A (en) | 2017-06-30 |
Family
ID=59200514
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510979443.4A Pending CN106910164A (en) | 2015-12-23 | 2015-12-23 | A kind of method and device being filtered to CT data for projection |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106910164A (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1207656A (en) * | 1996-01-17 | 1999-02-10 | 模拟技术公司 | Computed tomography streak suppression filter |
| CN103559729A (en) * | 2013-11-18 | 2014-02-05 | 首都师范大学 | Method for iterating and reconstructing double-energy-spectrum CT image |
| CN103645502A (en) * | 2013-12-11 | 2014-03-19 | 中国海洋石油总公司 | Seismic wave attenuation compensation method in curvelet domain |
| CN104867157A (en) * | 2015-06-01 | 2015-08-26 | 华润万东医疗装备股份有限公司 | Correction method of pixel response inconsistence of CT detector |
-
2015
- 2015-12-23 CN CN201510979443.4A patent/CN106910164A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1207656A (en) * | 1996-01-17 | 1999-02-10 | 模拟技术公司 | Computed tomography streak suppression filter |
| CN103559729A (en) * | 2013-11-18 | 2014-02-05 | 首都师范大学 | Method for iterating and reconstructing double-energy-spectrum CT image |
| CN103645502A (en) * | 2013-12-11 | 2014-03-19 | 中国海洋石油总公司 | Seismic wave attenuation compensation method in curvelet domain |
| CN104867157A (en) * | 2015-06-01 | 2015-08-26 | 华润万东医疗装备股份有限公司 | Correction method of pixel response inconsistence of CT detector |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103310432B (en) | Based on the metal artifacts reduction method of the CT image normalization of the full variation stream of quadravalence | |
| CN104992409B (en) | A kind of metal artifacts reduction method of CT images | |
| US8897527B2 (en) | Motion-blurred imaging enhancement method and system | |
| US8903150B2 (en) | Filtration imaging enhancement method and system | |
| CN103971387B (en) | CT image rebuilding methods | |
| CN106296615B (en) | The method and system of metal artifacts is corrected in CT images | |
| Naranjo et al. | Metal artifact reduction in dental CT images using polar mathematical morphology | |
| Rouf et al. | Glare encoding of high dynamic range images | |
| CN103116879A (en) | Neighborhood windowing based non-local mean value CT (Computed Tomography) imaging de-noising method | |
| CN106023204A (en) | Method and system for removing mosquito noise based on edge detection algorithm | |
| CN103186888A (en) | Method and device for removing CT (computed tomography) image noises | |
| CN104751432A (en) | Image reconstruction based visible light and infrared image fusion method | |
| CN103810735A (en) | Statistical iterative reconstructing method for low-dose X-ray CT image | |
| CN103810733A (en) | Statistical iteration reconstruction method of sparse-angle X-ray CT (electronic computer X-ray tomography technique) image | |
| CN107403419A (en) | A kind of low dose X-ray image de-noising method based on concatenated convolutional neutral net | |
| CN103106675A (en) | X-ray computed tomography (CT) metal artifact processing method | |
| CN103971349B (en) | Computed tomography images method for reconstructing and ct apparatus | |
| US8369643B2 (en) | Image processing method, image processing apparatus, and image processing program | |
| CN110070510A (en) | A kind of CNN medical image denoising method for extracting feature based on VGG-19 | |
| CN107194902A (en) | The method that wave filter parameter is automatically determined in image co-registration | |
| CN106910164A (en) | A kind of method and device being filtered to CT data for projection | |
| CN105184741A (en) | Three-dimensional CBCT (cone-beam computed tomography) image denoising method on the basis of improved nonlocal means | |
| DE112014005888T5 (en) | Method and apparatus for domain-specific MR-based attenuation correction of time domain filtered PET studies | |
| DE102011017710A1 (en) | A method for providing a 3D X-ray image data set to a moving object with highly absorbent material | |
| CN107274409A (en) | woven fabric defect segmentation method |
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 |