Biomedical subjects

Willi A Kalender

Publications and source records attributed to Willi A Kalender.

At least 19 recordsLinked to original sources

PubMed2006

Combining deterministic and Monte Carlo calculations for fast estimation of scatter intensities in CT.

A side effect of increased volume coverage by using multi-row and flat-panel detectors in computed tomography (CT) is the concurrently growing contribution of scattered radiation to the measured signal. In order to investigate the effect of scatter on x-ray projections used for CT imaging, our study aimed at the development of a simulation tool for fast calculation of primary and scatter intensities. We developed a deterministic method to assess the contribution of single-scatter events to the measured signal. The investigation of multiple scatter by Monte Carlo simulations showed that it results in a smooth signal as compared to single scatter. A hybrid method is proposed in order to optimize the performance of the scatter simulation: a fast and exact analytical calculation of the single-scatter intensity combined with a coarse Monte Carlo (MC) estimate of multiple scatter to reduce overall computational expenses, while assuring an acceptable signal quality. The results of the hybrid simulation of total scatter were in excellent agreement with the corresponding MC only simulations, thereby allowing us to reduce computational time by orders of magnitude. Estimates of two-dimensional scatter distributions for flat-panel CT imaging took about 30-40 s (per projection). The hybrid method provides a realistic simulation of x-ray scatter and offers a basis for scatter correction approaches.

Algorithms↗

PubMed2006

A hierarchical 3D segmentation method and the definition of vertebral body coordinate systems for QCT of the lumbar spine.

We have developed a new hierarchical 3D technique to segment the vertebral bodies in order to measure bone mineral density (BMD) with high trueness and precision in volumetric CT datasets. The hierarchical approach starts with a coarse separation of the individual vertebrae, applies a variety of techniques to segment the vertebral bodies with increasing detail and ends with the definition of an anatomic coordinate system for each vertebral body, relative to which up to 41 trabecular and cortical volumes of interest are positioned. In a pre-segmentation step constraints consisting of Boolean combinations of simple geometric shapes are determined that enclose each individual vertebral body. Bound by these constraints viscous deformable models are used to segment the main shape of the vertebral bodies. Volume growing and morphological operations then capture the fine details of the bone-soft tissue interface. In the volumes of interest bone mineral density and content are determined. In addition, in the segmented vertebral bodies geometric parameters such as volume or the length of the main axes of inertia can be measured. Intra- and inter-operator precision errors of the segmentation procedure were analyzed using existing clinical patient datasets. Results for segmented volume, BMD, and coordinate system position were below 2.0%, 0.6%, and 0.7%, respectively. Trueness was analyzed using phantom scans. The bias of the segmented volume was below 4%; for BMD it was below 1.5%. The long-term goal of this work is improved fracture prediction and patient monitoring in the field of osteoporosis. A true 3D segmentation also enables an accurate measurement of geometrical parameters that may augment the clinical value of a pure BMD analysis.

Absorptiometry, Photon↗

PubMed2006

X-ray computed tomography.

X-ray computed tomography (CT), introduced into clinical practice in 1972, was the first of the modern slice-imaging modalities. To reconstruct images mathematically from measured data and to display and to archive them in digital form was a novelty then and is commonplace today. CT has shown a steady upward trend with respect to technology, performance and clinical use independent of predictions and expert assessments which forecast in the 1980s that it would be completely replaced by magnetic resonance imaging. CT not only survived but exhibited a true renaissance due to the introduction of spiral scanning which meant the transition from slice-by-slice imaging to true volume imaging. Complemented by the introduction of array detector technology in the 1990s, CT today allows imaging of whole organs or the whole body in 5 to 20 s with sub-millimetre isotropic resolution. This review of CT will proceed in chronological order focussing on technology, image quality and clinical applications. In its final part it will also briefly allude to novel uses of CT such as dual-source CT, C-arm flat-panel-detector CT and micro-CT. At present CT possibly exhibits a higher innovation rate than ever before. In consequence the topical and most recent developments will receive the greatest attention.

Biotechnology↗

PubMed2006

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↗

PubMed2006

Contrast-enhanced coronary artery visualization by dual-source computed tomography--initial experience.

UNLABELLED: Multi-detector computed tomography (CT) scanners, by virtue of their high temporal and spatial resolution, permit imaging of the coronary arteries. However, motion artifacts, especially in patients with higher heart rates, can impair image quality. We thus evaluated the performance of a new dual-source CT (DSCT) with a heart rate independent temporal resolution of 83 ms for the visualization of the coronary arteries in 14 consecutive patients. METHODS: Fourteen patients (mean age 61 years, mean heart rate 71 min(-1)) were studied by DSCT. The system combines two arrays of an X-ray tube plus detector (64 slices) mounted on a single gantry at an angle of 90 degrees With a rotation speed of 330 ms, a temporal resolution of 83 ms (one-quarter rotation) can be achieved independent of heart rate. For data acquisition, intraveous contrast agent was injected at a rate of 5 ml/s. Images were reconstructed with 0.75 slice thickness and 0.5 mm increment. The data sets were evaluated concerning visibility of the coronary arteries and occurrence of motion artifact. RESULTS: Visualization of the coronary arteries was successful in all patients. Most frequently, image reconstruction at 70% of the cardiac cycle provided for optimal image quality (50% of patients). Of a total of 226 coronary artery segments, 222 (98%) were visualized free of motion artifact. In summary, DSCT constitutes a promising new concept for cardiac CT. High and heart rate independent temporal resolution permits imaging of the coronary arteries without motion artifacts in a substantially increased number of patients as compared to earlier scanner generations. Larger and appropriately designed studies will need to determine the method's accuracy for detection of coronary artery stenoses.

Artifacts↗

PubMed2006

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↗

PubMed2006

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↗

PubMed2005

Bone status in elite male runners.

The aim of our study was to compare long distance runners to body mass index (BMI)- and age-matched healthy controls with respect to bone parameters at all relevant loaded and nonloaded skeletal sites. Furthermore, we assessed the effect of running volume on bone parameters. Twenty elite male runners (21.1 km<1:15 h; volume >75 km/week/year) participated in the study (RG), 11 age- and BMI-matched male subjects (28+/-5 years) served as nontraining controls (CG). Subjects with any medication or illness affecting bone metabolism or with a family history of osteoporosis were not included. Bone parameters at various sites (total body, lumbar spine, femoral neck/hip, calcaneus) were measured by dual energy X-ray (DXA), quantitative computed tomography and quantitative ultrasound. Body composition was assessed via DXA and bioimpedance analysis; nutritional parameters were determined by 5-day dietary protocols. Training variables were assessed by questionnaires. Compared with nontraining controls runners had significantly higher BMD at all loaded sites (calcaneus, lower limbs, femoral neck, pelvis, and trabecular lumbar spine). BMD at nonloaded sites (ribs, upper limbs, and skull) was slightly but not significantly higher in the runners. We observed a low (r=0.30), nonsignificant association between training volume (km/week/year) and trabecular BMD of the femoral neck, which disappeared after adjusting for age, BMI, and body fat in this group of highly trained male runners. The effect of long distance running per se on bone parameters is not deleterious.

Adult↗

PubMed2005

An anatomic coordinate system of the femoral neck for highly reproducible BMD measurements using 3D QCT.

In this paper, a procedure for the determination of an anatomically oriented coordinate system of the femoral neck (NCS) in 3D spiral CT datasets is described. The origin of the NCS is centered in the minimal cross-sectional area of the neck. Its three axes are defined as follows: the so called neck axis is perpendicular to this area and points towards the femoral head, the second axis is the principal axis of the minimal cross-sectional area and the third axis is perpendicular to the other two. After a semi-automatic 3D segmentation of the proximal femur the NCS is automatically determined in a two-step minimization procedure. Relative to the coordinate system volumes of interest (VOIs) are positioned in which bone mineral density (BMD) and cortical thickness are analyzed. We investigated intra- and inter-operator precision of the position of the NCS, the BMD in cortical and trabecular VOIs, and cortical thickness in nine pelvic CT datasets obtained from clinical routine examinations. We further investigated the effect of increased noise by adding Gaussian distributed noise to measured projections before tomographic reconstruction. The mean precision error (averaged form the results of the nine datasets) of the NCS position was less than 0.5 mm and smaller than 2.25 degrees . There were no significant differences between inter- and intra-operator analyses. Precision errors in trabecular BMD were smaller than 3% in a stack of five 1 mm thin slices cut perpendicularly to the neck axis and smaller than 1% in a spherical VOI encompassing the neck. Relative precision errors for cortical BMD were smaller than 3% for both VOIs. An increase of noise up to a factor of 5 caused a maximal displacement of the NCS origin position by less than 1mm and a rotation by less than 2 degrees .

Bone Density↗

PubMed2005

A navigation system for minimally invasive CT-guided interventions.

The purpose of our project was to develop a novel navigation system for interventional radiology. Fields of application are minimally invasive percutaneous interventions performed under local anaesthesia. In order to reduce unintentional patient movements we used a patient vacuum immobilization device. Together with the vacuum fixation and a newly developed reference frame we achieved a fully automatic patient-to-image registration independent from the tracking system. The combination of the software and a novel designed needle holder allows for an adjustment of the needle within a few seconds. The complete system is adapted to the requirements of the radiologist and to the clinical work-flow. For evaluation of the navigation system we performed a phantom study with a perspex phantom and achieved an average needle positioning accuracy of less than 0.7 mm.

Algorithms↗

PubMed2005

Statistical face models for the rediction of soft-tissue deformations after orthognathic osteotomies.

This paper describes a technique to approximately predict the facial morphology after standardized orthognathic ostoetomies. The technique only relies on the outer facial morphology represented as a set of surface points and does not require computed tomography (CT) images as input. Surface points may either be taken from 3D surface scans or from 3D positions palpated on the face using a tracking system. The method is based on a statistical model generated from a set of pre- and postoperative 3D surface scans of patients that underwent the same standardized surgery. The model contains both the variability of preoperative facial morphologies and the corresponding postoperative deformations. After fitting the preoperative part to 3D data from a new patient the preoperative face is approximated by the model and the preiction of the postoperative morphology can be extracted at the same time. We built a model based on a set of 15 patient data sets and tested the predictive power in leave-one-out tests for a set of relevant cephalometric landmarks. The average prediction error was found to be between 0.3 and 1.2 mm at all important facial landmarks in the relevant areas of upper and lower jaw. Thus the technique provides an easy and powerful way of prediction which avoids time, cost and radiation required by other prediction techniques such as those based on CT scans.

Computer Simulation↗

PubMed2005

CT: the unexpected evolution of an imaging modality.

The evolution of CT can be characterized remarkably well by three phases of developments assigned to the 1970s, 1980s and 1990s. The first decade saw rapid innovations followed by a phase of consolidation and slow growth in the 1980s, and finally a third phase of very rapid development including spiral CT and multirow detectors in the 1990s. The evolution of technical developments during these first three decades is briefly reviewed. CT has reached a very high degree of maturity, mastering almost all clinical demands. The focus of this review is set on potential future developments and trends. Further doubling or multiplying of the numbers of slices acquired simultaneously cannot be expected to provide further essential innovations. New paradigms are required to advance the field. Future potential developments are outlined, including multisource, multidetector scanners for cardiac and dual-energy CT, new detector technologies, approaches to data handling and dose management.

Equipment Design↗

PubMed2005

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↗

PubMed2004

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↗

PubMed2004

A pragmatic approach to metal artifact reduction in CT: merging of metal artifact reduced images.

The purpose of this study was to improve metal artifact reduction (MAR) in X-ray computed tomography (CT) by the combination of two artifact reduction methods. The presented method constitutes an image-based weighted superposition of images processed with two known methods for MAR: linear interpolation of reprojected metal traces (LI) and multi-dimensional adaptive filtering of the raw data (MAF). Two weighting concepts were realized that take into account mean distances of image points from metal objects or additional directional components. Artifact reduction on patient data from the jaw and the hip region shows that although the application of only one of the MAR algorithms can already improve image quality, these methods have specific drawbacks. While MAF does not correct corrupted CT values, LI often introduces secondary artifacts. The corrective impact of the merging algorithm is almost always superior to the application of only one of the methods. The results obtained with directional weighting are equal to or in many cases better than those of the distance weighting scheme. Merging combines the advantages of two fundamentally different approaches to artifact reduction and can improve the quality of images that are affected by metal artifacts.

Algorithms↗

PubMed2004

Interactive 3D editing tools for image segmentation.

Segmentation is an important part of image processing, which often has a large impact on quantitative image analysis results. Fully automated operator independent segmentation procedures that successfully work in a population with a larger biological variation are extremely difficult to design and usually some kind of operator intervention is required, at least in pathological cases. We developed a variety of 3D editing tools that can be used to correct or improve results of initial automatic segmentation procedures. Specifically we will discuss and show examples for three types of editing tools that we termed: hole-filling (tool 1), point-bridging (tool 2), and surface-dragging (tool 3). Each tool comes in a number of flavors, all of which are implemented in a truly 3D manner. We describe the principles, evaluate efficiency and flexibility, and discuss advantages and disadvantages of each tool. We further demonstrate the superiority of the 3D approach over the time-consuming slice-by-slice editing of 3D datasets, which is still widely used in medical image processing today. We conclude that performance criteria for automatic segmentation algorithms may be eased significantly by including 3D editing tools early in the design process.

Algorithms↗

PubMed2004

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↗

PubMed2004

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↗