CN110097597B - Coordinate corresponding method for series X-ray images of target object - Google Patents
Coordinate corresponding method for series X-ray images of target object Download PDFInfo
- Publication number
- CN110097597B CN110097597B CN201910376664.0A CN201910376664A CN110097597B CN 110097597 B CN110097597 B CN 110097597B CN 201910376664 A CN201910376664 A CN 201910376664A CN 110097597 B CN110097597 B CN 110097597B
- Authority
- CN
- China
- Prior art keywords
- ray
- series
- ray image
- ray images
- 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.)
- Active
Links
Images
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 coordinate corresponding method of a series of X-ray images of a target object. The method comprises the following steps: a. recording a series of X-ray images of a target object, summing the X-ray intensities recorded by the recording units in each X-ray image to obtain the total intensity of each X-ray to obtain A; b. multiplying the X-ray intensity of each recording unit in each X-ray image by the corresponding two-dimensional coordinate, and summing to obtain B; c.C = B/A, and the light intensity barycentric coordinates of each X-ray image are obtained; d. and giving a two-dimensional coordinate corresponding relation of each X-ray image in the series of X-ray images by utilizing the attribute corresponding to the light intensity barycentric coordinate of each X-ray image in the series of X-ray images. The method can realize the accurate correspondence of the coordinates of each X-ray image in different observation directions, quantitatively analyze the time evolution of the coordinates of each X-ray image in different time periods, and ensure that the coordinate correspondence precision reaches the spatial resolution level of an imaging system, thereby having wide and important application prospect in inertial confinement fusion, laboratory celestial body physics or high-energy density physics.
Description
Technical Field
The invention belongs to the field of X-ray image data processing, and particularly relates to a coordinate corresponding method of a series of X-ray images of a target object.
Background
X-ray imaging is very widely used in inertial confinement fusion, laboratory celestial physics, and high energy density physical scientific research for diagnosing two-dimensional spatial, time-resolved information of three-dimensional targets, such as fusion shots, radiation sources, and substances under extreme conditions. In general, in two-dimensional spatially resolved diagnosis, one obtains a plurality of X-ray images of different observation positions, and the above scientific research needs to accurately correspond the images to one another through coordinates, because important information such as electron temperature, electron density, symmetry analysis and the like can be given through the correspondence. In time-resolved diagnosis, people can obtain a plurality of X-ray images which change along with time, and related scientific researches need to quantitatively analyze the coordinate time evolution of the images so as to research the problems of symmetry change and the like.
The existing X-ray image data processing has the following defects: 1. for the space resolution X-ray imaging diagnosis of different observation directions, according to the projection principle and the non-uniformity of target light emission, the obtained X-ray images of different observation directions are obviously different, which comprises the difference of space intensity distribution and appearance geometry, so that the X-ray images are difficult to realize the accurate correspondence of coordinates; 2. for time-resolved X-ray imaging diagnosis at different time periods, since the X-ray light state of the target object changes along with the time evolution, the obtained X-ray images at different time periods are also obviously different, which also includes the difference of space intensity distribution and shape geometry, but how to quantitatively analyze the coordinate time evolution of the images still lacks a corresponding coordinate corresponding method.
Disclosure of Invention
The invention aims to provide a coordinate corresponding method of a series of X-ray images of a target object.
The coordinate correspondence method of the series of X-ray images of the target object comprises the following steps:
a. recording a series of X-ray images of a target object, summing the X-ray intensities recorded by the recording units in each X-ray image to obtain the total intensity of each X-ray to obtain A;
b. multiplying the X-ray intensity of each recording unit in each X-ray image by the corresponding two-dimensional coordinate, and summing to obtain B;
c.C = B/A, and the light intensity barycentric coordinates of each X-ray image are obtained;
d. and giving a two-dimensional coordinate corresponding relation of each X-ray image in the series of X-ray images by utilizing the attribute corresponding to the light intensity barycentric coordinate of each X-ray image in the series of X-ray images.
The X-ray image is formed by one of a pinhole, a Kirkpatrick-Baez mirror or a curved crystal imaging device.
The series of X-ray images comprise X-ray images recorded at different positions aiming at the same target object, or X-ray images recorded at different time periods.
The coordinate corresponding method of the series X-ray images of the target object is applied to the data processing of the X-ray images of the target object in inertial confinement fusion, laboratory celestial body physics or high-energy density physics.
The coordinate corresponding method of the series X-ray images of the target object of the invention comprises the steps of obtaining the total intensity A of each X-ray, the sum B of the product of the X-ray intensity of each recording unit in each X-ray image and the corresponding two-dimensional coordinate, and the corresponding relation between the light intensity barycentric coordinate C of each X-ray image and the two-dimensional coordinate of each X-ray image through manual recording, software or programming.
The coordinate corresponding method of the series X-ray images of the target object has the following advantages:
1. the accurate correspondence of the coordinates of each X-ray image in different observation directions can be realized.
2. The time evolution of the coordinates of each X-ray image in different time periods can be quantitatively analyzed.
3. The coordinate corresponding precision reaches the spatial resolution level of the imaging system.
The coordinate correspondence method of the series X-ray images of the target object can realize the accurate correspondence and quantitative analysis of the coordinates of the series X-ray images, and has wide and important application prospects in inertial confinement fusion, laboratory celestial body physics or high-energy density physics.
Drawings
Fig. 1 is a projection coordinate diagram of a coordinate mapping method of a series of X-ray images of an object according to the present invention.
Detailed Description
The invention is described in detail below with reference to the figures and specific embodiments.
The coordinate correspondence method of the series of X-ray images of the target object comprises the following steps:
a. recording a series of X-ray images of a target object, summing the X-ray intensities recorded by the recording units in each X-ray image to obtain the total intensity of each X-ray to obtain A;
b. multiplying the X-ray intensity of each recording unit in each X-ray image by the corresponding two-dimensional coordinate, and summing to obtain B;
c.C = B/A, and the light intensity barycentric coordinates of each X-ray image are obtained;
d. and giving a two-dimensional coordinate corresponding relation of each X-ray image in the series of X-ray images by utilizing the attribute corresponding to the light intensity barycentric coordinate of each X-ray image in the series of X-ray images.
The X-ray image is formed by one of a pinhole, a Kirkpatrick-Baez mirror or a curved crystal imaging device.
The series of X-ray images comprise X-ray images recorded at different positions aiming at the same target object, or X-ray images recorded at different time periods.
The coordinate corresponding method of the series X-ray images of the target object is applied to the data processing of the X-ray images of the target object in inertial confinement fusion, laboratory celestial body physics or high-energy density physics.
Example 1
The X-ray image is formed by a pinhole. The series of X-ray images are recorded in different directions aiming at the same target object. The coordinate corresponding method of the series X-ray images of the target object is applied to the data processing of the X-ray images of the target object in inertial confinement fusion.
According to basic physical knowledge, a mass-integrated two-dimensional projection is carried out on a mass object to different directions, and the gravity centers of the mass object are fixed on each projection plane and correspond to each other one by one. I.e. the centre of gravity of the three-dimensional object must pass the centre of gravity of the two-dimensional projection in the projection direction. If the object mass is exchanged for the intensity of the X-ray source, the conclusion is consistent, i.e. the center of gravity of the intensity of the X-ray source of the object must pass through the center of gravity of the intensity of the two-dimensional projection in the projection direction (viewing direction).
The light intensity distribution of the X-ray source of the target object isSee FIG. 1, along the observation directionIs projected as
Wherein the content of the first and second substances,is thatIn thatProjection of a plane andthe included angle of the axes is set by the angle,is thatAndthe included angle of the plane is formed by the angle,is a rotation matrix of three-dimensional coordinates
Let the X-ray source of the object have the light intensity barycentric coordinate ofThe light intensity barycentric coordinate of each projected X-ray image is CAccording to the definition of the center of gravity
T is the intensity of the X-ray source of the targetIs the total intensity of the X-ray source of the three-dimensional object. The light intensity barycentric coordinate C of each X-ray image projected (observed)
Wherein the content of the first and second substances,the sum B of the products of the X-ray intensities of the recording elements in each X-ray image and the corresponding two-dimensional coordinates,for projecting (observing) the intensity of X-ray imagesI.e. the total intensity a of each X-ray projected (observed).The relationship with T is
WhereinR is a constant related to R for the solid angle occupied by the projection (observation) X-ray image to the target X-ray source.
Is thatObtained by two rotations, the relationship between the two is still determined by formula (2), i.e.
Where S is the matrix formed by the first and second rows of the matrix K, i.e.
Therefore, the total intensity A of each X-ray is extracted from each X-ray imageAnd then extracting the sum B of the product of the X-ray intensity of each recording unit in each X-ray image and the corresponding two-dimensional coordinate, thus obtaining the light intensity barycentric coordinate C of each X-ray image through the formula (5). And (4) giving a two-dimensional coordinate corresponding relation of each X-ray image in the series of X-ray images according to a formula (7), namely, by utilizing the attribute corresponding to the light intensity barycentric coordinate of each X-ray image in the series of X-ray images. The coordinate correspondence accuracy thus obtained is determined by the spatial resolution of the pinhole imaging system, and is of the order of 10 μm.
In addition, after the coordinate corresponding relation is determined, quantitative research of various X-ray images with energy spectrum resolution is carried out, so that the electron temperature and the electron density can be extracted according to the X-ray radiation property, such as a common bremsstrahlung mechanism. And evaluating the three-dimensional symmetry of the X-ray emission of the target object according to all coordinate corresponding conditions.
Example 2
This example is the same as the example 1 except that the X-ray image was formed by a Kirkpatrick-Baez mirror. The series of X-ray images are recorded aiming at the same target object in different time periods. The coordinate correspondence method of the series X-ray images of the target object is applied to the data processing of the X-ray images of the target object in high-energy density physics. After the light intensity gravity center coordinate C of each X-ray image is obtained through the formula (5), the time evolution of the coordinates of each X-ray image in different time periods can be quantitatively analyzed by comparing the distance change between the light intensity gravity center and the geometric center of each X-ray image, and therefore the symmetrical change trend of the target object is researched. The coordinate correspondence accuracy thus obtained is determined by the spatial resolution of the Kirkpatrick-Baez mirror imaging system, of the order of 2.5 μm.
The embodiment can also be applied to the X-ray image data processing of the target object in the physical treatment of the celestial body in a laboratory.
It will be appreciated by those of ordinary skill in the art that the embodiments described herein are intended to assist the reader in understanding the principles of the invention and are to be construed as being without limitation to such specifically recited embodiments and examples. Those skilled in the art can make various other specific changes and combinations based on the teachings of the present invention without departing from the spirit of the invention, and these changes and combinations are within the scope of the invention.
Claims (4)
1. A method for mapping coordinates of a series of X-ray images of an object, said method comprising the steps of:
a. recording a series of X-ray images of a target object, summing the X-ray intensities recorded by the recording units in each X-ray image to obtain the total intensity of each X-ray to obtain A;
b. multiplying the X-ray intensity of each recording unit in each X-ray image by the corresponding two-dimensional coordinate, and summing to obtain B;
c.C = B/A, and the light intensity barycentric coordinates of each X-ray image are obtained;
d. and giving a two-dimensional coordinate corresponding relation of each X-ray image in the series of X-ray images by utilizing the attribute corresponding to the light intensity barycentric coordinate of each X-ray image in the series of X-ray images.
2. The method of claim 1, wherein the X-ray image is an X-ray image formed by one of a pinhole, a Kirkpatrick-Baez mirror, or a curved crystal imaging device.
3. The method of claim 1, wherein the series of X-ray images includes X-ray images recorded at different orientations or time periods for the same object.
4. The coordinate mapping method of the series of X-ray images of the target object according to any one of claims 1 to 3, wherein the coordinate mapping method is applied to the data processing of the X-ray images of the target object in inertial confinement fusion, laboratory celestial physics or high energy density physics.
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 CN110097597A (en) | 2019-08-06 |
CN110097597B true 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 (15)
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 |
CN102530850A (en) * | 2012-03-14 | 2012-07-04 | 哈尔滨工业大学 | Method for millimeter-sized micro nanostructure nano carving and processing through adopting antifrictional metal (AFM) needle |
CN102542600A (en) * | 2011-12-14 | 2012-07-04 | 北京工业大学 | Simulated projection DRR( digitally reconstructed radiograph) generating method based on CUDA (compute unified device architecture) technology |
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 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8244017B2 (en) * | 2010-02-08 | 2012-08-14 | James Jiwen Chun | Constructing three dimensional images using panoramic images |
-
2019
- 2019-05-05 CN CN201910376664.0A patent/CN110097597B/en active Active
Patent Citations (15)
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 |
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 |
---|
"The Neutron star Interior Composition ExploreR (NICER): an Explorer mission of opportunity for soft x-ray timing spectroscopy";Keith C. Gendreau等;《SPIE》;20160722;正文第1-8页 * |
"Three-dimensional strain mapping using in situ X-ray synchrotron microtomography";H Toda1等;《SPECIAL ISSUE PAPER》;20110729;第549-561页 * |
"Understanding X-ray reflection emissivity profiles in AGN: locating the X-ray source";D. R. Wilkins等;《X-ray reflection emissivity profiles in AGN》;20121231;第1284–1296页 * |
"X 射线高空间分辨多色显微成像系统研制及应用";曹柱荣等;《强激光与粒子束》;20150331;第27卷(第3期);正文第1-7页 * |
"准分子激光靶面焦斑合束研究分析";薛全喜等;《中国激光》;20140331;第41卷(第3期);正文第1-15页 * |
"基于多通道 Kirkpatrick-Baez 显微镜的内爆热斑不对称性实验研究";董建军等;《核聚变与等离子体物理》;20180630;第38卷(第2期);第125-129页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110097597A (en) | 2019-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6683206B2 (en) | Oblique CT device | |
US10660580B2 (en) | Directed X-ray fields for tomosynthesis | |
US20180329225A1 (en) | Pattern Detection at Low Signal-To-Noise Ratio | |
CN107709968A (en) | Fluorescence biopsy sample imager and method | |
US8345817B2 (en) | Apparatus and method for non-rotational computer tomography | |
KR20150079560A (en) | Geometric characterization and calibration of a cone-beam computer tomography apparatus | |
KR100576228B1 (en) | A multiple view angle x-ray stereoscopic imaging method and system | |
WO1997034995A1 (en) | Analytical imaging system and process | |
WO2015129642A1 (en) | Flow velocimetry method for visualized fluid, and flow velocimetry system | |
Le Besnerais et al. | Fast and dense 2D and 3D displacement field estimation by a highly parallel image correlation algorithm | |
JP2014123124A (en) | Method for high-resolution 3d-localization microscopy | |
Stenz et al. | High-precision 3D object capturing with static and kinematic terrestrial laser scanning in industrial applications—approaches of quality assessment | |
Grubišić et al. | Active 3D scanning based 3D thermography system and medical applications | |
Huang et al. | Improved visual inspection through 3D image reconstruction of defects based on the photometric stereo technique | |
Sansoni et al. | Feasibility of contactless 3D optical measurement for the analysis of bone and soft tissue lesions: new technologies and perspectives in forensic sciences | |
CN110097597B (en) | Coordinate corresponding method for series X-ray images of target object | |
Jiang et al. | Calibration and uncertainty analysis of a combined tracking-based vision measurement system using Monte Carlo simulation | |
CN103126710A (en) | High time resolution three-dimensional imaging method based on framing camera | |
Nakasako | X-ray diffraction imaging of biological cells | |
Jiang et al. | Calibration method for binocular vision with large FOV based on normalized 1D homography | |
McGray et al. | Centroid precision and orientation precision of planar localization microscopy | |
US20050047637A1 (en) | Image buffers and access schedules for image reconstruction systems | |
Bureev et al. | Digital X-ray tomography | |
Blanton et al. | Extragalactic jets as probes of distant clusters of galaxies and the clusters occupied by bent radio AGN (COBRA) survey | |
JP7008325B2 (en) | Radiation fluoroscopy non-destructive inspection method and radiation fluoroscopy non-destructive inspection equipment |
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 |