CN115588011B - X-ray machine projection model parameter calibration method - Google Patents
X-ray machine projection model parameter calibration method Download PDFInfo
- Publication number
- CN115588011B CN115588011B CN202211587624.9A CN202211587624A CN115588011B CN 115588011 B CN115588011 B CN 115588011B CN 202211587624 A CN202211587624 A CN 202211587624A CN 115588011 B CN115588011 B CN 115588011B
- Authority
- CN
- China
- Prior art keywords
- image
- imaging
- points
- plane
- ray machine
- 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/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
- A61B6/584—Calibration using calibration phantoms determining position of components of the apparatus or device using images of the phantom
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- 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
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30204—Marker
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Biophysics (AREA)
- Quality & Reliability (AREA)
- High Energy & Nuclear Physics (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The invention discloses a method for calibrating parameters of a projection model of an X-ray machine, which comprises the following steps: placing a calibration phantom in an X-ray machine for perspective, wherein at least 3 asymmetric marking points are arranged on the calibration phantom; calculating to obtain a homography between the calibration phantom and the image plane according to the coordinates of the imaging points of the marking points on the obtained perspective image and the space coordinates of the marking points, and determining corresponding rays by combining a perspective imaging principle to further obtain the position of the X-ray source; and calculating to obtain the homography between the imaging plane and the image plane according to the homography between the calibration phantom and the image plane and the projection model of the X-ray machine, and accordingly obtaining the pose of the imaging plane, so as to calculate the position deviation of the X-ray source between the vertical point of the imaging plane and the central point of the perspective image, and further obtain the deviation between the imaging plane and the image plane. The invention can calibrate the projection model parameters of the current X-ray machine in real time.
Description
Technical Field
The invention relates to the technical field of X-ray machines, in particular to a method for calibrating parameters of a projection model of an X-ray machine.
Background
At present, with the wide application of medical surgical robots, a C-arm X-ray machine is widely used in image imaging in interventional radiology department orthopedic surgery, the C-arm X-ray machine is used as a medium for imaging and an actual patient, the imaging model is used for establishing mapping between the actual patient and the formed image, based on the mapping relation, a surgical robot system can process and analyze the patient image, so that the processing and analysis of the actual patient information can be completed, and the C-arm X-ray machine is a very key ring in a medical surgical robot, particularly an auxiliary surgical robot system.
With the long-term use of the X-ray machine, the rigidity of the structural part and the position of the internal key components are loosened, so that the parameters of the projection model are changed, and the deviation exists between the imaging plane and the image plane, therefore, the accurate calibration of the parameters of the projection model of the X-ray machine to obtain the deviation between the imaging plane and the image plane is a very critical work.
Disclosure of Invention
The purpose of the invention is as follows: aiming at the defects, the invention provides an efficient X-ray machine projection model parameter calibration method, which can calibrate the projection model parameters of the current X-ray machine in real time to obtain the deviation between an imaging plane and an image plane.
The technical scheme is as follows:
an X-ray machine projection model parameter calibration method comprises the following steps:
placing a calibration phantom in an X-ray machine for perspective, wherein at least 3 asymmetric marking points are arranged on the calibration phantom;
calculating to obtain a homography between the calibration phantom and the image plane according to the coordinates of the imaging points of the marking points on the obtained perspective image and the space coordinates of the marking points, and determining corresponding rays by combining a perspective imaging principle to further obtain the position of the X-ray source;
and calculating to obtain the homography between the imaging plane and the image plane according to the homography between the calibration phantom and the image plane and the projection model of the X-ray machine, and accordingly obtaining the pose of the imaging plane, so as to calculate the position deviation of the X-ray source between the vertical point of the imaging plane and the central point of the perspective image, and further obtain the deviation between the imaging plane and the image plane.
The homography between the imaging plane and the image plane obtained by calculation is specifically as follows:
and calculating the coordinates of the imaging points in the actual space according to the direction vectors of the rays between the marking points on the calibration phantom and the X-ray source and the coordinates of the imaging points of the marking points on the obtained perspective image in the image coordinate system, and further calculating according to the projection model of the X-ray machine to obtain the homography between the imaging plane and the image plane.
The obtained imaging plane pose specifically comprises the following steps:
and obtaining coordinates of four corner points on the perspective image in an image coordinate system according to the resolution of the perspective image, calculating the coordinates of the four corner points on the perspective image in the imaging plane coordinate system according to the homography between the imaging plane and the image plane, and calculating the pose of the imaging plane according to the coordinates.
And calculating the direction vectors of the x axis and the y axis of the coordinate system corresponding to the imaging plane according to the coordinates of the four corner points on the perspective image in the coordinate system of the imaging plane.
At least 3 coplanar mark points are embedded on two end faces of the calibration die body, which are used for penetrating through X rays, and the mark points on the two end faces are asymmetric.
The marking points are steel balls.
The method comprises the steps of respectively obtaining homography matrixes between two end surfaces and an image plane by identifying coordinates of imaging points of marking points on the two end surfaces on a perspective image in an image coordinate system and space coordinates of each marking point according to a perspective imaging principle, further calculating rays penetrating through each marking point, and obtaining an intersection point of each ray as an X-ray source.
Has the advantages that: the invention can calibrate the parameters of the X-ray machine projection model, accurately calibrate the X-ray machine projection model, prevent the equipment precision from being reduced after long-term use by calculating the deviation between the imaging plane and the image plane, and has important application scenes in the surgical robot system based on image analysis.
Drawings
FIG. 1 is a diagram of a projection model of an X-ray machine;
FIG. 2 is a schematic view of a calibration phantom with embedded markers;
FIG. 3 is a schematic view of the use of a calibration phantom;
FIG. 4 is a schematic view of a projection model with a calibration phantom positioned for fluoroscopy within the range of the X-ray machine;
FIG. 5 is a schematic diagram of a coordinate system of a fluoroscopic image obtained by fluoroscopy of an X-ray machine and a mark point therein;
FIG. 6 is a schematic view of projection model imaging;
FIG. 7 is a schematic diagram of a homography;
FIG. 8 is a schematic view of the projection imaging of the steel ball on the end Plane 1.
Detailed Description
The invention is further elucidated with reference to the drawings and the embodiments.
A projection model of an X-ray machine by which projection is made is shown in FIG. 1The phantom samples an object placed in a cone-shaped region formed by an X-ray source and a flat-bed receptor as a perspective image on the flat-bed receptor, wherein,P x-ray the X-ray source is a projection model and corresponds to the center of an X-ray emission tube of an actual X-ray machine;andrespectively setting direction vectors of an x axis and a y axis of an imaging plane coordinate system established on an imaging plane under a reference coordinate system, wherein the reference coordinate system is a robot system coordinate system;P imgCenter is the projection point of the X-ray source on the imaging plane, namely the central point of the obtained perspective image in an ideal state, wherein the pointP imgCenter And、the positions of the imaging plane of the X-ray projection model under the reference coordinate system are represented together, but with the long-term use of the X-ray machine, the rigidity of a structural part and the positions of internal key components are loosened, so that the parameters of the projection model are changed, namely, the deviation exists between the imaging plane and the image plane obtained through final projection, and the calibration of the parameters of the X-ray machine projection model is to calibrate the deviation.
The invention relates to a method for calibrating parameters of a projection model of an X-ray machine, which comprises the following steps:
(1) Placing the calibration phantom in the range of fluoroscopy of the X-ray machine for fluoroscopy, as shown in FIG. 3;
the structure of the calibration phantom is shown in fig. 2, the calibration phantom is in a cube shape and is prepared by adopting an X-ray transmitting material, at least 3 coplanar mark points are embedded on a Plane1 and a Plane2 which are used for transmitting X-rays, and the mark points on the two end surfaces are asymmetric; then each mark is known according to design parametersRelative position relation between points, and establishing a reference coordinate system C according to the relative position relation 0 Then, the transformation relation between the reference coordinate system and the reference coordinate system can be obtained through measuring or calibrating the installation position of the mold body, so that the projection model is associated with the robot system, and therefore the position and attitude information of the parameters in the projection model in the reference coordinate system is obtained through parameter calibration of the projection model;
in the invention, the number of the marking points on the two end surfaces Plane1 and Plane2 is 4, specifically steel balls, correspondingly, the serial numbers of the steel balls on the end surface Plane1 close to the ray source are sequentiallyp1~p4, the steel balls on the end surface Plane2 far away from the ray source are numbered sequentiallyp5~p8。
In the invention, at least 3 asymmetrical marking points can be arranged on the calibration phantom, but for increasing the accuracy, the invention preferably embeds 4 coplanar marking points on both end faces Plane1 and Plane2 of the calibration phantom for transmitting X-rays.
(2) Respectively calculating coordinates of each marking point on a perspective image obtained by perspective and coordinates of each marking point under a reference coordinate system to obtain homography matrixes between two end faces of a calibration mold body and an image plane, and further determining corresponding rays by a perspective imaging principle to obtain the position of an X-ray source;
the calibration phantom is placed in a projection model for fluoroscopy in the fluoroscopy range of the X-ray machine as shown in fig. 4, and a perspective image obtained through fluoroscopy is shown in fig. 5, wherein C 1 The image coordinate system of the output perspective image is a coordinate system established by taking the upper left corner point of the perspective image as an origin,p1′~p8' are imaging points of the marking points on the calibration phantom on the output perspective image respectively;
FIG. 6 is a schematic view of a projection model, as shown in FIG. 6, in which a single imaging point on a fluoroscopic image is imaged with a pointp1' imaging procedure as an example: x-ray sourceP x-ray With alignment of the upper face Plane1 of the phantomp1 the intersection of the component rays L1, L1 with the end Plane2 isp1', the intersection point with the image plane is the imaging pointp1', other marker imaging procedures can be obtained in the same way, i.e. coplanar markers on the end-face Plane1p1~p4 imaging onto image planep1′~p4', coplanar markers on the end Plane2p5~p8 imaging on the image planep5′~p8', a homography matrix between the two end faces and the image plane can be obtainedhm1 andhm2, wherein,hm1 andhm2 are all 3 x 3 matrices;
by identifying imaging points on the fluoroscopic imagep1′~p8' in the coordinate of image coordinate system and the coordinate of each mark point in reference coordinate system, the homography matrix between two end faces and image plane can be obtained respectively by perspective imaging principlehm1 andhm2:
further, the passing through mark point can be obtainedp1~p8, the intersection point of each ray is the X-ray source, and the X-ray source can be obtained only through the intersection point of any two rays;
in the invention, at least 3 asymmetric marking points can be arranged on the calibration phantom, so that a single mapping matrix between the calibration phantom and an image plane can be obtained by a perspective imaging principle according to the coordinates of the corresponding marking points on a perspective image and the coordinates of each marking point under a reference coordinate system, rays passing through each marking point can be further obtained, and the intersection point of each ray is an X-ray source.
(3) Solving the deviation between the imaging plane and the image plane;
(31) The coordinates of the imaging points under the reference coordinate system can be solved based on the projection model of the end surface Plane1 or Plane2, and the homography relation between the imaging Plane and the image Plane is obtained through calculation according to the coordinates;
the invention takes the projection model of the end Plane1 as an example, the marking points on the end Plane1p1~p4 coordinate is alreadyKnown as an X-ray sourceP x-ray And each mark point on the end surface Plane1p1~p4 form rays L1-L4, respectively, and the rays L1-L4 intersect the image plane atp1′~p4', then the direction vectors of the rays L1~ L4 can be obtained as follows:
imaging pointp1′~pThe coordinates of 4' in the reference coordinate system can be expressed as:
imaging pointp1′~p4' are known in the image coordinate system, and are respectively set asp1′(x1,y1)、p2′(x2,y2)、p3′(x3,y3)、p4′(x4,y4) Size of picture element on imaging planekIs an inherent property of an X-ray machine and is known; the imaging point is formed in the real spacep1′~pThe physical size coordinates of 4' in the image coordinate system can be expressed as follows:
from which the coefficients can be solveda、b、cAnddthereby obtaining an imaging pointp1′~p4' in a reference coordinate system;
then according to the imaging pointp1′~p4' in the image coordinate system, the transformation relation between the reference coordinate system and the image coordinate system can be obtained by calculation, so that the homography relation between the imaging plane and the image plane can be obtained according to the X-ray machine projection model, and the homography matrix is adoptedhm3 represents;
similarly, the homography between the imaging plane and the image plane can also be obtained by only adopting at least 3 asymmetric marking points arranged on the calibration phantom and the corresponding imaging points;
(32) Solving to obtain the pose of an imaging plane, and calculating to obtain the position transformation relation between the projection point of the X-ray source on the imaging plane and the image center point;
the resolution of the fluoroscopic image obtained by the X-ray machine is (u,v) Then the coordinates of four corner points on the perspective image under the image coordinate system are D1 (0, 0) respectively D2 (a)u,0)、D3(u,-v)、D4(0,-v) Then, the coordinates D1' = of the four corner points on the fluoroscopic image in the imaging plane coordinate system can be calculated according to the obtained homography between the imaging plane and the image planehm3*D1、D2′=hm3*D2、D3′=hm3*D3、D4′=hm3 x D4, then the coordinate system of the imaging plane can be obtainedx、yAxial direction vector:
wherein norm () represents a normalized vector function;
the position and the attitude of the imaging plane can be calculated through the coordinates of the four corner points in the imaging plane coordinate system, and then the X-ray source can be obtainedP x-ray Drop foot point on imaging planeP x According to the pointsP x Obtaining the position transformation relation between the imaging plane and the image plane by the position deviation between the imaging plane and the central point of the image;
(33) And (4) calculating the deviation between the imaging plane and the image plane according to the homography relation between the imaging plane and the image plane obtained in the step (31) and the position conversion relation between the imaging plane and the image plane obtained in the step (32).
In the invention, the perspective image is collected by the X-ray machine, and the perspective image is processed and analyzed by different robot systems according to the perspective image, so as to realize different purposes.
The invention calculates the homography between the imaging plane and the image plane by the coordinate of the marking point in the calibration die body on the perspective image and the actual space coordinate, and calculates the pose of the imaging plane on the basis of the homography so as to obtain the X-ray sourceP x-ray Drop foot point on imaging planeP x According to the pointsP x And the position deviation between the X-ray machine projection model parameter calibration device and the image center point obtains the position conversion relation between the imaging plane and the image plane, and finally the deviation between the imaging plane and the image plane is obtained through calculation so as to calibrate the X-ray machine projection model parameter and complete accurate calibration of the X-ray machine projection model parameter.
Although the preferred embodiments of the present invention have been described in detail, the present invention is not limited to the details of the foregoing embodiments, and various equivalent changes (such as number, shape, position, etc.) may be made to the technical solution of the present invention within the technical spirit of the present invention, and these equivalent changes are all within the protection scope of the present invention.
Claims (6)
1. A method for calibrating parameters of a projection model of an X-ray machine is characterized by comprising the following steps: the method comprises the following steps:
placing a calibration phantom in an X-ray machine for fluoroscopy, wherein at least 3 asymmetric marking points are arranged on the calibration phantom;
calculating to obtain a homography between a calibration phantom and an image plane according to the coordinates of imaging points of the marking points on the obtained perspective image and the space coordinates of the marking points, and determining corresponding rays by combining a perspective imaging principle to further obtain the position of an X-ray source;
calculating to obtain the homography between the imaging plane and the image plane according to the homography between the calibration phantom and the image plane and the projection model of the X-ray machine, and accordingly obtaining the pose of the imaging plane, so as to calculate to obtain the position deviation of the X-ray source between the vertical point of the imaging plane and the central point of the perspective image, and further obtain the deviation between the imaging plane and the image plane;
the homography between the imaging plane and the image plane obtained by calculation is specifically as follows:
and calculating the coordinates of the imaging points in the actual space according to the direction vectors of the rays between the marking points on the calibration phantom and the X-ray source and the coordinates of the imaging points of the marking points on the obtained perspective image in the image coordinate system, and further calculating according to the projection model of the X-ray machine to obtain the homography between the imaging plane and the image plane.
2. The method for calibrating the parameters of the projection model of the X-ray machine according to claim 1, wherein: the method for obtaining the pose of the imaging plane specifically comprises the following steps:
and obtaining coordinates of four corner points on the perspective image in an image coordinate system according to the resolution of the perspective image, calculating the coordinates of the four corner points on the perspective image in the imaging plane coordinate system according to the homography between the imaging plane and the image plane, and calculating the pose of the imaging plane according to the coordinates.
3. The method for calibrating the parameters of the projection model of the X-ray machine according to claim 2, characterized in that: and calculating the direction vectors of the x axis and the y axis of the coordinate system corresponding to the imaging plane according to the coordinates of the four corner points on the perspective image in the coordinate system of the imaging plane.
4. The method for calibrating the parameters of the projection model of the X-ray machine according to any one of claims 1 to 3, wherein the method comprises the following steps: at least 3 coplanar mark points are embedded on two end faces of the calibration die body, which are used for penetrating through X rays, and the mark points on the two end faces are asymmetric.
5. The method for calibrating the parameters of the projection model of the X-ray machine according to claim 4, wherein: the marking points are steel balls.
6. The method for calibrating the parameters of the projection model of the X-ray machine according to claim 4, wherein: the method comprises the steps of identifying coordinates of imaging points of marking points on two end surfaces in an image coordinate system of a perspective image and space coordinates of the marking points, respectively obtaining homography matrixes between the two end surfaces and an image plane according to a perspective imaging principle, further calculating rays passing through the marking points, and obtaining intersection points of the rays as an X-ray source.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211587624.9A CN115588011B (en) | 2022-12-12 | 2022-12-12 | X-ray machine projection model parameter calibration method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211587624.9A CN115588011B (en) | 2022-12-12 | 2022-12-12 | X-ray machine projection model parameter calibration method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115588011A CN115588011A (en) | 2023-01-10 |
CN115588011B true CN115588011B (en) | 2023-03-10 |
Family
ID=84783088
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211587624.9A Active CN115588011B (en) | 2022-12-12 | 2022-12-12 | X-ray machine projection model parameter calibration method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115588011B (en) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101467887B (en) * | 2007-12-29 | 2011-04-27 | 复旦大学 | X ray perspective view calibration method in operation navigation system |
JP6131606B2 (en) * | 2013-01-21 | 2017-05-24 | 株式会社島津製作所 | Radiation imaging apparatus and image processing method therefor |
CN113538259A (en) * | 2021-06-17 | 2021-10-22 | 中科超精(南京)科技有限公司 | Real-time geometric correction method for perspective imaging device |
CN113808211A (en) * | 2021-08-09 | 2021-12-17 | 中科超精(南京)科技有限公司 | Geometric correction die body for image guidance system and manufacturing method |
-
2022
- 2022-12-12 CN CN202211587624.9A patent/CN115588011B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN115588011A (en) | 2023-01-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3254621B1 (en) | 3d image special calibrator, surgical localizing system and method | |
US9560291B2 (en) | Automatic improvement of tracking data for intraoperative C-arm images in image guided surgery | |
CN112006779A (en) | Precision detection method for surgical navigation system | |
CN107049488A (en) | A kind of monoplane location of operation method and model | |
CN112168357B (en) | System and method for constructing spatial positioning model of C-arm machine | |
CN115568945A (en) | Position calibration method for surgical navigation reference mark point | |
JPH11160021A (en) | Wide area three-dimensional position measuring method and equipment | |
US20220054199A1 (en) | Robotic surgery systems and surgical guidance methods thereof | |
CN112190328A (en) | Holographic perspective positioning system and positioning method | |
WO2024027857A1 (en) | Method and system for registration of surgical robot coordinate system with ct scanner coordinate system | |
CN110645928A (en) | Space coordinate positioning method of three-dimensional scanner | |
CN115588011B (en) | X-ray machine projection model parameter calibration method | |
JP2006223868A (en) | Method and apparatus for determining geometrical parameters of imaging | |
US10991113B2 (en) | Gyroscope-based system and method for assisting in tracking heat source on mechanical arm | |
JP4187830B2 (en) | Medical image synthesizer | |
WO2022166182A1 (en) | Two-dimensional image guided intramedullary needle distal locking robot system and locking method therefor | |
CN115836875A (en) | Correction method and system | |
CN115690192A (en) | Imaging geometric relation online calibration method and related product | |
JP4409004B2 (en) | Surgical calibration system | |
CN117462267B (en) | Aiming method of robot end effector under perspective guidance | |
CN115908121B (en) | Endoscope registration method, device and calibration system | |
CN114176779B (en) | Surgical robot navigation positioning method and device | |
Cheryauka et al. | 3-D geometry calibration and markerless electromagnetic tracking with a mobile C-arm | |
JPH1080882A (en) | Coordinate transformation parameter measuring device for robot | |
CN116959675A (en) | External fixation three-dimensional reconstruction and electronic prescription generation method based on scribing |
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 |