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

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.