O(N3 log N) backprojection algorithm for the 3-D radon transform.
We present a novel backprojection algorithm for three-dimensional (3-D) radon transform data that requires O(N3 log2 N) operations for reconstruction of an N x N x N volume from O(N2) plane-integral projections. Our algorithm uses a hierarchical decomposition of the 3-D radon transform to recursively decompose the backprojection operation. Simulations are presented demonstrating reconstruction quality comparable to the standard filtered backprojection, which requires O(N5) computations under the same circumstances.
Basu
Samit
S
Coordinated Science Laboratory, University of Illinois at Urbana-Champaign, 61801, USA.
Bresler
Yoram
Y
D000465:Algorithms; D006257:Head / Q000033:anatomy & histology / Q000000981:diagnostic imaging; D006801:Humans; D021621:Imaging, Three-Dimensional / Q000295:instrumentation* / Q000379:methods*; D008962:Models, Theoretical; D019047:Phantoms, Imaging; D011856:Radiographic Image Enhancement / Q000379:methods*; D015203:Reproducibility of Results; D012680:Sensitivity and Specificity; D014057:Tomography, X-Ray Computed / Q000379:methods*