Modeling the performance characteristics of computed radiography (CR) systems

Maximizing microcalcification detectability in low-dose dedicated cone-beam breast CT: parallel cascades-based theoretical analysis

Purpose

We aim to determine the combination of X-ray spectrum and detector scintillator thickness that maximizes the detectability of microcalcification clusters in dedicated cone-beam breast CT.

Approach

A cascaded linear system analysis was implemented in the spatial frequency domain and was used to determine the detectability index using numerical observers for the imaging task of detecting a microcalcification cluster with 0.17 mm diameter calcium carbonate spheres. The analysis considered a thallium-doped cesium iodide scintillator coupled to a complementary metal-oxide semiconductor detector and an analytical filtered-back-projection reconstruction algorithm. Independent system parameters considered were the scintillator thickness, applied X-ray tube voltage, and X-ray beam filtration. The combination of these parameters that maximized the detectability index was considered optimal.

Results

Prewhitening, nonprewhitening, and nonprewhitening with eye filter numerical observers indicate that the combination of 0.525 to 0.6 mm thick scintillator, 70 kV, and 0.25 to 0.4 mm added copper filtration maximized the detectability index at a mean glandular dose (MGD) of 4.5 mGy.

Conclusion

Using parallel cascade systems' analysis, the combination of parameters that could maximize the detection of microcalcifications was identified. The analysis indicates that a harder beam than that used in current practice may be beneficial for the task of detecting microcalcifications at an MGD suitable for breast cancer screening.

Development of a computational tool for estimating computed tomography dose parameters

BACKGROUND

Computed Tomographic (CT) imaging procedures have been reported as the main source of radiation in diagnostic procedures compared to other modalities. To provide the optimal quality of CT images at the minimum radiation risk to the patient, periodic inspections and calibration tests for CT equipment are required. These tests involve a series of measurements that are time consuming and may require specific skills and highly-trained personnel.

OBJECTIVE

This study aims to develop a new computational tool to estimate the dose of CT radiation outputs and assist in the calibration of CT scanners. It may also provide an educational resource by which radiological practitioners can learn the influence of technique factors on both patient radiation dose and the produced image quality.

METHODS

The computational tool was developed using MATLAB in order to estimate the CT radiation dose parameters for different technique factors. The CT radiation dose parameters were estimated from the calibrated energy spectrum of the x-ray tube for a CT scanner.

RESULTS

The estimated dose parameters and the measured values utilising an Adult CT Head Dose Phantom showed linear correlations for different tube voltages (80 kVp, 100 kVp, 120 kVp, and 140 kVp), with R2 nearly equal to 1 (0.99). The maximum differences between the estimated and measured CTDIvol were under 5 %. For 80 kVp and low tube currents (50 mA, 100 mA), the maximum differences were under 10%.

CONCLUSIONS

The prototyped computational model provides a tool for the simulation of a machine-specific spectrum and CT dose parameters using a single dose measurement.

Readout models for BaFBr0.85I0.15:Eu image plates

The linearity of the photostimulated luminescence process makes repeated image-plate scanning a viable technique to extract a more dynamic range. In order to obtain a response estimate, two semi-empirical models for the readout fading of an image plate are introduced; they relate the depth distribution of activated photostimulated luminescence centers within an image plate to the recorded signal. Model parameters are estimated from image-plate scan series with BAS-MS image plates and the Typhoon FLA 7000 scanner for the hard x-ray image-plate diagnostic over a collection of experiments providing x-ray energy spectra whose approximate shape is a double exponential.

Emerging Breast Imaging Technologies on the Horizon

Early detection of breast cancers by mammography in conjunction with adjuvant therapy has contributed to reduction in breast cancer mortality. Mammography remains the "gold-standard" for breast cancer screening but is limited by tissue superposition. Digital breast tomosynthesis and more recently, dedicated breast computed tomography have been developed to alleviate the tissue superposition problem. However, all of these modalities rely upon x-ray attenuation contrast to provide anatomical images, and there are ongoing efforts to develop and clinically translate alternative modalities. These emerging modalities could provide for new contrast mechanisms and may potentially improve lesion detection and diagnosis. In this article, several of these emerging modalities are discussed with a focus on technologies that have advanced to the stage of in vivo clinical evaluation.

A comprehensive model for quantum noise characterization in digital mammography

A version of cascaded systems analysis was developed specifically with the aim of studying quantum noise propagation in x-ray detectors. Signal and quantum noise propagation was then modelled in four types of x-ray detectors used for digital mammography: four flat panel systems, one computed radiography and one slot-scan silicon wafer based photon counting device. As required inputs to the model, the two dimensional (2D) modulation transfer function (MTF), noise power spectra (NPS) and detective quantum efficiency (DQE) were measured for six mammography systems that utilized these different detectors. A new method to reconstruct anisotropic 2D presampling MTF matrices from 1D radial MTFs measured along different angular directions across the detector is described; an image of a sharp, circular disc was used for this purpose. The effective pixel fill factor for the FP systems was determined from the axial 1D presampling MTFs measured with a square sharp edge along the two orthogonal directions of the pixel lattice. Expectation MTFs were then calculated by averaging the radial MTFs over all possible phases and the 2D EMTF formed with the same reconstruction technique used for the 2D presampling MTF. The quantum NPS was then established by noise decomposition from homogenous images acquired as a function of detector air kerma. This was further decomposed into the correlated and uncorrelated quantum components by fitting the radially averaged quantum NPS with the radially averaged EMTF(2). This whole procedure allowed a detailed analysis of the influence of aliasing, signal and noise decorrelation, x-ray capture efficiency and global secondary gain on NPS and detector DQE. The influence of noise statistics, pixel fill factor and additional electronic and fixed pattern noises on the DQE was also studied. The 2D cascaded model and decompositions performed on the acquired images also enlightened the observed quantum NPS and DQE anisotropy.

Pixel Bleeding Correction in Laser Scanning Luminescence Imaging Demonstrated Using Optically Stimulated Luminescence

This paper describes and investigates the performance of an algorithm to correct for "pixel bleeding" caused by slow luminescence centers in laser scanning imaging (e.g., X-ray imaging using photostimulable phosphors and 2D dosimetry using optically stimulated luminescence). The algorithm is based on a deconvolution procedure that takes into account the lifetime of the slow luminescence center and is further constrained by the detection of fast and slow luminescence centers and combining rows scanned in opposite directions. The algorithm was tested using simulated data and demonstrated experimentally by applying it to image reconstruction of two types of Al2O3 X-ray detector films ( Al2O3:C and Al2O3 :C,Mg), whose use in 2D dosimetry in conjunction with laser-scanning readout has so far been prevented by slow luminescence centers (F-centers, 35 ms lifetime). We show that the algorithm allows the readout of Al2O3 film detectors 300-500 times faster than generally allowed considering the lifetime of the main luminescence centers. By relaxing the stringent requirements on the detector's luminescence lifetime, the algorithm opens the possibility of using new materials in 2D dosimetry as well as other laser scanning applications, such as X-ray imaging using storage phosphors and scanning confocal microscopy, although the effect of the noise introduced must be investigated for each specific application.

Optical absorption characteristics in the assessment of powder phosphor-based x-ray detectors: from nano- to micro-scale

X-ray phosphor-based detectors have enormously improved the quality of medical imaging examinations through the optimization of optical diffusion. In recent years, with the development of science and technology in the field of materials, improved powder phosphors require structural and optical properties that contribute to better optical signal propagation. The purpose of this paper was to provide a quantitative and qualitative understanding of the optical absorption characteristics in the assessment of powder phosphor-based detectors (from nano- scale up to micro-scale). Variations on the optical absorption parameters (i.e. the light extinction coefficient [Formula: see text] and the percentage probability of light absorption p%) were evaluated based on Mie calculations examining a wide range of light wavelengths, particle refractive indices and sizes. To model and assess the effects of the aforementioned parameters on optical diffusion, Monte Carlo simulation techniques were employed considering: (i) phosphors of different layer thickness, 100 μm (thin layer) and 300 μm (thick layer), respectively, (ii) light extinction coefficient values, 1, 3 and 6 μm(-1), and (iii) percentage probability of light absorption p% in the range 10(-4)-10(-2). Results showed that the [Formula: see text] coefficient is high for phosphor grains in the submicron scale and for low light wavelengths. At higher wavelengths (above 650 nm), optical quanta follow approximately similar depths until interaction for grain diameter 500 nm and 1 μm. Regarding the variability of the refractive index, high variations of the [Formula: see text] coefficient occurred above 1.6. Furthermore, results derived from Monte Carlo modeling showed that high spatial resolution phosphors can be accomplished by increasing the [Formula: see text] parameter. More specifically, the FWHM was found to decrease (i.e. higher resolution): (i) 4.8% at 100 μm and (ii) 9.5%, at 300 μm layer thickness. This study attempted to examine the role of the optical absorption parameters on optical diffusion studies. A significant outcome of the present investigation was that the improvement of phosphor spatial resolution without decreasing the light collection efficiency too much can be better achieved by increasing the parameter [Formula: see text] rather than the parameter p%.

Image simulation and a model of noise power spectra across a range of mammographic beam qualities

PURPOSE

The aim of this work is to create a model to predict the noise power spectra (NPS) for a range of mammographic radiographic factors. The noise model was necessary to degrade images acquired on one system to match the image quality of different systems for a range of beam qualities.

METHODS

Five detectors and x-ray systems [Hologic Selenia (ASEh), Carestream computed radiography CR900 (CRc), GE Essential (CSI), Carestream NIP (NIPc), and Siemens Inspiration (ASEs)] were characterized for this study. The signal transfer property was measured as the pixel value against absorbed energy per unit area (E) at a reference beam quality of 28 kV, Mo/Mo or 29 kV, W/Rh with 45 mm polymethyl methacrylate (PMMA) at the tube head. The contributions of the three noise sources (electronic, quantum, and structure) to the NPS were calculated by fitting a quadratic at each spatial frequency of the NPS against E. A quantum noise correction factor which was dependent on beam quality was quantified using a set of images acquired over a range of radiographic factors with different thicknesses of PMMA. The noise model was tested for images acquired at 26 kV, Mo/Mo with 20 mm PMMA and 34 kV, Mo/Rh with 70 mm PMMA for three detectors (ASEh, CRc, and CSI) over a range of exposures. The NPS were modeled with and without the noise correction factor and compared with the measured NPS. A previous method for adapting an image to appear as if acquired on a different system was modified to allow the reference beam quality to be different from the beam quality of the image. The method was validated by adapting the ASEh flat field images with two thicknesses of PMMA (20 and 70 mm) to appear with the imaging characteristics of the CSI and CRc systems.

RESULTS

The quantum noise correction factor rises with higher beam qualities, except for CR systems at high spatial frequencies, where a flat response was found against mean photon energy. This is due to the dominance of secondary quantum noise in CR. The use of the quantum noise correction factor reduced the difference from the model to the real NPS to generally within 4%. The use of the quantum noise correction improved the conversion of ASEh image to CRc image but had no difference for the conversion to CSI images.

CONCLUSIONS

A practical method for estimating the NPS at any dose and over a range of beam qualities for mammography has been demonstrated. The noise model was incorporated into a methodology for converting an image to appear as if acquired on a different detector. The method can now be extended to work for a wide range of beam qualities and can be applied to the conversion of mammograms.

Light wavelength effects in submicrometer phosphor materials using Mie scattering and Monte Carlo simulation

PURPOSE

Phosphor materials provide challenges to both fundamental research and breakthrough development of technologies in research areas. In recent years, with the development of science and technology in the field of materials, a number of physical or chemical synthesis methods have been developed and successfully used for the preparation of submicrometer-sized phosphors. The present paper provides a rigorous analysis of light diffusion capabilities of phosphor materials in submicrometer-scale investigating the effect of light wavelength.

METHODS

Mie scattering theory and Monte Carlo simulation techniques were used for the optical diffusion performance providing numerical calculations. The Monte Carlo model included: (i) phosphor layers composed of different thickness (200, 500, 1000 μm) and (ii) different light wavelength values (420, 545, 610 nm) corresponding to different types of activators, such as Ce, Tb, and Eu activators, respectively.

RESULTS

Based on Mie calculations, it was found that for low values of refractive index (e.g., 1.4) and for particle radius from 250 up to 500 nm no significant variations occurred on optical parameters. Monte Carlo simulations showed that the resolution increases as light wavelength decreases, respectively, however, this increase is more obvious at lower thickness values (i.e., at 200 μm). In particular, as light wavelength decreases from 610 down to 545 and 420 nm, the resolution increases 4.4% and 13.9%, respectively (at 200 μm layer thickness). In addition, as layer thickness increases from 200 up to 500 μm the resolution decreases 50.2% while an increase up to 1000 μm causes a decrease of 70.2% (at 420 nm light wavelength).

CONCLUSIONS

The goal of the author's study was to investigate the optical diffusion characteristics of submicrometer phosphor materials using Mie scattering theory and Monte Carlo simulation. The present investigation indicated that a key parameter on resolution improvement was the amount of light loss which depends on the choice of activator and affects the lateral spreading.