Choonik Lee

Hybrid computational phantoms of the male and female newborn patient : NURBS-based whole-body models

Anthropomorphic computational phantoms are computer models of the human body for use in the evaluation of dose distributions resulting from either internal or external radiation sources. Currently, two classes of computational phantoms have been developed and widely utilized for organ dose assessment: (1) stylized phantoms and (2) voxel phantoms which describe the human anatomy via mathematical surface equations or 3D voxel matrices, respectively. Although stylized phantoms based on mathematical equations can be very flexible in regard to making changes in organ position and geometrical shape, they are limited in their ability to fully capture the anatomic complexities of human internal anatomy. In turn, voxel phantoms have been developed through image-based segmentation and correspondingly provide much better anatomical realism in comparison to simpler stylized phantoms. However, they themselves are limited in defining organs presented in low contrast within either magnetic resonance or computed tomography images-the two major sources in voxel phantom construction. By definition, voxel phantoms are typically constructed via segmentation of transaxial images, and thus while fine anatomic features are seen in this viewing plane, slice-to-slice discontinuities become apparent in viewing the anatomy of voxel phantoms in the sagittal or coronal planes. This study introduces the concept of a hybrid computational newborn phantom that takes full advantage of the best features of both its stylized and voxel counterparts: flexibility in phantom alterations and anatomic realism. Non-uniform rational B-spline (NURBS) surfaces, a mathematical modeling tool traditionally applied to graphical animation studies, was adopted to replace the limited mathematical surface equations of stylized phantoms. A previously developed whole-body voxel phantom of the newborn female was utilized as a realistic anatomical framework for hybrid phantom construction. The construction of a hybrid phantom is performed in three steps: polygonization of the voxel phantom, organ modeling via NURBS surfaces and phantom voxelization. Two 3D graphic tools, 3D-DOCTOR and Rhinoceros, were utilized to polygonize the newborn voxel phantom and generate NURBS surfaces, while an in-house MATLAB code was used to voxelize the resulting NURBS model into a final computational phantom ready for use in Monte Carlo radiation transport calculations. A total of 126 anatomical organ and tissue models, including 38 skeletal sites and 31 cartilage sites, were described within the hybrid phantom using either NURBS or polygon surfaces. A male hybrid newborn phantom was constructed following the development of the female phantom through the replacement of female-specific organs with male-specific organs. The outer body contour and internal anatomy of the NURBS-based phantoms were adjusted to match anthropometric and reference newborn data reported by the International Commission on Radiological Protection in their Publication 89. The voxelization process was designed to accurately convert NURBS models to a voxel phantom with minimum volumetric change. A sensitivity study was additionally performed to better understand how the meshing tolerance and voxel resolution would affect volumetric changes between the hybrid-NURBS and hybrid-voxel phantoms. The male and female hybrid-NURBS phantoms were constructed in a manner so that all internal organs approached their ICRP reference masses to within 1%, with the exception of the skin (-6.5% relative error) and brain (-15.4% relative error). Both hybrid-voxel phantoms were constructed with an isotropic voxel resolution of 0.663 mm--equivalent to the ICRP 89 reference thickness of the newborn skin (dermis and epidermis). Hybrid-NURBS phantoms used to create their voxel counterpart retain the non-uniform scalability of stylized phantoms, while maintaining the anatomic realism of segmented voxel phantoms with respect to organ shape, depth and inter-organ positioning.

Assessment of parotid gland dose changes during head and neck cancer radiotherapy using daily megavoltage computed tomography and deformable image registration.

PURPOSE To analyze changes in parotid gland dose resulting from anatomic changes throughout a course of radiotherapy in a cohort of head-and-neck cancer patients. METHODS AND MATERIALS The study population consisted of 10 head-and-neck cancer patients treated definitively with intensity-modulated radiotherapy on a helical tomotherapy unit. A total of 330 daily megavoltage computed tomography images were retrospectively processed through a deformable image registration algorithm to be registered to the planning kilovoltage computed tomography images. The process resulted in deformed parotid contours and voxel mappings for both daily and accumulated dose-volume histogram calculations. The daily and cumulative dose deviations from the original treatment plan were analyzed. Correlations between dosimetric variations and anatomic changes were investigated. RESULTS The daily parotid mean dose of the 10 patients differed from the plan dose by an average of 15%. At the end of the treatment, 3 of the 10 patients were estimated to have received a greater than 10% higher mean parotid dose than in the original plan (range, 13-42%), whereas the remaining 7 patients received doses that differed by less than 10% (range, -6-8%). The dose difference was correlated with a migration of the parotids toward the high-dose region. CONCLUSIONS The use of deformable image registration techniques and daily megavoltage computed tomography imaging makes it possible to calculate daily and accumulated dose-volume histograms. Significant dose variations were observed as result of interfractional anatomic changes. These techniques enable the implementation of dose-adaptive radiotherapy.

Evaluation of geometric changes of parotid glands during head and neck cancer radiotherapy using daily MVCT and automatic deformable registration

BACKGROUND AND PURPOSE To assess and evaluate geometrical changes in parotid glands using deformable image registration and megavoltage CT (MVCT) images. METHODS A deformable registration algorithm was applied to 330 daily MVCT images (10 patients) to create deformed parotid contours. The accuracy and robustness of the algorithm was evaluated through visual review, comparison with manual contours, and precision analysis. Temporal changes in the parotid gland geometry were observed. RESULTS The deformed parotid contours were qualitatively judged to be acceptable. Compared with manual contours, the uncertainties of automatically deformed contours were similar with regard to geometry and dosimetric endpoint. The day-to-day variations (1 standard deviation of errors) in the center-of-mass distance and volume were 1.61mm and 4.36%, respectively. The volumes tended to decrease with a median total loss of 21.3% (6.7-31.5%) and a median change rate of 0.7%/day (0.4-1.3%/day). Parotids migrated toward the patient center with a median total distance change of -5.26mm (0.00 to -16.35mm) and a median change rate of -0.22mm/day (0.02 to -0.56mm/day). CONCLUSION The deformable image registration and daily MVCT images provide an efficient and reliable assessment of parotid changes over the course of a radiation therapy.

State of dose prescription and compliance to international standard (ICRU-83) in intensity modulated radiation therapy among academic institutions

PURPOSE The purpose of this study was to evaluate dose prescription and recording compliance to international standard (International Commission on Radiation Units & Measurements [ICRU]-83) in patients treated with intensity modulated radiation therapy (IMRT) among academic institutions. METHODS AND MATERIALS Ten institutions participated in this study to collect IMRT data to evaluate compliance to ICRU-83. Under institutional review board clearance, data from 5094 patients-including treatment site, technique, planner, physician, prescribed dose, target volume, monitor units, planning system, and dose calculation algorithm-were collected anonymously. The dose-volume histogram of each patient, as well as dose points, doses delivered to 100% (D100), 98% (D98), 95% (D95), 50% (D50), and 2% (D2), of sites was collected and sent to a central location for analysis. Homogeneity index (HI) as a measure of the steepness of target and is a measure of the shape of the dose-volume histogram was calculated for every patient and analyzed. RESULTS In general, ICRU recommendations for naming the target, reporting dose prescription, and achieving desired levels of dose to target were relatively poor. The nomenclature for the target in the dose prescription had large variations, having every permutation of name and number contrary to ICRU recommendations. There was statistically significant variability in D95, D50, and HI among institutions, tumor site, and technique with P values < .01. Nearly 95% of patients had D50 higher than 100% (103.5 ± 6.9) of prescribed dose and varied among institutions. On the other hand, D95 was close to 100% (97.1 ± 9.4) of prescribed dose. Liver and lung sites had a higher D50 compared with other sites. Pelvic sites had a lower variability indicated by HI (0.13 ± 1.21). Variability in D50 is 101.2 ± 8.5, 103.4 ± 6.8, 103.4 ± 8.2, and 109.5 ± 11.5 for IMRT, tomotherapy, volume modulated arc therapy, and stereotactic body radiation therapy with IMRT, respectively. CONCLUSIONS Nearly 95% of patient treatments deviated from the ICRU-83 recommended D50 prescription dose delivery. This variability is significant (P < .01) in terms of treatment site, technique, and institution. To reduce dosimetric and associated radiation outcome variability, dose prescription in every clinical trial should be unified with international guidelines.

Single or multi-channel vaginal cuff high-dose-rate brachytherapy: Is replanning necessary prior to each fraction?

PURPOSE Adjuvant high-dose-rate vaginal brachytherapy (VB) is commonly used in endometrial cancer. We evaluated the dosimetric and cost differences of using either a single plan or replan prior to each fraction for single- and multi-channel VB. METHODS AND MATERIALS We evaluated 84 fractions from 25 patients at our institution (16 single-channel patients each 3 fractions; 9 multi-channel patients each 4 fractions). All fractions were preceded by a computed tomographic (CT) simulation scan, after which a unique treatment plan was generated, dose points per International Commission on Radiation Units and Measurements (ICRU) 38. We calculated the dose to critical organs based on a decay-and-treat method utilizing the original catheter dwell-times for the initial fraction, and also the interfractional motion of the critical organ points between the initial and the subsequent CT scans. RESULTS The absolute mean dose difference was 14 cGy for bladder and 15 cGy for rectum between the replan and decay methods for single-channel, and 14 cGy for both organ points for the multi-channel cylinder. The bladder and rectum doses were not found to be significantly different between the replan and decay methods for either single-channel (bladder, P = .08; rectum, P = .19) or multi-channel cylinders (bladder, P = .85; rectum, P = .10). The mean interfractional displacement of the organ points between the initial and subsequent CT scans was 1.10 cm for the bladder and 0.67 cm for the rectum for single-channel, and 0.87 cm and 0.51 cm for multi-channel cylinders. The maximum interfractional motion was seen in the transverse plane for both organ points for both types of cylinders. At our institution, the decay method was 19% and 22% more cost-effective for single-channel and multi-channel cylinders, respectively. CONCLUSIONS Our data show no dosimetric advantage, but higher costs, associated with replanning prior to each fraction for both single- and multi-channel VB. Fractional replanning should not be utilized on a routine basis.

Improving treatment plan evaluation with automation

The goal of this work is to evaluate the effectiveness of Plan‐Checker Tool (PCT) which was created to improve first‐time plan quality, reduce patient delays, increase the efficiency of our electronic workflow, and standardize and automate the physics plan review in the treatment planning system (TPS). PCT uses an application programming interface to check and compare data from the TPS and treatment management system (TMS). PCT includes a comprehensive checklist of automated and manual checks that are documented when performed by the user as part of a plan readiness check for treatment. Prior to and during PCT development, errors identified during the physics review and causes of patient treatment start delays were tracked to prioritize which checks should be automated. Nineteen of 33 checklist items were automated, with data extracted with PCT. There was a 60% reduction in the number of patient delays in the six months after PCT release. PCT was successfully implemented for use on all external beam treatment plans in our clinic. While the number of errors found during the physics check did not decrease, automation of checks increased visibility of errors during the physics check, which led to decreased patient delays. The methods used here can be applied to any TMS and TPS that allows queries of the database. PACS number(s): 87.55.‐x, 87.55.N‐, 87.55.Qr, 87.55.tm, 89.20.BbThe goal of this work is to evaluate the effectiveness of Plan-Checker Tool (PCT) which was created to improve first-time plan quality, reduce patient delays, increase the efficiency of our electronic workflow, and standardize and automate the physics plan review in the treatment planning system (TPS). PCT uses an application programming interface to check and compare data from the TPS and treatment management system (TMS). PCT includes a comprehensive checklist of automated and manual checks that are documented when performed by the user as part of a plan readiness check for treatment. Prior to and during PCT development, errors identified during the physics review and causes of patient treatment start delays were tracked to prioritize which checks should be automated. Nineteen of 33 checklist items were automated, with data extracted with PCT. There was a 60% reduction in the number of patient delays in the six months after PCT release. PCT was successfully implemented for use on all external beam treatment plans in our clinic. While the number of errors found during the physics check did not decrease, automation of checks increased visibility of errors during the physics check, which led to decreased patient delays. The methods used here can be applied to any TMS and TPS that allows queries of the database. PACS number(s): 87.55.-x, 87.55.N-, 87.55.Qr, 87.55.tm, 89.20.Bb.

A Comparison of Soft-Tissue Implanted Markers and Bony Anatomy Alignments for Image-Guided Treatments of Head-and-Neck Cancers

PURPOSE To compare the geometric alignments of soft-tissue implanted markers to the traditional bony-based alignments in head-and-neck cancers, on the basis of daily image guidance. Dosimetric impact of the two alignment techniques on target coverage is presented. METHODS AND MATERIALS A total of 330 retrospective alignments (5 patients) were performed on daily megavoltage computed tomography (MVCT) image sets using both alignment techniques. Intermarker distances were tracked for all fractions to assess marker interfractional stability. Using a deformable image registration algorithm, target cumulative doses were calculated according to generated shifts on daily MVCT image sets. Target D95 was used as a dosimetric endpoint to evaluate each alignment technique. RESULTS Intermarker distances overall were stable, with a standard deviation of <1.5 mm for all fractions and no observed temporal trends. Differences in shift magnitudes between both alignment techniques were found to be statistically significant, with a maximum observed difference of 8 mm in a given direction. Evaluation of technique-specific dose coverage based on D95 of target clinical target volume and planning target volume shows small differences (within +/-5%) compared with the kilovoltage CT plan. CONCLUSION The use of daily MVCT imaging demonstrates that implanted markers in oral tongue and soft-palate cancers are stable localization surrogates. Alignments based on implanted markers generate shifts comparable overall to the traditional bony-based alignment, with no observed systematic difference in magnitude or direction. The cumulative dosimetric impact on target clinical target volume and planning target volume coverage was found to be similar, despite large observed differences in daily alignment shifts between the two techniques.

Dosimetry calculations for internal electron sources using a Korean reference adult stylised phantom

Absorbed fractions (AFs) and specific absorbed fractions (SAFs) for internally deposited electron were calculated using a Korean reference adult stylised phantom, where a total of 15 internal organ volumes and external body dimension were designed to match average Korean adult male. The walls of oesophagus, stomach, colon and urinary bladder were additionally divided into the mucosal layer and residual wall to accommodate dose calculation for weakly penetrating electron. The mucosal wall thicknesses were determined by the data reported in the International Commission on Radiological Protection Publication 89 and other literature resources and by direct measurements. The Monte Carlo transport code MCNPX (version 2.5.0) was employed to calculate the electron energy deposited. The SAFs and AFs for monoenergetic electrons with the energies ranging from 10 keV to 2 MeV were calculated. The results were compared with those of the revised Oak Ridge National Laboratory phantoms and showed considerable differences up to 150% in SAFs, whereas no substantial differences were observed in the AFs.

Revisiting fetal dose during radiation therapy: evaluating treatment techniques and a custom shield

To create a comprehensive dataset of peripheral dose (PD) measurements from a new generation of linear accelerators with and without the presence of a newly designed fetal shield, PD measurements were performed to evaluate the effects of depth, field size, distance from the field edge, collimator angle, and beam modifiers for common treatment protocols and modalities. A custom fetal lead shield was designed and made for our department that allows external beam treatments from multiple angles while minimizing the need to adjust the shield during patient treatments. PD measurements were acquired for a comprehensive series of static fields on stack of Solid Water. Additionally, PDs from various clinically relevant treatment scenarios for pregnant patients were measured using an anthropomorphic phantom that was abutted to a stack of Solid Water. As expected, the PD decreased as the distance from the field edge increased and the field size decreased. On average, a PD reduction was observed when a 90° collimator rotation was applied and/or when the tertiary MLCs and jaws defined the field aperture. However, the effect of the collimator rotation (90° versus 0°) in PD reduction was not found to be clinically significant when the tertiary MLCs were used to define the field aperture. In the presence of both the MLCs and the fetal shield, the PD was reduced by 58% at a distance of 10 cm from the field edge. The newly designed fetal shield may effectively reduce fetal dose and is relatively easy to setup. Due to its design, we are able to use a broad range of treatment techniques and beam angles. We believe the acquired comprehensive PD dataset collected with and without the fetal shield will be useful for treatment teams to estimate fetal dose and help guide decisions on treatment techniques without the need to perform pretreatment phantom measurements. PACS numbers: 87.53.Bn, 87.55.D-, 87.55.N.To create a comprehensive dataset of peripheral dose (PD) measurements from a new generation of linear accelerators with and without the presence of a newly designed fetal shield, PD measurements were performed to evaluate the effects of depth, field size, distance from the field edge, collimator angle, and beam modifiers for common treatment protocols and modalities. A custom fetal lead shield was designed and made for our department that allows external beam treatments from multiple angles while minimizing the need to adjust the shield during patient treatments. PD measurements were acquired for a comprehensive series of static fields on a stack of Solid Water. Additionally, PDs from various clinically relevant treatment scenarios for pregnant patients were measured using an anthropomorphic phantom that was abutted to a stack of Solid Water. As expected, the PD decreased as the distance from the field edge increased and the field size decreased. On average, a PD reduction was observed when a 90° collimator rotation was applied and/or when the tertiary MLCs and jaws defined the field aperture. However, the effect of the collimator rotation (90° versus 0°) in PD reduction was not found to be clinically significant when the tertiary MLCs were used to define the field aperture. In the presence of both the MLCs and the fetal shield, the PD was reduced by 58% at a distance of 10 cm from the field edge. The newly designed fetal shield may effectively reduce fetal dose and is relatively easy to setup. Due to its design, we are able to use a broad range of treatment techniques and beam angles. We believe the acquired comprehensive PD dataset collected with and without the fetal shield will be useful for treatment teams to estimate fetal dose and help guide decisions on treatment techniques without the need to perform pretreatment phantom measurements. PACS number(s): 87.53.Bn, 87.55.D‐, 87.55.N

SafetyNet: Streamlining and Automating QA in radiotherapy

Proper quality assurance (QA) of the radiotherapy process can be time‐consuming and expensive. Many QA efforts, such as data export and import, are inefficient when done by humans. Additionally, humans can be unreliable, lose attention, and fail to complete critical steps that are required for smooth operations. In our group we have sought to break down the QA tasks into separate steps and to automate those steps that are better done by software running autonomously or at the instigation of a human. A team of medical physicists and software engineers worked together to identify opportunities to streamline and automate QA. Development efforts follow a formal cycle of writing software requirements, developing software, testing and commissioning. The clinical release process is separated into clinical evaluation testing, training, and finally clinical release. We have improved six processes related to QA and safety. Steps that were previously performed by humans have been automated or streamlined to increase first‐time quality, reduce time spent by humans doing low‐level tasks, and expedite QA tests. Much of the gains were had by automating data transfer, implementing computer‐based checking and automation of systems with an event‐driven framework. These coordinated efforts by software engineers and clinical physicists have resulted in speed improvements in expediting patient‐sensitive QA tests. PACS number(s): 87.55.Ne, 87.55.Qr, 87.55.tg, 87.55.tmProper quality assurance (QA) of the radiotherapy process can be time-consuming and expensive. Many QA efforts, such as data export and import, are inefficient when done by humans. Additionally, humans can be unreliable, lose attention, and fail to complete critical steps that are required for smooth operations. In our group we have sought to break down the QA tasks into separate steps and to automate those steps that are better done by software running autonomously or at the instigation of a human. A team of medical physicists and software engineers worked together to identify opportunities to streamline and automate QA. Development efforts follow a formal cycle of writing software requirements, developing software, testing and commissioning. The clinical release process is separated into clinical evaluation testing, training, and finally clinical release. We have improved six processes related to QA and safety. Steps that were previously performed by humans have been automated or streamlined to increase first-time quality, reduce time spent by humans doing low-level tasks, and expedite QA tests. Much of the gains were had by automating data transfer, implementing computer-based checking and automation of systems with an event-driven framework. These coordinated efforts by software engineers and clinical physicists have resulted in speed improvements in expediting patient-sensitive QA tests. PACS number(s): 87.55.Ne, 87.55.Qr, 87.55.tg, 87.55.tm.