Aldo Badano

Accelerating Monte Carlo simulations of photon transport in a voxelized geometry using a massively parallel graphics processing unit

PURPOSE It is a known fact that Monte Carlo simulations of radiation transport are computationally intensive and may require long computing times. The authors introduce a new paradigm for the acceleration of Monte Carlo simulations: The use of a graphics processing unit (GPU) as the main computing device instead of a central processing unit (CPU). METHODS A GPU-based Monte Carlo code that simulates photon transport in a voxelized geometry with the accurate physics models from PENELOPE has been developed using the CUDATM programming model (NVIDIA Corporation, Santa Clara, CA). RESULTS An outline of the new code and a sample x-ray imaging simulation with an anthropomorphic phantom are presented. A remarkable 27-fold speed up factor was obtained using a GPU compared to a single core CPU. CONCLUSIONS The reported results show that GPUs are currently a good alternative to CPUs for the simulation of radiation transport. Since the performance of GPUs is currently increasing at a faster pace than that of CPUs, the advantages of GPU-based software are likely to be more pronounced in the future.

MANTIS: combined x-ray, electron and optical Monte Carlo simulations of indirect radiation imaging systems

We describe MANTIS (Monte carlo x-rAy electroN opTical Imaging Simulation), a tool for simulating imaging systems that tracks x-rays, electrons and optical photons in arbitrary materials and complex geometries. The x-ray and electron transport and involved physics models are from the PENELOPE package, and the optical transport and corresponding physics models are from DETECT-II and include Fresnel refraction and reflection at material boundaries, bulk absorption and scattering. Complex geometries can be handled with the aid of the geometry routines included in PENELOPE. When x-rays or electrons interact and deposit energy in the scintillator, the code generates a number of optical quanta according to a user-selected model for the conversion process. The optical photons are then tracked until they reach an absorption event, which in some cases contributes to the output signal, or escape from the geometry. We demonstrate the capabilities of this new tool with respect to the statistics of the optical signal detected and to the three-dimensional point-response functions corresponding to columnar phosphor screens.

Channelized Hotelling observers for the assessment of volumetric imaging data sets

Current clinical practice is rapidly moving in the direction of volumetric imaging. For two-dimensional (2D) images, task-based medical image quality is often assessed using numerical model observers. For three-dimensional (3D) images, however, these models have been little explored so far. In this work, first, two novel designs of a multislice channelized Hotelling observer (CHO) are proposed for the task of detecting 3D signals in 3D images. The novel designs are then compared and evaluated in a simulation study with five different CHO designs: a single-slice model, three multislice models, and a volumetric model. Four different random background statistics are considered, both gaussian (noncorrelated and correlated gaussian noise) and non-gaussian (lumpy and clustered lumpy backgrounds). Overall, the results show that the volumetric model outperforms the others, while the disparity between the models decreases for greater complexity of the detection task. Among the multislice models, the second proposed CHO could most closely approach the volumetric model, whereas the first new CHO seems to be least affected by the number of training samples.

Observer Variability in the Interpretation of HER2/neu Immunohistochemical Expression With Unaided and Computer-Aided Digital Microscopy

CONTEXT Observer variability in digital microscopy and the effect of computer-aided digital microscopy are underexamined areas in need of further research, considering the increasing use and future role of digital imaging in pathology. A reduction in observer variability using computer aids could enhance the statistical power of studies designed to determine the utility of new biomarkers and accelerate their incorporation in clinical practice. OBJECTIVES To quantify interobserver and intraobserver variability in immunohistochemical analysis of HER2/neu with digital microscopy and computer-aided digital microscopy, and to test the hypothesis that observer agreement in the quantitative assessment of HER2/neu immunohistochemical expression is increased with the use of computer-aided microscopy. DESIGN A set of 335 digital microscopy images extracted from 64 breast cancer tissue slides stained with a HER2 antibody, were read by 14 observers in 2 reading modes: the unaided mode and the computer-aided mode. In the unaided mode, HER2 images were displayed on a calibrated color monitor with no other information, whereas in the computer-aided mode, observers were shown a HER2 image along with a corresponding feature plot showing computer-extracted values of membrane staining intensity and membrane completeness for the particular image under examination and, at the same time, mean feature values of the different HER2 categories. In both modes, observers were asked to provide a continuous score of HER2 expression. RESULTS Agreement analysis performed on the output of the study showed significant improvement in both interobserver and intraobserver agreement when the computer-aided reading mode was used to evaluate preselected image fields. CONCLUSION The role of computer-aided digital microscopy in reducing observer variability in immunohistochemistry is promising.

A statistical, task-based evaluation method for three-dimensional x-ray breast imaging systems using variable-background phantoms

PURPOSE For the last few years, development and optimization of three-dimensional (3D) x-ray breast imaging systems, such as digital breast tomosynthesis (DBT) and computed tomography, have drawn much attention from the medical imaging community, either academia or industry. However, there is still much room for understanding how to best optimize and evaluate the devices over a large space of many different system parameters and geometries. Current evaluation methods, which work well for 2D systems, do not incorporate the depth information from the 3D imaging systems. Therefore, it is critical to develop a statistically sound evaluation method to investigate the usefulness of inclusion of depth and background-variability information into the assessment and optimization of the 3D systems. METHODS In this paper, we present a mathematical framework for a statistical assessment of planar and 3D x-ray breast imaging systems. Our method is based on statistical decision theory, in particular, making use of the ideal linear observer called the Hotelling observer. We also present a physical phantom that consists of spheres of different sizes and materials for producing an ensemble of randomly varying backgrounds to be imaged for a given patient class. Lastly, we demonstrate our evaluation method in comparing laboratory mammography and three-angle DBT systems for signal detection tasks using the phantoms projection data. We compare the variable phantom case to that of a phantom of the same dimensions filled with water, which we call the uniform phantom, based on the performance of the Hotelling observer as a function of signal size and intensity. RESULTS Detectability trends calculated using the variable and uniform phantom methods are different from each other for both mammography and DBT systems. CONCLUSIONS Our results indicate that measuring the systems detection performance with consideration of background variability may lead to differences in system performance estimates and comparisons. For the assessment of 3D systems, to accurately determine trade offs between image quality and radiation dose, it is critical to incorporate randomness arising from the imaging chain including background variability into system performance calculations.

Angular dependence of the luminance and contrast in medical monochrome liquid crystal displays

Active-matrix liquid crystal displays (AMLCDs) are light-modulating devices that generate images by differentially transmitting a nearly uniform luminous field provided by a backlight. While emissive displays exhibit a quasi-Lambertian emission with almost constant contrast at off-normal viewing, the anisotropy of the electro-optic effect that controls light transmission in AMLCDs causes a pixel luminance that varies, sometimes strongly, with viewing angle. These variations are not identical for all gray levels and can eventually cause grayscale inversions. In this paper, we measured the luminance emission of a monochrome medical AMLCD, a medical cathode-ray tube monitor, and a color desktop AMLCD, using a collimated photopic probe positioned on a manual rotation arm, and a research radiometer with automatic readout. The probe measures luminance with a small acceptance angle and provides optical shielding from emissions at other viewing directions that contaminate the readings. We obtained luminance response curves versus angle in the vertical, horizontal and at 45 degrees diagonal directions. The display systems were calibrated to reflect the DICOM Part 3.14 standard grayscale display function (GDF) when measured using the manufacturers probe and software tools. We analyzed the measurements at different viewing directions with respect to their departure from the GDF by computing the normalized contrast (deltaL/L) as a function of the DICOM just-noticeable difference index. Although cathode-ray tubes are known to be quasi-Lambertian emitters, the luminance at normal viewing is higher than the luminance observed at large angles. This decrease in luminance is however proportionally similar for all gray levels, resulting in a relatively flat contrast response for all angles. In addition to being more pronounced, the angular variation in AMLCDs does not follow the same profile at different intensities with the subsequent variation in the achieved display contrast. The changes due to off-normal viewing are substantial at large angles in the horizontal and vertical directions, and much worse in the diagonal viewing directions.

Anisotropic imaging performance in indirect x‐ray imaging detectors

We report on the variability in imaging system performance due to oblique x-ray incidence, and the associated transport of quanta (both x rays and optical photons) through the phosphor, in columnar indirect digital detectors. The analysis uses MANTIS, a combined x-ray, electron, and optical Monte Carlo transport code freely available. We describe the main features of the simulation method and provide some validation of the phosphor screen models considered in this work. We report x-ray and electron three-dimensional energy deposition distributions and point-response functions (PRFs), including optical spread in columnar phosphor screens of thickness 100 and 500 microm, for 19, 39, 59, and 79 keV monoenergetic x-ray beams incident at 0 degrees, 10 degrees, and 15 degrees. In addition, we present pulse-height spectra for the same phosphor thickness, x-ray energies, and angles of incidence. Our results suggest that the PRF due to the phosphor blur is highly nonsymmetrical, and that the resolution properties of a columnar screen in a tomographic, or tomosynthetic imaging system varies significantly with the angle of x-ray incidence. Moreover, we find that the noise due to the variability in the number of light photons detected per primary x-ray interaction, summarized in the information or Swank factor, is somewhat independent of thickness and incidence angle of the x-ray beam. Our results also suggest that the anisotropy in the PRF is not less in screens with absorptive backings, while the noise introduced by variations in the gain and optical transport is larger. Predictions from MANTIS, after additional validation, can provide the needed understanding of the extent of such variations, and eventually, lead to the incorporation of the changes in imaging performance with incidence angle into the reconstruction algorithms for volumetric x-ray imaging systems.

Image Quality Degradation by Light Scattering in Display Devices

Veiling glare and ambient light reflection can significantly degrade the quality of an image on a display device. Veiling glare is primarily associated with the diffuse spread of image signal caused by multiple light scattering in the emissive structure of the device. The glare ratio associated with a test image with a 1-cm-diameter black spot is reported as 555 for film, 89 for a monochrome monitor, and 25 for a color monitor. Diffuse light reflection results from ambient light entering the display surface and returning at random emission angles. The diffuse reflection coefficient (luminance/illuminance, 1/sr) is reported as 0.026 for film, 0.058 for a monochrome monitor, and 0.025 for a color monitor with an antireflective surface coating. Both processes increase the luminance in black regions and cause contrast reduction. Specular reflections interfere with detail in the displayed scene. The specular reflection coefficient (luminance/luminance) is reported as 0.011 for film, 0.041 for a monochrome monitor, and 0.021 for a color monitor with an antireflective coating.

penMesh —Monte Carlo Radiation Transport Simulation in a Triangle Mesh Geometry

We have developed a general-purpose Monte Carlo simulation code, called penMesh, that combines the accuracy of the radiation transport physics subroutines from PENELOPE and the flexibility of a geometry based on triangle meshes. While the geometric models implemented in most general-purpose codes-such as PENELOPEs quadric geometry-impose some limitations in the shape of the objects that can be simulated, triangle meshes can be used to describe any free-form (arbitrary) object. Triangle meshes are extensively used in computer-aided design and computer graphics. We took advantage of the sophisticated tools already developed in these fields, such as an octree structure and an efficient ray-triangle intersection algorithm, to significantly accelerate the triangle mesh ray-tracing. A detailed description of the new simulation code and its ray-tracing algorithm is provided in this paper. Furthermore, we show how it can be readily used in medical imaging applications thanks to the detailed anatomical phantoms already available. In particular, we present a whole body radiography simulation using a triangulated version of the anthropomorphic NCAT phantom. An example simulation of scatter fraction measurements using a standardized abdomen and lumbar spine phantom, and a benchmark of the triangle mesh and quadric geometries in the ray-tracing of a mathematical breast model, are also presented to show some of the capabilities of penMesh.

Anisotropic imaging performance in breast tomosynthesis.

We describe the anisotropy in imaging performance caused by oblique x-ray incidence in indirect detectors for breast tomosynthesis based on columnar scintillator screens. We use MANTIS, a freely available combined x-ray, electron, and optical Monte Carlo transport package which models the indirect detection processes in columnar screens, interaction by interaction. The code has been previously validated against published optical distributions. In this article, initial validation results are provided concerning the blur for particular designs of phosphor screens for which some details with respect to the columnar geometry are available from scanning electron microscopy. The polyenergetic x-ray spectrum utilized comes from a database of experimental data for three different anode/filter/kVp combinations: Mo/Mo at 28 kVp, Rh/Rh at 28 kVp, and W/Al at 42 kVp. The x-ray spectra were then filtered with breast tissue (3, 4, and 6 cm thickness), compression paddle, and support base, according to the oblique paths determined by the incidence angle. The composition of the breast tissue was 50%/50% adipose/glandular tissue mass ratio. Results are reported on the pulse-height statistics of the light output and on spatial blur, expressed as the response of the detector to a pencil beam with a certain incidence angle. Results suggest that the response is nonsymmetrical and that the resolution properties of a tomosynthesis system vary significantly with the angle of x-ray incidence. In contrast, it is found that the noise due to the variability in the number of light photons detected per primary x-ray interaction changes only a few percent. The anisotropy in the response is not less in screens with absorptive backings while the noise introduced by variations in the depth-dependent light output and optical transport is larger. The results suggest that anisotropic imaging performance across the detector area can be incorporated into reconstruction algorithms for improving the image quality of breast tomosynthesis. This study also demonstrates that the assessment of image quality of breast tomosynthesis systems requires a more complete description of the detector response beyond local, center measurements of resolution and noise that assume some degree of symmetry in the detector performance.