Jan Timmer

The gridding method for image reconstruction by Fourier transformation

The authors explore a computational method for reconstructing an n-dimensional signal f from a sampled version of its Fourier transform f;. The method involves a window function w; and proceeds in three steps. First, the convolution g;=w;*f; is computed numerically on a Cartesian grid, using the available samples of f;. Then, g=wf is computed via the inverse discrete Fourier transform, and finally f is obtained as g/w. Due to the smoothing effect of the convolution, evaluating w;*f; is much less error prone than merely interpolating f;. The method was originally devised for image reconstruction in radio astronomy, but is actually applicable to a broad range of reconstructive imaging methods, including magnetic resonance imaging and computed tomography. In particular, it provides a fast and accurate alternative to the filtered backprojection. The basic method has several variants with other applications, such as the equidistant resampling of arbitrarily sampled signals or the fast computation of the Radon (Hough) transform.

The n-PI-method for helical cone-beam CT

A new class of acquisition schemes for helical cone-beam computed tomography (CB-CT) scanning is introduced, and their effect on the reconstruction methods is analyzed. These acquisition schemes are based on a new detector shape that is bounded by the helix. It will be shown that the data acquired with these schemes are compatible with exact reconstruction methods, and the adaptation of exact reconstruction algorithms to the new acquisition geometry is described. At the same time, the so-called PI-sufficiency condition is fulfilled. Moreover, a good fit to the acquisition requirements of the various medical applications of cone-beam CT is achieved. In contrast to other helical cone-beam acquisition and reconstruction methods, the n-PI-method introduced in this publication allows for variable pitches of the acquisition helix. This additional feature will introduce a higher flexibility into the acquisition protocols of future medical cone-beam scanners. An approximative n-PI-filtered backprojection (n-PI-FBP) reconstruction method is presented and verified. It yields convincing image quality.

A rational approach to dose reduction in CT: individualized scan protocols

The aim of this study was to demonstrate that dose reduction and constant image quality can be achieved by adjusting X-ray dose to patient size. To establish the relation between patient size, image quality and dose we scanned 19 patients with reduced dose. Image noise was measured. Four radiologists scored image quality subjectively, whereby a higher score meant less image quality. A reference patient diameter was determined for which the dose was just sufficient. Then 22 patients were scanned with the X-ray dose adjusted to their size. Again, image noise was measured and subjective image quality was scored. The dose reduction compared with the standard protocol was calculated. In the first group the measured noise was correlated to the patient diameter (ρ=0.78). This correlation is lost in the second group (ρ=–0.13). The correlation between patient diameter and subjective image quality scores changes from ρ=0.60 (group 1) to ρ=–0.69 (group 2). Compared with the standard protocol, the dose was reduced (mean 28%, range 0–76%) in 19 of 22 patients (86%). Dose reduction and constant noise can be achieved when the X-ray dose is adjusted to the patient diameter. With constant image noise the subjective image quality increases with larger patients.

Performance of standard fluoroscopy antiscatter grids in flat-detector-based cone-beam CT

In this paper, the performance of focused lamellar anti-scatter grids, which are currently used in fluoroscopy, is studied in order to determine guidelines of grid usage for flat detector based cone beam CT. The investigation aims at obtaining the signal to noise ratio improvement factor by the use of anti-scatter grids. First, the results of detailed Monte Carlo simulations as well as measurements are presented. From these the general characteristics of the impinging field of scattered and primary photons are derived. Phantoms modeling the head, thorax and pelvis regions have been studied for various imaging geometries with varying phantom size, cone and fan angles and patient-detector distances. Second, simulation results are shown for ideally focused and vacuum spaced grids as best case approach as well as for grids with realistic spacing materials. The grid performance is evaluated by means of the primary and scatter transmission and the signal to noise ratio improvement factor as function of imaging geometry and grid parameters. For a typical flat detector cone beam CT setup, the grid selectivity and thus the performance of anti-scatter grids is much lower compared to setups where the grid is located directly behind the irradiated object. While for small object-to-grid distances a standard grid improves the SNR, the SNR for geometries as used in flat detector based cone beam CT is deteriorated by the use of an anti-scatter grid for many application scenarios. This holds even for the pelvic region. Standard fluoroscopy anti-scatter grids were found to decrease the SNR in many application scenarios of cone beam CT due to the large patient-detector distance and have, therefore, only a limited benefit in flat detector based cone beam CT.

Artefacts in spiral-CT images and their relation to pitch and subject morphology.

Abstract. This qualitative study is intended to create awareness of artefacts that are associated with spiral-CT imaging. A simple description of spiral-CT reconstruction is used to explain how these artefacts depend on the pitch and subject morphology, and shows when these artefacts are likely to impair the diagnostic value of the acquired images. We scanned a cone and rod phantom with pitch 2, and used the acquired images to demonstrate how spiral data acquisition and interpolation leads to artefacts in the reconstructed images. We then demonstrated the effects of various pitches in scans of a human cadaver, whereas the slice thickness was kept constant. Some patient studies are presented in order to show the possible clinical consequences. Spiral acquisition may cause geometric distortions and apparent inhomogeneity of homogeneous structures. We were able to link these artefacts to the way in the acquisitions were done, and the reconstructions were performed. We have shown how these artefacts can be anticipated in clinical studies. When areas of low contrast, surrounded by hypo- or hyperdense structures, are scanned with a large pitch and viewed with a narrow window, spiral artefacts may influence the diagnostic quality of the images. These effects should be considered when choosing the pitch.

Image quality of flat-panel cone beam CT

We present results on 3D image quality in terms of spatial resolution (MTF) and low contrast detectability, obtained on a flat dynamic X-ray detector (FD) based cone-beam CT (CB-CT) setup. Experiments have been performed on a high precision bench-top system with rotating object table, fixed X-ray tube and 176 x 176 mm2 active detector area (Trixell Pixium 4800). Several objects, including CT performance-, MTF- and pelvis phantoms, have been scanned under various conditions, including a high dose setup in order to explore the 3D performance limits. Under these optimal conditions, the system is capable of resolving less than 1% (~10 HU) contrast in a water background. Within a pelvis phantom, even inserts of muscle and fat equivalent are clearly distinguishable. This also holds for fast acquisitions of up to 40 fps. Focusing on the spatial resolution, we obtain an almost isotropic three-dimensional resolution of up to 30 lp/cm at 10% modulation.

Multiscale denoising and enhancement of 3D rotational X-ray imaging for percutaneous vertebroplasty

3D rotational X-ray (3DRX) imaging has been suggested as an alternative to CT scans for evaluating percutaneous vertebroplasty procedures. We present in this study preliminary work on denoising and enhancement of 3DRX data sets with dyadic wavelet thresholding. Thresholding operators were tuned to accommodate for spatial variability of the anatomical features and the noise energy at each level of decomposition. Results are presented for three clinical data sets with quantitative measurements of contrast to noise ratio (CNR) and qualitative evaluation with three- dimensional rendering of vertebral bodies having cement fillings. Wavelet denoising demonstrated the efficient on enhancement of subtle vertebral edges and clearer delineation of the cement filling contours which improves the accuracy in clinical practice and helps to make 3DRX more competitive when compared to CT.