Engineering topics

Lee, C

Publications and source records attributed to Lee, C.

At least 19 records

DOE OSTI2025

A comparative analysis of numerical approaches for the description of gas flow in clay-based repository systems: From a laboratory to a large-scale gas injection test

There is nowadays a consensus among many countries that geological disposal is a favourable solution for the long-term management. Although different host formations and different barrier systems are under consideration around the world, clay-based materials form an important component for waste isolation in most national programmes. Hence, a good comprehension of the effect of gas flow on the hydro-mechanical behaviour of clay-based soils is essential, both at laboratory and field scale. Task B under the international cooperative project DECOVALEX-2023 has recently shown that, after some enhancement, models can be employed to reproduce laboratory scale tests, even with different sample geometries 37 . However, further work is required to understand whether they can be applied to simulate a large-scale experiment. Up-scaling of models for the advective transport of gas through clay-based low permeable material presents a number of problems related to the difficulty in obtaining consistent hydrogeological parameters and constitutive relationships at both laboratory and field scale. Based on a unique dataset from a large-scale gas injection test (Lasgit) performed at the Äspö Hard Rock Laboratory (Sweden), Task B within DECOVALEX-2023 has explored the refinement of these numerical strategies applied to the simulation of gas flow. Work performed within the task reveals that codes do not need to be substantially modified from the laboratory models to reproduce full-scale tests: indeed, model parameters calibrated and validated at laboratory scale have been applied to predict field scale gas flow at Lasgit, including peak gas pressure and injected cumulative gas volume. By means of (1) the introduction of interfaces between blocks to reflect the experimental configuration and the (2) adjustment of some parameters (e.g., higher permeability), the updated models are able to represent most of the key features observed in the experimental data, even at a large scale.

Tamayo-Mas, E↗

DOE OSTI2016

SU-F-J-174: A Series of Computational Human Phantoms in DICOM-RT Format for Normal Tissue Dose Reconstruction in Epidemiological Studies

Purpose: Epidemiological studies of second cancer risk in radiotherapy patients often require individualized dose estimates of normal tissues. Prior to 3D conformal radiation therapy planning, patient anatomy information was mostly limited to 2D radiological images or not even available. Generic patient CT images are often used in commercial radiotherapy treatment planning system (TPS) to reconstruct normal tissue doses. The objective of the current work was to develop a series of reference size computational human phantoms in DICOM-RT format for direct use in dose reconstruction in TPS. Methods: Contours of 93 organs and tissues were extracted from a series of pediatric and adult hybrid computational human phantoms (newborn, 1-, 5-, 10-, 15-year-old, and adult males and females) using Rhinoceros software. A MATLAB script was created to convert the contours into the DICOM-RT structure format. The simulated CT images with the resolution of 1×1×3 mm3 were also generated from the binary phantom format and coupled with the DICOM-structure files. Accurate volumes of the organs were drawn in the format using precise delineation of the contours in converted format. Due to complex geometry of organs, higher resolution (1×1×1 mm3) was found to be more efficient in the conversion of newborn and 1-year-old phantoms. Results: Contour sets were efficiently converted into DICOM-RT structures in relatively short time (about 30 minutes for each phantom). A good agreement was observed in the volumes between the original phantoms and the converted contours for large organs (NRMSD<1.0%) and small organs (NRMSD<7.7%). Conclusion: A comprehensive series of computational human phantoms in DICOM-RT format was created to support epidemiological studies of second cancer risks in radiotherapy patients. We confirmed the DICOM-RT phantoms were successfully imported into the TPS programs of major vendors.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-T-111: Investigation of the Attila Deterministic Solver as a Supplement to Monte Carlo for Calculating Out-Of-Field Radiotherapy Dose

Purpose: To use the Attila deterministic solver as a supplement to Monte Carlo for calculating out-of-field organ dose in support of epidemiological studies looking at the risks of second cancers. Supplemental dosimetry tools are needed to speed up dose calculations for studies involving large-scale patient cohorts. Methods: Attila is a multi-group discrete ordinates code which can solve the 3D photon-electron coupled linear Boltzmann radiation transport equation on a finite-element mesh. Dose is computed by multiplying the calculated particle flux in each mesh element by a medium-specific energy deposition cross-section. The out-of-field dosimetry capability of Attila is investigated by comparing average organ dose to that which is calculated by Monte Carlo simulation. The test scenario consists of a 6 MV external beam treatment of a female patient with a tumor in the left breast. The patient is simulated by a whole-body adult reference female computational phantom. Monte Carlo simulations were performed using MCNP6 and XVMC. Attila can export a tetrahedral mesh for MCNP6, allowing for a direct comparison between the two codes. The Attila and Monte Carlo methods were also compared in terms of calculation speed and complexity of simulation setup. A key perquisite for this work was the modeling of a Varian Clinac 2100 linear accelerator. Results: The solid mesh of the torso part of the adult female phantom for the Attila calculation was prepared using the CAD software SpaceClaim. Preliminary calculations suggest that Attila is a user-friendly software which shows great promise for our intended application. Computational performance is related to the number of tetrahedral elements included in the Attila calculation. Conclusion: Attila is being explored as a supplement to the conventional Monte Carlo radiation transport approach for performing retrospective patient dosimetry. The goal is for the dosimetry to be sufficiently accurate for use in retrospective epidemiological investigations.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-T-114: A Novel Anatomically Predictive Extension Model of Computational Human Phantoms for Dose Reconstruction in Retrospective Epidemiological Studies of Second Cancer Risks in Radiotherapy Patients

Purpose: Recent advances in cancer treatments have greatly increased the likelihood of post-treatment patient survival. Secondary malignancies, however, have become a growing concern. Epidemiological studies determining secondary effects in radiotherapy patients require assessment of organ-specific dose both inside and outside the treatment field. An essential input for Monte Carlo modeling of particle transport is radiological images showing full patient anatomy. However, in retrospective studies it is typical to only have partial anatomy from CT scans used during treatment planning. In this study, we developed a multi-step method to extend such limited patient anatomy to full body anatomy for estimating dose to normal tissues located outside the CT scan coverage. Methods: The first step identified a phantom from a library of body size-dependent computational human phantoms by matching the height and weight of patients. Second, a Python algorithm matched the patient CT coverage location in relation to the whole body phantom. Third, an algorithm cut the whole body phantom and scaled them to match the size of the patient. Then, merged the two anatomies into one whole body. We entitled this new approach, Anatomically Predictive Extension (APE). Results: The APE method was examined by comparing the original chest-abdomen-pelvis CT images of the five patients with the APE phantoms developed from only the chest part of the CAP images and whole body phantoms. We achieved average percent differences of tissue volumes of 25.7%, 34.2%, 16.5%, 26.8%, and 31.6% with an average of 27% across all patients. Conclusion: Our APE method extends the limited CT patient anatomy to whole body anatomy by using image processing and computational human phantoms. Our ongoing work includes evaluating the accuracy of these APE phantoms by comparing normal tissue doses in the APE phantoms and doses calculated for the original full CAP images under generic radiotherapy simulations. This research was supported by the NIH Intramural Research Program.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-T-115: Uncertainty in the Esophagus Dose in Retrospective Epidemiological Study of Breast Cancer Radiotherapy Patients

Purpose: Epidemiological studies of second cancer risks in breast cancer radiotherapy patients often use generic patient anatomy to reconstruct normal tissue doses when CT images of patients are not available. To evaluate the uncertainty involved in the dosimetry approach, we evaluated the esophagus dose in five sample patients by simulating breast cancer treatments. Methods: We obtained the diagnostic CT images of five anonymized adult female patients in different Body Mass Index (BMI) categories (16– 36kg/m2) from National Institutes of Health Clinical Center. We contoured the esophagus on the CT images and imported them into a Treatment Planning System (TPS) to create treatment plans and calculate esophagus doses. Esophagus dose was calculated once again via experimentally-validated Monte Carlo (MC) transport code, XVMC under the same geometries. We compared the esophagus doses from TPS and the MC method. We also investigated the degree of variation in the esophagus dose across the five patients and also the relationship between the patient characteristics and the esophagus doses. Results: Eclipse TPS using Analytical Anisotropic Algorithm (AAA) significantly underestimates the esophagus dose in breast cancer radiotherapy compared to MC. In the worst case, the esophagus dose from AAA was only 40% of the MC dose. The Coefficient of Variation across the patients was 48%. We found that the maximum esophagus dose was up to 2.7 times greater than the minimum. We finally observed linear relationship (Dose = 0.0218 × BMI – 0.1, R2=0.54) between patient’s BMI and the esophagus doses. Conclusion: We quantified the degree of uncertainty in the esophagus dose in five sample breast radiotherapy patients. The results of the study underscore the importance of individualized dose reconstruction for the study cohort to avoid misclassification in the risk analysis of second cancer. We are currently extending the number of patients up to 30.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-T-117: A Pilot Study of Organ Dose Reconstruction for Wilms Tumor Patients Treated with Radiation Therapy

Purpose: To reconstruct major organ doses for the Wilms tumor pediatric patients treated with radiation therapy using pediatric computational phantoms, treatment planning system (TPS), and Monte Carlo (MC) dose calculation methods. Methods: A total of ten female and male pediatric patients (15–88 months old) were selected from the National Wilms Tumor Study cohort and ten pediatric computational phantoms corresponding to the patient’s height and weight were selected for the organ dose reconstruction. Treatment plans were reconstructed on the computational phantoms in a Pinnacle TPS (v9.10) referring to treatment records and exported into DICOM-RT files, which were then used to generate the input files for XVMC MC code. The mean doses to major organs and the dose received by 50% of the heart were calculated and compared between TPS and MC calculations. The same calculations were conducted by replacing the computational human phantoms with a series of diagnostic patient CT images selected by matching the height and weight of the patients to validate the anatomical accuracy of the computational phantoms. Results: Dose to organs located within the treatment fields from the computational phantoms and the diagnostic patient CT images agreed within 2% for all cases for both TPS and MC calculations. The maximum difference of organ doses was 55.9 % (thyroid), but the absolute dose difference in this case was 0.33 Gy which was 0.96% of the prescription dose. The doses to ovaries and testes from MC in out-of-field provided more discrepancy (the maximum difference of 13.2% and 50.8%, respectively). The maximum difference of the 50% heart volume dose between the phantoms and the patient CT images was 40.0%. Conclusion: This study showed the pediatric computational phantoms are applicable to organ doses reconstruction for the radiotherapy patients whose three-dimensional radiological images are not available.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-T-119: Development of Heart Prediction Model to Increase Accuracy of Dose Reconstruction for Radiotherapy Patients

Purpose: To assess individual variation in heart volume and location in order to develop a prediction model of the heart. This heart prediction model will be used to calculate individualized heart doses for radiotherapy patients in epidemiological studies. Methods: Chest CT images for 30 adult male and 30 adult female patients were obtained from NIH Clinical Center. Image-analysis computer programs were used to segment the whole heart and 8 sub-regions and to measure the volume of each sub- region and the dimension of the whole heart. An analytical dosimetry method was used for the 30 adult female patients to estimate mean heart dose during conventional left breast radiotherapy. Results: The average volumes of the whole heart were 803.37 cm{sup 3} (COV 18.8%) and 570.19 cm{sup 3} (COV 18.8%) for adult male and female patients, respectively, which are comparable with the international reference volumes of 807.69 cm{sup 3} for males and 596.15 cm{sup 3} for females. Some patient characteristics were strongly correlated (R{sup 2}>0.5) with heart volume and heart dimensions (e.g., Body Mass Index vs. heart depth in males: R{sup 2}=0.54; weight vs. heart width in the adult females: R{sup 2}=0.63). We found that the mean heart dose 3.805 Gy (assuming prescribed dose of 50 Gy) in the breast radiotherapy simulations of the 30 adult females could be an underestimate (up to 1.6-fold) or overestimate (up to 1.8-fold) of the patient-specific heart dose. Conclusion: The study showed the significant variation in patient heart volumes and dimensions, resulting in substantial dose errors when a single average heart model is used for retrospective dose reconstruction. We are completing a multivariate analysis to develop a prediction model of the heart. This model will increase accuracy in dose reconstruction for radiotherapy patients and allow us to individualize heart dose calculations for patients whose CT images are not available.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-T-138: Commissioning and Evaluating Dose Computation Models for a Dedicated Proton Line Scanning Beam Nozzle in Eclipse Treatment Planning System

Purpose: In this study, we present an effective method to derive low dose envelope of the proton in-air spot fluence at beam positions other than the isocenter to reduce amount of measurements required for planning commission. Also, we demonstrate commissioning and validation results of this method to the Eclipse treatment planning system (version 13.0.29) for a Sumitomo dedicated proton line scanning beam nozzle. Methods: The in-air spot profiles at five beam-axis positions (±200, ±100 and 0 mm) were obtained in trigger mode using a MP3 Water tank (PTW-Freiburg) and a pinpoint ionization chamber (model 31014, PTW-Freiburg). Low dose envelope (below 1% of the center dose) of the spot profile at isocenter was obtained by repeated point measurements to minimize dosimetry uncertainty. The double Gaussian (DG) model was used to fit and obtain optimal σ1, σ2 and their corresponding weightings through our in-house MATLAB (Mathworks) program. σ1, σ2 were assumed to expand linearly along the beam axis from a virtual source position calculated by back projecting fitted sigmas from the single Gaussian (SG) model. Absolute doses in water were validated using an Advanced Markus chamber at the depth of 2cm with Pristine Peak (BP) R90d ranging from 5–32 cm for 10×10 cm2 scanned fields. The field size factors were verified with square fields from 2 to 20 cm at 2cm and before BP depth. Results: The absolute dose outputs were found to be within ±3%. For field size factor, the agreement between calculated and measurement were within ±2% at 2cm and ±3% before BP, except for the field size below 2×2 cm2. Conclusion: The double Gaussian model was found to be sufficient for characterizing the Sumitomo dedicated proton line scanning nozzle. With our effective double Gaussian fitting method, we are able to save significant proton beam time with acceptable output accuracy.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-T-405: Development of a Rapid Cardiac Contouring Tool Using Landmark-Driven Modeling

Purpose: This study aims to develop a tool to rapidly delineate cardiac substructures for use in dosimetry for large-scale clinical trial or epidemiological investigations. The goal is to produce a system that can semi-automatically delineate nine cardiac structures to a reasonable accuracy within a couple of minutes. Methods: The cardiac contouring tool employs a Most Similar Atlas method, where a selection criterion is used to pre-select the most similar model to the patient from a library of pre-defined atlases. Sixty contrast-enhanced cardiac computed tomography angiography (CTA) scans (30 male and 30 female) were manually contoured to serve as the atlas library. For each CTA 12 structures were delineated. Kabsch algorithm was used to compute the optimum rotation and translation matrices between the patient and atlas. Minimum root mean squared distance between the patient and atlas after transformation was used to select the most-similar atlas. An initial study using 10 CTA sets was performed to assess system feasibility. Leave-one patient out method was performed, and fit criteria were calculated to evaluate the fit accuracy compared to manual contours. Results: For the pilot study, mean dice indices of .895 were achieved for the whole heart, .867 for the ventricles, and .802 for the atria. In addition, mean distance was measured via the chord length distribution (CLD) between ground truth and the atlas structures for the four coronary arteries. The mean CLD for all coronary arteries was below 14mm, with the left circumflex artery showing the best agreement (7.08mm). Conclusion: The cardiac contouring tool is able to delineate cardiac structures with reasonable accuracy in less than 90 seconds. Pilot data indicates that the system is able to delineate the whole heart and ventricles within a reasonable accuracy using even a limited library. We are extending the atlas sets to 60 adult males and females in total.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-G-IeP2-10: Lens Dose Reduction by Patient Position Modification During Neck CT Exams

Purpose: Irradiation of the lens during a neck CT may increase a patient’s risk of developing cataracts later in life. Radiologists and technologists at the National Institutes of Health Clinical Center (NIHCC) have developed new CT imaging protocols that include a reduction in scan range and modifying neck positioning using a head tilt. This study will evaluate the efficacy of this protocol in the reduction of lens dose. Methods: We retrieved CT images of five male patients who had two sets of CT images: before and after the implementation of the new protocol. The lens doses before the new protocol were calculated using an in-house CT dose calculator, National Cancer Institute dosimetry system for CT (NCICT), where computational human phantoms with no head tilt are included. We also calculated the lens dose for the patient CT conducted after the new protocol by using an adult male computational phantom with the neck position deformed to match the angle of the head tilt. We also calculated the doses to other radiosensitive organs including the globes of the eye, brain, pituitary gland and salivary glands before and after head tilt. Results: Our dose calculations demonstrated that modifying neck position reduced dose to the lens by 89% on average (range: 86–96%). Globe, brain, pituitary and salivary gland doses also decreased by an average of 65% (51–95%), 38% (−8–66%), 34% (−43–84%) and 14% (13–14%), respectively. The new protocol resulted in a nearly ten-fold decrease in lens dose. Conclusion: The use of a head tilt and scan range reduction is an easy and effective method to reduce radiation exposure to the lens and other radiosensitive organs, while still allowing for the inclusion of critical neck structures in the CT image. We are expanding our study to a total of 10 males and 10 females.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-G-TeP4-06: An Integrated Application for Radiation Therapy Treatment Plan Directives, Management, and Reporting

Purpose: With electronic medical records, patient information for the treatment planning process has become disseminated across multiple applications with limited quality control and many associated failure modes. We present the development of a single application with a centralized database to manage the planning process. Methods: The system was designed to replace current functionalities of (i) static directives representing the physician intent for the prescription and planning goals, localization information for delivery, and other information, (ii) planning objective reports, (iii) localization and image guidance documents and (iv) the official radiation therapy prescription in the medical record. Using the Eclipse Scripting Application Programming Interface, a plug-in script with an associated domain-specific SQL Server database was created to manage the information in (i)–(iv). The system’s user interface and database were designed by a team of physicians, clinical physicists, database experts, and software engineers to ensure usability and robustness for clinical use. Results: The resulting system has been fully integrated within the TPS via a custom script and database. Planning scenario templates, version control, approvals, and logic-based quality control allow this system to fully track and document the planning process as well as physician approval of tradeoffs while improving the consistency of the data. Multiple plans and prescriptions are supported along with non-traditional dose objectives and evaluation such as biologically corrected models, composite dose limits, and management of localization goals. User-specific custom views were developed for the attending physician review, physicist plan checks, treating therapists, and peer review in chart rounds. Conclusion: A method was developed to maintain cohesive information throughout the planning process within one integrated system by using a custom treatment planning management application that interfaces directly with the TPS. Future work includes quantifying the improvements in quality, safety and efficiency that are possible with the routine clinical use of this system. Supported in part by NIH-P01-CA-059827.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

MO-FG-CAMPUS-TeP1-01: An Efficient Method of 3D Patient Dose Reconstruction Based On EPID Measurements for Pre-Treatment Patient Specific QA

Purpose: To demonstrate an efficient and clinically relevant patient specific QA method by reconstructing 3D patient dose from 2D EPID images for IMRT plans. Also to determine the usefulness of 2D QA metrics when assessing 3D patient dose deviations. Methods: Using the method developed by King et al (Med Phys 39(5),2839–2847), EPID images of IMRT fields were acquired in air and converted to dose at 10 cm depth (SAD setup) in a flat virtual water phantom. Each EPID measured dose map was then divided by the corresponding treatment planning system (TPS) dose map calculated with an identical setup, to derive a 2D “error matrix”. For each field, the error matrix was used to adjust the doses along the respective ray lines in the original patient 3D dose. All field doses were combined to derive a reconstructed 3D patient dose for quantitative analysis. A software tool was developed to efficiently implement the entire process and was tested with a variety of IMRT plans for 2D (virtual flat phantom) and 3D (in-patient) QA analysis. Results: The method was tested on 60 IMRT plans. The mean (± standard deviation) 2D gamma (2%,2mm) pass rate (2D-GPR) was 97.4±3.0% and the mean 2D gamma index (2D-GI) was 0.35±0.06. The 3D PTV mean dose deviation was 1.8±0.8%. The analysis showed very weak correlations between both the 2D-GPR and 2D-GI when compared with PTV mean dose deviations (R2=0.3561 and 0.3632 respectively). Conclusion: Our method efficiently calculates 3D patient dose from 2D EPID images, utilising all of the advantages of an EPID-based dosimetry system. In this study, the 2D QA metrics did not predict the 3D patient dose deviation. This tool allows reporting of the 3D volumetric dose parameters thus providing more clinically relevant patient specific QA.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TU-H-CAMPUS-TeP1-04: Novel 3D Printed Plastic Cutouts Filled with Aluminum Oxide for Same Day Electron Radiotherapy

Purpose: Clinics that outsource electron cutout manufacturing may be unable to simulate and treat patients on the same day. To enable same day treatment, we investigate the use of 3D printed hollow cutouts filled with 30 grit Al{sub 2}O{sub 3} powder. We verified the dosimetric equivalence of such a cutout relative to an outsourced Copper cutout. Methods: Printing was performed using a Ultibots Kossel 250 V-Slot 3D printer and polylactic acid filament. Printing files were derived from an in-house 3D model designed to mate with a Varian 6 cm electron cone. Relative to conventional cutouts (Copper or Cerrobend), the height of the hollow plastic cutout was extended by 1.0 cm to increase attenuation. Measurements were performed for 6 MeV in solid water at dmax (1.4 cm) with Gafchromic™ EBT3 film; the cutout was kidney-shaped with a long and short axis of approximately 5 and 2 cm, respectively. The Copper cutout was based on an outline of the 3D printed cutout.A calibration film was exposed immediately after the electron irradiations. All films, including an un-irradiated one, were from the same batch. Films were scanned on an Epson 10000XL flatbed scanner. Film analysis was performed in DoseLab (MOBIUS Medical Systems, Houston, Tx). Results: Visual comparison of the physical cutouts revealed that the Copper cutout had a slightly smaller opening than the printed cutout. Line profiles through the registered films indicated agreement within 5% in the open section. 99.9% of pixels passed gamma analysis with 2% local percent difference, 2 mm DTA, and a 25% threshold. Conclusion: Same day simulation and treatment with electrons is feasible with 3D printing of a hollow cutout filled with Al{sub 2}O{sub 3}. Future work will include evaluations of additional cutout shapes at different depths for higher energies. Other printing materials, such as bismuth, are being tested.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TH-AB-202-06: BEST IN PHYSICS (JOINT IMAGING-THERAPY): A Real-Time Tumor Tracking Using Novel Scatter Imaging Modality During Lung SBRT

Purpose: Compared to traditional radiotherapy techniques, stereotactic body radiation therapy (SBRT) provides more favorable outcomes during the treatment of certain lung tumors. Despite advancements in image guidance, accurate target localization still remains a challenge. In this work, we expand our knowledge of a novel scatter imaging modality in order to develop a real-time tumor localization method using scattered photons from the patient during treatment. Methods: Images of the QUASAR™ Respiratory Motion Phantom were taken by irradiating it on a Varian TrueBeam accelerator. The scattered radiation was detected using a flat panel-based pinhole camera detection system. Two motion settings were investigated: static and dynamic. In the former, the lung tumor was manually shifted between imaging. In the latter, the lung tumor was set to move at a certain frequency and amplitude while the images were acquired continuously for one minute. The accuracy of tumor localization and the irradiation time required to distinguish the lung tumor were studied. Results: The comparison of measured and expected location of the lung tumor during static motion was shown to be under standard deviation (STD) of 0.064 with a mean STD of 0.031cm. The dynamic motion was taken at a rate of 1400 MU/min for one minute and the measured location of the lung tumor was then compared with the QUASAR phantom’s sinusoidal motion pattern and the agreement found to be at an average STD of 0.275cm. The location of the lung tumor was investigated using aggregate images consisting of 1 or 2 frames/image and the change was below STD of 0.30cm. The lung tumor also appeared to be blurrier in images consisting of two frames. Conclusion: Based on our preliminary results real-time image guidance using the scatter imaging modality to localize and track tumors during lung SBRT has the potential to become clinical reality.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TH-CD-207A-08: Simulated Real-Time Image Guidance for Lung SBRT Patients Using Scatter Imaging

Purpose: To develop a comprehensive Monte Carlo-based model for the acquisition of scatter images of patient anatomy in real-time, during lung SBRT treatment. Methods: During SBRT treatment, images of patient anatomy can be acquired from scattered radiation. To rigorously examine the utility of scatter images for image guidance, a model is developed using MCNP code to simulate scatter images of phantoms and lung cancer patients. The model is validated by comparing experimental and simulated images of phantoms of different complexity. The differentiation between tissue types is investigated by imaging objects of known compositions (water, lung, and bone equivalent). A lung tumor phantom, simulating materials and geometry encountered during lung SBRT treatments, is used to investigate image noise properties for various quantities of delivered radiation (monitor units(MU)). Patient scatter images are simulated using the validated simulation model. 4DCT patient data is converted to an MCNP input geometry accounting for different tissue composition and densities. Lung tumor phantom images acquired with decreasing imaging time (decreasing MU) are used to model the expected noise amplitude in patient scatter images, producing realistic simulated patient scatter images with varying temporal resolution. Results: Image intensity in simulated and experimental scatter images of tissue equivalent objects (water, lung, bone) match within the uncertainty (∼3%). Lung tumor phantom images agree as well. Specifically, tumor-to-lung contrast matches within the uncertainty. The addition of random noise approximating quantum noise in experimental images to simulated patient images shows that scatter images of lung tumors can provide images in as fast as 0.5 seconds with CNR∼2.7. Conclusions: A scatter imaging simulation model is developed and validated using experimental phantom scatter images. Following validation, lung cancer patient scatter images are simulated. These simulated patient images demonstrate the clinical utility of scatter imaging for real-time tumor tracking during lung SBRT.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

SU-D-BRD-03: Improving Plan Quality with Automation of Treatment Plan Checks

Purpose: To evaluate the effectiveness of an automated plan check tool to improve first-time plan quality as well as standardize and document performance of physics plan checks. Methods: The Plan Checker Tool (PCT) uses the Eclipse Scripting API to check and compare data from the treatment planning system (TPS) and treatment management system (TMS). PCT was created to improve first-time plan quality, reduce patient delays, increase efficiency of our electronic workflow, and to standardize and partially automate plan checks in the TPS. A framework was developed which can be configured with different reference values and types of checks. One example is the prescribed dose check where PCT flags the user when the planned dose and the prescribed dose disagree. PCT includes a comprehensive checklist of automated and manual checks that are documented when performed by the user. A PDF report is created and automatically uploaded into the TMS. Prior to and during PCT development, errors caught during plan checks and also patient delays were tracked in order to prioritize which checks should be automated. The most common and significant errors were determined. Results: Nineteen of 33 checklist items were automated with data extracted with the PCT. These include checks for prescription, reference point and machine scheduling errors which are three of the top six causes of patient delays related to physics and dosimetry. Since the clinical roll-out, no delays have been due to errors that are automatically flagged by the PCT. Development continues to automate the remaining checks. Conclusion: With PCT, 57% of the physics plan checklist has been partially or fully automated. Treatment delays have declined since release of the PCT for clinical use. By tracking delays and errors, we have been able to measure the effectiveness of automating checks and are using this information to prioritize future development. This project was supported in part by P01CA059827.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

SU-E-I-10: Investigation On Detectability of a Small Target for Different Slice Direction of a Volumetric Cone Beam CT Image

Purpose: To investigate the detectability of a small target for different slice direction of a volumetric cone beam CT image and its impact on dose reduction. Methods: Analytic projection data of a sphere object (1 mm diameter, 0.2/cm attenuation coefficient) were generated and reconstructed by FDK algorithm. In this work, we compared the detectability of the small target from four different backprojection Methods: hanning weighted ramp filter with linear interpolation (RECON 1), hanning weighted ramp filter with Fourier interpolation (RECON2), ramp filter with linear interpolation (RECON 3), and ramp filter with Fourier interpolation (RECON4), respectively. For noise simulation, 200 photons per measurement were used, and the noise only data were reconstructed using FDK algorithm. For each reconstructed volume, axial and coronal slice were extracted and detection-SNR was calculated using channelized Hotelling observer (CHO) with dense difference-of-Gaussian (D-DOG) channels. Results: Detection-SNR of coronal images varies for different backprojection methods, while axial images have a similar detection-SNR. Detection-SNR{sup 2} ratios of coronal and axial images in RECON1 and RECON2 are 1.33 and 1.15, implying that the coronal image has a better detectability than axial image. In other words, using coronal slices for the small target detection can reduce the patient dose about 33% and 15% compared to using axial slices in RECON 1 and RECON 2. Conclusion: In this work, we investigated slice direction dependent detectability of a volumetric cone beam CT image. RECON 1 and RECON 2 produced the highest detection-SNR, with better detectability in coronal slices. These results indicate that it is more beneficial to use coronal slice to improve detectability of a small target in a volumetric cone beam CT image. This research was supported by the MSIP (Ministry of Science, ICT and Future Planning), Korea, under the IT Consilience Creative Program (NIPA-2014-H0201-14-1002) supervised by the NIPA (National IT Industry Promotion Agency). Authors declares that s/he has no conflict of Interest in relation to the work in this abstract.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

SU-E-J-93: Development of Pre-Contoured Human Model Library in DICOM-RT Format for the Epidemiological Study of the Radiotherapy Patients

Purpose: Prior to 3D conformal radiation therapy planning, patient anatomy information was mostly limited to 2D beams-eye-view from the conventional simulator. To analyze the outcomes of such treatments for radiation late effects, 3D computational human models are often used in commercial treatment planning systems (TPSs). However, several underlying difficulties such as time-consuming manual delineation procedures of a large number of structures in the model have always limited its applications. Primary objective of this work was to develop a human model library for the epidemiological study by converting 3D-surface model organs to DICOM-RT format (DICOM-RT structure) using an in-house built software. We converted the ICRP reference human models to DICOM-RT models, which can be readily adopted for various dose calculations. Methods: MATLAB based code were utilized to convert the contour drawings extracted in text-format from the 3D graphic-tool, Rhinoceros into DICOM-RT structure format for 50 different organs of each model using a 16GB dual-core processor. The conversion periods were measured for each DICOM-RT models, and the reconstructed structure volumes were validated against the original 3D-surface models in the TPS. Ten reference hybrid whole-body models (8-pediatric and 2-adults) were automatically processed to create DICOM-RT computational human model library. Results: Mean contour conversion period was found to be 580 (N=2) and 394.5 (N=8) seconds for 50 organs in the adult and pediatric models respectively. A good agreement for large organs (NRMSD <1.0%) and small organs (NRMSD <7.7%) was also observed between the original volumes and corresponding DICOM-RT structure volumes of the organs. Conclusion: The ICRP reference human models were converted into DICOM-RT format to support the epidemiological study using a large cohort of conventional radiotherapy patients. Due to its DICOM-compatibility, the library may be implemented to many other different applications. We also expect to develop the library by including additional models in future.

60 APPLIED LIFE SCIENCES↗