Reduction of streak artifacts in low dose CT imaging through multi image compounding
View Patent ↗Disclosed is a method and system for constructing, from a computerized tomography (CT) scan, an image relating to a physical structure. Projection data associated with the image is obtained and divided into a plurality of subsets. Filtered back projection (FBP) is then applied to each subset in the plurality of subsets. The image is constructed based on the application of the FBP to each subset in the plurality of subsets.
1. A method for constructing, from a computerized tomography scan, an image relating to a physical structure, comprising:
obtaining projection data associated with the physical structure;
dividing said projection data into a plurality of subsets;
applying filtered back projection (FBP) to each subset in said plurality of subsets to generate a plurality of sub-reconstructions, wherein each sub-reconstruction is an entire reconstruction of said image based on a respective subset in said plurality of subsets; and
constructing said image based on the said plurality of sub-reconstructions generated based on said plurality of subsets.
2. The method of claim 1 wherein said applying FBP to each subset further comprises obtaining an estimate of streaking noise for each sub-reconstruction.
3. The method of claim 2 further comprising comparing said estimate of streaking noise for each sub-reconstruction with at least a portion of said constructed image.
4. The method of claim 3 further comprising analyzing dependency between estimates of streaking noise for different sub-reconstructions.
5. The method of claim 1 further comprising regularizing an independence constraint between sub-reconstructions in said plurality of sub-reconstructions.
6. The method of claim 5 wherein said regularizing an independence constraint further comprises including a regularization term defined in the following equation:
e i,j ( {circumflex over (N)} i ,{circumflex over (N)} j )=∥ E ( h 1 ( {circumflex over (N)} i ) h 2 ( {circumflex over (N)} j ))− E ( h 1 ( {circumflex over (N)} i )) E ( h 2 ( {circumflex over (N)} j ))∥ 2
where E is the expectation, {circumflex over (N)} i and {circumflex over (N)} j are random variables, and h 1 ( ) and h 2 ( ) are functions of {circumflex over (N)} i and {circumflex over (N)} j , respectively.
7. The method of claim 6 wherein said constructing is formed from an energy minimization framework.
8. The method of claim 1 wherein said obtaining said projection data further comprises solving integrals along parallel rays directed along and through said image.
9. The method of claim 1 wherein said dividing of said projection data into a plurality of subsets further comprises down-sampling at least a portion of said projection data by a predetermined factor.
10. A system for constructing, from a computerized tomography scan, an image relating to a physical structure, comprising:
means for obtaining projection data associated with said physical structure;
means for dividing said projection data into a plurality of subsets;
means for applying filtered back projection (FBP) to each subset in said plurality of subsets to generate a plurality of sub-reconstructions, wherein each sub-reconstruction is an entire reconstruction of said image based on a respective subset in said plurality of subsets; and
means for constructing said image based on the said plurality of sub-reconstructions generated based on said plurality of subsets.
11. The system of claim 10 wherein said means for applying FBP to each subset further comprises means for obtaining an estimate of streaking noise for each sub-reconstruction.
12. The system of claim 11 further comprising means for comparing said estimate of streaking noise for each sub-reconstruction with at least a portion of said image after construction.
13. The system of claim 12 further comprising means for analyzing dependency between estimates of streaking noise for different sub-reconstructions.
14. The system of claim 10 further comprising means for regularizing an independence constraint between sub-reconstructions in said plurality of sub-reconstructions.
15. The system of claim 14 wherein said means for regularizing an independence constraint further comprises means for including a regularization term defined in the following equation:
e i,j ( {circumflex over (N)} i ,{circumflex over (N)} j )=∥ E ( h 1 ( {circumflex over (N)} i ) h 2 ( {circumflex over (N)} j ))− E ( h 1 ( {circumflex over (N)} i )) E ( h 2 ( {circumflex over (N)} j ))∥ 2
where E is the expectation, {circumflex over (N)} i and {circumflex over (N)} j are random variables, and h 1 ( ) and h 2 ( ) are functions of {circumflex over (N)} i and {circumflex over (N)} j , respectively.
16. The system of claim 15 wherein said means for constructing is formed from an energy minimization framework.
17. The system of claim 10 wherein said means for obtaining said projection data further comprises means for solving integrals along parallel rays.
18. The system of claim 10 wherein said means for dividing of said projection data into a plurality of subsets further comprises means for down-sampling at least a portion of said projection data by a predetermined factor.
19. A computer readable medium comprising computer program instructions capable of being executed in a processor and defining the steps comprising:
obtaining projection data associated with a physical structure;
dividing said projection data into a plurality of subsets;
applying filtered back projection (FBP) to each subset in said plurality of subsets to generate a plurality of sub-reconstructions, wherein each sub-reconstruction is an entire reconstruction of said image based on a respective subset in said plurality of subsets; and
constructing an image based on the said plurality of sub-reconstructions generated based on said plurality of subsets.
20. The computer readable medium of claim 19 wherein said applying FBP to each subset further comprises the step of obtaining an estimate of streaking noise for each sub-reconstruction.
21. The computer readable medium of claim 20 further comprising the step of comparing said estimate of streaking noise for each sub-reconstruction with at least a portion of said constructed image.
22. The computer readable medium of claim 21 further comprising the step of analyzing dependency between estimates of streaking noise for different sub-reconstructions.
23. The computer readable medium of claim 19 further comprising the step of regularizing an independence constraint between sub-reconstructions in said plurality of sub-reconstructions.
24. The computer readable medium of claim 23 wherein said regularizing an independence constraint further comprises including a regularization term defined in the following equation:
e i,j ( {circumflex over (N)} i ,{circumflex over (N)} j )=∥ E ( h 1 ( {circumflex over (N)} i ) h 2 ( {circumflex over (N)} j ))− E ( h 1 ( {circumflex over (N)} i )) E ( h 2 ( {circumflex over (N)} j ))∥ 2
where E is the expectation, {circumflex over (N)} i and {circumflex over (N)} j are random variables, and h 1 ( ) and h 2 ( ) are functions of {circumflex over (N)} i and {circumflex over (N)} j , respectively.
25. The computer readable medium of claim 24 wherein said step of constructing is formed from an energy minimization framework.
26. The computer readable medium of claim 19 wherein said step of obtaining said projection data further comprises solving integrals along parallel rays directed along and through said image.
27. The computer readable medium of claim 19 wherein said dividing of said projection data into a plurality of subsets further comprises the step of down-sampling at least a portion of said projection data by a predetermined factor.