CN103822929A - Separated imaging method for superimposed targets based on multi-spectrum X rays - Google Patents
Separated imaging method for superimposed targets based on multi-spectrum X rays Download PDFInfo
- Publication number
- CN103822929A CN103822929A CN201410050519.0A CN201410050519A CN103822929A CN 103822929 A CN103822929 A CN 103822929A CN 201410050519 A CN201410050519 A CN 201410050519A CN 103822929 A CN103822929 A CN 103822929A
- Authority
- CN
- China
- Prior art keywords
- ray
- imaging
- decay
- voltage
- power spectrum
- 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
Images
Abstract
The invention provides a separated imaging method for superimposed targets based on multi-spectrum X rays. Imaging equipment sequentially comprises an X ray source, an imaging object, an X ray detector, a CCD (Charge Coupled Device), a data analyzing and processing device and terminal output equipment, which are put by an order. First, multiple X ray images are respectively shot for the imaging object in the same field of view through the change of the voltage of the X ray source. The imaging processes of the multiple X ray images are respectively described in combination of an X ray attenuation imaging formula and the gradation of the image so as to establish an X ray attenuation imaging equation set. The X ray attenuation imaging equation set is solved by an independent component analysis technology so as to obtain the attenuation thickness of the superimposed targets. The three-dimensional image of each superimposed target is reconstructed according to the thickness of the superimposed targets, that is, the superimposed targets in the original imaging object are separately imaged. Compared with a conventional X ray imaging method, the method has the advantages that the visual effect of the images is improved, and the medical diagnosis accuracy is improved.
Description
Technical field
The x-ray imaging that utilizes that the present invention relates to Medical Imaging and field of non destructive testing carries out diagnosis and detection, is specifically related to a kind of X ray three-dimensional imaging new technology.
Background technology
X ray is widely used in medical diagnosis, industrial flaw detection and Non-Destructive Testing.Typical x-ray imaging system is mainly made up of x-ray source, X-ray detector, data processing unit and image output device, penetrates the X ray intensity being attenuated by object form radioscopic image by record.Material attenuation coefficient changes with ray energy, and is determined by the chemical composition of material.X ray is obviously greater than the decay to low atomic number material to the decay of high atomic number material, causes the radioscopic image contrast of high atomic number material high, and the picture contrast of low atomic number object is lower.Therefore, in radioscopic image, low atomic number material is easy to superimposed material in the above and blocks.
In this invention, adopt X ray multi-power spectrum technology, object to be detected is taken to several not images of transmissive voltage, otherness according to different materials attenuation coefficient with change in voltage, set up x-ray imaging system of equations, utilize blind source analytic approach independent component analysis to solve system of equations, can try to achieve the transmission thickness of different materials, thereby different materials is separated from single X ray two dimensional image, reached three-dimensional imaging object.This invention is different from image-forming principle and the method for the classical technology CT of current X ray three-dimensional imaging, imaging is without multi-angled shooting and machinery rotation, in imaging process, only energy requirement switches, and therefore its unique distinction and technological difficulties are to apply multi-power spectrum imaging and independent component analysis separation of different materials.
Summary of the invention
Technical matters: the object of this invention is to provide a kind of stack target method for separate imaging based on multi-power spectrum X ray, adopt X ray multi-power spectrum technology, utilize different material to X ray decay different and same substance the X ray different character that decays is set up to multi-power spectrum X ray attenuation equation group, solve multidimensional nonlinear system of equations in conjunction with blind source separating treatment technology, reach overlapping target separate imaging object in object, strengthen image observation effect, improve image authentication rate.
Technical scheme: a kind of its imaging device of stack target method for separate imaging based on multi-power spectrum X ray of the present invention is followed successively by x-ray source, imaging object, X-ray detector, CCD, data analysis and processor and terminal output device according to placement order; First by changing x-ray source voltage, imaging object is taken respectively to several radioscopic images under same field; Respectively several radioscopic image imaging processes are described to set up X ray decay imaging equation group in conjunction with X ray decay imaging formula and gradation of image; Adopt Independent Component Analysis Technology to solve the X ray decay imaging equation group target decay thickness that obtains superposeing; Rebuilding each stack target 3-D view according to stack target thickness is by the target separate imaging that superposes in former imaging object.
Described X ray decay imaging equation is I=I
0× exp (μ T), wherein: I
0for initial ray intensity; I is the transmitted intensity after decay; μ is the attenuation coefficient that penetrates object; T is the thickness that penetrates object; Each radiographic source voltage U
icorresponding attenuation coefficient mu
i, therefore multiple radiographic source voltage U
icorresponding multiple X ray decay imaging equations, can form X ray decay imaging equation group jointly.
Described X ray decay imaging equation group is expressed as attenuation coefficient Matrix Multiplication and equals pixel grey scale matrix with decay thickness matrix, the decay thickness that can try to achieve each stack target by matrixing solving equation can obtain the superposeing high three-dimensional information of length and width of target.
Described matrixing solving equation group, adopt Independent Component Analysis Technology in the analytic approach of blind source to solve X ray decay imaging equation group, method is: according to the overlapping target of the interior of articles composition material separately of grasping, obtain corresponding attenuation coefficient μ under the different voltage of different materials
ij, subscript i distinguishes different voltage; Subscript j distinguishes different materials, brings X ray decay imaging equation group into and carries out blind separation, the thickness Matrix Calculating of stack target can be obtained.
Change x-ray source voltage and determine according to the attribute of studied object and thick and heavy degree, to penetrate object and imaging gray-level clearly as standard; X-ray source tube voltage selects to follow the application standard of conventional X ray aspect medical diagnosis, industrial flaw detection and Non-Destructive Testing; The family curve changing with x-ray source voltage U according to the X ray attenuation coefficient mu of material for voltage value interval is determined.
Several described radioscopic images adopt voltage automatic switchover technology or photon counter technology to obtain; Wherein, voltage automatic switchover technology is x-ray source voltage to be automatically adjusted according to selected voltage range and voltage spaces by voltage controller, takes several radioscopic images; Photon counting technique is to use photon counter to obtain multiband narrow band x-ray image simultaneously.
While using voltage automatic switchover technology, need to before imaging detector, lay and consider wiregrating to guarantee ray energy spectrum unicity.
Described X ray decay imaging equation group, setting up before system of equations, need to be to multiple image denoising, denoising method adopts Independent Component Analysis Technology, and noise and picture signal are considered as to two kinds of components and then separate, and retains picture signal.
Beneficial effect: adopt X ray multi-power spectrum technology, utilize different material to X ray decay different and same substance the X ray different character that decays is set up to multi-power spectrum X ray attenuation equation group, solve multidimensional nonlinear system of equations in conjunction with blind source separating treatment technology, reach overlapping target separate imaging object in object, can independent three-dimensional imaging after overlapping target is separated, overcome doubling of the image shortcoming, improve radioscopic image observation effect, conveniently according to image, imaging object details is judged.
Accompanying drawing explanation
Fig. 1 is the method that conventional X ray gathers imaging;
Fig. 2 is image acquisition, analysis and the disposal route that the present invention proposes.
Embodiment
A kind of its imaging device of stack target method for separate imaging based on multi-power spectrum X ray of the present invention is followed successively by x-ray source, imaging object, X-ray detector, CCD, data analysis and processor and terminal output device according to placement order; First by changing x-ray source voltage, imaging object is taken respectively to several radioscopic images under same field; Respectively several radioscopic image imaging processes are described to set up X ray decay imaging equation group in conjunction with X ray decay imaging formula and gradation of image; Adopt Independent Component Analysis Technology to solve the X ray decay imaging equation group target decay thickness that obtains superposeing; Rebuilding each stack target 3-D view according to stack target thickness is by the target separate imaging that superposes in former imaging object.
The concrete steps of embodiment of the present invention are:
Step 1: choose x-ray source according to object to be imaged attribute and thick and heavy degree, to penetrate object and imaging gray-level clearly as standard.
Step 2: build as shown in Figure 2 multi-power spectrum x-ray imaging system, system composition device is followed successively by x-ray source, imaging object, X-ray detector, CCD, data analysis and processor and terminal output device according to placement order.
Step 3: utilize the imaging respectively of different x-ray source voltage penetration imaging object, obtain several radioscopic images under same field.If use voltage automatic switchover technology, need to before imaging detector, lay and consider wiregrating to guarantee ray energy spectrum unicity.The family curve that voltage value interval changes with x-ray source voltage U according to the X ray attenuation coefficient mu of material is determined.
Step 4: several radioscopic images of taking are carried out to image denoising and process and image normalization processing.Denoising method adopts Independent Component Analysis Technology, and noise and picture signal are considered as to two kinds of components and then separate, and retains picture signal.
Step 5: respectively several radioscopic image imaging processes are described to set up X ray decay imaging equation group in conjunction with X ray decay imaging formula and radioscopic image gray scale.
Step 6: adopt Independent Component Analysis Technology (ICA) to carry out the separation of blind source, solve the X ray decay imaging equation group target decay thickness that obtains superposeing.
Step 7: rebuilding each stack target 3-D view according to stack target thickness is by the target separate imaging that superposes in former imaging object.
Step 8: the separate picture that step 7 is obtained carries out pseudo-colours processing, can obtain separate targets coloured image; Recycle image fusion technology, can rebuild the colored X ray stereo-picture of former imaging object.
Claims (8)
1. the stack target method for separate imaging based on multi-power spectrum X ray, is characterized in that its imaging device is followed successively by x-ray source, imaging object, X-ray detector, CCD, data analysis and processor and terminal output device according to placement order; First utilize different x-ray source voltage to take respectively several radioscopic images under same field to imaging object; Respectively several radioscopic image imaging processes are described to set up X ray decay imaging equation group in conjunction with X ray decay imaging formula and gradation of image; Adopt Independent Component Analysis Technology to solve the X ray decay imaging equation group target decay thickness that obtains superposeing; Rebuilding each stack target 3-D view according to stack target thickness is by the target separate imaging that superposes in former imaging object.
2. the stack target method for separate imaging based on multi-power spectrum X ray according to claim 1, is characterized in that described X ray decay imaging equation is I=I
0× exp (μ T), wherein: I
0for initial ray intensity; I is the transmitted intensity after decay; μ is the X ray attenuation coefficient that penetrates object; T is the thickness that penetrates object; Each radiographic source voltage U
icorresponding attenuation coefficient mu
i, therefore multiple radiographic source voltage U
icorresponding multiple X ray decay imaging equations, can form X ray decay imaging equation group jointly.
3. the stack target method for separate imaging based on multi-power spectrum X ray according to claim 1, it is characterized in that described X ray decay imaging equation group is expressed as attenuation coefficient Matrix Multiplication and equals pixel grey scale matrix with decay thickness matrix, the decay thickness that can try to achieve each stack target by matrixing solving equation can obtain the superposeing high three-dimensional information of length and width of target.
4. the stack target method for separate imaging based on multi-power spectrum X ray according to claim 3, it is characterized in that described matrixing solving equation group, adopt Independent Component Analysis Technology in the analytic approach of blind source to solve X ray decay imaging equation group, method is: according to the overlapping target of the interior of articles composition material separately of grasping, obtain corresponding attenuation coefficient μ under the different voltage of different materials
ij, subscript i distinguishes different voltage; Subscript j distinguishes different materials, brings X ray decay imaging equation group into and carries out blind separation, the thickness Matrix Calculating of stack target can be obtained.
5. the stack target method for separate imaging based on multi-power spectrum X ray according to claim 1, is characterized in that changing x-ray source voltage and determines according to the attribute of studied object and thick and heavy degree, to penetrate object and imaging gray-level clearly as standard; X-ray source tube voltage selects to follow the application standard of conventional X ray aspect medical diagnosis, industrial flaw detection and Non-Destructive Testing; The family curve changing with x-ray source voltage U according to the X ray attenuation coefficient mu of material for voltage value interval is determined.
6. the stack target method for separate imaging based on multi-power spectrum X ray according to claim 1, is characterized in that several described radioscopic images adopt voltage automatic switchover technology or photon counter technology to obtain; Wherein, voltage automatic switchover technology is x-ray source voltage to be automatically adjusted according to selected voltage range and voltage spaces by voltage controller, takes several radioscopic images; Photon counting technique is to use photon counter to obtain multiband narrow band x-ray image simultaneously.
7. the stack target method for separate imaging based on multi-power spectrum X ray according to claim 6 while it is characterized in that using voltage automatic switchover technology, need to be laid and consider wiregrating to guarantee ray energy spectrum unicity before imaging detector.
8. the stack target method for separate imaging based on multi-power spectrum X ray according to claim 1, it is characterized in that described X ray decay imaging equation group, setting up before system of equations, need to be to multiple image denoising, denoising method adopts Independent Component Analysis Technology, noise and picture signal are considered as to two kinds of components and then separate, retain picture signal.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410050519.0A CN103822929A (en) | 2014-02-13 | 2014-02-13 | Separated imaging method for superimposed targets based on multi-spectrum X rays |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410050519.0A CN103822929A (en) | 2014-02-13 | 2014-02-13 | Separated imaging method for superimposed targets based on multi-spectrum X rays |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN103822929A true CN103822929A (en) | 2014-05-28 |
Family
ID=50758096
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410050519.0A Pending CN103822929A (en) | 2014-02-13 | 2014-02-13 | Separated imaging method for superimposed targets based on multi-spectrum X rays |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103822929A (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104091356A (en) * | 2014-07-04 | 2014-10-08 | 南京邮电大学 | X-ray medical image objective reconstruction based on independent component analysis |
| CN104545960A (en) * | 2014-12-18 | 2015-04-29 | 沈阳东软医疗***有限公司 | Dual-energy CT (Computed Tomography) scanning voltage selection method and equipment |
| CN105654528A (en) * | 2016-01-04 | 2016-06-08 | 南京邮电大学 | Compressed sensing-based pleuripotent X-ray separate imaging method |
| CN108181326A (en) * | 2016-12-07 | 2018-06-19 | 同方威视技术股份有限公司 | Multi-power spectrum x-ray imaging system and the method for carrying out Object Classification to article to be measured using multi-power spectrum x-ray imaging system |
| CN109429418A (en) * | 2017-08-21 | 2019-03-05 | 上海影信息科技有限公司 | X-ray generator and the X-ray multipotency imaging system for using it |
| CN109997029A (en) * | 2016-11-29 | 2019-07-09 | 莱特拉姆有限责任公司 | The multi-energy X-ray detected to the foreign matter on conveyer absorbs imaging |
| CN112748132A (en) * | 2020-12-29 | 2021-05-04 | 四川赛康智能科技股份有限公司 | Interference removing method for aluminum sheath in X-ray detection cable imaging |
| CN113926737A (en) * | 2021-09-26 | 2022-01-14 | 河南中平自动化股份有限公司 | Coal gangue device utilizing dual-energy dual-machine X-ray and identification method thereof |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004008460A (en) * | 2002-06-06 | 2004-01-15 | Kawasaki Heavy Ind Ltd | X-ray energy analysis imaging apparatus |
| US20070053489A1 (en) * | 2005-04-25 | 2007-03-08 | The University Of North Carolina At Chapel Hill | X-ray imaging systems and methods using temporal digital signal processing for reducing noise and for obtaining multiple images simultaneously |
| CN100391407C (en) * | 2002-10-07 | 2008-06-04 | 通用电气公司 | Continuous scanning RAD chromatorgraphy X-ray lighting combination system and method |
-
2014
- 2014-02-13 CN CN201410050519.0A patent/CN103822929A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004008460A (en) * | 2002-06-06 | 2004-01-15 | Kawasaki Heavy Ind Ltd | X-ray energy analysis imaging apparatus |
| CN100391407C (en) * | 2002-10-07 | 2008-06-04 | 通用电气公司 | Continuous scanning RAD chromatorgraphy X-ray lighting combination system and method |
| US20070053489A1 (en) * | 2005-04-25 | 2007-03-08 | The University Of North Carolina At Chapel Hill | X-ray imaging systems and methods using temporal digital signal processing for reducing noise and for obtaining multiple images simultaneously |
Non-Patent Citations (3)
| Title |
|---|
| DONG-GOO KANG等: "Multiple object decomposition based on independent component analysis of multi-energy X-Ray projections", 《IMAGE PROCESSING（ICIP），2009 16TH IEEE INTERNATIONAL CONFERENCE ON》 * |
| SONGYUAN TANG等: "Application of ICA to X –ray coronary digital subtraction angiography", 《NEUROCOMPUTING》 * |
| 隋海洋等: "多能谱射线检测技术研究", 《2005年全国射线检测技术及加速器检测设备和应用技术交流会论文集》 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104091356A (en) * | 2014-07-04 | 2014-10-08 | 南京邮电大学 | X-ray medical image objective reconstruction based on independent component analysis |
| CN104545960A (en) * | 2014-12-18 | 2015-04-29 | 沈阳东软医疗***有限公司 | Dual-energy CT (Computed Tomography) scanning voltage selection method and equipment |
| US10105110B2 (en) | 2014-12-18 | 2018-10-23 | Shenyang Neusoft Medical Systems Co., Ltd. | Selecting scanning voltages for dual energy CT scanning |
| CN105654528A (en) * | 2016-01-04 | 2016-06-08 | 南京邮电大学 | Compressed sensing-based pleuripotent X-ray separate imaging method |
| CN105654528B (en) * | 2016-01-04 | 2018-11-27 | 南京邮电大学 | Compressed sensing based multipotency X-ray method for separate imaging |
| CN109997029A (en) * | 2016-11-29 | 2019-07-09 | 莱特拉姆有限责任公司 | The multi-energy X-ray detected to the foreign matter on conveyer absorbs imaging |
| CN108181326A (en) * | 2016-12-07 | 2018-06-19 | 同方威视技术股份有限公司 | Multi-power spectrum x-ray imaging system and the method for carrying out Object Classification to article to be measured using multi-power spectrum x-ray imaging system |
| CN108181327A (en) * | 2016-12-07 | 2018-06-19 | 同方威视技术股份有限公司 | Multi-power spectrum x-ray imaging system and the method for carrying out Object Classification to article to be measured using multi-power spectrum x-ray imaging system |
| CN108181327B (en) * | 2016-12-07 | 2021-02-05 | 同方威视技术股份有限公司 | Multi-energy spectrum X-ray imaging system and method for identifying substance of object to be detected by using multi-energy spectrum X-ray imaging system |
| CN109429418A (en) * | 2017-08-21 | 2019-03-05 | 上海影信息科技有限公司 | X-ray generator and the X-ray multipotency imaging system for using it |
| CN112748132A (en) * | 2020-12-29 | 2021-05-04 | 四川赛康智能科技股份有限公司 | Interference removing method for aluminum sheath in X-ray detection cable imaging |
| CN113926737A (en) * | 2021-09-26 | 2022-01-14 | 河南中平自动化股份有限公司 | Coal gangue device utilizing dual-energy dual-machine X-ray and identification method thereof |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103822929A (en) | Separated imaging method for superimposed targets based on multi-spectrum X rays | |
| US9625310B2 (en) | Systems and methods for sorting and summing signals from an imaging detector | |
| DE60215932T2 (en) | X INVESTIGATION SYSTEM | |
| US8306180B2 (en) | Image reconstruction method for high-energy, dual-energy CT system | |
| 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 | |
| EP3242126A1 (en) | Dual-energy ray imaging method and system | |
| ATE476142T1 (en) | METHOD AND SYSTEM FOR BINOCULAR STEREOSCOPIC IMAGING BY SCANNING RADIOGRAPHY | |
| US8254656B2 (en) | Methods and system for selective resolution improvement in computed tomography | |
| WO2010135901A1 (en) | Pseudo dual-energy under-sampling material discriminating system and method | |
| US8989346B2 (en) | Bone mineral density analysis method, bone mineral density analysis apparatus, and recording medium | |
| CN103808739B (en) | A kind of safety inspection device of transmission imaging and back scattering imaging integration | |
| WO2010135900A1 (en) | Dual-energy under-sampling material discriminating method and system | |
| CN105264361A (en) | High-resolution computed tomography | |
| JP2008268076A (en) | Non-destructive discrimination method, and device | |
| Paziresh et al. | Tomography of atomic number and density of materials using dual-energy imaging and the Alvarez and Macovski attenuation model | |
| CN103202707A (en) | X-ray Ct (computed Tomography) Device | |
| CN104274201A (en) | Method, system and equipment for tomography of mammary gland and image acquisition and processing method | |
| US20200187893A1 (en) | Beam hardening correction in x-ray dark-field imaging | |
| EP3388868A1 (en) | Inspection devices and methods for inspecting a container | |
| CN105784741A (en) | Multi-visual-angle X ray safety check device | |
| CN204882414U (en) | Various visual angles X ray safety check device | |
| Firsching et al. | Detection of enclosed diamonds using dual energy X-ray imaging | |
| Balaskó et al. | Neutron-, gamma-and X-ray three-dimensional computed tomography at the Budapest research reactor site | |
| CN102670231B (en) | Determine method and the computer tomograph of the distance of check object and focus | |
| CN105222730B (en) | A kind of industry CT physical dimension measuring method based on image restoration |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140528 |