## Anti-Aliased Three-Dimensional Cone-Beam Reconstruction of Low-Contrast Objects with Algebraic Methods (1999)

Venue: | IEEE Trans. Med. Imag |

Citations: | 10 - 1 self |

### BibTeX

@ARTICLE{Mueller99anti-aliasedthree-dimensional,

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

title = {Anti-Aliased Three-Dimensional Cone-Beam Reconstruction of Low-Contrast Objects with Algebraic Methods},

journal = {IEEE Trans. Med. Imag},

year = {1999},

volume = {18},

pages = {519--537}

}

### OpenURL

### Abstract

This paper examines the use of the algebraic reconstruction technique (ART) and related techniques to reconstruct 3-D objects from a relatively sparse set of cone-beam projections. Although ART has been widely used for cone-beam reconstruction of high-contrast objects, e.g., in computed angiography, the work presented here explores the more challenging low-contrast case which represents a little-investigated scenario for ART. Preliminary experiments indicate that for cone angles greater than 20 ffiffiffi , traditional ART produces reconstructions with strong aliasing artifacts. These artifacts are in addition to the usual off-midplane inaccuracies of cone-beam tomography with planar orbits. We find that the source of these artifacts is the nonuniform reconstruction grid sampling and correction by the cone-beam rays during the ART projection--backprojection procedure. A new method to compute the weights of the reconstruction matrix is devised, which replaces the usual constant-size in...

### Citations

479 |
The Mathematics of Computerized Tomography
- NATTERER
- 1986
(Show Context)
Citation Context ...ne-beam ART in recent years is in 3-D computed angiography. Here, one acquires a limited set of projection images of blood-perfused structures, such as vascular trees in the head or abdominal regions =-=[29]-=-, [34]. It should be noted, however, that the objects reconstructed in 3-D computed angiography are of rather high contrast, which poses the reconstruction problem as almost a binary one. Another nota... |

463 | Footprint evaluation for volume rendering
- Westover
- 1990
(Show Context)
Citation Context ...el and accumulating the interpolated values. Since accurate integration requires many sampling points, this is very time consuming. A more efficient way was proposed by many authors [13], [22], [23], =-=[41], -=-in the context of parallel-beam and fan-beam ART. It consists of reordering the ray integral so that each voxel’s contribution to the integral can be viewed isolated from the other voxels. To achiev... |

352 |
Principles of Computerized Tomographic Imaging
- Kak, Slaney
- 1988
(Show Context)
Citation Context ... photon scattering artifacts in the reconstruction procedure, FBP tends to either produce inferior results ors520 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 18, NO. 6, JUNE 1999 simply cannot be used =-=[20]-=-. We can think of a variety of clinical and industrial applications where one or more of these conditions arise, i.e., when projections are few, sparse, or missing. A clinical example for the latter a... |

215 |
Image Reconstruction from Projections: the Fundamentals of Computerized Tomography
- Herman
- 1980
(Show Context)
Citation Context ...rently gives rise to a significantly slower reconstruction speed and is the main reason for ART’s under-utilization in present clinical applications. While the research on general 2-D ART [11], [12]=-=, [14]-=-, [17] and 2-D fan-beam ART [13], [26] is numerous, the published literature on 3-D cone-beam reconstructors using ART-type algorithms is rather sparse. Exceptions are the early works by Colsher [9] a... |

172 |
Volume Visualization
- Kaufman
- 1991
(Show Context)
Citation Context ...tifacts that were observed with ray-iterative ART. (The same author also demonstrated SART’s strength in the limited angle problem [4].) The projection step of SART performs a summed volume renderin=-=g [18]-=- 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... |

146 |
Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and x-ray photography.”J. Theoret
- Gordon, Bender, et al.
- 1970
(Show Context)
Citation Context ...Under these conditions, an alternative reconstruction method has been shown to have a great margin of advantage. The algebraic reconstruction technique (ART), originated by Gordon, Bender, and Herman =-=[11] (-=-see [15] for a discussion of ART’s utility in practical cases where projections are limited and/or are noisy). In contrast to FBP, ART is an iterative procedure, i.e., it works by iteratively updati... |

118 |
Fourier Reconstruction of a Head Section
- Shepp, Logan
- 1974
(Show Context)
Citation Context ...entioned before, a 1 In this context, we define low-contrast objects as objects that have features of little variation in density (i.e., that have low contrast), such as the Shepp–Logan brain phanto=-=m [35]-=-, with a dynamic range of the main features of only 2.0%. This definition was also used by Tam [39]. scenario that has not been studied much in the past. We will see that the application of the standa... |

98 |
The Fourier Transform and Its
- Bracewell
- 2000
(Show Context)
Citation Context ...encies of up to but not greater than . Then, since is considered an ideal box in the frequency domain with bandwidth ,it removes all aliases of and we can write (9) as follows: using the relationship =-=[7]-=- comb (8) (9) (10) comb (11) In the parallel-beam case, for all and the aliases of in will not overlap. In that case, (attenuated by a nonideal ). However, in the cone-beam case there is a chance that... |

98 |
Angenäherte Auflösung von Systemen linearer Gleichungen. Bulletin de l’Académie Polonaise des Sciences et Lettres
- Kaczmarz
- 1937
(Show Context)
Citation Context ...ion or least squares methods. In addition, noise and sampling errors in the ART implementation normally do not provide for a consistent equation system. Thus, an iterative scheme proposed by Kaczmarz =-=[19]-=- is used. Starting from an initial guess for the volume vector we select at each iteration step , one of the equations in (1): say the one for .A value is measured which is the value of pixel computed... |

83 |
Simultaneous algebraic reconstruction technique (SART): a superior implementation of the art algorithm, Ultrason
- Andersen, Kak
- 1984
(Show Context)
Citation Context ...terpolation filter that is employed during volume projection and backprojection. In that respect, the quality of the interpolation filter determines the accuracy of the weights. While several authors =-=[5]-=-, [22], [23] have operated in this framework to determine the accurate weights for parallel-beam reconstruction, the scenario of diverging rays, as occurring in cone-beam reconstruction, imposes new c... |

59 |
Practical considerations for 3-D image reconstruction using spherically symmatric volume elements
- Matej, Lewitt
- 1996
(Show Context)
Citation Context ... filter that is employed during volume projection and backprojection. In that respect, the quality of the interpolation filter determines the accuracy of the weights. While several authors [5], [22], =-=[23]-=- have operated in this framework to determine the accurate weights for parallel-beam reconstruction, the scenario of diverging rays, as occurring in cone-beam reconstruction, imposes new constraints o... |

56 |
Alternatives to voxels for image representations in iterative reconstruction algorithms’, Phys
- Lewitt
- 1992
(Show Context)
Citation Context ...lation filter that is employed during volume projection and backprojection. In that respect, the quality of the interpolation filter determines the accuracy of the weights. While several authors [5], =-=[22]-=-, [23] have operated in this framework to determine the accurate weights for parallel-beam reconstruction, the scenario of diverging rays, as occurring in cone-beam reconstruction, imposes new constra... |

29 |
Evaluation of task-oriented performance of several fully 3d PET reconstruction algorithms
- Matej, Herman, et al.
- 1994
(Show Context)
Citation Context ...he objects reconstructed in 3-D computed angiography are of rather high contrast, which poses the reconstruction problem as almost a binary one. Another notable recent publication is the one by Matej =-=[25]-=-, whose studies indicate that ART also has significant merit for noisy projection data. Matej showed for PET that ART can produce quantitatively better reconstruction results than the more popular FBP... |

28 |
Algebraic reconstruction in ct from limited views
- Andersen
- 1989
(Show Context)
Citation Context ...by Andersen [5], which was shown to significantly reduce the noise artifacts that were observed with ray-iterative ART. (The same author also demonstrated SART’s strength in the limited angle proble=-=m [4]-=-.) The projection step of SART performs a summed volume rendering [18] of the reconstruction grid, then subtracts the rendered image from the acquired projection image, normalizes the result, and back... |

25 |
A general cone-beam reconstruction algorithm
- Wang, Lin, et al.
(Show Context)
Citation Context ... the past 15 years cone-beam researchers have mostly focused on reconstruction algorithms based on the filtered backprojection (FBP) approach (see Smith [37] for a review of these algorithms and Wang =-=[40] f-=-or a more recent paper on practical implementations of noncircular source orbits). This focus can be attributed to FBP’s convenient analytical formulation which enables fast computation of predictab... |

22 |
Local basis-function approach to computed tomography
- Hanson, Wecksung
- 1985
(Show Context)
Citation Context ...ntly slower reconstruction speed and is the main reason for ART’s under-utilization in present clinical applications. While the research on general 2-D ART [11], [12], [14], [17] and 2-D fan-beam AR=-=T [13]-=-, [26] is numerous, the published literature on 3-D cone-beam reconstructors using ART-type algorithms is rather sparse. Exceptions are the early works by Colsher [9] and Schlindwein [36]. However, th... |

19 |
ART: mathematics and applications (a report on the mathematical foundations and on the applicability to real data of the algebraic reconstruction techniques
- Herman, Lent, et al.
- 1973
(Show Context)
Citation Context ...e conditions, an alternative reconstruction method has been shown to have a great margin of advantage. The algebraic reconstruction technique (ART), originated by Gordon, Bender, and Herman [11] (see =-=[15] f-=-or a discussion of ART’s utility in practical cases where projections are limited and/or are noisy). In contrast to FBP, ART is an iterative procedure, i.e., it works by iteratively updating a recon... |

15 | Exact cone beam CT with a spiral scan
- Tam, Samarasekera, et al.
- 1998
(Show Context)
Citation Context ...of little variation in density (i.e., that have low contrast), such as the Shepp–Logan brain phantom [35], with a dynamic range of the main features of only 2.0%. This definition was also used by Ta=-=m [39]-=-. scenario that has not been studied much in the past. We will see that the application of the standard ART algorithm in the conebeam setting produces strong aliasing-related noise artifacts for cone ... |

14 |
An anti-aliasing technique for splatting
- Mueller, Möller, et al.
- 1997
(Show Context)
Citation Context ..., � is the Bessel filter described by Matej [23]. interpolation. An alternative method is to simply decrease the bandwidth of to which gives rise to a filter comb (13) This technique was also used i=-=n [38]-=- to achieve accurate perspective volume rendering with the splatting technique [41]. The frequency response of is shown in Fig. 5(a) for the slice at where . The frequency response of is also shown. I... |

13 |
Three dimensional reconstruction in nuclear medicine emission imaging
- Budinger, Gullberg
- 1974
(Show Context)
Citation Context ... Altschuler proposed two true cone-beam algorithms for the DSR, one using an analytic series expansion approach [1], [2] and one using an iterative Bayesian framework [3]. In unrelated work, Budinger =-=[8]-=- developed a solution based on least squares. It is apparent that over the past 15 years cone-beam researchers have mostly focused on reconstruction algorithms based on the filtered backprojection (FB... |

13 |
Cone-Beam Tomography: Recent Advances and a
- Smith
- 1990
(Show Context)
Citation Context ...ased on least squares. It is apparent that over the past 15 years cone-beam researchers have mostly focused on reconstruction algorithms based on the filtered backprojection (FBP) approach (see Smith =-=[37] f-=-or a review of these algorithms and Wang [40] for a more recent paper on practical implementations of noncircular source orbits). This focus can be attributed to FBP’s convenient analytical formulat... |

12 |
Computed tomography scanning with simultaneous patient translation, Med Phys 17
- Crawford, King
- 1990
(Show Context)
Citation Context ...y patient movement or breathing are much less of an issue than in slower forms of volumetric CT, such as the stacks-of-(fanbeam-)slices representation or the more recent helical/spiral CT (see, e.g., =-=[10]-=-). Likewise, imaging of dynamic structures such as the human heart is also greatly facilitated [32]. Cone-beam imaging received much attention with the construction of the dynamic spatial reconstructo... |

12 | A relaxation method for reconstructing objects from noisy X-rays - Herman |

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 ...rk currently gives rise to a significantly slower reconstruction speed and is the main reason for ART’s under-utilization in present clinical applications. While the research on general 2-D ART [11]=-=, [12]-=-, [14], [17] and 2-D fan-beam ART [13], [26] is numerous, the published literature on 3-D cone-beam reconstructors using ART-type algorithms is rather sparse. Exceptions are the early works by Colsher... |

11 |
The influence of a relaxation parameter on SPECT iterative reconstruction algorithms,” Phys
- Ros, Falcon, et al.
- 1996
(Show Context)
Citation Context ...a cone-beam reconstructor was not used, but the projection rays were rebinned, which simplified ART to the parallel-beam case. Another group of researchers has successfully applied ART for SPECT data =-=[33]-=-. That ART can produce superior results in the presence of noise was also demonstrated in an early paper by Herman [15]. However, this was found only to be true in the limited projection case. In anot... |

10 |
Interpolative algebraic reconstruction techniques without beam partitioning for computed tomography
- Mazur, Gordon
- 1995
(Show Context)
Citation Context ...lower reconstruction speed and is the main reason for ART’s under-utilization in present clinical applications. While the research on general 2-D ART [11], [12], [14], [17] and 2-D fan-beam ART [13]=-=, [26]-=- is numerous, the published literature on 3-D cone-beam reconstructors using ART-type algorithms is rather sparse. Exceptions are the early works by Colsher [9] and Schlindwein [36]. However, these im... |

9 | Cornhill, The Weighted Distance Scheme: A Globally Optimizing Projection Ordering Method for ART
- Mueller, Yagel, et al.
- 1997
(Show Context)
Citation Context ...taken.) Thus, in SART all rays of a projection are simultaneously processed, hence the name simultaneous ART. Apart from the correction algorithm and also from the projection access order (see, e.g., =-=[28]-=-), there are other important parameters that influence both reconstruction quality and speed of convergence. One of these factors is the relaxation coefficient . Natterer [29] gives exact bounds on th... |

9 |
In vivo evaluation of a new system for 3D computerized angiography
- Saint-Félix, Trousset, et al.
- 1994
(Show Context)
Citation Context ...m ART in recent years is in 3-D computed angiography. Here, one acquires a limited set of projection images of blood-perfused structures, such as vascular trees in the head or abdominal regions [29], =-=[34]-=-. It should be noted, however, that the objects reconstructed in 3-D computed angiography are of rather high contrast, which poses the reconstruction problem as almost a binary one. Another notable re... |

8 |
Direct Fourier methods in 3D reconstruction from cone-beam data
- Axelson
- 1994
(Show Context)
Citation Context ...artifacts. Henceforth, we will prescribe a number of modifications of ART’s projection and backproTABLE I THE DEFINITION OF A 3-D EXTENSION OF THE SHEPP–LOGAN PHANTOM [35], SIMILAR TO THE ONE USED=-= BY [6]. -=-THE ANGLES AND 0 ARE THE POLAR AND AZIMUTHAL ANGLES OF THE ELLIPSOID -AXIS. THE SCANNER ROTATES ABOUT THE ‰ -AXIS jection mechanisms with which accurate reconstructions can be obtained and which do ... |

8 |
Algebraic reconstruction can be made computationally efficient
- Herman, Meyer
- 1993
(Show Context)
Citation Context ... gives rise to a significantly slower reconstruction speed and is the main reason for ART’s under-utilization in present clinical applications. While the research on general 2-D ART [11], [12], [14]=-=, [17]-=- and 2-D fan-beam ART [13], [26] is numerous, the published literature on 3-D cone-beam reconstructors using ART-type algorithms is rather sparse. Exceptions are the early works by Colsher [9] and Sch... |

8 |
RA.The dynamic spatial reconstructor: an x-ray video-fluoroscopic CT scanner for dynamic volume imaging of moving organs
- Robb
- 1982
(Show Context)
Citation Context ...ch as the stacks-of-(fanbeam-)slices representation or the more recent helical/spiral CT (see, e.g., [10]). Likewise, imaging of dynamic structures such as the human heart is also greatly facilitated =-=[32]-=-. Cone-beam imaging received much attention with the construction of the dynamic spatial reconstructor (DSR) at the Mayo Clinic [32] for dynamic volume imaging of moving organs. However, since the con... |

6 | Image intensifier-based volume angiography imaging system: work in progress - Ning, Rooker - 1993 |

5 |
Demonstration of a software package for the reconstruction of a dynamically changing structure of the human heart from cone-beam X-ray projections
- Altschuler, Censor, et al.
- 1980
(Show Context)
Citation Context ...object in parallel layers. Later, Altschuler proposed two true cone-beam algorithms for the DSR, one using an analytic series expansion approach [1], [2] and one using an iterative Bayesian framework =-=[3]-=-. In unrelated work, Budinger [8] developed a solution based on least squares. It is apparent that over the past 15 years cone-beam researchers have mostly focused on reconstruction algorithms based o... |

4 |
Comparison of three-dimensional x-ray cone-beam reconstruction algorithms with circular source trajectories
- Rizo, Grangeat, et al.
- 1991
(Show Context)
Citation Context ... only. It is well known that this configuration gives rise to artifacts in object planes further off the midplane, due to incomplete coverage of the object’s 3-D Radon domain (see, for instance, Riz=-=o [31]-=-). Since this is an unavoidable issue with circular source orbits, our work will not eliminate these kind of artifacts. Only those artifacts related to aliasing during the ART reconstruction procedure... |

4 |
Iterative three-dimensional reconstruction from twincone beam projections
- Schlindwein
- 1978
(Show Context)
Citation Context ...fan-beam ART [13], [26] is numerous, the published literature on 3-D cone-beam reconstructors using ART-type algorithms is rather sparse. Exceptions are the early works by Colsher [9] and Schlindwein =-=[36]-=-. However, these implementations were rather inaccurate. The most popular use of 3-D cone-beam ART in recent years is in 3-D computed angiography. Here, one acquires a limited set of projection images... |

3 |
Iterative three-dimensional image reconstruction, from tomographic projections
- Colsher
- 1977
(Show Context)
Citation Context ... [14], [17] and 2-D fan-beam ART [13], [26] is numerous, the published literature on 3-D cone-beam reconstructors using ART-type algorithms is rather sparse. Exceptions are the early works by Colsher =-=[9]-=- and Schlindwein [36]. However, these implementations were rather inaccurate. The most popular use of 3-D cone-beam ART in recent years is in 3-D computed angiography. Here, one acquires a limited set... |

3 |
Multislice reconstruction from twin-cone beam scanning
- Kowalski
- 1979
(Show Context)
Citation Context ...t. This problem is generally called the teepee effect and is shown in Fig. 13(a). We can prevent the teepee effect by using noncircular source paths (see, e.g., [40]) or twin-cone source arrangements =-=[21]-=-, [36]. In the latter solution, two cone-beam sources rotate on noncoplanar coaxial congruent circular orbits around the patient and the reconstruction algorithm utilizes only rays, due to adjacent ha... |

3 |
3D grids for image reconstruction using sphericallysymmetric volume elements
- “Efficient
- 1995
(Show Context)
Citation Context ...s for SART this ratio can be brought down to 1.15. It should be mentioned that the ART methods outlined in this paper for cubic grids also fully extend to the dodecahedral grids that were proposed in =-=[24]-=-. These grids were shown to reduce the number of voxels to be processed by about 30%. We have demonstrated that our new variant of ART (and also SART) removes cone-beam related aliasing artifacts, eve... |

3 | Fast Implementations of Algebraic Methods for Three-Dimensional Reconstruction from Cone-Beam Data
- Mueller, Yagel, et al.
- 1999
(Show Context)
Citation Context ... then traverse all (scaled) voxel footprints for each projection ray and, in this way, accumulate so that the respective ray and weight sums in the nominator and denominator of (3), respectively (see =-=[27]-=- for more detail). Backprojection is performed in a similar way except that, here, the voxels receive (corrective) energy, scaled by their weight factors, instead of emitting it. The choice of varies ... |

2 |
Fully three-dimensional image reconstruction from cone-beam sources
- Altschuler, Herman, et al.
- 1978
(Show Context)
Citation Context ...g a traditional fan-beam algorithm, reconstructing the object in parallel layers. Later, Altschuler proposed two true cone-beam algorithms for the DSR, one using an analytic series expansion approach =-=[1]-=-, [2] and one using an iterative Bayesian framework [3]. In unrelated work, Budinger [8] developed a solution based on least squares. It is apparent that over the past 15 years cone-beam researchers h... |