CN103876772A - Multi-spectrum imaging method and device - Google Patents
Multi-spectrum imaging method and device Download PDFInfo
- Publication number
- CN103876772A CN103876772A CN201410105317.1A CN201410105317A CN103876772A CN 103876772 A CN103876772 A CN 103876772A CN 201410105317 A CN201410105317 A CN 201410105317A CN 103876772 A CN103876772 A CN 103876772A
- Authority
- CN
- China
- Prior art keywords
- component
- energy
- emax
- sigma
- emin
- 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
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
The invention discloses a multi-spectrum imaging method and device. The multi-spectrum imaging method includes the steps of determining an energy spectrum filtering range and an energy spectrum separating time number of a detected object, carrying out energy-variable imaging on the detected object to obtain a gradient energy projection sequence of energy spectrum separating, carrying out multi-spectrum CT reconstruction according to the gradient energy projection sequence to obtain attenuation images under different energies, and carrying out multi-spectrum fusing on the obtained attenuation images to obtain a required true-color CT image. By means of the scheme, the application range of a CT system can be expanded.
Description
Technical field
The present invention relates to computed tomography (CT, Computer Tomography) technology, particularly a kind of multispectral formation method and device.
Background technology
X ray CT imaging technique has been widely used in the fields such as medical science, industry, safety detection at present, and has brought into play the effect can not be substituted.
But along with the development of technology, be applicable to the existing Single energy X ray absorptionmetry CT imaging technique of flaw inspection, structural analysis etc., can not have met the functional imaging demand that in modern industry and modern medicine, component quantification detects.This is because existing CT system is all to realize the structure imaging of detected object based on monoenergetic hypothesis, cannot realize multispectral quantitative analysis and functional imaging, cannot realize multispectral imaging, thereby limit the scope of application of CT system.
Summary of the invention
In view of this, the invention provides a kind of multispectral formation method and device, can expand the scope of application of CT system.
In order to achieve the above object, technical scheme of the present invention is achieved in that
A kind of multispectral formation method, comprising:
The power spectrum filter range of determining detected object separates number of times with power spectrum;
According to definite result, described detected object is become to energy imaging, obtain the alternation energy projection sequence that power spectrum separates;
Carry out multispectral computer tomography CT reconstruction according to described alternation energy projection sequence, obtain the decay pattern picture under different-energy;
The each decay pattern obtaining is looked like to carry out multispectral fusion, obtain required very color CT image.
A kind of multispectral imaging device, comprising:
Ray filtering and separation module, separate number of times for the power spectrum filter range of determining detected object with power spectrum, and will determine that result sends to alternation energy imaging module;
Described alternation energy imaging module, for according to determining result, becomes energy imaging to described detected object, obtains the alternation energy projection sequence that power spectrum separates, and sends to multispectral computer tomography CT to rebuild module described alternation energy projection sequence;
Described multispectral CT rebuilds module, for carrying out multispectral CT reconstruction according to described alternation energy projection sequence, obtains the decay pattern picture under different-energy, and the each decay pattern obtaining is looked like to send to multispectral Fusion Module;
Described multispectral Fusion Module, for the each decay pattern obtaining is looked like to carry out multispectral fusion, obtains required very color CT image.
Visible, adopt scheme of the present invention, can be based on existing CT system, utilize ray filtering and become the technology such as energy imaging, realizing ray emission end power spectrum and separate, on this basis, can further utilize the technology such as multispectral CT reconstruction and multispectral fusion, obtain very color CT image, thereby realized multispectral imaging, and then expanded the scope of application of CT system.
Brief description of the drawings
Fig. 1 is the flow chart of the multispectral formation method embodiment of the present invention.
Fig. 2 is the composition structural representation of the multispectral imaging device embodiment of the present invention.
Detailed description of the invention
For make technical scheme of the present invention clearer, understand, referring to the accompanying drawing embodiment that develops simultaneously, scheme of the present invention is described in further detail.
Fig. 1 is the flow chart of the multispectral formation method embodiment of the present invention.As shown in Figure 1, comprise the following steps 11~14.
Step 11: the power spectrum filter range of determining detected object separates number of times with power spectrum.
Realize very color CT imaging, need to guarantee that the data for projection under each energy is monoenergetic or approximate monoenergetic data, and in existing CT system, due to bremsstrahlung, what obtain is projected as wide range projection, therefore cannot further realize power spectrum and separate.
For this reason, in this step, the first prior information such as component and overall dimensions based on detected object, the power spectrum filter range of determining detected object separates number of times with power spectrum, thereby provides energy parameter that foundation is set for follow-up obtaining of alternation energy projection sequence.
In addition, in x-ray imaging process, by the filter plate at the additional different materials of ray emission end and thickness, as aluminum, steel and tantalum etc., can make wider continuous spectrum become the narrow spectrum of approximate monoenergetic, thereby can overcome sclerosis artifact that continuous spectrum imaging occurs etc.; Meanwhile, the attenuation quotient of different materials, along with the rate of change difference of ray energy, under the condition of monoenergetic imaging, can reach the object that detected object component is distinguished according to the rate of change difference of attenuation quotient, thereby meet the demand of multispectral imaging; But, can increase the noise level of projection to ray filtering (energy spectral filter), thereby affect projection quality; , in scheme of the present invention, need to meet under the prerequisite of component discrimination for this reason, reduce as much as possible ray filtering degree, thereby ensure projection quality.
Correspondingly, in scheme of the present invention, for every kind of component i in detected object, can be respectively according to step 1)~3) shown in mode process.
1) according to the attenuation quotient rate of change model of component i, determine the ceiling capacity Emax that component i is corresponding;
The specific implementation of this step is prior art.
2) determine the span of the least energy Emin that this Emax is corresponding according to predetermined relationship;
Described predetermined relationship is:
Wherein, S (E) is illustrated in after energy spectral filter, and the number of photons that energy is E accounts for ray through the number of photons I before detected object
0ratio; I represents the component species number that detected object comprises, m
i(E) represent the mass attentuation coefficient of component i under energy E, ρ
irepresent the density of component i, d
ithe thickness that represents component i, χ represents the Limit of J-validity that detector can receive, and the value of above each parameter is all known, and Δ E represents the value interval of E, if value is 1;
In actual applications, need to ensure that in the span [Emin, Emax] of energy, the data that detector receives are effective, there will not be overexposure and under exposed situation after energy spectral filter.
3) determine respectively whether the each Emin that is positioned at span meets the following conditions:
If so, retain this Emin, otherwise, this Emin abandoned;
In order to ensure projection quality, using retain each Emin in value minimum as Emin corresponding to component i.
By the way, obtain the power spectrum filter range of detected object, afterwards, can add up the different value numbers that occur in Emax corresponding to each component, and statistical result is separated to number of times as power spectrum, i.e. power spectrum separation number of times can be less than or equal to the component species number that detected object comprises.
Illustrate: the value of supposing I is 5, in detected object, comprise altogether 5 component kinds, if Emax corresponding to various component is respectively: 10,15,20,25,30, power spectrum separation number of times is 5 so; If Emax corresponding to various components is respectively: 10,15,20,25,10, power spectrum separation number of times is 4 so.
In addition, in the time that Emax corresponding to N component is identical, N is greater than 1 positive integer, can using value minimum in Emax corresponding this N component as this N component and Emin corresponding to this Emax.
And, for the Emax of each different values, can, according to monoenergetic approximate condition, determine material and the thickness of the filter plate that this Emax is corresponding according to (its corresponding component is corresponding) Emin of its correspondence respectively, specific implementation is prior art.
Step 12: according to definite result, detected object is become to energy imaging, obtain the alternation energy projection sequence that power spectrum separates.
In this step, can separate number of times with power spectrum according to the power spectrum filter range of determining in step 11, carry out CT scan imaging, the alternation energy projection sequence separating to obtain power spectrum.
In the time carrying out CT scan, can adopt traditional CT circle track scanning pattern, angle stepping is 1 degree, and angular range is 0~359 degree, and difference is that energy parameter changes in CT scan process.
Specifically, can be according to order from small to large, Emax to variant value sorts, and according to clooating sequence, carry out the CT scan of alternation energy according to the Emax of variant value, separate number of times according to power spectrum and determine the number of times of alternation energy, and can arrange by the tube voltage parameter of adjusting CT scan, realize the change of energy; In addition, also need to add respectively at ray emission window place the filter plate corresponding with the Emax of variant value.
Step 13: carry out multispectral CT reconstruction according to alternation energy projection sequence, obtain the decay pattern picture under different-energy.
In this step, the two-dimentional alternation energy projection sequence getting is carried out to three-dimensional reconstruction, to obtain the internal data structure of detected object under each energy in step 12.Preferably, can adopt the multispectral CT algorithm for reconstructing of low dosage, alternation energy projection sequence is carried out to multispectral CT reconstruction, to obtain the image that under each energy, quality better, component discrimination is higher.
In CT system, it is a stochastic process that the factors such as the scattering of X ray, the fluctuation noise of detector make the acquisition process of data for projection, the research of nuclear physics simultaneously also shows already, the radiation of X-ray meets Poisson random distribution, therefore, can utilize stochastic signal processing method, set up the statistical model of CT imaging process, and the theory that application parameter is estimated is carried out image reconstruction.
Visit first reading d (y) for y for detector under each energy, suppose that it is the Poisson random distribution of g (y, μ) that its obedience is expected:
Wherein, I
0(y, E) is illustrated in while thering is no detected object, visits the number of photons that first energy receiving is E, a for y
xyrepresent that y is visited the length of ray corresponding to unit through detected object x position, μ (x, E) represents line attenuation coefficient, relevant with energy and position; X represents to rebuild x pixel position in the decay pattern picture that obtains, in scheme of the present invention, in the time of reconstruction by discrete the cross section of detected object be digital picture;
Owing to there being the prior information of component, μ (x, E) can be expressed as the linear combination of I kind component, that is:
Wherein, μ
i(E) represent the line attenuation coefficient of component i under energy E, c
i(x) represent that component i is at the shared proportion in x position;
According to the character of Poisson random distribution, the logarithm maximum likelihood function of the first reading of all spies is:
Correspondingly, what in scheme of the present invention, will do is exactly, and is meeting under formula (3) and (4) and non-negative constraint, obtains the c that makes formula (5) maximum
i(x).
Solve this restricted problem, can obtain reconstruction formula as follows:
Wherein,
K represents iterations, as 1,2,3 etc.; The span of i is from 1~I, and the span of E is from Emin~Emax, and Emin and Emax are Emin and the Emax that component i is corresponding, and the span of y is 1~M, and M represents the spy unit number that detector comprises; In addition, above-mentioned Z
i(x) be only preferably definition mode of one, also can adopt according to actual needs other definition mode, be not restricted; Above-mentioned μ
i(E), a
xy, I
0the isoparametric value of (y, E), d (y) is all known.
Correspondingly, in scheme of the present invention, can, for each energy, determine respectively decay pattern picture corresponding to every kind of component i under this energy;
Wherein, the mode of determining the decay pattern picture that every kind of component i under this energy is corresponding can be:
Carry out following calculating according to the mode of iteration:
When often obtaining a c
i k+1(x) afterwards, determine respectively c
i k+1and c (x)
i k(x) whether difference is less than predetermined threshold, if so, and by c
i k+1(x) as final result corresponding to x; The concrete value of described threshold value can be decided according to the actual requirements; Wherein, c
i k(x) initial value, value when k=1 can preset;
According to final result corresponding to variant x difference obtaining, generate decay pattern picture corresponding to component i, how to be generated as prior art.
Step 14: the each decay pattern obtaining is looked like to carry out multispectral fusion, obtain required very color CT image.
In this step, the each decay pattern getting in step 13 is looked like to carry out fusion treatment, to obtain can specification configuration distinguishing again the very color CT image of material.
Specifically, first, can utilize the multispectral blending algorithm based on principal component analysis, extract respectively relatively independent main constituent image in every width decay pattern picture, afterwards, can be according to extracting result, by adaptive weighted fusion, obtain fusion image, last, can be painted by fusion image is carried out, obtain required very color CT image; Specific implementation is prior art.
Based on above-mentioned introduction, Fig. 2 is the composition structural representation of the multispectral imaging device embodiment of the present invention.As shown in Figure 2, comprising:
Ray filtering and separation module, separate number of times for the power spectrum filter range of determining detected object with power spectrum, and will determine that result sends to alternation energy imaging module;
Alternation energy imaging module, for according to determining result, becomes energy imaging to detected object, obtains the alternation energy projection sequence that power spectrum separates, and sends to multispectral computer tomography CT to rebuild module alternation energy projection sequence;
Multispectral CT rebuilds module, for carrying out multispectral CT reconstruction according to alternation energy projection sequence, obtains the decay pattern picture under different-energy, and the each decay pattern obtaining is looked like to send to multispectral Fusion Module;
Multispectral Fusion Module, for the each decay pattern obtaining is looked like to carry out multispectral fusion, obtains required very color CT image.
Preferably,
Ray filtering and separation module can, for every kind of component i in detected object, according to the attenuation quotient rate of change model of component i, be determined the ceiling capacity Emax that component i is corresponding respectively; Determine the span of the least energy Emin that this Emax is corresponding according to predetermined relationship; Determine respectively whether the each Emin that is positioned at span meets predetermined condition; If so, retain this Emin, otherwise, this Emin abandoned; Using retain each Emin in value minimum as Emin corresponding to component i;
Wherein, described predetermined relationship is:
Described predetermined condition is:
S (E) is illustrated in after energy spectral filter, and the number of photons that energy is E accounts for ray through the number of photons I before detected object
0ratio; I represents the component species number that detected object comprises, m
i(E) represent the mass attentuation coefficient of component i under energy E, ρ
irepresent the density of component i, d
ithe thickness that represents component i, χ represents the Limit of J-validity that detector can receive, Δ E represents the value interval of E;
Ray filtering can add up from separation module the different value numbers that occur in Emax corresponding to each component, and statistical result is separated to number of times as power spectrum.
In addition,
Ray filtering and separation module also can be further used for, and in the time that Emax corresponding to N component is identical, N is greater than 1 positive integer, using value minimum in Emax corresponding this N component as this N component and Emin corresponding to this Emax; For the Emax of each different values, determine respectively material and the thickness of the filter plate that this Emax is corresponding according to its corresponding Emin, and notice is to alternation energy imaging module;
Correspondingly, alternation energy imaging module can, according to order from small to large, sort to the Emax of variant value; According to clooating sequence, carry out the CT scan of alternation energy according to the Emax of variant value; And, add respectively the filter plate corresponding with the Emax of variant value at ray emission window place.
Preferably,
Multispectral CT rebuilds module can, for each energy, determine respectively decay pattern picture corresponding to every kind of component i under this energy; Wherein, for every kind of component i under this energy, carry out following calculating according to the mode of iteration respectively:
When often obtaining a c
i k+1(x) afterwards, determine respectively c
i k+1and c (x)
i k(x) whether difference is less than predetermined threshold, if so, and by c
i k+1(x) as final result corresponding to x; According to final result corresponding to variant x difference obtaining, generate decay pattern picture corresponding to component i;
Wherein, μ
i(E) represent the line attenuation coefficient of component i under energy E, a
xyrepresent that y in detector is visited the length of ray corresponding to unit through detected object x position, I
0(y, E) is illustrated in while thering is no detected object, visits the number of photons that first energy receiving is E, c for y in detector
i(x) represent that component i is at the shared proportion in x position, d (y) represents that y in detector visits first reading, and x represents to rebuild x pixel position in the decay pattern picture obtaining, and k represents iterations; The span of i is from 1~I, and the span of E is from Emin~Emax, and Emin and Emax are Emin and the Emax that component i is corresponding, and the span of y is 1~M, and M represents the spy unit number that detector comprises.
Preferably,
The multispectral blending algorithm of multispectral Fusion Module utilization based on principal component analysis, extracts respectively relatively independent main constituent image in every width decay pattern picture; According to extracting result, by adaptive weighted fusion, obtain fusion image; Painted by fusion image is carried out, obtain very color CT image.
The specific works flow process of Fig. 2 shown device embodiment please refer to the respective description in preceding method embodiment, repeats no more herein.In addition, in actual applications, Fig. 2 shown device can be applicable in CT system.
In sum, these are only preferred embodiment of the present invention, be not intended to limit protection scope of the present invention.Within the spirit and principles in the present invention all, any amendment of doing, be equal to replacement, improvement etc., within all should being included in protection scope of the present invention.
Claims (10)
1. a multispectral formation method, is characterized in that, comprising:
The power spectrum filter range of determining detected object separates number of times with power spectrum;
According to definite result, described detected object is become to energy imaging, obtain the alternation energy projection sequence that power spectrum separates;
Carry out multispectral computer tomography CT reconstruction according to described alternation energy projection sequence, obtain the decay pattern picture under different-energy;
The each decay pattern obtaining is looked like to carry out multispectral fusion, obtain required very color CT image.
2. method according to claim 1, is characterized in that,
The power spectrum filter range of described definite detected object comprises:
For every kind of component i in described detected object, carry out respectively following processing:
According to the attenuation quotient rate of change model of component i, determine the ceiling capacity Emax that component i is corresponding;
Determine the span of the least energy Emin that this Emax is corresponding according to predetermined relationship;
Described predetermined relationship is:
Wherein, S (E) is illustrated in after energy spectral filter, and the number of photons that energy is E accounts for ray through the number of photons I before described detected object
0ratio; I represents the component species number that described detected object comprises, m
i(E) represent the mass attentuation coefficient of component i under energy E, ρ
irepresent the density of component i, d
ithe thickness that represents component i, χ represents the Limit of J-validity that detector can receive, Δ E represents the value interval of E;
Determine respectively whether the each Emin that is positioned at span meets the following conditions:
If so, retain this Emin, otherwise, this Emin abandoned;
Using retain each Emin in value minimum as Emin corresponding to component i;
The power spectrum of described definite detected object separates number of times and comprises: add up the different value numbers that occur in Emax corresponding to each component, statistical result is separated to number of times as described power spectrum.
3. method according to claim 2, is characterized in that,
The method further comprises: in the time that Emax corresponding to N component is identical, N is greater than 1 positive integer, using value minimum in Emax corresponding this N component as this N component and Emin corresponding to this Emax;
For the Emax of each different values, determine respectively material and the thickness of the filter plate that this Emax is corresponding according to its corresponding Emin;
Described described detected object is become to energy imaging, obtains the alternation energy projection sequence that power spectrum separates and comprise:
According to order from small to large, the Emax of variant value is sorted;
According to clooating sequence, carry out the CT scan of alternation energy according to the Emax of variant value; And, add respectively the filter plate corresponding with the Emax of variant value at ray emission window place.
4. method according to claim 3, is characterized in that,
Describedly carry out multispectral CT reconstruction according to described alternation energy projection sequence, the decay pattern obtaining under different-energy looks like to comprise:
For each energy, determine respectively decay pattern picture corresponding to every kind of component i under this energy;
Wherein, every kind of decay pattern corresponding to component i under described definite this energy looks like to comprise:
Carry out following calculating according to the mode of iteration:
Wherein, μ
i(E) represent the line attenuation coefficient of component i under energy E, a
xyrepresent that y in detector is visited the length of ray corresponding to unit through described detected object x position, I
0(y, E) is illustrated in while thering is no described detected object, visits the number of photons that first energy receiving is E, c for y in detector
i(x) represent that component i is at the shared proportion in x position, d (y) represents that y in detector visits first reading, and x represents to rebuild x pixel position in the decay pattern picture obtaining, and k represents iterations; The span of i is from 1~I, and the span of E is from Emin~Emax, and Emin and Emax are Emin and the Emax that component i is corresponding, and the span of y is 1~M, and M represents the spy unit number that detector comprises;
When often obtaining a c
i k+1(x) afterwards, determine respectively c
i k+1and c (x)
i k(x) whether difference is less than predetermined threshold, if so, and by c
i k+1(x) as final result corresponding to x;
According to final result corresponding to variant x difference obtaining, generate decay pattern picture corresponding to component i.
5. according to the method described in claim 1,2,3 or 4, it is characterized in that,
Described the each decay pattern obtaining is looked like to carry out multispectral fusion, obtains required very color CT image and comprise:
Utilize the multispectral blending algorithm based on principal component analysis, extract respectively relatively independent main constituent image in every width decay pattern picture;
According to extracting result, by adaptive weighted fusion, obtain fusion image;
Painted by described fusion image is carried out, obtain described very color CT image.
6. a multispectral imaging device, is characterized in that, comprising:
Ray filtering and separation module, separate number of times for the power spectrum filter range of determining detected object with power spectrum, and will determine that result sends to alternation energy imaging module;
Described alternation energy imaging module, for according to determining result, becomes energy imaging to described detected object, obtains the alternation energy projection sequence that power spectrum separates, and sends to multispectral computer tomography CT to rebuild module described alternation energy projection sequence;
Described multispectral CT rebuilds module, for carrying out multispectral CT reconstruction according to described alternation energy projection sequence, obtains the decay pattern picture under different-energy, and the each decay pattern obtaining is looked like to send to multispectral Fusion Module;
Described multispectral Fusion Module, for the each decay pattern obtaining is looked like to carry out multispectral fusion, obtains required very color CT image.
7. device according to claim 6, is characterized in that,
Described ray filtering and separation module, for every kind of component i in described detected object, according to the attenuation quotient rate of change model of component i, are determined the ceiling capacity Emax that component i is corresponding respectively; Determine the span of the least energy Emin that this Emax is corresponding according to predetermined relationship; Determine respectively whether the each Emin that is positioned at span meets predetermined condition; If so, retain this Emin, otherwise, this Emin abandoned; Using retain each Emin in value minimum as Emin corresponding to component i;
Wherein, described predetermined relationship is:
Described predetermined condition is:
S (E) is illustrated in after energy spectral filter, and the number of photons that energy is E accounts for ray through the number of photons I before described detected object
0ratio; I represents the component species number that described detected object comprises, m
i(E) represent the mass attentuation coefficient of component i under energy E, ρ
irepresent the density of component i, d
ithe thickness that represents component i, χ represents the Limit of J-validity that detector can receive, Δ E represents the value interval of E;
Described ray filtering adds up from separation module the different value numbers that occur in Emax corresponding to each component, and statistical result is separated to number of times as described power spectrum.
8. device according to claim 7, is characterized in that,
Described ray filtering and separation module are further used for, and in the time that Emax corresponding to N component is identical, N is greater than 1 positive integer, using value minimum in Emax corresponding this N component as this N component and Emin corresponding to this Emax; For the Emax of each different values, determine respectively material and the thickness of the filter plate that this Emax is corresponding according to its corresponding Emin, and notice is to described alternation energy imaging module;
Described alternation energy imaging module, according to order from small to large, sorts to the Emax of variant value; According to clooating sequence, carry out the CT scan of alternation energy according to the Emax of variant value; And, add respectively the filter plate corresponding with the Emax of variant value at ray emission window place.
9. device according to claim 8, is characterized in that,
Described multispectral CT rebuilds module for each energy, determines respectively decay pattern picture corresponding to every kind of component i under this energy; Wherein, for every kind of component i under this energy, carry out following calculating according to the mode of iteration respectively:
When often obtaining a c
i k+1(x) afterwards, determine respectively c
i k+1and c (x)
i k(x) whether difference is less than predetermined threshold, if so, and by c
i k+1(x) as final result corresponding to x; According to final result corresponding to variant x difference obtaining, generate decay pattern picture corresponding to component i;
Wherein, μ
i(E) represent the line attenuation coefficient of component i under energy E, a
xyrepresent that y in detector is visited the length of ray corresponding to unit through described detected object x position, I
0(y, E) is illustrated in while thering is no described detected object, visits the number of photons that first energy receiving is E, c for y in detector
i(x) represent that component i is at the shared proportion in x position, d (y) represents that y in detector visits first reading, and x represents to rebuild x pixel position in the decay pattern picture obtaining, and k represents iterations; The span of i is from 1~I, and the span of E is from Emin~Emax, and Emin and Emax are Emin and the Emax that component i is corresponding, and the span of y is 1~M, and M represents the spy unit number that detector comprises.
10. according to the device described in claim 6,7,8 or 9, it is characterized in that,
The multispectral blending algorithm of described multispectral Fusion Module utilization based on principal component analysis, extracts respectively relatively independent main constituent image in every width decay pattern picture; According to extracting result, by adaptive weighted fusion, obtain fusion image; Painted by described fusion image is carried out, obtain described very color CT image.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410105317.1A CN103876772B (en) | 2014-03-20 | 2014-03-20 | A kind of multispectral formation method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410105317.1A CN103876772B (en) | 2014-03-20 | 2014-03-20 | A kind of multispectral formation method and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103876772A true CN103876772A (en) | 2014-06-25 |
CN103876772B CN103876772B (en) | 2015-12-09 |
Family
ID=50946131
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410105317.1A Active CN103876772B (en) | 2014-03-20 | 2014-03-20 | A kind of multispectral formation method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103876772B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106534866A (en) * | 2016-11-25 | 2017-03-22 | 中北大学 | Image compression method and apparatus |
CN109813259A (en) * | 2019-01-18 | 2019-05-28 | 中北大学 | Merge CT imaging method, storage medium and device |
CN109875593A (en) * | 2019-02-02 | 2019-06-14 | 中北大学 | CT imaging method, storage medium and device |
CN111476856A (en) * | 2020-04-08 | 2020-07-31 | 中北大学 | Multispectral CT imaging method |
WO2021004157A1 (en) * | 2019-07-09 | 2021-01-14 | 上海联影医疗科技有限公司 | Medical scanning imaging method and apparatus, storage medium, and computer device |
CN112603345A (en) * | 2020-12-02 | 2021-04-06 | 赛诺威盛科技(北京)有限公司 | Model training method, multi-energy spectrum CT scanning method, device and electronic equipment |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030053597A1 (en) * | 2000-09-29 | 2003-03-20 | Thomas Flohr | X-ray computer tomograph |
CN101074937A (en) * | 2006-05-19 | 2007-11-21 | 清华大学 | Energy spectrum modulator, method and apparatus for discriminating material and image processing method |
CN101308102A (en) * | 2008-07-16 | 2008-11-19 | 中北大学 | Computer tomography scanned imagery apparatus and method |
CN101732057A (en) * | 2008-11-26 | 2010-06-16 | 通用电气公司 | Systems and methods for displaying multi-energy data |
WO2010070554A1 (en) * | 2008-12-17 | 2010-06-24 | Koninklijke Philips Electronics N.V. | X-ray examination apparatus and method |
CN103308537A (en) * | 2013-06-13 | 2013-09-18 | 中北大学 | Gradient-energy X-ray imaging image fusion method |
-
2014
- 2014-03-20 CN CN201410105317.1A patent/CN103876772B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030053597A1 (en) * | 2000-09-29 | 2003-03-20 | Thomas Flohr | X-ray computer tomograph |
CN101074937A (en) * | 2006-05-19 | 2007-11-21 | 清华大学 | Energy spectrum modulator, method and apparatus for discriminating material and image processing method |
CN101308102A (en) * | 2008-07-16 | 2008-11-19 | 中北大学 | Computer tomography scanned imagery apparatus and method |
CN101732057A (en) * | 2008-11-26 | 2010-06-16 | 通用电气公司 | Systems and methods for displaying multi-energy data |
WO2010070554A1 (en) * | 2008-12-17 | 2010-06-24 | Koninklijke Philips Electronics N.V. | X-ray examination apparatus and method |
CN103308537A (en) * | 2013-06-13 | 2013-09-18 | 中北大学 | Gradient-energy X-ray imaging image fusion method |
Non-Patent Citations (5)
Title |
---|
何鹏: "基于MARS系统的X射线能谱CT研究", 《中国优秀博士论文全文数据库》 * |
何鹏等: "Color_CT技术研究简述", 《全国射线数字成像与CT新技术》 * |
陈平等: "一种基于能谱滤波分离的彩色CT成像方法", 《全国射线数字成像与CT新技术》 * |
陈平等: "变能量X射线多谱成像方法研究", 《光谱学与光谱分析》 * |
魏交统等: "基于主成分分析的递变能量X射线图像融合", 《中国体视学与图像分析》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106534866A (en) * | 2016-11-25 | 2017-03-22 | 中北大学 | Image compression method and apparatus |
CN109813259A (en) * | 2019-01-18 | 2019-05-28 | 中北大学 | Merge CT imaging method, storage medium and device |
CN109875593A (en) * | 2019-02-02 | 2019-06-14 | 中北大学 | CT imaging method, storage medium and device |
CN109875593B (en) * | 2019-02-02 | 2023-03-21 | 中北大学 | X-ray multi-spectrum separation imaging method, storage medium and device |
WO2021004157A1 (en) * | 2019-07-09 | 2021-01-14 | 上海联影医疗科技有限公司 | Medical scanning imaging method and apparatus, storage medium, and computer device |
CN111476856A (en) * | 2020-04-08 | 2020-07-31 | 中北大学 | Multispectral CT imaging method |
CN111476856B (en) * | 2020-04-08 | 2023-06-06 | 中北大学 | Multispectral CT imaging method |
CN112603345A (en) * | 2020-12-02 | 2021-04-06 | 赛诺威盛科技(北京)有限公司 | Model training method, multi-energy spectrum CT scanning method, device and electronic equipment |
Also Published As
Publication number | Publication date |
---|---|
CN103876772B (en) | 2015-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103876772A (en) | Multi-spectrum imaging method and device | |
Zhao et al. | An extended algebraic reconstruction technique (E-ART) for dual spectral CT | |
Taguchi et al. | An analytical model of the effects of pulse pileup on the energy spectrum recorded by energy resolved photon counting x‐ray detectors | |
CN102834736B (en) | Improving image quality in photon counting-mode detector systems | |
US9808216B2 (en) | Material decomposition of multi-spectral x-ray projections using neural networks | |
Karimi et al. | Segmentation of artifacts and anatomy in CT metal artifact reduction | |
DE102017207125A1 (en) | A method of performing a material decomposition using a dual-energy X-ray CT and a corresponding dual-energy X-ray CT apparatus | |
Ding et al. | Image‐domain multimaterial decomposition for dual‐energy CT based on prior information of material images | |
Lee et al. | Quantitative material decomposition using spectral computed tomography with an energy-resolved photon-counting detector | |
Tremblay et al. | A theoretical comparison of tissue parameter extraction methods for dual energy computed tomography | |
Latief et al. | The effect of X‐ray micro computed tomography image resolution on flow properties of porous rocks | |
Teles et al. | Rock porosity quantification by dual-energy X-ray computed microtomography | |
Lee et al. | A Monte Carlo simulation study of the effect of energy windows in computed tomography images based on an energy-resolved photon counting detector | |
Zhang et al. | Experimental implementation of a joint statistical image reconstruction method for proton stopping power mapping from dual‐energy CT data | |
Hsieh et al. | Improving pulse detection in multibin photon-counting detectors | |
Schmidt et al. | The effects of extending the spectral information acquired by a photon-counting detector for spectral CT | |
Hao et al. | A novel image optimization method for dual-energy computed tomography | |
US10448904B2 (en) | Decomposition method and apparatus based on basis material combination | |
Gondzio et al. | Material-separating regularizer for multi-energy x-ray tomography | |
Liao et al. | Noise Estimation for Single‐Slice Sinogram of Low‐Dose X‐Ray Computed Tomography Using Homogenous Patch | |
Ewert et al. | Cross-sectional imaging of building elements by new non-linear tomosynthesis techniques using imaging plates and 60Co radiation | |
Boisson et al. | NEMA NU 4-2008 validation and applications of the PET-SORTEO Monte Carlo simulations platform for the geometry of the Inveon PET preclinical scanner | |
Dyakowski et al. | A dual modality tomography system for imaging gas/solids flows | |
Cao et al. | A simulation-based study on the influence of the x-ray spectrum on the performance of multi-material beam hardening correction algorithms | |
Persson et al. | Bias–variance tradeoff in anticorrelated noise reduction for spectral CT |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230427 Address after: 030006 High tech Power Port 1110-1112, No. 226 Changzhi Road, Taiyuan Xuefu Park, Shanxi Comprehensive Reform Demonstration Zone, Taiyuan City, Shanxi Province Patentee after: SHANXI TEAMWORK PHOTOELECTRIC INDUSTRIES CO.,LTD. Address before: No. 3, Xueyuan Road, Beijing, Haidian District Patentee before: NORTH University OF CHINA |
|
TR01 | Transfer of patent right |