CN110097597A - A kind of coordinate corresponding method of the series of X light image of object - Google Patents
A kind of coordinate corresponding method of the series of X light image of object Download PDFInfo
- Publication number
- CN110097597A CN110097597A CN201910376664.0A CN201910376664A CN110097597A CN 110097597 A CN110097597 A CN 110097597A CN 201910376664 A CN201910376664 A CN 201910376664A CN 110097597 A CN110097597 A CN 110097597A
- Authority
- CN
- China
- Prior art keywords
- actinogram
- width
- series
- light image
- coordinate
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/66—Analysis of geometric attributes of image moments or centre of gravity
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/68—Analysis of geometric attributes of symmetry
-
- 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/10116—X-ray image
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Multimedia (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Image Processing (AREA)
Abstract
The invention discloses a kind of coordinate corresponding methods of the series of X light image of object.This method includes the following steps: the series of X light image of a. record object, and the X-ray intensity of each recording unit records in every width actinogram is summed, every width X-ray overall strength is obtained, obtains A；B. the X-ray intensity of each recording unit in every width actinogram is multiplied with corresponding two-dimensional coordinate, and summed, obtain B；C.C=B/A obtains the light intensity barycentric coodinates of every width actinogram；D. using the corresponding attribute of light intensity barycentric coodinates of every width actinogram in series of X light image, the two-dimensional coordinate corresponding relationship of every width actinogram in series of X light image is provided.Each width actinogram coordinate that this method can be realized different observed bearings accurately corresponds to, quantitative analysis each width actinogram coordinate time in different time periods develops, coordinate pair answers precision to reach imaging system spatial discrimination level, has wide and important application prospect in inertial confinement fusion, Laboratory astrophysics or high-energy density physics.
Description
Technical field
The invention belongs to actinogram data processing fields, and in particular to a kind of coordinate counterparty of the series of X light image of object
Method.
Background technique
X-ray imaging is most widely used to inertial confinement fusion, Laboratory astrophysics and high-energy density physics
In scientific research, the two dimension for diagnosing substance etc. under objective object, such as target capsule of fusion, radiation source and extreme condition is empty
Between, time resolution information.Usually in the diagnosis that two-dimensional space is differentiated, people can obtain several actinograms of different observed bearings,
Above-mentioned scientific research needs by coordinate accurately to be corresponded these images, because can be provided by correspondence such as electric
The important informations such as sub- temperature, electron density, symmetrical analysis.In time-resolved diagnosis, people can be become at any time again
Several actinograms changed, related science research need the coordinate time of these images of quantitative analysis to develop to study symmetry change
The problems such as change.
And existing actinogram data processing has the disadvantage that 1, the imaging of different observed bearing spatial discrimination X-rays is examined
Disconnected, luminous according to projection theory and object inhomogeneities, each width actinogram of gained difference observed bearing is significantly different
, this includes the difference of space strength distribution and shape geometry, therefore these actinograms are to be difficult to realize coordinate accurately to correspond to
's；2, for different time sections time resolution X-ray imaging diagnosis, since object is with temporal evolution, X-ray luminance is
Variation, gained each width actinogram in different time periods be also it is visibly different, this also includes that space strength distribution and shape are several
What difference, but how the coordinate time evolution of quantitative analysis these images still lacks corresponding coordinate corresponding method.
Summary of the invention
Technical problem to be solved by the invention is to provide a kind of coordinate corresponding methods of the series of X light image of object.
The coordinate corresponding method of the series of X light image of object of the invention the following steps are included:
A. the series of X light image for recording object, the X-ray intensity of each recording unit records in every width actinogram is summed, and is obtained every
Width X-ray overall strength, obtains A；
B. the X-ray intensity of each recording unit in every width actinogram is multiplied with corresponding two-dimensional coordinate, and summed, obtain B；
C.C=B/A obtains the light intensity barycentric coodinates of every width actinogram；
D. it using the corresponding attribute of light intensity barycentric coodinates of every width actinogram in series of X light image, provides in series of X light image
The two-dimensional coordinate corresponding relationship of every width actinogram.
The actinogram is the X-ray as formed by one of pin hole, Kirkpatrick-Baez mirror or bent crystal image device
Picture.
The series of X light image includes the actinogram recorded for the same object in different direction, or when different
Between segment record actinogram.
The coordinate corresponding method of the series of X light image of object of the invention is applied to inertial confinement fusion, laboratory celestial body
Object actinogram data processing in physics or high-energy density physics.
In every width X-ray overall strength A, every width actinogram in the coordinate corresponding method of the series of X light image of object of the invention
X-ray intensity and the summation B of corresponding two-dimensional coordinate product of each recording unit, the light intensity barycentric coodinates C of every width actinogram and every
The two-dimensional coordinate corresponding relationship of width actinogram can be obtained by hand-kept, software or programming mode.
The coordinate corresponding method of the series of X light image of object of the invention has the advantage that
1. each width actinogram coordinate that can be realized different observed bearings accurately corresponds to.
2. can be realized quantitative analysis each width actinogram coordinate time in different time periods to develop.
3. it is horizontal that coordinate pair answers precision to reach imaging system spatial discrimination.
The coordinate that the coordinate corresponding method of the series of X light image of object of the invention can be realized series of X light image is accurately right
It answers and quantitative analysis, has in inertial confinement fusion, Laboratory astrophysics or high-energy density physics wide and again
Want application prospect.
Detailed description of the invention
Fig. 1 is projection coordinate's schematic diagram of the coordinate corresponding method of the series of X light image of object of the invention.
Specific embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.
The coordinate corresponding method of the series of X light image of object of the invention the following steps are included:
A. the series of X light image for recording object, the X-ray intensity of each recording unit records in every width actinogram is summed, and is obtained every
Width X-ray overall strength, obtains A；
B. the X-ray intensity of each recording unit in every width actinogram is multiplied with corresponding two-dimensional coordinate, and summed, obtain B；
C.C=B/A obtains the light intensity barycentric coodinates of every width actinogram；
D. it using the corresponding attribute of light intensity barycentric coodinates of every width actinogram in series of X light image, provides in series of X light image
The two-dimensional coordinate corresponding relationship of every width actinogram.
The actinogram is the X-ray as formed by one of pin hole, Kirkpatrick-Baez mirror or bent crystal image device
Picture.
The series of X light image includes the actinogram recorded for the same object in different direction, or when different
Between segment record actinogram.
The coordinate corresponding method of the series of X light image of object of the invention is applied to inertial confinement fusion, laboratory celestial body
Object actinogram data processing in physics or high-energy density physics.
Embodiment 1
The actinogram is the actinogram as formed by pin hole.The series of X light image is for the same object in not Tongfang
The actinogram of position record.The coordinate corresponding method of the series of X light image of object of the invention is applied to the mesh in inertial confinement fusion
Mark the data processing of object actinogram.
According to basic physical knowledge, there is the object of quality to carry out quality integral two-dimensional projection to different direction, center of gravity exists
It is fixed on each perspective plane, and is one-to-one.The center of gravity of namely three-dimension object centainly passes through along projecting direction
The center of gravity of two-dimensional projection.If changing mass of object into X source light intensity, conclusion is consistent, i.e. the light of object X source
Strong center of gravity centainly passes through the light intensity center of gravity of two-dimensional projection along projecting direction (observed bearing).
The light distribution for enabling object X source is, see Fig. 1, then along observed bearingBe projected as
(1)
(2)
Wherein,It is?The projection of plane withThe angle of axis,It isWithThe angle of plane,It is three-dimensional
The spin matrix of coordinate
(3)
The object X source light intensity barycentric coodinates are enabled to be, the light intensity barycentric coodinates of every width actinogram of projection are C, the definition of foundation center of gravity
(4)
T is object X source intensityVolume point i.e. objective object X source overall strength.Then project (observation)
Every width actinogram light intensity barycentric coodinates C
(5)
Wherein,For X-ray intensity and the corresponding two dimension of each recording unit in every width actinogram
The summation B of coordinate product,For projection (observation) actinogram light intensityLine Integral be project (observation) every width
X-ray overall strength A.Relationship between T is
(6)
WhereinIt is projection (observation) actinogram to solid angle shared by object X source, R is constant relevant to r.
It isAs obtained by rotating twice, relationship between the two is still determined by formula (2),
I.e.
(7)
Wherein S is the first row of matrix K and the matrix that the second row is constituted, i.e.,
(8)
Therefore, by extracting every width X-ray overall strength A from each width actinogram i.e., then extract each note in every width actinogram
X-ray intensity and the summation B of corresponding two-dimensional coordinate product for recording unit, can obtain the light intensity of every width actinogram by formula (5)
Barycentric coodinates C.The corresponding attribute of light intensity barycentric coodinates of every width actinogram in series of X light image is utilized according to formula (7),
Provide the two-dimensional coordinate corresponding relationship of every width actinogram in series of X light image.The coordinate pair obtained in this way answers precision by pin-hole imaging
The spatial discrimination of system determines, is 10 μm of magnitudes.
In addition, carry out the Quantitative study for each width actinogram that there is power spectrum to differentiate after coordinate correspondence relationship determines, thus
According to X-ray radiative property, such as common bremsstrahlung mechanism, electron temperature, electron density can extract.According to all coordinate pairs
Answer the three dimensional symmetry of situation assessment object X-ray transmitting.
Embodiment 2
The present embodiment is identical as the embodiment of embodiment 1, the difference is that, the actinogram is by Kirkpatrick-
Actinogram formed by Baez mirror.The series of X light image is the actinogram for the same object in different time segment record.This
The coordinate corresponding method of the series of X light image of the object of invention is applied to the object actinogram data in high-energy density physics
Processing.After the light intensity barycentric coodinates C for obtaining every width actinogram by formula (5), by compare each width actinogram light intensity center of gravity with
The variation of the distance between geometric center, can quantitative analysis each width actinogram coordinate time in different time periods develop, to grind
Study carefully object symmetry variation tendency.The coordinate pair obtained in this way answers precision by the sky of Kirkpatrick-Baez mirror imaging system
Between differentiate determine, be 2.5 μm of magnitudes.
The present embodiment applies also for the object actinogram data processing in Laboratory astrophysics.
Those of ordinary skill in the art will understand that the embodiments described herein, which is to help reader, understands this hair
Bright principle, it should be understood that protection scope of the present invention is not limited to such specific embodiments and embodiments.This field
Those of ordinary skill disclosed the technical disclosures can make according to the present invention and various not depart from the other each of essence of the invention
The specific variations and combinations of kind, these variations and combinations are still within the scope of the present invention.
Claims (4)
1. a kind of coordinate corresponding method of the series of X light image of object, which is characterized in that the method the following steps are included:
A. the series of X light image for recording object, the X-ray intensity of each recording unit records in every width actinogram is summed, and is obtained every
Width X-ray overall strength, obtains A；
B. the X-ray intensity of each recording unit in every width actinogram is multiplied with corresponding two-dimensional coordinate, and summed, obtain B；
C.C=B/A obtains the light intensity barycentric coodinates of every width actinogram；
D. it using the corresponding attribute of light intensity barycentric coodinates of every width actinogram in series of X light image, provides in series of X light image
The two-dimensional coordinate corresponding relationship of every width actinogram.
2. the coordinate corresponding method of the series of X light image of object according to claim 1, which is characterized in that the X-ray
As being the actinogram as formed by one of pin hole, Kirkpatrick-Baez mirror or bent crystal image device.
3. the coordinate corresponding method of the series of X light image of object according to claim 1, which is characterized in that described is
Column actinogram includes the actinogram recorded for the same object in different direction, or the actinogram in different time segment record.
4. the coordinate corresponding method of the series of X light image of any object according to claim 1 ~ 3, which is characterized in that institute
The coordinate corresponding method stated is applied to the object X in inertial confinement fusion, Laboratory astrophysics or high-energy density physics
Light image data processing.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910376664.0A CN110097597B (en) | 2019-05-05 | 2019-05-05 | Coordinate corresponding method for series X-ray images of target object |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910376664.0A CN110097597B (en) | 2019-05-05 | 2019-05-05 | Coordinate corresponding method for series X-ray images of target object |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110097597A true CN110097597A (en) | 2019-08-06 |
| CN110097597B CN110097597B (en) | 2022-02-11 |
Family
ID=67447192
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910376664.0A Active CN110097597B (en) | 2019-05-05 | 2019-05-05 | Coordinate corresponding method for series X-ray images of target object |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110097597B (en) |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH01262453A (en) * | 1988-04-13 | 1989-10-19 | Rigaku Corp | Method of reading diffracted x-ray intensity of x-ray diffraction device using accumulation type fluorescent plate |
| US5333065A (en) * | 1991-12-09 | 1994-07-26 | Agfa-Gevaert | Signal to density mapping with controlled contrast |
| JP2009249445A (en) * | 2008-04-03 | 2009-10-29 | Mitsubishi Chemicals Corp | Fluorescent substance and method for producing the same, fluorescent substance-containing composition, light-emitting device, lighting device and image display device |
| EP1834584A4 (en) * | 2004-12-21 | 2010-11-17 | Univ Gunma Nat Univ Corp | Region-in-object measuring system, computing device for measuring region in object, program for measuring region in object, computer-readable recording medium where the program is recorded |
| US20110194787A1 (en) * | 2010-02-08 | 2011-08-11 | James Jiwen Chun | Constructing Three Dimensional Images Using Panoramic Images |
| CN102542600A (en) * | 2011-12-14 | 2012-07-04 | 北京工业大学 | Simulated projection DRR( digitally reconstructed radiograph) generating method based on CUDA (compute unified device architecture) technology |
| CN102530850A (en) * | 2012-03-14 | 2012-07-04 | 哈尔滨工业大学 | Method for millimeter-sized micro nanostructure nano carving and processing through adopting antifrictional metal (AFM) needle |
| CN103018896A (en) * | 2012-12-19 | 2013-04-03 | 哈尔滨工业大学 | Three-point high-precision large-aperture electric reflector frame |
| WO2013164368A1 (en) * | 2012-05-01 | 2013-11-07 | Universität Bern | Image distortion correction and robust phantom detection |
| CN103512911A (en) * | 2012-06-18 | 2014-01-15 | 上海梅山钢铁股份有限公司 | Metallurgy miscellaneous material fast spectral analysis method |
| CN103815923A (en) * | 2012-11-16 | 2014-05-28 | 索尼公司 | Image processing apparatus, image processing method, and program |
| CN104807843A (en) * | 2015-04-13 | 2015-07-29 | 江阴市产品质量监督检验所 | Method for measuring sulfur and phosphorus in soldering flux with X-ray fluorescence spectrometry |
| CN107635468A (en) * | 2015-03-18 | 2018-01-26 | 韩国威泰有限公司 | For rebuilding the apparatus and method of medical image |
| CN107941830A (en) * | 2017-12-27 | 2018-04-20 | 钢研纳克检测技术股份有限公司 | The distributional analysis Image Acquisition and data handling system of Xray fluorescence spectrometer |
| CN109147049A (en) * | 2018-07-28 | 2019-01-04 | 天津大学 | Image reconstruction method for X-ray optical dynamic therapy |
| CN109637691A (en) * | 2018-12-13 | 2019-04-16 | 中国工程物理研究院激光聚变研究中心 | A kind of choosing of binaryzation X-ray can device and preparation method thereof |
-
2019
- 2019-05-05 CN CN201910376664.0A patent/CN110097597B/en active Active
Patent Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH01262453A (en) * | 1988-04-13 | 1989-10-19 | Rigaku Corp | Method of reading diffracted x-ray intensity of x-ray diffraction device using accumulation type fluorescent plate |
| US5333065A (en) * | 1991-12-09 | 1994-07-26 | Agfa-Gevaert | Signal to density mapping with controlled contrast |
| EP1834584A4 (en) * | 2004-12-21 | 2010-11-17 | Univ Gunma Nat Univ Corp | Region-in-object measuring system, computing device for measuring region in object, program for measuring region in object, computer-readable recording medium where the program is recorded |
| JP2009249445A (en) * | 2008-04-03 | 2009-10-29 | Mitsubishi Chemicals Corp | Fluorescent substance and method for producing the same, fluorescent substance-containing composition, light-emitting device, lighting device and image display device |
| US20110194787A1 (en) * | 2010-02-08 | 2011-08-11 | James Jiwen Chun | Constructing Three Dimensional Images Using Panoramic Images |
| CN102542600A (en) * | 2011-12-14 | 2012-07-04 | 北京工业大学 | Simulated projection DRR( digitally reconstructed radiograph) generating method based on CUDA (compute unified device architecture) technology |
| CN102530850A (en) * | 2012-03-14 | 2012-07-04 | 哈尔滨工业大学 | Method for millimeter-sized micro nanostructure nano carving and processing through adopting antifrictional metal (AFM) needle |
| WO2013164368A1 (en) * | 2012-05-01 | 2013-11-07 | Universität Bern | Image distortion correction and robust phantom detection |
| CN103512911A (en) * | 2012-06-18 | 2014-01-15 | 上海梅山钢铁股份有限公司 | Metallurgy miscellaneous material fast spectral analysis method |
| CN103815923A (en) * | 2012-11-16 | 2014-05-28 | 索尼公司 | Image processing apparatus, image processing method, and program |
| CN103018896A (en) * | 2012-12-19 | 2013-04-03 | 哈尔滨工业大学 | Three-point high-precision large-aperture electric reflector frame |
| CN107635468A (en) * | 2015-03-18 | 2018-01-26 | 韩国威泰有限公司 | For rebuilding the apparatus and method of medical image |
| CN104807843A (en) * | 2015-04-13 | 2015-07-29 | 江阴市产品质量监督检验所 | Method for measuring sulfur and phosphorus in soldering flux with X-ray fluorescence spectrometry |
| CN107941830A (en) * | 2017-12-27 | 2018-04-20 | 钢研纳克检测技术股份有限公司 | The distributional analysis Image Acquisition and data handling system of Xray fluorescence spectrometer |
| CN109147049A (en) * | 2018-07-28 | 2019-01-04 | 天津大学 | Image reconstruction method for X-ray optical dynamic therapy |
| CN109637691A (en) * | 2018-12-13 | 2019-04-16 | 中国工程物理研究院激光聚变研究中心 | A kind of choosing of binaryzation X-ray can device and preparation method thereof |
Non-Patent Citations (6)
| Title |
|---|
| D. R. WILKINS等: ""Understanding X-ray reflection emissivity profiles in AGN: locating the X-ray source"", 《X-RAY REFLECTION EMISSIVITY PROFILES IN AGN》 * |
| H TODA1等: ""Three-dimensional strain mapping using in situ X-ray synchrotron microtomography"", 《SPECIAL ISSUE PAPER》 * |
| KEITH C. GENDREAU等: ""The Neutron star Interior Composition ExploreR (NICER): an Explorer mission of opportunity for soft x-ray timing spectroscopy"", 《SPIE》 * |
| 曹柱荣等: ""Ｘ 射线高空间分辨多色显微成像系统研制及应用"", 《强激光与粒子束》 * |
| 董建军等: ""基于多通道 Kirkpatrick-Baez 显微镜的内爆热斑不对称性实验研究"", 《核聚变与等离子体物理》 * |
| 薛全喜等: ""准分子激光靶面焦斑合束研究分析"", 《中国激光》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110097597B (en) | 2022-02-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gralla et al. | Black hole shadows, photon rings, and lensing rings | |
| Limousin et al. | The three-dimensional shapes of galaxy clusters | |
| Psaltis et al. | A quantitative test of the no-hair theorem with Sgr A* using stars, pulsars, and the event horizon telescope | |
| Chan et al. | The power of imaging: constraining the plasma properties of GRMHD simulations using EHT observations of Sgr A | |
| Schleede et al. | Emphysema diagnosis using X-ray dark-field imaging at a laser-driven compact synchrotron light source | |
| Kingston et al. | Neutron ghost imaging | |
| CN103444194B (en) | Image processing system, image processing apparatus and image processing method | |
| Nguyen et al. | Scattered radiation emission imaging: principles and applications | |
| Shaviv et al. | Magnetic lensing near ultramagnetized neutron stars | |
| Elíasson et al. | Spatial tomography of individual atoms in a quantum gas microscope | |
| KR101268246B1 (en) | Method and apparatus for super-resolution using wobble motion and psf in positron emission tomography | |
| Deliyski et al. | Observing naked singularities by the present and next-generation Event Horizon Telescope | |
| CN110097597A (en) | A kind of coordinate corresponding method of the series of X light image of object | |
| Golovkin et al. | Simulation of X-ray scattering diagnostics in multi-dimensional plasma | |
| Rahman et al. | Task-based assessment of binned and list-mode SPECT systems | |
| Sharma et al. | Design and development of a high-energy gamma camera for use with NSECT imaging: feasibility for breast imaging | |
| Fleenor et al. | Correlated statistical uncertainties in coded-aperture imaging | |
| McNally | Computational optical-sectioning microscopy for 3D quantification of cell motion: results and challenges | |
| JP2018022376A (en) | Information processor, information processing method and program | |
| Caucci et al. | Radiance and photon noise: imaging in geometrical optics, physical optics, quantum optics and radiology | |
| Huang et al. | Images and flares of geodesic hot spots around a Kerr black hole | |
| Lorsakul et al. | 4D numerical observer for lesion detection in respiratory‐gated PET | |
| von Niederhäusern et al. | Augmenting camera images with gamma detector data: a novel approach to support sentinel lymph node biopsy | |
| Barros et al. | Monte Carlo studies and response assessment of a new Rotational Modulation Collimator (RMC) design | |
| JP2000135211A (en) | Radiographic instrument, radiography, memory medium recording program to realize radiography, and transmission medium to be transmitted |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |