PubMed Health⌕ Search

Biomedical subjects

Marc Kachelriess

Publications and source records attributed to Marc Kachelriess.

11 recordsLinked to original sources

Impact of the z-flying focal spot on resolution and artifact behavior for a 64-slice spiral CT scanner.

The effect of the z-flying focal spot (zFFS) technology was evaluated by simulations and measurements with respect to resolution and artifact behavior for a 64-slice spiral cone-beam computed tomography (CT) scanner. The zFFS alternates between two z-positions of the X-ray focal spot, acquiring two slices per detector row, which results in double sampling in the z-direction. We implemented a modified reconstruction that is able to obtain images as they would be without zFFS. A delta phantom equipped with a thin gold disc was used to measure slice sensitivity profiles (SSP), and a high-contrast bar phantom was used to quantify the resolution in the x/z-plane with and without zFFS. The zFFS decreases the full width at half maximum (FWHM) of the SSPs by a factor of about 1.4. The double z-sampling allows the separation of 0.4 mm bars in the z-direction compared with 0.6 mm in the case without zFFS. The zFFS effectively reduces windmill artifacts in the reconstructed images while maintaining the transverse resolution, even at the largest available pitch value of 1.5.

Artifacts↗

Empirical cupping correction: a first-order raw data precorrection for cone-beam computed tomography.

We propose an empirical cupping correction (ECC) algorithm to correct for CT cupping artifacts that are induced by nonlinearities in the projection data. The method is raw data based, empirical, and requires neither knowledge of the x-ray spectrum nor of the attenuation coefficients. It aims at linearizing the attenuation data using a precorrection function of polynomial form. The coefficients of the polynomial are determined once using a calibration scan of a homogeneous phantom. Computing the coefficients is done in image domain by fitting a series of basis images to a template image. The template image is obtained directly from the uncorrected phantom image and no assumptions on the phantom size or of its positioning are made. Raw data are precorrected by passing them through the once-determined polynomial. As an example we demonstrate how ECC can be used to perform water precorrection for an in vivo micro-CT scanner (TomoScope 30 s, VAMP GmbH, Erlangen, Germany). For this particular case, practical considerations regarding the definition of the template image are given. ECC strives to remove the cupping artifacts and to obtain well-calibrated CT values. Although ECC is a first-order correction and cannot compete with iterative higher-order beam hardening or scatter correction algorithms, our in vivo mouse images show a significant reduction of bone-induced artifacts as well. A combination of ECC with analytical techniques yielding a hybrid cupping correction method is possible and allows for channel-dependent correction functions.

Algorithms↗

Multithreaded cardiac CT.

Phase-correlated CT, as it is used for cardiac imaging, is the most popular and the most important but also the most demanding special CT application in the clinical routine, today. Basically, it fulfills the four-dimensional imaging task of depicting a quasiperiodically moving object at any desired motion phase with significantly reduced motion artifacts. Although image quality with phase-correlated reconstruction is far better than with standard reconstruction, there are motion artifacts remaining and improvements of temporal resolution are required. As a well-known alternative to simply decreasing rotation time, we consider a spiral cone-beam CT scanner that has G x-ray guns and detectors mounted. We call this a multisource or a multithreaded CT scanner. Aiming for improved temporal resolution the relative temporal resolution tau, which measures the fraction of a motion period that enters the image, is studied as a function of the motion rate (heart rate) and the degree of scan overlap (pitch value) for various configurations. The parameters to optimize are the number of threads G and the interthread parameters delta alpha and delta z, which are the angular and the longitudinal separation between adjacent threads, respectively. To demonstrate the improvements approximate image reconstruction of multithreaded raw data is performed by using a generalization of the extended parallel back projection cone-beam reconstruction algorithm [Med. Phys. 31(6), 1623-1641 (2004)] to the case of multithreaded CT. Reconstructions of a simulated cardiac motion phantom and of simulated semi-antropomorphic phantoms are presented for two and three threads and compared to the single-threaded case to demonstrate the potential of multithreaded cardiac CT. Patient data were acquired using a clinical double-threaded CT scanner to validate the theoretical results. The optimum angle delta alpha between the tubes is 90 degrees for a double-threaded system, and for triple-threaded scanners it is 60 degrees or 120 degrees. In all cases, delta z = 0 results as an optimum, which means that the threads should be mounted in the same transversal plane. However, the dependency of the temporal resolution on delta z is very weak and a longitudinal separation delta z not = 0 would not deteriorate image quality. The mean temporal resolution achievable with an optimized multithreaded CT scanner is a factor of G better than the mean temporal resolution obtained with a single-threaded scanner. The standard reconstructions showed decreased cone-beam artifacts with multithreaded CT compared to the single-threaded case. Our phase-correlated reconstructions demonstrate that temporal resolution is significantly improved with multithreaded CT. The clinical patient data confirm our results.

Algorithms↗

Presampling, algorithm factors, and noise: considerations for CT in particular and for medical imaging in general.

CT scanners acquire noisy data at discrete sample positions. Typically, a convention of how to continue these data from discrete integer positions to the continuous domain must be applied during processing. We study the properties of three typical one-dimensional spatial domain interpolation algorithms in terms of a cost or quality factor Q. This figure of merit Q is a function of spatial resolution, data noise, and dose and is used to optimize detector design. Spatial resolution R is defined as either mean square width delta or as the full width at half maximum W of the point spread function (PSF). Our results show that a trapezoidal interpolation algorithm is optimal for the high resolution domain (relative to the detector aperture size g) and should be replaced by a triangular or Gaussian interpolation function for spatial resolutions of about 1.3g or larger; these result in bell-shaped PSFs. Assuming such a hybrid algorithm we find a 1.5-fold increase of Q2-this is equivalent to 50% improved dose usage-when smoothing the data to a spatial resolution of 3g or more compared to a highest resolution reconstruction. Therefore it is advisable to use detectors of one-third of the size of the desired spatial resolution W and to compensate for the 1.5-fold increase in Q2 by reducing dose by 33%. Under the presence of moderately sized septa (e.g., 10% of the spatial resolution element size) the benefit of optimizing still lies in the order of 30% improved dose usage; in that case the detector size g should be on the order of W/2 and a dose reduction of 23% can be achieved. Again, bell-shaped PSFs show a better tradeoff between noise and resolution for a given dose than rectangular-shaped PSFs. The general interpretation of our results is that the degree of freedom of choosing the weighting or interpolation function for a given resolution is large for small detectors and small for large detectors. Thus systems with small g have a higher potential of optimization compared to systems with large g. Similarly, detector binning, which corresponds to replacing g by 2g, should be avoided. Note that the figures reported correspond to a one-dimensional interpolation. Two-dimensional detectors typically separate and resulting quality factors can be easily obtained by multiplication. Then, Q2 is expected to improve by a factor of 1.52 without septa and by a factor of 1.32 with septa. This indicates that dose can be reduced by about 56% and about 41%, respectively. Our findings are general and not restricted to CT. They can be readily applied to medical or nonmedical imaging devices and digital detectors and they may also turn out to be useful in other fields.

Algorithms↗

Improvement of image quality of multislice spiral CT scans of the head and neck region using a raw data-based multidimensional adaptive filtering (MAF) technique.

The purpose was to evaluate the potential of the multidimensional adaptive filtering (MAF) technique by investigating its effects on image noise and image quality in multislice spiral CT (MSCT) examinations of the head and neck region. Fifty patients with head and neck tumors were examined using MSCT with a high resolution protocol. Reconstructions were performed using dedicated reconstruction software with a standard algorithm both without and with MAF using different modification. In all reconstructions, we measured the noise in seven different anatomical structures. The image quality and image noise were rated on a five-point scale. There was a significant (P<0.05) reduction in mean pixel noise in the reconstructions using MAF in comparison to the standard reconstructions, but there was no significant difference between the different modification fractions. With MAF the mean reduction in noise level was 60%, depending upon body shape and anatomical region. Independently from the used modification fraction, MAF led to a significant (P<0.05) improvement of image quality. In direct comparison of the different filter strength, the optimal image quality was achieved in the investigations with 15% MAF. The use of MAF facilitates the distinction of anatomical and pathological structures from artifacts in the supraclavicular fossae and the upper mediastinum, whereas the image quality of the upper portions of the neck remained unchanged. MAF improved image quality by reducing the noise level and removing noise structures without loss of image sharpness. This technique offers new perspectives to reduce the patient dose.

Adipose Tissue↗

Extended parallel backprojection for standard three-dimensional and phase-correlated four-dimensional axial and spiral cone-beam CT with arbitrary pitch, arbitrary cone-angle, and 100% dose usage.

We have developed a new approximate Feldkamp-type algorithm that we call the extended parallel backprojection (EPBP). Its main features are a phase-weighted backprojection and a voxel-by-voxel 180 degrees normalization. The first feature ensures three-dimensional (3-D) and 4-D capabilities with one and the same algorithm; the second ensures 100% detector usage (each ray is accounted for). The algorithm was evaluated using simulated data of a thorax phantom and a cardiac motion phantom for scanners with up to 256 slices. Axial (circle and sequence) and spiral scan trajectories were investigated. The standard reconstructions (EPBPStd) are of high quality, even for as many as 256 slices. The cardiac reconstructions (EPBPCI) are of high quality as well and show no significant deterioration of objects even far off the center of rotation. Since EPBPCI uses the cardio interpolation (CI) phase weighting the temporal resolution is equivalent to that of the well-established single-slice and multislice cardiac approaches 180 degrees CI, 180 degrees MCI, and ASSRCI, respectively, and lies in the order of 50 to 100 ms for rotation times between 0.4 and 0.5 s. EPBP appears to fulfill all required demands. Especially the phase-correlated EPBP reconstruction of cardiac multiple circle scan data is of high interest, e.g., for dynamic perfusion studies of the heart.

Algorithms↗

Geometric misalignment and calibration in cone-beam tomography.

We present a new high-precision method for the geometric calibration in cone-beam computed tomography. It is based on a Fourier analysis of the projection-orbit data, recorded with a flat-panel area detector, of individual point-like objects. For circular scan trajectories the complete set of misalignment parameters which determine the deviation of the detector alignment from the ideal scan geometry are obtained from explicit analytic expressions. To derive these expressions we show how to disentangle the problems of calculating misalignment parameters and point coordinates. The calculation of the coordinates of the point objects inside the scanned volume, in units of the distance from the focal spot to the center of rotation, is then possible analytically likewise. We simulate point-projection data on a misaligned detector with various amounts of randomness added to mimic measurement uncertainties. This data is then employed in our calibration to validate the method by comparing the resulting misalignment parameters and point coordinates to the known true ones. We also present our implementation and results for the geometric calibration of micro-CT systems. The effectiveness of the corresponding misalignment correction in reducing image artifacts is exemplified by reconstructed micro-CT images.

Algorithms↗

Comparison of image quality in contrast-enhanced coronary-artery visualization by electron beam tomography and retrospectively electrocardiogram-gated multislice spiral computed tomography.

RATIONALE AND OBJECTIVES: To compare the image quality of electron beam tomography (EBT) and multislice spiral CT (MSCT) for coronary artery visualization. MATERIALS AND METHODS: Two groups of 30 patients without coronary stenoses were studied by MSCT (4 x 1 mm collimation) or EBT (3 mm slice thickness). Contrast-to-noise ratio (CNR), overall length of the visualized arteries and vessel length free of motion artifacts were measured. RESULTS: Length of visualized arteries was equal in MSCT and EBT. In EBT, longer segments were depicted free of motion artifacts (MSCT: 73%, EBT: 92% of visualized length, P< 0.001) and CNR was significantly higher than in MSCT (15.4 vs. 9.0; P< 0.001). In both modalities, vessel diameters correlated closely to quantitative coronary angiography. CONCLUSIONS: EBT and MSCT permit reliable coronary artery visualization and measurement of vessel diameters. For the used scan protocol, MSCT images had a lower CNR and were more frequently affected by motion.

Aged↗

Kymogram detection and kymogram-correlated image reconstruction from subsecond spiral computed tomography scans of the heart.

Subsecond single-slice, multi-slice or cone-beam spiral computed tomography (SSCT, MSCT, CBCT) offer great potential for improving heart imaging. Together with the newly developed phase-correlated cardiac reconstruction algorithms 180 degrees MCD and 180 degrees MCI [Med. Phys. 27, 1881-1902 (2000)] or related algorithms provided by the CT manufacturers, high image quality can be achieved. These algorithms require information about the cardiac motion, i.e., typically the simultaneously recorded electrocardiogram (ECG), to synchronize the reconstruction with the cardiac motion. Neither data acquired without ECG information (standard patients) nor acquisitions with corrupted ECG information can be handled adequately. We developed a method to extract the appropriate information about cardiac motion directly from the measured raw data (projection data). The so-called kymogram function is a measure of the cardiac motion as a function of time t or as a function of the projection angle alpha. In contrast to the ECG which is a global measure of the heart's electric excitation, the kymogram is a local measure of the heart motion at the z-position z(a) at projection angle a. The patient's local heart rate as well as the necessary synchronization information to be used with phase-correlated algorithms can be extracted from the kymogram by using a series of signal processing steps. The kymogram information is shown to be adequate to substitute the ECG information. Computer simulations with simulated ECG and patient measurements with simultaneously acquired ECG were carried out for a multislice scanner providing M = 4 slices to evaluate these new approaches. Both the ECG function and the kymogram function were used for reconstruction. Both were highly correlated regarding the periodicity information used for reconstruction. In 21 out of 25 consecutive cases the kymogram approach was equivalent to the ECG-correlated reconstruction; only minor differences in image quality between both methods were observed. For one patient the synchronization information detected by the ECG monitor turned out to be wrong; here, the kymogram constituted the only approach that provided useful reconstructions. Patient studies with 12 and 16 slices indicate the usefulness of our approach for cone-beam CT scans. Kymogram-correlated reconstructions also appear to have the potential to improve imaging of pericardial lung areas in general.

Algorithms↗

An extensible software-based platform for reconstruction and evaluation of CT images.

A collection of software applications dedicated to clinical x-ray computed tomography (CT) has been developed. These tools are designed for routine and scientific work and represent a complete platform, which is called syngo Explorer. The user can reconstruct, process, and view CT images in a personal computer environment independently of specialized hardware. To manage patient data and to allow full database access, the system uses the syngo platform (Siemens Medical Solutions, Erlangen, Germany), which serves as a framework for the management and visualization of Digital Imaging and Communications in Medicine (DICOM) images. Part of the philosophy behind syngo Explorer is to quickly provide specialized solutions and to be able to respond to user requests at once. Thus, syngo Explorer was built by choosing a modular setup that comprises standard and specific reconstruction algorithms as well as various artifact reduction techniques; it also helps one attain better insights into dose issues. Various raw data formats of commercial CT scanners are supported. Reconstruction algorithms equivalent to those on the manufacturer's scanner console are provided. A non-syngo Windows XP (Microsoft, Redmond, Wash) version of the application is available.

Algorithms↗