X-ray CT device and method for correcting scattered X-rays
View Patent ↗It is an object to prevent a poor X-ray CT image due to scattered X-rays. An X-ray CT scanner ( 100 ) uses an X-ray detector ( 4 ), etc., to scan (F 1 ) an object ( 3 ) to obtain X-ray transmission image data, and estimate (F 3 ) an X-ray absorption coefficient distribution inside the object ( 3 ) based on the scanned X-ray transmission image data. Next, the X-ray CT scanner ( 100 ) performs a Monte Carlo simulation on a simulated object ( 3 ) having the X-ray absorption coefficient distribution estimated, and estimates (F 4 and F 7 ) a point spread function or scattered X-ray distribution derived from the object. Then, the X-ray CT scanner ( 100 ) corrects (F 5 and F 8 ), based on the point spread function or scattered X-ray distribution estimated, the X-ray transmission image data, and finally creates (F 6 and F 9 ) an image of the X-ray absorption coefficient distribution of the object ( 3 ).
1. An X-ray CT scanner comprising:
a scanning section that scans an object to obtain X-ray transmission image data of the object from a plurality of projection directions, the scanning section comprising an X-ray source generating X-rays from an X-ray focus and an X-ray detector having X-ray detection elements in a two-dimensional array so as to detect the X-rays, wherein the X-ray source and the X-ray detector rotate around the object while facing each other and the object is interposed therebetween;
an internal-distribution-estimating section that estimates an X-ray absorption coefficient distribution inside the object, based on the X-ray transmission image data scanned by the scanning section;
a point-spread-function-estimating section that estimates a point spread function of a scatter derived from the object by performing a Monte Carlo simulation to simulate a physical interaction of the X-rays in a simulated object having the X-ray absorption coefficient distribution estimated by the internal-distribution-estimating section;
a correction section that corrects the X-ray transmission image data by processing the point spread function estimated by the point-spread-function-estimating section and the X-ray transmission image date according to a deconvolution integral method; and
an image-creating section that creates an image of the X-ray absorption coefficient distribution of the object by using the X-ray transmission image data corrected by the correction section.
2. The X-ray CT scanner according to claim 1 , wherein the point-spread-function-estimating section performs the Monte Carlo simulation with respect to representative projection directions, the number of which is less than the number of the plurality of projection directions; the point spread functions obtained by the Monte Carlo simulation are used for interpolation between the representative projection directions to estimate point spread functions regarding all the plurality of projection directions.
3. The X-ray CT scanner according to claim 1 , wherein the point-spread-function-estimating section performs the Monte Carlo simulation with respect to representative element positions, the number of which is predetermined in view of the X-ray detection elements in the two-dimensional array; the point spread functions obtained by the Monte Carlo simulation are used for interpolation between the representative element positions to estimate point spread functions regarding all the positions of the X-ray detection elements in the two-dimensional array.
4. The X-ray CT scanner according to claim 1 , wherein the internal-distribution-estimating section simulates the object as a mixture of a plurality of substances having a known CT number to estimate the X-ray absorption coefficient distribution.
5. The X-ray CT scanner according to claim 1 , further comprising:
a user interface allowing for at least one operation selected from the group consisting of: a setting of X-ray statistics for the Monte Carlo simulation; a setting of a calculation time for execution of the Monte Carlo simulation; and an interruption at an any desired timing during the Monte Carlo simulation.
6. The X-ray CT scanner according to claim 1 , wherein the internal-distribution-estimating section estimates the X-ray absorption coefficient distribution inside the object from an X-ray absorption coefficient distribution of a simulated object that simulates the object instead of estimating, based on the X-ray transmission image data scanned by the scanning section, the X-ray absorption coefficient distribution inside the object.
7. An X-ray CT scanner comprising:
a scanning section that scans an object to obtain X-ray transmission image data of the object from a plurality of projection directions, the scanning section comprising an X-ray source generating X-rays from an X-ray focus and an X-ray detector having X-ray detection elements in a two-dimensional array so as to detect the X-rays, wherein the X-ray source and the X-ray detector rotate around the object while facing each other and the object is interposed therebetween;
an internal-distribution-estimating section that estimates an X-ray absorption coefficient distribution inside the object, based on the X-ray transmission image data scanned by the scanning section;
an X-ray-distribution-estimating section that estimates a distribution of scattered X rays derived from the object by performing a Monte Carlo simulation to simulate a physical interaction of the X-rays in a simulated object having the X-ray absorption coefficient distribution estimated by the internal-distribution-estimating section;
a correction section that removes components of the scattered X-rays from the X-ray transmission image data, based on the distribution estimated by the X-ray-distribution-estimating section; and
an image-creating section that creates an image of the X-ray absorption coefficient distribution of the object by using the X-ray transmission image data in which the components of the scattered X-rays have been removed by the correction section.
8. The X-ray CT scanner according to claim 7 , wherein the X-ray-distribution-estimating section performs the Monte Carlo simulation with respect to representative projection directions, the number of which is less than the number of the plurality of projection directions; the distributions obtained by the Monte Carlo simulation are used for interpolation between the representative projection directions to estimate distributions regarding all the plurality of projection directions.
9. The X-ray CT scanner according to claim 8 , wherein the X-ray-distribution-estimating section estimates a projection direction having a maximal change in an amount of the scattered X-rays from the distributions of the representative projection directions and repeats a process for determining a scattered X-ray distribution of the projection direction having the maximal change by performing the Monte Carlo simulation; and all the representative projection directions are sequentially selected.
10. A method for scatter correction, comprising:
a scanning step of scanning an object to obtain X-ray transmission image data of the object from a plurality of projection directions, a scanning section including an X-ray source generating X-rays from an X-ray focus and an X-ray detector having X-ray detection elements in a two-dimensional array so as to detect the X-rays, wherein the X-ray source and the X-ray detector rotate around the object while facing each other and the object is interposed therebetween;
an internal-distribution-estimating step of estimating an X-ray absorption coefficient distribution inside the object, based on the X-ray transmission image data scanned in the scanning step;
a point-spread-function-estimating step of estimating a point spread function of a scatter derived from the object by performing a Monte Carlo simulation to simulate a physical interaction of the X-rays in a simulated object having the X-ray absorption coefficient distribution estimated in the internal-distribution-estimating step;
a correction step of correcting the X-ray transmission image data by processing the point spread function estimated in the point-spread-function-estimating step and the X-ray transmission image date according to a deconvolution integral method; and
an image-creating step of creating an image of the X-ray absorption coefficient distribution of the object by using the X-ray transmission image data corrected in the correction step.
11. A method for scatter correction, comprising:
a scanning step of scanning an object to obtain X-ray transmission image data of the object from a plurality of projection directions, a scanning section including an X-ray source generating X-rays from an X-ray focus and an X-ray detector having X-ray detection elements in a two-dimensional array so as to detect the X-rays, wherein the X-ray source and the X-ray detector rotate around the object while facing each other and the object is interposed therebetween;
an internal-distribution-estimating step of estimating an X-ray absorption coefficient distribution inside the object, based on the X-ray transmission image data scanned in the scanning step;
an X-ray-distribution-estimating step of estimating a distribution of scattered X rays derived from the object by performing a Monte Carlo simulation to simulate a physical interaction of the X-rays in a simulated object having the X-ray absorption coefficient distribution estimated in the internal-distribution-estimating step;
a correction step of removing components of the scattered X-rays from the X-ray transmission image data, based on the distribution estimated in the X-ray-distribution-estimating step; and
an image-creating step of creating an image of the X-ray absorption coefficient distribution of the object by using the X-ray transmission image data in which the components of the scattered X-rays have been removed in the correction step.