Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Radiological Data”

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

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

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Characterization of a spectroscopic detector for application in x-ray computed tomography

Recent advances in cadmium telluride (CdTe) energy-discriminating pixelated detectors have enabled the possibility of Multi-Spectral X-ray Computed Tomography (MSXCT) to incorporate spectroscopic information into CT. MultiX ME 100 V2 is a CdTe-based spectroscopic x-ray detector array capable of recording energies from 20 to 160 keV in 1.1 keV energy bin increments. Hardware and software have been designed to perform radiographic and computed tomography tasks with this spectroscopic detector. Energy calibration is examined using the end-point energy of a bremsstrahlung spectrum and radioisotope spectral lines. When measuring the spectrum from Am-241 across 500 detector elements, the standard deviation of the peak-location and FWHM measurements are ±0.4 and ±0.6 keV, respectively. As these values are within the energy bin size (1.1 keV), detector elements are consistent with each other. The count rate is characterized, using a nonparalyzable model with a dead time of 64 ± 5 ns. This is consistent with the manufacturer’s quoted per detector-element linear-deviation at 2 Mpps (million photons per sec) of 8.9% (typical) and 12% (max). When comparing measured and simulated spectra, a low-energy tail is visible in the measured data due to the spectral response of the detector. If no valid photon detections are expected in the low-energy tail, then a background subtraction may be applied to allow for a possible first-order correction. If photons are expected in the low-energy tail, a detailed model must be implemented. A radiograph of an aluminum step wedge with a maximum height of about 20 mm shows an underestimation of attenuation by about 10% at 60 keV. This error is due to partial energy deposition from higher-energy (> 60 keV) photons into a lower-energy (~60 keV) bin, reducing the apparent attenuation. A radiograph of a PTFE cylinder taken using a bremsstrahlung spectrum from an x-ray voltage of 100 kV filtered by 1.3 mm Cu is reconstructed using Abel inversion. As no counts are expected in the low energy tail, a first order background correction is applied to the spectrum. The measured linear attenuation coefficient (LAC) is within 10% of the expected value in the 60 to 100 keV range. Below 60 keV, low counts in the corrected spectrum and partial energy deposition from incident photons of energy greater than 60 keV into energy bins below 60 keV impact the LAC measurements. This report ends with a demonstration of the tomographic capability of the system. The quantitative understanding of the detector developed in this report will enable further study in evaluating the system for characterization of an object’s chemical make-up for industrial and security purposes.

36 MATERIALS SCIENCE↗
DOE OSTI2016

SU-F-R-18: Updates to the Computational Environment for Radiological Research for Image Analysis

Purpose: To present new tools in CERR for Texture Analysis and Visualization. Method: (1) Quantitative Image Analysis: We added the ability to compute Haralick texture features based on local neighbourhood. The Texture features depend on many parameters used in their derivation. For example: (a) directionality, (b) quantization of image, (c) patch-size for the neighborhood, (d) handling of the edge voxels within the region of interest, (e) Averaging co-occurance matrix vs texture features for different directions etc. A graphical user interface was built to set these parameters and then visualize their impact on the resulting texture maps. The entire functionality was written in Matlab. Array indexing was used to speed up the texture calculation. The computation speed is very competitive with the ITK library. Moreover, our implementation works with multiple CPUs and the computation time can be further reduced by using multiple processor threads. In order to reduce the Haralick texture maps into scalar features, we propose the use of Texture Volume Histograms. This lets users make use of the entire distribution of texture values within the region of interest rather than using just the mean and the standard deviations. (2) Qualitative/Visualization tools: The derived texture maps are stored as a new scan (derived) within CERR’s planC data structure. A display that compares various scans was built to show the raw image and the derived texture maps side-by-side. These images are positionally linked and can be navigated together. CERR’s graphics handling was updated and sped-up to be compatible with the newer Matlab versions. As a result, the users can use (a) different window levels and colormaps for different viewports, (b) click-and-drag or use mouse scroll-wheel to navigate slices. Results: The new features and updates are available via https://www.github.com/adityaapte/cerr . Conclusion: Features added to CERR increase its utility in Radiomics and Outcomes modeling.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Low Cerebral Blood Volume Identifies Poor Outcome in Stent Retriever Thrombectomy

BackgroundMechanical thrombectomy (MT) is an efficient treatment of acute stroke caused by large-vessel occlusion. We evaluated the factors predicting poor clinical outcome (3-month modified Rankin Scale, mRS >2) although MT performed with modern stent retrievers.MethodsWe prospectively collected the clinical and imaging data of 105 consecutive anterior circulation stroke patients who underwent MT after multimodal CT imaging. Patients with occlusion of the internal carotid artery and/or middle cerebral artery up to the M2 segment were included. We recorded baseline clinical, procedural and imaging variables, technical outcome, 24-h imaging outcome and the clinical outcome. Differences between the groups were studied with appropriate statistical tests and binary logistic regression analysis.ResultsLow cerebral blood volume Alberta stroke program early CT score (CBV-ASPECTS) was associated with poor clinical outcome (median 7 vs. 9, p = 0.01). Lower collateral score (CS) significantly predicted poor outcome in regression modelling with CS = 0 increasing the odds of poor outcome 4.4-fold compared to CS = 3 (95% CI 1.27–15.5, p = 0.02). Lower CBV-ASPECTS significantly predicted poor clinical outcome among those with moderate or severe stroke (OR 0.82, 95% CI 0.68–1, p = 0.05) or poor collateral circulation (CS 0–1, OR 0.66, 95% CI 0.48–0.90, p = 0.009) but not among those with mild strokes or good collaterals.ConclusionsCBV-ASPECTS estimating infarct core is a significant predictor of poor clinical outcome among anterior circulation stroke patients treated with MT, especially in the setting of poor collateral circulation and/or moderate or severe stroke.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2011

Synthetic CT: Simulating low dose single and dual energy protocols from a dual energy scan

Purpose: The choice of CT protocol can greatly impact patient dose and image quality. Since acquiring multiple scans at different techniques on a given patient is undesirable, the ability to predict image quality changes starting from a high quality exam can be quite useful. While existing methods allow one to generate simulated images of lower exposure (mAs) from an acquired CT exam, the authors present and validate a new method called synthetic CT that can generate realistic images of a patient at arbitrary low dose protocols (kVp, mAs, and filtration) for both single and dual energy scans. Methods: The synthetic CT algorithm is derived by carefully ensuring that the expected signal and noise are accurate for the simulated protocol. The method relies on the observation that the material decomposition from a dual energy CT scan allows the transmission of an arbitrary spectrum to be predicted. It requires an initial dual energy scan of the patient to either synthesize raw projections of a single energy scan or synthesize the material decompositions of a dual energy scan. The initial dual energy scan contributes inherent noise to the synthesized projections that must be accounted for before adding more noise to simulate low dose protocols. Therefore, synthetic CT is subject to the constraint that the synthesized data have noise greater than the inherent noise. The authors experimentally validated the synthetic CT algorithm across a range of protocols using a dual energy scan of an acrylic phantom with solutions of different iodine concentrations. An initial 80/140 kVp dual energy scan of the phantom provided the material decomposition necessary to synthesize images at 100 kVp and at 120 kVp, across a range of mAs values. They compared these synthesized single energy scans of the phantom to actual scans at the same protocols. Furthermore, material decompositions of a 100/120 kVp dual energy scan are synthesized by adding correlated noise to the initial material decompositions. The aforementioned noise constraint also allows us to compute feasible mAs values that can be synthesized for each kVp. Results: The single energy synthesized and actual reconstructed images exhibit identical signal and noise properties at 100 kVp and at 120 kVp, and across a range of mAs values. For example, the noise in both the synthesized and actual images at 100 kVp increases by {radical}(2) when the mAs is halved. The synthesized and actual material decompositions of a dual energy protocol show excellent agreement when the decomposition images are linearly weighted to form monoenergetic images at energies from 40 to 100 keV. For simulated single energy protocols with kVp between 80 and 140, the highest feasible mAs exceeds that of either initial scan. Conclusions: This work describes and validates the synthetic CT theory and algorithm by comparing its results to actual scans. Synthetic CT is a powerful new tool that allows users to realistically see how protocol selection affects CT images and enables radiologists to retrospectively identify the lowest dose protocol achievable that provides diagnostic quality images on real patients.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2012

Photon energy-modulated radiotherapy: Monte Carlo simulation and treatment planning study

Purpose: To demonstrate the feasibility of photon energy-modulated radiotherapy during beam-on time. Methods: A cylindrical device made of aluminum was conceptually proposed as an energy modulator. The frame of the device was connected with 20 tubes through which mercury could be injected or drained to adjust the thickness of mercury along the beam axis. In Monte Carlo (MC) simulations, a flattening filter of 6 or 10 MV linac was replaced with the device. The thickness of mercury inside the device varied from 0 to 40 mm at the field sizes of 5 x 5 cm{sup 2} (FS5), 10 x 10 cm{sup 2} (FS10), and 20 x 20 cm{sup 2} (FS20). At least 5 billion histories were followed for each simulation to create phase space files at 100 cm source to surface distance (SSD). In-water beam data were acquired by additional MC simulations using the above phase space files. A treatment planning system (TPS) was commissioned to generate a virtual machine using the MC-generated beam data. Intensity modulated radiation therapy (IMRT) plans for six clinical cases were generated using conventional 6 MV, 6 MV flattening filter free, and energy-modulated photon beams of the virtual machine. Results: As increasing the thickness of mercury, Percentage depth doses (PDD) of modulated 6 and 10 MV after the depth of dose maximum were continuously increased. The amount of PDD increase at the depth of 10 and 20 cm for modulated 6 MV was 4.8% and 5.2% at FS5, 3.9% and 5.0% at FS10 and 3.2%-4.9% at FS20 as increasing the thickness of mercury from 0 to 20 mm. The same for modulated 10 MV was 4.5% and 5.0% at FS5, 3.8% and 4.7% at FS10 and 4.1% and 4.8% at FS20 as increasing the thickness of mercury from 0 to 25 mm. The outputs of modulated 6 MV with 20 mm mercury and of modulated 10 MV with 25 mm mercury were reduced into 30%, and 56% of conventional linac, respectively. The energy-modulated IMRT plans had less integral doses than 6 MV IMRT or 6 MV flattening filter free plans for tumors located in the periphery while maintaining the similar quality of target coverage, homogeneity, and conformity. Conclusions: The MC study for the designed energy modulator demonstrated the feasibility of energy-modulated photon beams available during beam-on time. The planning study showed an advantage of energy-and intensity modulated radiotherapy in terms of integral dose without sacrificing any quality of IMRT plan.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2011

A Bayesian approach to real-time 3D tumor localization via monoscopic x-ray imaging during treatment delivery

Purpose: Monoscopic x-ray imaging with on-board kV devices is an attractive approach for real-time image guidance in modern radiation therapy such as VMAT or IMRT, but it falls short in providing reliable information along the direction of imaging x-ray. By effectively taking consideration of projection data at prior times and/or angles through a Bayesian formalism, the authors develop an algorithm for real-time and full 3D tumor localization with a single x-ray imager during treatment delivery. Methods: First, a prior probability density function is constructed using the 2D tumor locations on the projection images acquired during patient setup. Whenever an x-ray image is acquired during the treatment delivery, the corresponding 2D tumor location on the imager is used to update the likelihood function. The unresolved third dimension is obtained by maximizing the posterior probability distribution. The algorithm can also be used in a retrospective fashion when all the projection images during the treatment delivery are used for 3D localization purposes. The algorithm does not involve complex optimization of any model parameter and therefore can be used in a ''plug-and-play'' fashion. The authors validated the algorithm using (1) simulated 3D linear and elliptic motion and (2) 3D tumor motion trajectories of a lung and a pancreas patient reproduced by a physical phantom. Continuous kV images were acquired over a full gantry rotation with the Varian TrueBeam on-board imaging system. Three scenarios were considered: fluoroscopic setup, cone beam CT setup, and retrospective analysis. Results: For the simulation study, the RMS 3D localization error is 1.2 and 2.4 mm for the linear and elliptic motions, respectively. For the phantom experiments, the 3D localization error is < 1 mm on average and < 1.5 mm at 95th percentile in the lung and pancreas cases for all three scenarios. The difference in 3D localization error for different scenarios is small and is not statistically significant. Conclusions: The proposed algorithm eliminates the need for any population based model parameters in monoscopic image guided radiotherapy and allows accurate and real-time 3D tumor localization on current standard LINACs with a single x-ray imager.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

A limited-angle intrafraction verification (LIVE) system for radiation therapy

Purpose: Currently, no 3D or 4D volumetric x-ray imaging techniques are available for intrafraction verification of target position during actual treatment delivery or in-between treatment beams, which is critical for stereotactic radiosurgery (SRS) and stereotactic body radiation therapy (SBRT) treatments. This study aims to develop a limited-angle intrafraction verification (LIVE) system to use prior information, deformation models, and limited angle kV-MV projections to verify target position intrafractionally. Methods: The LIVE system acquires limited-angle kV projections simultaneously during arc treatment delivery or in-between static 3D/IMRT treatment beams as the gantry moves from one beam to the next. Orthogonal limited-angle MV projections are acquired from the beam's eye view (BEV) exit fluence of arc treatment beam or in-between static beams to provide additional anatomical information. MV projections are converted to kV projections using a linear conversion function. Patient prior planning CT at one phase is used as the prior information, and the on-board patient volume is considered as a deformation of the prior images. The deformation field is solved using the data fidelity constraint, a breathing motion model extracted from the planning 4D-CT based on principal component analysis (PCA) and a free-form deformation (FD) model. LIVE was evaluated using a 4D digital extended cardiac torso phantom (XCAT) and a CIRS 008A dynamic thoracic phantom. In the XCAT study, patient breathing pattern and tumor size changes were simulated from CT to treatment position. In the CIRS phantom study, the artificial target in the lung region experienced both size change and position shift from CT to treatment position. Varian Truebeam research mode was used to acquire kV and MV projections simultaneously during the delivery of a dynamic conformal arc plan. The reconstruction accuracy was evaluated by calculating the 3D volume percentage difference (VPD) and the center of mass (COM) difference of the tumor in the true on-board images and reconstructed images. Results: In both simulation and phantom studies, LIVE achieved substantially better reconstruction accuracy than reconstruction using PCA or FD deformation model alone. In the XCAT study, the average VPD and COM differences among different patient scenarios for LIVE system using orthogonal 30° scan angles were 4.3% and 0.3 mm when using kV+BEV MV. Reducing scan angle to 15° increased the average VPD and COM differences to 15.1% and 1.7 mm. In the CIRS phantom study, the VPD and COM differences for the LIVE system using orthogonal 30° scan angles were 6.4% and 1.4 mm. Reducing scan angle to 15° increased the VPD and COM differences to 51.9% and 3.8 mm. Conclusions: The LIVE system has the potential to substantially improve intrafraction target localization accuracy by providing volumetric verification of tumor position simultaneously during arc treatment delivery or in-between static treatment beams. With this improvement, LIVE opens up a new avenue for margin reduction and dose escalation in both fractionated treatments and SRS and SBRT treatments.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

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

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

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

A Biodistribution and Toxicity Study of Cobalt Dichloride-N-Acetyl Cysteine in an Implantable MRI Marker for Prostate Cancer Treatment

Purpose: C4, a cobalt dichloride-N-acetyl cysteine complex, is being developed as a positive-signal magnetic resonance imaging (MRI) marker to localize implanted radioactive seeds in prostate brachytherapy. We evaluated the toxicity and biodistribution of C4 in rats with the goal of simulating the systemic effects of potential leakage from C4 MRI markers within the prostate. Methods and Materials: 9-μL doses (equivalent to leakage from 120 markers in a human) of control solution (0.9% sodium chloride), 1% (proposed for clinical use), and 10% C4 solution were injected into the prostates of male Sprague-Dawley rats via laparotomy. Organ toxicity and cobalt disposition in plasma, tissues, feces, and urine were evaluated. Results: No C4-related morbidity or mortality was observed in the biodistribution arm (60 rats). Biodistribution was measurable after 10% C4 injection: cobalt was cleared rapidly from periprostatic tissue; mean concentrations in prostate were 163 μg/g and 268 μg/g at 5 and 30 minutes but were undetectable by 60 minutes. Expected dual renal-hepatic elimination was observed, with percentages of injected dose recovered in tissues of 39.0 ± 5.6% (liver), >11.8 ± 6.5% (prostate), and >5.3 ± 0.9% (kidney), with low plasma concentrations detected up to 1 hour (1.40 μg/mL at 5-60 minutes). Excretion in urine was 13.1 ± 4.6%, with 3.1 ± 0.54% recovered in feces by 24 hours. In the toxicity arm, 3 animals died in the control group and 1 each in the 1% and 10% groups from surgical or anesthesia-related complications; all others survived to scheduled termination at 14 days. No C4-related adverse clinical signs or organ toxicity were observed. Conclusion: C4-related toxicity was not observed at exposures at least 10-fold the exposure proposed for use in humans. These data demonstrating lack of systemic toxicity with dual routes of elimination in the event of in situ rupture suggest that C4 warrants further investigation as an MRI marker for prostate brachytherapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS1999

Three-Dimensional, Transgenic Cell Models to Quantify Space Genotoxic Effects

The space environment contains radiation and chemical agents known to be mutagenic and carcinogenic to humans. Additionally, microgravity is a complicating factor that may modify or synergize induced genotoxic effects. Most in vitro models fail to use human cells (making risk extrapolation to humans more difficult), overlook the dynamic effect of tissue intercellular interactions on genotoxic damage, and lack the sensitivity required to measure low-dose effects. Currently a need exists for a model test system that simulates cellular interactions present in tissue, and can be used to quantify genotoxic damage induced by low levels of radiation and chemicals, and extrapolate assessed risk to humans. A state-of-the-art, three-dimensional, multicellular tissue equivalent cell culture model will be presented. It consists of mammalian cells genetically engineered to contain multiple copies of defined target genes for genotoxic assessment,. NASA-designed bioreactors were used to coculture mammalian cells into spheroids, The cells used were human mammary epithelial cells (H184135) and Stratagene's (Austin, Texas) Big Blue(TM) Rat 2 lambda fibroblasts. The fibroblasts were genetically engineered to contain -a high-density target gene for mutagenesis (60 copies of lacl/LacZ per cell). Tissue equivalent spheroids were routinely produced by inoculation of 2 to 7 X 10(exp 5) fibroblasts with Cytodex 3 beads (150 micrometers in diameter). at a 20:1 cell:bead ratio, into 50-ml HARV bioreactors (Synthecon, Inc.). Fibroblasts were cultured for 5 days, an equivalent number of epithelial cells added, and the fibroblast/epithelial cell coculture continued for 21 days. Three-dimensional spheroids with diameters ranging from 400 to 600 micrometers were obtained. Histological and immunohistochemical Characterization revealed i) both cell types present in the spheroids, with fibroblasts located primarily in the center, surrounded by epithelial cells; ii) synthesis of extracellular matrix; and iii,, mitotic cells located throughout the spheroids. Spheroidal integrity and cell viability were retained for the 30-day test period after removal of spheroids from the bioreactor. Potential utility of this three-dimensional, transgenic model for genotoxicity was initially assessed by exposure of spheroids to 0-2 Gy neon at dose rates of 0.3 to 1.5 Gy/min (National Institute of Radiological Sciences, Chiba, Japan). Quantification of mutation at the lacl gene revealed a linear dose response for mutation induction. Limited sequencing analysis of mutant clones revealed higher frequencies of deletions and multiple base sequence changes with increasing dose. These results suggest that our three-dimensional, transgenic model is applicable to a wide variety of studies involving the quantification, identification, and characterization of genotoxicity incurred in space and on Earth. This model uniquely allows investigation of the interaction of relevant factors, namely cell-to-cell interactions and the mechanistic interaction of microgravity with radiation insults and DNA repair. Using this three-dimensional model will allow us to obtain dual genotoxic information (i.e., mutation rate plus chromosome aberration data) from the same system so that one endpoint can be used to reference the other, thereby increasing the fidelity of the data set. Moreover, the tissue-equivalent nature of the three-dimensional model provides high confidence for relevance of risk assessment, i.e., the establishment of quality factors directly applicable to the microgravity environment.

Gonda, S. R.↗
DOE OSTI2014

SU-E-J-67: Evaluation of Breathing Patterns for Respiratory-Gated Radiation Therapy Using Respiration Regularity Index

Purpose: Despite the importance of accurately estimating the respiration regularity of a patient in motion compensation treatment, an effective and simply applicable method has rarely been reported. The authors propose a simple respiration regularity index based on parameters derived from a correspondingly simplified respiration model. Methods: In order to simplify a patient's breathing pattern while preserving the data's intrinsic properties, we defined a respiration model as a power of cosine form with a baseline drift. According to this respiration formula, breathing-pattern fluctuation could be explained using four factors: sample standard deviation of respiration period, sample standard deviation of amplitude and the results of simple regression of the baseline drift (slope and standard deviation of residuals of a respiration signal. Overall irregularity (δ) was defined as a Euclidean norm of newly derived variable using principal component analysis (PCA) for the four fluctuation parameters. Finally, the proposed respiration regularity index was defined as ρ=ln(1+(1/ δ))/2, a higher ρ indicating a more regular breathing pattern. Subsequently, we applied it to simulated and clinical respiration signals from real-time position management (RPM; Varian Medical Systems, Palo Alto, CA) and investigated respiration regularity. Moreover, correlations between the regularity of the first session and the remaining fractions were investigated using Pearson's correlation coefficient. Results: The respiration regularity was determined based on ρ; patients with ρ<0.3 showed worse regularity than the others, whereas ρ>0.7 was suitable for respiratory-gated radiation therapy (RGRT). Fluctuations in breathing cycle and amplitude were especially determinative of ρ. If the respiration regularity of a patient's first session was known, it could be estimated through subsequent sessions. Conclusions: Respiration regularity could be objectively determined using a respiration regularity index, ρ. Such single-index testing of respiration regularity can facilitate determination of RGRT availability in clinical settings, especially for free-breathing cases. This work was supported by a Korea Science and Engineering Foundation (KOSEF) grant funded by the Korean Ministry of Science, ICT and Future Planning (No. 2013043498)

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

A Comparison of Lumpectomy Cavity Delineations Between Use of Magnetic Resonance Imaging and Computed Tomography Acquired With Patient in Prone Position for Radiation Therapy Planning of Breast Cancer

Purpose: To compare lumpectomy cavity (LC) and planning target volume (PTV) delineated with the use of magnetic resonance imaging (MRI) and computed tomography (CT) and to examine the possibility of replacing CT with MRI for radiation therapy (RT) planning for breast cancer. Methods and Materials: MRI and CT data were acquired for 15 patients with early-stage breast cancer undergoing lumpectomy during RT simulation in prone positions, the same as their RT treatment positions. The LCs were delineated manually on both CT (LC-CT) and MRI acquired with 4 sequences: T1, T2, STIR, and DCE. Various PTVs were created by expanding a 15-mm margin from the corresponding LCs and from the union of the LCs for the 4 MRI sequences (PTV-MRI). Differences were measured in terms of cavity visualization score (CVS) and dice coefficient (DC). Results: The mean CVSs for T1, T2, STIR, DCE, and CT defined LCs were 3.47, 3.47, 3.87, 3.50. and 2.60, respectively, implying that the LC is mostly visible with a STIR sequence. The mean reductions of LCs from those for CT were 22%, 43%, 36%, and 17% for T1, T2, STIR, and DCE, respectively. In 14 of 15 cases, MRI (union of T1, T2, STIR, and DCE) defined LC included extra regions that would not be visible from CT. The DCs between CT and MRI (union of T1, T2, STIR, and DCE) defined volumes were 0.65 ± 0.20 for LCs and 0.85 ± 0.06 for PTVs. There was no obvious difference between the volumes of PTV-MRI and PTV-CT, and the average PTV-STIR/PTV-CT volume ratio was 0.83 ± 0.23. Conclusions: The use of MRI improves the visibility of LC in comparison with CT. The volumes of LC and PTV generated based on a MRI sequence are substantially smaller than those based on CT, and the PTV-MRI volumes, defined by the union of T1, T2, STIR, and DCE, were comparable with those of PTV-CT for most of the cases studied.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

Dosimetric effect of intrafraction tumor motion in phase gated lung stereotactic body radiotherapy

Purpose: A major concern for lung intensity modulated radiation therapy delivery is the deviation of actually delivered dose distribution from the planned one due to simultaneous movements of multileaf collimator (MLC) leaves and tumor. For gated lung stereotactic body radiotherapy treatment (SBRT), the situation becomes even more complicated because of SBRT's characteristics such as fewer fractions, smaller target volume, higher dose rate, and extended fractional treatment time. The purpose of this work is to investigate the dosimetric effect of intrafraction tumor motion during gated lung SBRT delivery by reconstructing the delivered dose distribution with real-time tumor motion considered. Methods: The tumor motion data were retrieved from six lung patients. Each of them received three fractions of stereotactic radiotherapy treatments with Cyberknife Synchrony (Accuray, Sunnyvale, CA). Phase gating through an external surrogate was simulated with a gating window of 5 mm. The resulting residual tumor motion curves during gating (beam-on) were retrieved. Planning target volume (PTV) was defined as physician-contoured clinical target volume (CTV) surrounded by an isotropic 5 mm margin. Each patient was prescribed with 60 Gy/3 fractions. The authors developed an algorithm to reconstruct the delivered dose with tumor motion. The DMLC segments, mainly leaf position and segment weighting factor, were recalculated according to the probability density function of tumor motion curve. The new DMLC sequence file was imported back to treatment planning system to reconstruct the dose distribution. Results: Half of the patients in the study group experienced PTV D95% deviation up to 26% for fractional dose and 14% for total dose. CTV mean dose dropped by 1% with tumor motion. Although CTV is almost covered by prescribed dose with 5 mm margin, qualitative comparison on the dose distributions reveals that CTV is on the verge of underdose. The discrepancy happens due to tumor excursion outside of the gating window, which, for our study group, is mainly caused by baseline shift, i.e., the change in general trend of the motion curve during extended period of treatment time. Conclusions: The dose deviation in PTV and CTV due to target motion is not always negligible in gated SBRT. Although CTVs are covered sufficiently with prescribed dose in most cases, some are on the verge of underdose due to large tumor excursion caused by factors such as baseline shift.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Limited-angle effect compensation for respiratory binned cardiac SPECT

Purpose: In cardiac single photon emission computed tomography (SPECT), respiratory-binned study is used to combat the motion blur associated with respiratory motion. However, owing to the variability in respiratory patterns during data acquisition, the acquired data counts can vary significantly both among respiratory bins and among projection angles within individual bins. If not properly accounted for, such variation could lead to artifacts similar to limited-angle effect in image reconstruction. In this work, the authors aim to investigate several reconstruction strategies for compensating the limited-angle effect in respiratory binned data for the purpose of reducing the image artifacts. Methods: The authors first consider a model based correction approach, in which the variation in acquisition time is directly incorporated into the imaging model, such that the data statistics are accurately described among both the projection angles and respiratory bins. Afterward, the authors consider an approximation approach, in which the acquired data are rescaled to accommodate the variation in acquisition time among different projection angles while the imaging model is kept unchanged. In addition, the authors also consider the use of a smoothing prior in reconstruction for suppressing the artifacts associated with limited-angle effect. In our evaluation study, the authors first used Monte Carlo simulated imaging with 4D NCAT phantom wherein the ground truth is known for quantitative comparison. The authors evaluated the accuracy of the reconstructed myocardium using a number of metrics, including regional and overall accuracy of the myocardium, uniformity and spatial resolution of the left ventricle (LV) wall, and detectability of perfusion defect using a channelized Hotelling observer. As a preliminary demonstration, the authors also tested the different approaches on five sets of clinical acquisitions. Results: The quantitative evaluation results show that the three compensation methods could all, but to different extents, reduce the reconstruction artifacts over no compensation. In particular, the model based approach reduced the mean-squared-error of the reconstructed myocardium by as much as 40%. Compared to the approach of data rescaling, the model based approach further improved both the overall and regional accuracy of the myocardium; it also further improved the lesion detectability and the uniformity of the LV wall. When ML reconstruction was used, the model based approach was notably more effective for improving the LV wall; when MAP reconstruction was used, the smoothing prior could reduce the noise level and artifacts with little or no increase in bias, but at the cost of a slight resolution loss of the LV wall. The improvements in image quality by the different compensation methods were also observed in the clinical acquisitions. Conclusions: Compensating for the uneven distribution of acquisition time among both projection angles and respiratory bins can effectively reduce the limited-angle artifacts in respiratory-binned cardiac SPECT reconstruction. Direct incorporation of the time variation into the imaging model together with a smoothing prior in reconstruction can lead to the most improvement in the accuracy of the reconstructed myocardium.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Quantitative modeling of the accuracy in registering preoperative patient-specific anatomic models into left atrial cardiac ablation procedures

Purpose: In cardiac ablation therapy, accurate anatomic guidance is necessary to create effective tissue lesions for elimination of left atrial fibrillation. While fluoroscopy, ultrasound, and electroanatomic maps are important guidance tools, they lack information regarding detailed patient anatomy which can be obtained from high resolution imaging techniques. For this reason, there has been significant effort in incorporating detailed, patient-specific models generated from preoperative imaging datasets into the procedure. Both clinical and animal studies have investigated registration and targeting accuracy when using preoperative models; however, the effect of various error sources on registration accuracy has not been quantitatively evaluated. Methods: Data from phantom, canine, and patient studies are used to model and evaluate registration accuracy. In the phantom studies, data are collected using a magnetically tracked catheter on a static phantom model. Monte Carlo simulation studies were run to evaluate both baseline errors as well as the effect of different sources of error that would be present in a dynamicin vivo setting. Error is simulated by varying the variance parameters on the landmark fiducial, physical target, and surface point locations in the phantom simulation studies. In vivo validation studies were undertaken in six canines in which metal clips were placed in the left atrium to serve as ground truth points. A small clinical evaluation was completed in three patients. Landmark-based and combined landmark and surface-based registration algorithms were evaluated in all studies. In the phantom and canine studies, both target registration error and point-to-surface error are used to assess accuracy. In the patient studies, no ground truth is available and registration accuracy is quantified using point-to-surface error only. Results: The phantom simulation studies demonstrated that combined landmark and surface-based registration improved landmark-only registration provided the noise in the surface points is not excessively high. Increased variability on the landmark fiducials resulted in increased registration errors; however, refinement of the initial landmark registration by the surface-based algorithm can compensate for small initial misalignments. The surface-based registration algorithm is quite robust to noise on the surface points and continues to improve landmark registration even at high levels of noise on the surface points. Both the canine and patient studies also demonstrate that combined landmark and surface registration has lower errors than landmark registration alone. Conclusions: In this work, we describe a model for evaluating the impact of noise variability on the input parameters of a registration algorithm in the context of cardiac ablation therapy. The model can be used to predict both registration error as well as assess which inputs have the largest effect on registration accuracy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

Model-based dose calculations for {sup 125}I lung brachytherapy

Purpose: Model-baseddose calculations (MBDCs) are performed using patient computed tomography (CT) data for patients treated with intraoperative {sup 125}I lung brachytherapy at the Mayo Clinic Rochester. Various metallic artifact correction and tissue assignment schemes are considered and their effects on dose distributions are studied. Dose distributions are compared to those calculated under TG-43 assumptions. Methods: Dose distributions for six patients are calculated using phantoms derived from patient CT data and the EGSnrc user-code BrachyDose. {sup 125}I (GE Healthcare/Oncura model 6711) seeds are fully modeled. Four metallic artifact correction schemes are applied to the CT data phantoms: (1) no correction, (2) a filtered back-projection on a modified virtual sinogram, (3) the reassignment of CT numbers above a threshold in the vicinity of the seeds, and (4) a combination of (2) and (3). Tissue assignment is based on voxel CT number and mass density is assigned using a CT number to mass density calibration. Three tissue assignment schemes with varying levels of detail (20, 11, and 5 tissues) are applied to metallic artifact corrected phantoms. Simulations are also performed under TG-43 assumptions, i.e., seeds in homogeneous water with no interseed attenuation. Results: Significant dose differences (up to 40% for D{sub 90}) are observed between uncorrected and metallic artifact corrected phantoms. For phantoms created with metallic artifact correction schemes (3) and (4), dose volume metrics are generally in good agreement (less than 2% differences for all patients) although there are significant local dose differences. The application of the three tissue assignment schemes results in differences of up to 8% for D{sub 90}; these differences vary between patients. Significant dose differences are seen between fully modeled and TG-43 calculations with TG-43 underestimating the dose (up to 36% in D{sub 90}) for larger volumes containing higher proportions of healthy lung tissue. Conclusions: Metallic artifact correction is necessary for accurate application of MBDCs for lung brachytherapy; simpler threshold replacement methods may be sufficient for early adopters concerned with clinical dose metrics. Rigorous determination of voxel tissue parameters and tissue assignment is required for accurate dose calculations as different tissue assignment schemes can result in significantly different dose distributions. Significant differences are seen between MBDCs and TG-43 dose distributions with TG-43 underestimating dose in volumes containing healthy lung tissue.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

Comparison of TG-43 and TG-186 in breast irradiation using a low energy electronic brachytherapy source

Purpose: The recently updated guidelines for dosimetry in brachytherapy in TG-186 have recommended the use of model-based dosimetry calculations as a replacement for TG-43. TG-186 highlights shortcomings in the water-based approach in TG-43, particularly for low energy brachytherapy sources. The Xoft Axxent is a low energy (<50 kV) brachytherapy system used in accelerated partial breast irradiation (APBI). Breast tissue is a heterogeneous tissue in terms of density and composition. Dosimetric calculations of seven APBI patients treated with Axxent were made using a model-based Monte Carlo platform for a number of tissue models and dose reporting methods and compared to TG-43 based plans. Methods: A model of the Axxent source, the S700, was created and validated against experimental data. CT scans of the patients were used to create realistic multi-tissue/heterogeneous models with breast tissue segmented using a published technique. Alternative water models were used to isolate the influence of tissue heterogeneity and backscatter on the dose distribution. Dose calculations were performed using Geant4 according to the original treatment parameters. The effect of the Axxent balloon applicator used in APBI which could not be modeled in the CT-based model, was modeled using a novel technique that utilizes CAD-based geometries. These techniques were validated experimentally. Results were calculated using two dose reporting methods, dose to water (D{sub w,m}) and dose to medium (D{sub m,m}), for the heterogeneous simulations. All results were compared against TG-43-based dose distributions and evaluated using dose ratio maps and DVH metrics. Changes in skin and PTV dose were highlighted. Results: All simulated heterogeneous models showed a reduced dose to the DVH metrics that is dependent on the method of dose reporting and patient geometry. Based on a prescription dose of 34 Gy, the average D{sub 90} to PTV was reduced by between ∼4% and ∼40%, depending on the scoring method, compared to the TG-43 result. Peak skin dose is also reduced by 10%–15% due to the absence of backscatter not accounted for in TG-43. The balloon applicator also contributed to the reduced dose. Other ROIs showed a difference depending on the method of dose reporting. Conclusions: TG-186-based calculations produce results that are different from TG-43 for the Axxent source. The differences depend strongly on the method of dose reporting. This study highlights the importance of backscatter to peak skin dose. Tissue heterogeneities, applicator, and patient geometries demonstrate the need for a more robust dose calculation method for low energy brachytherapy sources.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

Combined Recipe for Clinical Target Volume and Planning Target Volume Margins

Purpose: To develop a combined recipe for clinical target volume (CTV) and planning target volume (PTV) margins. Methods and Materials: A widely accepted PTV margin recipe is M{sub geo} = aΣ{sub geo} + bσ{sub geo}, with Σ{sub geo} and σ{sub geo} standard deviations (SDs) representing systematic and random geometric uncertainties, respectively. On the basis of histopathology data of breast and lung tumors, we suggest describing the distribution of microscopic islets around the gross tumor volume (GTV) by a half-Gaussian with SD Σ{sub micro}, yielding as possible CTV margin recipe: M{sub micro} = ƒ(N{sub i}) × Σ{sub micro}, with N{sub i} the average number of microscopic islets per patient. To determine ƒ(N{sub i}), a computer model was developed that simulated radiation therapy of a spherical GTV with isotropic distribution of microscopic disease in a large group of virtual patients. The minimal margin that yielded D{sub min} <95% in maximally 10% of patients was calculated for various Σ{sub micro} and N{sub i}. Because Σ{sub micro} is independent of Σ{sub geo}, we propose they should be added quadratically, yielding for a combined GTV-to-PTV margin recipe: M{sub GTV-PTV} = √([aΣ{sub geo}]{sup 2} + [ƒ(N{sub i})Σ{sub micro}]{sup 2}) + bσ{sub geo}. This was validated by the computer model through numerous simultaneous simulations of microscopic and geometric uncertainties. Results: The margin factor ƒ(N{sub i}) in a relevant range of Σ{sub micro} and N{sub i} can be given by: ƒ(N{sub i}) = 1.4 + 0.8log(N{sub i}). Filling in the other factors found in our simulations (a = 2.1 and b = 0.8) yields for the combined recipe: M{sub GTV-PTV} = √((2.1Σ{sub geo}){sup 2} + ([1.4 + 0.8log(N{sub i})] × Σ{sub micro}){sup 2}) + 0.8σ{sub geo}. The average margin difference between the simultaneous simulations and the above recipe was 0.2 ± 0.8 mm (1 SD). Calculating M{sub geo} and M{sub micro} separately and adding them linearly overestimated PTVs by on average 5 mm. Margin recipes based on tumor control probability (TCP) instead of D{sub min} criteria yielded similar results. Conclusions: A general recipe for GTV-to-PTV margins is proposed, which shows that CTV and PTV margins should be added in quadrature instead of linearly.

62 RADIOLOGY AND NUCLEAR MEDICINE↗