Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Biomedical diagnostics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

477 records · Page 27

Investigation of energy weighting using an energy discriminating photon counting detector for breast CT

Purpose: Breast CT is an emerging imaging technique that can portray the breast in 3D and improve visualization of important diagnostic features. Early clinical studies have suggested that breast CT has sufficient spatial and contrast resolution for accurate detection of masses and microcalcifications in the breast, reducing structural overlap that is often a limiting factor in reading mammographic images. For a number of reasons, image quality in breast CT may be improved by use of an energy resolving photon counting detector. In this study, the authors investigate the improvements in image quality obtained when using energy weighting with an energy resolving photon counting detector as compared to that with a conventional energy integrating detector.Methods: Using computer simulation, realistic CT images of multiple breast phantoms were generated. The simulation modeled a prototype breast CT system using an amorphous silicon (a-Si), CsI based energy integrating detector with different x-ray spectra, and a hypothetical, ideal CZT based photon counting detector with capability of energy discrimination. Three biological signals of interest were modeled as spherical lesions and inserted into breast phantoms; hydroxyapatite (HA) to represent microcalcification, infiltrating ductal carcinoma (IDC), and iodine enhanced infiltrating ductal carcinoma (IIDC). Signal-to-noise ratio (SNR) of these three lesions was measured from the CT reconstructions. In addition, a psychophysical study was conducted to evaluate observer performance in detecting microcalcifications embedded into a realistic anthropomorphic breast phantom.Results: In the energy range tested, improvements in SNR with a photon counting detector using energy weighting was higher (than the energy integrating detector method) by 30%–63% and 4%–34%, for HA and IDC lesions and 12%–30% (with Al filtration) and 32%–38% (with Ce filtration) for the IIDC lesion, respectively. The average area under the receiver operating characteristic curve (AUC) for detection of microcalcifications was higher by greater than 19% (for the different energy weighting methods tested) as compared to the AUC obtained with an energy integrating detector.Conclusions: This study showed that breast CT with a CZT photon counting detector using energy weighting can provide improvements in pixel SNR, and detectability of microcalcifications as compared to that with a conventional energy integrating detector. Since a number of degrading physical factors were not modeled into the photon counting detector, this improvement should be considered as an upper bound on achievable performance.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Energy dependent calibration of XR-QA2 radiochromic film with monochromatic and polychromatic x-ray beams

Purpose: This work investigates the energy response and dose-response curve determinations for XR-QA2 radiochromic film dosimetry system used for synchrotron radiation work and for quality assurance in diagnostic radiology, in the range of effective energies 18–46.5 keV. Methods: Pieces of XR-QA2 films were irradiated, in a plane transverse to the beam axis, with a monochromatic beam of energy in the range 18–40 keV at the ELETTRA synchrotron radiation facility (Trieste, Italy) and with a polychromatic beam from a laboratory x-ray tube operated at 80, 100, and 120 kV. The film calibration curve was expressed as air kerma (measured free-in-air with an ionization chamber) versus the net optical reflectance change (netΔR) derived from the red channel of the RGB scanned film image. Four functional relationships (rational, linear exponential, power, and logarithm) were tested to evaluate the best curve for fitting the calibration data. The adequacy of the various fitting functions was tested by using the uncertainty analysis and by assessing the average of the absolute air kerma error calculated as the difference between calculated and delivered air kerma. The sensitivity of the film was evaluated as the ratio of the change in net reflectance to the corresponding air kerma. Results: The sensitivity of XR-QA2 films increased in the energy range 18–39 keV, with a maximum variation of about 170%, and decreased in the energy range 38–46.5 keV. The present results confirmed and extended previous findings by this and other groups, as regards the dose response of the radiochromic film XR-QA2 to monochromatic and polychromatic x-ray beams, respectively. Conclusions: The XR-QA2 radiochromic film response showed a strong dependence on beam energy for both monochromatic and polychromatic beams in the range of half value layer values from 0.55 to 6.1 mm Al and corresponding effective energies from 18 to 46.5 keV. In this range, the film response varied by 170%, from a minimum sensitivity of 0.0127 to a maximum sensitivity of 0.0219 at 10 mGy air kerma in air. The more suitable function for air kerma calibration of the XR-QA2 radiochromic film was the power function. A significant batch-to-batch variation, up to 55%, in film response at 120 kV (46.5 keV effective energy) was observed in comparison with published data.

61 RADIATION PROTECTION AND DOSIMETRY↗

Denoising of polychromatic CT images based on their own noise properties

Purpose: Because of high diagnostic accuracy and fast scan time, computed tomography (CT) has been widely used in various clinical applications. Since the CT scan introduces radiation exposure to patients, however, dose reduction has recently been recognized as an important issue in CT imaging. However, low-dose CT causes an increase of noise in the image and thereby deteriorates the accuracy of diagnosis. In this paper, the authors develop an efficient denoising algorithm for low-dose CT images obtained using a polychromatic x-ray source. The algorithm is based on two steps: (i) estimation of space variant noise statistics, which are uniquely determined according to the system geometry and scanned object, and (ii) subsequent novel conversion of the estimated noise to Gaussian noise so that an existing high performance Gaussian noise filtering algorithm can be directly applied to CT images with non-Gaussian noise. Methods: For efficient polychromatic CT image denoising, the authors first reconstruct an image with the iterative maximum-likelihood polychromatic algorithm for CT to alleviate the beam-hardening problem. We then estimate the space-variant noise variance distribution on the image domain. Since there are many high performance denoising algorithms available for the Gaussian noise, image denoising can become much more efficient if they can be used. Hence, the authors propose a novel conversion scheme to transform the estimated space-variant noise to near Gaussian noise. In the suggested scheme, the authors first convert the image so that its mean and variance can have a linear relationship, and then produce a Gaussian image via variance stabilizing transform. The authors then apply a block matching 4D algorithm that is optimized for noise reduction of the Gaussian image, and reconvert the result to obtain a final denoised image. To examine the performance of the proposed method, an XCAT phantom simulation and a physical phantom experiment were conducted. Results: Both simulation and experimental results show that, unlike the existing denoising algorithms, the proposed algorithm can effectively reduce the noise over the whole region of CT images while preventing degradation of image resolution. Conclusions: To effectively denoise polychromatic low-dose CT images, a novel denoising algorithm is proposed. Because this algorithm is based on the noise statistics of a reconstructed polychromatic CT image, the spatially varying noise on the image is effectively reduced so that the denoised image will have homogeneous quality over the image domain. Through a simulation and a real experiment, it is verified that the proposed algorithm can deliver considerably better performance compared to the existing denoising algorithms.

60 APPLIED LIFE SCIENCES↗

TU-C-18C-01: Medical Physics 1.0 to 2.0: Introduction and Panel Discussion

Medical Physics 2.0, a new frontier in clinical imaging physics: Diagnostic imaging has always been a technological highlight of modern medicine. Imaging systems, with their ever-expanding advancement in terms of technology and application, increasingly require skilled expertise to understand the delicacy of their operation, monitor their performance, design their effective use, and ensure their overall quality and safety, scientifically and in quantitative terms. Physicists can play a crucial role in that process. But that role has largely remained a severely untapped resource. Many imaging centers fail to appreciate this potential, with medical physics groups either nonexistent or highly understaffed and their services poorly integrated into the patient care process. As a field, we have yet to define and enact how the clinical physicist can engage as an active, effective, and integral member of the clinical team, and how the services that she/he provides can be financially accounted for. Physicists do and will always contribute to research and development. However, their indispensible contribution to clinical imaging operations is something that has not been adequately established. That, in conjunction with new realities of healthcare practice, indicates a growing need to establish an updated approach to clinical medical imaging physics. This presentation aims to describe a vision as how clinical imaging physics can expand beyond traditional insular models of inspection and acceptance testing, oriented toward compliance, towards team-based models of operational engagement addressing topics such as new non-classical challenges of new technologies, quantitative imaging, and operational optimization. The Medical Physics 2.0 paradigm extends clinical medical physics from isolated characterization of inherent properties of the equipment to effective use of the equipment and to retrospective evaluation of clinical performance. This is an existential transition of the field that speaks to the new paradigms of value-based and evidence-based medicine, comparative effectiveness, and meaningful use. The panel discussion that follows includes prominent practitioners, thinkers, and leaders that would lead the discussion on how Medical Physics 2.0 can be actualized. Topics of discussion will include the administrative, financial, regulatory, and accreditation requirements of the new paradigm, effective models of practice, and the steps that we need to take to make MP 2.0 a reality. Learning Objectives: To understand the new paradigm of clinical medical physics practice extending from traditional insular models of compliance towards teambased models of operational engagement. To understand how clinical physics can most effectively contribute to clinical care. Learn to identify strengths and weaknesses in studies designed to measure the effect of low doses of ionizing radiation To recognize the impediments to Medical Physics 2.0 paradigm.

60 APPLIED LIFE SCIENCES↗

WE-A-BRF-01: Dual-Energy CT Imaging in Diagnostic Imaging and Radiation Therapy

The quantification accuracy of dual-energy imaging is influenced by the fundamentals of x-ray physics, system geometry, data acquisition hardware/protocol, system calibration, and image processing technique. This symposium will provide updates on the following advanced application areas: Mammography. Volumetric breast density techniques based on standard mammograms require estimation of breast thickness, which is difficult to accurately measure. By comparison, calculation of breast density using dual energy mammography does not require measurement of breast thickness. Dual energy mammography has been implemented using both energy integrating flat panel detectors in conjunction with beam energy switching and energy resolved photon counting detectors. These techniques have been optimized using simulation studies and validated using physical phantoms and postmortem breasts. Chemical decomposition was used as the gold standard for volumetric breast density measurement in postmortem breasts. Breast density measurements have also been compared with results from four-category BI-RADS density rankings, standard image thresholding and Fuzzy k-mean clustering techniques. These studies indicate that dual energy mammography can be used to accurately measure volumetric breast density. Cardiovascular CT. The predicative accuracy of risk models for recurrent stroke and cardiac arrest depends heavily on accurate differentiation of thrombus or calcium from iodine in left atrial appendage or coronary arteries. The amount of energy separation is constrained by image noise; therefore, optimal kVp, beam filtration, and balanced flux are essential for the quantification accuracy of iodine and calcium. The basis materials are combined linearly to generate monochromatic energy images, where CT# accuracy and CNR are energy dependent. With optimal monochromatic energy, the mean iodine concentration for the thrombus, circulatory stasis, and control groups are significantly different. Risk classification based on calcium scores shows excellent agreement with classification on the basis of conventional coronary artery calcium scoring. These studies demonstrate dual-energy cardiovascular CT can potentially be a noninvasive and sensitive modality in high risk patients. On-board KV/MV Imaging. To enhance soft tissue contrast and reduce metal artifacts, we have developed a dual-energy CBCT technique and a novel on-board kV/MV imaging technique based on hardware available on modern linear accelerators. We have also evaluated the feasibility of these two techniques in various phantom studies. Optimal techniques (energy, beam filtration, # of overlapping projections, etc) have been investigated with unique calibration procedures, which leads to successful decomposition of imaged material into acrylic-aluminum basis material pair. This enables the synthesis of virtual monochromatic (VM) CBCT images that demonstrate much less beam hardening, significantly reduced metal artifacts, and/or higher soft tissue CNR compared to single-energy CBCT. Adaptive Radiation Therapy. DECT could actually contribute to the area of Dose-Guided Radiation Therapy (or Adaptive Therapy). The application of DECT imaging using 80kV and 140 kV combinations could potentially increase the image quality by reducing the bone or high density material artifacts and also increase the soft tissue contrast by a light contrast agent. The result of this higher contrast / quality images is beneficial for deformable image registration / segmentation algorithm to improve its accuracy hence to make adaptive therapy less time consuming in its recontouring process. The real time re-planning prior to per treatment fraction could become more realistic with this improvement especially in hypofractional SBRT cases. Learning Objectives: Learn recent developments of dual-energy imaging in diagnosis and radiation therapy; Understand the unique clinical problem and required quantification accuracy in each application; Understand the different approaches to optimize dual-energy imaging techniques for different applications.

60 APPLIED LIFE SCIENCES↗

Method for measuring the focal spot size of an x-ray tube using a coded aperture mask and a digital detector

Purpose: The goal of this study is to evaluate a new method based on a coded aperture mask combined with a digital x-ray imaging detector for measurements of the focal spot sizes of diagnostic x-ray tubes. Common techniques for focal spot size measurements employ a pinhole camera, a slit camera, or a star resolution pattern. The coded aperture mask is a radiation collimator consisting of a large number of apertures disposed on a predetermined grid in an array, through which the radiation source is imaged onto a digital x-ray detector. The method of the coded mask camera allows one to obtain a one-shot accurate and direct measurement of the two dimensions of the focal spot (like that for a pinhole camera) but at a low tube loading (like that for a slit camera). A large number of small apertures in the coded mask operate as a ''multipinhole'' with greater efficiency than a single pinhole, but keeping the resolution of a single pinhole. Methods: X-ray images result from the multiplexed output on the detector image plane of such a multiple aperture array, and the image of the source is digitally reconstructed with a deconvolution algorithm. Images of the focal spot of a laboratory x-ray tube (W anode: 35-80 kVp; focal spot size of 0.04 mm) were acquired at different geometrical magnifications with two different types of digital detector (a photon counting hybrid silicon pixel detector with 0.055 mm pitch and a flat panel CMOS digital detector with 0.05 mm pitch) using a high resolution coded mask (type no-two-holes-touching modified uniformly redundant array) with 480 0.07 mm apertures, designed for imaging at energies below 35 keV. Measurements with a slit camera were performed for comparison. A test with a pinhole camera and with the coded mask on a computed radiography mammography unit with 0.3 mm focal spot was also carried out. Results: The full width at half maximum focal spot sizes were obtained from the line profiles of the decoded images, showing a focal spot of 0.120 mmx0.105 mm at 35 kVp and M=6.1, with a detector entrance exposure as low as 1.82 mR (0.125 mA s tube load). The slit camera indicated a focal spot of 0.112 mmx0.104 mm at 35 kVp and M=3.15, with an exposure at the detector of 72 mR. Focal spot measurements with the coded mask could be performed up to 80 kVp. Tolerance to angular misalignment with the reference beam up to 7 deg. in in-plane rotations and 1 deg. deg in out-of-plane rotations was observed. The axial distance of the focal spot from the coded mask could also be determined. It is possible to determine the beam intensity via measurement of the intensity of the decoded image of the focal spot and via a calibration procedure. Conclusions: Coded aperture masks coupled to a digital area detector produce precise determinations of the focal spot of an x-ray tube with reduced tube loading and measurement time, coupled to a large tolerance in the alignment of the mask.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Coronary artery analysis: Computer-assisted selection of best-quality segments in multiple-phase coronary CT angiography

Purpose: The authors are developing an automated method to identify the best-quality coronary arterial segment from multiple-phase coronary CT angiography (cCTA) acquisitions, which may be used by either interpreting physicians or computer-aided detection systems to optimally and efficiently utilize the diagnostic information available in multiple-phase cCTA for the detection of coronary artery disease. Methods: After initialization with a manually identified seed point, each coronary artery tree is automatically extracted from multiple cCTA phases using our multiscale coronary artery response enhancement and 3D rolling balloon region growing vessel segmentation and tracking method. The coronary artery trees from multiple phases are then aligned by a global registration using an affine transformation with quadratic terms and nonlinear simplex optimization, followed by a local registration using a cubic B-spline method with fast localized optimization. The corresponding coronary arteries among the available phases are identified using a recursive coronary segment matching method. Each of the identified vessel segments is transformed by the curved planar reformation (CPR) method. Four features are extracted from each corresponding segment as quality indicators in the original computed tomography volume and the straightened CPR volume, and each quality indicator is used as a voting classifier for the arterial segment. A weighted voting ensemble (WVE) classifier is designed to combine the votes of the four voting classifiers for each corresponding segment. The segment with the highest WVE vote is then selected as the best-quality segment. In this study, the training and test sets consisted of 6 and 20 cCTA cases, respectively, each with 6 phases, containing a total of 156 cCTA volumes and 312 coronary artery trees. An observer preference study was also conducted with one expert cardiothoracic radiologist and four nonradiologist readers to visually rank vessel segment quality. The performance of our automated method was evaluated by comparing the automatically identified best-quality segments identified by the computer to those selected by the observers. Results: For the 20 test cases, 254 groups of corresponding vessel segments were identified after multiple phase registration and recursive matching. The AI-BQ segments agreed with the radiologist’s top 2 ranked segments in 78.3% of the 254 groups (Cohen’s kappa 0.60), and with the 4 nonradiologist observers in 76.8%, 84.3%, 83.9%, and 85.8% of the 254 groups. In addition, 89.4% of the AI-BQ segments agreed with at least two observers’ top 2 rankings, and 96.5% agreed with at least one observer’s top 2 rankings. In comparison, agreement between the four observers’ top ranked segment and the radiologist’s top 2 ranked segments were 79.9%, 80.7%, 82.3%, and 76.8%, respectively, with kappa values ranging from 0.56 to 0.68. Conclusions: The performance of our automated method for selecting the best-quality coronary segments from a multiple-phase cCTA acquisition was comparable to the selection made by human observers. This study demonstrates the potential usefulness of the automated method in clinical practice, enabling interpreting physicians to fully utilize the best available information in cCTA for diagnosis of coronary disease, without requiring manual search through the multiple phases and minimizing the variability in image phase selection for evaluation of coronary artery segments across the diversity of human readers with variations in expertise.

61 RADIATION PROTECTION AND DOSIMETRY↗

MABAL: a Novel Deep-Learning Architecture for Machine-Assisted Bone Age Labeling

Bone age assessment (BAA) is a commonly performed diagnostic study in pediatric radiology to assess skeletal maturity. The most commonly utilized method for assessment of BAA is the Greulich and Pyle method (Pediatr Radiol 46.9:1269–1274, 2016; Arch Dis Child 81.2:172–173, 1999) atlas. The evaluation of BAA can be a tedious and time-consuming process for the radiologist. As such, several computer-assisted detection/diagnosis (CAD) methods have been proposed for automation of BAA. Classical CAD tools have traditionally relied on hard-coded algorithmic features for BAA which suffer from a variety of drawbacks. Recently, the advent and proliferation of convolutional neural networks (CNNs) has shown promise in a variety of medical imaging applications. There have been at least two published applications of using deep learning for evaluation of bone age (Med Image Anal 36:41–51, 2017; JDI 1–5, 2017). However, current implementations are limited by a combination of both architecture design and relatively small datasets. The purpose of this study is to demonstrate the benefits of a customized neural network algorithm carefully calibrated to the evaluation of bone age utilizing a relatively large institutional dataset. In doing so, this study will aim to show that advanced architectures can be successfully trained from scratch in the medical imaging domain and can generate results that outperform any existing proposed algorithm. The training data consisted of 10,289 images of different skeletal age examinations, 8909 from the hospital Picture Archiving and Communication System at our institution and 1383 from the public Digital Hand Atlas Database. The data was separated into four cohorts, one each for male and female children above the age of 8, and one each for male and female children below the age of 10. The testing set consisted of 20 radiographs of each 1-year-age cohort from 0 to 1 years to 14–15+ years, half male and half female. The testing set included left-hand radiographs done for bone age assessment, trauma evaluation without significant findings, and skeletal surveys. A 14 hidden layer-customized neural network was designed for this study. The network included several state of the art techniques including residual-style connections, inception layers, and spatial transformer layers. Data augmentation was applied to the network inputs to prevent overfitting. A linear regression output was utilized. Mean square error was used as the network loss function and mean absolute error (MAE) was utilized as the primary performance metric. MAE accuracies on the validation and test sets for young females were 0.654 and 0.561 respectively. For older females, validation and test accuracies were 0.662 and 0.497 respectively. For young males, validation and test accuracies were 0.649 and 0.585 respectively. Finally, for older males, validation and test set accuracies were 0.581 and 0.501 respectively. The female cohorts were trained for 900 epochs each and the male cohorts were trained for 600 epochs. An eightfold cross-validation set was employed for hyperparameter tuning. Test error was obtained after training on a full data set with the selected hyperparameters. Using our proposed customized neural network architecture on our large available data, we achieved an aggregate validation and test set mean absolute errors of 0.637 and 0.536 respectively. To date, this is the best published performance on utilizing deep learning for bone age assessment. Our results support our initial hypothesis that customized, purpose-built neural networks provide improved performance over networks derived from pre-trained imaging data sets. We build on that initial work by showing that the addition of state-of-the-art techniques such as residual connections and inception architecture further improves prediction accuracy. This is important because the current assumption for use of residual and/or inception architectures is that a large pre-trained network is required for successful implementation given the relatively small datasets in medical imaging. Instead we show that a small, customized architecture incorporating advanced CNN strategies can indeed be trained from scratch, yielding significant improvements in algorithm accuracy. It should be noted that for all four cohorts, testing error outperformed validation error. One reason for this is that our ground truth for our test set was obtained by averaging two pediatric radiologist reads compared to our training data for which only a single read was used. This suggests that despite relatively noisy training data, the algorithm could successfully model the variation between observers and generate estimates that are close to the expected ground truth.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

TH-E-18A-01: Developments in Monte Carlo Methods for Medical Imaging

Monte Carlo simulation methods are widely used in medical physics research and are starting to be implemented in clinical applications such as radiation therapy planning systems. Monte Carlo simulations offer the capability to accurately estimate quantities of interest that are challenging to measure experimentally while taking into account the realistic anatomy of an individual patient. Traditionally, practical application of Monte Carlo simulation codes in diagnostic imaging was limited by the need for large computational resources or long execution times. However, recent advancements in high-performance computing hardware, combined with a new generation of Monte Carlo simulation algorithms and novel postprocessing methods, are allowing for the computation of relevant imaging parameters of interest such as patient organ doses and scatter-to-primaryratios in radiographic projections in just a few seconds using affordable computational resources. Programmable Graphics Processing Units (GPUs), for example, provide a convenient, affordable platform for parallelized Monte Carlo executions that yield simulation times on the order of 10{sup 7} xray/ s. Even with GPU acceleration, however, Monte Carlo simulation times can be prohibitive for routine clinical practice. To reduce simulation times further, variance reduction techniques can be used to alter the probabilistic models underlying the x-ray tracking process, resulting in lower variance in the results without biasing the estimates. Other complementary strategies for further reductions in computation time are denoising of the Monte Carlo estimates and estimating (scoring) the quantity of interest at a sparse set of sampling locations (e.g. at a small number of detector pixels in a scatter simulation) followed by interpolation. Beyond reduction of the computational resources required for performing Monte Carlo simulations in medical imaging, the use of accurate representations of patient anatomy is crucial to the virtual generation of medical images and accurate estimation of radiation dose and other imaging parameters. For this, detailed computational phantoms of the patient anatomy must be utilized and implemented within the radiation transport code. Computational phantoms presently come in one of three format types, and in one of four morphometric categories. Format types include stylized (mathematical equation-based), voxel (segmented CT/MR images), and hybrid (NURBS and polygon mesh surfaces). Morphometric categories include reference (small library of phantoms by age at 50th height/weight percentile), patient-dependent (larger library of phantoms at various combinations of height/weight percentiles), patient-sculpted (phantoms altered to match the patient's unique outer body contour), and finally, patient-specific (an exact representation of the patient with respect to both body contour and internal anatomy). The existence and availability of these phantoms represents a very important advance for the simulation of realistic medical imaging applications using Monte Carlo methods. New Monte Carlo simulation codes need to be thoroughly validated before they can be used to perform novel research. Ideally, the validation process would involve comparison of results with those of an experimental measurement, but accurate replication of experimental conditions can be very challenging. It is very common to validate new Monte Carlo simulations by replicating previously published simulation results of similar experiments. This process, however, is commonly problematic due to the lack of sufficient information in the published reports of previous work so as to be able to replicate the simulation in detail. To aid in this process, the AAPM Task Group 195 prepared a report in which six different imaging research experiments commonly performed using Monte Carlo simulations are described and their results provided. The simulation conditions of all six cases are provided in full detail, with all necessary data on material composition, source, geometry, scoring and other parameters provided. The results of these simulations when performed with the four most common publicly available Monte Carlo packages are also provided in tabular form. The Task Group 195 Report will be useful for researchers needing to validate their Monte Carlo work, and for trainees needing to learn Monte Carlo simulation methods. In this symposium we will review the recent advancements in highperformance computing hardware enabling the reduction in computational resources needed for Monte Carlo simulations in medical imaging. We will review variance reduction techniques commonly applied in Monte Carlo simulations of medical imaging systems and present implementation strategies for efficient combination of these techniques with GPU acceleration. Trade-offs involved in Monte Carlo acceleration by means of denoising and “sparse sampling” will be discussed. A method for rapid scatter correction in cone-beam CT (<5 min/scan) will be presented as an illustration of the simulation speeds achievable with optimized Monte Carlo simulations. We will also discuss the development, availability, and capability of the various combinations of computational phantoms for Monte Carlo simulation of medical imaging systems. Finally, we will review some examples of experimental validation of Monte Carlo simulations and will present the AAPM Task Group 195 Report. Learning Objectives: Describe the advances in hardware available for performing Monte Carlo simulations in high performance computing environments. Explain variance reduction, denoising and sparse sampling techniques available for reduction of computational time needed for Monte Carlo simulations of medical imaging. List and compare the computational anthropomorphic phantoms currently available for more accurate assessment of medical imaging parameters in Monte Carlo simulations. Describe experimental methods used for validation of Monte Carlo simulations in medical imaging. Describe the AAPM Task Group 195 Report and its use for validation and teaching of Monte Carlo simulations in medical imaging.

60 APPLIED LIFE SCIENCES↗