CN115024740A - Virtual radiation field display method for common X-ray photography - Google Patents
Virtual radiation field display method for common X-ray photography Download PDFInfo
- Publication number
- CN115024740A CN115024740A CN202210962947.5A CN202210962947A CN115024740A CN 115024740 A CN115024740 A CN 115024740A CN 202210962947 A CN202210962947 A CN 202210962947A CN 115024740 A CN115024740 A CN 115024740A
- Authority
- CN
- China
- Prior art keywords
- beam limiter
- plane
- ray
- point
- field
- 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
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000005855 radiation Effects 0.000 title claims abstract description 22
- 230000002452 interceptive effect Effects 0.000 claims abstract description 25
- 239000011159 matrix material Substances 0.000 claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000003287 optical effect Effects 0.000 claims description 5
- 238000005096 rolling process Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 238000002601 radiography Methods 0.000 claims description 4
- 238000013329 compounding Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 abstract description 4
- 238000003384 imaging method Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 4
- 230000002238 attenuated effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000004804 winding Methods 0.000 description 2
- 206010010356 Congenital anomaly Diseases 0.000 description 1
- 101100410811 Mus musculus Pxt1 gene Proteins 0.000 description 1
- 230000003190 augmentative effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000005865 ionizing radiation Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/54—Control of apparatus or devices for radiation diagnosis
- A61B6/542—Control of apparatus or devices for radiation diagnosis involving control of exposure
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/54—Control of apparatus or devices for radiation diagnosis
- A61B6/547—Control of apparatus or devices for radiation diagnosis involving tracking of position of the device or parts of the device
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The application discloses a virtual radiation field display method for common X-ray photography, which comprises an interactive panel terminal device for displaying a shot picture and a virtual radiation field, a depth camera for shooting a real-time picture, interactive software and back-end processing software which are installed on the interactive panel terminal device, wherein the interactive panel terminal device comprises the step of acquiring any point of a plane where a known beam limiter opening is locatedv n And converting the coordinate into a portal plane pixel coordinate, and displaying through the interactive panel terminal equipment. According to the invention, the homogeneous conversion matrix from the camera coordinate to the image coordinate and the homogeneous conversion matrix from the world coordinate to the camera coordinate are obtained by combining the internal reference and the external reference of the depth camera, so that the size and the shape of the radiation field of each X-ray passing through the beam limiter opening on the plane where the SOD is located can be accurately calculated on the premise of knowing the shape, the size and the SID of the beam limiter opening, the relative position relation between the radiation field and the shot object is displayed in real time through the interactive panel terminal equipment, and the aim of accurate alignment shooting is finally achieved.
Description
Technical Field
The invention relates to the field of computers, in particular to an image data processing method based on electric digital data processing, which belongs to the augmented reality technology, and specifically relates to a virtual field display method for common X-ray photography.
Background
In a common radiography system, a beam limiter is positioned between an X-ray source and an irradiated object, and has the main functions of: 1. the range and the shape of the X-ray are limited, and unnecessary ionizing radiation is prevented from being applied to an irradiated object; 2. blocking scattered rays; 3. the size, shape and center of the X-ray irradiation field (hereinafter referred to as the radiation field) are marked by using LED projection light field (hereinafter referred to as the light field) and black cross projection.
In an actual X-ray photographing process, an operator (hereinafter, referred to as a technician) of the X-ray photographing system needs to adjust the size of the opening of the beam limiter before photographing each time so as to adapt to different examination items. Since the SID (distance of X-ray source to image receptor) of each X-ray photograph is uncertain, the height and thinness of each patient are uncertain, and the technician can usually determine the approximate field of view only by observing the LED field of the beam limiter, which is a time-consuming and labor-consuming process that requires the technician to adjust the beam limiter opening back and forth in and out of the studio; due to the technical limitation of the LED light source, the edge of the light field may become blurred as the shooting distance increases, the brightness of the light field may also be attenuated, and the black cross line may also become thick and blurred, so that the technician may misjudge the size and the center of the current light field.
The virtual radiation field display method is provided for accurately judging the shooting position and the radiation field position and accurately shooting the X-ray film of the interest area, so that the part of the shot object can be synchronously displayed with the virtual radiation field in real time, a technician does not need to repeatedly go between an observation object and an X-ray machine for adjustment, the precision is higher, the convenience is better, and the shooting efficiency is higher.
Disclosure of Invention
In order to solve the technical problems that the irradiation field is visually indicated by an LED light source in the prior X-ray photography technology, such as the irradiation field cannot be clearly identified due to the influence of factors such as the change of the shooting distance and the like, and the irradiation field needs to be adjusted repeatedly, the application provides a common X-ray photography virtual irradiation field display party which is used for displaying the shot object and the irradiation field through an interactive panel terminal in real time, avoids the trouble that a technician needs to repeatedly check the relative position of a light spot of the LED light source on the shot object to determine the position of the irradiation field, can completely cancel the arrangement of the LED light source, avoids the desynchrony of the light field displayed by the LED light source and the actual irradiation field and the congenital error caused by the desynchrony, and further ensures that the judgment of the technician before the shooting is more consistent with the actual shooting condition and is accurate.
In order to achieve the purpose, the application adopts the following specific technical scheme:
the virtual radiation field display method for common X-ray photography comprises interactive panel terminal equipment for displaying a photographed picture and a virtual radiation field, a depth camera for photographing a real-time picture, interactive software and back-end processing software which are installed on the interactive panel terminal equipment, and comprises the following steps:
STP100, establishing a world coordinate system with the X-ray source as an origin O through back-end processing software, and acquiring any point on the opening plane of the beam limiterv n Has world coordinates of
Passing through any one pointv n Spatial rotation matrix occurring relative to originTCalculating to obtain the rotated point
World coordinate
Step STP200, decision point
If the current position is within the known opening range of the beam limiter, if the judgment result is yes, the step STP300 is carried out, if the judgment result is no, the step STP100 is repeated, and the next point is judgedv n+1 Until all points within the range of the known beam limiter opening are traversedv n To therein, whereinn≥1；
Step STP300, according to the points
Obtaining rays
Intersection point with field planep v World coordinate ofp vw 
；
Step STP400, by reading the external reference torque of the depth cameraK 2 Calculating to obtain an intersection pointp v Camera coordinates ofp vc And then obtaining the internal reference torque of the depth camera through back-end processing softwareK 1 Calculating to obtain an intersection pointp v Image coordinates ofp vI Torque of external referenceK 2 The torque is expressed as follows:
wherein,
is a matrix of rotations of the optical system,
is a translation vector;
internal reference torqueK 1 The torque is expressed as follows:
wherein,
is the focal length of the camera and,
and
respectively the ratio of the actual object distance to the pixel in both directions,
is a translation vector from the origin of coordinates of the depth camera to the origin of coordinates of the image;
step STP500, using interactive software to coordinate the image obtained in step STP400p vI Finally displayed on the interactive panel terminal,
wherein SED is the distance from the X-ray source to the opening plane of the beam limiter;
the opening plane of the beam limiter is a plane on the beam limiter, which is perpendicular to the X central ray;
the opening range of the beam limiter refers to the maximum cross-sectional range of the beam limiter for allowing X-rays to pass through;
the field plane is a plane passing through the field center and parallel to the image receiving surface;
the field is the area formed by the intersection of all the X-rays passing through the beam limiter opening and the field plane.
Furthermore, the simplest method for obtaining the rotation matrix R is to establish a trigonometric function relationship of an X-ray shooting scene through data of an IMU accelerometer, which is also the most common calculation method, but this method is not the optimal and most efficient method in practical application; in order to optimize the calculation accuracy of the rotation matrix R, the rotation matrix R in step STP100 is obtained by:
step STP110, passing any point
The spatial rotation generated by the X-ray is divided into a roll angle alpha, a pitch angle beta and a yaw angle gamma;
STP120, respectively calculating a rolling angle alpha, a pitch angle beta and a yaw angle gamma;
based on the X-ray photography principle, if the change of the roll angle alpha is not involved, the roll angle alpha is 0;
the pitch angle β is calculated as follows:
the formula for calculating the yaw angle γ is as follows:
step STP130, the rotation matrices for pitch and yaw are calculated separately as follows:
step STP140, compounding an X-ray rotation matrix R by a pitch angle beta and a yaw angle gamma as follows:
wherein rolling means winding
The rotary movement of the axis, pitch being about
The rotary movement of the shaft, yaw being about
The rotational movement of the shaft is such that,
、
、
respectively, three components of an IMU accelerometer mounted within the beam limiter.
Still further, the camera coordinates in step STP400p vc The calculation method of (c) is as follows:
p
vc
image coordinatesp vI The calculation method is as follows:
p vI 
。
has the advantages that:
according to the invention, the homogeneous conversion matrix from the camera coordinate to the image coordinate and the homogeneous conversion matrix from the world coordinate to the camera coordinate are obtained by combining the internal reference and the external reference of the depth camera, so that the size and the shape of the radiation field of each X-ray passing through the beam limiter opening on the plane where the SOD is located can be accurately calculated on the premise of knowing the shape, the size and the SID of the beam limiter opening, the relative position relation between the radiation field and the shot object is displayed in real time through the interactive panel terminal equipment, and the problem of accurate alignment shooting is finally achieved.
According to the invention, the existing LED light source is not required to be used as visible light indication, so that the error influence of the existing judgment by replacing the light field indicated by the LED light source for the light field can be fundamentally eliminated in the application, and the error introduction link is further reduced; meanwhile, due to the fact that visible light indication is omitted, a technician can directly judge the visual shot through the interactive panel terminal, and therefore the efficiency, the convenience and the accuracy of the visual shot work are improved.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments of the present application, and for those skilled in the art, other drawings can be obtained according to the drawings without inventive exercise.
Fig. 1 is a schematic diagram of the present invention in a shooting scenario.
Figure 2 is a schematic view of the pitch and yaw of an X-ray source.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are some embodiments of the present application, but not all embodiments. The components of the embodiments of the present application, generally described and illustrated in the figures herein, can be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the present application, presented in the accompanying drawings, is not intended to limit the scope of the claimed application, but is merely representative of selected embodiments of the application. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
The preferred embodiment:
since the present application has not been applied in the field of X-ray photography, and belongs to a completely new invention creation, and since the present invention belongs to the field of computer electrical digital data processing, many terms of art will be involved, in order to show the present invention more thoroughly and clearly, and to understand the present invention better, first, the applicant makes the following description of the terms of art and the necessary principles involved in the present invention, and aims to provide a better reading environment for facilitating quick and accurate understanding of the present invention.
Before understanding the present invention, first, the working principle of X-ray is clarified, and in combination with the principle of the present invention shown in fig. 1, the X-ray source emits X-ray toward the detector for receiving X-ray and imaging to form a macroscopic X-ray film, and the detector is located at the imaging plane position as shown in the figure. The most popular image receiver is now the flat panel detector, and therefore the flat panel detector is also discussed herein as an image receiver. The flat panel detector can be seen as a square in space, which is always perpendicular to the horizontal plane in X-ray photography, and has X-ray sensing components thereon for receiving all the X-rays attenuated by the irradiated object and converting them into digital images. It is known that, because X-rays have radiation damage to living bodies, the radiation range of X-rays is controllable and strictly limited in X-ray detection, however, a device for limiting the X-ray range is called a beam limiter and is installed between an X-ray source and an imaging plane, only the position where the beam limiter allows the X-rays to pass through can the X-rays reach the imaging plane smoothly for imaging, and other positions where the X-rays are not allowed can be completely blocked to avoid unnecessary radiation damage, and the blocked material is generally an existing material such as a lead plate. The cross section of the beam limiter that allows X-rays to pass through is called the beam limiter opening.
Then, after understanding the basic principles of X-ray operation, the following definitions of the relevant terms are needed.
Definition 1: starting from the X-ray source, the X-ray starting right along the inclination angle direction of the X-ray source is taken as a central ray. If the direction of the X-ray source is exactly perpendicular to the vertically mounted detector, the X-ray in the horizontal direction and perpendicular to the vertical plane is the central ray.
Definition 2: the plane where the image receiver (i.e. flat panel detector) is located is the image receiving surface.
The central symmetry point of the image receiver is referred to as the center point of the image receiver.
In radiography, all intersection points of X-rays with the image receiving surface are on the image receiver, and the central ray intersects the image receiving surface at a midpoint of the image receiver.
Definition 3: the projection of the central ray on the surface of the irradiated object is the center of the field.
Definition 4: and a plane passing through the center of the field and parallel to the image receiving surface is a field plane.
Definition 5: and all X-rays pass through the beam limiter opening and form a region with the intersection point on the field plane, which is the field.
Definition 6: in any plane perpendicular to the central ray, a region can be chosen to equivalently describe the effect of the beam limiter on limiting the field, which is called the beam limiter opening.
During X-ray photography, a large number of X-rays are emitted from the X-ray source in a conical shape along the inclination angle direction of the X-ray source, and the X-rays just pass through points in the opening range of the beam limiter and are projected on a radiation field plane, and each single X-ray propagation path is the ray taking the X-ray source as an end point. It is to be noted that the plane of the beam limiter opening is the plane of the beam limiter opening, and in the initial position, the plane of the beam limiter opening is selected as the plane z = SED, which is the distance from the X-ray source to the plane of the beam limiter opening.
Let the X-ray source be the origin O of the world coordinate system.
Then, assuming a point p in the beam limiter opening plane, the intersection of the ray Op and the field plane is within the field if and only if the point p is within the beam limiter opening range.
That is, any point p of the beam limiter aperture plane, which is on some X-ray path if it is within the aperture range, has its intersection point with the field plane in the field according to definition 5; if the point p is not within the aperture, it is not necessarily on the optical path of any X-ray, and its intersection with the field plane is not in the field according to definition 5.
Definition 7: the distance of the X-ray source to the imaging plane is SID.
In a certain X-ray photograph, the midpoint of the image receiver is denoted as point I, and the distance OI between the X-ray source and the midpoint of the image receiver is:
OI=D=SID/(cosβ∙cosγ )
the image receiving surface is a plane where the SID is located.
Definition 8: the distance from the X-ray source to the radiation field plane is SOD.
If SOD is obtained, the ejection field plane is the plane where SOD is located.
Definition 9: the distance from the X-ray source to the optical center plane of the depth camera is SCD.
The SCD is determined after the depth camera installation is completed.
After the above definitions are clarified, the general X-ray photography virtual portal display method shown in fig. 1 and fig. 2 in conjunction with the description includes an interactive panel terminal device for displaying a photographed picture and a virtual portal, a depth camera for photographing a real-time picture, and interactive software and back-end processing software installed on the interactive panel terminal device, including the following steps:
STP100, establishing a world coordinate system with the X-ray source as an origin O through back-end processing software, and obtaining any point on the opening plane of the beam limiterv n Has world coordinates of
Pass through any pointv n Spatial rotation matrix occurring relative to originTCalculating to obtain the rotated point
World coordinate
Further, as a preferred embodiment of the present application, the method for obtaining the rotation matrix R in the present embodiment is not obtained by using the simplest existing method, that is, the trigonometric function relationship of the X-ray shooting scene is established by using data of the IMU accelerometer, which is the most common calculation method, but this method is not the optimal and most efficient method in practical application; in order to optimize the calculation accuracy of the rotation matrix R, the rotation matrix R in step STP100 in this implementation is obtained by the following method:
step STP110, passing any point
The spatial rotation generated by the X-ray is divided into a roll angle alpha, a pitch angle beta and a yaw angle gamma;
STP120, respectively calculating a rolling angle alpha, a pitch angle beta and a yaw angle gamma;
based on the X-ray photography principle, if the change of the roll angle alpha is not involved, the roll angle alpha is 0;
the pitch angle β is calculated as follows:
the formula for calculating the yaw angle γ is as follows:
step STP130, the rotation matrices for pitch and yaw are calculated separately as follows:
step STP140, compounding an X-ray rotation matrix R by a pitch angle beta and a yaw angle gamma as follows:
wherein rolling means winding
The rotary movement of the axis, pitch being about
The rotary movement of the shaft, yaw being about
The rotational movement of the shaft is such that,
、
、
respectively, three components of an IMU accelerometer mounted within the beam limiter.
Step STP200, decision point
If the current point is within the known opening range of the beam limiter, if the judgment result is yes, the step STP300 is carried out, if the judgment result is no, the step STP100 is repeated and the next point is judgedv n+1 Until all points within the range of the known beam limiter opening are traversedv n To whichn≥1；
Step STP300, according to the points
Obtaining rays
Intersection point with field planep v World coordinate ofp vw 
；
Specifically, the position of the portal plane is z =d+SCD，p vw 
The calculation method is as follows:
the principle of linear light propagation includes:
wherein
Is a similar ratio, so:
namely:
step STP400, by reading the external reference torque of the depth cameraK 2 Calculating to obtain an intersection pointp v Camera coordinates ofp vc And then acquiring internal reference torque of the depth camera through back-end processing softwareK 1 Calculating to obtain an intersection pointp v Image coordinates ofp vI Torque of external referenceK 2 The expression is as follows:
wherein,
is a matrix of rotations of the optical system,
is a translation vector;
internal reference torqueK 1 The expression is as follows:
wherein,
is the focal length of the camera and,
and
respectively the ratio of the actual object distance to the pixel in both directions,
is a translation vector from the origin of coordinates of the depth camera to the origin of coordinates of the image; further, in this embodiment, the camera coordinates are described in STP400 of this stepp vc The calculation of (c) is as follows:
p
vc
image coordinatesp vI The calculation method is as follows:
p vI 
。
step STP500, using interactive software to coordinate the image obtained in step STP400p vI Finally displayed on the interactive panel terminal,
wherein SED is the distance from the X-ray source to the opening plane of the beam limiter;
the opening plane of the beam limiter is a plane on the beam limiter, which is perpendicular to the X central ray;
the opening range of the beam limiter refers to the maximum cross-sectional range of the beam limiter for allowing X-rays to pass through;
the field plane is a plane passing through the field center and parallel to the image receiving surface;
the field is the area formed by the intersection of all the X-rays passing through the beam limiter opening and the field plane.
It should be noted that the beam limiter opening is usually rectangular, but the method of the present application includes, but is not limited to, only rectangular, and the application scenario includes, but is not limited to, rectangular cross section as long as it can be shown in a closed cross section shape, such as a circle, a triangle, a polygon, an ellipse, and other closed irregular cross sections, since the present application can convert any point world coordinate on a specific known plane (such as a plane where the SED is known) into a pixel coordinate on a plane where another specific distance (such as a plane where the SOD is located). Of course, in order to simplify the application, only the boundary line can be displayed on the interactive panel terminal, and the virtual field can be indicated by using the highlighted boundary line, so that the purpose of visually indicating the field to assist the technician in quickly positioning the shooting part can be achieved.
From the above conclusion, as another extended application of the present invention, when the field shape, size and SID of the illuminated object are known, the shape and size of the beam limiter opening can be obtained. The specific principle is as follows:
the world coordinate of any point p of the portal plane is
。
Since the size and shape of the field are known, it can be determined whether the point p is in the field. If the point p is in the field, according to conclusion 1, the intersection of Op and the plane of the beam limiter opening at that moment is within the beam limiter opening. At this time, the world coordinates of the intersection point v are easily solved. According to a point p in the opening plane of the beam limiter, if the intersection point of the ray Op and the field plane is in the field and only if the point p is in the range of the opening of the beam limiter, finding out the points in all the fields, and repeating the algorithm to obtain the size and the shape of the opening of the beam limiter.
In practical implementation, since the field can be approximated by a quadrilateral and the beam limiter opening is rectangular, the above process can be simplified as follows: the world coordinates of the four vertexes of the field quadrilateral are selected, so that the world coordinates of the four vertexes of the opening rectangle of the beam limiter can be determined, and the size and the shape of the opening of the beam limiter can be determined.
It follows that the invention can be applied to two opposite scenarios, namely:
the first situation is as follows: when the shape, the size and the SID of the opening of the beam limiter are known, the shape, the size and the center of the radiation field can be obtained;
case two: knowing the desired field shape, size and SID of the illuminated object, the desired shape and size of the beam limiter opening can be obtained. Although the two scenarios are different, the actual calculation is not different, and in practical application, the technical problems to be solved by the two scenarios can be solved in the present invention.
The concrete brief description is as follows:
suppose any point on the portal plane
Has world coordinates of
。
Since the size and shape of the field are known, the point can be identified
Whether it is in the field. If point
In the field of emission, the
The intersection point with the plane of the beam limiter opening at this moment is within the beam limiter opening. The intersection is solved below
The world coordinates of (2) can be solved first
Corresponding point when the X-ray is not rotated
，
There is also a corresponding point when the X-ray is not rotated
：
Intersection point
Corresponding point when the X-ray is not rotated
Can be regarded as
With the field plane when the X-ray is not rotated
The intersection of (a), in this case, by the nature of a similar triangle, is:
wherein
Is a similar ratio, so:
namely:
the size and shape of the beam limiter opening can be obtained by repeating the above algorithm after finding out the points in all fields of view.
In practical implementation, since the field can be approximated by a quadrilateral and the beam limiter opening is rectangular, the above process can be simplified as follows: the world coordinates of the four vertexes of the field quadrilateral are selected, so that the world coordinates of the four vertexes of the opening rectangle of the beam limiter can be determined, and the size and the shape of the opening of the beam limiter can be determined.
The above description is only a preferred embodiment of the present application and is not intended to limit the present application, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present application shall be included in the protection scope of the present application.
Claims (3)
1. The virtual radiation field display method for common X-ray photography comprises interactive panel terminal equipment for displaying a photographed picture and a virtual radiation field, a depth camera for photographing a real-time picture, interactive software and back-end processing software which are installed on the interactive panel terminal equipment, and is characterized in that: the method comprises the following steps:
STP100, establishing a world coordinate system with the X-ray source as an origin O through back-end processing software, and acquiring any point on the opening plane of the beam limiterv n Has world coordinates of
Pass through any pointv n Spatial rotation matrix occurring relative to originTCalculating to obtain the rotated point
World coordinate
STP200, judging point
If the current point is within the known opening range of the beam limiter, the step STP300 is carried out, if the current point is not within the known opening range of the beam limiter, the step STP100 is repeated, and the next point is judgedv n+1 Until all points within the range of the known beam limiter opening are traversedv n To therein, whereinn≥1；
Step STP300, according to the points
Obtaining rays
Point of intersection with field planep v World coordinate ofp vw 
；
Step STP400, by reading the external reference torque of the depth cameraK 2 Calculating to obtain an intersection pointp v Camera coordinates ofp vc And then acquiring internal reference torque of the depth camera through back-end processing softwareK 1 Calculating to obtain an intersection pointp v Image coordinates ofp vI Torque of external referenceK 2 The expression is as follows:
wherein,Ris a matrix of rotations of the optical system,
is a translation vector;
internal reference torqueK 1 The expression is as follows:
wherein,
is the focal length of the camera and,
and
respectively the ratio of the actual object distance to the pixel in both directions,
is a translation vector from the origin of coordinates of the depth camera to the origin of coordinates of the image;
step STP500, using interactive software to coordinate the image obtained in step STP400p vI Finally displayed on the interactive panel terminal,
wherein SED is the distance from the X-ray source to the opening plane of the beam limiter;
the opening plane of the beam limiter is a plane on the beam limiter, which is perpendicular to the X central ray;
the opening range of the beam limiter refers to the maximum cross-sectional range of the beam limiter for allowing X-rays to pass through;
the field plane is a plane passing through the field center and parallel to the image receiving surface;
the field is the area formed by the intersection of all the X-rays passing through the beam limiter opening and the field plane.
2. The method of displaying a virtual portal for general radiography according to claim 1, wherein: the rotation matrix R in step STP100 is obtained by the following method:
step STP110, passing any point
The spatial rotation generated by the X-ray is divided into a roll angle alpha, a pitch angle beta and a yaw angle gamma;
STP120, respectively calculating a rolling angle alpha, a pitch angle beta and a yaw angle gamma;
based on the X-ray photography principle, if the change of the roll angle alpha is not involved, the roll angle alpha is 0;
the pitch angle β is calculated as follows:
the formula for calculating the yaw angle γ is as follows:
step STP130, the rotation matrices for pitch and yaw are calculated separately as follows:
step STP140, compounding an X-ray rotation matrix R by a pitch angle beta and a yaw angle gamma as follows:
wherein roll refers to rotational motion about the z-axis, pitch refers to rotational motion about the x-axis, yaw refers to rotational motion about the y-axis,
、
、
respectively, three components of an IMU accelerometer mounted within the beam limiter.
3. The method of displaying a virtual portal for general radiography according to claim 1, wherein: the camera coordinates in step STP400p vc The calculation of (c) is as follows:
p
vc
image coordinatesp vI The calculation method is as follows:
p vI 
。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210962947.5A CN115024740B (en) | 2022-08-11 | 2022-08-11 | Virtual radiation field display method for common X-ray photography |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210962947.5A CN115024740B (en) | 2022-08-11 | 2022-08-11 | Virtual radiation field display method for common X-ray photography |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115024740A true CN115024740A (en) | 2022-09-09 |
| CN115024740B CN115024740B (en) | 2022-10-25 |
Family
ID=83130505
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210962947.5A Active CN115024740B (en) | 2022-08-11 | 2022-08-11 | Virtual radiation field display method for common X-ray photography |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115024740B (en) |
Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6385287B1 (en) * | 2000-12-06 | 2002-05-07 | Siemens Medical Solutions Usa, Inc. | Method and system for providing virtual grid for portal imaging in a radiotherapy system |
| CN101219054A (en) * | 2007-12-26 | 2008-07-16 | 东软飞利浦医疗设备系统有限责任公司 | Method for displaying image variation in X-ray irradiation area with synchronous simulation |
| CN104605872A (en) * | 2014-12-31 | 2015-05-13 | 沈阳东软医疗系统有限公司 | Radiation field control method and device for medical X-ray system |
| US20160296195A1 (en) * | 2015-04-09 | 2016-10-13 | Kabushiki Kaisha Toshiba | X-ray diagnostic apparatus, and method of adjusting irradiation range of x-ray diagnostic apparatus |
| CN107242880A (en) * | 2017-06-06 | 2017-10-13 | 深圳市创谷科技发展有限公司 | X-ray beam limiter and X-ray photography system |
| EP3311749A1 (en) * | 2015-06-22 | 2018-04-25 | FUJIFILM Corporation | Radiation-irradiating device, control method for radiation-irradiating device, and program |
| CN109999371A (en) * | 2017-12-21 | 2019-07-12 | 瓦里安医疗系统国际股份公司 | For the virtual beam weights optimization imaging of patient's setting in radiotherapy |
| WO2019228359A1 (en) * | 2018-05-28 | 2019-12-05 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for taking x-ray images |
| CN110740687A (en) * | 2017-06-15 | 2020-01-31 | 皇家飞利浦有限公司 | X-ray imaging apparatus |
| CN110811654A (en) * | 2019-11-12 | 2020-02-21 | 飞瑞医疗器械(嘉兴)有限公司 | X-ray exposure control system and control method thereof |
| CN110859641A (en) * | 2019-11-14 | 2020-03-06 | 深圳市博视科技有限公司 | Light radiation field visual detection method and light radiation field visual detection equipment |
| CN111540022A (en) * | 2020-05-14 | 2020-08-14 | 深圳市艾为智能有限公司 | Image uniformization method based on virtual camera |
| CN111870268A (en) * | 2020-07-30 | 2020-11-03 | 上海联影医疗科技有限公司 | Method and system for determining target position information of beam limiting device |
| CN112102413A (en) * | 2020-07-22 | 2020-12-18 | 西安交通大学 | Virtual lane line-based automatic calibration method for vehicle-mounted camera |
| CN114343689A (en) * | 2022-03-17 | 2022-04-15 | 晓智未来(成都)科技有限公司 | Method for measuring opening area of beam limiter based on photogrammetry and application |
-
2022
- 2022-08-11 CN CN202210962947.5A patent/CN115024740B/en active Active
Patent Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6385287B1 (en) * | 2000-12-06 | 2002-05-07 | Siemens Medical Solutions Usa, Inc. | Method and system for providing virtual grid for portal imaging in a radiotherapy system |
| CN101219054A (en) * | 2007-12-26 | 2008-07-16 | 东软飞利浦医疗设备系统有限责任公司 | Method for displaying image variation in X-ray irradiation area with synchronous simulation |
| CN104605872A (en) * | 2014-12-31 | 2015-05-13 | 沈阳东软医疗系统有限公司 | Radiation field control method and device for medical X-ray system |
| US20160296195A1 (en) * | 2015-04-09 | 2016-10-13 | Kabushiki Kaisha Toshiba | X-ray diagnostic apparatus, and method of adjusting irradiation range of x-ray diagnostic apparatus |
| EP3311749A1 (en) * | 2015-06-22 | 2018-04-25 | FUJIFILM Corporation | Radiation-irradiating device, control method for radiation-irradiating device, and program |
| CN107242880A (en) * | 2017-06-06 | 2017-10-13 | 深圳市创谷科技发展有限公司 | X-ray beam limiter and X-ray photography system |
| CN110740687A (en) * | 2017-06-15 | 2020-01-31 | 皇家飞利浦有限公司 | X-ray imaging apparatus |
| CN109999371A (en) * | 2017-12-21 | 2019-07-12 | 瓦里安医疗系统国际股份公司 | For the virtual beam weights optimization imaging of patient's setting in radiotherapy |
| WO2019228359A1 (en) * | 2018-05-28 | 2019-12-05 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for taking x-ray images |
| CN110811654A (en) * | 2019-11-12 | 2020-02-21 | 飞瑞医疗器械(嘉兴)有限公司 | X-ray exposure control system and control method thereof |
| CN110859641A (en) * | 2019-11-14 | 2020-03-06 | 深圳市博视科技有限公司 | Light radiation field visual detection method and light radiation field visual detection equipment |
| CN111540022A (en) * | 2020-05-14 | 2020-08-14 | 深圳市艾为智能有限公司 | Image uniformization method based on virtual camera |
| CN112102413A (en) * | 2020-07-22 | 2020-12-18 | 西安交通大学 | Virtual lane line-based automatic calibration method for vehicle-mounted camera |
| CN111870268A (en) * | 2020-07-30 | 2020-11-03 | 上海联影医疗科技有限公司 | Method and system for determining target position information of beam limiting device |
| CN114343689A (en) * | 2022-03-17 | 2022-04-15 | 晓智未来(成都)科技有限公司 | Method for measuring opening area of beam limiter based on photogrammetry and application |
Non-Patent Citations (3)
| Title |
|---|
| CHRISTOPHER HAHNE等: "Baseline of virtual cameras acquired by a standard plenoptic camera setup", 《3DTV-CON 2014》 * |
| 任冲等: "一种扇面激光X线束中心线指示装置", 《医疗卫生装备》 * |
| 谈友恒: "基于非晶硅平板探测器的电子射野影像系统的研究", 《中国优秀硕士学位论文全文数据库医药卫生科技辑》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115024740B (en) | 2022-10-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chakraborty | Image intensifier distortion correction | |
| US7085430B2 (en) | Correcting geometric distortion in a digitally captured image | |
| WO2022193558A1 (en) | Projector correction method and system, and storage medium and electronic device | |
| US7697126B2 (en) | Three dimensional spatial imaging system and method | |
| US9795021B2 (en) | X-ray equipment and alignment method of X-ray equipment | |
| US20150116691A1 (en) | Indoor surveying apparatus and method | |
| US20150216498A1 (en) | Geometric Characterization and Calibration of a Cone-Beam Computer Tomography Apparatus | |
| CN105486702B (en) | A kind of target defect detecting system based on X ray | |
| JP2011019707A (en) | X-ray imaging device, control method for x-ray imaging device, and program | |
| EP2702449B1 (en) | System and method for correction of geometric distortion of multi-camera flat panel x-ray detectors | |
| KR20160089762A (en) | An radiography apparatus, a method for controlling the radiography apparatus | |
| CN116019474A (en) | Multi-source imaging device and method | |
| CN117717367A (en) | Auxiliary positioning system and method for standing position computer tomography | |
| JP2002310943A (en) | Computed tomograph | |
| CN115024740B (en) | Virtual radiation field display method for common X-ray photography | |
| US10610177B2 (en) | Method for imaging by means of an X-ray device and X-ray device | |
| JP2019196991A (en) | Nuclear medicine diagnosis device and position correction method | |
| JP5684670B2 (en) | Gamma camera system | |
| CN116678322A (en) | Crack width measurement method and system considering parallel laser beam inclination angle | |
| US11948290B2 (en) | Individual channel characterization of collimator | |
| JP2021173716A (en) | Inspection device and inspection method | |
| WO2020209313A1 (en) | Image processing device, and image processing method | |
| CN108052213A (en) | Position indication method, device and system | |
| KR20170107893A (en) | Computed tomography apparatus | |
| JPWO2020209312A1 (en) | Inspection equipment and inspection method |
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 |