Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HISTOGRAM”

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 163 records · Page 9

The first clinical implementation of electromagnetic transponder-guided MLC tracking

Purpose: We report on the clinical process, quality assurance, and geometric and dosimetric results of the first clinical implementation of electromagnetic transponder-guided MLC tracking which occurred on 28 November 2013 at the Northern Sydney Cancer Centre. Methods: An electromagnetic transponder-based positioning system (Calypso) was modified to send the target position output to in-house-developed MLC tracking code, which adjusts the leaf positions to optimally align the treatment beam with the real-time target position. Clinical process and quality assurance procedures were developed and performed. The first clinical implementation of electromagnetic transponder-guided MLC tracking was for a prostate cancer patient being treated with dual-arc VMAT (RapidArc). For the first fraction of the first patient treatment of electromagnetic transponder-guided MLC tracking we recorded the in-room time and transponder positions, and performed dose reconstruction to estimate the delivered dose and also the dose received had MLC tracking not been used. Results: The total in-room time was 21 min with 2 min of beam delivery. No additional time was needed for MLC tracking and there were no beam holds. The average prostate position from the initial setup was 1.2 mm, mostly an anterior shift. Dose reconstruction analysis of the delivered dose with MLC tracking showed similar isodose and target dose volume histograms to the planned treatment and a 4.6% increase in the fractional rectal V{sub 60}. Dose reconstruction without motion compensation showed a 30% increase in the fractional rectal V{sub 60} from that planned, even for the small motion. Conclusions: The real-time beam-target correction method, electromagnetic transponder-guided MLC tracking, has been translated to the clinic. This achievement represents a milestone in improving geometric and dosimetric accuracy, and by inference treatment outcomes, in cancer radiotherapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Impact of dose size in single fraction spatially fractionated (grid) radiotherapy for melanoma

Purpose: To evaluate the impact of dose size in single fraction, spatially fractionated (grid) radiotherapy for selectively killing infiltrated melanoma cancer cells of different tumor sizes, using different radiobiological models. Methods: A Monte Carlo technique was employed to calculate the 3D dose distribution of a commercially available megavoltage grid collimator in a 6 MV beam. The linear-quadratic (LQ) and modified linear quadratic (MLQ) models were used separately to evaluate the therapeutic outcome of a series of single fraction regimens that employed grid therapy to treat both acute and late responding melanomas of varying sizes. The dose prescription point was at the center of the tumor volume. Dose sizes ranging from 1 to 30 Gy at 100% dose line were modeled. Tumors were either touching the skin surface or having their centers at a depth of 3 cm. The equivalent uniform dose (EUD) to the melanoma cells and the therapeutic ratio (TR) were defined by comparing grid therapy with the traditional open debulking field. The clinical outcomes from recent reports were used to verify the authors’ model. Results: Dose profiles at different depths and 3D dose distributions in a series of 3D melanomas treated with grid therapy were obtained. The EUDs and TRs for all sizes of 3D tumors involved at different doses were derived through the LQ and MLQ models, and a practical equation was derived. The EUD was only one fifth of the prescribed dose. The TR was dependent on the prescribed dose and on the LQ parameters of both the interspersed cancer and normal tissue cells. The results from the LQ model were consistent with those of the MLQ model. At 20 Gy, the EUD and TR by the LQ model were 2.8% higher and 1% lower than by the MLQ, while at 10 Gy, the EUD and TR as defined by the LQ model were only 1.4% higher and 0.8% lower, respectively. The dose volume histograms of grid therapy for a 10 cm tumor showed different dosimetric characteristics from those of conventional radiotherapy. A significant portion of the tumor volume received a very large dose in grid therapy, which ensures significant tumor cell killing in these regions. Conversely, some areas received a relatively small dose, thereby sparing interspersed normal cells and increasing radiation tolerance. The radiobiology modeling results indicated that grid therapy could be useful for treating acutely responding melanomas infiltrating radiosensitive normal tissues. The theoretical model predictions were supported by the clinical outcomes. Conclusions: Grid therapy functions by selectively killing infiltrating tumor cells and concomitantly sparing interspersed normal cells. The TR depends on the radiosensitivity of the cell population, dose, tumor size, and location. Because the volumes of very high dose regions are small, the LQ model can be used safely to predict the clinical outcomes of grid therapy. When treating melanomas with a dose of 15 Gy or higher, single fraction grid therapy is clearly advantageous for sparing interspersed normal cells. The existence of a threshold fraction dose, which was found in the authors’ theoretical simulations, was confirmed by clinical observations.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Dose impact in radiographic lung injury following lung SBRT: Statistical analysis and geometric interpretation

Purpose: To demonstrate a new method of evaluating dose response of treatment-induced lung radiographic injury post-SBRT (stereotactic body radiotherapy) treatment and the discovery of bimodal dose behavior within clinically identified injury volumes. Methods: Follow-up CT scans at 3, 6, and 12 months were acquired from 24 patients treated with SBRT for stage-1 primary lung cancers or oligometastic lesions. Injury regions in these scans were propagated to the planning CT coordinates by performing deformable registration of the follow-ups to the planning CTs. A bimodal behavior was repeatedly observed from the probability distribution for dose values within the deformed injury regions. Based on a mixture-Gaussian assumption, an Expectation-Maximization (EM) algorithm was used to obtain characteristic parameters for such distribution. Geometric analysis was performed to interpret such parameters and infer the critical dose level that is potentially inductive of post-SBRT lung injury. Results: The Gaussian mixture obtained from the EM algorithm closely approximates the empirical dose histogram within the injury volume with good consistency. The average Kullback-Leibler divergence values between the empirical differential dose volume histogram and the EM-obtained Gaussian mixture distribution were calculated to be 0.069, 0.063, and 0.092 for the 3, 6, and 12 month follow-up groups, respectively. The lower Gaussian component was located at approximately 70% prescription dose (35 Gy) for all three follow-up time points. The higher Gaussian component, contributed by the dose received by planning target volume, was located at around 107% of the prescription dose. Geometrical analysis suggests the mean of the lower Gaussian component, located at 35 Gy, as a possible indicator for a critical dose that induces lung injury after SBRT. Conclusions: An innovative and improved method for analyzing the correspondence between lung radiographic injury and SBRT treatment dose has been demonstrated. Bimodal behavior was observed in the dose distribution of lung injury after SBRT. Novel statistical and geometrical analysis has shown that the systematically quantified low-dose peak at approximately 35 Gy, or 70% prescription dose, is a good indication of a critical dose for injury. The determined critical dose of 35 Gy resembles the critical dose volume limit of 30 Gy for ipsilateral bronchus in RTOG 0618 and results from previous studies. The authors seek to further extend this improved analysis method to a larger cohort to better understand the interpatient variation in radiographic lung injury dose response post-SBRT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SU-E-J-93: Parametrisation of Dose to the Mucosa of the Anterior Rectal Wall in Transrectal Ultrasound Guided High-Dose-Rate Brachytherapy of the Prostate

Purpose: In high-dose-rate (HDR) brachytherapy of the prostate, radiation is delivered from a number of radioactive sources which are inserted via catheter into the target volume. The rectal mucosa also receives dose during the treatment, which may lead to late toxicity effects. To allow possible links between rectal dose and toxicity to be investigated, suitable methods of parametrising the rectal dose are needed. Methods: During treatment of a series of 95 patients, anatomy and catheter locations were monitored by transrectal ultrasound, and target volume positions were contoured on the ultrasound scan by the therapist. The anterior rectal mucosal wall was identified by contouring the transrectal ultrasound balloon within the ultrasound scan. Source positions and dwell times, along with the dose delivered to the patient were computed using the Oncentra Prostate treatment planning system (TPS). Data for the series of patients were exported from the TPS in Dicom format, and a series of parametrisation methods were developed in a Matlab environment to assess the rectal dose. Results: Contours of the anterior rectal mucosa were voxelised within Matlab to allow the dose to the rectal mucosa to be analysed directly from the 3D dose grid. Dose parametrisations based on dose-surface (DSH) and dose-line (DLH) histograms were obtained. Both lateral and longitudinal extents of the mucosal dose were parametrised using dose-line histograms in the relevant directions. Conclusion: We have developed a series of dose parametrisations for quantifying the dose to the rectal mucosa during HDR prostate brachytherapy which are suitable for future studies investigating potential associations between mucosal dose and late toxicity effects. The geometry of the transrectal probe standardises the rectal anatomy, making this treatment technique particularly suited to studies of this nature.

60 APPLIED LIFE SCIENCES↗

SU-E-T-275: Radiobiological Evaluation of Intensity Modulated Radiotherapy Treatment for Locally Advanced Head and Neck Squamous Cell Carcinomas

Purpose: To evaluate the radiobiological outcome of Intensity Modulated Radiotherapy Treatment (IMRT) for locally advanced head and neck squamous cell carcinomas using HART (Histogram Analysis in Radiation Therapy; J Appl Clin Med Phys 11(1): 137–157, 2010) program and compare with the clinical outcomes. Methods: We have treated 20 patients of stage III and IV HNSCC Oropharynx and hypopharynx with accelerated IMRT technique and concurrent chemotherapy. Delineation of tumor and normal tissues were done using Danish Head and Neck Cancer Group (DAHANCA) contouring guidelines and radiotherapy was delivered to a dose of 70Gy in 35 fractions to the primary and involved lymph nodes, 63Gy to intermediate risk areas and 56 Gy to lower risk areas, Monday to Saturday, 6 Days/week using 6 MV Photons with an expected overall treatment time of 6 weeks. The TCP and NTCP's were calculated from the dose-volume histogram (DVH) statistics using the Poisson Statistics (PS) and JT Lyman models respectively and the Resultwas correlated with clinical outcomes of the patients with mean follow up of 24 months. Results: Using HART program, the TCP (0.89± 0.01) of primary tumor and the NTCP for parotids (0.20±0.12), spinal cord (0.05±0.01), esophagus (0.30±0.2), mandible (0.35±0.21), Oral cavity (0.37±0.18), Larynx (0.30±0.15) were estimated and correlated with clinical outcome of the patients. Conclusion: Accelerated IMRT with Chemotherapy is a clinical feasible option in the treatment of locally advanced HNSCC with encouraging initial tumour response and acceptable acute toxicities. The correlation between the clinical outcomes and radiobiological model estimated parameters using HART programs are found to be satisfactory.

60 APPLIED LIFE SCIENCES↗

SU-E-T-375: Passive Scattering to Pencil-Beam-Scanning Comparison for Medulloblastoma Proton Therapy: LET Distributions and Radiobiological Implications

Purpose: To compare the linear energy transfer (LET) distributions between passive scattering and pencil beam scanning proton radiation therapy techniques for medulloblastoma patients and study the potential radiobiological implications. Methods: A group of medulloblastoma patients, previously treated with passive scattering (PS) proton craniospinal irradiation followed by prosterior fossa or involved field boost, were selected from the patient database of our institution. Using the beam geometry and planning computed tomography (CT) image sets of the original treatment plans, pencil beam scanning (PBS) treatment plans were generated for the cranial treatment for each patient, with average beam spot size of 8mm (sigma in air at isocenter). 3-dimensional dose and LET distributions were calculated by Monte Carlo methods (TOPAS) both for the original passive scattering and new pencil beam scanning treatment plans. LET volume histograms were calculated for the target and OARs and compared for the two delivery methods. Variable RBE weighted dose distributions and volume histograms were also calculated using a variable dose and LET-based model. Results: Better dose conformity was achieved with PBS planning compared to PS, leading to increased dose coverage for the boost target area and decreased average dose to the structures adjacent to it and critical structures outside the whole brain treatment field. LET values for the target were lower for PBS plans. Elevated LET values for OARs close to the boosted target areas were noticed, due to end of range of proton beams falling inside these structures, resulting in higher RBE weighted dose for these structures compared to the clinical RBE value of 1.1. Conclusion: Transitioning from passive scattering to pencil beam scanning proton radiation treatment can be dosimetrically beneficial for medulloblastoma patients. LET–guided treatment planning could contribute to better decision making for these cases, especially for critical structures at close proximity to the boosted target area.

60 APPLIED LIFE SCIENCES↗

SU-E-T-604: Dosimetric Dependence On the Collimator Angle in Prostate Volumetric Modulated Arc Therapy

Purpose: The purpose of this study is to investigate the dose-volume variations of planning target volume (PTV) and organs-at-risk (OARs) in prostate volumetric modulated arc therapy (VMAT) when using different collimator angles. It is because collimator angle awareness is essential for planner to produce an optimal prostate VMAT plan in a rational time. Methods: Single-arc VMAT plans at different collimator angles (0o, 15o, 30o, 45o, 60o, 75o and 90o) were created systematically using a Harold heterogeneous pelvis phantom. For each change of collimator angle, a new plan was re-optimized for that angle. The prescription dose was 78 Gy per 39 fractions. Conformity index (CI), homogeneity index (HI), gradient index, machine monitor unit, dose-volume histogram, the mean and maximum doses of the PTV were calculated and analyzed. On the other hand, dose-volume histogram, the mean and maximum doses of the OARs such as bladder, rectum and femoral heads for different collimator angles were determined from the plans. Results: There was no significance difference, based on the plan dose-volume evaluation criteria, found in the VMAT optimizations for all studied collimator angles. Higher CI and lower HI were found for the 45o collimator angle. In addition, the 15o collimator angle provided lower HI similar to the 45o collimator angle. The 75o and 90o collimator angle were found good for the rectum sparing, and the 75o and 30o collimator angle were found good for the right and left femur sparing, respectively. The PTV dose coverage for each plan was comparatively independent of the collimator angle. Conclusion: The dosimetric results in this study are useful to the planner to select different collimator angles to improve the PTV coverage and OAR sparing in prostate VMAT.

60 APPLIED LIFE SCIENCES↗

SU-E-T-632: A Dosimetric Comparison of the 3D-CRT Planning of Chest Wall in Post-Mastectomy Breast Cancer Patients, with and Without Breast Board Setup

Purpose: Breast boards are used in breast radiation which increases normal lung and heart doses, when supraclavicular field is included. Therefore, in this study through dose volume histogram (DVHs), lung and heart doses comparison was done between two different setups i.e. with and without breast board, for the treatment of left chest wall and supraclavicular fossa in postmastectomy left breast cancer. Methods: In this study, CT-Simulation scans of ten breast cancer patients were done with and without breast board, at Shifa International Hospitals Islamabad, to investigate the differences between the two different setups of the irradiation of left chest wall in terms of lung and heart doses. For immobilization, support under the neck, shoulders and arms was used. Precise PLAN 2.15 treatment planning system (TPS) was used for 3D-CRT planning. The total prescribed dose for both the plans was 5000 cGy/25 fractions. The chest wall was treated with a pair of tangential photon fields and the upper supraclavicular nodal regions were treated with an anterior photon field. A mono-isocentric technique was used to match the tangential fields with the anterior field at the isocentre. The dose volume histogram was used to compare the doses of heart and ipsilateral lung. Results: Both the plans of each patient were generated and compared. DVH results showed that for the same PTV dose coverage, plans without breast board resulted in a reduction of lung and heart doses compared with the plans with breast board. There was significant reductions in V20, V<25 and mean doses for lung and V<9 and mean doses for heart. Conclusion: In comparison of both the plans, setup without breast board significantly reduced the dose-volume of the ipsilateral lung and heart in left chest wall patients. Waived registration request has been submitted.

60 APPLIED LIFE SCIENCES↗

SU-E-T-634: Analysis of Volume Based GYN HDR Brachytherapy Plans for Dose Calculation to Organs At Risk(OAR)

Purpose: We have analyzed the dose volume histogram of 140 CT based HDR brachytherapy plans and evaluated the dose received to OAR ; rectum, bladder and sigmoid colon based on recommendations from ICRU and Image guided brachytherapy working group for cervical cancer . Methods: Our treatment protocol consist of XRT to whole pelvis with 45 Gy at 1.8Gy/fraction followed by 30 Gy at 6 Gy per fraction by HDR brachytherapy in 2 weeks . The CT compatible tandem and ovoid applicators were used and stabilized with radio opaque packing material. The patient was stabilized using special re-locatable implant table and stirrups for reproducibility of the geometry during treatment. The CT scan images were taken at 3mm slice thickness and exported to the treatment planning computer. The OAR structures, bladder, rectum and sigmoid colon were outlined on the images along with the applicators. The prescription dose was targeted to A left and A right as defined in Manchester system and optimized on geometry . The dosimetry was compared on all plans using the parameter Ci.sec.cGy-1 . Using the Dose Volume Histogram (DVH) obtained from the plans the doses to rectum, sigmoid colon and bladder for ICRU defined points and 2cc volume were analyzed and reported. The following criteria were used for limiting the tolerance dose by volume (D2cc) were calculated. The rectum and sigmoid colon doses were limited to <75Gy. The bladder dose was limited to < 90Gy from both XRT and HDR brachytherapy. Results: The average total (XRT+HDRBT) BED values to prescription volume was 120 Gy. Dose 2cc to rectum was 70Gy +/− 17Gy, dose to 2cc bladder was 82+/−32 Gy. The average Ci.sec.cGy-1 calculated for the HDR plans was 6.99 +/− 0.5 Conclusion: The image based treatment planning enabled to evaluati volume based dose to critical structures for clinical interpretation.

60 APPLIED LIFE SCIENCES↗

SU-F-BRD-11: Prediction of Dosimetric Endpoints From Patient Geometry Using Neural Nets

Purpose: The previously-published overlap volume histogram (OVH) technique lends itself naturally to prediction of the dose received by a given volume of tissue (e.g. D90) in intensity-modulated radiotherapy (IMRT) treatment plans. Here we extend the OVH technique using artificial neural networks in order to predict the volume of tissue receiving a given dose (e.g. V90) in both prostate IMRT and conventional breast radiotherapy. Methods: Twenty-nine prostate treatment plans and forty-three breast treatment plans were analyzed. The spatial relationships between the prostate and rectum and between the breast and ipsilateral lung were characterized using OVHs. The OVH is a cumulative histogram representing the fractional volume of the risk organ overlapped by a series of isotropic expansions of the planning target volume (PTV). Seven cases were identified as outliers and replanned. OVH points were used as inputs to a one hidden layer feed forward artificial neural network with quality parameters of the corresponding plan, such as the rectum V50, as targets. A 3-fold cross-validation was used to estimate the prediction error. Results: The root mean square (RMS) error between the predicted rectum V50s and the planned values was 2.3, which was 35% of the standard deviation of V50 for the twenty-nine plans. The RMS error of prediction of V20 of the ipsilateral lung in breast cases was 3.9, which was 90% of the standard deviation of the V20 values in the breast plan database. Conclusion: This study demonstrates that artificial neural nets can be used to extend the OVH technique to predict dosimetric endpoints taking the form of a volume receiving a given dose, rather than the minimum dose received by a given volume. Prediction of ipsilateral lung dose in breast radiotherapy using the OVH technique remains a work in progress.

60 APPLIED LIFE SCIENCES↗

MO-G-BRE-03: Automated Continuous Monitoring of Patient Setup with Second-Check Independent Image Registration

Purpose: To create a non-supervised quality assurance program to monitor image-based patient setup. The system acts a secondary check by independently computing shifts and rotations and interfaces with Varian's database to verify therapist's work and warn against sub-optimal setups. Methods: Temporary digitally-reconstructed radiographs (DRRs) and OBI radiographic image files created by Varian's treatment console during patient setup are intercepted and used as input in an independent registration module customized for accuracy that determines the optimal rotations and shifts. To deal with the poor quality of OBI images, a histogram equalization of the live images to the DDR counterparts is performed as a pre-processing step. A search for the most sensitive metric was performed by plotting search spaces subject to various translations and convergence analysis was applied to ensure the optimizer finds the global minima. Final system configuration uses the NCC metric with 150 histogram bins and a one plus one optimizer running for 2000 iterations with customized scales for translations and rotations in a multi-stage optimization setup that first corrects and translations and subsequently rotations. Results: The system was installed clinically to monitor and provide almost real-time feedback on patient positioning. On a 2 month-basis uncorrected pitch values were of a mean 0.016° with standard deviation of 1.692°, and couch rotations of − 0.090°± 1.547°. The couch shifts were −0.157°±0.466° cm for the vertical, 0.045°±0.286 laterally and 0.084°± 0.501° longitudinally. Uncorrected pitch angles were the most common source of discrepancies. Large variations in the pitch angles were correlated with patient motion inside the mask. Conclusion: A system for automated quality assurance of therapist's registration was designed and tested in clinical practice. The approach complements the clinical software's automated registration in terms of algorithm configuration and performance and constitutes a practical approach to implement safe and cost-effective radiotherapy.

60 APPLIED LIFE SCIENCES↗

WE-E-17A-05: Complementary Prognostic Value of CT and 18F-FDG PET Non-Small Cell Lung Cancer Tumor Heterogeneity Features Quantified Through Texture Analysis

Purpose: Previous studies have shown that CT or 18F-FDG PET intratumor heterogeneity features computed using texture analysis may have prognostic value in Non-Small Cell Lung Cancer (NSCLC), but have been mostly investigated separately. The purpose of this study was to evaluate the potential added value with respect to prognosis regarding the combination of non-enhanced CT and 18F-FDG PET heterogeneity textural features on primary NSCLC tumors. Methods: One hundred patients with non-metastatic NSCLC (stage I–III), treated with surgery and/or (chemo)radiotherapy, that underwent staging 18F-FDG PET/CT images, were retrospectively included. Morphological tumor volumes were semi-automatically delineated on non-enhanced CT using 3D SlicerTM. Metabolically active tumor volumes (MATV) were automatically delineated on PET using the Fuzzy Locally Adaptive Bayesian (FLAB) method. Intratumoral tissue density and FDG uptake heterogeneities were quantified using texture parameters calculated from co-occurrence, difference, and run-length matrices. In addition to these textural features, first order histogram-derived metrics were computed on the whole morphological CT tumor volume, as well as on sub-volumes corresponding to fine, medium or coarse textures determined through various levels of LoG-filtering. Association with survival regarding all extracted features was assessed using Cox regression for both univariate and multivariate analysis. Results: Several PET and CT heterogeneity features were prognostic factors of overall survival in the univariate analysis. CT histogram-derived kurtosis and uniformity, as well as Low Grey-level High Run Emphasis (LGHRE), and PET local entropy were independent prognostic factors. Combined with stage and MATV, they led to a powerful prognostic model (p<0.0001), with median survival of 49 vs. 12.6 months and a hazard ratio of 3.5. Conclusion: Intratumoral heterogeneity quantified through textural features extracted from both CT and FDG PET images have complementary and independent prognostic value in NSCLC.

60 APPLIED LIFE SCIENCES↗

WE-G-16A-01: Evolution of Radiation Treatment Planning

Welcome and Introduction - Lawrence N. Rothenberg This symposium is one a continuing series of presentations at AAPM Annual Meetings on the historical aspects of medical physics, radiology, and radiation oncology that have been organized by the AAPM History Committee. Information on previous presentations including “Early Developments in Teletherapy” (Indianapolis 2013), “Historical Aspects of Cross-Sectional Imaging” (Charlotte 2012), “Historical Aspects of Brachytherapy” (Vancouver 2011), “50 Years of Women in Medical Physics” (Houston 2008), and “Roentgen's Early Investigations” (Minneapolis 2007) can be found in the Education Section of the AAPM Website. The Austin 2014 History Symposium will be on “Evolution of Radiation Treatment Planning.” Overview - Radhe Mohan Treatment planning is one of the most critical components in the chain of radiation therapy of cancers. Treatment plans of today contain a wide variety of sophisticated information conveying the potential clinical effectiveness of the designed treatment to practitioners. Examples of such information include dose distributions superimposed on three- or even four-dimensional anatomic images; dose volume histograms, dose, dose-volume and dose-response indices for anatomic structures of interest; etc. These data are used for evaluating treatment plans and for making treatment decisions. The current state-of-the-art has evolved from the 1940s era when the dose to the tumor and normal tissues was estimated approximately by manual means. However, the symposium will cover the history of the field from the late-1950's, when computers were first introduced for treatment planning, to the present state involving the use of high performance computing and advanced multi-dimensional anatomic, functional and biological imaging, focusing only on external beam treatment planning. The symposium will start with a general overview of the treatment planning process including imaging, structure delineation, assignment of dose requirements, consideration of uncertainties, selection of beam configurations and shaping of beams, and calculations, optimization and evaluation of dose distributions. This will be followed by three presentations covering the evolution of treatment planning, which parallels the evolution of computers, availability of advanced volumetric imaging and the development of novel technologies such as dynamic multi-leaf collimators and online image guidance. This evolution will be divided over three distinct periods - prior to 1970's, the 2D era; from 1980 to the mid-1990's, the 3D era; and from the mid 1990's to today, the IMRT era. When the World was Flat: The Two-Dimensional Radiation Therapy Era” - Jacob Van Dyk In the 2D era, anatomy was defined with the aid of solder wires, special contouring devices and projection x-rays. Dose distributions were calculated manually from single field, flat surface isodoses on transparencies. Precalculated atlases of generic dose distributions were produced by the International Atomic Energy Agency. Massive time-shared main frames and mini-computers were used to compute doses at individual points or dose distributions in a single plane. Beam shapes were generally rectangular, with wedges, missing tissue compensators and occasional blocks to shield critical structures. Dose calculations were measurement-based or they used primary and scatter calculations based on scatter-air ratio methodologies. Dose distributions were displayed on line printers as alpha-numeric character maps or isodose patterns made with pen plotters. More than Pretty Pictures: 3D Treatment Planning and Conformal Therapy - Benedick A. Fraass The introduction of computed tomography allowed the delineation of anatomy three-dimensionally and, supported partly by contracts from the National Cancer Institute, made possible the introduction and clinical use of 3D treatment planning, leading to development and use of 3D conformal therapy in the 1980's. 3D computer graphics and 3D anatomical structure definitions made possible Beam's Eye View (BEV) displays, making conformal beam shaping and much more sophisticated beam arrangements possible. These conformal plans significantly improved target dose coverage as well as normal tissue sparing. The use of dose volume histograms, gross/clinical/planning target volumes, MRI and PET imaging, multileaf collimators, and computer-controlled treatment delivery made sophisticated planning approaches practical. The significant improvements in dose distributions and analysis achievable with 3D conformal therapy made possible formal dose escalation and normal tissue tolerance clinical studies that set new and improved expectations for improved local control and decreasing complications in many clinical sites. From the Art to the State of the Art: Inverse Planning and IMRT - Thomas R. Bortfeld While the potential of intensity modulation was recognized in the mid- 1980's, intensity-modulated radiotherapy (IMRT) did not become a reality until the mid-1990's. Broad beams of photons could be sub-divided into narrow beamlets whose intensities could be determined using sophisticated optimization algorithms to appropriately balance tumor dose with normal tissue sparing. The development of dynamic multi-leaf collimators (on conventional linear accelerators as well as in helical delivery devices) enabled the efficient delivery of IMRT. The evolution of IMRT planning is continuing in the form of Volumetric Modulated Arc Therapy (VMAT) and through advanced optimization tools, such as multi-criteria optimization, automated IMRT planning, and robust optimization to protect dose distributions against uncertainties. IMRT also facilitates “dose painting” in which different sub-volumes of the target are prescribed different doses. Clearly, these advancements are being made possible by the increasing power and lower cost of computers and developments in other fields such as imaging and operations research. Summary - Radhe Mohan The history does not end here. The advancement of treatment planning is expected to continue, leading to further automation and improvements in conformality and robustness of dose distributions, particularly in the area of particle therapy. Radiobiological modeling will gain emphasis as part of the planning process. Learning Objectives: The scope of changes in technology and the capabilities of radiation treatment planning The impact of these changes in the quality of treatment plans and optimality of dose distributions The impact of development in other fields (imaging, computers, operations research, etc.) on the evolution of radiation treatment planning.

60 APPLIED LIFE SCIENCES↗

Effect of dosimeter type for commissioning small photon beams on calculated dose distribution in stereotactic radiosurgery

Purpose: To assess the impact of the detector used to commission small photon beams on the calculated dose distribution in stereotactic radiosurgery (SRS). Methods: In this study, six types of detectors were used to characterize small photon beams: three diodes [a silicon stereotactic field diode SFD, a silicon diode SRS, and a silicon diode E], an ionization chamber CC01, and two types of radiochromic film models EBT and EBT2. These detectors were used to characterize circular collimated beams that were generated by a Novalis linear accelerator. This study was conducted in two parts. First, the following dosimetric data, which are of particular interest in SRS, were compared for the different detectors: the total scatter factor (TSF), the tissue phantom ratios (TPRs), and the off-axis ratios (OARs). Second, the commissioned data sets were incorporated into the treatment planning system (TPS) to compare the calculated dose distributions and the dose volume histograms (DVHs) that were obtained using the different detectors. Results: The TSFs data measured by all of the detectors were in good agreement with each other within the respective statistical uncertainties: two exceptions, where the data were systematically below those obtained for the other detectors, were the CC01 results for all of the circular collimators and the EBT2 film results for circular collimators with diameters below 10.0 mm. The OAR results obtained for all of the detectors were in excellent agreement for all of the circular collimators. This observation was supported by the gamma-index test. The largest difference in the TPR data was found for the 4.0 mm circular collimator, followed by the 10.0 and 20.0 mm circular collimators. The results for the calculated dose distributions showed that all of the detectors passed the gamma-index test at 100% for the 3 mm/3% criteria. The aforementioned observation was true regardless of the size of the calculation grid for all of the circular collimators. Finally, the dose volume histogram results were independent of the size of the calculation grid used. Conclusions: The results of this study showed that all of the studied detectors produced similar commissioned data sets for the TPS dose calculations. However, this result only validated the dose distribution calculation in the TPS and could not be used to assess the dose delivery to the target in which the TFS data were used to calculate the monitor units (the TFS data were not used in the TPS dose distribution calculation). Therefore, this study could not be used to determine the most accurate detector commissioning data set; however, all of the detectors exhibited superior performance for the relative dosimetry of small photon beams.

07 ISOTOPE AND RADIATION SOURCES↗

Assessing and accounting for the impact of respiratory motion on FDG uptake and viable volume for liver lesions in free-breathing PET using respiration-suspended PET images as reference

Purpose: To assess and account for the impact of respiratory motion on the variability of activity and volume determination of liver tumor in positron emission tomography (PET) through a comparison between free-breathing (FB) and respiration-suspended (RS) PET images. Methods: As part of a PET/computed tomography (CT) guided percutaneous liver ablation procedure performed on a PET/CT scanner, a patient's breathing is suspended on a ventilator, allowing the acquisition of a near-motionless PET and CT reference images of the liver. In this study, baseline RS and FB PET/CT images of 20 patients undergoing thermal ablation were acquired. The RS PET provides near-motionless reference in a human study, and thereby allows a quantitative evaluation of the effect of respiratory motion on PET images obtained under FB conditions. Two methods were applied to calculate tumor activity and volume: (1) threshold-based segmentation (TBS), estimating the total lesion glycolysis (TLG) and the segmented volume and (2) histogram-based estimation (HBE), yielding the background-subtracted lesion (BSL) activity and associated volume. The TBS method employs 50% of the maximum standardized uptake value (SUV{sub max}) as the threshold for tumors with SUV{sub max} ≥ 2× SUV{sub liver-bkg}, and tumor activity above this threshold yields TLG{sub 50%}. The HBE method determines local PET background based on a Gaussian fit of the low SUV peak in a SUV-volume histogram, which is generated within a user-defined and optimized volume of interest containing both local background and lesion uptakes. Voxels with PET intensity above the fitted background were considered to have originated from the tumor and used to calculate the BSL activity and its associated lesion volume. Results: Respiratory motion caused SUV{sub max} to decrease from RS to FB by −15% ± 11% (p = 0.01). Using TBS method, there was also a decrease in SUV{sub mean} (−18% ± 9%, p = 0.01), but an increase in TLG{sub 50%} (18% ± 36%) and in the segmented volume (47% ± 52%, p = 0.01) from RS to FB PET images. The background uptake in normal liver was stable, 1% ± 9%. In contrast, using the HBE method, the differences in both BSL activity and BSL volume from RS to FB were −8% ± 10% (p = 0.005) and 0% ± 16% (p = 0.94), respectively. Conclusions: This is the first time that almost motion-free PET images of the human liver were acquired and compared to free-breathing PET. The BSL method's results are more consistent, for the calculation of both tumor activity and volume in RS and FB PET images, than those using conventional TBS. This suggests that the BSL method might be less sensitive to motion blurring and provides an improved estimation of tumor activity and volume in the presence of respiratory motion.

60 APPLIED LIFE SCIENCES↗

A novel approach for evaluation of prostate deformation and associated dosimetric implications in IGRT of the prostate

Purpose: Prostate deformation is assumed to be a secondary correction and is typically ignored in the planning target volume (PTV) margin calculations. This assumption needs to be tested, especially when planning margins are reduced with daily image-guidance. In this study, deformation characteristics of the prostate and seminal vesicles were determined, and the dosimetric impact on treatment plans with different PTV margins was investigated. Methods: Ten prostate cancer patients were retrospectively selected for the study, each with three fiducial markers implanted in the prostate. Two hundred CBCT images were registered to respective planning CT images using a B-spline-based deformable image registration (DIR) software. A manual bony anatomy-based match was first applied based on the alignment of the pelvic bones and fiducial landmarks. DIR was then performed. For each registration, deformation vector fields (DVFs) of the prostate and seminal vesicles (SVs) were quantified using deformation-volume histograms. In addition, prostate rotation was evaluated and compared with prostate deformation. For a patient demonstrating small and large prostate deformations, target coverage degradation was analyzed in each of three treatment plans with PTV margins of 10 mm (6 mm at the prostate/rectum interface), as well as 5, and 3 mm uniformly. Results: Deformation of the prostate was most significant in the anterior direction. Maximum prostate deformation of greater than 10, 5, and 3 mm occurred in 1%, 17%, and 76% of the cases, respectively. Based on DVF-histograms, DVF magnitudes greater than 5 and 3 mm occurred in 2% and 27% of the cases, respectively. Deformation of the SVs was most significant in the posterior direction, and it was greater than 5 and 3 mm in 7.5% and 44.9% of the cases, respectively. Prostate deformation was found to be poorly correlated with rotation. Fifty percent of the cases showed rotation with negligible deformation and 7% of the cases showed significant deformation with minimal rotation (<3°). Average differences in the D{sub 95} dose to the prostate + SVs between the planning CT and CBCT images was 0.4% ± 0.5%, 3.0% ± 2.8%, and 6.6% ± 6.1%, respectively, for the plans with 10/6, 5, and 3 mm margins. For the case with both a large degree of prostate deformation (≈10% of the prostate volume) and rotation (≈8°), D{sub 95} was reduced by 0.5% ± 0.1%, 6.8% ± 0.6%, and 20.9% ± 1.6% for 10/6, 5, and 3 mm margin plans, respectively. For the case with large prostate deformation but negligible rotation (<1°), D{sub 95} was reduced by 0.4 ± 0.3, 3.9 ± 1.0, and 11.5 ± 2.5 for 10/6, 5, and 3 mm margin plans, respectively. Conclusions: Prostate deformation over a course of fractionated prostate radiotherapy may not be insignificant and may need to be accounted for in the planning margin design. A consequence of these results is that use of highly reduced planning margins must be viewed with caution.

60 APPLIED LIFE SCIENCES↗

Poster — Thur Eve — 69: Computational Study of DVH-guided Cancer Treatment Planning Optimization Methods

Purpose: To compare methods to incorporate the Dose Volume Histogram (DVH) curves into the treatment planning optimization. Method: The performance of three methods, namely, the conventional Mixed Integer Programming (MIP) model, a convex moment-based constrained optimization approach, and an unconstrained convex moment-based penalty approach, is compared using anonymized data of a prostate cancer patient. Three plans we generated using the corresponding optimization models. Four Organs at Risk (OARs) and one Tumor were involved in the treatment planning. The OARs and Tumor were discretized into total of 50,221 voxels. The number of beamlets was 943. We used commercially available optimization software Gurobi and Matlab to solve the models. Plan comparison was done by recording the model runtime followed by visual inspection of the resulting dose volume histograms. Conclusion: We demonstrate the effectiveness of the moment-based approaches to replicate the set of prescribed DVH curves. The unconstrained convex moment-based penalty approach is concluded to have the greatest potential to reduce the computational effort and holds a promise of substantial computational speed up.

60 APPLIED LIFE SCIENCES↗

Toward adaptive radiotherapy for head and neck patients: Uncertainties in dose warping due to the choice of deformable registration algorithm

Purpose: The aims of this work were to evaluate the performance of several deformable image registration (DIR) algorithms implemented in our in-house software (NiftyReg) and the uncertainties inherent to using different algorithms for dose warping. Methods: The authors describe a DIR based adaptive radiotherapy workflow, using CT and cone-beam CT (CBCT) imaging. The transformations that mapped the anatomy between the two time points were obtained using four different DIR approaches available in NiftyReg. These included a standard unidirectional algorithm and more sophisticated bidirectional ones that encourage or ensure inverse consistency. The forward (CT-to-CBCT) deformation vector fields (DVFs) were used to propagate the CT Hounsfield units and structures to the daily geometry for “dose of the day” calculations, while the backward (CBCT-to-CT) DVFs were used to remap the dose of the day onto the planning CT (pCT). Data from five head and neck patients were used to evaluate the performance of each implementation based on geometrical matching, physical properties of the DVFs, and similarity between warped dose distributions. Geometrical matching was verified in terms of dice similarity coefficient (DSC), distance transform, false positives, and false negatives. The physical properties of the DVFs were assessed calculating the harmonic energy, determinant of the Jacobian, and inverse consistency error of the transformations. Dose distributions were displayed on the pCT dose space and compared using dose difference (DD), distance to dose difference, and dose volume histograms. Results: All the DIR algorithms gave similar results in terms of geometrical matching, with an average DSC of 0.85 ± 0.08, but the underlying properties of the DVFs varied in terms of smoothness and inverse consistency. When comparing the doses warped by different algorithms, we found a root mean square DD of 1.9% ± 0.8% of the prescribed dose (pD) and that an average of 9% ± 4% of voxels within the treated volume failed a 2%pD DD-test (DD{sub 2%-pp}). Larger DD{sub 2%-pp} was found within the high dose gradient (21% ± 6%) and regions where the CBCT quality was poorer (28% ± 9%). The differences when estimating the mean and maximum dose delivered to organs-at-risk were up to 2.0%pD and 2.8%pD, respectively. Conclusions: The authors evaluated several DIR algorithms for CT-to-CBCT registrations. In spite of all methods resulting in comparable geometrical matching, the choice of DIR implementation leads to uncertainties in dose warped, particularly in regions of high gradient and/or poor imaging quality.

60 APPLIED LIFE SCIENCES↗