Thomas Mertelmeier

Optimizing filtered backprojection reconstruction for a breast tomosynthesis prototype device

Digital breast tomosynthesis is a new technique intended to overcome the limitations of conventional projection mammography by reconstructing slices through the breast from projection views acquired from different angles with respect to the breast. We formulate a general theory of filtered backprojection reconstruction for linear tomosynthesis. The filtering step consists of an MTF inversion filter, a spectral filter, and a slice thickness filter. In this paper the method is applied first to simulated data to understand the basic effects of the various filtering steps. We then demonstrate the impact of the filter functions with simulated projections and with clinical data acquired with a research breast tomosynthesis system.** With this reconstruction method the image quality can be controlled regarding noise and spatial resolution. In a wide range of spatial frequencies the slice thickness can be kept constant and artifacts caused by the incompleteness of the data can be suppressed.

X-ray spectrum optimization of full-field digital mammography : Simulation and phantom study

In contrast to conventional analog screen-film mammography new flat detectors have a high dynamic range and a linear characteristic curve. Hence, the radiographic technique can be optimized independently of the receptor exposure. It can be exclusively focused on the improvement of the image quality and the reduction of the patient dose. In this paper we measure the image quality by a physical quantity, the signal difference-to-noise ratio (SDNR), and the patient risk by the average glandular dose (AGD). Using these quantities, we compare the following different setups through simulations and phantom studies regarding the detection of microcalcifications and tumors for different breast thicknesses and breast compositions: Monochromatic radiation, three different anode/filter combinations: Molybdenum/molybdenum (Mo/Mo), molybdenum/rhodium (Mo/Rh), and tungsten/rhodium (W/Rh), different filter thicknesses, use of anti-scatter grids, and different tube voltages. For a digital mammography system based on an amorphous selenium detector it turned out that, first, the W/Rh combination is the best choice for all detection tasks studied. Second, monochromatic radiation can further reduce the AGD by a factor of up to 2.3, maintaining the image quality in comparison with a real polychromatic spectrum of an x-ray tube. And, third, the use of an anti-scatter grid is only advantageous for breast thicknesses larger than approximately 5 cm.

Experimental validation of a three-dimensional linear system model for breast tomosynthesis

A three-dimensional (3D) linear model for digital breast tomosynthesis (DBT) was developed to investigate the effects of different imaging system parameters on the reconstructed image quality. In the present work, experimental validation of the model was performed on a prototype DBT system equipped with an amorphous selenium (a-Se) digital mammography detector and filtered back-projection (FBP) reconstruction methods. The detector can be operated in either full resolution with 85 microm pixel size or 2 x 1 pixel binning mode to reduce acquisition time. Twenty-five projection images were acquired with a nominal angular range of +/- 20 degrees. The images were reconstructed using a slice thickness of 1 mm with 0.085 x 0.085 mm in-plane pixel dimension. The imaging performance was characterized by spatial frequency-dependent parameters including a 3D noise power spectrum (NPS) and in-plane modulation transfer function (MTF). Scatter-free uniform x-ray images were acquired at four different exposure levels for noise analysis. An aluminum (Al) edge phantom with 0.2 mm thickness was imaged to measure the in-plane presampling MTF. The measured in-plane MTF and 3D NPS were both in good agreement with the model. The dependence of DBT image quality on reconstruction filters was investigated. It was found that the slice thickness (ST) filter, a Hanning window to limit the high-frequency components in the slice thickness direction, reduces noise aliasing and improves 3D DQE. An ACR phantom was imaged to investigate the effects of angular range and detector operational modes on reconstructed image quality. It was found that increasing the angular range improves the MTF at low frequencies, resulting in better detection of large-area, low-contrast mass lesions in the phantom. There is a trade-off between noise and resolution for pixel binning and full resolution modes, and the choice of detector mode will depend on radiation dose and the targeted lesion.

Improving 3D image quality of x-ray C-arm imaging systems by using properly designed pose determination systems for calibrating the projection geometry

C-arm volume reconstruction has become increasingly popular over the last years. These imaging systems generate 3D data sets for various interventional procedures such as endovascular treatment of aneurysms or orthopedic applications. Due to their open design and mechanical instability, C-arm imaging systems acquire projections along non-ideal scan trajectories. Volume reconstruction from filtered 2D X-ray projections requires a very precise knowledge of the imaging geometry. We show that the 3D image quality of C-arm cone beam imaging devices can be improved by proper design of the calibration phantom.

Why and how is soft copy reading possible in clinical practice

The properties of the human visual system (HVS) relevant to the diagnostic process are described after a brief introduction on the general problems and advantages of using soft copy for primary radiology interpretations. At various spatial and temporal frequencies the contrast sensitivity defines the spatial resolution of the eye-brain system and the sensitivity to flicker. The adaptation to the displayed radiological scene and the ambient illumination determine the dynamic range for the operation of the HVS. Although image display devices are determined mainly by state-of-the-art technology, analysis of the HVS may suggest technical characteristics for electronic displays that will help to optimize the display to the operation of the HVS. These include display size, spatial resolution, contrast resolution, luminance range, and noise, from which further consequences for the technical components of a monitor follow. It is emphasized that routine monitor quality control must be available in clinical practice. These image quality measures must be simple enough to be applied as part of the daily routine. These test instructions might also serve as elements of technical acceptance and constancy tests.

Optimization of detector operation and imaging geometry for breast tomosynthesis

In breast tomosynthesis there are tradeoffs between resolution, noise and acquisition speed for a given glandular dose. The purpose of the present work is to investigate the dependence of tomosynthesis imaging performance on system configuration, which includes detector operational modes and image acquisition geometry. A prototype Siemens breast tomosynthesis system with maximum angular range of +/- 25 degrees was used in our investigation. The system was equipped with an amorphous selenium (a-Se) full field digital mammography detector with pixel size of 85µm. The detector can be read out with full resolution or 2x1 binning (binning in the tube travel direction), which increases the image readout rate and decreases the degradation effect of electronic noise. The total number of views can be varied from 11 to 49, and filtered back projection (FBP) method was used to reconstruct the tomosynthesis images. We investigated the effects of detector operational modes (binning) and imaging geometry (view angle and number) on temporal performance and spatial resolution of the projection images. The focal spot blur due to continuous tube travel was measured for different acquisition geometry, and its effect on in-plane presampling modulation transfer function (MTF) was compared to that due to pixel binning. A three-dimensional cascaded linear system model was developed for tomosynthesis to predict the 3D MTF, NPS and DQE. The results were compared with experimental measurements, and reasonable agreement was achieved. The understanding of the relationship between the 3D and projection image quality will lead to optimization of the x-ray spectrum, imaging geometry and reconstruction filters for digital breast tomosynthesis.

High-temporal-resolution volume heart imaging with multirow computed tomography

Functional cardiac imaging with 3rd generation CT scanners is challenging, because the temporal resolution seems to be limited to approximately 2/3 of the rotation time of the gantry. We propose a new method for high temporal resolution volume heart imaging with multirow detectors based on a retrospective electrocardiogram-gated rebinning procedure. The limited time resolution is overcome using time consistent projection data retrieved from more than one cardiac cycle. In principle the method provides volume heart imaging with adjustable time resolution at arbitrary cardiac phases. It can be applied both for spiral and axial scan imaging. The presented study is based on computer simulations incorporating a model of the human heart taking into account anatomy, motion and heart rate variability. For multirow detectors we were able to show that good image quality can be obtained even during systole with temporal resolution which even exceeds that provided by an Electron Beam Scanner in standard mode of operation. Using an area detector with detector height > 3 cm (center of rotation) the total measurement time is within one breathhold for complete volume imaging of the heart. Furthermore, freezing motion of the coronary arteries during enddiastole allows high quality 3D display of coronary anatomy.

A Novel Approach for Filtered Backprojection in Tomosynthesis Based on Filter Kernels Determined by Iterative Reconstruction Techniques

Breast tomosynthesis is a new 3D imaging technique providing 3D slices of the breast. The computation of 3D slices from a set of 2D projections is performed with a reconstruction algorithm --- either a filtered backprojection (FBP) or an algebraic iterative reconstruction method. Both approaches yield different image characteristics of the reconstructed object. The algebraic method has the main disadvantage of a very long processing time. We experimentally developed a set of filter kernels for FBP, determined by iterative reconstruction providing similar image characteristics and quality as an algebraic reconstruction. Additionally we showed that it is possible to approximate these kernels by a polynomial function of order 4. Using clinical data sets we demonstrate how the image quality and the image impression can be varied by using different reconstruction methods.

Optimization of Tomosynthesis Acquisition Parameters: Angular Range and Number of Projections

In breast tomosynthesis several projection views of the object are acquired over a limited angular range from which slice images are reconstructed. The image quality of the reconstructed slices, however, depends on the geometric acquisitions parameters and on the reconstruction algorithm. We investigated the impact of the angular range of data acquisition and number of projection views on image quality by performing 3-D cascaded linear system modeling and phantom measurements. The results show that a larger angular range not only increases resolution perpendicular to the slices but also improves the visibility of low-frequency objects since the sampling area in Fourier space is mainly limited at low-frequencies. The effect on contrast resolution is demonstrated also with clinical data sets.

Impact of phosphor luminance noise on the specification of high-resolution CRT displays for medical imaging

We studied the impact of CRT spot size, phosphor luminance noise and image noise on the specification of high- resolution CRT displays that address the critical needs of general chest radiography. Using Argus CRT simulation software, the design of high-resolution CRTs for the display of adult chest radiographs was studied. The simulated images were printed on a laser printer and evaluated by a board- certified radiologist, RMS. The validity of the Argus simulation was assessed by modeling a 1k X 1k pixels CRT, whose technical parameters were sufficiently well known. Comments from the observer are presented comparing the simulated 2k display and a size-matched replicate of the original screen/film image. Critical parameters like phosphor luminance efficiency and its impact on electron beam size and phosphor luminance noise and its impact on radiographic image noise are discussed. We conclude that Argus CRT simulation software can successfully model the performance of CRTs intended to display medical images permitting consideration of critical parameters without costly manufacturing trials. Based on the 2k CRT simulation results, we suggest that a low luminance noise phosphor such as type p45 be used to ensure that specifying a small spot size would yield the anticipated sharpness improvements.