CN110702708A - X-ray detection transillumination geometric parameter measuring method - Google Patents
X-ray detection transillumination geometric parameter measuring method Download PDFInfo
- Publication number
- CN110702708A CN110702708A CN201911065064.9A CN201911065064A CN110702708A CN 110702708 A CN110702708 A CN 110702708A CN 201911065064 A CN201911065064 A CN 201911065064A CN 110702708 A CN110702708 A CN 110702708A
- Authority
- CN
- China
- Prior art keywords
- plane
- ray machine
- measuring
- transillumination
- imaging plate
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pulmonology (AREA)
- Radiology & Medical Imaging (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Length-Measuring Devices Using Wave Or Particle Radiation (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The embodiment of the application discloses a method for measuring transillumination geometric parameters in X-ray detection, which comprises the steps of determining the space coordinates of an X-ray machine, detected equipment and an imaging device in a global coordinate system by establishing the global coordinate system, calculating the transillumination geometric parameters of the X-ray machine, the detected equipment and the imaging device during detection according to the space coordinates, and providing a foundation for the application of a CT (computed tomography) technology in the field detection of power equipment.
Description
Technical Field
The invention relates to the field of power equipment detection, in particular to the field of transillumination geometric parameter measurement in an X-ray detection process.
Background
At present, an X-ray digital imaging technology is widely applied to nondestructive testing of power equipment, can visually and conveniently detect common fault defects in the power equipment, and plays a vital role in ensuring the safe and stable operation of the power equipment. The CT technology based on the X-ray detection technology can realize three-dimensional imaging of the internal structure of the equipment, more intuitively and accurately present the detection result, and has important significance in fault and defect diagnosis of the power equipment. In the CT technique, the transillumination geometric parameters are an important basis for the three-dimensional reconstruction of images. Due to the fact that the field power equipment is complex in structure and large in size, transillumination geometric parameters among the X-ray machine, the detected equipment and the imaging device are difficult to accurately measure through conventional means in the detection process.
The invention provides a method for measuring transillumination geometric parameters of X-ray detection, which determines the space coordinates of an X-ray machine, detected equipment and an imaging device in a global coordinate system by establishing the global coordinate system, calculates the transillumination geometric parameters of the X-ray machine, the detected equipment and the imaging device during detection according to the space coordinates, and provides a foundation for the application of a CT (computed tomography) technology in the field detection of power equipment.
Disclosure of Invention
Aiming at the problems, the invention provides a measuring method of X-ray detection transillumination geometric parameters, which determines the space coordinates of an X-ray machine, detected equipment and an imaging device in a global coordinate system by establishing the global coordinate system, calculates the transillumination geometric parameters of the X-ray machine, the detected equipment and the imaging device during detection according to the space coordinates, and provides a foundation for the application of the CT technology in the field detection of power equipment.
The invention provides an X-ray detection transillumination geometric parameter measuring method which comprises an origin positioning device, a distance measurement signal emitting device, a signal receiving device of an imaging plate, a signal receiving device of an X-ray machine, a signal receiving device of an object to be detected and a data analysis processing system.
The measuring steps are as follows:
the method comprises the following steps: selecting a point in a detection field space as an origin of global coordinates, and establishing a global coordinate system;
step two: the imaging plate is a plane model, the distances from three points which are not positioned on a straight line on the plane of the imaging plate to the original point are measured, and the plane position of the imaging plate is determined;
step three: the window of the X-ray machine is a plane model, the distances from three points, which are not positioned on a straight line, of the window of the X-ray machine to the origin are measured, and the position of the window of the X-ray machine is determined;
step four: the electric power equipment is generally a cylinder, the position of the electric power equipment can be determined by determining the axis of the electric power equipment in combination with a drawing of the electric power equipment, the center of a circle can be determined by the axis through three points on an equal altitude surface, and a normal line passing through the center of the circle is the axis.
Step five: after the position equations of the imaging plate, the X-ray machine and the equipment are determined, the geometric parameters during transillumination can be obtained.
Furthermore, the signal emitting device and the signal receiver are connected with a data processing terminal in a limited or wireless mode, and the data processing terminal can convert the signals into distance data and calculate the positions of the imaging plate, the X-ray machine and the detected object according to the distance data.
Further, the specific implementation method for determining the plane position of the imaging plate in the second step includes:
a, B, C three-point coordinates which are not on the same straight line are measured, and are respectively:
A(X1，Y1，Z1)
B(X2，Y2，Z2)
C(X3，Y3，Z3)
normal vector of plane
And in-plane vector ofAnd
is perpendicular, therefore
a＝(y2-y1)*(z3-z1)-(y3-y1)*(z2-z1)
b＝(z2-z1)*(x3-x1)-(z3-z1)*(x2-x1)
c＝(x2-x1)*(y3-y1)-(x3-x1)*(y2-y1)
The plane equation is: a (x-x1) + b (y-y1) + c (z-z1) ═ 0.
Further, the determination of the axis of the cylinder in the fourth step specifically includes:
measuring three-point coordinates A, B, C on a circle with the same height, which are respectively as follows:
A(X1，Y1，Z1)
B(X2，Y2，Z2)
C(X3，Y3，Z3)
normal vector of plane
And in-plane vector of
Andperpendicular, therefore:
a＝(y2-y1)*(z3-z1)-(y3-y1)*(z2-z1)
b＝(z2-z1)*(x3-x1)-(z3-z1)*(x2-x1)
c＝(x2-x1)*(y3-y1)-(x3-x1)*(y2-y1)
let the center coordinate be (X)0，Y0，Z0)；
The three formulas are mutually equal to each other to obtain,
the plane equation is:
wherein:
A1＝y1·z2-y1·z3-z1·y2+z1·y3+y2·z3-y3·z2
B1＝-x1·z2+x1·z3+z1·x2-z1·x3-x2·z3+x3·z2
C1＝x1·y2-x1·y3-y1·x2+y1·x3+x2·y3-x3·y2
D1＝-x1·y2·z3+x1·y3·z2+x2·y1·z3-x3·y1·z2-x2·y3·z1+x3·y2·z1
the simultaneous (4), (5) and (6) can obtain:
the circle center coordinate is obtained by solving:
the equation for the cylinder axis is:
(X-X0)(X1-X0)+(Y-Y0)(Y1-Y0)+(Z-Z0)(Z1-Z0)＝0。
the beneficial effect of this application: the invention provides a method for measuring transillumination geometric parameters of X-ray detection, which determines the space coordinates of an X-ray machine, detected equipment and an imaging device in a global coordinate system by establishing the global coordinate system, calculates the transillumination geometric parameters of the X-ray machine, the detected equipment and the imaging device during detection according to the space coordinates, and provides a foundation for the application of a CT (computed tomography) technology in the field detection of power equipment.
Drawings
In order to more clearly explain the technical solution of the present application, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious to those skilled in the art that other drawings can be obtained according to the drawings without any creative effort.
FIG. 1 is a schematic diagram of a distance measurement principle;
FIG. 2 is a schematic view of a planar position measurement;
FIG. 3 is a schematic view of cylinder axis determination;
1: a ranging signal transmitting device; 2: an imaging plate; 2-1/2-2/2-3: a signal receiver integrated on the imaging plate; 3: a test object; 3-1/3-2/3-3: a signal receiver integrated on the object to be detected; 4: an X-ray machine; 4-1/4-2/4-3: and a signal receiver integrated on the imaging plate.
Detailed Description
The invention provides a method for measuring transillumination geometric parameters of X-ray detection, which determines the space coordinates of an X-ray machine, detected equipment and an imaging device in a global coordinate system by establishing the global coordinate system, calculates the transillumination geometric parameters of the X-ray machine, the detected equipment and the imaging device during detection according to the space coordinates, and provides a foundation for the application of a CT (computed tomography) technology in the field detection of power equipment.
The invention provides a method for measuring X-ray detection transillumination geometric parameters, which comprises an origin positioning device, a distance measurement signal emitting device, a signal receiving device at an imaging plate, a signal receiving device of an X-ray machine, a signal receiving device of a detected object and a data analysis processing system, wherein the origin positioning device, the distance measurement signal emitting device, the signal receiving device at the imaging plate are arranged on the imaging plate.
The measuring steps are as follows:
the method comprises the following steps: selecting a point in a detection field space as an origin of global coordinates, and establishing a global coordinate system;
step two: the imaging plate is a plane model, the distances from three points which are not positioned on a straight line on the plane of the imaging plate to the original point are measured, and the plane position of the imaging plate is determined;
step three: the window of the X-ray machine is a plane model, the distances from three points, which are not positioned on a straight line, of the window of the X-ray machine to the origin are measured, and the position of the window of the X-ray machine is determined;
step four: the electric power equipment is generally a cylinder, the position of the electric power equipment can be determined by determining the axis of the electric power equipment in combination with a drawing of the electric power equipment, the center of a circle can be determined by the axis through three points on an equal altitude surface, and a normal line passing through the center of the circle is the axis.
Step five: after the position equations of the imaging plate, the X-ray machine and the equipment are determined, the geometric parameters during transillumination can be obtained.
Furthermore, the signal emitting device and the signal receiver are connected with a data processing terminal in a limited or wireless mode, and the data processing terminal can convert the signals into distance data and calculate the positions of the imaging plate, the X-ray machine and the detected object according to the distance data.
Further, as shown in fig. 2, the specific implementation method for determining the plane position of the imaging plate in the second step includes:
a, B, C three-point coordinates which are not on the same straight line are measured, and are respectively:
A(X1，Y1，Z1)
B(X2，Y2，Z2)
C(X3，Y3，Z3)
normal vector of plane
And in-plane vector ofAnd
is perpendicular, therefore
a＝(y2-y1)*(z3-z1)-(y3-y1)*(z2-z1)
b＝(z2-z1)*(x3-x1)-(z3-z1)*(x2-x1)
c＝(x2-x1)*(y3-y1)-(x3-x1)*(y2-y1)
The plane equation is: a (x-x1) + b (y-y1) + c (z-z1) ═ 0.
Further, as shown in fig. 3, the determination of the cylinder axis in the fourth step specifically includes:
measuring three-point coordinates A, B, C on a circle with the same height, which are respectively as follows:
A(X1，Y1，Z1)
B(X2，Y2，Z2)
C(X3，Y3，Z3)
normal vector of plane
And in-plane vector of
And
perpendicular, therefore:
a＝(y2-y1)*(z3-z1)-(y3-y1)*(z2-z1)
b＝(z2-z1)*(x3-x1)-(z3-z1)*(x2-x1)
c＝(x2-x1)*(y3-y1)-(x3-x1)*(y2-y1)
let the center coordinate be (X)0，Y0，Z0)；
The three formulas are mutually equal to each other to obtain,
the plane equation is:
wherein:
A1＝y1·z2-y1·z3-z1·y2+z1·y3+y2·z3-y3·z2
B1＝-x1·z2+x1·z3+z1·x2-z1·x3-x2·z3+x3·z2
C1＝x1·y2-x1·y3-y1·x2+y1·x3+x2·y3-x3·y2
D1＝-x1·y2·z3+x1·y3·z2+x2·y1·z3-x3·y1·z2-x2·y3·z1+x3·y2·z1
the simultaneous (4), (5) and (6) can obtain:
the circle center coordinate is obtained by solving:
the equation for the cylinder axis is:
(X-X0)(X1-X0)+(Y-Y0)(Y1-Y0)+(Z-Z0)(Z1-Z0)＝0。
other embodiments of the invention will be apparent to those skilled in the art from consideration of the specification and practice of the invention disclosed herein. This application is intended to cover any variations, uses, or adaptations of the invention following, in general, the principles of the invention and including such departures from the present disclosure as come within known or customary practice within the art to which the invention pertains. It is intended that the specification and examples be considered as exemplary only, with a true scope and spirit of the invention being indicated by the following claims.
It will be understood that the invention is not limited to the precise arrangements described above and shown in the drawings and that various modifications and changes may be made without departing from the scope thereof. The scope of the invention is limited only by the appended claims.
Claims (4)
1. A method for measuring X-ray detection transillumination geometric parameters is characterized by comprising an origin positioning device, a distance measurement signal transmitting device, a signal receiving device and a data analysis processing system;
the measuring steps are as follows:
the method comprises the following steps: the origin positioning device selects one point of a detection field space as a global coordinate origin and establishes a global coordinate system;
step two: the distance measuring signal transmitting device and the distance measuring signal receiving device measure the distance from three points which are not on a straight line on the plane of the imaging plate to the original point, and the plane position of the imaging plate is determined;
step three: the window of the X-ray machine is a plane model, the distance measuring signal emitting device and the distance measuring signal receiving device measure the distance from three points on a straight line, which is not positioned on the window of the X-ray machine, to the origin, and the position of the window of the X-ray machine is determined;
step four: the electric equipment is generally a cylinder, the position of the electric equipment can be determined by determining the axis of the electric equipment in combination with a drawing of the electric equipment, the center of a circle can be determined by the axis of the electric equipment through three points on an equal altitude surface, and a normal line passing through the center of the circle is the axis;
step five: after the data analysis processing system determines the position equation of the imaging plate, the X-ray machine and the equipment, the geometric parameters during transillumination can be obtained.
2. The method according to claim 1, wherein the signal transmitter and the signal receiver are connected with a data processing terminal in a limited or wireless manner, and the data processing terminal can convert the signals into distance data and calculate the positions of the imaging plate, the X-ray machine and the object to be detected.
3. The method for measuring X-ray detection transillumination geometric parameters according to claim 1, wherein the specific implementation method for determining the planar position of the imaging plate in the second step is as follows:
a, B, C three-point coordinates which are not on the same straight line are measured, and are respectively:
A(X1，Y1，Z1)
B(X2，Y2，Z2)
C(X3，Y3，Z3)
normal vector of plane
And in-plane vector of
And
is perpendicular, therefore
a＝(y2-y1)*(z3-z1)-(y3-y1)*(z2-z1)
b＝(z2-z1)*(x3-x1)-(z3-z1)*(x2-x1)
c＝(x2-x1)*(y3-y1)-(x3-x1)*(y2-y1)
The plane equation is: a (x-x1) + b (y-y1) + c (z-z1) ═ 0.
4. The method for measuring trans-illumination geometric parameters in X-ray detection according to claim 1, wherein the determination of the cylinder axis in the fourth step is specifically as follows:
measuring three-point coordinates A, B, C on a circle with the same height, which are respectively as follows:
A(X1，Y1，Z1)
B(X2，Y2，Z2)
C(X3，Y3，Z3)
normal vector of planeAnd in-plane vector of
And
perpendicular, therefore:
a＝(y2-y1)*(z3-z1)-(y3-y1)*(z2-z1)
b＝(z2-z1)*(x3-x1)-(z3-z1)*(x2-x1)
c＝(x2-x1)*(y3-y1)-(x3-x1)*(y2-y1)
let the center coordinate be (X)0，Y0，Z0)；
The three formulas are mutually equal to each other to obtain,
is recorded as: a. the2x+B2y+C2z+D2＝0 (5)
Is recorded as: a. the3 x+B3 y+C3 z+D3＝0 (6)
The plane equation is:
wherein:
A1＝y1·z2-y1·z3-z1·y2+z1·y3+y2·z3-y3·z2
B1＝-x1·z2+x1·z3+z1·x2-z1·x3-x2·z3+x3·z2
C1＝x1·y2-x1·y3-y1·x2+y1·x3+x2·y3-x3·y2
D1＝-x1·y2·z3+x1·y3·z2+x2·y1·z3-x3·y1·z2-x2·y3·z1+x3·y2·z1
the simultaneous (4), (5) and (6) can obtain:
the circle center coordinate is obtained by solving:
the equation for the cylinder axis is:
(X-X0)(X1-X0)+(Y-Y0)(Y1-Y0)+(Z-Z0)(Z1-Z0)＝0。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911065064.9A CN110702708A (en) | 2019-11-04 | 2019-11-04 | X-ray detection transillumination geometric parameter measuring method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911065064.9A CN110702708A (en) | 2019-11-04 | 2019-11-04 | X-ray detection transillumination geometric parameter measuring method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110702708A true CN110702708A (en) | 2020-01-17 |
Family
ID=69204125
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911065064.9A Pending CN110702708A (en) | 2019-11-04 | 2019-11-04 | X-ray detection transillumination geometric parameter measuring method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110702708A (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6163589A (en) * | 1998-06-13 | 2000-12-19 | General Electric Company | Monte Carlo scatter correction method for computed tomography of general object geometries |
| CN102743184A (en) * | 2012-05-14 | 2012-10-24 | 清华大学 | Geometrical parameter calibration method of X-ray cone beam computed tomography system |
| CN104280013A (en) * | 2014-10-30 | 2015-01-14 | 中国电建集团成都勘测设计研究院有限公司 | Method for determining attitude of rock mass structural plane based on measurement coordinates |
| CN104665862A (en) * | 2015-02-16 | 2015-06-03 | 清华大学 | Method for eliminating geometrical artifacts from CBCT (Cone Beam Computed Tomography) and CBCT system using method for eliminating geometrical artifacts from CBCT |
| CN105783811A (en) * | 2016-04-26 | 2016-07-20 | 中国海洋石油总公司 | Detection method of circular tube end part center three-dimensional coordinate |
| DE102015207102A1 (en) * | 2015-04-20 | 2016-10-20 | Siemens Aktiengesellschaft | X-ray computed tomography method |
| CN107389705A (en) * | 2017-07-20 | 2017-11-24 | 云南电网有限责任公司电力科学研究院 | A kind of power equipment X ray Non-Destructive Testing three-D imaging method |
-
2019
- 2019-11-04 CN CN201911065064.9A patent/CN110702708A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6163589A (en) * | 1998-06-13 | 2000-12-19 | General Electric Company | Monte Carlo scatter correction method for computed tomography of general object geometries |
| CN102743184A (en) * | 2012-05-14 | 2012-10-24 | 清华大学 | Geometrical parameter calibration method of X-ray cone beam computed tomography system |
| CN104280013A (en) * | 2014-10-30 | 2015-01-14 | 中国电建集团成都勘测设计研究院有限公司 | Method for determining attitude of rock mass structural plane based on measurement coordinates |
| CN104665862A (en) * | 2015-02-16 | 2015-06-03 | 清华大学 | Method for eliminating geometrical artifacts from CBCT (Cone Beam Computed Tomography) and CBCT system using method for eliminating geometrical artifacts from CBCT |
| DE102015207102A1 (en) * | 2015-04-20 | 2016-10-20 | Siemens Aktiengesellschaft | X-ray computed tomography method |
| CN105783811A (en) * | 2016-04-26 | 2016-07-20 | 中国海洋石油总公司 | Detection method of circular tube end part center three-dimensional coordinate |
| CN107389705A (en) * | 2017-07-20 | 2017-11-24 | 云南电网有限责任公司电力科学研究院 | A kind of power equipment X ray Non-Destructive Testing three-D imaging method |
Non-Patent Citations (1)
| Title |
|---|
| 连保胜 等: "《新编高等数学》", 31 August 2019 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6193657B1 (en) | Image based probe position and orientation detection | |
| JP3871747B2 (en) | Ultrasonic diagnostic equipment | |
| US9123096B2 (en) | Information processing apparatus and control method thereof | |
| CN108430334B (en) | Ultrasound imaging apparatus and ultrasound imaging method for examining a volume of an object | |
| US8179132B2 (en) | Method and system for integrating eddy current inspection with a coordinate measuring device | |
| US10222352B2 (en) | Method and device for improving the SAFT analysis when measuring irregularities | |
| CN102626324A (en) | Ultrasound measuring apparatus and control method thereof | |
| CN104655658A (en) | Large-sized high-temperature blade internal defect three-dimensional nondestructive detection method | |
| CN107255673A (en) | High temperature blade internal flaw three dimensional lossless detection method based on ultrasonic phase array | |
| CN108267741A (en) | A kind of ultrasonic probe caliberating device and the method for demarcating ultrasonic probe | |
| CN105954359B (en) | Complicated shape inside parts defect distribution formula ultrasonic no damage detection device and method | |
| CN107389705A (en) | A kind of power equipment X ray Non-Destructive Testing three-D imaging method | |
| CN106842241A (en) | A kind of high dynamic satellite navigation receiver performance estimating method | |
| KR100859434B1 (en) | Image processing system and method for providing contour editing function using multiple sectional plane images | |
| WO2020113787A1 (en) | Ultrasound probe calibration method | |
| CN110702708A (en) | X-ray detection transillumination geometric parameter measuring method | |
| CN218443999U (en) | Ultrasonic flow detection system | |
| US10578459B2 (en) | Method for operating an x-ray device | |
| CN116576806B (en) | Precision control system for thickness detection equipment based on visual analysis | |
| JP4099196B2 (en) | Ultrasonic diagnostic equipment | |
| CN107951512A (en) | A kind of method and apparatus of generation fetal weight for ultrasound scanning device | |
| JP4141466B2 (en) | Ultrasonic diagnostic equipment | |
| JPS595950A (en) | Ultrasonic diagnostic method | |
| Tamashevich | Advanced OST System for the Secondsound Test of Fully Dressed Cavities | |
| KR20230102590A (en) | Ultrasound scanner and method of guiding aim |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200117 |