Method and arrangement for medical X-ray imaging
First Claim
1. A method for producing three-dimensional information of an object (4) in medical X-ray imaging, characterized in thatthe object is X-radiated from at least two different directions and the said X-radiation is detected to form projection data of the object (4)the object is modelled mathematically utilizing the projection data to solve the imaging geometry and/or the motion of the object, where the solving concerns either some or all parts of the imaging geometry and/or the motion of the objectand said projection data and said mathematical modelling of the object are utilized in Bayesian inversion based on Bayes'"'"' formula
-
( x , θ | m ) = p pr ( θ ) p pr ( x ) p ( m | x , θ ) p ( m )
to produce three-dimensional information of the object, the prior distribution ppr(θ
) representing the prior knowledge of the imaging geometry and/or the motion of the object (4), the prior distribution ppr(x) representing mathematical modelling of the object, x representing the object image vector, which comprises values of the X-ray attenuation coefficient inside the object, θ
representing the parameter vector of the imaging geometry and/or the motion of the object (4), m representing projection data, the likelihood distribution p(m|x,θ
) representing the X-radiation attenuation model between the object image vector x, geometry parameter vector θ
, and projection data m, p(m) being a normalization constant and the posteriori distribution p(x,θ
|m) representing the three-dimensional information of the object (4) and the imaging geometry including the motion of the object.
1 Assignment
0 Petitions
Accused Products
Abstract
A medical X-ray device 5 arrangement for producing three-dimensional information of an object 4 in a medical X-ray imaging comprises an X-ray source 2 for X-radiating the object from at least two different directions; a detector 6 for detecting the X-radiation to form projection data of the object 4; a computational device 15 for modelling the object 4 mathematically utilizing the projection data to solve the imaging geometry and/or the motion of the object, where the solving concerns either some or all parts of the imaging geometry and/or the motion of the object. The computational device 15 utilizes said projection data and said mathematical modelling of the object in Bayesian inversion based on Bayes'"'"' formula
to produce three-dimensional information of the object.
16 Citations
36 Claims
-
1. A method for producing three-dimensional information of an object (4) in medical X-ray imaging, characterized in that
the object is X-radiated from at least two different directions and the said X-radiation is detected to form projection data of the object (4) the object is modelled mathematically utilizing the projection data to solve the imaging geometry and/or the motion of the object, where the solving concerns either some or all parts of the imaging geometry and/or the motion of the object and said projection data and said mathematical modelling of the object are utilized in Bayesian inversion based on Bayes'"'"' formula -
( x , θ | m ) = p pr ( θ ) p pr ( x ) p ( m | x , θ ) p ( m )
to produce three-dimensional information of the object, the prior distribution ppr(θ
) representing the prior knowledge of the imaging geometry and/or the motion of the object (4), the prior distribution ppr(x) representing mathematical modelling of the object, x representing the object image vector, which comprises values of the X-ray attenuation coefficient inside the object, θ
representing the parameter vector of the imaging geometry and/or the motion of the object (4), m representing projection data, the likelihood distribution p(m|x,θ
) representing the X-radiation attenuation model between the object image vector x, geometry parameter vector θ
, and projection data m, p(m) being a normalization constant and the posteriori distribution p(x,θ
|m) representing the three-dimensional information of the object (4) and the imaging geometry including the motion of the object. - View Dependent Claims (2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18)
-
-
19. A medical X-ray device (5) arrangement for producing three-dimensional information of an object (4) in a medical X-ray imaging, characterized in that the medical X-ray device (5) arrangement comprises:
-
an X-ray source (2) for X-radiating the object from at least two different directions a detector (6) for detecting the X-radiation to form projection data of the object (4) means (15) for modelling the object (4) mathematically utilizing the projection data to solve the imaging geometry and/or the motion of the object, where the solving concerns either some or all parts of the imaging geometry and/or the motion of the object and the medical X-ray device (5) arrangement includes means (15) for utilizing said projection data and said mathematical modelling of the object in Bayesian inversion based on Bayes'"'"' formula
to produce three-dimensional information of the object, the prior distribution ppr(θ
) representing the prior knowledge of the imaging geometry and/or the motion of the object (4), the prior distribution ppr(x) representing mathematical modelling of the object, x representing the object image vector, which comprises values of the X-ray attenuation coefficient inside the object, θ
representing the parameter vector of the imaging geometry and/or the motion of the object (4), m representing projection data, the likelihood distribution p(m|x,θ
) representing the X-radiation attenuation model between the object image vector x, geometry parameter vector θ
, and projection data m, p(m) being a normalization constant and the posteriori distribution p(x,θ
|m) representing the three-dimensional information of the object (4) and the imaging geometry including the motion of the object.- View Dependent Claims (20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36)
where matrix Aθ
contains the weights describing how much each voxel contributes to the X-ray attenuation along the X-ray paths and the noise e is independent of object image vector x and the geometry parameter vector θ
leading to the likelihood distribution
p(m|x,θ
)=pnoise(m−
Aθ
x).
-
-
23. A medical X-ray device (5) arrangement according to claim 19 characterized in that the medical X-ray device arrangement comprises means (15) for modelling the object (4) mathematically so that X-radiation attenuates when passing the object (4), which means that every image voxel is nonnegative.
-
24. A medical X-ray device (5) arrangement according to claim 19, characterized in that the medical X-ray device arrangement comprises means (15) for modelling the object (4) mathematically by the formula:
-
where the sum is taken over a collection of 3D neighbourhoods N and the value UN(x) depends only on the values of voxels belonging to the neighborhood N, and α
is a positive regularization parameter used to tune the width of the prior distribution.
-
-
25. A medical x-ray device (5) arrangement according to claim 19, characterized in that the 3D tomographic problem is divided into a stack of 2D tomographic problems and on every such 2D problem, the medical X-ray device arrangement comprises means (15) for modelling the object (4) mathematically by the formula:
-
where the sum is taken over a collection of 2D neighbourhoods N and the value UN(x) depends only on the values of pixels belonging to the neighborhood N, and α
is a positive regularization parameter used to tune the width of the prior distribution, and the 2D tomographic problems are related to each other by the formula
pr3D(x(j))=exp(−
γ
Σ
Σ
|x(j)[k,q]−
x(j−
1)[k,q]|),where the sums are taken over all pixels (k=1, . . . , K, q=1, . . . , Q) and γ
>
0 is another regularization parameter.
-
-
26. A medical X-ray device (5) arrangement according to claim 25, characterized in that the neighborhood systems comprise two neighboring pixels xj, xk or voxels xj, xk and UN(x) calculates a power of the
-
( x ( j ) ) = exp ( - α ( ∑ k = 1 K - 1 ∑ q = 1 Q x ( j ) [ k , q ] - x ( j ) [ k + 1 , q ] 8 ++ ∑ k = 1 K ∑ q = 1 Q - 1 x ( j ) [ k , q ] - x ( j ) [ k , q + 1 ] 8 ) ) absolute value of the difference, leading to the formula where s is a positive real number and α
is a regularization parameter used to tune the width of the prior distribution.
-
-
27. A medical X-ray device (5) arrangement according to claim 26, characterized in that the medical X-ray device arrangement comprises means (15) for modelling the object (4) mathematically by setting s=1 corresponding to total variation (TV) prior describing objects consisting of different regions with well-defined boundaries.
-
28. A medical X-ray device (5) arrangement according to claim 19, characterized in that the medical X-ray device arrangement comprises means (15) for modelling the object (4) mathematically by assuming that mathematical modelling is qualitative structural information of the target where the structural information is encoded in prior distributions that are concentrated around image vectors x that correspond to the physiological structures of the target.
-
29. A medical X-ray device (5) arrangement according to claim 19, characterized in that the medical X-ray device arrangement comprises means (15) for modelling the object (4) mathematically by assuming that mathematical modelling comprises a list of possible attenuation coefficient values in the object.
-
30. A medical X-ray device (5) arrangement according to claim 19, characterized in that the medical X-ray device arrangement comprises means (15) for modelling the object (4) mathematically by assuming that the X-ray imaging geometry, such as X-ray source position, has unknown error modelled in the distribution p(m|x,θ
- ).
-
31. A medical X-ray device (5) arrangement according to claim 19, characterized in that the medical X-ray device arrangement comprises means (15) for modelling the object (4) mathematically by assuming that X-radiation measurement noise is Poisson distributed.
-
32. A medical X-ray device (5) arrangement according to claim 19, characterized in that the medical X-ray imaging is dental radiography.
-
33. A medical X-ray device (5) arrangement according to claim 19, characterized in that the medical X-ray imaging is surgical C-arm imaging.
-
34. A medical X-ray device (5) arrangement according to claim 19, characterized in that the medical X-ray imaging is mammography.
-
35. A medical X-ray device (5) arrangement according to claim 19, characterized in that the medical X-ray device arrangement comprises means (15) for producing three-dimensional information of the object (4) on the basis of the maximum a posteriori estimator (MAP), which is calculated by the equation:
-
p(yMAP|m)=maxp(y|m)m representing projection data and y=(x,θ
) representing the vector composed of the object image vector and the geometry parameter vector, and where the maximum on the right hand side of the equation is taken over all y.
-
-
36. A medical X-ray device (5) arrangement according to claim 19, characterized in that the medical X-ray device arrangement comprises means (15) for producing three-dimensional information of the object (4) on the basis of the conditional mean estimator (CM), which is calculated by the equation
yCM=∫- yp(y|m)dy
where m represents projection data and y=(x,θ
) represents the vector composed of the object image vector and the geometry parameter vector.
- yp(y|m)dy
Specification