Computed tomography apparatus and empirical pre-weighting method for decreasing image noise nonuniformity

Patent metadata

Field	Value
Publication number	US-10055860-B2
Application number	US-201615167521-A
Country	US
Kind code	B2
Filing date	May 27, 2016
Priority date	May 27, 2016
Publication date	Aug 21, 2018
Grant date	Aug 21, 2018

Abstract

A computed tomography (CT) method and apparatus including a radiation source configured to produce radiation directed to an object space, and a plurality of detector elements configured to detect the radiation produced from the radiation source through the object space and generate projection data. A rotation mount is configured to rotate the radiation source around the object space. Processing circuitry is configured to cause the rotation mount to rotate the radiation source, and to receive the projection data. The projection data includes a plurality of projection data sets. The processing circuitry calculates a set of weights based on the projection data sets, calculates a set of pre-weights based on the weights, and minimizes a penalized weighted least-squares cost function to produce a reconstructed image. The cost function is a sum of a weighted least-squares term, weighted using the weights, and a penalty term weighted using the pre-weights.

First claim

The invention claimed is:

1. A method for computed tomography (CT) image reconstruction, comprising:
acquiring a plurality of projection data sets generated by a plurality of detector elements that detect X-ray radiation directed to an object by an X-ray source;
calculating a set of weights based on the acquired projection data sets;
calculating a set of pre-weights based on the calculated weights; and
minimizing a penalized weighted least-squares cost function to produce a reconstructed image, the cost function being a sum of a weighted least-squares term and a penalty term,
wherein the weighted least-squares term is weighted using the calculated set of weights,
wherein the penalty term is weighted using the calculated set of pre-weights, and
wherein the penalty term includes a penalty function defined for each voxel of the reconstructed image.

2. The method according to claim 1 , wherein the step of calculating the set of weights includes calculating statistical information from the plurality of projection data sets.

3. The method according to claim 1 , wherein the penalty term includes a sum of the penalty function evaluated at each of the voxels of the reconstructed image.

4. The method according to claim 1 , wherein
the step of calculating the set of weights includes calculating elements of a square weighting matrix, and
the step of calculating the set of pre-weights includes applying a backprojection operator to a product of a function of diagonal elements of the weighting matrix and X-ray penetration terms calculated by summing the rows of a forward projection operator.

5. The method according to claim 4 , wherein the function of the diagonal elements of the weighting matrix is a square root of the diagonal elements of the weighting matrix.

6. A method of computed tomography (CT) image reconstruction, comprising:
acquiring a plurality of projection data sets generated by a plurality of detector elements that detect X-ray radiation directed to an object by an X-ray source;
calculating a set of weights based on the acquired projection data sets;
calculating a set of pre-weights based on the projection data sets; and
minimizing a penalized weighted least-squares cost function to produce a reconstructed image, the cost function being a sum of a weighted least-squares term and a penalty term,
wherein the weighted least-squares term is weighted using the calculated set of weights,
wherein the penalty term is weighted using the calculated set of pre-weights, and
wherein the penalty term includes a penalty function defined for each voxel of the reconstructed image.

7. The method according to claim 6 , wherein the step of calculating the set of weights includes calculating statistical information from the plurality of projection data sets.

8. The method according to claim 6 , wherein the penalty term includes a sum of the penalty function evaluated at each of the voxels of the reconstructed image.

9. The method according to claim 6 , wherein
the step of calculating the set of weights includes calculating matrix elements of a square weighting matrix,
the step of calculating the set of pre-weights includes estimating a variance of volume elements of the reconstructed image, and
the pre-weights are a function of the variance of the volume elements.

10. The method according to claim 9 , wherein the pre-weights are a square root of the variance of the volume elements.

11. A computed tomography (CT) apparatus, comprising:
a radiation source configured to produce radiation directed to an object space;
a plurality of detector elements configured to detect the radiation produced by the radiation source, the plurality of detector elements being configured to generate projection data;
a rotation mount configured to rotate the radiation source around the object space, the radiation source being fixedly connected to the rotation mount; and
processing circuitry configured to
receive the projection data from the plurality of detector elements, the projection data including a plurality of projection data sets;
calculate a set of weights based on the received projection data sets;
calculate a set of pre-weights based on the calculated weights; and
minimize a penalized weighted least-squares cost function to produce a reconstructed image, the cost function being a sum of a weighted least-squares term and a penalty term,
wherein the weighted least-squares term is weighted using the calculated set of weights,
wherein the sum in the penalty term is weighted using the calculated set of pre-weights, and
wherein the penalty term includes a penalty function defined for each voxel of the reconstructed image.

12. The CT apparatus according to claim 11 , wherein the processing circuitry is further configured to calculate the set of weights using statistical information calculated from the plurality of projection data sets.

13. The CT apparatus according to claim 11 , wherein the processing circuitry is further configured to calculate the penalty term using a sum of the penalty function evaluated at each of the voxels of the reconstructed image.

14. The CT apparatus according to claim 11 , wherein
the processing circuitry is further configured to
calculate the set of weights using calculated elements of a square weighting matrix, and
calculate the set of pre-weights by applying a backprojection operator to a product of a function of diagonal elements of the weighting matrix and X-ray penetration terms calculated by summing the rows of a forward projection operator.

15. The CT apparatus according to claim 14 , wherein the function of the diagonal elements of the weighting matrix is a square root of the diagonal elements of the weighting matrix.

Assignees

• Toshiba Medical Sys Corp

Inventors

• Shi Daxin
• Zou Yu

Classifications

• G06T12/20 Primary

  Inverse problem, i.e. transformations from projection space into object space · CPC title

• A61B6/5205 Primary

  involving processing of raw data to produce diagnostic data · CPC title

• A61B6/032

  Transmission computed tomography [CT] · CPC title

• G06T2207/20056

  Discrete and fast Fourier transform, [DFT, FFT] · CPC title

• G06T2207/10081

  Computed x-ray tomography [CT] · CPC title