The 100 most-cited articles in Medical Physics, 2016

The increasing complexity of modern radiation therapy planning and delivery challenges traditional prescriptive quality management (QM) methods, such as many of those included in guidelines published by organizations such as the AAPM, ASTRO, ACR, ESTRO, and IAEA. These prescriptive guidelines have traditionally focused on monitoring all aspects of the functional performance of radiotherapy (RT) equipment by comparing parameters against tolerances set at strict but achievable values. Many errors that occur in radiation oncology are not due to failures in devices and software; rather they are failures in workflow and process. A systematic understanding of the likelihood and clinical impact of possible failures throughout a course of radiotherapy is needed to direct limit QM resources efficiently to produce maximum safety and quality of patient care. Task Group 100 of the AAPM has taken a broad view of these issues and has developed a framework for designing QM activities, based on estimates of the probability of identified failures and their clinical outcome through the RT planning…

PURPOSE: Accurate segmentation of organs-at-risks (OARs) is the key step for efficient planning of radiation therapy for head and neck (HaN) cancer treatment. In the work, we proposed the first deep learning-based algorithm, for segmentation of OARs in HaN CT images, and compared its performance against state-of-the-art automated segmentation algorithms, commercial software, and interobserver variability. METHODS: Convolutional neural networks (CNNs)-a concept from the field of deep learning-were used to study consistent intensity patterns of OARs from training CT images and to segment the OAR in a previously unseen test CT image. For CNN training, we extracted a representative number of positive intensity patches around voxels that belong to the OAR of interest in training CT images, and negative intensity patches around voxels that belong to the surrounding structures. These patches then passed through a sequence of CNN layers that captured local image features such as corners, end-points, and edges, and combined them into more complex high-order features that can efficiently…

PURPOSE: Develop a computer-aided detection (CAD) system for masses in digital breast tomosynthesis (DBT) volume using a deep convolutional neural network (DCNN) with transfer learning from mammograms. METHODS: A data set containing 2282 digitized film and digital mammograms and 324 DBT volumes were collected with IRB approval. The mass of interest on the images was marked by an experienced breast radiologist as reference standard. The data set was partitioned into a training set (2282 mammograms with 2461 masses and 230 DBT views with 228 masses) and an independent test set (94 DBT views with 89 masses). For DCNN training, the region of interest (ROI) containing the mass (true positive) was extracted from each image. False positive (FP) ROIs were identified at prescreening by their previously developed CAD systems. After data augmentation, a total of 45 072 mammographic ROIs and 37 450 DBT ROIs were obtained. Data normalization and reduction of non-uniformity in the ROIs across heterogeneous data was achieved using a background correction method applied to each ROI. A DCNN with…

PURPOSE: A computational toolkit (spektr 3.0) has been developed to calculate x-ray spectra based on the tungsten anode spectral model using interpolating cubic splines (TASMICS) algorithm, updating previous work based on the tungsten anode spectral model using interpolating polynomials (TASMIP) spectral model. The toolkit includes a matlab (The Mathworks, Natick, MA) function library and improved user interface (UI) along with an optimization algorithm to match calculated beam quality with measurements. METHODS: The spektr code generates x-ray spectra (photons/mm(2)/mAs at 100 cm from the source) using TASMICS as default (with TASMIP as an option) in 1 keV energy bins over beam energies 20-150 kV, extensible to 640 kV using the TASMICS spectra. An optimization tool was implemented to compute the added filtration (Al and W) that provides a best match between calculated and measured x-ray tube output (mGy/mAs or mR/mAs) for individual x-ray tubes that may differ from that assumed in TASMICS or TASMIP and to account for factors such as anode angle. RESULTS: The median percent…

PURPOSE: The GE SIGNA PET/MR is a new whole body integrated time-of-flight (ToF)-PET/MR scanner from GE Healthcare. The system is capable of simultaneous PET and MR image acquisition with sub-400 ps coincidence time resolution. Simultaneous PET/MR holds great potential as a method of interrogating molecular, functional, and anatomical parameters in clinical disease in one study. Despite the complementary imaging capabilities of PET and MRI, their respective hardware tends to be incompatible due to mutual interference. In this work, the GE SIGNA PET/MR is evaluated in terms of PET performance and the potential effects of interference from MRI operation. METHODS: The NEMA NU 2-2012 protocol was followed to measure PET performance parameters including spatial resolution, noise equivalent count rate, sensitivity, accuracy, and image quality. Each of these tests was performed both with the MR subsystem idle and with continuous MR pulsing for the duration of the PET data acquisition. Most measurements were repeated at three separate test sites where the system is installed. RESULTS: The…

PURPOSE: The authors are developing a computerized system for bladder segmentation in CT urography (CTU) as a critical component for computer-aided detection of bladder cancer. METHODS: A deep-learning convolutional neural network (DL-CNN) was trained to distinguish between the inside and the outside of the bladder using 160 000 regions of interest (ROI) from CTU images. The trained DL-CNN was used to estimate the likelihood of an ROI being inside the bladder for ROIs centered at each voxel in a CTU case, resulting in a likelihood map. Thresholding and hole-filling were applied to the map to generate the initial contour for the bladder, which was then refined by 3D and 2D level sets. The segmentation performance was evaluated using 173 cases: 81 cases in the training set (42 lesions, 21 wall thickenings, and 18 normal bladders) and 92 cases in the test set (43 lesions, 36 wall thickenings, and 13 normal bladders). The computerized segmentation accuracy using the DL likelihood map was compared to that using a likelihood map generated by Haar features and a random forest classifier,…

PURPOSE: Automated detection of solitary pulmonary nodules using positron emission tomography (PET) and computed tomography (CT) images shows good sensitivity; however, it is difficult to detect nodules in contact with normal organs, and additional efforts are needed so that the number of false positives (FPs) can be further reduced. In this paper, the authors propose an improved FP-reduction method for the detection of pulmonary nodules in PET/CT images by means of convolutional neural networks (CNNs). METHODS: The overall scheme detects pulmonary nodules using both CT and PET images. In the CT images, a massive region is first detected using an active contour filter, which is a type of contrast enhancement filter that has a deformable kernel shape. Subsequently, high-uptake regions detected by the PET images are merged with the regions detected by the CT images. FP candidates are eliminated using an ensemble method; it consists of two feature extractions, one by shape/metabolic feature analysis and the other by a CNN, followed by a two-step classifier, one step being rule based…

PURPOSE: Automated segmentation of breast and fibroglandular tissue (FGT) is required for various computer-aided applications of breast MRI. Traditional image analysis and computer vision techniques, such atlas, template matching, or, edge and surface detection, have been applied to solve this task. However, applicability of these methods is usually limited by the characteristics of the images used in the study datasets, while breast MRI varies with respect to the different MRI protocols used, in addition to the variability in breast shapes. All this variability, in addition to various MRI artifacts, makes it a challenging task to develop a robust breast and FGT segmentation method using traditional approaches. Therefore, in this study, we investigated the use of a deep-learning approach known as "U-net." MATERIALS AND METHODS: We used a dataset of 66 breast MRI's randomly selected from our scientific archive, which includes five different MRI acquisition protocols and breasts from four breast density categories in a balanced distribution. To prepare reference segmentations, we…

Purpose The aim of this study was to assess the suitability of Gafchromic EBT3 films for reference dose measurements in the beam of a prototype high dose‐per‐pulse linear accelerator (linac), capable of delivering electron beams with a mean dose‐rate (Ḋm) ranging from 0.07 to 3000 Gy/s and a dose‐rate in pulse (Ḋp) of up to 8 × 106 Gy/s. To do this, we evaluated the overall uncertainties in EBT3 film dosimetry as well as the energy and dose‐rate dependence of their response. Material and methods Our dosimetric system was composed of EBT3 Gafchromic films in combination with a flatbed scanner and was calibrated against an ionization chamber traceable to primary standard. All sources of uncertainties in EBT3 dosimetry were carefully analyzed using irradiations at a clinical radiotherapy linac. Energy dependence was investigated with the same machine by acquiring and comparing calibration curves for three different beam energies (4, 8 and 12 MeV), for doses between 0.25 and 30 Gy. Ḋm dependence was studied at the clinical linac by changing the pulse repetition frequency (f) of the beam…

PURPOSE: It is common practice to perform patient-specific pretreatment verifications to the clinical delivery of IMRT. This process can be time-consuming and not altogether instructive due to the myriad sources that may produce a failing result. The purpose of this study was to develop an algorithm capable of predicting IMRT QA passing rates a priori. METHODS: From all treatment, 498 IMRT plans sites were planned in eclipse version 11 and delivered using a dynamic sliding window technique on Clinac iX or TrueBeam Linacs. 3%/3 mm local dose/distance-to-agreement (DTA) was recorded using a commercial 2D diode array. Each plan was characterized by 78 metrics that describe different aspects of their complexity that could lead to disagreements between the calculated and measured dose. A Poisson regression with Lasso regularization was trained to learn the relation between the plan characteristics and each passing rate. RESULTS: Passing rates 3%/3 mm local dose/DTA can be predicted with an error smaller than 3% for all plans analyzed. The most important metrics to describe the passing…

Peer-reviewed journals and conference proceedings publish hundreds of papers that describe new medical imaging algorithms, including, for example, techniques for computer-aided diagnosis or characterization, segmentation, image registration, image reconstruction, and radiomics. It is difficult, if not impossible, to fairly compare the performance of these algorithms as investigators must either use different data sets, or if using the same data, use different implementations of competing algorithms. Grand Challenges facilitate the fair comparison of algorithms by providing a common data set to all participants and by having each participant be responsible for implementation of their own algorithm. The dissemination of findings from Grand Challenges provides important information to the scientific community and helps to determine which approaches have the greatest promise for successful translation to clinical practice. In this session we will review the outcomes and lessons learned from the 2015 SPIE-AAPM-NCI Lung Nodule Classification Challenge. We will then turn to the 2016…

PURPOSE: Magnetic resonance imaging-guided radiotherapy (MRIgRT) provides superior soft-tissue contrast and real-time imaging compared with standard image-guided RT, which uses x-ray based imaging. Several groups are developing integrated MRIgRT machines. Reference dosimetry with these new machines requires accounting for the effects of the magnetic field on the response of the ionization chambers used for dose calibration. Here, the authors propose a formalism for reference dosimetry with integrated MRIgRT devices. The authors also examined the suitability of the TPR10 (20) and %dd(10)x beam quality specifiers in the presence of magnetic fields and calculated detector correction factors to account for the effects of the magnetic field for a range of detectors. METHODS: The authors used full-head and point-source Monte Carlo models of an MR-linac along with detailed detector models of an Exradin A19, an NE2571, and several PTW Farmer chambers to calculate magnetic field correction factors for six commercial ionization chambers in three chamber configurations. Calculations of…

PURPOSE: Imaging biomarker research focuses on discovering relationships between radiological features and histological findings. In glioblastoma patients, methylation of the O(6)-methylguanine methyltransferase (MGMT) gene promoter is positively correlated with an increased effectiveness of current standard of care. In this paper, the authors investigate texture features as potential imaging biomarkers for capturing the MGMT methylation status of glioblastoma multiforme (GBM) tumors when combined with supervised classification schemes. METHODS: A retrospective study of 155 GBM patients with known MGMT methylation status was conducted. Co-occurrence and run length texture features were calculated, and both support vector machines (SVMs) and random forest classifiers were used to predict MGMT methylation status. RESULTS: The best classification system (an SVM-based classifier) had a maximum area under the receiver-operating characteristic (ROC) curve of 0.85 (95% CI: 0.78-0.91) using four texture features (correlation, energy, entropy, and local intensity) originating from the…

PURPOSE: This paper describes the design, fabrication, and characterization of multilayered tissue mimicking skin and vessel phantoms with tunable mechanical, optical, and acoustic properties. The phantoms comprise epidermis, dermis, and hypodermis skin layers, blood vessels, and blood mimicking fluid. Each tissue component may be individually tailored to a range of physiological and demographic conditions. METHODS: The skin layers were constructed from varying concentrations of gelatin and agar. Synthetic melanin, India ink, absorbing dyes, and Intralipid were added to provide optical absorption and scattering in the skin layers. Bovine serum albumin was used to increase acoustic attenuation, and 40 μm diameter silica microspheres were used to induce acoustic backscatter. Phantom vessels consisting of thin-walled polydimethylsiloxane tubing were embedded at depths of 2-6 mm beneath the skin, and blood mimicking fluid was passed through the vessels. The phantoms were characterized through uniaxial compression and tension experiments, rheological frequency sweep studies, diffuse…

Attenuation correction is an essential component of the long chain of data correction techniques required to achieve the full potential of quantitative positron emission tomography (PET) imaging. The development of combined PET/magnetic resonance imaging (MRI) systems mandated the widespread interest in developing novel strategies for deriving accurate attenuation maps with the aim to improve the quantitative accuracy of these emerging hybrid imaging systems. The attenuation map in PET/MRI should ideally be derived from anatomical MR images; however, MRI intensities reflect proton density and relaxation time properties of biological tissues rather than their electron density and photon attenuation properties. Therefore, in contrast to PET/computed tomography, there is a lack of standardized global mapping between the intensities of MRI signal and linear attenuation coefficients at 511 keV. Moreover, in standard MRI sequences, bones and lung tissues do not produce measurable signals owing to their low proton density and short transverse relaxation times. MR images are also inevitably…

PURPOSE: Ultrasound neuromodulation is a promising noninvasive technique for controlling neural activity. Previous small animal studies suffered from low targeting specificity because of the low ultrasound frequencies (<690 kHz) used. In this study, the authors demonstrated the capability of focused ultrasound (FUS) neuromodulation in the megahertz-range to achieve superior targeting specificity in the murine brain as well as demonstrate modulation of both motor and sensory responses. METHODS: FUS sonications were carried out at 1.9 MHz with 50% duty cycle, pulse repetition frequency of 1 kHz, and duration of 1 s. The robustness of the FUS neuromodulation was assessed first in sensorimotor cortex, where elicited motor activities were observed and recorded on videos and electromyography. Deeper brain regions were then targeted where pupillary dilation served as an indicative of successful modulation of subcortical brain structures. RESULTS: Contralateral and ipsilateral movements of the hind limbs were repeatedly observed when the FUS was targeted at the sensorimotor cortex. Induced…

PURPOSE: Accuracy in proton therapy treatment planning can be improved using Monte Carlo (MC) simulations. However the long computation time of such methods hinders their use in clinical routine. This work aims to develop a fast multipurpose Monte Carlo simulation tool for proton therapy using massively parallel central processing unit (CPU) architectures. METHODS: A new Monte Carlo, called MCsquare (many-core Monte Carlo), has been designed and optimized for the last generation of Intel Xeon processors and Intel Xeon Phi coprocessors. These massively parallel architectures offer the flexibility and the computational power suitable to MC methods. The class-II condensed history algorithm of MCsquare provides a fast and yet accurate method of simulating heavy charged particles such as protons, deuterons, and alphas inside voxelized geometries. Hard ionizations, with energy losses above a user-specified threshold, are simulated individually while soft events are regrouped in a multiple scattering theory. Elastic and inelastic nuclear interactions are sampled from ICRU 63 differential…

PURPOSE: To develop a semiautomated computer-aided diagnosis (cad) system for thyroid cancer using two-dimensional ultrasound images that can be used to yield a second opinion in the clinic to differentiate malignant and benign lesions. METHODS: A total of 118 ultrasound images that included axial and longitudinal images from patients with biopsy-confirmed malignant (n = 30) and benign (n = 29) nodules were collected. Thyroid cad software was developed to extract quantitative features from these images based on thyroid nodule segmentation in which adaptive diffusion flow for active contours was used. Various features, including histogram, intensity differences, elliptical fit, gray-level co-occurrence matrixes, and gray-level run-length matrixes, were evaluated for each region imaged. Based on these imaging features, a support vector machine (SVM) classifier was used to differentiate benign and malignant nodules. Leave-one-out cross-validation with sequential forward feature selection was performed to evaluate the overall accuracy of this method. Additionally, analyses with…

Purpose This work aims at investigating intensity corrected cone‐beam x‐ray computed tomography (CBCT) images for accurate dose calculation in adaptive intensity modulated proton therapy (IMPT) for prostate and head and neck (H&N) cancer. A deformable image registration (DIR)‐based method and a scatter correction approach using the image data obtained from DIR as prior are characterized and compared on the basis of the same clinical patient cohort for the first time. Methods Planning CT (pCT) and daily CBCT data (reconstructed images and measured projections) of four H&N and four prostate cancer patients have been considered in this study. A previously validated Morphons algorithm was used for DIR of the planning CT to the current CBCT image, yielding a so‐called virtual CT (vCT). For the first time, this approach was translated from H&N to prostate cancer cases in the scope of proton therapy. The warped pCT images were also used as prior for scatter correction of the CBCT projections for both tumor sites. Single field uniform dose and IMPT (only for H&N cases) treatment plans have…

PURPOSE: Detection (diagnosis) of diabetic retinopathy (DR) in optical coherence tomography (OCT) images for patients with type 2 diabetes, but almost clinically normal retina appearances. METHODS: The proposed computer-aided diagnostic (CAD) system detects the DR in three steps: (a) localizing and segmenting 12 distinct retinal layers on the OCT image; (b) deriving features of the segmented layers, and (c) learning most discriminative features and classifying each subject as normal or diabetic. To localise and segment the retinal layers, signals (intensities) of the OCT image are described with a joint Markov-Gibbs random field (MGRF) model of intensities and shape descriptors. Each segmented layer is characterized with cumulative probability distribution functions (CDF) of its locally extracted features, such as reflectivity, curvature, and thickness. A multistage deep fusion classification network (DFCN) with a stack of non-negativity-constrained autoencoders (NCAE) is trained to select the most discriminative retinal layers' features and use their CDFs for detecting the DR. A…

PURPOSE: To describe in detail a dataset consisting of serial four-dimensional computed tomography (4DCT) and 4D cone beam CT (4DCBCT) images acquired during chemoradiotherapy of 20 locally advanced, nonsmall cell lung cancer patients we have collected at our institution and shared publicly with the research community. ACQUISITION AND VALIDATION METHODS: As part of an NCI-sponsored research study 82 4DCT and 507 4DCBCT images were acquired in a population of 20 locally advanced nonsmall cell lung cancer patients undergoing radiation therapy. All subjects underwent concurrent radiochemotherapy to a total dose of 59.4-70.2 Gy using daily 1.8 or 2 Gy fractions. Audio-visual biofeedback was used to minimize breathing irregularity during all fractions, including acquisition of all 4DCT and 4DCBCT acquisitions in all subjects. Target, organs at risk, and implanted fiducial markers were delineated by a physician in the 4DCT images. Image coordinate system origins between 4DCT and 4DCBCT were manipulated in such a way that the images can be used to simulate initial patient setup in the…

PURPOSE: Modern facilities for actively scanned ion beam radiotherapy allow in principle the use of helium beams, which could present specific advantages, especially for pediatric tumors. In order to assess the potential use of these beams for radiotherapy, i.e., to create realistic treatment plans, the authors set up a dedicated (4)He beam model, providing base data for their treatment planning system TRiP98, and they have reported that in this work together with its physical and biological validations. METHODS: A semiempirical beam model for the physical depth dose deposition and the production of nuclear fragments was developed and introduced in TRiP98. For the biological effect calculations the last version of the local effect model was used. The model predictions were experimentally verified at the HIT facility. The primary beam attenuation and the characteristics of secondary charged particles at various depth in water were investigated using (4)He ion beams of 200 MeV/u. The nuclear charge of secondary fragments was identified using a ΔE/E telescope. 3D absorbed dose…

PURPOSE: To demonstrate the preliminary feasibility of a longitudinal diffusion magnetic resonance imaging (MRI) strategy for assessing patient response to radiotherapy at 0.35 T using an MRI-guided radiotherapy system (ViewRay). METHODS: Six patients (three head and neck cancer, three sarcoma) who underwent fractionated radiotherapy were enrolled in this study. A 2D multislice spin echo single-shot echo planar imaging diffusion pulse sequence was implemented on the ViewRay system and tested in phantom studies. The same pulse sequence was used to acquire longitudinal diffusion data (every 2-5 fractions) on the six patients throughout the entire course of radiotherapy. The reproducibility of the apparent diffusion coefficient (ADC) measurements was assessed using reference regions and the temporal variations of the tumor ADC values were evaluated. RESULTS: In diffusion phantom studies, the ADC values measured on the ViewRay system matched well with reference ADC values with <5% error for a range of ground truth diffusion coefficients of 0.4-1.1 × 10(-3) mm(2)/s. The remote reference…

PURPOSE: Dual energy CT (DECT) has recently been proposed as an improvement over single energy CT (SECT) for stopping power ratio (SPR) estimation for proton therapy treatment planning (TP), thereby potentially reducing range uncertainties. Published literature investigated phantoms. This study aims at performing proton therapy TP on SECT and DECT head images of the same patients and at evaluating whether the reported improved DECT SPR accuracy translates into clinically relevant range shifts in clinical head treatment scenarios. METHODS: Two phantoms were scanned at a last generation dual source DECT scanner at 90 and 150 kVp with Sn filtration. The first phantom (Gammex phantom) was used to calibrate the scanner in terms of SPR while the second served as evaluation (CIRS phantom). DECT images of five head trauma patients were used as surrogate cancer patient images for TP of proton therapy. Pencil beam algorithm based TP was performed on SECT and DECT images and the dose distributions corresponding to the optimized proton plans were calculated using a Monte Carlo (MC) simulation…

PURPOSE: This paper provides a comparison between a fast, commercial, in-patient Monte Carlo dose calculation algorithm (GPUMCD) and geant4. It also evaluates the dosimetric impact of the application of an external 1.5 T magnetic field. METHODS: A stand-alone version of the Elekta™ GPUMCD algorithm, to be used within the Monaco treatment planning system to model dose for the Elekta™ magnetic resonance imaging (MRI) Linac, was compared against GEANT4 (v10.1). This was done in the presence or absence of a 1.5 T static magnetic field directed orthogonally to the radiation beam axis. Phantoms with material compositions of water, ICRU lung, ICRU compact-bone, and titanium were used for this purpose. Beams with 2 MeV monoenergetic photons as well as a 7 MV histogrammed spectrum representing the MRI Linac spectrum were emitted from a point source using a nominal source-to-surface distance of 142.5 cm. Field sizes ranged from 1.5 × 1.5 to 10 × 10 cm(2). Dose scoring was performed using a 3D grid comprising 1 mm(3) voxels. The production thresholds were equivalent for both codes. Results…

The problem of attenuation correction (AC) for quantitative positron emission tomography (PET) had been considered solved to a large extent after the commercial availability of devices combining PET with computed tomography (CT) in 2001; single photon emission computed tomography (SPECT) has seen a similar development. However, stimulated in particular by technical advances toward clinical systems combining PET and magnetic resonance imaging (MRI), research interest in alternative approaches for PET AC has grown substantially in the last years. In this comprehensive literature review, the authors first present theoretical results with relevance to simultaneous reconstruction of attenuation and activity. The authors then look back at the early history of this research area especially in PET; since this history is closely interwoven with that of similar approaches in SPECT, these will also be covered. We then review algorithmic advances in PET, including analytic and iterative algorithms. The analytic approaches are either based on the Helgason-Ludwig data consistency conditions of…

PURPOSE: The amount of coronary artery calcification (CAC) is a strong and independent predictor of cardiovascular disease (CVD) events. In clinical practice, CAC is manually identified and automatically quantified in cardiac CT using commercially available software. This is a tedious and time-consuming process in large-scale studies. Therefore, a number of automatic methods that require no interaction and semiautomatic methods that require very limited interaction for the identification of CAC in cardiac CT have been proposed. Thus far, a comparison of their performance has been lacking. The objective of this study was to perform an independent evaluation of (semi)automatic methods for CAC scoring in cardiac CT using a publicly available standardized framework. METHODS: Cardiac CT exams of 72 patients distributed over four CVD risk categories were provided for (semi)automatic CAC scoring. Each exam consisted of a noncontrast-enhanced calcium scoring CT (CSCT) and a corresponding coronary CT angiography (CCTA) scan. The exams were acquired in four different hospitals using…

Purpose: To assess the performance of using transferred features from pre-trained deep convolutional networks (CNNs) in the task of classifying cancer in breast ultrasound images, and to compare this method of transfer learning with previous methods involving human-designed features. Methods: A breast ultrasound dataset consisting of 1125 cases and 2393 regions of interest (ROIs) was used. Each ROI was labeled as cystic, benign, or malignant. Features were extracted from each ROI using pre-trained CNNs and used to train support vector machine (SVM) classifiers in the tasks of distinguishing non-malignant (benign+cystic) vs malignant lesions and benign vs malignant lesions. For a baseline comparison, classifiers were also trained on prior analytically-extracted tumor features. Five-fold cross-validation (by case) was conducted with the area under the receiver operating characteristic curve (AUC) as the performance metric. Results: Classifiers trained on CNN-extracted features were comparable to classifiers trained on human-designed features. In the non-malignant vs malignant task,…

PURPOSE: Patient motion is a common problem during dynamic positron emission tomography (PET) scans for quantification of myocardial blood flow (MBF). The purpose of this study was to quantify the prevalence of body motion in a clinical setting and evaluate with realistic phantoms the effects of motion on blood flow quantification, including CT attenuation correction (CTAC) artifacts that result from PET-CT misalignment. METHODS: A cohort of 236 sequential patients was analyzed for patient motion under resting and peak stress conditions by two independent observers. The presence of motion, affected time-frames, and direction of motion was recorded; discrepancy between observers was resolved by consensus review. Based on these results, patient body motion effects on MBF quantification were characterized using the digital NURBS-based cardiac-torso phantom, with characteristic time activity curves (TACs) assigned to the heart wall (myocardium) and blood regions. Simulated projection data were corrected for attenuation and reconstructed using filtered back-projection. All simulations…

PURPOSE: Magnetic resonance imaging (MRI) is increasingly used for radiotherapy target delineation, image guidance, and treatment response monitoring. Recent studies have shown that an entire external x-ray radiotherapy treatment planning (RTP) workflow for brain tumor or prostate cancer patients based only on MRI reference images is feasible. This study aims to show that a MRI-only based RTP workflow is also feasible for proton beam therapy plans generated in MRI-based substitute computed tomography (sCT) images of the head and the pelvis. METHODS: The sCTs were constructed for ten prostate cancer and ten brain tumor patients primarily by transforming the intensity values of in-phase MR images to Hounsfield units (HUs) with a dual model HU conversion technique to enable heterogeneous tissue representation. HU conversion models for the pelvis were adopted from previous studies, further extended in this study also for head MRI by generating anatomical site-specific conversion models (a new training data set of ten other brain patients). This study also evaluated two other types of…

PURPOSE: Modern radiotherapy modalities often use small or nonstandard fields to ensure highly localized and precise dose delivery, challenging conventional clinical dosimetry protocols. The emergence of preclinical spatially fractionated synchrotron radiotherapies with high dose-rate, sub-millimetric parallel kilovoltage x-ray beams, has pushed clinical dosimetry to its limit. A commercially available synthetic single crystal diamond detector designed for small field dosimetry has been characterized to assess its potential as a dosimeter for synchrotron microbeam and minibeam radiotherapy. METHODS: Experiments were carried out using a synthetic diamond detector on the imaging and medical beamline (IMBL) at the Australian Synchrotron. The energy dependence of the detector was characterized by cross-referencing with a calibrated ionization chamber in monoenergetic beams in the energy range 30-120 keV. The dose-rate dependence was measured in the range 1-700 Gy/s. Dosimetric quantities were measured in filtered white beams, with a weighted mean energy of 95 keV, in broadbeam and…

PURPOSE: In radiotherapy based only on magnetic resonance imaging (MRI), knowledge about tissue electron densities must be derived from the MRI. This can be achieved by converting the MRI scan to the so-called pseudo-computed tomography (pCT). An obstacle is that the voxel intensities in conventional MRI scans are not uniquely related to electron density. The authors previously demonstrated that a patch-based method could produce accurate pCTs of the brain using conventional T1-weighted MRI scans. The method was driven mainly by local patch similarities and relied on simple affine registrations between an atlas database of the co-registered MRI/CT scan pairs and the MRI scan to be converted. In this study, the authors investigate the applicability of the patch-based approach in the pelvis. This region is challenging for a method based on local similarities due to the greater inter-patient variation. The authors benchmark the method against a baseline pCT strategy where all voxels inside the body contour are assigned a water-equivalent bulk density. Furthermore, the authors implement…

Purpose Physical phantoms are central to the evaluation of 2D and 3D breast‐imaging systems. Currently, available physical phantoms have limitations including unrealistic uniform background structure, large expense, or excessive fabrication time. The purpose of this work is to outline a method for rapidly creating realistic, inexpensive physical anthropomorphic phantoms for use in full‐field digital mammography (FFDM) and digital breast tomosynthesis (DBT). Methods The phantom was first modeled using analytical expressions and then discretized into voxels of a specified size. The interior of the breast was divided into glandular and adipose tissue classes using Voronoi segmentation, and additional structures like blood vessels, chest muscle, and ligaments were added. The physical phantom was then fabricated from the virtual model in a slice by slice fashion through inkjet printing, using parchment paper and a radiopaque ink containing 33% (I33%) or 25% (I25%) iohexol by volume. Three types of parchment paper (P1, P2, and P3) were examined. The phantom materials were characterized in…

Purpose To determine whether individual liver tumor patients can be safely treated with pencil beam scanning proton therapy. This study reports a planning preparation workflow that can be used for beam angle selection and the evaluation of the efficacy of abdominal compression (AC) to mitigate motion. Methods Four‐dimensional computed tomography scans (4DCT) with and without AC were available from 10 liver tumor patients with fluoroscopy‐proven motion reduction by AC, previously treated using photons. For each scan, the motion amplitudes and the motion‐induced variation of water‐equivalent thickness (ΔWET) in each voxel of the target volume were evaluated during treatment plan preparation. Optimal proton beam angles were selected after volume analysis of the respective beam‐specific planning target volume (BSPTV). M⊥80 and ΔWET80 derived from the 80th percentiles of perpendicular motion amplitude (M⊥) and ΔWET were compared with and without AC. Proton plans were created on the average CT to achieve target coverage similar to that of the conventional photon treatments. 4D dynamic…

Purpose To develop and validate an automated technique for evaluating the spatial resolution characteristics of clinical computed tomography (CT) images. Methods Twenty one chest and abdominopelvic clinical CT datasets were examined in this study. An algorithm was developed to extract a CT resolution index (RI) analogous to the modulation transfer function from clinical CT images by measuring the edge‐spread function (ESF) across the patient's skin. A polygon mesh of the air‐skin boundary was created. The faces of the mesh were then used to measure the ESF across the air‐skin interface. The ESF was differentiated to obtain the line‐spread function (LSF), and the LSF was Fourier transformed to obtain the RI. The algorithm's ability to detect the radial dependence of the RI was investigated. RIs measured with the proposed method were compared with a conventional phantom‐based method across two reconstruction algorithms (FBP and iterative) using the spatial frequency at 50% RI, f50, as the metric for comparison. Three reconstruction kernels were investigated for each reconstruction…

PURPOSE: Computational models of microwave ablation (MWA) are widely used during the design optimization of novel devices and are under consideration for patient-specific treatment planning. The objective of this study was to assess the sensitivity of computational models of MWA to tissue biophysical properties. METHODS: The Morris method was employed to assess the global sensitivity of the coupled electromagnetic-thermal model, which was implemented with the finite element method (FEM). The FEM model incorporated temperature dependencies of tissue physical properties. The variability of the model was studied using six different outputs to characterize the size and shape of the ablation zone, as well as impedance matching of the ablation antenna. Furthermore, the sensitivity results were statistically analyzed and absolute influence of each input parameter was quantified. A framework for systematically incorporating model uncertainties for treatment planning was suggested. RESULTS: A total of 1221 simulations, incorporating 111 randomly sampled starting points, were performed.…

PURPOSE: To develop methods for evaluation and mitigation of dosimetric impact due to respiratory and diaphragmatic motion during free breathing in treatment of distal esophageal cancers using intensity-modulated proton therapy (IMPT). METHODS: This was a retrospective study on 11 patients with distal esophageal cancer. For each patient, four-dimensional computed tomography (4D CT) data were acquired, and a nominal dose was calculated on the average phase of the 4D CT. The changes of water equivalent thickness (ΔWET) to cover the treatment volume from the peak of inspiration to the valley of expiration were calculated for a full range of beam angle rotation. Two IMPT plans were calculated: one at beam angles corresponding to small ΔWET and one at beam angles corresponding to large ΔWET. Four patients were selected for the calculation of 4D-robustness-optimized IMPT plans due to large motion-induced dose errors generated in conventional IMPT. To quantitatively evaluate motion-induced dose deviation, the authors calculated the lowest dose received by 95% (D95) of the internal clinical…

PURPOSE: The purpose of this work was to develop, validate, and compare a highly computer-aided method for the segmentation of hot lesions in head and neck 18F-FDG PET scans. METHODS: A semiautomated segmentation method was developed, which transforms the segmentation problem into a graph-based optimization problem. For this purpose, a graph structure around a user-provided approximate lesion centerpoint is constructed and a suitable cost function is derived based on local image statistics. To handle frequently occurring situations that are ambiguous (e.g., lesions adjacent to each other versus lesion with inhomogeneous uptake), several segmentation modes are introduced that adapt the behavior of the base algorithm accordingly. In addition, the authors present approaches for the efficient interactive local and global refinement of initial segmentations that are based on the "just-enough-interaction" principle. For method validation, 60 PET/CT scans from 59 different subjects with 230 head and neck lesions were utilized. All patients had squamous cell carcinoma of the head and neck.…

PURPOSE: The accuracy of proton dose calculation is dependent on the ability to correctly characterize patient tissues with medical imaging. The most common method is to correlate computed tomography (CT) numbers obtained via single-energy CT (SECT) with proton stopping power ratio (SPR). CT numbers, however, cannot discriminate between a change in mass density and change in chemical composition of patient tissues. This limitation can have consequences on SPR calibration accuracy. Dual-energy CT (DECT) is receiving increasing interest as an alternative imaging modality for proton therapy treatment planning due to its ability to discriminate between changes in patient density and chemical composition. In the current work we use a phantom of known composition to demonstrate the dosimetric advantages of proton therapy treatment planning with DECT over SECT. METHODS: A phantom of known composition was scanned with a clinical SECT radiotherapy CT-simulator. The phantom was rescanned at a lower X-ray tube potential to generate a complimentary DECT image set. A set of reference materials…

PURPOSE: White matter hyperintensities (WMH) are seen on FLAIR-MRI in several neurological disorders, including multiple sclerosis, dementia, Parkinsonism, stroke and cerebral small vessel disease (SVD). WMHs are often used as biomarkers for prognosis or disease progression in these diseases, and additionally longitudinal quantification of WMHs is used to evaluate therapeutic strategies. Human readers show considerable disagreement and inconsistency on detection of small lesions. A multitude of automated detection algorithms for WMHs exists, but since most of the current automated approaches are tuned to optimize segmentation performance according to Jaccard or Dice scores, smaller WMHs often go undetected in these approaches. In this paper, the authors propose a method to accurately detect all WMHs, large as well as small. METHODS: A two-stage learning approach was used to discriminate WMHs from normal brain tissue. Since small and larger WMHs have quite a different appearance, the authors have trained two probabilistic classifiers: one for the small WMHs (⩽3 mm effective diameter)…

PURPOSE: To accurately and efficiently implement charged particle transport in a magnetic field in EGSnrc and validate the code for the use in phantom and ion chamber simulations. METHODS: The effect of the magnetic field on the particle motion and position is determined using one- and three-point numerical integrations of the Lorentz force on the charged particle and is added to the condensed history calculation performed by the EGSnrc PRESTA-II algorithm. The code is tested with a Fano test adapted for the presence of magnetic fields. The code is compatible with all EGSnrc based applications, including egs++. Ion chamber calculations are compared to experimental measurements and the effect of the code on the efficiency and timing is determined. RESULTS: Agreement with the Fano test's theoretical value is obtained at the 0.1% level for large step-sizes and in magnetic fields as strong as 5 T. The NE2571 dose calculations achieve agreement with the experiment within 0.5% up to 1 T beyond which deviations up to 1.2% are observed. Uniform air gaps of 0.5 and 1 mm and a misalignment of…

PURPOSE: Due to the increased axial coverage of multislice computed tomography (CT) and the introduction of flat detectors, the size of x-ray illumination fields has grown dramatically, causing an increase in scatter radiation. For CT imaging, scatter is a significant issue that introduces shading artifact, streaks, as well as reduced contrast and Hounsfield Units (HU) accuracy. The purpose of this work is to provide a fast and accurate scatter artifacts correction algorithm for cone beam CT (CBCT) imaging. METHODS: The method starts with an estimation of coarse scatter profiles for a set of CBCT data in either image domain or projection domain. A denoising algorithm designed specifically for Poisson signals is then applied to derive the final scatter distribution. Qualitative and quantitative evaluations using thorax and abdomen phantoms with Monte Carlo (MC) simulations, experimental Catphan phantom data, and in vivo human data acquired for a clinical image guided radiation therapy were performed. Scatter correction in both projection domain and image domain was conducted and the…

Purpose: The authors previously reported successful luminescence imaging of water during proton irradiation and its application to range estimation. However, since the feasibility of this approach for carbon-ion irradiation remained unclear, the authors conducted luminescence imaging during carbon-ion irradiation and estimated the ranges. Methods: The authors placed a pure-water phantom on the patient couch of a carbon-ion therapy system and measured the luminescence images with a high-sensitivity, cooled charge-coupled device camera during carbon-ion irradiation. The authors also carried out imaging of three types of phantoms (tap-water, an acrylic block, and a plastic scintillator) and compared their intensities and distributions with those of a phantom containing pure-water. Results: The luminescence images of pure-water phantoms during carbon-ion irradiation showed clear Bragg peaks, and the measured carbon-ion ranges from the images were almost the same as those obtained by simulation. The image of the tap-water phantom showed almost the same distribution as that of the…

PURPOSE: The FDA approved the use of digital breast tomosynthesis (DBT) in 2011 as an adjunct to 2D full field digital mammography (FFDM) with the constraint that all DBT acquisitions must be paired with a 2D image to assure adequate interpretative information is provided. Recently manufacturers have developed methods to provide a synthesized 2D image generated from the DBT data with the hope of sparing patients the radiation exposure from the FFDM acquisition. While this much needed alternative effectively reduces the total radiation burden, differences in image quality must also be considered. The goal of this study was to compare the intrinsic image quality of synthesized 2D c-view and 2D FFDM images in terms of resolution, contrast, and noise. METHODS: Two phantoms were utilized in this study: the American College of Radiology mammography accreditation phantom (ACR phantom) and a novel 3D printed anthropomorphic breast phantom. Both phantoms were imaged using a Hologic Selenia Dimensions 3D system. Analysis of the ACR phantom includes both visual inspection and objective…

PURPOSE: To allow for a purely image-based motion estimation and compensation in weight-bearing cone-beam computed tomography of the knee joint. METHODS: Weight-bearing imaging of the knee joint in a standing position poses additional requirements for the image reconstruction algorithm. In contrast to supine scans, patient motion needs to be estimated and compensated. The authors propose a method that is based on 2D/3D registration of left and right femur and tibia segmented from a prior, motion-free reconstruction acquired in supine position. Each segmented bone is first roughly aligned to the motion-corrupted reconstruction of a scan in standing or squatting position. Subsequently, a rigid 2D/3D registration is performed for each bone to each of K projection images, estimating 6 × 4 × K motion parameters. The motion of individual bones is combined into global motion fields using thin-plate-spline extrapolation. These can be incorporated into a motion-compensated reconstruction in the backprojection step. The authors performed visual and quantitative comparisons between a…

PURPOSE: The new radiochromic film, GAFChromic EBT-XD, contains the same active material, lithium-10,12-pentacosadiynoate, as GAFChromic EBT3, but the crystalline form is different. This work investigates the effect of this change on the well-known lateral response artifact when EBT-XD film is digitized on a flatbed scanner. METHODS: The dose response of a single production lot of EBT-XD was characterized by scanning an unexposed film plus a set of films exposed to doses between 2.5 and 50 Gy using 6 MV photons. To characterize the lateral response artifact, the authors used the unexposed film plus a subset of samples exposed to doses between 20 and 50 Gy. Digital images of these films were acquired at seven discrete lateral locations perpendicular to the scan direction on three Epson 10000XL scanners. Using measurements at the discrete lateral positions, the scanner responses were determined as a function of the lateral position of the film. From the data for each scanner, a set of coefficients were derived whereby measured response values could be corrected to remove the effects…

PURPOSE: Proton computed tomography (pCT) will enable accurate prediction of proton and ion range in a patient while providing the benefit of lower radiation exposure than in x-ray CT. The accuracy of the range prediction is essential for treatment planning in proton or ion therapy and depends upon the detector used to evaluate the water-equivalent path length (WEPL) of a proton passing through the object. A novel approach is presented for an inexpensive WEPL detector for pCT and proton radiography. METHODS: A novel multistage detector with an aperture of 10 × 37.5 cm was designed to optimize the accuracy of the WEPL measurements while simplifying detector construction and the performance requirements of its components. The design of the five-stage detector was optimized through simulations based on the geant4 detector simulation toolkit, and the fabricated prototype was calibrated in water-equivalent millimeters with 200 MeV protons in the research beam line of the clinical proton synchrotron at Loma Linda University Medical Center. A special polystyrene step phantom was designed…

PURPOSE: To examine the abilities of a traditional failure mode and effects analysis (FMEA) and modified healthcare FMEA (m-HFMEA) scoring methods by comparing the degree of congruence in identifying high risk failures. METHODS: The authors applied two prospective methods of the quality management to surface image guided, linac-based radiosurgery (SIG-RS). For the traditional FMEA, decisions on how to improve an operation were based on the risk priority number (RPN). The RPN is a product of three indices: occurrence, severity, and detectability. The m-HFMEA approach utilized two indices, severity and frequency. A risk inventory matrix was divided into four categories: very low, low, high, and very high. For high risk events, an additional evaluation was performed. Based upon the criticality of the process, it was decided if additional safety measures were needed and what they comprise. RESULTS: The two methods were independently compared to determine if the results and rated risks matched. The authors' results showed an agreement of 85% between FMEA and m-HFMEA approaches for top 20…

PURPOSE: A protocol for reference dosimetry for the MR-linac is under development. The 1.5 T magnetic field changes the mean path length of electrons in an air-filled ionization chamber but has little effect on the electron trajectories in a surrounding phantom. It is therefore necessary to correct the response of an ionization chamber for the influence of the magnetic field. Solid phantoms are used for dosimetry measurements on the MR-linac, but air is present between the chamber wall and phantom insert. This study aimed to determine if this air influences the ion chamber measurements on the MR-linac. The absolute response of the chamber and reproducibility of dosimetry measurements were assessed on an MR-linac in solid and water phantoms. The sensitivity of the chamber response to the distribution of air around the chamber was also investigated. METHODS: Measurements were performed on an MR-linac and replicated on a conventional linac for five chambers. The response of three waterproof chambers was measured with air and with water between the chamber and the insert to measure the…

PURPOSE: To evaluate the accuracy and robustness of a simple, linear, separable, two-parameter model (basis vector model, BVM) in mapping proton stopping powers via dual energy computed tomography (DECT) imaging. METHODS: The BVM assumes that photon cross sections (attenuation coefficients) of unknown materials are linear combinations of the corresponding radiological quantities of dissimilar basis substances (i.e., polystyrene, CaCl2 aqueous solution, and water). The authors have extended this approach to the estimation of electron density and mean excitation energy, which are required parameters for computing proton stopping powers via the Bethe-Bloch equation. The authors compared the stopping power estimation accuracy of the BVM with that of a nonlinear, nonseparable photon cross section Torikoshi parametric fit model (VCU tPFM) as implemented by the authors and by Yang et al. ["Theoretical variance analysis of single- and dual-energy computed tomography methods for calculating proton stopping power ratios of biological tissues," Phys. Med. Biol. 55, 1343-1362 (2010)]. Using an…

PURPOSE: The Leksell Gamma Knife (LGK) Icon has been recently introduced to provide Gamma Knife technology with frameless stereotactic treatments which use an additional cone-beam CT (CBCT) imaging system and a motion tracking system (IFMM, Intra-Fraction Motion Management). The system was commissioned for the treatment unit itself as well as the imaging system. METHODS: The LGK Icon was calibrated using an A1SL ionization chamber. EBT3 radiochromic films were employed to independently check the machine calibration, to measure the relative output factors (ROFs) and to collect dose distributions. Coincidence between CBCT isocenter and radiological focus was evaluated by means of EBT3 films. CBCT image quality was investigated in terms of spatial resolution, contrast-to-noise ratio (CNR), and uniformity for the two presets available (low dose and high dose). Computed Tomography Dose Index (CTDI) was also measured for both presets. RESULTS: The absolute dose rate of the LGK Icon was 3.86 ± 0.09 Gy/min. This result was confirmed by EBT3 readings. ROF were found to be 0.887 ± 0.035 and…

Purpose: Magnetic resonance imaging (MRI)‐only radiotherapy treatment planning requires accurate pseudo‐CT (pCT) images for precise dose calculation. The current work introduced an atlas‐based method combined with MR intensity information. pCT analyses and Monte Carlo dose calculations for intracranial stereotactic treatments were performed. Methods: Twenty‐two patients, representing 35 tumor targets, were scanned using a 3D T 1‐weighted MRI sequence according to the clinical protocol. The MR atlas image was registered to the MR patient image using a deformable algorithm, and the deformation was then applied to the atlas CT. Two methods were applied. The first method (MRdef) was based on deformations only, while the second (MRint) also used the actual MR intensities. pCT analysis was performed using the mean (absolute) error, as well as an in‐house tool based on a gamma index. Dose differences between pCT and true CT were analyzed using dose–volume histogram (DVH) parameters, statistical tests, the gamma index, and probability density functions. An unusual case, where the patient…

PURPOSE: To analyze the most recent results of the Imaging and Radiation Oncology Core Houston Quality Assurance Center's (IROC-H) anthropomorphic head and neck (H&N) phantom to determine the nature of failing irradiations and the feasibility of altering credentialing criteria. METHODS: IROC-H's H&N phantom, used for intensity-modulated radiation therapy credentialing for National Cancer Institute-sponsored clinical trials, requires that an institution's treatment plan agrees within ±7% of measured thermoluminescent dosimeter (TLD) doses; it also requires that ≥85% of pixels pass ±4 mm distance to agreement (7%/4 mm gamma analysis for film). The authors re-evaluated 156 phantom irradiations (November 1, 2014-October 31, 2015) according to the following tighter criteria: (1) 5% TLD and 5%/4 mm, (2) 5% TLD and 5%/3 mm, (3) 4% TLD and 4%/4 mm, and (4) 3% TLD and 3%/3 mm. Failure rates were evaluated with respect to individual film and TLD performance by location in the phantom. Overall poor phantom results were characterized qualitatively as systematic errors (correct shape and…

PURPOSE: Mild hyperthermia can be used as an adjuvant therapy to enhance radiation therapy or chemotherapy of cancer. However, administering mild hyperthermia is technically challenging due to the high accuracy required of the temperature control. MR guided high-intensity focused ultrasound (MR-HIFU) is a technology that can address this challenge. In this work, accurate and spatially uniform mild hyperthermia is demonstrated for deep-seated clinically relevant heating volumes using a HIFU system under MR guidance. METHODS: Mild hyperthermia heating was evaluated for temperature accuracy and spatial uniformity in 11 in vivo porcine leg experiments. Hyperthermia was induced with a commercial Philips Sonalleve MR-HIFU system embedded in a 1.5T Ingenia MR scanner. The operating software was modified to allow extended duration mild hyperthermia. Heating time varied from 10 min up to 60 min and the assigned target temperature was 42.5 °C. Electronic focal point steering, mechanical transducer movement, and dynamic transducer element switch-off were exploited to enlarge the heated volume…

PURPOSE: To use novel voxel-based 3D printed textured phantoms in order to compare low-contrast detectability between two reconstruction algorithms, FBP (filtered-backprojection) and SAFIRE (sinogram affirmed iterative reconstruction) and determine what impact background texture (i.e., anatomical noise) has on estimating the dose reduction potential of SAFIRE. METHODS: Liver volumes were segmented from 23 abdominal CT cases. The volumes were characterized in terms of texture features from gray-level co-occurrence and run-length matrices. Using a 3D clustered lumpy background (CLB) model, a fitting technique based on a genetic optimization algorithm was used to find CLB textures that were reflective of the liver textures, accounting for CT system factors of spatial blurring and noise. With the modeled background texture as a guide, four cylindrical phantoms (Textures A-C and uniform, 165 mm in diameter, and 30 mm height) were designed, each containing 20 low-contrast spherical signals (6 mm diameter at nominal contrast levels of ∼3.2, 5.2, 7.2, 10, and 14 HU with four repeats per…

PURPOSE: Both humans and software are notoriously challenging to account for in traditional hazard analysis models. The purpose of this work is to investigate and demonstrate the application of a new, extended accident causality model, called systems theoretic accident model and processes (STAMP), to radiation oncology. Specifically, a hazard analysis technique based on STAMP, system-theoretic process analysis (STPA), is used to perform a hazard analysis. METHODS: The STPA procedure starts with the definition of high-level accidents for radiation oncology at the medical center and the hazards leading to those accidents. From there, the hierarchical safety control structure of the radiation oncology clinic is modeled, i.e., the controls that are used to prevent accidents and provide effective treatment. Using STPA, unsafe control actions (behaviors) are identified that can lead to the hazards as well as causal scenarios that can lead to the identified unsafe control. This information can be used to eliminate or mitigate potential hazards. The STPA procedure is demonstrated on a new…

PURPOSE: To develop a clinically feasible method for dose calculations on cone beam CT (CBCT) images of two different vendors, and to determine the accuracy of these dose calculations for lung cancer patients. METHODS: Lung cancer patients with CBCT imaging (n = 10 for Elekta, n = 6 for Varian) and a repeated planning CT scan on the same day were selected. For CBCT dose calculations, an adapted Hounsfield units-to-mass density table (HU table) was used which was obtained by comparing CT values of corresponding points on the CBCT and the repeated planning CT scan. Dose calculations with three different HU tables were compared: a patient-specific, a general thorax-CBCT, and the standard CT HU table. Planning CT data were used to compensate for the limited field of view (FOV) (Elekta) or scan length (Varian) of the CBCT. For evaluation, clinically relevant dose metrics were compared between the repeated CT and CBCT to assess the accuracy of dose calculations on CBCT for both vendors. RESULTS: For both vendors, isodose lines and dose volume histograms were very similar between dose…

PURPOSE: Small field output correction factors have been studied by several research groups for the PTW 60019 microDiamond (MD) dosimeter, by comparing the response of such a device with both reference dosimeters and Monte Carlo simulations. A general good agreement is observed for field sizes down to about 1 cm. However, evident inconsistencies can be noticed when comparing some experimental results and Monte Carlo simulations obtained for smaller irradiation fields. This issue was tentatively attributed by some authors to unintentional large variations of the MD active surface area. The aim of the present study is a nondestructive experimental determination of the MD active surface area and active volume. METHODS: Ten MD dosimeters, one MD prototype, and three synthetic diamond samples were investigated in the present work. 2D maps of the MD response were recorded under scanned soft x-ray microbeam irradiation, leading to an experimental determination of the device active surface area. Profiles of the device responses were measured as well. In order to evaluate the MD active…

Although a multicenter, Phase III, prospective, randomized trial is the gold standard for evidence-based medicine, it is rarely used in the evaluation of innovative devices because of many practical and ethical reasons. It is usually sufficient to compare the dose distributions and dose rates for determining the equivalence of the innovative treatment modality to an existing one. Thus, quantitative evaluation of the dosimetric characteristics of innovative radiotherapy devices or applications is a critical part in which physicists should be actively involved. The physicist's role, along with physician colleagues, in this process is highlighted for innovative brachytherapy devices and applications and includes evaluation of (1) dosimetric considerations for clinical implementation (including calibrations, dose calculations, and radiobiological aspects) to comply with existing societal dosimetric prerequisites for sources in routine clinical use, (2) risks and benefits from a regulatory and safety perspective, and (3) resource assessment and preparedness. Further, it is suggested that…

PURPOSE: Proton beam dosimetry using bare plastic optical fibers has emerged as a simple approach to proton beam dosimetry. The source of the signal in this method has been attributed to Čerenkov radiation. The aim of this work was a phenomenological study of the nature of the visible light responsible for the signal in bare fiber optic dosimetry of proton therapy beams. METHODS: Plastic fiber optic probes embedded in solid water phantoms were irradiated with proton beams of energies 100, 180, and 225 MeV produced by a proton therapy cyclotron. Luminescence spectroscopy was performed by a CCD-coupled spectrometer. The spectra were acquired at various depths in phantom to measure the percentage depth dose (PDD) for each beam energy. For comparison, the PDD curves were acquired using a standard multilayer ion chamber device. In order to further analyze the contribution of the Čerenkov radiation in the spectra, Monte Carlo simulation was performed using fluka Monte Carlo code to stochastically simulate radiation transport, ionizing radiation dose deposition, and optical emission of…

PURPOSE: The objective of this study is to measure through simulation the impact of (1) heterogeneity of biophysical parameters in tumor vs healthy tissue, (2) applicator placement relative to the tumor, and (3) proximity to large blood vessels on microwave ablation (MWA) treatment effect area. This will help identify the biophysical properties that have the greatest impact on improving clinical modeling of MWA procedures. METHODS: The authors' approach was to develop two-compartment models with variable tissue properties and simulate MWA procedures performed in liver with Perseon Medical's 915 MHz short-tip applicator. Input parameters for the dielectric and thermal properties considered in this study were based on measurements for healthy and malignant (primary or metastatic) liver tissue previously reported in the literature. Compartment 1 (C1) represented normal, fatty, or cirrhotic liver, and compartment 2 (C2) represented a primary hepatocellular carcinoma tumor sample embedded within C1. To evaluate the sensitivity to tissue parameters, a range of clinically relevant tissue…

PURPOSE: To determine how training set size affects the accuracy of knowledge-based treatment planning (KBP) models. METHODS: The authors selected four models from three classes of KBP approaches, corresponding to three distinct quantities that KBP models may predict: dose-volume histogram (DVH) points, DVH curves, and objective function weights. DVH point prediction is done using the best plan from a database of similar clinical plans; DVH curve prediction employs principal component analysis and multiple linear regression; and objective function weights uses either logistic regression or K-nearest neighbors. The authors trained each KBP model using training sets of sizes n = 10, 20, 30, 50, 75, 100, 150, and 200. The authors set aside 100 randomly selected patients from their cohort of 315 prostate cancer patients from Princess Margaret Cancer Center to serve as a validation set for all experiments. For each value of n, the authors randomly selected 100 different training sets with replacement from the remaining 215 patients. Each of the 100 training sets was used to train a model…

PURPOSE: This paper represents the initial phase of a proposed operator-friendly semiautomatic method for positioning and directing an intravascular three-magnet tip catheter in the human heart using an electromagnetic system. METHODS: A predictive computer algorithm based on a comprehensive mathematical model is developed, which accurately calculates the magnetic field generated by the electromagnet system as well as the magnetic torques and forces exerted on a three-magnet tip catheter, and generates the necessary electromagnet currents for arbitrary displacement and deflection of the catheter tip within a workspace of 128 × 128 × 128 mm. RESULTS: We demonstrate the ability of the developed mathematical model to accurately position a three-magnet tip catheter within the 128 × 128 × 128 mm workspace of a 3D eight-electromagnet system. CONCLUSIONS: The ability of the developed mathematical model in predicting the displacement, direction, and deflection of the catheter tip as a function of the electromagnet current values has been verified through experimental results.

PURPOSE: , and mean excitation energy, I, are avoided. METHODS: The SPR parametrization proposed in this study is a purely empirical fit based on the theoretical SPR values for a list of 34 reference human tissues. To investigate the SPR estimation made with this new method the authors performed a calibration and an evaluation with the method. The authors initially calculated CT numbers using CT energy spectrum characterization parameters obtained from calibration based on a Gammex 467 electron density calibration phantom. These CT numbers were fitted to the theoretical SPR for the reference human tissues using the new SPR parametrization presented in this study. The method was evaluated based on theoretical CT numbers for the reference human tissues. The root-mean-square error (RMSE) of the SPR and the proton range error from the continuous slowing down approximation were calculated for the reference human tissues. To test the stability of the parametrization the authors varied the density and elemental composition of the reference human tissues and calculated their new SPR…

PURPOSE: To investigate the effects of dose level and reconstruction method on density and texture based features computed from CT lung nodules. METHODS: This study had two major components. In the first component, a uniform water phantom was scanned at three dose levels and images were reconstructed using four conventional filtered backprojection (FBP) and four iterative reconstruction (IR) methods for a total of 24 different combinations of acquisition and reconstruction conditions. In the second component, raw projection (sinogram) data were obtained for 33 lung nodules from patients scanned as a part of their clinical practice, where low dose acquisitions were simulated by adding noise to sinograms acquired at clinical dose levels (a total of four dose levels) and reconstructed using one FBP kernel and two IR kernels for a total of 12 conditions. For the water phantom, spherical regions of interest (ROIs) were created at multiple locations within the water phantom on one reference image obtained at a reference condition. For the lung nodule cases, the ROI of each nodule was…

PURPOSE: To develop a practical background compensation (BC) technique to improve quantitative (90)Y-bremsstrahlung single-photon emission computed tomography (SPECT)/computed tomography (CT) using a commercially available imaging system. METHODS: All images were acquired using medium-energy collimation in six energy windows (EWs), ranging from 70 to 410 keV. The EWs were determined based on the signal-to-background ratio in planar images of an acrylic phantom of different thicknesses (2-16 cm) positioned below a (90)Y source and set at different distances (15-35 cm) from a gamma camera. The authors adapted the widely used EW-based scatter-correction technique by modeling the BC as scaled images. The BC EW was determined empirically in SPECT/CT studies using an IEC phantom based on the sphere activity recovery and residual activity in the cold lung insert. The scaling factor was calculated from 20 clinical planar (90)Y images. Reconstruction parameters were optimized in the same SPECT images for improved image quantification and contrast. A count-to-activity calibration factor was…

PURPOSE: To present a clinical case in which real-time intratreatment imaging identified an multileaf collimator (MLC) leaf to be consistently deviating from its programmed and logged position by >1 mm. METHODS: An EPID-based exit-fluence dosimetry system designed to prevent gross delivery errors was used to capture cine during treatment images. The author serendipitously visually identified a suspected MLC leaf displacement that was not otherwise detected. The leaf position as recorded on the EPID images was measured and log-files were analyzed for the treatment in question, the prior day's treatment, and for daily MLC test patterns acquired on those treatment days. Additional standard test patterns were used to quantify the leaf position. RESULTS: Whereas the log-file reported no difference between planned and recorded positions, image-based measurements showed the leaf to be 1.3 ± 0.1 mm medial from the planned position. This offset was confirmed with the test pattern irradiations. CONCLUSIONS: It has been clinically observed that log-file derived leaf positions can differ from…

PURPOSE: The gamma comparison is widely used to evaluate the agreement between measurements and treatment planning system calculations in patient-specific intensity modulated radiation therapy (IMRT) quality assurance (QA). However, recent publications have raised concerns about the lack of sensitivity when employing commonly used gamma criteria. Understanding the actual sensitivity of a wide range of different gamma criteria may allow the definition of more meaningful gamma criteria and tolerance limits in IMRT QA. We present a method that allows the quantitative determination of gamma criteria sensitivity to induced errors which can be applied to any unique combination of device, delivery technique, and software utilized in a specific clinic. METHODS: A total of 21 DMLC IMRT QA measurements (ArcCHECK®, Sun Nuclear) were compared to QA plan calculations with induced errors. Three scenarios were studied: MU errors, multi-leaf collimator (MLC) errors, and the sensitivity of the gamma comparison to changes in penumbra width. Gamma comparisons were performed between measurements and…

PURPOSE: To propose a formalism for the reference dosimetry of scanned light-ion beams consistent with IAEA TRS-398 and Alfonso et al. [Med. Phys. 35, 5179-5186 (2008)]. To identify machine-specific reference (msr) fields and plan-class specific reference (pcsr) fields consistent with the definitions given by Alfonso et al. To review the literature of beam monitor calibration in scanned beams using three different methods in terms of this common formalism. METHODS: Four types of msr fields are identified as those that are meant to calibrate the beam monitor for scanned beams with particular energies. Two types of pcsr fields are identified as those that are meant to apply one or more tuning factors to the entire delivery chain. RESULTS: The formalism establishes the energy-dependent relation between the number of particles incident on the phantom surface and the beam monitor reading and distinguishes three routes to determine the beam monitor calibration function: (i) the use of a calibrated reference ionization chamber in a single-layer scanned beam, (ii) the use of a…

Purpose The mechanical and imaging properties of polyvinyl chloride (PVC) can be adjusted to meet the needs of researchers as a tissue‐mimicking material. For instance, the hardness can be adjusted by changing the ratio of softener to PVC polymer, mineral oil can be added for lubrication in needle insertion, and glass beads can be added to scatter acoustic energy similar to biological tissue. Through this research, the authors sought to develop a regression model to design formulations of PVC with targeted mechanical and multimodal medical imaging properties. Methods The design of experiment was conducted by varying three factors—(1) the ratio of softener to PVC polymer, (2) the mass fraction of mineral oil, and (3) the mass fraction of glass beads—and measuring the mechanical properties (elastic modulus, hardness, viscoelastic relaxation time constant, and needle insertion friction force) and the medical imaging properties [speed of sound, acoustic attenuation coefficient, magnetic resonance imaging time constants T 1 and T 2 , and the transmittance of the visible light at…

Purpose: An x‐ray photon interacts with photon counting detectors (PCDs) and generates an electron charge cloud or multiple clouds. The clouds (thus, the photon energy) may be split between two adjacent PCD pixels when the interaction occurs near pixel boundaries, producing a count at both of the pixels. This is called double‐counting with charge sharing. (A photoelectric effect with K‐shell fluorescence x‐ray emission would result in double‐counting as well). As a result, PCD data are spatially and energetically correlated, although the output of individual PCD pixels is Poisson distributed. Major problems include the lack of a detector noise model for the spatio‐energetic cross talk and lack of a computationally efficient simulation tool for generating correlated Poisson data. A Monte Carlo (MC) simulation can accurately simulate these phenomena and produce noisy data; however, it is not computationally efficient. Methods: In this study, the authors developed a new detector model and implemented it in an efficient software simulator that uses a Poisson random number generator to…

PURPOSE: Statistical modeling of an ultrasound backscattered echo envelope is used for tissue characterization. However, in the presence of complex structures within the analyzed area, estimation of parameters is disturbed and unreliable, e.g., in the case of breast tumor classification. In order to improve the differentiation of breast lesions, the authors proposed a method based on the segmentation of homodyned K distribution parameter maps. Regions within lesions of different scattering properties were extracted and analyzed. In order to improve the classification, the best-performing features were selected from various regions and then combined. METHODS: A radio-frequency data set consisting of 103 breast lesions was used in the authors' analysis. Maps of homodyned K distribution parameters were created using an algorithm based on signal-to-noise ratio, kurtosis, and skewness of fractional-order envelope moments. A Markov random field model was used to segment parametric maps. Features of different segments were extracted and evaluated based on bootstrapping and the receiver…

PURPOSE: The pretreatment physics plan review is a standard tool for ensuring treatment quality. Studies have shown that the majority of errors in radiation oncology originate in treatment planning, which underscores the importance of the pretreatment physics plan review. This quality assurance measure is fundamentally important and central to the safety of patients and the quality of care that they receive. However, little is known about its effectiveness. The purpose of this study was to analyze reported incidents to quantify the effectiveness of the pretreatment physics plan review with the goal of improving it. METHODS: This study analyzed 522 potentially severe or critical near-miss events within an institutional incident learning system collected over a three-year period. Of these 522 events, 356 originated at a workflow point that was prior to the pretreatment physics plan review. The remaining 166 events originated after the pretreatment physics plan review and were not considered in the study. The applicable 356 events were classified into one of the three categories: (1)…

PURPOSE: Multiatlas based segmentation is largely used in many clinical and research applications. Due to its good performances, it has recently been included in some commercial platforms for radiotherapy planning and surgery guidance. Anyway, to date, a software with no restrictions about the anatomical district and image modality is still missing. In this paper we introduce plastimatch mabs, an open source software that can be used with any image modality for automatic segmentation. METHODS: plastimatch mabs workflow consists of two main parts: (1) an offline phase, where optimal registration and voting parameters are tuned and (2) an online phase, where a new patient is labeled from scratch by using the same parameters as identified in the former phase. Several registration strategies, as well as different voting criteria can be selected. A flexible atlas selection scheme is also available. To prove the effectiveness of the proposed software across anatomical districts and image modalities, it was tested on two very different scenarios: head and neck (H&N) CT segmentation for…

PURPOSE: The integration of MRI with a linear accelerator (MR-linac) offers great potential for high-precision delivery of radiation therapy (RT). However, the electron deflection resulting from the presence of a transverse magnetic field (TMF) can affect the dose distribution, particularly the electron return effect (ERE) at tissue interfaces. The purpose of the study is to investigate the dose effects of ERE at air-tissue and lung-tissue interfaces during intensity-modulated radiation therapy (IMRT) planning. METHODS: IMRT and volumetric modulated arc therapy (VMAT) plans for representative pancreas, lung, breast, and head and neck (HN) cases were generated following commonly used clinical dose volume (DV) criteria. In each case, three types of plans were generated: (1) the original plan generated without a TMF; (2) the reconstructed plan generated by recalculating the original plan with the presence of a TMF of 1.5 T (no optimization); and (3) the optimized plan generated by a full optimization with TMF = 1.5 T. These plans were compared using a variety of DV parameters,…

PURPOSE: The clinical commissioning of IMRT subject to a magnetic field is challenging. The purpose of this work is to develop a GPU-accelerated Monte Carlo dose calculation platform based on penelope and then use the platform to validate a vendor-provided MRIdian head model toward quality assurance of clinical IMRT treatment plans subject to a 0.35 T magnetic field. METHODS: penelope was first translated from fortran to c++ and the result was confirmed to produce equivalent results to the original code. The c++ code was then adapted to cuda in a workflow optimized for GPU architecture. The original code was expanded to include voxelized transport with Woodcock tracking, faster electron/positron propagation in a magnetic field, and several features that make gpenelope highly user-friendly. Moreover, the vendor-provided MRIdian head model was incorporated into the code in an effort to apply gpenelope as both an accurate and rapid dose validation system. A set of experimental measurements were performed on the MRIdian system to examine the accuracy of both the head model and…

PURPOSE: Phantom surrogates were developed to allow multimodal [computed tomography (CT), magnetic resonance imaging (MRI), and teletherapy] and anthropomorphic tissue simulation as well as materials and methods to construct deformable organ shapes and anthropomorphic bone models. METHODS: Agarose gels of variable concentrations and loadings were investigated to simulate various soft tissue types. Oils, fats, and Vaseline were investigated as surrogates for adipose tissue and bone marrow. Anthropomorphic shapes of bone and organs were realized using 3D-printing techniques based on segmentations of patient CT-scans. All materials were characterized in dual energy CT and MRI to adapt CT numbers, electron density, effective atomic number, as well as T1- and T2-relaxation times to patient and literature values. RESULTS: Soft tissue simulation could be achieved with agarose gels in combination with a gadolinium-based contrast agent and NaF to simulate muscle, prostate, and tumor tissues. Vegetable oils were shown to be a good representation for adipose tissue in all modalities. Inner…

PURPOSE: To develop an automated system for mediastinal lymph node detection and station mapping for chest CT. METHODS: The contextual organs, trachea, lungs, and spine are first automatically identified to locate the region of interest (ROI) (mediastinum). The authors employ shape features derived from Hessian analysis, local object scale, and circular transformation that are computed per voxel in the ROI. Eight more anatomical structures are simultaneously segmented by multiatlas label fusion. Spatial priors are defined as the relative multidimensional distance vectors corresponding to each structure. Intensity, shape, and spatial prior features are integrated and parsed by a random forest classifier for lymph node detection. The detected candidates are then segmented by the following curve evolution process. Texture features are computed on the segmented lymph nodes and a support vector machine committee is used for final classification. For lymph node station labeling, based on the segmentation results of the above anatomical structures, the textual definitions of mediastinal…

PURPOSE: Delivery errors during radiotherapy may lead to medical harm and reduced life expectancy for patients. Such serious incidents can be avoided by performing dose verification online, i.e., while the patient is being irradiated, creating the possibility of halting the linac in case of a large overdosage or underdosage. The offline EPID-based 3D in vivo dosimetry system clinically employed at our institute is in principle suited for online treatment verification, provided the system is able to complete 3D dose reconstruction and verification within 420 ms, the present acquisition time of a single EPID frame. It is the aim of this study to show that our EPID-based dosimetry system can be made fast enough to achieve online 3D in vivo dose verification. METHODS: The current dose verification system was sped up in two ways. First, a new software package was developed to perform all computations that are not dependent on portal image acquisition separately, thus removing the need for doing these calculations in real time. Second, the 3D dose reconstruction algorithm was sped up via…

PURPOSE: Intravascular optical coherence tomography (IV-OCT) is a high-resolution imaging method used to visualize the microstructure of arterial walls in vivo. IV-OCT enables the clinician to clearly observe and accurately measure stent apposition and neointimal coverage of coronary stents, which are associated with side effects such as in-stent thrombosis. In this study, the authors present an algorithm for quantifying stent apposition and neointimal coverage by automatically detecting lumen contours and stent struts in IV-OCT images. METHODS: The algorithm utilizes OCT intensity images and their first and second gradient images along the axial direction to detect lumen contours and stent strut candidates. These stent strut candidates are classified into true and false stent struts based on their features, using an artificial neural network with one hidden layer and ten nodes. After segmentation, either the protrusion distance (PD) or neointimal thickness (NT) for each strut is measured automatically. In randomly selected image sets covering a large variety of clinical scenarios,…

PURPOSE: Cortical thickness and density are critical components in determining the strength of bony structures. Computed tomography (CT) is one possible modality for analyzing the cortex in 3D. In this paper, a model-based approach for measuring the cortical bone thickness and density from clinical CT images is proposed. METHODS: Density variations across the cortex were modeled as a function of the cortical thickness and density, location of the cortex, density of surrounding tissues, and imaging blur. High resolution micro-CT data of cadaver proximal femurs were analyzed to determine a relationship between cortical thickness and density. This thickness-density relationship was used as prior information to be incorporated in the model to obtain accurate measurements of cortical thickness and density from clinical CT volumes. The method was validated using micro-CT scans of 23 cadaver proximal femurs. Simulated clinical CT images with different voxel sizes were generated from the micro-CT data. Cortical thickness and density were estimated from the simulated images using the…

PURPOSE: A method was developed utilizing Cherenkov imaging for rapid and thorough determination of the two gantry angles that produce the most uniform treatment plane during dual-field total skin electron beam therapy (TSET). METHODS: Cherenkov imaging was implemented to gather 2D measurements of relative surface dose from 6 MeV electron beams on a white polyethylene sheet. An intensified charge-coupled device camera time-gated to the Linac was used for Cherenkov emission imaging at sixty-two different gantry angles (1° increments, from 239.5° to 300.5°). Following a modified Stanford TSET technique, which uses two fields per patient position for full body coverage, composite images were created as the sum of two beam images on the sheet; each angle pair was evaluated for minimum variation across the patient region of interest. Cherenkov versus dose correlation was verified with ionization chamber measurements. The process was repeated at source to surface distance (SSD) = 441, 370.5, and 300 cm to determine optimal angle spread for varying room geometries. In addition, three…

PURPOSE: In many clinical MRI scenarios, existing imaging information can be used to significantly shorten acquisition time or to improve Signal to Noise Ratio (SNR). In this paper the authors present a framework, referred to as FASTMER, for fast MRI by exploiting a reference image. METHODS: The proposed approach utilizes the possible similarity of the reference image to the acquired image, which exists in many clinical MRI imaging scenarios. Examples include similarity between adjacent slices in high resolution MRI, similarity between various contrasts in the same scan and similarity between different scans of the same patient. To account for the fact that the reference image may exhibit low similarity with the acquired image the authors develop an iterative weighted reconstruction approach, which tunes the weights according to the degree of similarity. RESULTS: Experimental results demonstrate the performance of the method in three different clinical MRI scenarios: The first example demonstrates SNR improvement in high resolution brain MRI, the second scenario exploits similarity…

PURPOSE: Although current computed tomography (CT) systems can scan the head in a very short time, patient motion sometimes still induces artifacts. If motion occurs, one has to repeat the scan; to avoid motion, sedation or anesthesia is sometimes applied. METHODS: The authors propose a method to iteratively estimate and compensate this motion during the reconstruction. In every iteration, the rigid motion was estimated view-by-view and then used to update the system matrix. A multiresolution scheme was used to speed up the convergence of this joint estimation of the image and the motion of the subject. A final iterative reconstruction was performed with the last motion estimate. RESULTS: The method was evaluated on simulations, patient scans, and a phantom study. The quality of the reconstructed images was improved substantially after the compensation. In simulation and phantom studies, root-mean-square error was reduced and mean structural similarity was increased. In the patient studies, most of the motion blurring in the reconstructed images disappeared after the compensation.…

PURPOSE: The purpose of this study was to evaluate, over a wide range of phantom sizes, CT number stability achieved using two techniques for generating dual-energy computed tomography (DECT) virtual monoenergetic images. METHODS: Water phantoms ranging in lateral diameter from 15 to 50 cm and containing a CT number test object were scanned on a DSCT scanner using both single-energy (SE) and dual-energy (DE) techniques. The SE tube potentials were 70, 80, 90, 100, 110, 120, 130, 140, and 150 kV; the DE tube potential pairs were 80/140, 70/150Sn, 80/150Sn, 90/150Sn, and 100/150Sn kV (Sn denotes that the 150 kV beam was filtered with a 0.6 mm tin filter). Virtual monoenergetic images at energies ranging from 40 to 140 keV were produced from the DECT data using two algorithms, monoenergetic (mono) and monoenergetic plus (mono+). Particularly in large phantoms, water CT number errors and/or artifacts were observed; thus, datasets with water CT numbers outside ±10 HU or with noticeable artifacts were excluded from the study. CT numbers were measured to determine CT number stability…

PURPOSE: The aim of this work was to experimentally demonstrate the feasibility of x-ray acoustic computed tomography (XACT) as a dosimetry tool in a clinical radiotherapy environment. METHODS: The acoustic waves induced following a single pulse of linear accelerator irradiation in a water tank were detected with an immersion ultrasound transducer. By rotating the collimator and keeping the transducer stationary, acoustic signals at varying angles surrounding the field were detected and reconstructed to form an XACT image. Simulated XACT images were obtained using a previously developed simulation workflow. Profiles extracted from experimental and simulated XACT images were compared to profiles measured with an ion chamber. A variety of radiation field sizes and shapes were investigated. RESULTS: XACT images resembling the geometry of the delivered radiation field were obtained for fields ranging from simple squares to more complex shapes. When comparing profiles extracted from simulated and experimental XACT images of a 4 cm × 4 cm field, 97% of points were found to pass a 3%/3 mm…

PURPOSE: Radioembolization with yttrium-90 microspheres may be optimized with patient-specific pretherapy treatment planning. Dose verification and validation of treatment planning methods require quantitative imaging of the post-therapy distribution of yttrium-90 (Y-90). Methods for quantitative imaging of Y-90 using both bremsstrahlung SPECT and PET have previously been described. The purpose of this study was to compare the two modalities quantitatively in humans. METHODS: Calibration correction factors for both quantitative Y-90 bremsstrahlung SPECT and a non-time-of-flight PET system without compensation for prompt coincidences were developed by imaging three phantoms. The consistency of these calibration correction factors for the different phantoms was evaluated. Post-therapy images from both modalities were obtained from 15 patients with hepatocellular carcinoma who underwent hepatic radioembolization using Y-90 glass microspheres. Quantitative SPECT and PET images were rigidly registered and the total liver activities and activity distributions estimated for each modality…

PURPOSE: The highest photon fluence rate that a computed tomography (CT) detector must be able to measure is an important parameter. The authors calculate the maximum transmitted fluence rate in a commercial CT scanner as a function of patient size for standard head, chest, and abdomen protocols. METHODS: The authors scanned an anthropomorphic phantom (Kyoto Kagaku PBU-60) with the reference CT protocols provided by AAPM on a GE LightSpeed VCT scanner and noted the tube current applied with the tube current modulation (TCM) system. By rescaling this tube current using published measurements on the tube current modulation of a GE scanner [N. Keat, "CT scanner automatic exposure control systems," MHRA Evaluation Report 05016, ImPACT, London, UK, 2005], the authors could estimate the tube current that these protocols would have resulted in for other patient sizes. An ECG gated chest protocol was also simulated. Using measured dose rate profiles along the bowtie filters, the authors simulated imaging of anonymized patient images with a range of sizes on a GE VCT scanner and calculated…

PURPOSE: The pursuit of real-time image guided radiotherapy using optimal tissue contrast has seen the development of several hybrid magnetic resonance imaging (MRI)-treatment systems, high field and low field, and inline and perpendicular configurations. As part of a new MRI-linac program, an MRI scanner was integrated with a linear accelerator to enable investigations of a coupled inline MRI-linac system. This work describes results from a prototype experimental system to demonstrate the feasibility of a high field inline MR-linac. METHODS: The magnet is a 1.5 T MRI system (Sonata, Siemens Healthcare) was located in a purpose built radiofrequency (RF) cage enabling shielding from and close proximity to a linear accelerator with inline (and future perpendicular) orientation. A portable linear accelerator (Linatron, Varian) was installed together with a multileaf collimator (Millennium, Varian) to provide dynamic field collimation and the whole assembly built onto a stainless-steel rail system. A series of MRI-linac experiments was performed to investigate (1) image quality with…

PURPOSE: Dual-energy CT (DECT) expands applications of CT imaging in its capability to decompose CT images into material images. However, decomposition via direct matrix inversion leads to large noise amplification and limits quantitative use of DECT. Their group has previously developed a noise suppression algorithm via penalized weighted least-square optimization with edge-preservation regularization (PWLS-EPR). In this paper, the authors improve method performance using the same framework of penalized weighted least-square optimization but with similarity-based regularization (PWLS-SBR), which substantially enhances the quality of decomposed images by retaining a more uniform noise power spectrum (NPS). METHODS: The design of PWLS-SBR is based on the fact that averaging pixels of similar materials gives a low-noise image. For each pixel, the authors calculate the similarity to other pixels in its neighborhood by comparing CT values. Using an empirical Gaussian model, the authors assign high/low similarity value to one neighboring pixel if its CT value is close/far to the CT value…

PURPOSE: Neurovascular interventional procedures using biplane fluoroscopic imaging systems can lead to increased risk of radiation-induced skin injuries. The authors developed a biplane dose tracking system (Biplane-DTS) to calculate the cumulative skin dose distribution from the frontal and lateral x-ray tubes and display it in real-time as a color-coded map on a 3D graphic of the patient for immediate feedback to the physician. The agreement of the calculated values with the dose measured on phantoms was evaluated. METHODS: The Biplane-DTS consists of multiple components including 3D graphic models of the imaging system and patient, an interactive graphical user interface, a data acquisition module to collect geometry and exposure parameters, the computer graphics processing unit, and functions for determining which parts of the patient graphic skin surface are within the beam and for calculating dose. The dose is calculated to individual points on the patient graphic using premeasured calibration files of entrance skin dose per mAs including backscatter; corrections are applied…

PURPOSE: Although a number of studies have focused on finding anatomical regions in which iron concentrations are high, no study has been conducted to examine the overall variations in susceptibility maps of Alzheimer's disease (AD). The objective of this study, therefore, was to differentiate AD from cognitive normal (CN) and mild cognitive impairment (MCI) using a texture analysis of quantitative susceptibility maps (QSMs). METHODS: The study was approved by the local institutional review board, and informed consent was obtained from all subjects. In each participant group-CN, MCI, and AD-18 elderly subjects were enrolled. A fully first-order flow-compensated 3D gradient-echo sequence was run to obtain axial magnitudes and phase images and to produce QSM data. Sagittal structural 3D T1-weighted (3DT1W) images were also obtained with the magnetization-prepared rapid acquisition of gradient-echo sequence to obtain brain tissue images. The first- and second-order texture parameters of the QSMs and 3DT1W images were obtained to evaluate group differences using a one-way analysis of…

PURPOSE: The authors have developed a new 3D breast image reconstruction technique that utilizes the soft tissue spatial resolution of magnetic resonance imaging (MRI) and integrates the dielectric property differentiation from microwave imaging to produce a dual modality approach with the goal of augmenting the specificity of MR imaging, possibly without the need for nonspecific contrast agents. The integration is performed through the application of a soft prior regularization which imports segmented geometric meshes generated from MR exams and uses it to constrain the microwave tomography algorithm to recover nearly uniform property distributions within segmented regions with sharp delineation between these internal subzones. METHODS: Previous investigations have demonstrated that this approach is effective in 2D simulation and phantom experiments and also in clinical exams. The current study extends the algorithm to 3D and provides a thorough analysis of the sensitivity and robustness to misalignment errors in size and location between the spatial prior information and the…

PURPOSE: Dose calculation errors near metal implants are caused by limitations of the dose calculation algorithm in modeling tissue/metal interface effects as well as density assignment errors caused by imaging artifacts. The purpose of this study was to investigate two strategies for reducing dose calculation errors near metal implants: implementation of metal-based energy deposition kernels in the convolution/superposition (C/S) dose calculation method and use of metal artifact reduction methods for computed tomography (CT) imaging. METHODS: Both error reduction strategies were investigated using a simple geometric slab phantom with a rectangular metal insert (composed of titanium or Cerrobend), as well as two anthropomorphic phantoms (one with spinal hardware and one with dental fillings), designed to mimic relevant clinical scenarios. To assess the dosimetric impact of metal kernels, the authors implemented titanium and silver kernels in a commercial collapsed cone C/S algorithm. To assess the impact of CT metal artifact reduction methods, the authors performed dose calculations…

PURPOSE: Compton cameras (CCs) use electronic collimation to reconstruct the images of activity distribution. Although this approach can greatly improve imaging efficiency, due to complex geometry of the CC principle, image reconstruction with the standard iterative algorithms, such as ordered subset expectation maximization (OSEM), can be very time-consuming, even more so if resolution recovery (RR) is implemented. We have previously shown that the origin ensemble (OE) algorithm can be used for the reconstruction of the CC data. Here we propose a method of extending our OE algorithm to include RR. METHODS: To validate the proposed algorithm we used Monte Carlo simulations of a CC composed of multiple layers of pixelated CZT detectors and designed for imaging small animals. A series of CC acquisitions of small hot spheres and the Derenzo phantom placed in air were simulated. Images obtained from (a) the exact data, (b) blurred data but reconstructed without resolution recovery, and (c) blurred and reconstructed with resolution recovery were compared. Furthermore, the reconstructed…

PURPOSE: Biplane angiography systems provide time resolved 2D fluoroscopic images from two different angles, which can be used for the positioning of interventional devices such as guidewires and catheters. The purpose of this work is to provide a novel algorithm framework, which allows the 3D reconstruction of these curvilinear devices from the 2D projection images for each time frame. This would allow creating virtual projection images from arbitrary view angles without changing the position of the gantries, as well as virtual endoscopic 3D renderings. METHODS: The first frame of each time sequence is registered to and subtracted from the following frame using an elastic grid registration technique. The images are then preprocessed by a noise reduction algorithm using directional adaptive filter kernels and a ridgeness filter that emphasizes curvilinear structures. A threshold based segmentation of the device is then performed, followed by a flux driven topology preserving thinning algorithm to extract the segments of the device centerline. The exact device path is determined…

PURPOSE: Correction factors for the relative dosimetry of narrow megavoltage photon beams have recently been determined in several publications. These corrections are required because of the several small-field effects generally thought to be caused by the lack of lateral charged particle equilibrium (LCPE) in narrow beams. Correction factors for relative dosimetry are ultimately necessary to account for the fluence perturbation caused by the detector. For most small field detectors the perturbation depends on field size, resulting in large correction factors when the field size is decreased. In this work, electron and photon fluence differential in energy will be calculated within the radiation sensitive volume of a number of small field detectors for 6 MV linear accelerator beams. The calculated electron spectra will be used to determine electron fluence perturbation as a function of field size and its implication on small field dosimetry analyzed. METHODS: Fluence spectra were calculated with the user code PenEasy, based on the PENELOPE Monte Carlo system. The detectors simulated…

PURPOSE: To introduce the heterogeneous multiscale (HetMS) model for Monte Carlo simulations of gold nanoparticle dose-enhanced radiation therapy (GNPT), a model characterized by its varying levels of detail on different length scales within a single phantom; to apply the HetMS model in two different scenarios relevant for GNPT and to compare computed results with others published. METHODS: The HetMS model is implemented using an extended version of the EGSnrc user-code egs_chamber; the extended code is tested and verified via comparisons with recently published data from independent GNP simulations. Two distinct scenarios for the HetMS model are then considered: (a) monoenergetic photon beams (20 keV to 1 MeV) incident on a cylinder (1 cm radius, 3 cm length); (b) isotropic point source (brachytherapy source spectra) at the center of a 2.5 cm radius sphere with gold nanoparticles (GNPs) diffusing outwards from the center. Dose enhancement factors (DEFs) are compared for different source energies, depths in phantom, gold concentrations, GNP sizes, and modeling assumptions, as well…

PURPOSE: Bronchiectasis and airway wall thickening are commonly assessed in computed tomography (CT) by comparing the airway size with the size of the accompanying artery. Thus, in order to automate the quantification of bronchiectasis and wall thickening following a similar principle, there is a need for methods that automatically segment the airway and vascular trees, measure their size, and pair each airway branch with its accompanying artery. METHODS: This paper combines and extends existing techniques to present a fully automated pipeline that, given a thoracic chest CT, segments, measures, and pairs airway branches with the accompanying artery, then quantifies airway wall thickening and bronchiectasis by measuring the wall-artery ratio (WAR) and lumen and outer wall airway-artery ratio (AAR). Measurements that do not use the artery size for normalization are also extracted, including wall area percentage (WAP), wall thickness ratio (WTR), and airway diameters. RESULTS: The method was thoroughly evaluated using 8000 manual annotations of airway-artery pairs from 24…

Purpose Small field measurements are challenging, due to the physical characteristics coming from the lack of charged particle equilibrium, the partial occlusion of the finite radiation source, and to the detector response. These characteristics can be modeled in the dose calculations in the treatment planning systems. Aim of the present work is to evaluate the MU calculation accuracy for small fields, defined by jaw or MLC, for anisotropic analytical algorithm (AAA) and Acuros XB algorithms, relative to output measurements on the beam central axis. Methods Single point output factor measurement was acquired with a PTW microDiamond detector for 6 MV, 6 and 10 MV unflattened beams generated by a Varian TrueBeam STx equipped with high definition‐MLC. Fields defined by jaw or MLC apertures were set; jaw‐defined: 0.6 × 0.6, 0.8 × 0.8, 1 × 1, 2 × 2, 3 × 3, 4 × 4, 5 × 5, and 10 × 10 cm2; MLC‐defined: 0.5 × 0.5 cm2 to the maximum field defined by the jaw, with 0.5 cm stepping, and jaws set to: 2 × 2, 3 × 3, 4 × 4, 5 × 5, and 10 × 10 cm2. MU calculation was obtained with 1 mm grid in a…