WO2018133090A1 - Parameter calibration method, device, and system for x-ray machine - Google Patents
Parameter calibration method, device, and system for x-ray machine Download PDFInfo
- Publication number
- WO2018133090A1 WO2018133090A1 PCT/CN2017/072177 CN2017072177W WO2018133090A1 WO 2018133090 A1 WO2018133090 A1 WO 2018133090A1 CN 2017072177 W CN2017072177 W CN 2017072177W WO 2018133090 A1 WO2018133090 A1 WO 2018133090A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- rotating arm
- projection
- rotation
- standard phantom
- target
- 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.)
- Ceased
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
Definitions
- the present disclosure relates to signal processing techniques, such as a method, apparatus, and system for parameter calibration of an X-ray machine.
- parameter calibration for the difference between system parameters and design goals in static state can be used.
- the present disclosure provides a method, device and system for parameter calibration of an X-ray machine, which can calibrate errors generated in the movement of the X-ray machine.
- a parameter calibration method for an X-ray machine may include:
- the error parameter is calibrated according to the projection result of the standard phantom to obtain a calibration result of the error parameter
- a dynamic correction matrix for correcting the actual projection result of the X-ray machine is generated based on the mechanical structural motion state of the rotating arm at each set rotational position and the calibration result of the error parameter.
- a parameter calibration device for an X-ray machine which can be used to perform parameter calibration of the X-ray machine described in this embodiment
- the device can include:
- the projection result acquisition module is configured to control the X-ray machine rotating arm to perform a rotation scan around the standard phantom, and collect the projection results of the standard phantom at a plurality of set rotation positions;
- a motion state recording module configured to record a mechanical structural motion state of the rotating arm at each set rotational position
- the error parameter calibration module is configured to, at each set rotation position, calibrate the error parameter according to a projection result of the standard phantom to obtain a calibration result of the error parameter;
- a dynamic correction matrix generating module configured to generate a correction for correcting an actual projection result of the X-ray machine according to a mechanical structural motion state of the rotating arm at each set rotational position and a calibration result of the error parameter Dynamic correction matrix.
- a parameter calibration system for an X-ray machine may include:
- a rotating arm an X-ray source disposed opposite the rotating arm, and a detector, a processor, an image acquisition device, a motion state acquisition component, and a standard phantom;
- the rotating arm is configured to drive the X-ray light source and the detector to rotate a plurality of set rotation angles around the standard phantom according to a control signal sent by the processor;
- the X-ray source is configured to send an X-ray scan signal to the standard phantom
- the detector is configured to collect projection data of the X-ray scan signal on the standard phantom
- the image acquisition device is configured to generate a projection result of a standard phantom according to the projection data collected by the detector;
- the motion state collecting component is configured to acquire a mechanical structural motion state of the rotating arm
- the processor is configured to control the rotating arm to perform a rotational scan around the standard phantom to obtain a projection result at a plurality of set rotational positions; at each set rotational position, record the mechanical structural motion of the rotating arm a state; at each set rotational position, the error parameter is calibrated according to a projection result of the standard phantom to obtain a calibration result of the error parameter; according to the mechanism of the rotating arm at each set rotational position
- the structural motion state and the calibration result of the error parameter generate a dynamic correction matrix for correcting the actual projection result of the X-ray machine.
- the parameter calibration method, device and system for the X-ray machine provided by the embodiment are controlled by rotating the rotating arm of the X-ray machine around the standard phantom, and the projection result of the standard phantom at a plurality of rotational positions can be realized.
- the error parameters of the different mechanical structures of the rotating arm are calibrated, so that the method of calibrating the structural parameters of the system is usually calibrated according to the difference between the system parameters and the design target in the static state, and cannot be applied to the system.
- the problem of calibration and correction of the error caused by motion can be used to calibrate the error parameters generated during the movement of the X-ray machine.
- the number calibration provides a new idea. By using the error parameter calibration result to correct the actual projection result, the quality of the image after 3D reconstruction can be improved.
- FIG. 1 is a flow chart of a parameter calibration method of an X-ray machine according to Embodiment 1.
- 2a is a flow chart of a method for parameter calibration of an X-ray machine according to the second embodiment.
- FIG. 2b is a schematic diagram of a plurality of parameters for calculating a tangential displacement error according to the second embodiment.
- 2c is a schematic diagram of a plurality of parameters for calculating a radial displacement error provided by the second embodiment.
- FIG. 3a is a flowchart of a parameter calibration method of an X-ray machine according to Embodiment 3.
- Figure 3b is a top plan view of a standard phantom provided in the third embodiment.
- Figure 3c is a left side view of a standard phantom provided in the third embodiment.
- Figure 3d is a front elevational view of a standard phantom provided in the third embodiment.
- FIG. 3e is a schematic diagram of calculating a plurality of parameters of an acquisition angle in a projected image of a standard phantom according to the third embodiment.
- FIG. 4 is a flow chart of a method for parameter calibration of an X-ray machine according to Embodiment 4.
- FIG. 5 is a flowchart of a parameter calibration method of an X-ray machine according to Embodiment 5.
- FIG. 6 is a structural diagram of a parameter calibration apparatus for an X-ray machine according to Embodiment 6.
- FIG. 7 is a structural diagram of a parameter calibration system of an X-ray machine according to Embodiment 7.
- FIG. 8 is a flowchart of a calibration method corresponding to a parameter calibration system of an X-ray machine in an application scenario provided in the seventh embodiment.
- FIG. 9 is a structural diagram of a parameter calibration system of an X-ray machine in an application scenario provided in the seventh embodiment.
- This embodiment can be applied to calibrate parameters of an X-ray machine.
- X-ray machines can be used as an instrument for capturing images of human bodies.
- Parameter calibration is a method to reduce the system error. It can be that the pointer calibrates the error of the structural parameters of the system, and then a reasonable compensation amount can be set based on the error amount of the calibration to reduce the measurement error of the system. purpose.
- a standard human body model (hereinafter referred to as a standard phantom) is constructed, and the rotating arm of the X-ray machine is controlled to rotate around the standard phantom, and the projection result and ideal according to the standard phantom at a plurality of rotational positions are performed.
- the difference between the projection results can be used to calibrate the error parameters of the different mechanical structures of the rotating arm, so that the error caused by the X-ray machine motion can be calibrated.
- Embodiment 1 is a flow chart of a parameter calibration method of an X-ray machine according to Embodiment 1.
- the method of this embodiment can be applied to the case of calibrating an error generated during the rotation of the X-ray machine rotating arm.
- the method can be performed by a parameter calibration device of an X-ray machine, which can be implemented by software and/or hardware, and can generally be integrated into an X-ray machine.
- the method of this embodiment may include S110-S140.
- the X-ray machine rotating arm is controlled to perform a rotational scan around the standard phantom, and the projection results of the standard phantoms at the plurality of set rotational positions are acquired.
- the standard phantom may refer to a pre-configured, alternative to the actual scanning object whose shape conforms to the set rule, to calibrate the amount of error generated during the actual movement of the X-ray machine. Wherein, considering that the projection diameter of the spherical shape is the same in all directions, the standard phantom may be selected as a spherical phantom.
- the rotating arm of the X-ray machine may be a C-arm or an O-arm.
- the rotating arm can reach different rotation positions by controlling the rotation angle of the X-ray machine rotating arm.
- the rotation angle can be set to ⁇ /3, 5 ⁇ /6, 4 ⁇ /3, and 11 ⁇ /6, respectively. Wait.
- the projection result may refer to a projection image of the standard phantom obtained by the X-ray source rotating the X-ray machine rotating to the set rotation position and emitting X-rays by the X-ray source.
- the projection position of the X-ray source of the X-ray machine in the detector may be adjusted before the rotation of the X-ray machine rotating arm is controlled around the standard phantom. Overlaps the center of the detector and adjusts the geometric center of the standard phantom to overlap the design center of rotation of the X-ray machine's rotating arm.
- the detector can be a flat panel detector.
- the mechanical structure motion state may refer to a motion distance, a motion direction, or a motion speed corresponding to the one or more mechanical structures at a rotational position.
- the mechanical structure motion state can be represented by a vector.
- the mechanical structure motion state at one rotational position i can be represented by a vector s i
- the mechanical structural motion states corresponding to all the set rotational positions constitute a matrix S.
- the error parameter is calibrated according to the projection result of the standard phantom, and a calibration result of the error parameter is obtained.
- the projection result of the standard phantom should have a desired expected value.
- the error parameter in the rotation process may be determined by comparing the actual projection result of the standard phantom obtained during the actual rotation scanning process with the ideal expected value, wherein the error parameter may be actual.
- the difference between the projection result and the ideal expected value for example, the theoretical expected value of the center of rotation of the X-ray machine's rotating arm is F at a set rotation angle, and the actual projection center calculated by the actual projection result is H, Then (FH) can be used as an error parameter; the error parameter can also be a deviation value corresponding to the ideal expected value, for example, setting the ideal rotation angle of the X-ray machine rotating arm to 30°, determined by the actual projection result. If the actual rotation angle of the rotating arm is 31°, 31° can be used as an error parameter.
- the error parameters of the rotating arm in different mechanical structure motion states can be calibrated according to the difference between the actual projection result and the ideal expected value. Further, after the actual scanning object is placed on the X-ray machine, the projected image of the scanned object is corrected according to the mechanical structure motion state of the X-ray machine and the pre-calibrated error parameter.
- a dynamic correction matrix for correcting the actual projection result of the X-ray machine is generated based on the mechanical structure motion state of the rotating arm at each set rotational position and the calibration result of the error parameter.
- the set rotation angle of the rotating arm is limited. Therefore, it is impossible to calculate the corresponding error parameters in different mechanical structure states by exhausting the mechanical motion state of all the rotating arms.
- the rotating arm of the X-ray machine may be rotated to any angle, and the error parameter of the mechanical structure motion state corresponding to the angle may not be pre-calibrated. At this time, the X-ray machine's rotating arm cannot be corrected based on the pre-calibration result.
- an optimal estimation matrix can be constructed based on the error parameters of the motion state of the plurality of mechanical structures that have been calibrated, that is, the dynamic correction matrix is used to realize the error of the motion state of any one of the mechanical structures through the finite number of calibrated mechanical structure motion states.
- the parameters are estimated.
- the parameter calibration method of the X-ray machine provided by the embodiment is controlled by rotating the rotating arm of the X-ray machine around the standard phantom, and the movement state of the different mechanical structures of the rotating arm according to the projection result of the standard phantom at the plurality of rotating positions
- the calibration of the error parameters can solve the traditional calibration of the system structural parameters.
- the calibration method is usually based on the difference between the system parameters and the design objectives in the static state.
- the errors generated in the system motion cannot be calibrated and
- the revised problem can calibrate the error parameters generated in the X-ray machine motion, and provide a new idea for the parameter calibration of the X-ray machine.
- the error parameter calibration result By using the error parameter calibration result to calibrate the actual projection result, the three-dimensional can also be improved. The quality of the image after reconstruction.
- 2a is a flowchart of a method for calibrating a parameter of an X-ray machine according to the second embodiment.
- the standard body is selected as a spherical phantom, and the standard phantom is diametrically oriented.
- Two orthogonal alignment through holes may be included, and correspondingly, before the rotation of the rotating arm of the X-ray machine is controlled around the standard phantom to obtain a projection result of the standard phantom at a plurality of set rotation positions, Determining that a projection position of the X-ray source of the X-ray machine in the detector overlaps with a center of the detector, and a center of the standard phantom overlaps with a designed rotation center of the X-ray machine rotating arm;
- the error parameter may include: a tangential displacement error of the rotation center in the rotation plane and a radial displacement error, and based on the type of the error parameter, a calculation process of the error parameter is given.
- the method of this embodiment may include the following steps.
- the standard phantom is a spherical phantom, and the center of the standard phantom is adjusted such that the center of the standard phantom overlaps with the designed rotation center of the X-ray machine rotating arm, the standard The phantom may include two orthogonally aligned through holes in the diametrical direction.
- the standard phantom is placed in the C-arm imaging range, and the position of the standard phantom along a diameter direction is adjusted.
- the laser light emitted by the X-ray source passes through a through hole in the diameter direction of the standard phantom.
- the rotating arm is rotated by ⁇ /2 to adjust the position of the standard phantom along the other diameter direction, and the laser passes through the standard body.
- Another through hole in the diameter direction of the mold wherein the two through holes in the diameter direction are orthogonally aligned through holes. At this time, it can be determined that the center of the standard phantom is located at the design rotation center of the rotating arm.
- the X-ray machine rotating arm is controlled to perform a rotational scan around the standard phantom, and a plurality of projection results of the standard phantom at the set rotational position are acquired.
- a rotating scan period of the rotating arm can be set to 2 ⁇ .
- the set rotation angle can be multiple, so that the comprehensiveness of the set rotation angle can be ensured, and a relatively rich target projection result can be obtained.
- the standard phantom is a spherical phantom, so after the rotation of the rotating arm of the X-ray machine, the diameter of the standard phantom corresponds to the projected diameter.
- the projection center of the standard phantom has a position difference with the image center of the target projection result.
- the target tangential displacement error of the rotation center in the tangential direction of the rotation plane is calculated according to the projection diameter of the standard phantom, the actual diameter of the standard phantom, and the source plate distance of the rotating arm.
- Figure 2b is a schematic diagram of various parameters for calculating the target tangential displacement error, wherein the standard phantom has a projection diameter d 1 , the standard phantom has an actual diameter d 2 , and the rotating arm X-ray center A to the flat panel detector
- the distance between the distances (hereinafter referred to as the source plate distance of the rotating arm) is L.
- the target tangential displacement error of the projection point D of the rotating center C of the rotating arm in the tangential direction of the rotating plane is y i .
- point A is the center of the X-ray source
- B is the center of the standard phantom
- C is the center of rotation of the rotating arm
- D is the projection point of the center of rotation in the tangential direction of the plane of rotation
- point E is the DC connection and An intersection of the projection ray of the X-ray source
- F is the projection center of the standard phantom on the flat panel detector
- G is the intersection of the outermost ray tangential to the standard phantom and the flat panel detector of the rotating arm. Since the distance from the center of rotation C of the rotating arm to the center B of the standard phantom is small, that is, the value of y i is small, the length of DE is approximately d 2 /2.
- the center B of the standard phantom is placed on the rotating center C of the rotating arm.
- the distance from the center B of the standard phantom to the center A of the X source, that is, the length of AB, and the distance L from the source plate of the rotating arm is k, and the length of AD is kL+y i . According to the similarity principle of triangles, ⁇ AFG ⁇ ADE can be obtained. which is therefore
- FIG. 2c is a schematic diagram of a plurality of parameters for calculating a target radial displacement error, wherein the standard phantom has a projection diameter of d 1 , the standard phantom has an actual diameter of d 2 , and the rotating arm has a source plate distance of L, the target machine In the motion state, the target radial displacement error of the center of rotation in the tangential direction of the plane of rotation is x i .
- point A is the center of the X-ray source
- B is the center of the standard phantom
- C is the center of rotation.
- D1 is the projection point of the center of rotation in the radial direction of the plane of rotation
- D is the projection point of the center of rotation in the tangential direction of the plane of rotation
- point H is the intersection of the extension line of the AC line and the flat panel detector of the rotating arm, ie At the center of the target projection result
- F is the projection center of the standard phantom
- G is the intersection of the outermost ray tangential to the standard phantom and the flat panel detector.
- the length of the FH is the position difference X between the standard phantom projection center and the target projection result center. According to the similarity principle of the triangle, ⁇ AFH ⁇ ADC can be obtained. which is Combine therefore
- the technical solution of the embodiment obtains a dynamic correction matrix for correcting the rotation center of the actual projection result of the X-ray machine by calibrating the tangential displacement error and the radial displacement error of the rotation center of the rotating arm under different mechanical structural motion states,
- the dynamic correction matrix By applying the dynamic correction matrix, the tangential displacement error and the radial displacement error of the actual rotation center of the projection result can be corrected, and the image quality after three-dimensional imaging can be improved.
- FIG. 3a is a flowchart of a method for calibrating a parameter of an X-ray machine according to Embodiment 3.
- This embodiment is based on the above embodiment, and two standard intersecting with any one of the aligned through holes are disposed on the standard phantom.
- the marking through holes are symmetrical with respect to the center of the standard phantom and perpendicular to the plane defined by the two aligned through holes.
- the parameter calibration method of the X-ray machine may further include: In each setting At a fixed rotational position, the rotational angle of the rotating arm is recorded; the error parameter may be an acquisition angle included in a different mechanical structural motion state of the rotating arm, and the error parameter is calculated based on the type of the error parameter. process.
- the method of this embodiment may include the following steps.
- the standard phantom is a spherical phantom
- the standard phantom includes two orthogonally aligned through holes in a diameter direction, optionally, and is disposed on the standard phantom
- Two mark through holes intersecting the through holes are aligned, the two mark through holes being symmetric with respect to the center of the standard phantom and perpendicular to the plane defined by the two aligned through holes, respectively.
- the distance from the center of the sphere can be set to be 1/2 radius from the center of the sphere.
- Figure 3b shows a top view of a standard phantom
- Figure 3c shows a left side view of a standard phantom
- Figure 3d shows a front view of a standard phantom with alignment through holes 31 and markings Through hole 32.
- the X-ray machine rotating arm is controlled to perform a rotational scan around the standard phantom to obtain a projection result of a plurality of standard phantoms at the set rotational position.
- the rotation angle of the rotating arm may be a set rotation angle input to a central processing unit of the X-ray machine, and based on the set rotation angle, the X-ray machine can control the rotating arm to perform corresponding mechanical motion. In order to adjust the rotating arm to the set rotation angle, and then after the rotation of the rotating arm is finished, the corresponding mechanical structure motion parameters are recorded.
- a target acquisition angle in a moving state of the target mechanical structure is calculated according to a projection distance between the two marked through holes of the standard phantom, the target rotation angle, and a projection diameter of the standard phantom.
- 3e is a schematic diagram of calculating a plurality of parameters of a target acquisition angle in a projection image of a standard phantom, wherein a projection distance between two standard through-holes 32 of the standard phantom is a, an acquisition angle value is ⁇ , a standard phantom
- the projection diameter is d 1 .
- the rotation angle values corresponding to the four quadrants are ⁇ , ⁇ + ⁇ /2, ⁇ , respectively. + ⁇ , ⁇ + 3 ⁇ /2, the set of four angles as the acquisition angle.
- the acquisition angle ⁇ corresponding to the rotation angle ⁇ in one rotation scan period is:
- ⁇ ⁇ cos -1 (2a/d 1 ), (cos -1 (2a/d 1 ))+ ⁇ /2), (cos -1 (2a/d 1 ))+ ⁇ ), (cos -1 ( (2a/d 1 )) +3 ⁇ /2 ⁇ .
- the acquisition angle that minimizes the difference between the acquisition angle set and the target rotation angle is used as the calibration result of the target acquisition angle.
- the target rotation angle is When the target acquisition angle is corrected, That is, the calibration result of the target acquisition angle is the acquisition angle value with the smallest difference from the target rotation angle, thereby ensuring that the obtained acquisition angle value is closest to the target rotation angle, that is, the target acquisition angle is consistent with the actual rotation angle.
- a dynamic correction matrix for correcting the actual projection result of the X-ray machine is generated according to the calibration result of the target acquisition angle and the mechanical structure motion state of the rotating arm.
- the calibration result ⁇ i according to the target acquisition angle constitutes an error matrix P, which is a matrix composed of ⁇ i .
- the mechanical structure motion state s i at all set rotational positions is made up of a mechanical structure matrix S.
- Perform singular value decomposition on the error matrix to obtain P EW, W is the weight matrix, and E is the eigenvalue matrix.
- the technical solution of the embodiment calculates the target acquisition angle under the motion state of the target mechanical structure by the projection distance between the two calibrated through holes of the standard phantom, the target rotation angle and the projection diameter of the standard phantom, and the calibration according to the target acquisition angle.
- a dynamic correction matrix for correcting the actual projection result of the X-ray machine is generated, and the actual correction result can be corrected by using the dynamic correction matrix when the corresponding error parameter is used for the acquisition angle.
- the quality of the image is a dynamic correction matrix for correcting the actual projection result of the X-ray machine is generated, and the actual correction result can be corrected by using the dynamic correction matrix when the corresponding error parameter is used for the acquisition angle.
- the dynamic correction matrix for correcting the actual projection result of the X-ray machine may further comprise: forming a mechanical structure motion state at a plurality of set rotation positions into a mechanical structure matrix S; and moving the plurality of mechanical structures in a state of motion Error parameters constitute an error matrix P;
- the method in this embodiment may include:
- the X-ray machine rotating arm is controlled to perform a rotational scan around the standard phantom to obtain projection results of a plurality of standard phantoms at the set rotational position.
- the error parameter is calibrated according to the projection result of the standard phantom, and a calibration result of the error parameter is obtained.
- error parameters may be a radial displacement error x i, y i tangential displacement error calibration results and the acquisition angle ⁇ i
- the configuration E(WS T (SS T ) -1 ) is constructed as a dynamic correction matrix for correcting the actual projection result of the X-ray machine.
- S is a matrix of motion states s i of all mechanical structures.
- the technical solution of the embodiment constructs a dynamic correction matrix by using error parameters, which can correct the actual projection result and improve the quality of the image after three-dimensional imaging.
- FIG. 5 is a flowchart of a parameter calibration method for an X-ray machine according to Embodiment 5, on the basis of the foregoing embodiment, the parameter calibration method of the X-ray machine can further increase a dynamic repair matrix generated according to the calibration, and the actual The process of correcting the projection result.
- the method of this embodiment may include the following steps.
- the X-ray machine rotating arm is controlled to perform a rotational scan around the standard phantom, and a plurality of projection results of the standard phantom at the set rotational position are acquired.
- the error parameter is calibrated according to the projection result of the standard phantom, and a calibration result of the error parameter is obtained.
- the X-ray machine rotating arm is controlled to perform a rotation scan around the object to be measured, and the projection result to be corrected at the actual rotation position is obtained.
- an actual error parameter corresponding to the motion state of the measured mechanical structure is generated according to the measured mechanical structure motion state and the dynamic correction matrix.
- the actual error parameter may include: a measured target tangential displacement error y j of the rotation center on the rotation plane and a measured target radial displacement error x j under the measured mechanical structure motion state s j of the rotating arm;
- the actual error parameter, and performing image correction on the projection result to be corrected may include the following steps.
- the target correction target projection result is subjected to a scaling process of a set ratio.
- the measured target tangential displacement error wherein, the projection diameter d 1 of the standard phantom, the actual diameter d 2 of the standard phantom, and the source plate distance L, k of the rotating arm are the distance from the center of the sphere of the standard phantom to the light source, that is, the length of the AB, and the distance from the source plate.
- the ratio of L, the ratio of the actual image of the standard phantom to the projected image make Among them, ideally Therefore, the projection result of the target to be corrected can be multiplied as The scaling operation to eliminate scaling errors.
- k 1/2
- the target radial displacement error Where X is the position difference between the standard phantom projection center and the target projection result center, the projection diameter d 1 of the standard phantom, and the actual diameter d 2 of the standard phantom, the difference between the standard phantom projection center and the target projection result center Therefore, by performing the projection result of the target to be corrected Pan, you can eliminate the translation error.
- the actual error parameter may include an actual acquisition angle in a state of motion of the measured mechanical structure of the rotating arm.
- performing image correction on the to-be-corrected projection result according to the actual error parameter may include: performing three-dimensional reconstruction on each to-be-corrected projection result according to the actual acquisition angle in different measured mechanical structure motion states.
- the three-dimensional reconstruction is performed according to the corrected image and the corresponding angle ⁇ j to improve the image quality after the three-dimensional reconstruction.
- the target tangential displacement error is used to perform a scaling process on the target to be corrected projection result; according to the target radial displacement error, the target distance to be corrected is subjected to a set distance translation, according to different actual measurements.
- the actual acquisition angle of the mechanical structure in motion state, three-dimensional reconstruction of each projection result to be corrected, can achieve image correction, and achieve the effect of improving the image quality after three-dimensional reconstruction.
- FIG. 6 is a structural diagram of a parameter calibration apparatus for an X-ray machine according to Embodiment 6.
- the parameter calibration device of the X-ray machine provided by the embodiment of the present disclosure can be applied to the parameter calibration system of the X-ray machine described in the embodiment of the present disclosure.
- the apparatus includes a projection result acquisition module 610, a motion state recording module 620, an error parameter calibration module 630, and a dynamic correction matrix generation module 640.
- the projection result obtaining module 610 is configured to control the X-ray machine rotating arm to perform a rotation scan around the standard phantom, and collect the projection results at the plurality of set rotation positions.
- the motion state recording module 620 is configured to record the mechanical structural motion state of the rotating arm at each of the set rotational positions.
- the error parameter calibration module 630 is configured to calibrate the error parameter according to the projection result of the standard phantom at each set rotation position to obtain a calibration result of the error parameter.
- the dynamic correction matrix generating module 640 is configured to generate a correction for the actual projection result of the X-ray machine according to the mechanical structure motion state of the rotating arm at each set rotational position and the calibration result of the error parameter Dynamic correction matrix.
- the parameter calibration device of the X-ray machine provided by the embodiment is controlled to rotate around the standard phantom by controlling the rotating arm of the X-ray machine, and the mechanical structure of the rotating arm is moved according to the projection result of the standard phantom at a plurality of rotating positions.
- the calibration of the error parameters can solve the problem that the calibration of the system structural parameters is usually based on the difference between the system parameters and the design objectives in the static state, and the error caused by the system motion cannot be calibrated and corrected. It can realize the calibration of the error parameters generated in the X-ray machine motion, and provide a new idea for the parameter calibration of the X-ray machine. By using the error parameter calibration result to correct the actual projection result, the image quality after 3D reconstruction is improved.
- the standard phantom is a spherical phantom, and the standard phantom includes two orthogonally aligned through holes in the diametrical direction.
- the parameter calibration device of the X-ray machine may further include: a projection position determining module configured to perform a rotational scan around the standard phantom of the rotating arm of the X-ray machine to obtain a standard phantom at a plurality of set rotational positions; Before the projection result, determining that the projection position of the X-ray source of the X-ray machine in the detector overlaps with the center of the detector, the spherical center of the standard phantom and the design rotation of the X-ray machine rotating arm The center overlaps.
- the error parameter may include a tangential displacement error of the center of rotation on the plane of rotation and a radial displacement error in each mechanical structural motion state of the rotating arm.
- the error parameter calibration module 630 can be configured to: acquire a target projection result corresponding to the set rotation position i and a target mechanical structure motion state, i ⁇ [1, N], N is a rotation scan of the rotating arm a total number of set rotation angles through which the cycle passes; in the target projection result, a projection diameter of the standard phantom, and a position difference between a projection center of the standard phantom and an image center of the target projection result; a projection diameter of the standard phantom, an actual diameter of the standard phantom, and a source plate distance of the rotating arm, and calculating a target tangential displacement error of the rotation center in a tangential direction of the rotation plane; projection according to the standard phantom The diameter, the actual diameter of the standard phantom, and the position difference between the projection center of the standard phantom and the image center of the target projection result, and the target radial displacement error of the rotation center in the radial direction of the rotation plane is calculated.
- two standard through holes intersecting any one of the aligned through holes are disposed on the standard phantom, the two labeled through holes being symmetric with respect to a center of the standard phantom and perpendicular to the two alignments respectively The plane defined by the through hole.
- the motion state recording module 620 can also be configured to record the rotation angle of the rotating arm at each set rotation position.
- the error parameter may include an acquisition angle in a different mechanical structural motion state of the rotating arm.
- the error parameter calibration module 630 can be configured to: acquire a target projection result corresponding to the set rotation position i, a target mechanical structure motion state, and a target rotation angle of the rotating arm, 1 ⁇ [1, N], N is a total number of set rotation angles of a rotation scan period of the rotating arm; in the target projection result, a projection diameter of the standard phantom is acquired, and between the two mark through holes of the standard phantom Projection distance; calculating a target acquisition angle in a moving state of the target mechanical structure according to a projection distance between the two marked through holes of the standard phantom, the target rotation angle, and a projection diameter of the standard phantom.
- the dynamic correction matrix generating module 640 may be configured to: compose a mechanical structure motion state at a plurality of set rotation positions into a mechanical structure matrix S; and combine error parameters in a plurality of mechanical structure motion states into an error matrix P
- the parameter calibration device of the X-ray machine may further include: a projection result obtaining module to be corrected, configured to be calibrated according to a mechanical mechanism motion state of the rotating arm at each set position and the error parameter
- a projection result obtaining module to be corrected configured to be calibrated according to a mechanical mechanism motion state of the rotating arm at each set position and the error parameter
- the measured mechanical structure motion state recording module is configured to generate an actual projection result for the X-ray machine based on a calibration result of a mechanical mechanism motion state of the rotating arm at each set position and the error parameter After the modified dynamic correction matrix, the measured mechanical motion state of the rotating arm is recorded at the actual rotational position.
- the actual error parameter production module is configured to generate an actual error parameter corresponding to the motion state of the measured mechanical structure according to the measured mechanical structure motion state and the dynamic correction matrix.
- the projection result correction module to be corrected is configured to perform image correction on the projection result to be corrected according to the actual error parameter.
- the actual error parameter may include: an actual tangential displacement error of the center of rotation of the rotating arm in the plane of rotation and an actual radial displacement error at the actual rotational position;
- the actual error parameter production module may be configured to: acquire a target to be corrected projection result corresponding to the actual rotation position j, the actual mechanical structure motion state of the nuclear target, j ⁇ [1, M], M is a rotation scan period of the rotating arm The total number of actual rotation angles passed; according to the dynamic correction matrix, acquiring a target tangential displacement error and a target radial displacement error corresponding to the moving state of the target measured mechanical structure; and waiting for the target according to the target tangential displacement error Correcting the projection result to perform scaling processing of the set ratio; And shifting the set distance to the target to be corrected projection result according to the target radial displacement error.
- the actual error parameter may include an actual acquisition angle in a state of motion of the measured mechanical structure of the rotating arm.
- the actual error parameter production module may be further configured to perform three-dimensional reconstruction on each of the projections to be corrected according to the actual acquisition angles in different motion states of the mechanical structure.
- the parameter calibration device of the X-ray machine provided in this embodiment can be used to perform the parameter calibration method of the X-ray machine provided by any of the above embodiments, and has a corresponding function module.
- FIG. 7 is a structural diagram of a parameter calibration system for an X-ray machine according to Embodiment 7. As shown in FIG. 7, the system may include: a rotating arm 710, an X-ray source 720 disposed opposite the rotating arm, and a detector 730, a processor 740, an image capturing device 750, a motion state collecting component 760, and Standard phantom 770.
- the rotating arm 710 is configured to drive the X-ray source and the flat panel detector to rotate a plurality of set rotation angles around the standard phantom according to a control signal sent by the processor.
- the X-ray source 720 is configured to transmit an X-ray scan signal to the standard phantom.
- the detector 730 is configured to collect projection data of the X-ray scan signal on the standard phantom.
- the image acquisition device 750 is configured to generate a projection result of the standard phantom according to the projection data collected by the detector.
- the motion state collecting component 760 is configured to acquire a mechanical structural motion state of the rotating arm.
- the processor 740 is configured to control the rotating arm to perform a rotational scan around the standard phantom to acquire a projection result of setting a standard phantom at a plurality of rotational positions; at each set rotational position, the rotating arm is recorded The mechanical structure motion state; at each set rotation position, the error parameter is proofed according to the projection result of the standard phantom, and the proofreading result of the error parameter is obtained; according to each set rotation position A mechanical correction state of the rotating arm and a calibration result of the error parameter generate a dynamic correction matrix for correcting the actual projection result of the X-ray machine.
- the standard phantom 770 is a spherical phantom, and the standard phantom 770 includes two orthogonally aligned through holes in a diameter direction.
- the processor 740 may be further configured to: determine an X-ray source of the X-ray machine before controlling a rotation scan of the X-ray machine around the standard phantom 770 to obtain a projection result of the standard phantom at different rotational positions
- the projected position of 720 in detector 730 overlaps the center of detector 730
- the center of the standard phantom 770 overlaps the designed center of rotation of the X-ray machine's rotating arm.
- two standard through holes intersecting any one of the aligned through holes are disposed on the standard phantom 770, and the two labeled through holes are symmetric with respect to a center of a standard phantom, and are perpendicular to the two pairs respectively.
- the plane defined by the through hole is disposed on the standard phantom 770, and the two labeled through holes are symmetric with respect to a center of a standard phantom, and are perpendicular to the two pairs respectively. The plane defined by the through hole.
- the processor 740 may be further configured to record a rotation angle of the rotating arm at each set rotational position.
- FIG. 8 is a flow chart of a calibration method corresponding to a parameter calibration system of an X-ray machine, which may include:
- Initial alignment may refer to aligning the center of the standard phantom and X-ray source with the center of the flat panel detector.
- Exposure imaging can refer to a rotational scan of a standard phantom.
- the actual projection result of the measured object is corrected to obtain a more realistic high quality image.
- FIG. 9 is a structural diagram of a parameter calibration system of an X-ray machine.
- the parameter calibration system of the X-ray machine may include: an initial alignment module 910, which is set as an initial alignment of a standard phantom; and an image acquisition module 920. And set to collect the actual projection result of the object to be corrected; the system calibration module 930 is configured as an error parameter of the calibration system; and the image correction module 940 is configured to correct the actual projection result of the object to be corrected according to the error parameter;
- the human-computer interaction module 950 is set to human-computer interaction.
- the parameter calibration system of the X-ray machine provided by the embodiment is controlled by rotating the rotating arm of the X-ray machine around the standard phantom, and the different mechanical parts of the rotating arm according to the projection result of the standard phantom at the plurality of rotating positions Calibration of the error parameters in the structural motion state can be solved.
- the method of calibrating the system structural parameters is usually calibrated according to the difference between the system parameters and the design target in the static state, and the error caused by the system motion cannot be calibrated.
- the correction problem can calibrate the error parameters generated in the X-ray machine motion, and provide a new idea for the parameter calibration of the X-ray machine.
- the present disclosure provides a parameter calibration method, device and system for an X-ray machine, which can realize X-ray
- the error parameters generated during the motion of the machine are calibrated.
- the quality of the image after three-dimensional reconstruction can be improved.
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (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)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Description
本公开涉及信号处理技术,例如涉及一种X光机的参数校准方法、装置及系统。The present disclosure relates to signal processing techniques, such as a method, apparatus, and system for parameter calibration of an X-ray machine.
在X光三维图像重建中,C型臂光源和探测器的运行轨迹是否符合预设轨迹会影响重建图像的质量。C型臂系统运行过程中机械结构的抖动会使射线源至探测器的距离、射线源至旋转中心的距离、旋转中心位置等重要参数发生变化,导致投影图像出现缩放或是平移;或是出现投影采集角度与预设角度不同的情况,最终导致重建图像中可能出现伪影,影响图像重建质量。而机械结构的抖动缺乏观测手段,使得对这些因素的校准和修正较为困难。In the X-ray three-dimensional image reconstruction, whether the running trajectory of the C-arm light source and the detector conforms to the preset trajectory will affect the quality of the reconstructed image. The vibration of the mechanical structure during the operation of the C-arm system changes the important parameters such as the distance from the source to the detector, the distance from the source to the center of rotation, and the position of the center of rotation, resulting in scaling or translation of the projected image; The difference between the projection angle and the preset angle may eventually lead to artifacts in the reconstructed image, which may affect the image reconstruction quality. The lack of observation of the vibration of the mechanical structure makes it difficult to calibrate and correct these factors.
针对系统参数与设计目标之间的差异,可以采用针对静态状态下系统参数与设计目标之间的差异而进行的参数校准。For the difference between system parameters and design goals, parameter calibration for the difference between system parameters and design goals in static state can be used.
发明人发现对系统结构参数进行校准的方法是根据静态状态下系统参数与设计目标之间的差异而进行的校准,不能对系统运动产生的误差进行校准和修正。The inventors found that the method of calibrating the structural parameters of the system is based on the difference between the system parameters and the design goals in the static state, and the errors caused by the motion of the system cannot be calibrated and corrected.
发明内容Summary of the invention
本公开提供一种X光机的参数校准方法、装置及系统,可以实现对X光机运动中产生的误差进行校准。The present disclosure provides a method, device and system for parameter calibration of an X-ray machine, which can calibrate errors generated in the movement of the X-ray machine.
一种X光机的参数校准方法,该方法可以包括:A parameter calibration method for an X-ray machine, the method may include:
控制X光机旋转臂围绕标准体模进行旋转扫描,采集多个设定旋转位置处的标准体模的投影结果;Controlling the X-ray machine rotating arm to perform a rotation scan around the standard phantom, and collecting projection results of a plurality of standard phantoms at the set rotation position;
在每个设定旋转位置处,记录所述旋转臂的机械结构运动状态;Recording a mechanical structural motion state of the rotating arm at each set rotational position;
在每个设定旋转位置处,根据所述标准体模的投影结果,对所述误差参数进行校准,得到误差参数的校准结果;以及At each set rotation position, the error parameter is calibrated according to the projection result of the standard phantom to obtain a calibration result of the error parameter;
根据在每个设定旋转位置处的旋转臂的机械结构运动状态和所述误差参数的校准结果,生成用于对所述X光机的实际投影结果进行修正的动态修正矩阵。A dynamic correction matrix for correcting the actual projection result of the X-ray machine is generated based on the mechanical structural motion state of the rotating arm at each set rotational position and the calibration result of the error parameter.
一种X光机的参数校准装置,可用于执行本实施例所述的X光机的参数校 准方法,该装置可以包括:A parameter calibration device for an X-ray machine, which can be used to perform parameter calibration of the X-ray machine described in this embodiment A quasi-method, the device can include:
投影结果采集模块,设置为控制X光机旋转臂围绕标准体模进行旋转扫描,采集多个设定旋转位置处的标准体模的投影结果;The projection result acquisition module is configured to control the X-ray machine rotating arm to perform a rotation scan around the standard phantom, and collect the projection results of the standard phantom at a plurality of set rotation positions;
运动状态记录模块,设置为在每个设定旋转位置处,记录所述旋转臂的机械结构运动状态;a motion state recording module configured to record a mechanical structural motion state of the rotating arm at each set rotational position;
误差参数校准模块,设置为在每个设定旋转位置处,根据所述标准体模的投影结果,对所述误差参数进行校准,得到误差参数的校准结果;以及The error parameter calibration module is configured to, at each set rotation position, calibrate the error parameter according to a projection result of the standard phantom to obtain a calibration result of the error parameter;
动态修正矩阵生成模块,设置为根据在每个设定旋转位置处的旋转臂的机械结构运动状态和所述误差参数的校准结果,生成用于对所述X光机的实际投影结果进行修正的动态修正矩阵。a dynamic correction matrix generating module configured to generate a correction for correcting an actual projection result of the X-ray machine according to a mechanical structural motion state of the rotating arm at each set rotational position and a calibration result of the error parameter Dynamic correction matrix.
一种X光机的参数校准系统,该系统可以包括:A parameter calibration system for an X-ray machine, the system may include:
旋转臂、相对设置于所述旋转臂之上的X光光源以及探测器、处理器、图像采集器件、运动状态采集部件以及标准体模;a rotating arm, an X-ray source disposed opposite the rotating arm, and a detector, a processor, an image acquisition device, a motion state acquisition component, and a standard phantom;
所述旋转臂,设置为根据所述处理器发送的控制信号带动所述X光光源以及探测器围绕所述标准体模进行多个设定旋转角度的旋转;The rotating arm is configured to drive the X-ray light source and the detector to rotate a plurality of set rotation angles around the standard phantom according to a control signal sent by the processor;
所述X光光源,设置为向所述标准体模发送X光扫描信号;The X-ray source is configured to send an X-ray scan signal to the standard phantom;
所述探测器,设置为采集所述X光扫描信号在所述标准体模上的投影数据;The detector is configured to collect projection data of the X-ray scan signal on the standard phantom;
所述图像采集器件,设置为根据所述探测器采集的投影数据,生成标准体模的投影结果;The image acquisition device is configured to generate a projection result of a standard phantom according to the projection data collected by the detector;
所述运动状态采集部件,设置为获取所述旋转臂的机械结构运动状态;The motion state collecting component is configured to acquire a mechanical structural motion state of the rotating arm;
所述处理器,设置为控制所述旋转臂围绕标准体模进行旋转扫描,获取多个设定旋转位置处的投影结果;在每个设定旋转位置处,记录所述旋转臂的机械结构运动状态;在每个设定旋转位置处,根据所述标准体模的投影结果,对所述误差参数进行校准,得到误差参数的校准结果;根据在每个设定旋转位置处的旋转臂的机械结构运动状态和所述误差参数的校准结果,生成用于对所述X光机的实际投影结果进行修正的动态修正矩阵。The processor is configured to control the rotating arm to perform a rotational scan around the standard phantom to obtain a projection result at a plurality of set rotational positions; at each set rotational position, record the mechanical structural motion of the rotating arm a state; at each set rotational position, the error parameter is calibrated according to a projection result of the standard phantom to obtain a calibration result of the error parameter; according to the mechanism of the rotating arm at each set rotational position The structural motion state and the calibration result of the error parameter generate a dynamic correction matrix for correcting the actual projection result of the X-ray machine.
本实施例提供的X光机的参数校准方法、装置及系统,通过控制X光机旋转臂围绕标准体模进行旋转扫描,根据标准体模在多个旋转位置处的的投影结果,可以实现对旋转臂的不同机械结构运动状态下的误差参数进行校准,从而可以解决在对系统结构参数进行校准的方法通常是根据静态状态下系统参数与设计目标之间的差异而进行的校准,不能对系统运动产生的误差进行校准和修正的问题,可以实现对X光机运动中产生的误差参数进行校准,为X光机的参 数校准提供一种新的思路,通过利用误差参数校准结果对实际投影结果进行修正,可以提高三维重建后图像的质量。The parameter calibration method, device and system for the X-ray machine provided by the embodiment are controlled by rotating the rotating arm of the X-ray machine around the standard phantom, and the projection result of the standard phantom at a plurality of rotational positions can be realized. The error parameters of the different mechanical structures of the rotating arm are calibrated, so that the method of calibrating the structural parameters of the system is usually calibrated according to the difference between the system parameters and the design target in the static state, and cannot be applied to the system. The problem of calibration and correction of the error caused by motion can be used to calibrate the error parameters generated during the movement of the X-ray machine. The number calibration provides a new idea. By using the error parameter calibration result to correct the actual projection result, the quality of the image after 3D reconstruction can be improved.
图1是本实施例一提供的一种X光机的参数校准方法的流程图。1 is a flow chart of a parameter calibration method of an X-ray machine according to Embodiment 1.
图2a是本实施例二提供的一种X光机的参数校准方法的流程图。2a is a flow chart of a method for parameter calibration of an X-ray machine according to the second embodiment.
图2b是本实施例二提供的计算切向位移误差的多个参数的示意图。FIG. 2b is a schematic diagram of a plurality of parameters for calculating a tangential displacement error according to the second embodiment.
图2c是本实施例二提供的计算径向位移误差的多个参数的示意图。2c is a schematic diagram of a plurality of parameters for calculating a radial displacement error provided by the second embodiment.
图3a是本实施例三提供的一种X光机的参数校准方法的流程图。FIG. 3a is a flowchart of a parameter calibration method of an X-ray machine according to Embodiment 3.
图3b是本实施例三提供的一种标准体模的俯视图。Figure 3b is a top plan view of a standard phantom provided in the third embodiment.
图3c是本实施例三提供的一种标准体模的左视图。Figure 3c is a left side view of a standard phantom provided in the third embodiment.
图3d是本实施例三提供的一种标准体模的主视图。Figure 3d is a front elevational view of a standard phantom provided in the third embodiment.
图3e是本实施例三提供的在标准体模的投影图像中计算采集角度的多个参数的示意图。FIG. 3e is a schematic diagram of calculating a plurality of parameters of an acquisition angle in a projected image of a standard phantom according to the third embodiment.
图4是本实施例四提供的一种X光机的参数校准方法的流程图。4 is a flow chart of a method for parameter calibration of an X-ray machine according to Embodiment 4.
图5是本实施例五提供的一种X光机的参数校准方法的流程图。FIG. 5 is a flowchart of a parameter calibration method of an X-ray machine according to Embodiment 5.
图6是本实施例六提供的一种X光机的参数校准装置的结构图。FIG. 6 is a structural diagram of a parameter calibration apparatus for an X-ray machine according to Embodiment 6.
图7是本实施例七提供的一种X光机的参数校准系统的结构图。FIG. 7 is a structural diagram of a parameter calibration system of an X-ray machine according to Embodiment 7.
图8是本实施例七提供的一种应用场景中X光机的参数校准系统对应的校准方法的流程图。FIG. 8 is a flowchart of a calibration method corresponding to a parameter calibration system of an X-ray machine in an application scenario provided in the seventh embodiment.
图9是本实施例七提供的一种应用场景中X光机的参数校准系统的结构图。FIG. 9 is a structural diagram of a parameter calibration system of an X-ray machine in an application scenario provided in the seventh embodiment.
下面结合附图和实施例对本公开作相关说明。可以理解的是,此处所描述的实施例仅仅用于解释本公开,而非对本公开的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本公开相关的部分而非全部结构。在不冲突的情况下,以下实施例和实施例中的特征可以相互组合。The present disclosure will be described in conjunction with the accompanying drawings and embodiments. It is understood that the embodiments described herein are merely illustrative of the disclosure and are not intended to be limiting. In addition, it should be noted that, for the convenience of description, only some but not all of the structures related to the present disclosure are shown in the drawings. The features of the following embodiments and embodiments may be combined with each other without conflict.
另外还需要说明的是,为了便于描述,附图中仅示出了与本公开相关的部分而非全部内容。在讨论示例性实施例之前应当提到的是,一些示例性实施例被描述成作为流程图描绘的处理或方法。虽然流程图将操作(或步骤)描述成顺序的处理,但是许多操作可以被并行地、并发地或者同时实施。此外,操作的顺序可以被重新安排。当其操作完成时所述处理可以被终止,但是还可以具 有未包括在附图中的附加步骤。所述处理可以对应于方法、函数、规程、子例程、子程序等等。It should also be noted that, for the convenience of description, only some, but not all, of the present disclosure are shown in the drawings. It should be noted that prior to discussing the exemplary embodiments, some exemplary embodiments are described as a process or method depicted as a flowchart. Although a flowchart depicts operations (or steps) as a sequential process, many of the operations can be performed in parallel, concurrently or concurrently. In addition, the order of operations can be rearranged. The process can be terminated when its operation is completed, but can also have There are additional steps not included in the figures. The processing may correspond to methods, functions, procedures, subroutines, subroutines, and the like.
本实施例可以应用于对X光机的参数进行校准。在医学中,X光机可以作为一种采集人体图像的仪器。参数校准是一种减小系统误差的方法,可以是指针对系统的结构参数对误差量进行校准,进而可以基于该校准的误差量设定合理的补偿量,以达到减小系统的测量误差的目的。This embodiment can be applied to calibrate parameters of an X-ray machine. In medicine, X-ray machines can be used as an instrument for capturing images of human bodies. Parameter calibration is a method to reduce the system error. It can be that the pointer calibrates the error of the structural parameters of the system, and then a reasonable compensation amount can be set based on the error amount of the calibration to reduce the measurement error of the system. purpose.
发明人通过研究发现:参数校准的方法通常是在静态系统下对结构参数进行校准,不能实现对系统运动产生的误差进行校准。The inventors found through research that the method of parameter calibration is usually to calibrate the structural parameters in a static system, and the error caused by the motion of the system cannot be calibrated.
在本实施例中,构造了一个标准人体模型(以下简称标准体模),控制X光机旋转臂围绕该标准体模进行旋转扫描,根据标准体模在多个旋转位置处的投影结果与理想的投影结果之间的差异,可以实现对旋转臂的不同机械结构运动状态下的误差参数进行校准,从而可以实现对X光机运动产生的误差进行校准。In this embodiment, a standard human body model (hereinafter referred to as a standard phantom) is constructed, and the rotating arm of the X-ray machine is controlled to rotate around the standard phantom, and the projection result and ideal according to the standard phantom at a plurality of rotational positions are performed. The difference between the projection results can be used to calibrate the error parameters of the different mechanical structures of the rotating arm, so that the error caused by the X-ray machine motion can be calibrated.
实施例一Embodiment 1
图1是本实施例一提供的一种X光机的参数校准方法的流程图。本实施例的方法可以适用于对X光机旋转臂旋转过程中产生的误差进行校准的情况。该方法可以由X光机的参数校准装置执行,该装置可由软件和/或硬件的方式实现,并一般可集成于X光机中。如图1所示,本实施例的方法可以包括S110-S140。1 is a flow chart of a parameter calibration method of an X-ray machine according to Embodiment 1. The method of this embodiment can be applied to the case of calibrating an error generated during the rotation of the X-ray machine rotating arm. The method can be performed by a parameter calibration device of an X-ray machine, which can be implemented by software and/or hardware, and can generally be integrated into an X-ray machine. As shown in FIG. 1, the method of this embodiment may include S110-S140.
在S110中,控制X光机旋转臂围绕标准体模进行旋转扫描,采集多个设定旋转位置处的标准体模的投影结果。In S110, the X-ray machine rotating arm is controlled to perform a rotational scan around the standard phantom, and the projection results of the standard phantoms at the plurality of set rotational positions are acquired.
所述标准体模可以是指预先构造的,形状符合设定规则的实际扫描对象的替代物,用以校准X光机实际运动过程中产生的误差量。其中,考虑到球形在所有方向上的投影直径均相同,所述标准体模可选为球形体模。The standard phantom may refer to a pre-configured, alternative to the actual scanning object whose shape conforms to the set rule, to calibrate the amount of error generated during the actual movement of the X-ray machine. Wherein, considering that the projection diameter of the spherical shape is the same in all directions, the standard phantom may be selected as a spherical phantom.
可选地,X光机的旋转臂可以是C型臂,也可以是O型臂。其中,可以通过控制X光机旋转臂的旋转角度来使旋转臂到达不同的旋转位置,示例性的,旋转角度可以分别设定为π/3、5π/6、4π/3、以及11π/6等。Alternatively, the rotating arm of the X-ray machine may be a C-arm or an O-arm. Wherein, the rotating arm can reach different rotation positions by controlling the rotation angle of the X-ray machine rotating arm. For example, the rotation angle can be set to π/3, 5π/6, 4π/3, and 11π/6, respectively. Wait.
所述投影结果可以是指所述X光机旋转臂旋转至设定旋转位置处,通过X光源发射X光获取的该标准体模的投影图像。The projection result may refer to a projection image of the standard phantom obtained by the X-ray source rotating the X-ray machine rotating to the set rotation position and emitting X-rays by the X-ray source.
可选地,为了保证后续误差参数校准的准确性,在控制X光机旋转臂围绕标准体模进行旋转扫描之前,可以将所述X光机的X光光源在探测器中的投影位置调整至与所述探测器的中心重叠,并将标准体模的几何中心调整至与所述X光机旋转臂的设计旋转中心重叠。可选地,该探测器可为平板探测器。 Optionally, in order to ensure the accuracy of the subsequent error parameter calibration, the projection position of the X-ray source of the X-ray machine in the detector may be adjusted before the rotation of the X-ray machine rotating arm is controlled around the standard phantom. Overlaps the center of the detector and adjusts the geometric center of the standard phantom to overlap the design center of rotation of the X-ray machine's rotating arm. Alternatively, the detector can be a flat panel detector.
在S120中,在每个设定旋转位置处,记录所述旋转臂的机械结构运动状态。In S120, the mechanical structural motion state of the rotating arm is recorded at each set rotational position.
为了实现将旋转臂调整至设定的旋转位置处,需要调整旋转臂内部或者外部的一个或者多个机械结构进行相应的机械运动。所述机械结构运动状态可以是指在一个旋转位置处,与上述一个或者多个机械结构对应的运动距离、运动方向或者运动速度等参数。该机械结构运动状态可以通过一个向量来表示。示例性的,可以用向量si表示一个旋转位置i处的机械结构运动状态,所有设定旋转位置对应的机械结构运动状态构成矩阵S。In order to achieve adjustment of the swivel arm to the set rotational position, it is necessary to adjust one or more mechanical structures inside or outside the swivel arm for corresponding mechanical movement. The mechanical structure motion state may refer to a motion distance, a motion direction, or a motion speed corresponding to the one or more mechanical structures at a rotational position. The mechanical structure motion state can be represented by a vector. Illustratively, the mechanical structure motion state at one rotational position i can be represented by a vector s i , and the mechanical structural motion states corresponding to all the set rotational positions constitute a matrix S.
在S130中,在每个设定旋转位置处,根据所述标准体模的投影结果,对所述误差参数进行校准,得到误差参数的校准结果。In S130, at each set rotation position, the error parameter is calibrated according to the projection result of the standard phantom, and a calibration result of the error parameter is obtained.
可以理解的是,如果在旋转臂的旋转过程中,不发生任何机械抖动以及偏移,标准体模的投影结果应该有一个理想预期值。然而,在实际旋转扫描过程中,由于机械抖动或者机械偏移等因素,导致实际投影结果与该理想预期值之间存在误差。在本实施例中,可以通过将实际旋转扫描过程中获取的标准体模在的实际投影结果与这个理想预期值进行比较,确定上述旋转过程中的误差参数,其中,所述误差参数可以是实际投影结果与理想预期值之间的差值,例如:在一个设定旋转角度下,X光机旋转臂的旋转中心的理论预期值为F,通过实际投影结果计算得到的实际投影中心为H,则可以将(F-H)作为误差参数;所述误差参数还可以是一个与理想预期值对应的偏差值,例如,设定X光机旋转臂的理想旋转角度为30°,由实际投影结果确定的旋转臂的实际旋转角度为31°,则可以将31°作为误差参数。It can be understood that if no mechanical jitter and offset occur during the rotation of the rotating arm, the projection result of the standard phantom should have a desired expected value. However, during the actual rotational scanning process, there is an error between the actual projection result and the ideal expected value due to factors such as mechanical jitter or mechanical offset. In this embodiment, the error parameter in the rotation process may be determined by comparing the actual projection result of the standard phantom obtained during the actual rotation scanning process with the ideal expected value, wherein the error parameter may be actual. The difference between the projection result and the ideal expected value, for example, the theoretical expected value of the center of rotation of the X-ray machine's rotating arm is F at a set rotation angle, and the actual projection center calculated by the actual projection result is H, Then (FH) can be used as an error parameter; the error parameter can also be a deviation value corresponding to the ideal expected value, for example, setting the ideal rotation angle of the X-ray machine rotating arm to 30°, determined by the actual projection result. If the actual rotation angle of the rotating arm is 31°, 31° can be used as an error parameter.
相应的,可以根据实际投影结果与理想预期值之间的差异,对所述旋转臂在不同机械结构运动状态下的误差参数进行校准。进而可以实现当实际的扫描对象放置于X光机后,根据X光机的机械结构运动状态,以及预先校准的误差参数,对该扫描对象的投影图像进行修正。Correspondingly, the error parameters of the rotating arm in different mechanical structure motion states can be calibrated according to the difference between the actual projection result and the ideal expected value. Further, after the actual scanning object is placed on the X-ray machine, the projected image of the scanned object is corrected according to the mechanical structure motion state of the X-ray machine and the pre-calibrated error parameter.
在S140中,根据在每个设定旋转位置处的旋转臂的机械结构运动状态和所述误差参数的校准结果,生成用于对所述X光机的实际投影结果进行修正的动态修正矩阵。In S140, a dynamic correction matrix for correcting the actual projection result of the X-ray machine is generated based on the mechanical structure motion state of the rotating arm at each set rotational position and the calibration result of the error parameter.
可以理解的是,对误差参数进行校准时,所设定的旋转臂的旋转角度是有限的。因此,不可能通过穷举所有的旋转臂的机械结构运动状态来计算不同机械结构状态下对应的误差参数。但是在实际应用中,X光机的旋转臂则可能旋转至任意角度,该角度对应的机械结构运动状态的误差参数可能是未被预先校准的。此时,无法根据预先校准结果对X光机的旋转臂进行修正。为了解决上述 问题,可以基于已经校准的多个机械结构运动状态的误差参数构造一个最优估计矩阵,即:动态修正矩阵,以实现通过有限个已校准的机械结构运动状态对任意一个机械结构运动状态的误差参数进行估计。It can be understood that when the error parameter is calibrated, the set rotation angle of the rotating arm is limited. Therefore, it is impossible to calculate the corresponding error parameters in different mechanical structure states by exhausting the mechanical motion state of all the rotating arms. However, in practical applications, the rotating arm of the X-ray machine may be rotated to any angle, and the error parameter of the mechanical structure motion state corresponding to the angle may not be pre-calibrated. At this time, the X-ray machine's rotating arm cannot be corrected based on the pre-calibration result. In order to solve the above The problem is that an optimal estimation matrix can be constructed based on the error parameters of the motion state of the plurality of mechanical structures that have been calibrated, that is, the dynamic correction matrix is used to realize the error of the motion state of any one of the mechanical structures through the finite number of calibrated mechanical structure motion states. The parameters are estimated.
本实施例提供的X光机的参数校准方法,通过控制X光机旋转臂围绕标准体模进行旋转扫描,根据标准体模在多个旋转位置处的投影结果对旋转臂的不同机械结构运动状态下的误差参数进行校准,可以解决传统的对系统结构参数进行校准的方法通常是根据静态状态下系统参数与设计目标之间的差异而进行的校准,不能对系统运动中产生的误差进行校准和修正的问题,可以实现对X光机运动中产生的误差参数进行校准,为X光机的参数校准提供一种新的思路,通过利用误差参数校准结果对实际投影结果进行校准,还可以提高三维重建后图像的质量。The parameter calibration method of the X-ray machine provided by the embodiment is controlled by rotating the rotating arm of the X-ray machine around the standard phantom, and the movement state of the different mechanical structures of the rotating arm according to the projection result of the standard phantom at the plurality of rotating positions The calibration of the error parameters can solve the traditional calibration of the system structural parameters. The calibration method is usually based on the difference between the system parameters and the design objectives in the static state. The errors generated in the system motion cannot be calibrated and The revised problem can calibrate the error parameters generated in the X-ray machine motion, and provide a new idea for the parameter calibration of the X-ray machine. By using the error parameter calibration result to calibrate the actual projection result, the three-dimensional can also be improved. The quality of the image after reconstruction.
实施例二
图2a是本实施例二提供的一种X光机的参数校准方法的流程图,在上述实施例的基础上,将标准体模选为球形体模,且所述标准体模沿直径方向上可以包括两个正交对准通孔,相应的,在控制X光机的旋转臂围绕标准体模进行旋转扫描,获取多个设定旋转位置处的标准体模的投影结果之前,还可以包括:确定所述X光机的X光光源在探测器中的投影位置与所述探测器的中心重叠,所述标准体模的球心与所述X光机旋转臂的设计旋转中心重叠;可选地,所述误差参数可以包括:旋转中心在旋转平面内的切向位移误差以及径向位移误差,并基于该误差参数的类型,给出了误差参数的计算过程。相应的,本实施例的方法可以包括如下步骤。2a is a flowchart of a method for calibrating a parameter of an X-ray machine according to the second embodiment. On the basis of the above embodiment, the standard body is selected as a spherical phantom, and the standard phantom is diametrically oriented. Two orthogonal alignment through holes may be included, and correspondingly, before the rotation of the rotating arm of the X-ray machine is controlled around the standard phantom to obtain a projection result of the standard phantom at a plurality of set rotation positions, Determining that a projection position of the X-ray source of the X-ray machine in the detector overlaps with a center of the detector, and a center of the standard phantom overlaps with a designed rotation center of the X-ray machine rotating arm; Optionally, the error parameter may include: a tangential displacement error of the rotation center in the rotation plane and a radial displacement error, and based on the type of the error parameter, a calculation process of the error parameter is given. Correspondingly, the method of this embodiment may include the following steps.
在S210中,确定X光机的X光光源在探测器中的投影位置与所述探测器的中心重叠,所述标准体模的球心与所述X光机旋转臂的设计旋转中心重叠。In S210, it is determined that the projection position of the X-ray source of the X-ray machine in the detector overlaps with the center of the detector, and the center of the standard phantom overlaps with the designed rotation center of the X-ray machine rotating arm.
在本实施例中,所述标准体模为球形体模,调整该标准体模的球心,使得该标准体模的球心与所述X光机旋转臂的设计旋转中心重叠,所述标准体模沿直径方向上可以包括两个正交对准通孔。In this embodiment, the standard phantom is a spherical phantom, and the center of the standard phantom is adjusted such that the center of the standard phantom overlaps with the designed rotation center of the X-ray machine rotating arm, the standard The phantom may include two orthogonally aligned through holes in the diametrical direction.
相应的,通过在标准体模沿直径方向上设置两个正交对准通孔,可以方便将标准体模的球心设置在所述X光机旋转臂的设计旋转中心上。Correspondingly, by arranging two orthogonally aligned through holes in the diameter direction of the standard phantom, it is convenient to set the center of the standard phantom on the design rotation center of the X-ray machine rotating arm.
将标准体模置于C型臂成像范围内,调整该标准体模沿一个直径方向的位置,今X光源发射的激光穿过该标准体模直径方向上的一个通孔。此后,将旋转臂旋转π/2,调整标准体模沿另一个直径方向的位置,今激光穿过该标准体 模直径方向上的另一个通孔,其中,这两个直径方向上的通孔为正交对准的通孔。此时,可以确定该标准体模的球心位于旋转臂的设计旋转中心上。The standard phantom is placed in the C-arm imaging range, and the position of the standard phantom along a diameter direction is adjusted. The laser light emitted by the X-ray source passes through a through hole in the diameter direction of the standard phantom. Thereafter, the rotating arm is rotated by π/2 to adjust the position of the standard phantom along the other diameter direction, and the laser passes through the standard body. Another through hole in the diameter direction of the mold, wherein the two through holes in the diameter direction are orthogonally aligned through holes. At this time, it can be determined that the center of the standard phantom is located at the design rotation center of the rotating arm.
在S220中,控制X光机旋转臂围绕所述标准体模进行旋转扫描,采集多个设定旋转位置处标准体模的投影结果。In S220, the X-ray machine rotating arm is controlled to perform a rotational scan around the standard phantom, and a plurality of projection results of the standard phantom at the set rotational position are acquired.
在S230中,在每个设定旋转位置处,记录所述旋转臂的机械结构运动状态。In S230, the mechanical structural motion state of the rotating arm is recorded at each set rotational position.
在S240中,获取与设定旋转位置i对应的目标投影结果和目标机械结构运动状态,i∈[1,N],N为所述旋转臂一个旋转扫描周期经过的设定旋转角度总数。In S240, the target projection result corresponding to the set rotation position i and the target mechanical structure motion state are acquired, i ∈ [1, N], N is the total number of set rotation angles through which one rotation scan period of the rotary arm passes.
旋转臂一个旋转扫描周期可以设定为2π。在一个旋转扫描周期中,设定旋转角度可以为多个,从而可以保证设定旋转角度的全面性,获得比较丰富的目标投影结果。A rotating scan period of the rotating arm can be set to 2π. In a rotating scanning cycle, the set rotation angle can be multiple, so that the comprehensiveness of the set rotation angle can be ensured, and a relatively rich target projection result can be obtained.
在S250中,在所述目标投影结果中,获取所述标准体模的投影直径,以及所述标准体模的投影中心与所述目标投影结果的图像中心的位置差。In S250, in the target projection result, a projection diameter of the standard phantom, and a position difference between a projection center of the standard phantom and an image center of the target projection result are acquired.
标准体模是球型体模,所以经过X光机的旋转臂旋转扫描后,标准体模的直径对应有投影直径。另外,由于X光机旋转臂在旋转的过程中存在抖动的问题,使得标准体模的投影中心与目标投影结果的图像中心存在位置差。The standard phantom is a spherical phantom, so after the rotation of the rotating arm of the X-ray machine, the diameter of the standard phantom corresponds to the projected diameter. In addition, due to the problem of jitter in the rotating arm of the X-ray machine, the projection center of the standard phantom has a position difference with the image center of the target projection result.
在S260中,根据所述标准体模的投影直径、所述标准体模的实际直径和旋转臂的源板距离,计算旋转中心在旋转平面切向方向上的目标切向位移误差。In S260, the target tangential displacement error of the rotation center in the tangential direction of the rotation plane is calculated according to the projection diameter of the standard phantom, the actual diameter of the standard phantom, and the source plate distance of the rotating arm.
图2b为计算目标切向位移误差多种参数的示意图,其中,标准体模的投影直径为d1,标准体模的实际直径为d2,以及旋转臂的X光源中心A至平板探测器之间的距离(以下简称旋转臂的源板距离)为L,在目标机械结构运动状态下,旋转臂的旋转中心C在旋转平面切向方向上的投影点D的目标切向位移误差为yi。图2b中,A点为X光源中心,B为标准体模的球心,C为旋转臂的旋转中心,D为旋转中心在旋转平面切向方向上的投影点,E点为DC连线与X光源的投影射线的一个交点,F为标准体模在平板探测器上的投影中心,G为与标准体模相切的最外侧射线和旋转臂的平板探测器的交点。由于旋转臂的旋转中心C偏离标准体模的球心B的距离很小,即yi的值很小,因此,DE的长度近似为d2/2。Figure 2b is a schematic diagram of various parameters for calculating the target tangential displacement error, wherein the standard phantom has a projection diameter d 1 , the standard phantom has an actual diameter d 2 , and the rotating arm X-ray center A to the flat panel detector The distance between the distances (hereinafter referred to as the source plate distance of the rotating arm) is L. In the moving state of the target mechanical structure, the target tangential displacement error of the projection point D of the rotating center C of the rotating arm in the tangential direction of the rotating plane is y i . In Fig. 2b, point A is the center of the X-ray source, B is the center of the standard phantom, C is the center of rotation of the rotating arm, D is the projection point of the center of rotation in the tangential direction of the plane of rotation, and point E is the DC connection and An intersection of the projection ray of the X-ray source, F is the projection center of the standard phantom on the flat panel detector, and G is the intersection of the outermost ray tangential to the standard phantom and the flat panel detector of the rotating arm. Since the distance from the center of rotation C of the rotating arm to the center B of the standard phantom is small, that is, the value of y i is small, the length of DE is approximately d 2 /2.
今所述旋转臂旋转中心C至X光源中心A的距离与旋转臂源板距离L之间的比值为k,则将标准体模的球心B放置于所述旋转臂旋转中心C上后,标准体模的球心B至X光源中心A的距离,即AB的长度,与旋转臂的源板距离L的比值为k,则AD的长度为kL+yi。根据三角形的相似原理,ΔAFG~ΔADE,可以得到即因此 When the ratio of the distance between the rotation center C of the rotating arm C to the center A of the X light source and the distance L of the rotating arm source plate is k, the center B of the standard phantom is placed on the rotating center C of the rotating arm. The distance from the center B of the standard phantom to the center A of the X source, that is, the length of AB, and the distance L from the source plate of the rotating arm is k, and the length of AD is kL+y i . According to the similarity principle of triangles, ΔAFG~ΔADE can be obtained. which is therefore
在S270中,根据所述标准体模的投影直径、所述标准体模的实际直径和所述标准体模的投影中心与所述目标投影结果的图像中心的位置差,计算旋转中心在旋转平面径向方向上的目标径向位移误差。In S270, calculating a rotation center in a rotation plane according to a projection diameter of the standard phantom, an actual diameter of the standard phantom, and a position difference between a projection center of the standard phantom and an image center of the target projection result Target radial displacement error in the radial direction.
图2c为计算目标径向位移误差的多个参数的示意图,其中,标准体模的投影直径为d1,标准体模的实际直径为d2,以及旋转臂的源板距离为L,目标机械结构在运动状态下,旋转中心在旋转平面切向方向上的目标径向位移误差为xi,图2c中,A点为X光源中心,B为标准体模的球心,C为旋转中心,D1为旋转中心在旋转平面径向方向上的投影点,D为旋转中心在旋转平面切向方向上的投影点,H点为AC连线的延长线与旋转臂的平板探测器的交点,即目标投影结果的中心,F为标准体模的投影中心,G为与标准体模相切的最外侧射线与平板探测器的交点。FH的长度为标准体模投影中心与目标投影结果中心的位置差X。根据三角形的相似原理,ΔAFH~ΔADC,可以得到即结合因此 2c is a schematic diagram of a plurality of parameters for calculating a target radial displacement error, wherein the standard phantom has a projection diameter of d 1 , the standard phantom has an actual diameter of d 2 , and the rotating arm has a source plate distance of L, the target machine In the motion state, the target radial displacement error of the center of rotation in the tangential direction of the plane of rotation is x i . In Figure 2c, point A is the center of the X-ray source, B is the center of the standard phantom, and C is the center of rotation. D1 is the projection point of the center of rotation in the radial direction of the plane of rotation, D is the projection point of the center of rotation in the tangential direction of the plane of rotation, and point H is the intersection of the extension line of the AC line and the flat panel detector of the rotating arm, ie At the center of the target projection result, F is the projection center of the standard phantom, and G is the intersection of the outermost ray tangential to the standard phantom and the flat panel detector. The length of the FH is the position difference X between the standard phantom projection center and the target projection result center. According to the similarity principle of the triangle, ΔAFH~ΔADC can be obtained. which is Combine therefore
在S280中,根据所述目标切向位移误差及目标径向位移误差的校准结果,以及所述旋转臂的机械结构运动状态,生成用于对所述X光机的实际投影结果进行修正的动态修正矩阵。In S280, according to the calibration result of the target tangential displacement error and the target radial displacement error, and the mechanical structural motion state of the rotating arm, generating a dynamic for correcting the actual projection result of the X-ray machine Correct the matrix.
根据目标切向位移误差yi及目标径向位移误差xi,构成误差矩阵P,其中误差矩阵P是由pi组成的矩阵,pi=(xi,yi)。将所有设定旋转位置处的机械结构运动状态si组成机械结构矩阵S。对误差矩阵P进行奇异值分解,得到P=EW,W为权值矩阵,E为特征值矩阵。构造E(WST(SST)-1)作为动态修正矩阵。The target displacement error y i tangential and radial displacement of the target error x i, constitutes the error matrix P, where P is the error matrix is a matrix composed of p i, p i = (x i, y i). The mechanical structure motion state s i at all set rotational positions is made up of a mechanical structure matrix S. Singular value decomposition is performed on the error matrix P to obtain P=EW, W is a weight matrix, and E is a eigenvalue matrix. Construct E(WS T (SS T ) -1 ) as the dynamic correction matrix.
本实施例的技术方案通过校准旋转臂在不同机械结构运动状态下旋转中心的切向位移误差和径向位移误差,获得对X光机的实际投影结果的旋转中心进行修正的动态修正矩阵,通过应用该动态修正矩阵,可以实现对实际投影结果旋转中心的切向位移误差和径向位移误差进行修正的目的,进而可以提高三维成像后图像的质量。The technical solution of the embodiment obtains a dynamic correction matrix for correcting the rotation center of the actual projection result of the X-ray machine by calibrating the tangential displacement error and the radial displacement error of the rotation center of the rotating arm under different mechanical structural motion states, By applying the dynamic correction matrix, the tangential displacement error and the radial displacement error of the actual rotation center of the projection result can be corrected, and the image quality after three-dimensional imaging can be improved.
实施例三Embodiment 3
图3a是本实施例三提供的一种X光机的参数校准方法的流程图,本实施例以上述实施例为基础,在所述标准体模上设置与任意一个对准通孔相交的两个标记通孔,该两个标记通孔关于标准体模的球心对称,且分别垂直于两个对准通孔所确定的平面,相应的,该X光机的参数校准方法还可以包括:在每个设 定旋转位置处,记录所述旋转臂的转动角度;误差参数可为包括在所述旋转臂的不同机械结构运动状态下的采集角度,并基于该误差参数的类型,给出了误差参数的计算过程。相应的,本实施例的方法可以包括如下步骤。FIG. 3a is a flowchart of a method for calibrating a parameter of an X-ray machine according to Embodiment 3. This embodiment is based on the above embodiment, and two standard intersecting with any one of the aligned through holes are disposed on the standard phantom. The marking through holes are symmetrical with respect to the center of the standard phantom and perpendicular to the plane defined by the two aligned through holes. Correspondingly, the parameter calibration method of the X-ray machine may further include: In each setting At a fixed rotational position, the rotational angle of the rotating arm is recorded; the error parameter may be an acquisition angle included in a different mechanical structural motion state of the rotating arm, and the error parameter is calculated based on the type of the error parameter. process. Correspondingly, the method of this embodiment may include the following steps.
在S310中,确定X光机的X光光源在探测器中的投影位置与所述探测器的中心重叠,所述标准体模的球心与所述X光机旋转臂的设计旋转中心重叠。In S310, it is determined that the projection position of the X-ray source of the X-ray machine in the detector overlaps with the center of the detector, and the center of the standard phantom overlaps with the designed rotation center of the X-ray machine rotating arm.
在本实施例中,所述标准体模为球形体模,所述标准体模沿直径方向上包括两个正交对准通孔,可选地,在所述标准体模上设置与任意一个对准通孔相交的两个标记通孔,该两个标记通孔关于标准体模的球心对称,且分别垂直于两个对准通孔所确定的平面。In this embodiment, the standard phantom is a spherical phantom, and the standard phantom includes two orthogonally aligned through holes in a diameter direction, optionally, and is disposed on the standard phantom Two mark through holes intersecting the through holes are aligned, the two mark through holes being symmetric with respect to the center of the standard phantom and perpendicular to the plane defined by the two aligned through holes, respectively.
其中,为了后续误差参数方便计算,距球心等距离位置可以设置为距球心1/2半径处。Among them, in order to facilitate the calculation of the subsequent error parameters, the distance from the center of the sphere can be set to be 1/2 radius from the center of the sphere.
图3b示出了一种标准体模的俯视图,图3c示出了一种标准体模的左视图,图3d示出了一种标准体模的主视图,其中包含对准通孔31和标记通孔32。Figure 3b shows a top view of a standard phantom, Figure 3c shows a left side view of a standard phantom, and Figure 3d shows a front view of a standard phantom with alignment through
在S320中,控制X光机旋转臂围绕所述标准体模进行旋转扫描,获取多个设定旋转位置处标准体模的投影结果。In S320, the X-ray machine rotating arm is controlled to perform a rotational scan around the standard phantom to obtain a projection result of a plurality of standard phantoms at the set rotational position.
在S330中,在每个设定旋转位置处,记录所述旋转臂的机械结构运动状态和所述旋转臂的转动角度。In S330, at each set rotational position, the mechanical structural motion state of the rotating arm and the rotational angle of the rotating arm are recorded.
在本实施例中,所述旋转臂的转动角度可以是向X光机的中央处理器中输入的设定旋转角度,基于该设定旋转角度,X光机可以控制旋转臂进行相应的机械运动,以将旋转臂调整至该设定旋转角度下,并进而在旋转臂旋转结束后,记录相应的机械结构运动参数。In this embodiment, the rotation angle of the rotating arm may be a set rotation angle input to a central processing unit of the X-ray machine, and based on the set rotation angle, the X-ray machine can control the rotating arm to perform corresponding mechanical motion. In order to adjust the rotating arm to the set rotation angle, and then after the rotation of the rotating arm is finished, the corresponding mechanical structure motion parameters are recorded.
在S340中,获取与设定旋转位置i对应的目标投影结果、目标机械结构运动状态和所述旋转臂的目标转动角度,i∈[1,N],N为所述旋转臂一个旋转扫描周期经过的设定旋转角度总数。In S340, acquiring a target projection result corresponding to the set rotation position i, a target mechanical structure motion state, and a target rotation angle of the rotating arm, i ∈ [1, N], N is a rotation scan period of the rotating arm The total number of rotation angles set.
在S350中,在所述目标投影结果中,测量所述标准体模的投影直径,以及所述标准体模两个标记通孔之间的投影距离。In S350, in the target projection result, a projection diameter of the standard phantom and a projection distance between the two mark through holes of the standard phantom are measured.
在S360中,根据所述标准体模两个标记通孔之间的投影距离、所述目标转动角度以及所述标准体模的投影直径,计算目标机械结构运动状态下的目标采集角度。In S360, a target acquisition angle in a moving state of the target mechanical structure is calculated according to a projection distance between the two marked through holes of the standard phantom, the target rotation angle, and a projection diameter of the standard phantom.
图3e为在标准体模的投影图像中计算目标采集角度的多个参数的示意图,其中,标准体模两个标记通孔32之间的投影距离为a,采集角度值为Θ,标准体模的投影直径为d1。
3e is a schematic diagram of calculating a plurality of parameters of a target acquisition angle in a projection image of a standard phantom, wherein a projection distance between two standard through-
示例性的,当对应一个旋转臂旋转位置的旋转角度为θ,且θ=cos-1(2a/d1),则对应四个象限的旋转角度值分别为θ,θ+π/2,θ+π,θ+3π/2,该四个角度的集合作为采集角度。在一个旋转扫描周期中对应旋转角度θ的采集角度Θ为:Illustratively, when the rotation angle corresponding to the rotational position of one rotating arm is θ, and θ=cos −1 (2a/d 1 ), the rotation angle values corresponding to the four quadrants are θ, θ+π/2, θ, respectively. +π, θ + 3π/2, the set of four angles as the acquisition angle. The acquisition angle 对应 corresponding to the rotation angle θ in one rotation scan period is:
Θ={cos-1(2a/d1),(cos-1(2a/d1))+π/2),(cos-1(2a/d1))+π),(cos-1((2a/d1))+3π/2}。Θ={cos -1 (2a/d 1 ), (cos -1 (2a/d 1 ))+π/2), (cos -1 (2a/d 1 ))+π), (cos -1 ( (2a/d 1 )) +3π/2}.
将采集角度集合中与目标转动角度的差值最小的采集角度作为目标采集角度的校准结果。当目标转动角度为时,目标采集角度的校准结果为即目标采集角度的校准结果是与目标转动角度差值最小的采集角度值,由此可以保证得到的采集角度值与目标转动角度最接近,即该目标采集角度与实际转动角度相符合。The acquisition angle that minimizes the difference between the acquisition angle set and the target rotation angle is used as the calibration result of the target acquisition angle. When the target rotation angle is When the target acquisition angle is corrected, That is, the calibration result of the target acquisition angle is the acquisition angle value with the smallest difference from the target rotation angle, thereby ensuring that the obtained acquisition angle value is closest to the target rotation angle, that is, the target acquisition angle is consistent with the actual rotation angle.
在S370中,根据所述目标采集角度的校准结果,以及所述旋转臂的机械结构运动状态,生成用于对所述X光机的实际投影结果进行修正的动态修正矩阵。In S370, a dynamic correction matrix for correcting the actual projection result of the X-ray machine is generated according to the calibration result of the target acquisition angle and the mechanical structure motion state of the rotating arm.
根据目标采集角度的校准结果θi构成误差矩阵P,误差矩阵P是由θi组成的矩阵。将所有设定旋转位置处的机械结构运动状态si组成机械结构矩阵S。对误差矩阵进行奇异值分解,得到P=EW,W为权值矩阵,E为特征值矩阵。构造E(WST(SST)-1)作为动态修正矩阵。The calibration result θ i according to the target acquisition angle constitutes an error matrix P, which is a matrix composed of θ i . The mechanical structure motion state s i at all set rotational positions is made up of a mechanical structure matrix S. Perform singular value decomposition on the error matrix to obtain P=EW, W is the weight matrix, and E is the eigenvalue matrix. Construct E(WS T (SS T ) -1 ) as the dynamic correction matrix.
本实施例的技术方案通过标准体模两个校准通孔之间的投影距离,目标转动角度以及标准体模的投影直径,计算目标机械结构运动状态下的目标采集角度,根据目标采集角度的校准结果以及旋转臂的机械结构运动状态,生成用于对X光机的实际投影结果进行修正的动态修正矩阵,可以实现利用对应误差参数为采集角度时的动态修正矩阵进行实际投影结果的修正,提高图像的质量。The technical solution of the embodiment calculates the target acquisition angle under the motion state of the target mechanical structure by the projection distance between the two calibrated through holes of the standard phantom, the target rotation angle and the projection diameter of the standard phantom, and the calibration according to the target acquisition angle. As a result, as well as the mechanical structure motion state of the rotating arm, a dynamic correction matrix for correcting the actual projection result of the X-ray machine is generated, and the actual correction result can be corrected by using the dynamic correction matrix when the corresponding error parameter is used for the acquisition angle. The quality of the image.
实施例四Embodiment 4
图4是本实施例四提供的一种X光机的参数校准方法的流程图,在上述实施例基础上,根据所述误差参数的校准结果,以及所述旋转臂的机械结构运动状态,生成用于对所述X光机的实际投影结果进行修正的动态修正矩阵还可以包括:将多个设定旋转位置处的机械结构运动状态组成机械结构矩阵S;将多个机械结构运动状态下的误差参数组成误差矩阵P;4 is a flowchart of a method for calibrating a parameter of an X-ray machine according to the fourth embodiment. Based on the above embodiment, a calibration result of the error parameter and a motion state of the mechanical structure of the rotating arm are generated. The dynamic correction matrix for correcting the actual projection result of the X-ray machine may further comprise: forming a mechanical structure motion state at a plurality of set rotation positions into a mechanical structure matrix S; and moving the plurality of mechanical structures in a state of motion Error parameters constitute an error matrix P;
对误差矩阵P进行奇异值分解,得到P=EW,其中,W为权值矩阵,E为特征 值矩阵;构造E(WST(SST)-1)作为所述动态修正矩阵。相应的,本实施例的方法可以包括:Singular value decomposition is performed on the error matrix P to obtain P = EW, where W is a weight matrix, E is a feature value matrix, and E (WS T (SS T ) -1 ) is constructed as the dynamic correction matrix. Correspondingly, the method in this embodiment may include:
在S410中,控制X光机旋转臂围绕标准体模进行旋转扫描,获取多个设定旋转位置处标准体模的投影结果。In S410, the X-ray machine rotating arm is controlled to perform a rotational scan around the standard phantom to obtain projection results of a plurality of standard phantoms at the set rotational position.
在S420中,在每个设定旋转位置处,记录所述旋转臂的机械结构运动状态。In S420, the mechanical structural motion state of the rotating arm is recorded at each set rotational position.
在S430中,在每个设定旋转位置处,根据所述标准体模的投影结果,对所述误差参数进行校准,得到误差参数的校准结果。In S430, at each set rotation position, the error parameter is calibrated according to the projection result of the standard phantom, and a calibration result of the error parameter is obtained.
在S440中,将多个设定旋转位置处的机械结构运动状态组成机械结构矩阵S。In S440, the mechanical structural motion states at the plurality of set rotational positions are made up to form the mechanical structure matrix S.
在S450中,将多个机械结构运动状态下的误差参数组成误差矩阵P。In S450, error parameters in a plurality of mechanical structural motion states are formed into an error matrix P.
示例性的,误差参数可以是径向位移误差xi、切向位移误差yi和采集角度校准结果θi,误差矩阵P是由pi组成的矩阵,pi=(xi,yi,θi)。Exemplary, error parameters may be a radial displacement error x i, y i tangential displacement error calibration results and the acquisition angle θ i, the error matrix P is a matrix composed of p i, p i = (x i, y i, θ i ).
在S460中,对误差矩阵P进行奇异值分解,得到P=EW,其中,W为权值矩阵,E为特征值矩阵。In S460, the error matrix P is subjected to singular value decomposition to obtain P=EW, where W is a weight matrix and E is a feature value matrix.
在S470中,构造E(WST(SST)-1)作为动态修正矩阵,用于对所述X光机的实际投影结果进行修正。In S470, the configuration E(WS T (SS T ) -1 ) is constructed as a dynamic correction matrix for correcting the actual projection result of the X-ray machine.
S为所有机械结构的运动状态si组成的矩阵。S is a matrix of motion states s i of all mechanical structures.
本实施例的技术方案通过误差参数构造动态修正矩阵,可以实现对实际投影结果的校正,提高三维成像后图像的质量。The technical solution of the embodiment constructs a dynamic correction matrix by using error parameters, which can correct the actual projection result and improve the quality of the image after three-dimensional imaging.
实施例五Embodiment 5
图5为本实施例五提供的一种X光机的参数校准方法的流程图,在上述实施例基础上,该X光机的参数校准方法还可以增加根据校准生成的动态修成矩阵,对实际投影结果进行修正的过程。相应的,本实施例的方法可以包括如下步骤。FIG. 5 is a flowchart of a parameter calibration method for an X-ray machine according to Embodiment 5, on the basis of the foregoing embodiment, the parameter calibration method of the X-ray machine can further increase a dynamic repair matrix generated according to the calibration, and the actual The process of correcting the projection result. Correspondingly, the method of this embodiment may include the following steps.
在S510中,控制X光机旋转臂围绕标准体模进行旋转扫描,采集多个设定旋转位置处标准体模的投影结果。In S510, the X-ray machine rotating arm is controlled to perform a rotational scan around the standard phantom, and a plurality of projection results of the standard phantom at the set rotational position are acquired.
在S520中,在每个设定旋转位置处,记录所述旋转臂的机械结构运动状态。In S520, the mechanical structural motion state of the rotating arm is recorded at each set rotational position.
在S530中,在每个设定旋转位置处,根据所述标准体模的投影结果,对所述误差参数进行校准,得到误差参数的校准结果。In S530, at each set rotation position, the error parameter is calibrated according to the projection result of the standard phantom, and a calibration result of the error parameter is obtained.
在S540中,根据在每个而设定旋转位置处的旋转臂的机械结构运动状态和所述误差参数的校准结果,生成用于对所述X光机的实际投影结果进行修正的 动态修正矩阵。In S540, generating a correction result for correcting the actual projection result of the X-ray machine according to the mechanical structure motion state of the rotating arm at each of the set rotational positions and the calibration result of the error parameter Dynamic correction matrix.
在S550中,控制所述X光机旋转臂围绕被测对象进行旋转扫描,获取实际旋转位置处的待修正投影结果。In S550, the X-ray machine rotating arm is controlled to perform a rotation scan around the object to be measured, and the projection result to be corrected at the actual rotation position is obtained.
在S560中,在实际旋转位置处,记录所述旋转臂的实测机械结构运动状态。In S560, at the actual rotational position, the measured mechanical structure motion state of the rotating arm is recorded.
在S570中,根据所述实测机械结构运动状态以及所述动态修正矩阵,生成与所述实测机械结构运动状态对应的实际误差参数。In S570, an actual error parameter corresponding to the motion state of the measured mechanical structure is generated according to the measured mechanical structure motion state and the dynamic correction matrix.
在实际旋转过程中,对任意实际机械结构运动状态sj,对应误差矩阵P的最佳估计为E(WST(SST)-1)sj,由此可实现对实际投影结果的修正。In the actual rotation process, for any actual mechanical structure motion state s j , the best estimate of the corresponding error matrix P is E(WS T (SS T ) -1 ) s j , thereby realizing the correction of the actual projection result.
在S580中,根据所述实际误差参数,对所述待修正投影结果进行图像修正。In S580, image correction is performed on the projection result to be corrected according to the actual error parameter.
所述实际误差参数可以包括:在所述旋转臂的实测机械结构运动状态sj下,旋转中心在旋转平面上的实测目标切向位移误差yj以及实测目标径向位移误差xj;根据所述实际误差参数,对所述待修正投影结果进行图像修正可以包括如下步骤。The actual error parameter may include: a measured target tangential displacement error y j of the rotation center on the rotation plane and a measured target radial displacement error x j under the measured mechanical structure motion state s j of the rotating arm; The actual error parameter, and performing image correction on the projection result to be corrected may include the following steps.
获取与实际旋转位置j对应的目标待修正投影结果,目标实际机械结构运动状态,j∈[1,M],M为所述旋转臂一个旋转扫描周期经过的实际旋转角度总数;获取与目标实测机械结构运动状态对应的目标切向位移误差、目标径向位移误差。Obtaining a target to be corrected projection result corresponding to the actual rotation position j, the target actual mechanical structure motion state, j ∈ [1, M], M is the total number of actual rotation angles of a rotation scan period of the rotating arm; acquisition and target measurement The target tangential displacement error and the target radial displacement error corresponding to the motion state of the mechanical structure.
根据所述目标切向位移误差,对所述目标待修正投影结果进行设定比例的缩放处理。And according to the target tangential displacement error, the target correction target projection result is subjected to a scaling process of a set ratio.
由上述实施例可知,实测目标切向位移误差其中,标准体模的投影直径d1,标准体模的实际直径d2,以及旋转臂的源板距离L,k为标准体模的球心至光源距离,即AB的长度,与源板距离L的比值,则标准体模的实际图像与投影图像的比值令其中,理想状态下因此可以今目标待修正投影结果进行倍数为的缩放运算,以消除缩放误差。可选地,当k=1/2时, It can be known from the above embodiment that the measured target tangential displacement error Wherein, the projection diameter d 1 of the standard phantom, the actual diameter d 2 of the standard phantom, and the source plate distance L, k of the rotating arm are the distance from the center of the sphere of the standard phantom to the light source, that is, the length of the AB, and the distance from the source plate. The ratio of L, the ratio of the actual image of the standard phantom to the projected image make Among them, ideally Therefore, the projection result of the target to be corrected can be multiplied as The scaling operation to eliminate scaling errors. Optionally, when k=1/2,
根据所述目标径向位移误差,对所述目标待修正投影结果进行设定距离的平移。And shifting the set distance to the target to be corrected projection result according to the target radial displacement error.
由上述实施例可知,目标径向位移误差其中,X为标准体模投影中心与目标投影结果中心的位置差,标准体模的投影直径d1,标准体模的实际 直径d2,则标准体模投影中心与目标投影结果中心的位置差因此,通过对目标待修正投影结果进行平移,可以消除平移误差。According to the above embodiment, the target radial displacement error Where X is the position difference between the standard phantom projection center and the target projection result center, the projection diameter d 1 of the standard phantom, and the actual diameter d 2 of the standard phantom, the difference between the standard phantom projection center and the target projection result center Therefore, by performing the projection result of the target to be corrected Pan, you can eliminate the translation error.
所述实际误差参数可以包括:在所述旋转臂的实测机械结构运动状态下的实际采集角度。The actual error parameter may include an actual acquisition angle in a state of motion of the measured mechanical structure of the rotating arm.
根据所述实际误差参数,对所述待修正投影结果进行图像修正可以包括:根据不同实测机械结构运动状态下的所述实际采集角度,对每个待修正投影结果进行三维重建。And performing image correction on the to-be-corrected projection result according to the actual error parameter may include: performing three-dimensional reconstruction on each to-be-corrected projection result according to the actual acquisition angle in different measured mechanical structure motion states.
根据所述目标径向位移误差,对所述目标待修正投影结果进行设定距离的平移,根据所述目标切向位移误差,对所述目标待修正投影结果进行设定比例的缩放处理,达到对待修正投影结果的校正的目的。根据校正后的图像和对应角度θj进行三维重建,提高三维重建后的图像质量。And performing a set distance translation on the target to be corrected projection result according to the target radial displacement error, and performing scaling processing on the target to be corrected projection result according to the target tangential displacement error, to achieve The purpose of the correction of the corrected projection result. The three-dimensional reconstruction is performed according to the corrected image and the corresponding angle θ j to improve the image quality after the three-dimensional reconstruction.
本实施例的技术方案根据目标切向位移误差,对目标待修正投影结果进行设定比例的缩放处理;根据目标径向位移误差,对目标待修正投影结果进行设定距离的平移,根据不同实测机械结构运动状态下的实际采集角度,对每个待修正投影结果进行三维重建,可以实现对图像的校正,达到提高三维重建后的图像质量的效果。According to the technical solution of the embodiment, the target tangential displacement error is used to perform a scaling process on the target to be corrected projection result; according to the target radial displacement error, the target distance to be corrected is subjected to a set distance translation, according to different actual measurements. The actual acquisition angle of the mechanical structure in motion state, three-dimensional reconstruction of each projection result to be corrected, can achieve image correction, and achieve the effect of improving the image quality after three-dimensional reconstruction.
实施例六Embodiment 6
图6为本实施例六提供的一种X光机的参数校准装置的结构图。本公开实施例提供的X光机的参数校准装置可以应用于本公开实施例所述的X光机的参数校准系统中。如图6所示,所述装置包括:投影结果获取模块610、运动状态记录模块620、误差参数校准模块630和动态修正矩阵生成模块640。FIG. 6 is a structural diagram of a parameter calibration apparatus for an X-ray machine according to Embodiment 6. The parameter calibration device of the X-ray machine provided by the embodiment of the present disclosure can be applied to the parameter calibration system of the X-ray machine described in the embodiment of the present disclosure. As shown in FIG. 6, the apparatus includes a projection result acquisition module 610, a motion state recording module 620, an error parameter calibration module 630, and a dynamic correction matrix generation module 640.
投影结果获取模块610,设置为控制X光机旋转臂围绕标准体模进行旋转扫描,采集多个设定旋转位置处的投影结果。The projection result obtaining module 610 is configured to control the X-ray machine rotating arm to perform a rotation scan around the standard phantom, and collect the projection results at the plurality of set rotation positions.
运动状态记录模块620,设置为在每个设定旋转位置处,记录所述旋转臂的机械结构运动状态。The motion state recording module 620 is configured to record the mechanical structural motion state of the rotating arm at each of the set rotational positions.
误差参数校准模块630,设置为在每个设定旋转位置处,根据所述标准体模的投影结果,对所述误差参数进行校准,得到误差参数的校准结果。The error parameter calibration module 630 is configured to calibrate the error parameter according to the projection result of the standard phantom at each set rotation position to obtain a calibration result of the error parameter.
动态修正矩阵生成模块640,设置为根据在每个设定旋转位置处的旋转臂的机械结构运动状态和所述误差参数的校准结果,生成用于对所述X光机的实际投影结果进行修正的动态修正矩阵。 The dynamic correction matrix generating module 640 is configured to generate a correction for the actual projection result of the X-ray machine according to the mechanical structure motion state of the rotating arm at each set rotational position and the calibration result of the error parameter Dynamic correction matrix.
本实施例提供的X光机的参数校准装置,通过控制X光机旋转臂围绕标准体模进行旋转扫描,根据标准体模在多个旋转位置处的投影结果对旋转臂的机械结构运动状态下的误差参数进行校准,可以解决在对系统结构参数进行校准的方法通常为根据静态状态下系统参数与设计目标之间的差异而进行的校准,不能对系统运动产生的误差进行校准和修正的问题,可以实现对X光机运动中产生的误差参数进行校准,为X光机的参数校准提供一种新的思路,通过利用误差参数校准结果对实际投影结果进行修正,提高三维重建后图像的质量。在上述实施例的基础上,所述标准体模为球形体模,且所述标准体模沿直径方向上包括两个正交对准通孔。The parameter calibration device of the X-ray machine provided by the embodiment is controlled to rotate around the standard phantom by controlling the rotating arm of the X-ray machine, and the mechanical structure of the rotating arm is moved according to the projection result of the standard phantom at a plurality of rotating positions. The calibration of the error parameters can solve the problem that the calibration of the system structural parameters is usually based on the difference between the system parameters and the design objectives in the static state, and the error caused by the system motion cannot be calibrated and corrected. It can realize the calibration of the error parameters generated in the X-ray machine motion, and provide a new idea for the parameter calibration of the X-ray machine. By using the error parameter calibration result to correct the actual projection result, the image quality after 3D reconstruction is improved. . Based on the above embodiment, the standard phantom is a spherical phantom, and the standard phantom includes two orthogonally aligned through holes in the diametrical direction.
相应的,所述X光机的参数校准装置还可以包括:投影位置确定模块,设置为在控制X光机的旋转臂围绕标准体模进行旋转扫描,获取多个设定旋转位置处标准体模的投影结果之前,确定所述X光机的X光光源在探测器中的投影位置与所述探测器的中心重叠,所述标准体模的球心与所述X光机旋转臂的设计旋转中心重叠。Correspondingly, the parameter calibration device of the X-ray machine may further include: a projection position determining module configured to perform a rotational scan around the standard phantom of the rotating arm of the X-ray machine to obtain a standard phantom at a plurality of set rotational positions; Before the projection result, determining that the projection position of the X-ray source of the X-ray machine in the detector overlaps with the center of the detector, the spherical center of the standard phantom and the design rotation of the X-ray machine rotating arm The center overlaps.
可选地,所述误差参数可以包括:在所述旋转臂的每个机械结构运动状态下,旋转中心在旋转平面上的切向位移误差以及径向位移误差。Optionally, the error parameter may include a tangential displacement error of the center of rotation on the plane of rotation and a radial displacement error in each mechanical structural motion state of the rotating arm.
相应的,所述误差参数校准模块630可以设置为:获取与设定旋转位置i对应的目标投影结果和目标机械结构运动状态,i∈[1,N],N为所述旋转臂一个旋转扫描周期经过的设定旋转角度总数;在所述目标投影结果中,获取所述标准体模的投影直径,以及所述标准体模的投影中心与所述目标投影结果的图像中心的位置差;根据所述标准体模的投影直径、所述标准体模的实际直径和旋转臂的源板距离,计算旋转中心在旋转平面切向方向上的目标切向位移误差;根据所述标准体模的投影直径、所述标准体模的实际直径,以及所述标准体模的投影中心与所述目标投影结果的图像中心的位置差,计算旋转中心在旋转平面径向方向上的目标径向位移误差。Correspondingly, the error parameter calibration module 630 can be configured to: acquire a target projection result corresponding to the set rotation position i and a target mechanical structure motion state, i ∈ [1, N], N is a rotation scan of the rotating arm a total number of set rotation angles through which the cycle passes; in the target projection result, a projection diameter of the standard phantom, and a position difference between a projection center of the standard phantom and an image center of the target projection result; a projection diameter of the standard phantom, an actual diameter of the standard phantom, and a source plate distance of the rotating arm, and calculating a target tangential displacement error of the rotation center in a tangential direction of the rotation plane; projection according to the standard phantom The diameter, the actual diameter of the standard phantom, and the position difference between the projection center of the standard phantom and the image center of the target projection result, and the target radial displacement error of the rotation center in the radial direction of the rotation plane is calculated.
可选地,在所述标准体模上设置与任意一个对准通孔相交的两个标记通孔,该两个标记通孔关于标准体模的球心对称,且分别垂直于两个对准通孔所确定的平面。Optionally, two standard through holes intersecting any one of the aligned through holes are disposed on the standard phantom, the two labeled through holes being symmetric with respect to a center of the standard phantom and perpendicular to the two alignments respectively The plane defined by the through hole.
相应的,所述运动状态记录模块620还可以设置为:在每个设定旋转位置处,记录所述旋转臂的转动角度。Correspondingly, the motion state recording module 620 can also be configured to record the rotation angle of the rotating arm at each set rotation position.
可选地,所述误差参数可以包括:在所述旋转臂的不同机械结构运动状态下的采集角度。 Optionally, the error parameter may include an acquisition angle in a different mechanical structural motion state of the rotating arm.
相应的,所述误差参数校准模块630可以设置为:获取与设定旋转位置i对应的目标投影结果、目标机械结构运动状态和所述旋转臂的目标转动角度,1∈[1,N],N为所述旋转臂一个旋转扫描周期经过的设定旋转角度总数;在所述目标投影结果中,获取所述标准体模的投影直径,以及所述标准体模两个标记通孔之间的投影距离;根据所述标准体模两个标记通孔之间的投影距离、所述目标转动角度以及所述标准体模的投影直径,计算目标机械结构运动状态下的目标采集角度。Correspondingly, the error parameter calibration module 630 can be configured to: acquire a target projection result corresponding to the set rotation position i, a target mechanical structure motion state, and a target rotation angle of the rotating arm, 1 ∈ [1, N], N is a total number of set rotation angles of a rotation scan period of the rotating arm; in the target projection result, a projection diameter of the standard phantom is acquired, and between the two mark through holes of the standard phantom Projection distance; calculating a target acquisition angle in a moving state of the target mechanical structure according to a projection distance between the two marked through holes of the standard phantom, the target rotation angle, and a projection diameter of the standard phantom.
可选地,所述动态修正矩阵生成模块640可以设置为:将多个设定旋转位置处的机械结构运动状态组成机械结构矩阵S;将多个机械结构运动状态下的误差参数组成误差矩阵P;对误差矩阵P进行奇异值分解,得到P=EW,其中,W为权值矩阵,E为特征值矩阵;构造E(WST(SST)-1)作为所述动态修正矩阵。Optionally, the dynamic correction matrix generating module 640 may be configured to: compose a mechanical structure motion state at a plurality of set rotation positions into a mechanical structure matrix S; and combine error parameters in a plurality of mechanical structure motion states into an error matrix P The singular value decomposition is performed on the error matrix P to obtain P=EW, where W is a weight matrix, E is an eigenvalue matrix, and the structure E(WS T (SS T ) -1 ) is used as the dynamic correction matrix.
可选地,所述X光机的参数校准装置还可以包括:待修正投影结果获取模块,设置为在根据在每个设定位置处的旋转臂的机械机构运动状态和所述误差参数的校准结果,生成用于对所述X光机的实际投影结果进行修正的动态修正矩阵之后,控制所述X光机旋转臂围绕被测对象进行旋转扫描,采集实际旋转位置处的待修正投影结果。Optionally, the parameter calibration device of the X-ray machine may further include: a projection result obtaining module to be corrected, configured to be calibrated according to a mechanical mechanism motion state of the rotating arm at each set position and the error parameter As a result, after generating a dynamic correction matrix for correcting the actual projection result of the X-ray machine, the X-ray machine rotating arm is controlled to perform a rotation scan around the object to be measured, and the projection result to be corrected at the actual rotation position is acquired.
实测机械结构运动状态记录模块,设置为在根据在每个设定位置处的旋转臂的机械机构运动状态和所述误差参数的校准结果,生成用于对所述X光机的实际投影结果进行修正的动态修正矩阵之后,在实际旋转位置处,记录所述旋转臂的实测机械结构运动状态。The measured mechanical structure motion state recording module is configured to generate an actual projection result for the X-ray machine based on a calibration result of a mechanical mechanism motion state of the rotating arm at each set position and the error parameter After the modified dynamic correction matrix, the measured mechanical motion state of the rotating arm is recorded at the actual rotational position.
实际误差参数生产模块,设置为根据所述实测机械结构运动状态以及所述动态修正矩阵,生成与所述实测机械结构运动状态对应的实际误差参数。The actual error parameter production module is configured to generate an actual error parameter corresponding to the motion state of the measured mechanical structure according to the measured mechanical structure motion state and the dynamic correction matrix.
待修正投影结果修正模块,设置为根据所述实际误差参数,对所述待修正投影结果进行图像修正。The projection result correction module to be corrected is configured to perform image correction on the projection result to be corrected according to the actual error parameter.
可选地,所述实际误差参数可以包括:在实际旋转位置处,旋转臂的旋转中心在旋转平面内的实际切向位移误差以及实际径向位移误差;Optionally, the actual error parameter may include: an actual tangential displacement error of the center of rotation of the rotating arm in the plane of rotation and an actual radial displacement error at the actual rotational position;
所述实际误差参数生产模块可以设置为:获取与实际旋转位置j对应的目标待修正投影结果核目标实际机械结构运动状态,j∈[1,M],M为所述旋转臂一个旋转扫描周期经过的实际旋转角度总数;根据所述动态修正矩阵,获取与目标实测机械结构运动状态对应的目标切向位移误差、目标径向位移误差;根据所述目标切向位移误差,对所述目标待修正投影结果进行设定比例的缩放处理; 根据所述目标径向位移误差,对所述目标待修正投影结果进行设定距离的平移。The actual error parameter production module may be configured to: acquire a target to be corrected projection result corresponding to the actual rotation position j, the actual mechanical structure motion state of the nuclear target, j ∈ [1, M], M is a rotation scan period of the rotating arm The total number of actual rotation angles passed; according to the dynamic correction matrix, acquiring a target tangential displacement error and a target radial displacement error corresponding to the moving state of the target measured mechanical structure; and waiting for the target according to the target tangential displacement error Correcting the projection result to perform scaling processing of the set ratio; And shifting the set distance to the target to be corrected projection result according to the target radial displacement error.
可选地,所述实际误差参数可以包括:在所述旋转臂的实测机械结构运动状态下的实际采集角度。Optionally, the actual error parameter may include an actual acquisition angle in a state of motion of the measured mechanical structure of the rotating arm.
所述实际误差参数生产模块还可以设置为:根据不同实测机械结构运动状态下的所述实际采集角度,对每个待修正投影结果进行三维重建。The actual error parameter production module may be further configured to perform three-dimensional reconstruction on each of the projections to be corrected according to the actual acquisition angles in different motion states of the mechanical structure.
本实施例所提供的X光机的参数校准装置可用于执行上述任意实施例提供的X光机的参数校准方法,具备相应的功能模块。The parameter calibration device of the X-ray machine provided in this embodiment can be used to perform the parameter calibration method of the X-ray machine provided by any of the above embodiments, and has a corresponding function module.
实施例七Example 7
图7为本实施例七提供的一种X光机的参数校准系统的结构图。如图7所示,所述系统可以包括:旋转臂710、相对设置于所述旋转臂之上的X光光源720以及探测器730、处理器740、图像采集器件750、运动状态采集部件760以及标准体模770。FIG. 7 is a structural diagram of a parameter calibration system for an X-ray machine according to Embodiment 7. As shown in FIG. 7, the system may include: a
所述旋转臂710,设置为根据所述处理器发送的控制信号带动所述X光光源以及平板探测器围绕所述标准体模进行多个设定旋转角度的旋转。The
所述X光光源720,设置为向所述标准体模发送X光扫描信号。The
所述探测器730,设置为采集所述X光扫描信号在所述标准体模上的投影数据。The
所述图像采集器件750,设置为根据所述探测器采集的投影数据,生成标准体模的投影结果。The
所述运动状态采集部件760,设置为获取所述旋转臂的机械结构运动状态。The motion
所述处理器740,设置为控制所述旋转臂围绕标准体模进行旋转扫描,获取设定多个旋转位置处标准体模的投影结果;在每个设定旋转位置处,记录所述旋转臂的机械结构运动状态;在每个设定旋转位置处,根据所述标准体模的投影结果,对所述误差参数进行校对,得到误差参数的校对结果;根据在每个设定旋转位置处的旋转臂的机械结构运动状态和所述误差参数的校准结果,生成用于对所述X光机的实际投影结果进行修正的动态修正矩阵。The
可选的,所述标准体模770为球形体模,且所述标准体模770沿直径方向上包括两个正交对准通孔。Optionally, the
所述处理器740还可以设置为:在控制X光机的旋转臂围绕标准体模770进行旋转扫描,获取不同旋转位置处标准体模的投影结果之前,确定所述X光机的X光光源720在探测器730中的投影位置与所述探测器730的中心重叠,
所述标准体模770的球心与所述X光机旋转臂的设计旋转中心重叠。The
可选的,在所述标准体模770上设置与任意一个对准通孔相交的两个标记通孔,该两个标记通孔关于标准体模的球心对称,且分别垂直于两个对准通孔所确定的平面。Optionally, two standard through holes intersecting any one of the aligned through holes are disposed on the
所述处理器740还可以设置为:在每个设定旋转位置处,记录所述旋转臂的转动角度。The
图8为一种X光机的参数校准系统对应的校准方法的流程图,可以包括:8 is a flow chart of a calibration method corresponding to a parameter calibration system of an X-ray machine, which may include:
在S810中,初始对准。In S810, initial alignment.
初始对准可以是指对准标准体模和X光光源中心与平板探测器的中心。Initial alignment may refer to aligning the center of the standard phantom and X-ray source with the center of the flat panel detector.
在S820中,曝光成像。In S820, the exposure is imaged.
曝光成像可以是指对标准体模进行旋转扫描。Exposure imaging can refer to a rotational scan of a standard phantom.
在S830中,图像分析与校准。In S830, image analysis and calibration.
对扫描获得的投影结果进行图像分析,计算误差参数。Perform image analysis on the projection results obtained by scanning and calculate error parameters.
在S840中,被测物体的实际投影结果校正。In S840, the actual projection result of the measured object is corrected.
根据误差参数对被测物体的实际投影结果进行修正,得到更加符合实际的高质量图像。According to the error parameter, the actual projection result of the measured object is corrected to obtain a more realistic high quality image.
相应的,图9为一种X光机的参数校准系统的结构图,该X光机的参数校准系统可以包括:初始对准模块910,设置为标准体模的初始对准;图像采集模块920,设置为采集待校正被测物体实际的投影结果;系统校准模块930,设置为校准系统的误差参数;图像校正模块940,设置为根据误差参数对待校正的被测物体的实际投影结果进行修正;人机交互模块950,设置为人机交互。Correspondingly, FIG. 9 is a structural diagram of a parameter calibration system of an X-ray machine. The parameter calibration system of the X-ray machine may include: an
本实施例提供的一种X光机的参数校准系统,通过控制X光机旋转臂围绕标准体模进行旋转扫描,根据标准体模在多个旋转位置处的的投影结果对旋转臂的不同机械结构运动状态下的误差参数进行校准,可以解决在对系统结构参数进行校准的方法通常是根据静态状态下系统参数与设计目标之间的差异而进行的校准,不能对系统运动产生的误差进行校准和修正的问题,可以实现对X光机运动中产生的误差参数进行校准,为X光机的参数校准提供一种新的思路,通过利用误差参数校准结果对实际投影结果进行校准,可以提高三维重建后图像的质量。The parameter calibration system of the X-ray machine provided by the embodiment is controlled by rotating the rotating arm of the X-ray machine around the standard phantom, and the different mechanical parts of the rotating arm according to the projection result of the standard phantom at the plurality of rotating positions Calibration of the error parameters in the structural motion state can be solved. The method of calibrating the system structural parameters is usually calibrated according to the difference between the system parameters and the design target in the static state, and the error caused by the system motion cannot be calibrated. And the correction problem can calibrate the error parameters generated in the X-ray machine motion, and provide a new idea for the parameter calibration of the X-ray machine. By using the error parameter calibration result to calibrate the actual projection result, the three-dimensional can be improved. The quality of the image after reconstruction.
本公开提供了一种X光机的参数校准方法、装置及系统,可以实现对X光 机运动中产生的误差参数进行校准,通过利用误差参数校准结果对实际投影结果进行校正,可以提高三维重建后图像的质量。 The present disclosure provides a parameter calibration method, device and system for an X-ray machine, which can realize X-ray The error parameters generated during the motion of the machine are calibrated. By correcting the actual projection result by using the error parameter calibration result, the quality of the image after three-dimensional reconstruction can be improved.
Claims (13)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/CN2017/072177 WO2018133090A1 (en) | 2017-01-23 | 2017-01-23 | Parameter calibration method, device, and system for x-ray machine |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/CN2017/072177 WO2018133090A1 (en) | 2017-01-23 | 2017-01-23 | Parameter calibration method, device, and system for x-ray machine |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2018133090A1 true WO2018133090A1 (en) | 2018-07-26 |
Family
ID=62907666
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/CN2017/072177 Ceased WO2018133090A1 (en) | 2017-01-23 | 2017-01-23 | Parameter calibration method, device, and system for x-ray machine |
Country Status (1)
| Country | Link |
|---|---|
| WO (1) | WO2018133090A1 (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102144927A (en) * | 2010-02-10 | 2011-08-10 | 清华大学 | Motion-compensation-based computed tomography (CT) equipment and method |
| CN102652674A (en) * | 2011-03-04 | 2012-09-05 | 首都师范大学 | Method and system for eliminating geometrical artifacts in CT (Computerized Tomography) image |
| CN102692421A (en) * | 2012-06-18 | 2012-09-26 | 东营市三英精密工程研究中心 | High-precision x-ray microscope sample scanning table with metering rotary shaft |
| CN103759679A (en) * | 2013-11-04 | 2014-04-30 | 杭州长庚测量技术有限公司 | Measuring method for angular deviation of cone-beam CT system |
| CN204520733U (en) * | 2015-02-16 | 2015-08-05 | 清华大学 | A kind of die body and use the CBCT system of this die body |
| WO2016055167A1 (en) * | 2014-10-10 | 2016-04-14 | Werth Messtechnik Gmbh | Method and device for the computer tomography of a workpiece |
| CN105849537A (en) * | 2013-11-28 | 2016-08-10 | 尼康计量公众有限公司 | Calibration apparatus and method for computed tomography |
-
2017
- 2017-01-23 WO PCT/CN2017/072177 patent/WO2018133090A1/en not_active Ceased
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102144927A (en) * | 2010-02-10 | 2011-08-10 | 清华大学 | Motion-compensation-based computed tomography (CT) equipment and method |
| CN102652674A (en) * | 2011-03-04 | 2012-09-05 | 首都师范大学 | Method and system for eliminating geometrical artifacts in CT (Computerized Tomography) image |
| CN102692421A (en) * | 2012-06-18 | 2012-09-26 | 东营市三英精密工程研究中心 | High-precision x-ray microscope sample scanning table with metering rotary shaft |
| CN103759679A (en) * | 2013-11-04 | 2014-04-30 | 杭州长庚测量技术有限公司 | Measuring method for angular deviation of cone-beam CT system |
| CN105849537A (en) * | 2013-11-28 | 2016-08-10 | 尼康计量公众有限公司 | Calibration apparatus and method for computed tomography |
| WO2016055167A1 (en) * | 2014-10-10 | 2016-04-14 | Werth Messtechnik Gmbh | Method and device for the computer tomography of a workpiece |
| CN204520733U (en) * | 2015-02-16 | 2015-08-05 | 清华大学 | A kind of die body and use the CBCT system of this die body |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106821405B (en) | parameter calibration method, device and system of X-ray machine | |
| CN102743184B (en) | Geometrical parameter calibration method of X-ray cone beam computed tomography system | |
| CN108122203B (en) | Geometric parameter correction method, device, equipment and system | |
| CN102652674B (en) | A method and system for removing geometric artifacts in CT images | |
| CN103800032B (en) | The corrective system corrected for cone-beam CT system geometric position and bearing calibration thereof | |
| CN108201447A (en) | A kind of static state CT system geometric parameter bearing calibration | |
| CN108430334B (en) | Ultrasound imaging apparatus and ultrasound imaging method for examining volume of object | |
| CN103784160B (en) | A kind of correcting unit of cone-beam CT system geometric position and bearing calibration thereof | |
| CN100581479C (en) | Reconstruction Method of 3D Ultrasound Image | |
| CN107684435A (en) | Cone-beam CT system geometric calibration method and its calibrating installation | |
| CN103735282A (en) | Device and method for geometric correction of detector of cone-beam CT (computed tomography) system | |
| CN104665862B (en) | Method for eliminating geometrical artifacts from CBCT (Cone Beam Computed Tomography) and CBCT system using method for eliminating geometrical artifacts from CBCT | |
| CN202049120U (en) | System for eliminating geometric artifacts in CT (computed tomography) image | |
| CN104764756A (en) | Calibration method for cone beam CT imaging system | |
| JP2019158534A (en) | X-ray ct apparatus for measurement and method for generating fault image | |
| TW561241B (en) | Method and apparatus for calibrating laser three-dimensional digitizing sensor | |
| CN203776924U (en) | Calibration device for geometric position of cone-beam CT (computed tomography) system | |
| CN102253061B (en) | Vertical cone beam CT (Computed Tomography) imaging calibration system and method applying same | |
| CN106725556A (en) | X-ray machine, X-ray machine pivot method of adjustment and device | |
| Deli et al. | Automated landmark extraction for orthodontic measurement of faces using the 3-camera photogrammetry methodology | |
| JP2012010891A (en) | Radiographic apparatus and method and program | |
| CN110084855B (en) | Improved CBCT geometric parameter calibration method | |
| US10856829B2 (en) | Reconstructing a three-dimensional image data record recorded with a biplanar x-ray device | |
| CN209032406U (en) | A geometric calibration device for cone beam CT system | |
| WO2015146166A1 (en) | Radiographic image capturing device, method, and program |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 17892761 Country of ref document: EP Kind code of ref document: A1 |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| 122 | Ep: pct application non-entry in european phase |
Ref document number: 17892761 Country of ref document: EP Kind code of ref document: A1 |