## Fast Implementations of Algebraic Methods for 3D Reconstruction from Cone-Beam Data (1998)

Venue: | IEEE Transactions on Medical Imaging |

Citations: | 4 - 3 self |

### BibTeX

@ARTICLE{Mueller98fastimplementations,

author = {Klaus Mueller and Roni Yagel and John J. Wheller},

title = {Fast Implementations of Algebraic Methods for 3D Reconstruction from Cone-Beam Data},

journal = {IEEE Transactions on Medical Imaging},

year = {1998},

volume = {18},

pages = {538--547}

}

### OpenURL

### Abstract

The prime motivation of this work is to devise techniques that make the Algebraic Reconstruction Technique (ART) and related methods more efficient for routine clinical use, while not compromising their accuracy. In particular, we strive to push the overall cost for a ART reconstruction as close as possible to the theoretical cost for a reconstruction obtained with Filtered Backprojection (FBP). While we focus mostly on fast implementations of ART-type methods in the context of 3D cone-beam reconstruction, different parts of the material presented here is also applicable to speed up reconstruction from fan-beam and parallel-beam data. It was shown in previous research that three iterations are sufficient to obtain a high quality reconstruction for low-contrast cone-beam. Based on the observation that ART typically only requires only half the projections of FBP, we conclude that if the overall cost for ART's projection-backprojection operations could be cut in half, then one ...

### Citations

442 | Footprint Evaluation for Volume Rendering
- Westover
- 1990
(Show Context)
Citation Context ...o new, highly accurate projection algorithms, one voxel-driven and one ray-driven, and analyze their efficiency in both the parallel-beam and cone-beam setting. Although other voxel-driven projectors =-=[28]-=- and ray-driven projectors [17][18] have been described, these algorithms are only efficient for the parallel-beam case or do not allow the stretched interpolation kernels prescribed in [21] as necess... |

296 |
Principles of Computerized Tomographic Imaging
- Kak, Slaney
- 2001
(Show Context)
Citation Context ...does not have a large set of projections available, or when the projections are sparse or missing at certain orientations, it was found that ART produces reconstructions of better quality than FBP [1]=-=[16]-=-. More recently, the importance of ART was reinforced by Matej et. al. [20] who found, for the reconstruction from noisy PET data, that ART produces quantitatively better results than the more popular... |

248 |
Summed-Area Tables for Texture Mapping
- Crow
- 1984
(Show Context)
Citation Context ...intersecting each of the four plane pairs with the footprint table we obtain four table indexes. The indexed table values are then combined according to the well-known summed area table (SAT) formula =-=[7]-=- such that a voxel weight w ij is given by (see Fig. 6b): (6) It must be added, however, that the method bears slight inaccuracies. This is due to the fact that for most pyramid orientations one or mo... |

165 |
Volume visualization
- Kaufman
- 1990
(Show Context)
Citation Context ...ects the volume on an image-basis. This was shown to significantly reduce the noise artifacts that were observed with ray-iterative ART. The projection step of SART performs a summed volume rendering =-=[14]-=- of the reconstruction grid, then subtracts the rendered image from the acquired projection image, normalizes the result, and backprojects the image in an inverse volume rendering process. More formal... |

147 | Practical Cone-beam Algorithm
- Feldkamp, David, et al.
- 1984
(Show Context)
Citation Context ...However, much of our discussion is also valid for the 2D and 3D reconstruction from parallel-beam and fan-beam data. Although a variety of general cone-beam algorithms based on FBP have been proposed =-=[8]-=-[9][25], there are, as of today, no clinical 3D cone-beam scanners. The present literature for 3D cone-beam reconstruction with ART is mostly restricted to the reconstruction of high-contrast objects,... |

124 |
Algebraic Reconstruction Techniques (ART) for three Dimensional Electron Micoscopy and X-ray Photography
- Bender, R, et al.
- 1970
(Show Context)
Citation Context ...ed projections onto a reconstruction grid [5]. The first and still most prevalent representative of the iterative methods is the Algebraic Reconstruction Technique (ART), attributed to Gordon et. al. =-=[10]-=-, while another well-known ART-type method is Simultaneous ART (SART), 2 proposed by Andersen and Kak [2]. In both methods, a reconstruction grid is iteratively updated by a projection-backprojection ... |

114 | Texture Splats for 3D Scalar and Vector Field Visualization
- Crawfis, Max
- 1993
(Show Context)
Citation Context ...e sampling frequency, where the signal's aliases are largest. 3. An AccurateVoxel-Driven Splatting Algorithm for Cone-Beam ART Let us first introduce some terminology. As suggested by Crawfis and Max =-=[6]-=-, we can think of the interpolation kernel footprint as a polygon with a superimposed texture map that is placed in object (volume) space. Here, the texture map is given by the projected kernel functi... |

97 |
The Fourier Reconstruction of a Head Section
- Shepp, Logan
- 1974
(Show Context)
Citation Context ...2] contrasted various previously published projection access schemes with a new scheme, the Weighted Distance Scheme. It was found that for low-contrast objects, such as the Shepp-Logan brain phantom =-=[27]-=-, usually 3 to 4 iterations are sufficient for good reconstruction quality. The choice of the relaxation coefficient l is another important parameter, which was studied by Herman and Meyer [13]. The e... |

87 |
The Fourier transform
- Bracewell
- 1989
(Show Context)
Citation Context ...roduction The field of 2D and 3D reconstruction methods can be roughly divided into two main categories. On one side there is the domain of direct methods that capitalize on the Fourier Slice Theorem =-=[4]-=-, while on the other side lies the domain of iterative methods that seek to solve the reconstruction problem by solving a system of simultaneous linear equations. The most prominent member of the form... |

72 |
Simultaneous algebraic reconstruction technique (SART): a superior implementation of the art algorithm
- Andersen, Kak
- 1984
(Show Context)
Citation Context ...erative methods is the Algebraic Reconstruction Technique (ART), attributed to Gordon et. al. [10], while another well-known ART-type method is Simultaneous ART (SART), 2 proposed by Andersen and Kak =-=[2]-=-. In both methods, a reconstruction grid is iteratively updated by a projection-backprojection procedure until a convergence criterion is satisfied. Two circumstances have led to the preference of FBP... |

51 |
Practical considerations for 3-D image reconstruction using spherically-symmetric volume elements
- Matej, Lewitt
- 1996
(Show Context)
Citation Context ...lgorithms, one voxel-driven and one ray-driven, and analyze their efficiency in both the parallel-beam and cone-beam setting. Although other voxel-driven projectors [28] and ray-driven projectors [17]=-=[18]-=- have been described, these algorithms are only efficient for the parallel-beam case or do not allow the stretched interpolation kernels prescribed in [21] as necessary for accurate cone-beam reconstr... |

49 |
Alternatives to voxels for image representation in iterative reconstruction algorithms
- Lewitt
- 1992
(Show Context)
Citation Context ...on algorithms, one voxel-driven and one ray-driven, and analyze their efficiency in both the parallel-beam and cone-beam setting. Although other voxel-driven projectors [28] and ray-driven projectors =-=[17]-=-[18] have been described, these algorithms are only efficient for the parallel-beam case or do not allow the stretched interpolation kernels prescribed in [21] as necessary for accurate cone-beam reco... |

35 |
Mathematical Framework of Cone Beam 3-D Reconstruction via the First Derivative of the Radon Transform
- Grangeat
- 1991
(Show Context)
Citation Context ...ever, much of our discussion is also valid for the 2D and 3D reconstruction from parallel-beam and fan-beam data. Although a variety of general cone-beam algorithms based on FBP have been proposed [8]=-=[9]-=-[25], there are, as of today, no clinical 3D cone-beam scanners. The present literature for 3D cone-beam reconstruction with ART is mostly restricted to the reconstruction of high-contrast objects, su... |

31 |
Efficient 3D grids for image reconstruction using sphericallysymmetric volume elements
- Matej, Lewitt
- 1995
(Show Context)
Citation Context .... This approach, however, is not appropriate for cone-beam reconstruction, as it does not allow independent footprint stretching in the two ray sheet directions. In another approach, Matej and Lewitt =-=[19]-=- decompose the voxel grid into a set of 2D slices. Here the orientation of the slices is that orientation most parallel to the image plane. Recall that a footprint is the pre-integrated kernel functio... |

27 |
Evaluation of task-oriented performance of several fully 3D PET reconstruction algorithms, Phys
- Matej, Herman, et al.
- 1994
(Show Context)
Citation Context ... are sparse or missing at certain orientations, it was found that ART produces reconstructions of better quality than FBP [1][16]. More recently, the importance of ART was reinforced by Matej et. al. =-=[20]-=- who found, for the reconstruction from noisy PET data, that ART produces quantitatively better results than the more popular FBP and MLE (Maximum Likelihood Estimation) methods. Finally, since the AR... |

25 |
Algebraic reconstruction in CT from limited views
- Andersen
- 1989
(Show Context)
Citation Context ...ne does not have a large set of projections available, or when the projections are sparse or missing at certain orientations, it was found that ART produces reconstructions of better quality than FBP =-=[1]-=-[16]. More recently, the importance of ART was reinforced by Matej et. al. [20] who found, for the reconstruction from noisy PET data, that ART produces quantitatively better results than the more pop... |

24 |
Image Reconstruction from Cone-beam Projections: Necessary and Sufficient Conditions and reconstruction methods
- Smith
- 1985
(Show Context)
Citation Context ...r, much of our discussion is also valid for the 2D and 3D reconstruction from parallel-beam and fan-beam data. Although a variety of general cone-beam algorithms based on FBP have been proposed [8][9]=-=[25]-=-, there are, as of today, no clinical 3D cone-beam scanners. The present literature for 3D cone-beam reconstruction with ART is mostly restricted to the reconstruction of high-contrast objects, such a... |

20 |
Local basisfunction approach to computed tomography
- Hanson, Wecksung
- 1985
(Show Context)
Citation Context ...ion algorithm is considerably more accurate than the one described by Westover [28]. Our ray-driven algorithm, on the other hand, is a 3D extension of the 2D algorithm proposed by Hanson and Wecksung =-=[12]-=-. However, a fast projection algorithm is not enough. We must reduce the actual complexity of the overall projection-backprojection framework. Ideally, we only want to do the computational equivalent ... |

11 |
Computed tomography using algebraic reconstruction techniques (ARTs) with different projection access schemes : a comparison study under practical situations
- Guan, Gordon
- 1996
(Show Context)
Citation Context ...ction-backprojection operations. Particularly encouraging in this respect is the observation that ART-type methods theoretically only require about half the number of projections than the FBP methods =-=[11]-=-. Thus we can allows the number of iterations for ART to grow larger than one and still maintain the same complexity as FBP. Several groups of researchers have worked on reducing the number of iterati... |

11 |
The influence of a relaxation parameter on SPECT iterative reconstruction algorithms,” Phys
- Ros, Falcon, et al.
- 1996
(Show Context)
Citation Context ...ART is mostly restricted to the reconstruction of high-contrast objects, such as encountered in computed angiography [23][26](small 3 number of projections) and reconstruction from PET [20] and SPECT =-=[24]-=- data (noisy projections). Since high-contrast objects usually require a significantly smaller number of iterations than lowcontrast object, it may not be fair to use the outcome of these studies for ... |

8 |
Direct Fourier methods in 3D reconstruction from cone-beam data
- Axelson
- 1994
(Show Context)
Citation Context ... scheme goes well with the variable-size voxel kernels since here the slab width is constant for z>z c . 6. Results Fig. 7a shows a slice from the 3D extension of the Shepp-Logan phantom (see [21] or =-=[3]-=- for details), reconstructed from 60 cone-beam data with ART using the traditional constant-size interpolation kernel. The aliasing artifacts caused by the diverging rays in cone-beam are easily notic... |

7 |
Algebraic reconstruction can be made computationally efficient
- Herman, Meyer
- 1993
(Show Context)
Citation Context ...hantom [27], usually 3 to 4 iterations are sufficient for good reconstruction quality. The choice of the relaxation coefficient l is another important parameter, which was studied by Herman and Meyer =-=[13]-=-. The effect of l was also examined by us in [21] where we also studied the effect of various other ART parameters, such as grid initialization and correction algorithm, for the cone-beam reconstructi... |

7 | The weighted distance scheme: a globally optimizing projection ordering method for the Algebraic Reconstruction Technique (ART - Mueller, Yagel, et al. - 1997 |

6 | Angenherte auflsung von systemen linearer gleichungen - Kaczmarz - 1937 |

6 |
Image intensifier-based volume angiography imaging system: work in progress
- Ning, Rooker
- 1993
(Show Context)
Citation Context ...3D cone-beam scanners. The present literature for 3D cone-beam reconstruction with ART is mostly restricted to the reconstruction of high-contrast objects, such as encountered in computed angiography =-=[23]-=-[26](small 3 number of projections) and reconstruction from PET [20] and SPECT [24] data (noisy projections). Since high-contrast objects usually require a significantly smaller number of iterations t... |

3 | Accurate low-contrast 3D cone-beam reconstruction with algebraic methods” in review
- Mueller, Yagel, et al.
- 1998
(Show Context)
Citation Context ...ient for good reconstruction quality. The choice of the relaxation coefficient l is another important parameter, which was studied by Herman and Meyer [13]. The effect of l was also examined by us in =-=[21]-=- where we also studied the effect of various other ART parameters, such as grid initialization and correction algorithm, for the cone-beam reconstruction of the 3D extension of the Shepp-Logan phantom... |

3 |
Rougee, "In vivo evaluation of a new system for 3D computerized angiography," Phys
- Saint-Felix, Trousset, et al.
- 1994
(Show Context)
Citation Context ...one-beam scanners. The present literature for 3D cone-beam reconstruction with ART is mostly restricted to the reconstruction of high-contrast objects, such as encountered in computed angiography [23]=-=[26]-=-(small 3 number of projections) and reconstruction from PET [20] and SPECT [24] data (noisy projections). Since high-contrast objects usually require a significantly smaller number of iterations than ... |