CN103006251A - Calibration phantom, calibration device and calibration method for calibrating geometric parameters in CT (Computed Tomography) system - Google Patents

Abstract

The invention relates to a calibration phantom for calibrating geometric parameters in a CT (Computed Tomography) system. The calibration phantom comprises two same calibration plates which are arranged in parallel, and pellets, wherein the calibration plates are correspondingly provided with a plurality of same small holes respectively; and at least four pellets are arranged in the small holes in the calibration plates respectively. The invention relates to a calibration device for calibrating geometric parameters in a CT system. The calibration device comprises a CT system and the calibration phantom, wherein the CT system consists of a light source, a rotary platform and a detector; and the calibration phantom is vertically placed on the rotary platform between the light source and the detector. Furthermore, the invention further provides a calibration method for calibrating geometric parameters in a CT system. According to the calibration device and the calibration method, all seven geometric parameters of the CT system can be calibrated simultaneously by only performing one-time projection on the calibration phantom, so that operation is easy, and subsequent image correction of the CT system is met greatly.

Description

Be used for the CT system and demarcate demarcation phantom, caliberating device and the scaling method of geometric parameter

Technical field

The present invention relates to the demarcation of geometric parameter in the CT system, relate in particular to a kind of demarcation phantom, caliberating device and scaling method for CT system geometrical parameter calibration.

Background technology

Computer tomography (CT) is a kind of important imaging means that obtains internal structure of body information by lossless manner, and it has the inferior many merits of high-resolution, high sensitivity and multilamellar, is widely used in each medical clinical examination field.Geometrical parameter calibration is the important component part of micro-CT system debug, also is the precondition that obtains accurate reconstruction CT image.

The geometrical parameter calibration method of system " a kind of X ray gated cone-beam computed tomography " that Chinese patent application numbers 201210148432.8 proposes.The method reappears geometry site between x-ray source, flat panel detector and the rotating shaft in the cone-beam CT system by the camera calibration technology, thereby geometric parameter and error thereof to Cone-Beam CT are carried out direct solution, it is a kind of systematized measurement method for solving, can with Cone-Beam CT abstract be basic camera system model, thereby be mutually related simultaneously a plurality of geometric parameters in the solving system.

But the geometrical parameter calibration method of above-mentioned disclosed CT system can not be demarcated whole 7 geometric parameters of CT system simultaneously; And the CT system that can not be applicable to simultaneously exsomatize (demarcate the phantom rotation, light source and detector are motionless) and vivo CT system (demarcate phantom and fix, light source and detector rotation); In addition, the parameter calibration process is complicated, need to carry out multiple projections to demarcating phantom, is difficult to repeatedly repeat realize.

Summary of the invention

Based on this, be necessary the defective that the calibration system for geometric parameter in the above-mentioned CT system exists, a kind of demarcation phantom that can demarcate simultaneously geometric parameter in the CT system is provided.

A kind of demarcation phantom for CT system geometrical parameter calibration, comprise: comprise two the identical scaling board that be arranged in parallel and beads, correspondingly respectively on described two scaling boards offer some identical apertures, described bead is at least four and be installed on respectively in the aperture on described two scaling boards.

The present invention also provides a kind of caliberating device for CT system geometrical parameter calibration, comprise the CT system that is formed by light source, turntable and detector, described caliberating device also comprises the demarcation phantom on the turntable that vertically is placed between described light source and the detector, described demarcation phantom comprises two the identical scaling board that be arranged in parallel and beads, correspondingly respectively on described two scaling boards offer some identical apertures, described bead is at least four and be installed on respectively in the aperture on described two scaling boards; Described demarcation phantom rotates with described turntable, and described detector is used for the center point coordinate of the projected image of each bead of collection; And image processing module, be used for making up the first projection matrix according to the center point coordinate of described projected image, described image processing module also is used for calculating respectively each bead at the locus of described CT system coordinate system coordinate according to geometry and the described demarcation phantom geometry of described CT system, and make up the second projection matrix according to described locus coordinate, calculate the projection centre point coordinates of detector plane according to described the first projection matrix and described the second projection matrix, light source is to the distance of detector, the torsion angle of detector, the inclination angle of detector, the anglec of rotation of detector and light source are to the distance of turntable rotating shaft.

In addition, the present invention also provides the scaling method of geometric parameter in a kind of CT system, described CT system comprises light source, turntable and detector, described scaling method comprises the steps: that vertical the placement demarcated phantom on the turntable between described light source and the detector, described demarcation phantom comprises two the identical scaling board that be arranged in parallel and beads, correspondingly respectively on described two scaling boards offer some identical apertures, described bead is at least four and be installed on respectively in the aperture on described two scaling boards; Described demarcation phantom is rotated with described turntable, gather abscissa and the vertical coordinate of central point of the projected image of each bead by described detector; Center point coordinate according to described projected image makes up the first projection matrix; Calculate respectively the locus coordinate of each bead in described CT system coordinate system according to the geometric relationship of described CT system and the geometric relationship of demarcating phantom; Make up the second projection matrix according to described locus coordinate; According to described the first projection matrix and described the second projection matrix calculate the projection centre point coordinates, light source of detector plane to the inclination angle of the torsion angle of the distance of detector, detector, detector, the anglec of rotation of detector and light source be to the distance of turntable rotating shaft.

Above-mentioned caliberating device and scaling method for CT system geometric parameter, adopt the simple phantom of demarcating, rotate at turntable, demarcate the center point coordinate of the projected image of each bead on the phantom by the detector collection, center point coordinate according to projected image makes up the first projection matrix, geometry and demarcation phantom geometry according to the CT system calculate respectively the locus coordinate of each bead in described caliberating device coordinate system, and make up the second projection matrix according to described locus coordinate, calculate the projection centre point coordinates of detector plane according to the first projection matrix and the second projection matrix, light source is to the distance of detector, the torsion angle of detector, the inclination angle of detector, the anglec of rotation of detector and light source are to the distance of turntable rotating shaft.

Above-mentioned caliberating device and scaling method only need carry out whole 7 geometric parameters that 1 projection just can be demarcated the CT system simultaneously to demarcating phantom, and be simple to operate, greatly satisfied the image rectification of follow-up CT system; In addition, above-mentioned caliberating device and scaling method can be applicable to exsomatize CT system and vivo CT system, wide adaptability simultaneously; Simultaneously, above-mentioned nominal volume mode structure is simple, preparation cost is cheap, is easy to obtain.

Description of drawings

The sketch map of the demarcation phantom that is used for CT system geometrical parameter calibration that Fig. 1 provides for the embodiment of the invention one.

The structural representation of the aperture of offering on the scaling board that Fig. 2 provides for embodiment one provided by the invention.

The structural representation of the caliberating device that is used for CT system geometrical parameter calibration that Fig. 3 provides for the embodiment of the invention two.

The flow chart of the scaling method of geometric parameter in the CT system that Fig. 4 provides for the embodiment of the invention three.

Fig. 5 passes through the flow chart of steps of center point coordinate that detector gathers the projected image of each bead for what embodiment three provided by the invention provided.

Fig. 6 calculates respectively the locus coordinate flow chart of steps of each bead in the CT system coordinate system for what the embodiment of the invention three provided according to the geometric relationship of CT system and the geometric relationship of demarcation phantom.

The projected image that when embedding 18 prills on the demarcation phantom, obtains that Fig. 7 provides for the embodiment of the invention four.

The specific embodiment

In order to make purpose of the present invention, technical scheme and advantage clearer, below in conjunction with drawings and the specific embodiments, the present invention is further elaborated.Should be appreciated that specific embodiment described herein only in order to explain the present invention, is not intended to limit the present invention.

Embodiment one:

See also Fig. 1 and Fig. 2, the sketch map of the demarcation phantom that is used for CT system geometrical parameter calibration that provides for the embodiment of the invention one.

Demarcate phantom 100 and comprise scaling board 110 and bead 120.

The quantity of scaling board 110 is that two and structure are identical.Two scaling boards 110 over against and be arranged in parallel.In embodiment provided by the invention, the spacing of two scaling boards is preferably 46mm, and the width of each scaling board is preferably 150mm, thickness is preferably 2mm.Correspondingly respectively on two scaling boards 110 offer some identical apertures 111.Aperture 111 is the square formation formula and is distributed on the scaling board.See also Fig. 2, the structural representation of the aperture 111 of offering on the scaling board 110 that provides for embodiment one provided by the invention, the square formation that aperture 111 is 7*7 is distributed on the scaling board 110, and the diameter of aperture 111 is preferably 1mm, hole depth is preferably 0.5mm, and spacing is preferably 20mm between the adjacent holes.

Bead 120 is at least four.Bead 120 is installed on respectively in two apertures 111 on the scaling board 110, namely on any one scaling board 110 bead 120 is installed all.The quantity that is appreciated that aperture 120 can also be 6 or other more quantity.

In embodiment provided by the invention, the material of scaling board 110 is preferably vinyon.Bead 120 is preferably the rustless steel bead.

Embodiment two:

See also Fig. 3, the structural representation of the caliberating device that is used for CT system geometrical parameter calibration that Fig. 3 provides for the embodiment of the invention two.

Caliberating device 200 comprises CT system, demarcation phantom 100 and the image processing module (not shown) that is comprised of light source 210, turntable (not shown) and detector 230.

Be provided with turntable between light source 210 and the detector 230.Demarcating phantom 100 is vertically installed on the turntable.Demarcating phantom 100 rotates with turntable.Detector 230 is used for gathering the center point coordinate of the projected image of demarcating each bead 120 in the aperture 111 of offering on the phantom 100.The center point coordinate that is appreciated that the projected image of bead 120 on imaging plane does not all overlap.

Image processing module is used for making up the first projection matrix according to the center point coordinate of projected image.Image processing module also is used for according to the geometry of CT system and demarcates phantom 100 geometries and calculate respectively each bead 120 at the locus of CT system coordinate system coordinate, and makes up the second projection matrix according to the locus coordinate.Image processing module calculates the projection centre point coordinates on detector 230 planes, the distance that light source 210 arrives detector 230, the torsion angle of detector 230, the inclination angle of detector 230, the anglec of rotation of detector 230 and the distance that light source 210 arrives the turntable rotating shaft according to the first projection matrix and the second projection matrix.

Embodiment three:

See also Fig. 3 and Fig. 4, the flow chart of the scaling method of geometric parameter in the CT system that provides for the embodiment of the invention three of Fig. 4 wherein, wherein the CT system comprises light source 210, turntable and detector 230, specifically comprises the steps:

Step S310: the vertical placement demarcated phantom 100 on the turntable between light source 210 and the detector 230.In embodiment provided by the invention, it is identical with the demarcation phantom that the embodiment of the invention one provides to demarcate phantom 100.Be appreciated that bead 120 is at least four and be installed on respectively in two apertures 111 on the scaling board 110, namely demarcate on two scaling boards 110 of phantom 100 and all will place bead 120.

Step S320: will demarcate phantom 100 with turntable rotation, and gather abscissa and the vertical coordinate of central point of the projected image of each beads 120 by detector 230.

See also Fig. 5, pass through the flow chart of steps of center point coordinate that detector 230 gathers the projected image of each beads for what embodiment three provided by the invention provided, it comprises the steps:

Step S321: the plane at definition detector 230 places is as coordinate plane.Wherein take detector 230 lower left corners as zero, orientation U, the V of pixel is respectively axis of abscissas and axis of ordinates on the detector 230.

Step S322: will demarcate phantom 100 and rotate with turntable.Be appreciated that the bead 120 of demarcating in the aperture 111 of offering on the phantom 100 at the uniform velocity rotates with turntable, carries out projection by 210 pairs of beads of light source 120, and obtains projected image.In embodiment provided by the invention, the number of times of projection is preferably 1 time, is appreciated that the number of times of projection can also be for repeatedly, and when the number of times of projection was more, the geometric parameter of the CT system that finally obtains was more accurate.

Step S323: detector 230 gathers the projected image of each bead 120.Be appreciated that each projection bead all obtains 1 projected image at detector 230, namely bead is all obtaining projected image take detector 230 place coordinate planes as imaging plane.

Step S324: the central point that obtains successively each bead 120 projected images by the edge extracting processing method is with the plane at detector 230 places abscissa and the vertical coordinate as coordinate plane.The center point coordinate that is appreciated that the projected image of bead 120 on imaging plane does not all overlap.

Step S330: central point abscissa and vertical coordinate according to projected image make up the first projection matrix.In embodiment provided by the invention, it is specific as follows to make up the first projection matrix according to the central point abscissa of projected image and vertical coordinate:

With the abscissa of the central point of the projected image of each bead 120 and vertical coordinate substitution matrix [u, v, 1] respectively ^TIn, make up the first projection matrix [u, v, 1] ^T, wherein, u and v represent respectively abscissa and the vertical coordinate of bead 120 projected centre point, [u, v, 1] ^TBe the matrix of 3xN, N is the bead number.

Step S340: calculate respectively the locus coordinate of each bead in the CT system coordinate system according to the geometric relationship of CT system and the geometric relationship of demarcation phantom.

See also Fig. 6, calculate respectively the locus coordinate flow chart of steps of each bead 120 in the CT system coordinate system for what the embodiment of the invention three provided according to the geometric relationship of CT system and the geometric relationship of demarcation phantom 100, specifically comprise the steps:

Step S341: definition coordinate system (X, Y, Z), wherein Y-axis is the rotating shaft of turntable, Z axis is that light source 210 is to the ray of rotating shaft.

Step S342: the axis at the central point place of demarcating phantom 100 is overlapped with rotating shaft.In another embodiment of demarcation phantom 100 provided by the invention, the axis that the central point place of phantom 100 is demarcated on the edge is provided with a metal cylinder, and this metal cylinder overlaps with rotating shaft, and this metal cylinder can rotate around rotating shaft.

Step S343: the rotation turntable makes the plane parallel at two scaling boards 110 and detector 230 places, can obtain the locus coordinate (x, y, z) of each bead on coordinate system (X, Y, Z).Be appreciated that because to demarcate the structure of phantom 100 be known, demarcate like this relative tertiary location of each bead 120 on the phantom 100 and determine, will demarcate phantom 100 and be fixed on the turntable, and make that metal cylinder overlaps with the rotating shaft Y-axis in the demarcation phantom 100.The rotation turntable, so that the plane at two scaling board 110 places of demarcation phantom 100 and the plane parallel at detector 230 places, light source 210 is respectively the center on two planes at the vertical point of two scaling boards 110 demarcating phantom 100, can obtain rustless steel bead 120 at coordinate system (X, Y, Z) the locus coordinate (x, y, z) on.

Step S350: make up the second projection matrix according to the locus coordinate.In embodiment provided by the invention, wherein, according to rustless steel bead 120 at coordinate system (X, Y, Z) the locus coordinate (x, y, z) on makes up the second projection matrix, specifically comprise the steps: the locus coordinate (x with each bead, y, z) difference substitution matrix [x, y, z, 1] ^TIn, make up the second projection matrix [x, y, z, 1] ^TWherein, x, y, z are the locus coordinate (x, y, z) of each bead.

Step S360: according to the first projection matrix and the second projection matrix calculate the projection centre point coordinates, light source of detector plane to the inclination angle of the torsion angle of the distance of detector, detector, detector, the anglec of rotation of detector and light source be to the distance of turntable rotating shaft.

According to the first projection matrix and the second projection matrix calculate the projection centre point coordinates, light source of detector plane to the inclination angle of the torsion angle of the distance of detector, detector, detector, the anglec of rotation of detector and light source be to the flow chart of steps of the distance of turntable rotating shaft, specifically comprise the steps:

Definition coordinate system (X, Y, Z), wherein Y-axis is the rotating shaft of turntable, Z axis is that light source 210 is to the ray of rotating shaft.

Definition detector 230 lower left corners are zero, and U, V are the orientation of pixel on the detector 230, and λ is the size, (u of detector pixel ₀, v ₀) be the projection centre point coordinates on detector 230 planes, D is the distance that light source 210 arrives detector 230, and α is that the torsion angle of detector 230, inclination angle, the γ that β is detector 230 are the anglec of rotation of detector 230, and R is that light source 210 is to the distance of rotating shaft.

According to the first projection matrix [u, v, 1] ^TAnd the second projection matrix [x, y, z, 1] ^TMake up formula a[u, v, 1] ^T=A[x, y, z, 1] ^T, wherein, A is sytem matrix, a is coefficient.Be appreciated that a can rule of thumb choose.

According to formula A=Φ [ζ | ω] and calculate Φ, ζ, wherein,

$\mathrm{\φ}=\left[\begin{array}{ccc}\frac{D}{\mathrm{\λ}}& 0& u_0\\ 0& \frac{D}{\mathrm{\λ}}& v_0\\ 0& 0& 1\\ & & \end{array}\right]$

$\mathrm{\ζ}=\left[\begin{array}{ccc}\cos \mathrm{\β}\cos \mathrm{\γ}& \sin \mathrm{\α}\sin \mathrm{\β}\cos \mathrm{\γ}-\cos \mathrm{\α}\sin \mathrm{\γ}& \sin \mathrm{\α}\sin \mathrm{\γ}+\cos \mathrm{\α}\sin \mathrm{\β}\cos \mathrm{\γ}\\ \cos \mathrm{\β}\sin \mathrm{\γ}& \cos \mathrm{\α}\cos \mathrm{\γ}+\sin \mathrm{\α}\sin \mathrm{\β}\sin \mathrm{\γ}& \cos \mathrm{\α}\sin \mathrm{\β}\sin \mathrm{\γ}-\sin \mathrm{\α}\cos \mathrm{\γ}\\ -\sin \mathrm{\β}& \sin \mathrm{\α}\cos \mathrm{\β}& \cos \mathrm{\α}\cos \mathrm{\β}\end{array}\right]$

$\mathrm{\ω}=\left[\begin{array}{c}{\mathrm{\ω}}_1\\ {\mathrm{\ω}}_2\\ {\mathrm{\ω}}_3\end{array}\right]$

ω ₁＝A ₃₄、ω ₂＝(A ₂₄-Φ ₂₃A ₃₄)/Φ ₂₃、ω ₂＝(A ₂₄-Φ ₂₃A ₃₄)/Φ ₂₃、ω ₂＝(A ₁₄-Φ ₁₃A ₃₄-ω ₂A ₁₂)/Φ ₁₁。In the above-mentioned formula, D is that light source is to the distance of detector, (u ₀, v ₀) be light source 210 at the projection centre point coordinates on detector 230 planes, λ is the size of detector 230 pixels, α is that the torsion angle of detector 230, inclination angle, the γ that β is detector 230 are the anglec of rotation of detector 230.

Be appreciated that and obtain first the locus coordinate (x of N bead 120 on coordinate system (X, Y, Z) _i, y _i, z _i), i=1,2,3 ... N and with the plane at detector 230 places abscissa and the vertical coordinate (u as coordinate plane _i, v _i), i=1,2,3 ... N.By formula a[u, v, 1] ^T=A[x, y, z, 1] ^T, following equation is arranged:

u _ia _i＝A ₁₁x _i+A ₁₂y _i+A ₁₃z _i+A ₁₄

v _ia _i＝A ₂₁x _i+A ₂₂y _i+A ₂₃z _i+A ₂₄

ω _i＝A ₃₁x _i+A ₃₂y _i+A ₃₃z _i+A ₃₄

Coefficient a disappears _i, obtain following two equations:

A ₁₁x _i+A ₁₂y _i+A ₁₃z _i+A ₁₄-u _i(A ₃₁x _i+A ₃₂y _i+A ₃₃z _i+A ₃₄)＝0

A ₂₁x _i+A ₂₂y _i+A ₂₃z _i+A ₂₄-v _i(A ₃₁x _i+A ₃₂y _i+A ₃₃z _i+A ₃₄)＝0

Write above equation as matrix form, obtained: HA*=0, wherein,

$H=\left[\begin{array}{cccccccccccc}{x}_i& y_i& z_i& 1& 0& 0& 0& 0& -u_i{x}_i& {-u}_i{y}_i& -u_i{z}_i& -u_i\\ 0& 0& 0& 0& x_i& y_i& z_i& 1& -v_i{x}_i& -v_i{y}_i& -v_i{z}_i& -v_i\end{array}\right]$

$A^{*}=\left[\begin{array}{c}{A}_{11}\\ A_{12}\\ A_{13}\\ A_{14}\\ A_{21}\\ A_{22}\\ A_{23}\\ A_{24}\\ A_{31}\\ A_{32}\\ A_{33}\\ A_{34}\end{array}\right]$

For equation HA*=0, can pass through singular value decomposition method, matrix H is decomposed, obtain A*, get final product to get sytem matrix A.

According to formula A=Φ [ζ | ω], can obtain A _{3 * 3}=Φ ζ (namely having removed the last string of the 3x4 of matrix A), because Φ is upper triangular matrix, can be to matrix A _{3 * 3}Carry out " QR decomposition " and obtain Φ, ζ.

According to formula o=[o _x, o _y, o _z] ^T=-ζ ^Tω, calculating XYZ coordinate is center point coordinate o _x, o _yAnd o _z

Wherein, in the above-mentioned formula, ζ can be obtained by above-mentioned steps,

$\mathrm{\ω}=\left[\begin{array}{c}{\mathrm{\ω}}_1\\ {\mathrm{\ω}}_2\\ {\mathrm{\ω}}_3\end{array}\right]$

ω ₁＝A ₃₄、ω ₂＝(A ₂₄-Φ ₂₃A ₃₄)/Φ ₂₃、ω ₂＝(A ₂₄-Φ ₂₃A ₃₄)/Φ ₂₃、ω ₂＝(A ₁₄-Φ ₁₃A ₃₄-ω ₂A ₁₂)/Φ ₁₁。

Projection centre point coordinates (the u on calculating detector 230 planes ₀, v ₀), light source 210 arrives the distance B of detector 230, the anglec of rotation γ of the torsion angle α of detector 230, the angle of inclination beta of detector 230, detector 230, and light source 210 is to the distance R of rotating shaft.

Be appreciated that ζ and XYZ coordinate are center point coordinate o according to Φ obtained above _x, o _yAnd o _z, can calculate the projection centre point coordinates (u on detector 230 planes ₀, v ₀), light source 210 arrives the distance B of detector 230, the anglec of rotation γ of the torsion angle α of detector 230, the angle of inclination beta of detector 230, detector 230, and light source 210 is to the distance R of rotating shaft.

Wherein, the projection centre point coordinates (u on detector 230 planes ₀, v ₀) computing formula be: u ₀=Φ ₁₃, v ₀=Φ ₂₃, Φ ₁₃The 1st row of representing matrix Φ and element corresponding to the 3rd row, Φ ₂₃The 1st row of representing matrix Φ and element corresponding to the 3rd row.

Light source 210 to the computing formula of the distance B of detector 230 is: D=Φ ₁₁λ, Φ ₁₁The 1st row of representing matrix Φ and element corresponding to the 1st row.

The computing formula of the torsion angle α of detector 230 is: α=Arctan2 (ζ ₃₂, ζ ₃₃), ζ ₃₂The 3rd row of representing matrix ζ and element corresponding to the 2nd row, ζ ₃₃The 3rd row of representing matrix ζ and element corresponding to the 3rd row.

The computing formula of the angle of inclination beta of detector 230 is: β=sin ^-1(ζ ₃₁), ζ ₃₁The 3rd row of representing matrix ζ and element corresponding to the 1st row.

The computing formula of the anglec of rotation γ of detector 230 is: γ=Arctan2 (ζ ₁₂, ζ ₃₁), ζ _{, 12}The 1st row of representing matrix ζ and element corresponding to the 2nd row, ζ ₃₁The of representing matrix ζ, 3 row and element corresponding to the 1st row.

Light source 210 to the computing formula of the distance R of rotating shaft is: o _x, o _y, o _zFor XYZ coordinate is center point coordinate.

Embodiment four:

See also Fig. 7, for the embodiment of the invention four provide on demarcating phantom, embed 18 prills the time projected image that obtains.Situation when demarcating 18 prills of phantom 100 embeddings.Each correspondence is put into 9 beads, the projected image that obtains on two scaling boards up and down respectively.For the effectiveness of verification algorithm, press extreme case, only carried out 1 projection (equaling 1).In theory, when value was larger, the result can be more close to actual value.System's geometric parameter analogue value that the above-mentioned scaling method Calculation Simulation of foundation obtains and actual value are as shown in Table 1.As can be seen from the results, the method precision is very high.

Form 1

?	R	D	α	β	γ	u 0	v 0
								Actual value	570	1005	-1.5708	-1.5708	0.00	257	257
Value of calculation	571.7760	1008.5622	-1.5708	-1.5705	0.00	257	255.0620
								Error	0.312％	0.354％	0％	0.019％	0％	0％	0.754％

The scaling method of geometric parameter and caliberating device adopt and demarcate the phantom rotation in the above-mentioned CT system, and light source and detector are motionless, thereby realize the demarcation of geometric parameter in the CT system; Be appreciated that the demarcation that can realize equally geometric parameter in the CT system for vivo CT system (demarcate phantom and fix, light source and detector rotation).

Above-mentioned caliberating device and scaling method for CT system geometric parameter, adopt the simple phantom of demarcating, rotate at turntable, demarcate the center point coordinate of the projected image of each bead on the phantom by the detector collection, center point coordinate according to projected image makes up the first projection matrix, geometry and demarcation phantom geometry according to the CT system calculate respectively the locus coordinate of each bead in described caliberating device coordinate system, and make up the second projection matrix according to described locus coordinate, calculate the projection centre point coordinates of detector plane according to the first projection matrix and the second projection matrix, light source is to the distance of detector, the torsion angle of detector, the inclination angle of detector, the anglec of rotation of detector and light source are to the distance of turntable rotating shaft.

Above-mentioned caliberating device and scaling method only need carry out whole 7 geometric parameters that 1 projection just can be demarcated the CT system simultaneously to demarcating phantom, and be simple to operate, greatly satisfied the image rectification of follow-up CT system; In addition, above-mentioned caliberating device and scaling method can be applicable to exsomatize CT system and vivo CT system, wide adaptability simultaneously; Simultaneously, above-mentioned nominal volume mode structure is simple, preparation cost is cheap, is easy to obtain.

The above, it only is preferred embodiment of the present invention, be not that the present invention is done any pro forma restriction, although the present invention discloses as above with preferred embodiment, yet be not to limit the present invention, any those skilled in the art, within not breaking away from the technical solution of the present invention scope, when the technology contents that can utilize above-mentioned announcement is made a little change or is modified to the equivalent embodiment of equivalent variations, in every case be not break away from the technical solution of the present invention content, any simple modification that foundation technical spirit of the present invention is done above embodiment, equivalent variations and modification all still belong in the scope of technical solution of the present invention.

Claims (13)

1. one kind is used for the demarcation phantom that geometric parameter is demarcated by the CT system, it is characterized in that, described demarcation phantom comprises two the identical scaling board that be arranged in parallel and beads, correspondingly respectively on described two scaling boards offer some identical apertures, described bead is at least four and be installed on respectively in the aperture on described two scaling boards.

2. the demarcation phantom for CT system demarcation geometric parameter according to claim 1 is characterized in that the spacing of described two scaling boards is 46mm, and the width of described scaling board is that 150mm, thickness are 2mm.

3. the demarcation phantom for CT system demarcation geometric parameter according to claim 1 is characterized in that described aperture is the square formation formula and is distributed on the described scaling board.

4. the demarcation phantom of demarcating geometric parameter for the CT system according to claim 3, it is characterized in that, the square formation that described aperture is 7*7 is distributed on the described scaling board, and the diameter of described aperture is that 1mm, hole depth are 0.5mm, and spacing is 20mm between the adjacent holes.

5. the demarcation phantom for CT system demarcation geometric parameter according to claim 1 is characterized in that described scaling board is vinyon, and described bead is the rustless steel bead.

6. the caliberating device of geometric parameter in the CT system, comprise the CT system that is formed by light source, turntable and detector, it is characterized in that, described caliberating device also comprises the demarcation phantom on the turntable that vertically is placed between described light source and the detector, described demarcation phantom comprises two the identical scaling board that be arranged in parallel and beads, correspondingly respectively on described two scaling boards offer some identical apertures, described bead is at least four and be installed on respectively in the aperture on described two scaling boards; Described demarcation phantom rotates with described turntable, and described detector is used for the center point coordinate of the projected image of each bead of collection; And

Image processing module, be used for making up the first projection matrix according to the center point coordinate of described projected image, described image processing module also is used for calculating respectively each bead at the locus of described CT system coordinate system coordinate according to geometry and the described demarcation phantom geometry of described CT system, and make up the second projection matrix according to described locus coordinate, calculate the projection centre point coordinates of detector plane according to described the first projection matrix and described the second projection matrix, light source is to the distance of detector, the torsion angle of detector, the inclination angle of detector, the anglec of rotation of detector and light source are to the distance of turntable rotating shaft.

7. the caliberating device of geometric parameter in the CT according to claim 6 system is characterized in that the center point coordinate of the projected image of described bead on imaging plane does not all overlap.

8. the scaling method of geometric parameter in the CT system, described CT system comprises light source, turntable and detector, it is characterized in that described scaling method comprises the steps:

The vertical placement demarcated phantom on the turntable between described light source and the detector, described demarcation phantom comprises two the identical scaling board that be arranged in parallel and beads, correspondingly respectively on described two scaling boards offer some identical apertures, described bead is at least four and be installed on respectively in the aperture on described two scaling boards;

Described demarcation phantom is rotated with described turntable, gather abscissa and the vertical coordinate of central point of the projected image of each bead by described detector;

Central point abscissa and vertical coordinate according to described projected image make up the first projection matrix;

Calculate respectively the locus coordinate of each bead in described CT system coordinate system according to the geometric relationship of described CT system and the geometric relationship of demarcating phantom;

Make up the second projection matrix according to described locus coordinate;

According to described the first projection matrix and described the second projection matrix calculate the projection centre point coordinates, light source of detector plane to the inclination angle of the torsion angle of the distance of detector, detector, detector, the anglec of rotation of detector and light source be to the distance of turntable rotating shaft.

9. the scaling method of geometric parameter in the CT according to claim 8 system is characterized in that, wherein, the step of center point coordinate that gathers the projected image of each bead by described detector specifically comprises:

Define the plane at described detector place as coordinate plane; Wherein the detector lower left corner is zero, and orientation U, the V of pixel is respectively axis of abscissas and axis of ordinates on the described detector;

Described demarcation phantom is rotated with described turntable;

Described detector gathers the projected image of each bead;

The central point that obtains successively each bead projected image by the edge extracting processing method is with the plane at described detector place abscissa and the vertical coordinate as coordinate plane.

10. the scaling method of geometric parameter according to claim 8 or in the 9 described CT systems is characterized in that, wherein, the step that makes up the first projection matrix according to the center point coordinate of described projected image specifically comprises:

With the abscissa of the central point of the projected image of each bead and vertical coordinate substitution matrix [u, v, 1] respectively ^TIn, make up the first projection matrix [u, v, 1] ^T, wherein, u and v represent respectively abscissa and the vertical coordinate of bead projected centre point.

11. the scaling method of geometric parameter in the CT according to claim 8 system, it is characterized in that, wherein, calculate respectively the locus coordinate of each bead according to the geometric relationship of described CT system and the geometric relationship of demarcating phantom, comprise the steps:

Definition coordinate system (X, Y, Z), wherein Y-axis is the rotating shaft of described turntable, Z axis is that described light source is to the ray of described rotating shaft;

The axis at the central point place of described demarcation phantom is overlapped with described rotating shaft;

The rotation turntable makes the plane parallel at described two scaling boards and described detector place, can obtain the locus coordinate (x, y, z) of each bead on coordinate system (X, Y, Z).

12. the scaling method of geometric parameter is characterized in that according to claim 8 or in the 11 described CT systems, wherein, makes up the second projection matrix according to described locus coordinate, comprises the steps:

Locus coordinate (x, y, z) difference substitution matrix [x, y, z, 1] with each bead ^TIn, make up the second projection matrix [x, y, z, 1] ^T

13. the scaling method of geometric parameter is characterized in that in the CT according to claim 8 system, wherein, and according to described the first projection matrix [u, v, 1] ^TAnd described the second projection matrix [x, y, z, 1] ^TCalculate the projection centre point coordinates, light source of detector plane to the inclination angle of the torsion angle of the distance of detector, detector, detector, the anglec of rotation of detector and light source be to the distance of turntable rotating shaft, comprise the steps:

Definition coordinate system (X, Y, Z), wherein Y-axis is the rotating shaft of described turntable, Z axis is that described light source is to the ray of described rotating shaft;

Defining the described detector lower left corner is zero, and U, V are the orientation of pixel on the described detector, and λ is the size, (u of detector pixel ₀, v ₀) be the projection centre point coordinates of described detector plane, D be light source to the distance of detector, α is that the torsion angle of detector, inclination angle, the γ that β is detector are the anglec of rotation of detector, R is that light source is to the distance of rotating shaft;

According to described the first projection matrix [u, v, 1] ^TAnd described the second projection matrix [x, y, z, 1] ^TMake up formula a[u, v, 1] ^T=A[x, y, z, 1] ^T, wherein, A is sytem matrix, a is coefficient;

According to formula A=Φ [ζ | ω] and calculate Φ, ζ, wherein, $\mathrm{\φ}=\left[\begin{array}{ccc}\frac{D}{\mathrm{\λ}}& 0& u_0\\ 0& \frac{D}{\mathrm{\λ}}& v_0\\ 0& 0& 1\\ & & \end{array}\right]$

$\mathrm{\ζ}=\left[\begin{array}{ccc}\cos \mathrm{\β}\cos \mathrm{\γ}& \sin \mathrm{\α}\sin \mathrm{\β}\cos \mathrm{\γ}-\cos \mathrm{\α}\sin \mathrm{\γ}& \sin \mathrm{\α}\sin \mathrm{\γ}+\cos \mathrm{\α}\sin \mathrm{\β}\cos \mathrm{\γ}\\ \cos \mathrm{\β}\sin \mathrm{\γ}& \cos \mathrm{\α}\cos \mathrm{\γ}+\sin \mathrm{\α}\sin \mathrm{\β}\sin \mathrm{\γ}& \cos \mathrm{\α}\sin \mathrm{\β}\sin \mathrm{\γ}-\sin \mathrm{\α}\cos \mathrm{\γ}\\ -\sin \mathrm{\β}& \sin \mathrm{\α}\cos \mathrm{\β}& \cos \mathrm{\α}\cos \mathrm{\β}\end{array}\right]$

$\mathrm{\ω}=\left[\begin{array}{c}{\mathrm{\ω}}_1\\ {\mathrm{\ω}}_2\\ {\mathrm{\ω}}_3\end{array}\right]$

ω ₁＝A ₃₄、ω ₂＝(A ₂₄-Φ ₂₃A ₃₄)/Φ ₂₃、ω ₂＝(A ₂₄-Φ ₂₃A ₃₄)/Φ ₂₃、ω ₂＝(A ₁₄-Φ ₁₃A ₃₄-ω ₂A ₁₂)/Φ ₁₁；

According to formula o=[o _x, o _y, o _z] ^T=-ζ ^Tω, calculating XYZ coordinate is center point coordinate o _x, o _yAnd o _z

Calculate the projection centre point coordinates (u of described detector plane ₀, v ₀), light source is to the distance B of detector, the torsion angle α of detector, the angle of inclination beta of detector, the anglec of rotation γ of detector, and light source is to the distance R of rotating shaft;

Wherein, the projection centre point coordinates (u of detector plane ₀, v ₀) computing formula be: u ₀=Φ ₁₃, v ₀=Φ ₂₃, Φ ₁₃The 1st row of representing matrix Φ and element corresponding to the 3rd row, Φ ₂₃The 1st row of representing matrix Φ and element corresponding to the 3rd row,

Light source to the computing formula of the distance B of detector is: D=Φ ₁₁λ, Φ ₁₁The 1st row of representing matrix Φ and element corresponding to the 1st row,

The computing formula of the torsion angle α of detector is: α=Arctan2 (ζ ₃₂, ζ ₃₃), ζ ₃₂The 3rd row of representing matrix ζ and element corresponding to the 2nd row, ζ ₃₃The 3rd row of representing matrix ζ and element corresponding to the 3rd row,

The computing formula of the angle of inclination beta of detector is: β=sin ^-1(ζ ₃₁), ζ ₃₁The 3rd row of representing matrix ζ and element corresponding to the 1st row,

The computing formula of the anglec of rotation γ of detector is: γ=Arctan2 (ζ ₁₂, ζ ₃₁), ζ _{, 12}The 1st row of representing matrix ζ and element corresponding to the 2nd row, ζ ₃₁The of representing matrix ζ, 3 row and element corresponding to the 1st row,

Light source to the computing formula of the distance R of rotating shaft is:

o _x, o _y, o _zFor XYZ coordinate is center point coordinate.

Images