Search

Background and purpose Magnetic resonance imaging (MRI)-only treatment planning is gaining in popularity in radiation oncology, with various methods available to generate a synthetic computed tomography (sCT) for this purpose. The aim of this study was to validate a sCT generation software for MRI-only radiotherapy planning of male and female pelvic cancers. The secondary aim of this study was to improve dose agreement by applying a derived relative electron and mass density (RED) curve to the sCT. Method and materials Computed tomography (CT) and MRI scans of forty patients with pelvic neoplasms were used in the study. Treatment plans were copied from the CT scan to the sCT scan for dose comparison. Dose difference at reference point, 3D gamma comparison and dose volume histogram analysis was used to validate the dose impact of the sCT. The RED values were optimised to improve dose agreement by using a linear plot. Results The average percentage dose difference at isocentre was 1.2% and the mean 3D gamma comparison with a criteria of 1%/1 mm was 84.0% ± 9.7%. The results indicate an inherent systematic difference in the dosimetry of the sCT plans, deriving from the tissue densities. With the adapted REDmod table, the average percentage dose difference was reduced to −0.1% and the mean 3D gamma analysis improved to 92.9% ± 5.7% at 1%/1 mm. Conclusions CT generation software is a viable solution for MRI-only radiotherapy planning. The option makes it relatively easy for departments to implement a MRI-only planning workflow for cancers of male and female pelvic anatomy.

The introduction of MRI linear accelerators (MR-linacs) and the increased use of MR imaging in radiotherapy, requires improved approaches to MRI-only radiotherapy. MRI provides excellent soft tissue visualisation but does not provide any electron density information required for radiotherapy dose calculation, instead MRI is registered to CT images to enable dose calculations. MRI-only radiotherapy eliminates registration errors and reduces patient discomfort, workload and cost. Electron density requirements may be addressed in different ways, from manually applying bulk density corrections, to more computationally intensive methods to produce substitute CT datasets (sCT), requiring additional sequences, increasing overall imaging time. Reducing MR imaging time would reduce potential artefacts from intrafraction motion and patient discomfort. The aim of this study was to assess the impact of reducing MR imaging time on a hybrid atlas-voxel sCT conversion for prostate MRI-only treatment planning, considering both anatomical and dosimetric parameters. 10 volunteers were scanned on a Siemens Skyra 3T MRI. Sequences included the 3D T2-weighted (T2-w) SPACE sequence used for sCT conversion as previously validated against CT, along with variations to this sequence in repetition time (TR), turbo factor, and combinations of these to reduce the imaging time. All scans were converted to sCT and were compared to the sCT from the original SPACE sequence, evaluating for anatomical changes and dosimetric differences for a standard prostate VMAT plan. Compared to the previously validated T2-w SPACE sequence, scan times were reduced by up to 80%. The external volume and bony anatomy were compared, with all but one sequence meeting a DICE coefficient of 0.9 or better, with the largest variations occurring at the edges of the external body volume. The generated sCT agreed with the gold standard sCT within an isocentre dose of 1% and a gamma pass rate of 99% for a 1%/1 mm gamma tolerance for all but one sequence. This study demonstrates that the MR imaging sequence time was able to be reduced by approximately 80% with similar dosimetric results.

Radiotherapy planning for lung cancer typically requires both 3D and 4D Computed Tomography (CT) to account for respiratory related movement. 4D Magnetic Resonance Imaging (MRI) with self-navigation offers a potential alternative with greater reliability in patients with irregular breathing patterns and improved soft tissue contrast. In this study 4D-CT and a 4D-MRI Radial Volumetric Interpolated Breath-hold Examination (VIBE) sequence was evaluated with a 4D phantom and 13 patient respiratory patterns, simulating tumour motion. Quantification of motion related tumour displacement in 4D-MRI and 4D-CT found no statistically significant difference in mean motion range. The results demonstrated the potential viability of 4D-MRI for lung cancer treatment planning.

PurposeStereotactic Ablative Radiotherapy (SABR) has recently emerged as a favourable treatment option for prostate cancer patients. With higher doses delivered over fewer fractions, motion adaptation is a requirement for accurate delivery of SABR. This study compared the efficacy of multileaf collimator (MLC) tracking vs. gating as a real-time motion adaptation strategy for prostate SABR patients enrolled in a clinical trial.MethodsForty-four prostate cancer patients treated over five fractions in the TROG 15.01 SPARK trial were analysed in this study. Forty-nine fractions were treated using MLC tracking and 166 fractions were treated using beam gating and couch shifts. A time-resolved motion-encoded dose reconstruction method was used to evaluate the dose delivered using each motion adaptation strategy and compared to an estimation of what would have been delivered with no motion adaptation strategy implemented.ResultsMLC tracking and gating both delivered doses closer to the plan compared to when no motion adaptation strategy was used. Differences between MLC tracking and gating were small with differences in the mean discrepancy from the plan of -0.3% (CTV D98%), 1.4% (CTV D2%), 0.4% (PTV D95%), 0.2% (rectum V30Gy) and 0.0% (bladder V30Gy). On average, 0.5 couch shifts were required per gated fractions with a mean interruption duration of 1.8 ± 2.6 min per fraction treated using gating.ConclusionBoth MLC tracking and gating were effective strategies at improving the accuracy of the dose delivered to the target and organs at risk. While dosimetric performance was comparable, gating resulted in interruptions to treatment.Clinical trial registration numberNCT02397317.

Purpose Deep learning methods (DLMs) have recently been proposed to generate pseudo-CT (pCT) for magnetic resonance imaging (MRI) based dose planning. This study aims to evaluate and compare DLMs (U-Net and generative adversarial network [GAN]) using various loss functions (L2, single-scale perceptual loss [PL], multiscale PL, weighted multiscale PL) and a patch-based method (PBM). Methods and Materials Thirty-nine patients received a volumetric modulated arc therapy for prostate cancer (78 Gy). T2-weighted MRIs were acquired in addition to planning CTs. The pCTs were generated from the MRIs using 7 configurations: 4 GANs (L2, single-scale PL, multiscale PL, weighted multiscale PL), 2 U-Net (L2 and single-scale PL), and the PBM. The imaging endpoints were mean absolute error and mean error, in Hounsfield units, between the reference CT (CTref) and the pCT. Dose uncertainties were quantified as mean absolute differences between the dose volume histograms (DVHs) calculated from the CTref and pCT obtained by each method. Three-dimensional gamma indexes were analyzed. Results Considering the image uncertainties in the whole pelvis, GAN L2 and U-Net L2 showed the lowest mean absolute error (≤34.4 Hounsfield units). The mean errors were not different than 0 (P ≤ .05). The PBM provided the highest uncertainties. Very few DVH points differed when comparing GAN L2 or U-Net L2 DVHs and CTref DVHs (P ≤ .05). Their dose uncertainties were ≤0.6% for the prostate planning target Volume V95%, ≤0.5% for the rectum V70Gy, and ≤0.1% for the bladder V50Gy. The PBM, U-Net PL, and GAN PL presented the highest systematic dose uncertainties. The gamma pass rates were >99% for all DLMs. The mean calculation time to generate 1 pCT was 15 s for the DLMs and 62 min for the PBM. Conclusions Generating pCT for MRI dose planning with DLMs and PBM provided low-dose uncertainties. In particular, the GAN L2 and U-Net L2 provided the lowest dose uncertainties together with a low computation time.

Purpose Methods have been recently developed to generate pseudo–computed tomography (pCT) for dose calculation in magnetic resonance imaging (MRI)–only radiation therapy. This study aimed to propose an original nonlocal mean patch-based method (PBM) and to compare this PBM to an atlas-based method (ABM) and to a bulk density method (BDM) for prostate MRI-only radiation therapy. Materials and Methods Thirty-nine patients received a volumetric modulated arc therapy for prostate cancer. In addition to the planning computed tomography (CT) scans, T2-weighted MRI scans were acquired. pCTs were generated from MRIs using 3 methods: an original nonlocal mean PBM, ABM, and BDM. The PBM was performed using feature extraction and approximate nearest neighbor search in a training cohort. The PBM accuracy was evaluated in a validation cohort by using imaging and dosimetric endpoints. Imaging endpoints included mean absolute error and mean error between Hounsfield units of the pCT and the reference CT (CTref). Dosimetric endpoints were based on dose-volume histograms calculated from the CTref and the pCTs for various volumes of interest and on 3-dimensional gamma analyses. The PBM uncertainties were compared with those of the ABM and BDM. Results The mean absolute error and mean error obtained from the PBM were 41.1 and –1.1 Hounsfield units. The PBM dose-volume histogram differences were 0.7% for prostate planning target volume V95%, 0.5% for rectum V70Gy, and 0.2% for bladder V50Gy. Compared with ABM and BDM, PBM provided significantly lower dose uncertainties for the prostate planning target volume (70-78 Gy), the rectum (8.5-29 Gy, 40-48 Gy, and 61-73 Gy), and the bladder (12-78 Gy). The PBM mean gamma pass rate (99.5%) was significantly higher than that of ABM (94.9%) or BDM (96.1%). Conclusions The proposed PBM provides low uncertainties with dose planned on CTref. These uncertainties were smaller than those of ABM and BDM and are unlikely to be clinically significant.

Magnetic resonance images (MRI) require intensity standardisation if they are used for the purpose of quantitative analysis as inherent variations in image intensity levels between different image sets are manifest due to technical factors. One approach is to standardise the image intensity values using a statistically applied biological reference tissue. The aim of this study is to compare the performance of differing candidate biological reference tissues for standardising T2WI intensity distributions. Fifty-one prostate cancer patients across two centres with different scanners were evaluated using the percentage interpatient coefficient of variation (%interCV) for four different biological references; femoral bone marrow, ischioanal fossa, obturator-internus muscle and bladder urine. The tissue with the highest reproducibility (lowest %interCV) in both centres was used for intensity standardisation of prostate T2WI using three different statistical measures (mean, Z-score, median + Interquartile Range). The performance of different standardisation methods was evaluated from the assessment of image intensity histograms and the percentage normalised root mean square error (%NRSME) of the healthy peripheral zone tissue. Ischioanal fossa as a reference tissue demonstrated the highest reproducibility with %interCV of 18.9 for centre1 and 11.2 for centre2. Using ischioanal fossa for statistical intensity standardisation and the median + Interquartile Range method demonstrated the lowest %NRMSE across centres for healthy peripheral zone tissues. This study demonstrates ischioanal fossa as a preferred reference tissue for standardising intensity values from T2WI of the prostate. Subsequent image standardisation using the median + Interquartile Range intensity of the reference tissue demonstrated a robust and reliable standardisation method for quantitative image assessment.

Introduction This study quantified inter-observer contouring variations for multiple male pelvic structures, many of which are of emerging relevance for prostate cancer radiotherapy progression and toxicity response studies. Methods Five prostate cancer patient datasets (CT and T2-weighted MR) were distributed to 13 observers for contouring. CT structures contoured included the clinical target volume (CTV), seminal vesicles, rectum, colon, bowel bag, bladder and peri-rectal space (PRS). MR contours included CTV, trigone, membranous urethra, penile bulb, neurovascular bundle and multiple pelvic floor muscles. Contouring variations were assessed using the intraclass correlation coefficient (ICC), Dice similarity coefficient (DSC), and multiple additional metrics. Results Clinical target volume (CT and MR), bladder, rectum and PRS contours showed excellent inter-observer agreement (median ICC = 0.97; 0.99; 1.00; 0.95; 0.90, DSC = 0.83 ± 0.05; 0.88 ± 0.05; 0.93 ± 0.03; 0.81 ± 0.07; 0.80 ± 0.06, respectively). Seminal vesicle contours were more variable (ICC = 0.75, DSC = 0.73 ± 0.14), while colon and bowel bag contoured volumes were consistent (ICC = 0.97; 0.97), but displayed poor overlap (DSC = 0.58 ± 0.22; 0.67 ± 0.21). Smaller MR structures showed significant inter-observer variations, with poor overlap for trigone, membranous urethra, penile bulb, and left and right neurovascular bundles (DSC = 0.44 ± 0.22; 0.41 ± 0.21; 0.66 ± 0.21; 0.16 ± 0.17; 0.15 ± 0.15). Pelvic floor muscles recorded moderate to strong inter-observer agreement (ICC = 0.50-0.97), although large outlier variations were observed. Conclusions Inter-observer contouring variation was significant for multiple pelvic structures contoured on MR.

Objectiu múltiple; SRS; Fantasma virtual Objetivo múltiple; SRS; Fantasma virtual Multiple‐target; SRS; Virtual phantom The aim of this study was to benchmark the accuracy of the VIrtual Phantom Epid dose Reconstruction (VIPER) software for pre-treatment dosimetric verification of multiple-target stereotactic radiosurgery (SRS). VIPER is an EPID-based method to reconstruct a 3D dose distribution in a virtual phantom from in-air portal images. Validation of the VIPER dose calculation was assessed using several MLC-defined fields for a 6 MV photon beam. Central axis percent depth doses (PDDs) and output factors were measured with an ionization chamber in a water tank, while dose planes at a depth of 10 cm in a solid flat phantom were acquired with radiochromic films. The accuracy of VIPER for multiple-target SRS plan verification was benchmarked against Monte Carlo simulations. Eighteen multiple-target SRS plans designed with the Eclipse treatment planning system were mapped to a cylindrical water phantom. For each plan, the 3D dose distribution reconstructed by VIPER within the phantom was compared with the Monte Carlo simulation, using a 3D gamma analysis. Dose differences (VIPER vs. measurements) generally within 2% were found for the MLC-defined fields, while film dosimetry revealed gamma passing rates (GPRs) ≥95% for a 3%/1 mm criteria. For the 18 multiple-target SRS plans, average 3D GPRs greater than 93% and 98% for the 3%/2 mm and 5%/2 mm criteria, respectively. Our results validate the use of VIPER as a dosimetric verification tool for pre-treatment QA of single-isocenter multiple-target SRS plans. The method requires no setup time on the linac and results in an accurate 3D characterization of the delivered dose.

Background Current multiparametric MRI (mp-MRI) in routine clinical practice has poor-to-moderate diagnostic performance for transition zone prostate cancer. The aim of this study was to evaluate the potential diagnostic performance of novel 1H magnetic resonance spectroscopic imaging (MRSI) using a semi-localized adiabatic selective refocusing (sLASER) sequence with gradient offset independent adiabaticity (GOIA) pulses in addition to the routine mp-MRI, including T2-weighted imaging (T2WI), diffusion-weighted imaging (DWI) and quantitative dynamic contrast enhancement (DCE) for transition zone prostate cancer detection, localization and grading. Methods Forty-one transition zone prostate cancer patients underwent mp-MRI with an external phased-array coil. Normal and cancer regions were delineated by two radiologists and divided into low-risk, intermediate-risk, and high-risk categories based on TRUS guided biopsy results. Support vector machine models were built using different clinically applicable combinations of T2WI, DWI, DCE, and MRSI. The diagnostic performance of each model in cancer detection was evaluated using the area under curve (AUC) of the receiver operating characteristic diagram. Then accuracy, sensitivity and specificity of each model were calculated. Furthermore, the correlation of mp-MRI parameters with low-risk, intermediate-risk and high-risk cancers were calculated using the Spearman correlation coefficient. Results The addition of MRSI to T2WI + DWI and T2WI + DWI + DCE improved the accuracy, sensitivity and specificity for cancer detection. The best performance was achieved with T2WI + DWI + MRSI where the addition of MRSI improved the AUC, accuracy, sensitivity and specificity from 0.86 to 0.99, 0.83 to 0.96, 0.80 to 0.95, and 0.85 to 0.97 respectively. The (choline + spermine + creatine)/citrate ratio of MRSI showed the highest correlation with cancer risk groups (r = 0.64, p Conclusion The inclusion of GOIA-sLASER MRSI into conventional mp-MRI significantly improves the diagnostic accuracy of the detection and aggressiveness assessment of transition zone prostate cancer.

Purpose. The purpose of this work was to report on the geometric uncertainty for patients treated with multi-leaf collimator (MLC) tracking for lung SABR to verify the accuracy of the system. Methods. Seventeen patients were treated as part of the MLC tracking for lung SABR clinical trial using electromagnetic beacons implanted around the tumor acting as a surrogate for target motion. Sources of uncertainties evaluated in the study included the surrogate-target positional uncertainty, the beam-surrogate tracking uncertainty, the surrogate localization uncertainty, and the target delineation uncertainty. Probability density functions (PDFs) for each source of uncertainty were constructed for the cohort and each patient. The total PDFs was computed using a convolution approach. The 95% confidence interval (CI) was used to quantify these uncertainties. Results. For the cohort, the surrogate-target positional uncertainty 95% CIs were ±2.5 mm (−2.0/3.0 mm) in left-right (LR), ±3.0 mm (−1.6/4.5 mm) in superior–inferior (SI) and ±2.0 mm (−1.8/2.1 mm) in anterior–posterior (AP). The beam-surrogate tracking uncertainty 95% CIs were ±2.1 mm (−2.1/2.1 mm) in LR, ±2.8 mm (−2.8/2.7 mm) in SI and ±2.1 mm (−2.1/2.0 mm) in AP directions. The surrogate localization uncertainty minimally impacted the total PDF with a width of ±0.6 mm. The target delineation uncertainty distribution 95% CIs were ±5.4 mm. For the total PDF, the 95% CIs were ±5.9 mm (−5.8/6.0 mm) in LR, ±6.7 mm (−5.8/7.5 mm) in SI and ±6.0 mm (−5.5/6.5 mm) in AP. Conclusion. This work reports the geometric uncertainty of MLC tracking for lung SABR by accounting for the main sources of uncertainties that occurred during treatment. The overall geometric uncertainty is within ±6.0 mm in LR and AP directions and ±6.7 mm in SI. The dominant uncertainty was the target delineation uncertainty. This geometric analysis helps put into context the range of uncertainties that may be expected during MLC tracking for lung SABR (ClinicalTrials.gov registration number: NCT02514512).

Background and purpose This study aimed to identify anatomically-localised regions where planned radiotherapy dose is associated with gastrointestinal toxicities in healthy tissues throughout the pelvic anatomy. Materials and methods Planned dose distributions for up to 657 patients of the Trans Tasman Radiation Oncology Group 03.04 RADAR trial were deformably registered onto a single exemplar computed tomography dataset. Voxel-based multiple comparison permutation dose difference testing, Cox regression modelling and LASSO feature selection were used to identify regions where dose-increase was associated with grade ≥2 rectal bleeding (RB) or tenesmus, according to the LENT/SOMA scale. This was externally validated by registering dose distributions from the RT01 (n = 388) and CHHiP (n = 241) trials onto the same exemplar and repeating the tests on each of these data sets, and on all three datasets combined. Results Voxel-based Cox regression and permutation dose difference testing revealed regions where increased dose was correlated with gastrointestinal toxicity. Grade ≥2 RB was associated with posteriorly extended lateral beams that manifested high doses (>55 Gy) in a small rectal volume adjacent to the clinical target volume. A correlation was found between grade ≥2 tenesmus and increased low-intermediate dose (∼25 Gy) at the posterior beam region, including the posterior rectum and perirectal fat space (PRFS). Conclusions The serial response of the rectum with respect to RB has been demonstrated in patients with posteriorly extended lateral beams. Similarly, the parallel response of the PRFS with respect to tenesmus has been demonstrated in patients treated with the posterior beam.

There is increasing interest in MR-only radiotherapy planning since it provides superb soft-tissue contrast without the registration uncertainties inherent in a CT-MR registration. However, MR images cannot readily provide the electron density information necessary for radiotherapy dose calculation. An algorithm which generates synthetic CTs for dose calculations from MR images of the prostate using an atlas of 3 T MR images has been previously reported by two of the authors. This paper aimed to evaluate this algorithm using MR data acquired at a different field strength and a different centre to the algorithm atlas. Twenty-one prostate patients received planning 1.5 T MR and CT scans with routine immobilisation devices on a flat-top couch set-up using external lasers. The MR receive coils were supported by a coil bridge. Synthetic CTs were generated from the planning MR images with ([Formula: see text]) and without (sCT) a one voxel body contour expansion included in the algorithm. This was to test whether this expansion was required for 1.5 T images. Both synthetic CTs were rigidly registered to the planning CT (pCT). A 6 MV volumetric modulated arc therapy plan was created on the pCT and recalculated on the sCT and [Formula: see text]. The synthetic CTs' dose distributions were compared to the dose distribution calculated on the pCT. The percentage dose difference at isocentre without the body contour expansion (sCT-pCT) was [Formula: see text] and with ([Formula: see text]-pCT) was [Formula: see text] (mean ± one standard deviation). The [Formula: see text] result was within one standard deviation of zero and agreed with the result reported previously using 3 T MR data. The sCT dose difference only agreed within two standard deviations. The mean ± one standard deviation gamma pass rate was [Formula: see text] for the sCT and [Formula: see text] for the [Formula: see text] (with [Formula: see text] global dose difference and [Formula: see text] distance to agreement gamma criteria). The one voxel body contour expansion improves the synthetic CT accuracy for MR images acquired at 1.5 T but requires the MR voxel size to be similar to the atlas MR voxel size. This study suggests that the atlas-based algorithm can be generalised to MR data acquired using a different field strength at a different centre.

In MR only radiation therapy planning, generation of the tissue specific HU map directly from the MRI would eliminate the need of CT image acquisition and may improve radiation therapy planning. The aim of this work is to generate and validate substitute CT (sCT) scans generated from standard T2 weighted MR pelvic scans in prostate radiation therapy dose planning. A Siemens Skyra 3T MRI scanner with laser bridge, flat couch and pelvic coil mounts was used to scan 39 patients scheduled for external beam radiation therapy for localized prostate cancer. For sCT generation a whole pelvis MRI (1.6 mm 3D isotropic T2w SPACE sequence) was acquired. Patients received a routine planning CT scan. Co-registered whole pelvis CT and T2w MRI pairs were used as training images. Advanced tissue specific non-linear regression models to predict HU for the fat, muscle, bladder and air were created from co-registered CT-MRI image pairs. On a test case T2w MRI, the bones and bladder were automatically segmented using a novel statistical shape and appearance model, while other soft tissues were separated using an Expectation-Maximization based clustering model. The CT bone in the training database that was most 'similar' to the segmented bone was then transformed with deformable registration to create the sCT component of the test case T2w MRI bone tissue. Predictions for the bone, air and soft tissue from the separate regression models were successively combined to generate a whole pelvis sCT. The change in monitor units between the sCT-based plans relative to the gold standard CT plan for the same IMRT dose plan was found to be [Formula: see text] (mean ± standard deviation) for 39 patients. The 3D Gamma pass rate was [Formula: see text] (2 mm/2%). The novel hybrid model is computationally efficient, generating an sCT in 20 min from standard T2w images for prostate cancer radiation therapy dose planning and DRR generation.

Background Kilovoltage Intrafraction Monitoring (KIM) is a method which determines the three-dimensional position of the prostate from two-dimensional kilovoltage (kV) projections taken during linac based radiotherapy treatment with real-time feedback. Rectal displacement devices (RDDs) allow for improved rectal dosimetry during prostate cancer treatment. This study used KIM to perform a preliminary investigation of prostate intrafraction motion observed in patients with an RDD in place. Methods Ten patients with intermediate to high-risk prostate cancer were treated with a Rectafix RDD in place during two boost fractions of 9.5–10 Gy delivered using volumetric modulated arc therapy (VMAT) on Clinac iX and Truebeam linacs. Two-dimensional kV projections were acquired during treatment. KIM software was used following treatment to determine the displacement of the prostate over time. The displacement results were analyzed to determine the percentage of treatment time the prostate spent within 1 mm, between 1 and 2 mm, between 2 and 3 mm and greater than 3 mm from its initial position. Results KIM successfully measured displacement for 19 prostate stereotactic boost fractions. The prostate was within 1 mm of its initial position for 84.8%, 1–2 mm for 14%, 2–3 mm 1.2% and ≥3 mm only 0.4% of the treatment time. Conclusions In this preliminary study using KIM, KIM was successfully used to measure prostate intrafraction motion, which was found to be small in the presence of a rectal displacement device. Trial registration The Hunter New England Human Research Ethics Committee reference number is 14/08/20/3.01.

Purpose The Trans-Tasman Radiation Oncology Group (TROG) 15.01 Stereotactic Prostate Adaptive Radiotherapy utilizing Kilovoltage intrafraction monitoring (SPARK) trial is a multicenter trial using Kilovoltage Intrafraction Monitoring (KIM) to monitor prostate position during the delivery of prostate radiation therapy. KIM increases the accuracy of prostate radiation therapy treatments and allows for hypofractionation. However, an additional imaging dose is delivered to the patient. A standardized procedure to determine the imaging dose per frame delivered using KIM was developed and applied at four radiation therapy centers on three different types of linear accelerator. Methods Dose per frame for kilovoltage imaging in fluoroscopy mode was measured in air at isocenter using an ion chamber. Beam quality and dose were determined for a Varian Clinac iX linear accelerator, a Varian Trilogy, four Varian Truebeams and one Elekta Synergy at four different radiation therapy centers. The imaging parameters used on the Varian machines were 125 kV, 80 mA, and 13 ms. The Elekta machine was measured at 120 kV, 80 mA, and 12 ms. Absorbed doses to the skin and the prostate for a typical SBRT prostate treatment length were estimated according to the IPEMB protocol. Results The average dose per kV frame to the skin was 0.24 ± 0.03 mGy. The average estimated absorbed dose to the prostate for all five treatment fractions across all machines measured was 39.9 ± 2.6 mGy for 1 Hz imaging, 199.7 ± 13.2 mGy for 5 Hz imaging and 439.3 ± 29.0 mGy for 11 Hz imaging. Conclusions All machines measured agreed to within 20%. Additional dose to the prostate from using KIM is at most 1.3% of the prescribed dose of 36.25 Gy in five fractions delivered during the trial.

Purpose We present the first clinical implementation of a real-time six-degree of freedom (6DoF) Kilovoltage Intrafraction Monitoring (KIM) system which tracks the cancer target translational and rotational motions during treatment. The method was applied to measure and correct for target motion during stereotactic body radiotherapy (SBRT) for prostate cancer. Methods Patient: A patient with prostate adenocarcinoma undergoing SBRT with 36.25 Gy, delivered in 5 fractions was enrolled in the study. 6DoF KIM technology: 2D positions of three implanted gold markers in each of the kV images (125 kV, 10 mA at 11 Hz) were acquired continuously during treatment. The 2D → 3D target position estimation was based on a probability distribution function. The 3D → 6DoF target rotation was calculated using an iterative closest point algorithm. The accuracy and precision of the KIM method was measured by comparing the real-time results with kV-MV triangulation. Results Of the five treatment fractions, KIM was utilised successfully in four fractions. The intrafraction prostate motion resulted in three couch shifts in two fractions when the prostate motion exceeded the pre-set action threshold of 2 mm for more than 5 s. KIM translational accuracy and precision were 0.3 ± 0.6 mm, −0.2 ± 0.3 mm and 0.2 ± 0.7 mm in the Left-Right (LR), Superior-Inferior (SI) and Anterior-Posterior (AP) directions, respectively. The KIM rotational accuracy and precision were 0.8° ± 2.0°, −0.5° ± 3.3° and 0.3° ± 1.6° in the roll, pitch and yaw directions, respectively. Conclusion This treatment represents, to the best of our knowledge, the first time a cancer patient’s tumour position and rotation have been monitored in real-time during treatment. The 6 DoF KIM system has sub-millimetre accuracy and precision in all three translational axes, and less than 1° accuracy and 4° precision in all three rotational axes.

Prostate cancer treatment planning can be performed using magnetic resonance imaging (MRI) only with sCT scans. However, sCT scans are computer generated from MRI data and therefore robust, efficient, and accurate patient-specific quality assurance methods for dosimetric verification are required. Bulk anatomical density (BAD) maps can be generated based on anatomical contours derived from the MRI image. This study investigates and optimizes the BAD map approach for sCT quality assurance with a large patient CT and MRI dataset. 3D T2-weighted MRI and full density CT images of 54 patients were used to create BAD maps with different tissue class combinations. Mean Hounsfield units (HU) of Fat (F: below −30 HU), the entire Tissue [T: excluding bone (B)], and Muscle (M: excluding bone and fat) were derived from the CT scans. CT based BAD maps (BADBT,CT and BADBMF,CT) and a conventional bone and water bulk-density method (BADBW,CT) were compared to full CT calculations with bone assignments to 366 HU (measured) and 288 HU (obtained from literature). Optimal bulk densities of Tissue for BADBT,CT and Bone for BADBMF,CT were derived to provide zero mean isocenter dose agreement to the CT plan. Using the optimal densities, the dose agreement of BADBT,CT and BADBMF,CT to CT was redetermined. These maps were then created for the MRI dataset using auto-generated contours and dose calculations compared to CT. The average mean density of Bone, Fat, Muscle, and Tissue were 365.5 ± 62.2, −109.5 ± 12.9, 23.3 ± 9.7, and −46.3 ± 15.2 HU, respectively. Comparing to other bulk-density maps, BADBMF,CT maps provided the closest dose to CT. Calculated optimal mean densities of Tissue and Bone were −32.7 and 323.7 HU, respectively. The isocenter dose agreement of the optimal density assigned BADBT,CT and BADBMF,CT to full density CT were 0.10 ± 0.65% and 0.01 ± 0.45%, respectively. The isocenter dose agreement of MRI generated BADBT,MR and BADBMF,MR to full density CT were −0.15 ± 0.90% and −0.16 ± 0.65%, respectively. The BAD method with optimal bulk densities can provide robust, accurate and efficient patient-specific quality assurance for dose calculations in MRI-only radiotherapy.

© 2019 American Association of Physicists in Medicine Purpose: Kilovoltage intrafraction monitoring (KIM) allows for real-time image guidance for tracking tumor motion in six-degrees-of-freedom (6DoF) on a standard linear accelerator. This study assessed the geometric accuracy and precision of KIM used to guide patient treatments in the TROG 15.01 multi-institutional Stereotactic Prostate Ablative Radiotherapy with KIM trial and investigated factors affecting accuracy and precision. Methods: Fractions from 44 patients with prostate cancer treated using KIM-guided SBRT were analyzed across four institutions, on two different linear accelerator models and two different beam models (6 MV and 10 MV FFF). The geometric accuracy and precision of KIM was assessed from over 33 000 images (translation) and over 9000 images (rotation) by comparing the real-time measured motion to retrospective kV/MV triangulation. Factors potentially affecting accuracy, including contrast-to-noise ratio (CNR) of kV images and incorrect marker segmentation, were also investigated. Results: The geometric accuracy and precision did not depend on treatment institution, beam model or motion magnitude, but was correlated with gantry angle. The centroid geometric accuracy and precision of the KIM system for SABR prostate treatments was 0.0 ± 0.5, 0.0 ± 0.4 and 0.1 ± 0.3 mm for translation, and −0.1 ± 0.6°, −0.1 ± 1.4° and −0.1 ± 1.0° for rotation in the AP, LR and SI directions respectively. Centroid geometric error exceeded 2 mm for 0.05% of this dataset. No significant relationship was found between large geometric error and CNR or marker segmentation correlation. Conclusions: This study demonstrated the ability of KIM to locate the prostate with accuracy below other uncertainties in radiotherapy treatments, and the feasibility for KIM to be implemented across multiple institutions.

BACKGROUND Due to the histological heterogeneity of the central gland, accurate detection of central gland prostate cancer remains a challenge. PURPOSE To evaluate the efficacy of in vivo 3D 1 H MR spectroscopic imaging (3D 1 H MRSI) with a semi-localized adiabatic selective refocusing (sLASER) sequence and gradient-modulated offset-independent adiabatic (GOIA) pulses for detection of central gland prostate cancer. Additionally four risk models were developed to differentiate 1) normal vs. cancer, 2) low- vs. high-risk cancer, 3) low- vs. intermediate-risk cancer, and 4) intermediate- vs. high-risk cancer voxels. STUDY TYPE Prospective. SUBJECTS Thirty-six patients with biopsy-proven central gland prostate cancer. FIELD STRENGTH/SEQUENCE 3T MRI / 3D 1 H MRSI using GOIA-sLASER. ASSESSMENT Cancer and normal regions of interest (ROIs) were selected by an experienced radiologist and 1 H MRSI voxels were placed within the ROIs to calculate seven metabolite signal ratios. Voxels were split into two subsets, 80% for model training and 20% for testing. STATISTICAL TESTS Four support vector machine (SVM) models were built using the training dataset. The accuracy, sensitivity, and specificity for each model were calculated for the testing dataset. RESULTS High-quality MR spectra were obtained for the whole central gland of the prostate. The normal vs. cancer diagnostic model achieved the highest predictive performance with an accuracy, sensitivity, and specificity of 96.2%, 95.8%, and 93.1%, respectively. The accuracy, sensitivity, and specificity of the low- vs. high-risk cancer and low- vs. intermediate-risk cancer models were 82.5%, 89.2%, 70.2%, and 73.0%, 84.7%, 60.8%, respectively. The intermediate- vs. high-risk cancer model yielded an accuracy, sensitivity, and specificity lower than 55%. DATA CONCLUSION The GOIA-sLASER sequence with an external phased-array coil allows for fast assessment of central gland prostate cancer. The classification offers a promising diagnostic tool for discriminating normal vs. cancer, low- vs. high-risk cancer, and low- vs. intermediate-risk cancer. LEVEL OF EVIDENCE 2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;50:1926-1936.

Purpose: The dynamic keyhole is a new MRimage reconstruction method for thoracic and abdominal MRimaging. To date, this method has not been investigated with cancer patient magnetic resonance imaging (MRI) data. The goal of this study was to assess the dynamic keyhole method for the task of lungtumor localization using cine-MR imagesreconstructed in the presence of respiratory motion. Methods: The dynamic keyhole method utilizes a previously acquired a library of peripheral k-space datasets at similar displacement and phase (where phase is simply used to determine whether the breathing is inhale to exhale or exhale to inhale) respiratory bins in conjunction with central k-space datasets (keyhole) acquired. External respiratory signals drive the process of sorting, matching, and combining the two k-space streams for each respiratory bin, thereby achieving faster image acquisition without substantial motion artifacts. This study was the first that investigates the impact of k-space undersampling on lungtumor motion and area assessment across clinically available techniques (zero-filling and conventional keyhole). In this study, the dynamic keyhole, conventional keyhole and zero-filling methods were compared to full k-space dataset acquisition by quantifying (1) the keyhole size required for central k-space datasets for constant image quality across sixty four cine-MRI datasets from nine lungcancer patients, (2) the intensity difference between the original and reconstructed images in a constant keyhole size, and (3) the accuracy of tumor motion and area directly measured by tumor autocontouring. Results: For constant image quality, the dynamic keyhole method, conventional keyhole, and zero-filling methods required 22%, 34%, and 49% of the keyhole size (P < 0.0001), respectively, compared to the full k-space image acquisition method. Compared to the conventional keyhole and zero-filling reconstructed images with the keyhole size utilized in the dynamic keyhole method, an average intensity difference of the dynamic keyhole reconstructed images (P < 0.0001) was minimal, and resulted in the accuracy of tumor motion within 99.6% (P < 0.0001) and the accuracy of tumor area within 98.0% (P < 0.0001) for lungtumor monitoring applications. Conclusions: This study demonstrates that the dynamic keyhole method is a promising technique for clinical applications such as image-guided radiation therapy requiring the MR monitoring of thoracic tumors. Based on the results from this study, the dynamic keyhole method could increase the imaging frequency by up to a factor of five compared with full k-space methods for real-time lungtumor MRI.

Accurate radiation dose calculation is a major challenge in magnetic resonance imaging (MRI)-only radiation therapy (RT) treatment planning as the required electron density map is not readily provided by this modality. In this work, a number of state-of-the-art synthetic-CT (sCT) generation methods, exhibited promising results in the literature, were evaluated based on common quantitative metrics and patients dataset. This includes four atlas-based approaches, specifically median of atlas images (A-Median) [1], atlas-based voxel-wise weighting (A-VW) [2], bone enhanced atlas-based voxel-wise weighting (A-Bone) [3], iterative atlas-based voxel-wise weighting (A-Iter) [4], and a method based on deep learning convolutional neural network (DL-CNN) [5]. Automatic organ delineation was performed for bladder, rectum and bone. Overall, A-VW, A-Bone, A-Iter and A-VWexhibited comparable performance while DL-CNN showed slightly better segmentation performance resulting in Dice metrics of 0.93, 0.90, and 0.93, respectively. The dosimetric evaluation demonstrated that A-Median, A-VW, A-Bone, A-Iter and DL-CNN resulted in comparable mean dose errors within organs at risk and target volumes showing less than 1% dose difference against the CT-based RT planning. The two-dimensional gamma analysis performed at 1%/1 mm criterion demonstrated comparable pass rates of 94.99±5.15%, 94.59±5.65%, 93.68±5.53% and 93.10±5.99% for A-Bone, DL-CNN, A-Median and A-Iter, respectively. Whereas A-VW and water-only resulted in pass rates of 86.91±13.50% and 80.77±12.10%, respectively. DL-CNN and advanced atlas-based approaches showed promising dosimetric and segmentation accuracy (DL-CNN is slightly better) suggesting that these methods are able to resolve the challenge of synthetic-CT generation from MR images with clinically acceptable errors.

This study is aimed at evaluating the potential role of quantitative magnetic resonance diffusion tensor imaging (DTI) and tractography parameters in the detection and characterization of peripheral zone prostate cancer with a particular attention for fiber tract density.DTI was acquired from eleven high risk, transrectal ultrasound (TRUS)-guided biopsy proven prostate cancers with perineural invasion (histological Gleason score ≥ 7) on a 3 T magnet. Twenty parameters derived from DTI were quantified in cancer and healthy regions of the prostate. In addition, fiber tract density in normal versus cancer tissues was also calculated using DTI tractography. Support vector machine with a radial basis function kernel and area under receiver operator characteristic (ROC) were used to describe and compare the diagnostic performance of combined fractional anisotropy (FA) and mean diffusivity (MD) and other statistically significant DTI parameters. Spearman correlation analysis between DTI parameters and Gleason scores was conducted.Eighteen DTI parameters yielded statistically significant differences between cancer and healthy regions (p-value0.05). The ROC curve of all statistically significant DTI parameters between cancer and healthy regions was higher than the area under ROC curve using FA + MD alone (95% confidence interval = 0.988, range = 0.975-1.00) vs (95% confidence interval = 0.935, range = 0.898-0.999), respectively (p-value0.05). Fiber tract density was also found to be higher in cancer than in healthy tissues (+38.22%, p-value = 0.010) and may be related to the increase in nerve and vascular density reported in prostate cancer. The linear and relative anisotropy were highly correlated with Gleason score (Spearman correlation factor r = 0.655, p-value = 0.001 and r = 0.667, p-value0.001, respectively).DTI has the potential to provide imaging biomarkers in the detection and characterization of prostate cancer. Novel quantitative parameters derived from DTI and DTI tractography, including fiber tract density, support the use of DTI in the assessment of high grade prostate cancer.

PURPOSE: Breathing management can reduce breath-to-breath (intrafraction) and day-by-day (interfraction) variability in breathing motion while utilizing the respiratory motion of internal and external surrogates for respiratory guidance. Audiovisual (AV) biofeedback, an interactive personalized breathing motion management system, has been developed to improve reproducibility of intra- and interfraction breathing motion. However, the assumption of the correlation of respiratory motion between surrogates and tumors is not always verified during medical imaging and radiation treatment. Therefore, the aim of the study was to test the hypothesis that the correlation of respiratory motion between surrogates and tumors is the same under free breathing without guidance (FB) and with AV biofeedback guidance for voluntary motion management. METHODS: For 13 lung cancer patients receiving radiotherapy, 2D coronal and sagittal cine-MR images were acquired across two MRI sessions (pre- and mid-treatment) with two breathing conditions: (a) FB and (b) AV biofeedback, totaling 88 patient measurements. Simultaneously, the external respiratory motion of the abdomen was measured. The internal respiratory motion of the diaphragm and lung tumor was retrospectively measured from 2D coronal and sagittal cine-MR images. The correlation of respiratory motion between surrogates and tumors was calculated using Pearson's correlation coefficient for: (a) abdomen to tumor (abdomen-tumor) and (b) diaphragm to tumor (diaphragm-tumor). The correlations were compared between FB and AV biofeedback using several metrics: abdomen-tumor and diaphragm-tumor correlations with/without ≥5 mm tumor motion range and with/without adjusting for phase shifts between the signals. RESULTS: Compared to FB, AV biofeedback improved abdomen-tumor correlation by 11% (p = 0.12) from 0.53 to 0.59 and diaphragm-tumor correlation by 13% (p = 0.02) from 0.55 to 0.62. Compared to FB, AV biofeedback improved abdomen-tumor correlation by 17% (p = 0.01) and diaphragm-tumor correlation by 15% (p < 0.01) while correcting 0.3 s (p = 0.54) and 0.2 s (p = 0.19) phase shifts, respectively. In addition, AV biofeedback with ≥5 mm tumor motion range, compared to FB improved abdomen-tumor correlation by 14% (p = 0.18) and diaphragm-tumor correlation by 17% (p = 0.01). The highest abdomen-tumor and diaphragm-tumor correlations were found using ≥5 mm tumor motion range and phase shifts, resulting in a 12% improvement in AV biofeedback. CONCLUSIONS: Our results demonstrated that AV biofeedback improves the correlation of respiratory motion between surrogates and the tumor. This suggests a need for AV biofeedback for respiratory guidance utilizing respiratory surrogates during image-guided and MRI-guided radiotherapy in thoracic regions.

Many similarity metrics exist for inter-observer contouring variation studies, however no correlation between metric choice and prostate cancer radiotherapy dosimetry has been explored. These correlations were investigated in this study. Two separate trials were undertaken, the first a thirty-five patient cohort with three observers, the second a five patient dataset with ten observers. Clinical and planning target volumes (CTV and PTV), rectum, and bladder were independently contoured by all observers in each trial. Structures were contoured on T2-weighted MRI and transferred onto CT following rigid registration for treatment planning in the first trial. Structures were contoured directly on CT in the second trial. STAPLE and majority voting volumes were generated as reference gold standard volumes for each structure for the two trials respectively. VMAT treatment plans (78 Gy to PTV) were simulated for observer and gold standard volumes, and dosimetry assessed using multiple radiobiological metrics. Correlations between contouring similarity metrics and dosimetry were calculated using Spearman's rank correlation coefficient. No correlations were observed between contouring similarity metrics and dosimetry for CTV within either trial. Volume similarity correlated most strongly with radiobiological metrics for PTV in both trials, including TCPPoisson (ρ = 0.57, 0.65), TCPLogit (ρ = 0.39, 0.62), and EUD (ρ = 0.43, 0.61) for each respective trial. Rectum and bladder metric correlations displayed no consistency for the two trials. PTV volume similarity was found to significantly correlate with rectum normal tissue complication probability (ρ = 0.33, 0.48). Minimal to no correlations with dosimetry were observed for overlap or boundary contouring metrics. Future inter-observer contouring variation studies for prostate cancer should incorporate volume similarity to provide additional insights into dosimetry during analysis.

Purpose: Radiation treatments are trending toward delivering higher doses per fraction under stereotactic radiosurgery and hypofractionated treatment regimens. There is a need for accurate 3D in vivo patient dose verification using electronic portal imaging device (EPID) measurements. This work presents a model-based technique to compute full three-dimensional patient dose reconstructed from on-treatment EPID portal images (i.e., transmission images). Methods: EPID dose is converted to incident fluence entering the patient using a series of steps which include converting measured EPID dose to fluence at the detector plane and then back-projecting the primary source component of the EPID fluence upstream of the patient. Incident fluence is then recombined with predicted extra-focal fluence and used to calculate 3D patient dose via a collapsed-cone convolution method. This method is implemented in an iterative manner, although in practice it provides accurate results in a single iteration. The robustness of the dose reconstruction technique is demonstrated with several simple slab phantom and nine anthropomorphic phantom cases. Prostate, head and neck, and lung treatments are all included as well as a range of delivery techniques including VMAT and dynamic intensity modulated radiation therapy (IMRT). Results: Results indicate that the patient dose reconstruction algorithm compares well with treatmentmore » planning system computed doses for controlled test situations. For simple phantom and square field tests, agreement was excellent with a 2%/2 mm 3D chi pass rate ≥98.9%. On anthropomorphic phantoms, the 2%/2 mm 3D chi pass rates ranged from 79.9% to 99.9% in the planning target volume (PTV) region and 96.5% to 100% in the low dose region (>20% of prescription, excluding PTV and skin build-up region). Conclusions: An algorithm to reconstruct delivered patient 3D doses from EPID exit dosimetry measurements was presented. The method was applied to phantom and patient data sets, as well as for dynamic IMRT and VMAT delivery techniques. Results indicate that the EPID dose reconstruction algorithm presented in this work is suitable for clinical implementation.« less

Purpose Gantry-mounted megavoltage electronic portal imaging devices (EPIDs) have become ubiquitous on linear accelerators. WatchDog is a novel application of EPIDs, in which the image frames acquired during treatment are used to monitor treatment delivery in real time. We report on the preliminary use of WatchDog in a prospective study of cancer patients undergoing intensity modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT) and identify the challenges of clinical adoption. Methods and Materials At the time of submission, 28 cancer patients (head and neck, pelvis, and prostate) undergoing fractionated external beam radiation therapy (24 IMRT, 4 VMAT) had ≥1 treatment fraction verified in real time (131 fractions or 881 fields). EPID images acquired continuously during treatment were synchronized and compared with model-generated transit EPID images within a frame time (∼0.1 s). A χ comparison was performed to cumulative frames to gauge the overall delivery quality, and the resulting pass rates were reported graphically during treatment delivery. Every frame acquired (500-1500 per fraction) was saved for postprocessing and analysis. Results The system reported the mean ± standard deviation in real time χ 91.1% ± 11.5% (83.6% ± 13.2%) for cumulative frame χ analysis with 4%, 4 mm (3%, 3 mm) criteria, global over the integrated image. Conclusions A real-time EPID-based radiation delivery verification system for IMRT and VMAT has been demonstrated that aims to prevent major mistreatments in radiation therapy.

To clinically implement MRI simulation or MRI-alone treatment planning requires comprehensive end-to-end testing to ensure an accurate process. The purpose of this study was to design and build a geometric phantom simulating a human male pelvis that is suitable for both CT and MRI scanning and use it to test geometric and dosimetric aspects of MRI simulation including treatment planning and digitally reconstructed radiograph (DRR) generation.A liquid filled pelvic shaped phantom with simulated pelvic organs was scanned in a 3T MRI simulator with dedicated radiotherapy couch-top, laser bridge and pelvic coil mounts. A second phantom with the same external shape but with an internal distortion grid was used to quantify the distortion of the MR image. Both phantoms were also CT scanned as the gold-standard for both geometry and dosimetry. Deformable image registration was used to quantify the MR distortion. Dose comparison was made using a seven-field IMRT plan developed on the CT scan with the fluences copied to the MR image and recalculated using bulk electron densities. Without correction the maximum distortion of the MR compared with the CT scan was 7.5 mm across the pelvis, while this was reduced to 2.6 and 1.7 mm by the vendor's 2D and 3D correction algorithms, respectively. Within the locations of the internal organs of interest, the distortion was

The purpose of this study was to investigate performance of the couch and coil mounts designed for MR‐simulation prostate scanning using data from ten volunteers. Volunteers were scanned using the standard MR scanning protocol with the MR coil directly strapped on the external body and the volunteer lying on the original scanner table. They also were scanned using a MR‐simulation table top and pelvic coil mounts. MR images from both setups were compared in terms of body contour variation and image quality effects within particular organs of interest. Six‐field conformal plans were generated on the two images with assigned bulk density for dose calculation. With the MR‐simulation devices, the anterior skin deformation was reduced by up to 1.7 cm. The hard tabletop minimizes the posterior body deformation which can be up to 2.3 cm on the standard table, depending on the weight of volunteer. The image signal‐to‐noise ratio reduced by 14% and 25% on large field of view (FOV) and small FOV images, respectively, after using the coil mount; the prostate volume contoured on two images showed difference of 1.05±0.66 cm3. The external body deformation caused a mean dose reduction of 0.6±0.3 Gy, while the coverage reduced by 22%±13% and 27%±6% in V98 and V100, respectively. A dedicated MR simulation setup for prostate radiotherapy is essential to ensure the agreement between planning anatomy and treatment anatomy. The image signal was reduced after applying the coil mount, but no significant effect was found on prostate contouring. PACS numbers: 87.55.D‐, 87.61.‐c, 87.57.C‐

Purpose: A sensitivity analysis of the effect of variations in electron density data (ED) on dose calculation accuracy for MRI based cervical cancertreatment planning.Methods: Five cervical cancer patients were analysed in this work. Planning CT scans represented gold standard ED data. Standard four field 3DCRT plans (prescription 45Gy) were designed on these CT scans. The CT data was then manipulated to simulate the following methods of assigning ED to MRI; (1) homogenous bulk density corrections, (2) Bulk density correction to bones, (3) rigid image registration of CT to MR, and (4) regression analysis based pseudo CT. Plans were then generated on the manipulated data sets, and compared to the plans generated on the original. Dose was analysed using Chi analysis and equivalent uniform dose (EUD). Data was analysed to quantify (A) the effect on plan design (called optimisation error), and (B) the effect on dose calculation accuracy (systematic error). Results:Analysis of the averaged patient results showed that for 3DCRT, the use of imperfect electron density data had minimal impact on plan design for all tested data sets.Analysis of systematic error showed minimal errors for cases (1), (2) and (3), where average errors of less than 0.3 Gy in EUD were recorded and Chi analysis showed that over 95% of points within the high dose region (D>36Gy) were within 2% or 2mm of the original dose. For case (4), errors greater than .5 Gy in EUD were recorded; these were not considered acceptable errors. Conclusions: Using imperfect electron density data for 3DCRT treatment planning for cervical cancer patients is feasible for appropriately considered choices of electron density assignment. Further analysis is needed to test this result for IMRT, and is ongoing.

Purpose Respiratory variation can increase the variability of tumor position and volume, accounting for larger treatment margins and longer treatment times. Audiovisual biofeedback as a breath-hold technique could be used to improve the reproducibility of lung tumor positions at inhalation and exhalation for the radiation therapy of mobile lung tumors. This study aimed to assess the impact of audiovisual biofeedback breath-hold (AVBH) on interfraction lung tumor position reproducibility and volume consistency for respiratory-gated lung cancer radiation therapy. Methods Lung tumor position and volume were investigated in 9 patients with lung cancer who underwent a breath-hold training session with AVBH before 2 magnetic resonance imaging (MRI) sessions. During the first MRI session (before treatment), inhalation and exhalation breath-hold 3-dimensional MRI scans with conventional breath-hold (CBH) using audio instructions alone and AVBH were acquired. The second MRI session (midtreatment) was repeated within 6 weeks after the first session. Gross tumor volumes (GTVs) were contoured on each dataset. CBH and AVBH were compared in terms of tumor position reproducibility as assessed by GTV centroid position and position range (defined as the distance of GTV centroid position between inhalation and exhalation) and tumor volume consistency as assessed by GTV between inhalation and exhalation. Results Compared with CBH, AVBH improved the reproducibility of interfraction GTV centroid position by 46% ( P = .009) from 8.8 mm to 4.8 mm and GTV position range by 69% ( P = .052) from 7.4 mm to 2.3 mm. Compared with CBH, AVBH also improved the consistency of intrafraction GTVs by 70% ( P = .023) from 7.8 cm 3 to 2.5 cm 3 . Conclusions This study demonstrated that audiovisual biofeedback can be used to improve the reproducibility and consistency of breath-hold lung tumor position and volume, respectively. These results may provide a pathway to achieve more accurate lung cancer radiation treatment in addition to improving various medical imaging and treatments by using breath-hold procedures.

Background and purpose: Recent studies have shown that it is possible to conduct entire radiotherapy treatment planning (RTP) workflow using only MR images. This study aims to develop a generalized intensity-based method to generate synthetic CT (sCT) images from standard T2-weighted (T2(W)) MR images of the pelvis.Materials and methods: This study developed a generalized dual model HU conversion method to convert standard T2(W) MR image intensity values to synthetic HU values, separately inside and outside of atlas-segmented bone volume contour. The method was developed and evaluated with 20 and 35 prostate cancer patients, respectively. MR images with scanning sequences in clinical use were acquired with four different MR scanners of three vendors.Results: For the generated synthetic CT (sCT) images of the 35 prostate patients, the mean (and maximal) HU differences in soft and bony tissue volumes were 16 +/- 6 HUs (34 HUs) and -46 +/- 56 HUs (181 HUs), respectively, against the true CT images. The average of the PTV mean dose difference in sCTs compared to those in true CTs was -0.6 +/- 0.4% (-1.3%).Conclusions: The study provides a generalized method for sCT creation from standard T2(W) images of the pelvis. The method produced clinically acceptable dose calculation results for all the included scanners and MR sequences. (c) 2017 Elsevier B.V. All rights reserved.

Introduction Several methods have been developed to generate pseudo-CT (pCT) from MRI for dose calculation. The objective of this study was to compare an original non-local mean patch-basatlas-based method (ABM) and a bulk density method (BDM) (Fig.1). Methods Thirty-nine patients received VMAT for prostate cancer. T2-weighted MR images were acquired in addition to the planning CT images. pCT were generated from MRI by three methods: an original non-local mean PBM, ABM [1] (leave-one-out strategy) and BDM (water–air-bone density assignment). The PBM was carried out with feature extraction and approximate nearest neighbor search in a training cohort of 25 patients. The methods were compared in a validation cohort of 14 patients. Imaging endpoints were mean absolute error (MAE) and mean error (ME) of Hounsfield unit (HU) from voxel-wise comparisons between pCT and reference CT. Dosimetric endpoints consisted on the absolute mean DVH differences between the reference CT and the pCTs, for the organs at risk and PTV. 3D gamma index analyses (local, 1%/1 mm) were also performed. Results The PBM shows significantly less difference in MAE and ME compared to the ABM and BDM. The differences between pCT from PBM and CT were 0.7% for prostate PTV V 95 % , 0.5% for rectum V 70 Gy and 0.2% for bladder V 50 Gy . Compared to ABM, PBM showed significantly lower DVH difference: for the prostate PTV (from 70 to 78 Gy), the rectum (from 8.5 to 29 Gy, 40 to 48 Gy and 61 to 73 Gy) and the bladder (from 12 to 78 Gy). All the gamma passrate values were significantly higher with the PBM, than those of ABM and BDM. Conclusions In order to generate pCT from MRI for dose calculation, our non-local mean PBM provides lower imaging and dosimetric uncertainties than ABM and BDM. Download : Download high-res image (495KB) Download : Download full-size image

We proposed to perform a basic dosimetry commissioning on a new imager system, the Varian aS1200 electronic portal imaging device (EPID) and TrueBeam 2.0 linear accelerator for flattened (FF) and flattening filter‐free (FFF) beams, then to develop an image‐based quality assurance (QA) model for verification of the system delivery accuracy for intensity‐modulated radiation therapy (IMRT) treatments. For dosimetry testing, linearity of dose response with MU, imager lag, and effectiveness of backscatter shielding were investigated. Then, an image‐based model was developed to convert images to planar dose onto a virtual water phantom. The model parameters were identified using energy fluence of the Acuros treatment planning system (TPS) and, reference dose profiles and output factors measured at depths of 5, 10, 15, and 20 cm in water phantom for square fields. To validate the model, its calculated dose was compared to measured dose from MapCHECK 2 diode arrays for 36 IMRT fields at 10 cm depth delivered with 6X, 6XFFF, 10X, and 10XFFF energies. An in‐house gamma function was used to compare planar doses pixel‐by‐pixel. Finally, the method was applied to the same IMRT fields to verify their pretreatment delivery dose compared with Eclipse TPS dose. For the EPID commissioning, dose linearity was within 0.4% above 5 MU and ∼1% above 2 MU, measured lag was smaller than the previous EPIDs, and profile symmetry was improved. The model was validated with mean gamma pass rates (standard deviation) of 99.0% (0.4%), 99.5% (0.6%), 99.3% (0.4%), and 98.0% (0.8%) at 3%/3 mm for respectively 6X, 6XFFF, 10X, and 10XFFF beams. Using the same comparison criteria, the beam deliveries were verified with mean pass rates of 100% (0.0%), 99.6% (0.3%), 99.9% (0.1%), and 98.7% (1.4%). Improvements were observed in dosimetric response of the aS1200 imager compared to previous EPID models, and the model was successfully developed for the new system and delivery energies of 6 and 10 MV, FF, and FFF modes. PACS number(s): 87.53.Oq, 87.53.Xd

Purpose: In advanced radiotherapy treatments such as intensity modulated radiation therapy(IMRT) and volumetric modulated arc therapy (VMAT), verification of the performance of the multileaf collimator(MLC) is an essential part of the linac QA program. The purpose of this study is to use the existing measurement methods for geometric QA of the MLCs and extend them to more comprehensive evaluation techniques, and to develop dedicated robust algorithms to quantitatively investigate the MLC performance in a fast, accurate, and efficient manner. Methods: The behavior of leaves was investigated in the step-and-shoot mode by the analysis of integrated electronic portal imaging device(EPID)images acquired during picket fence tests at fixed gantry angles and arc delivery. The MLC was also studied in dynamic mode by the analysis of cine EPIDimages of a sliding gap pattern delivered in a variety of conditions including different leaf speeds, deliveries at fixed gantry angles or in arc mode, and changing the direction of leaf motion. The accuracy of the method was tested by detection of the intentionally inserted errors in the delivery patterns. Results: The algorithm developed for the picket fence analysis was able to find each individual leaf position, gap width, and leaf bank skewness in addition to the deviations from expected leaf positions with respect to the beam central axis with sub-pixel accuracy. For the three tested linacs over a period of 5 months, the maximum change in the gap width was 0.5 mm, the maximum deviation from the expected leaf positions was 0.1 mm and the MLC skewness was up to 0.2°. The algorithm developed for the sliding gap analysis could determine the velocity and acceleration/deceleration of each individual leaf as well as the gap width. There was a slight decrease in the accuracy of leaf performance with increasing leaf speeds. The analysis results were presented through several graphs. The accuracy of the method was assessed as 0.01 mm for both the gap size and peak position determination. Conclusions: This study provides fast, easy, and accurate test methods for routine QA of the MLC performance and helps in faster troubleshooting of MLC problems in both IMRT and VMAT treatments.

Dosimetry measurements with Varian amorphous silicon electronic portal imaging devices (a-Si EPIDs) are affected by the backscattered radiation from the EPID support arm. In this study, the nonuniform backscatter from an E-type support arm was reduced by fixing a thick (12.2 × 10.5 × 0.5 cm(3)) piece of lead on top of the arm, and the remaining backscatter was modeled and included in an existing dose prediction algorithm. The applied backscatter kernel was the average of kernels on different regions of the EPID over the arm. The lead-shielded arm reduced the nonuniform backscatter component by about 50% for field sizes ranging from 3 × 3 to 30 × 30 cm(2) and the field symmetry improved for medium to large fields up to 3%. Gamma evaluation of the measured and modeled doses (2%, 2-mm criteria) showed that using the lead-shielded arm in the model increased the number of points with Gamma index1 by 5.7% and decreased the mean Gamma by 0.201. Even using the lead alone (no modeling) could increase the number of points with Gamma index1 by 4.7% and decrease the mean Gamma by 0.153. This is a simple and easy method to decrease the nonuniform arm backscatter and improve the accuracy of dosimetry measurements with the existing EPIDs used for clinical applications.

Purpose : Electronic portal imaging devices (EPIDs) have been studied and used for pretreatment and in-vivo dosimetry applications for many years. The application of EPIDs for dosimetry in arc treatments requires accurate characterization of the mechanical sag of the EPID and gantry during rotation. Several studies have investigated the effects of gravity on the sag of these systems but each have limitations. In this study, an easy experiment setup and accurate algorithm have been introduced to characterize and correct for the effect of EPID and gantry sag during arc delivery. Methods : Three metallic ball bearings were used as markers in the beam: two of them fixed to the gantry head and the third positioned at the isocenter. EPID images were acquired during a 360° gantry rotation in cine imaging mode. The markers were tracked in EPID images and a robust in-house developedMATLAB code was used to analyse the images and find the EPID sag in three directions as well as the EPID + gantry sag by comparison to the reference gantry zero image. The algorithm results were then tested against independent methods. The method was applied to compare the effect in clockwise and counter clockwise gantry rotations and different source-to-detector distances (SDDs). The results were monitored for one linear accelerator over a course of 15 months and six other linear-accelerators from two treatment centers were also investigated using this method. The generalized shift patterns were derived from the data and used in an image registration algorithm to correct for the effect of the mechanical sag in the system. The Gamma evaluation (3%, 3 mm) technique was used to investigate the improvement in alignment of cine EPID images of a fixed field, by comparing both individual images and the sum of images in a series with the reference gantry zero image. Results : The mechanical sag during gantry rotation was dependent on the gantry angle and was larger in the in-plane direction, although the patterns were not identical for various linear-accelerators. The reproducibility of measurements was within 0.2 mm over a period of 15 months. The direction of gantry rotation and SDD did not affect the results by more than 0.3 mm. Results of independent tests agreed with the algorithm within the accuracy of the measurement tools. When comparing summed images, the percentage of points with Gamma index

Purpose: Electronic portal imaging devices (EPIDs) are high resolution systems that produce electronic dose maps with minimal time required for equipment setup, and therefore potentially present a time-saving alternative for intensity modulated radiation therapy (IMRT) pretreatment verification. A modified commercial EPID was investigated operated with an opaque sheet blocking the optical signal produced in the phosphor layer as a precursor to a switched mode dual dosimetry-imaging EPID system. The purpose of this study was to investigate the feasibility of using this system for direct dose to water dosimetry for pretreatment IMRT verification. Methods: A Varian amorphous silicon EPID was modified by placing an opaque sheet between the Gd{sub 2}S{sub 2}O:Tb phosphor layer and the photodiode array to block the optical photons. The EPID was thus converted to a direct-detecting system (dEPID), in which the high energy radiation deposits energy directly in the photodiode array. The copper build-up was replaced with d{sub max} solid water. Sixty-one IMRT beams of varying complexity were delivered to the EPID, to EDR2 dosimetric film and to a 2D ion chamber array (MapCheck). EPID data was compared to film and MapCheck data using gamma analysis with 3%, 3mm pass criteria. Results: The fraction of points thatmore » passed the gamma test was on average 98.1% and 98.6%, for the EPID versus film and EPID versus MapCheck comparisons, respectively. In the case of comparison with film, the majority of observed discrepancies were associated with problems related to film sensitivity or processing. Conclusions: The very close agreement between EPID and both film and MapCheck data demonstrates that the modified EPID is suitable for direct dose to water measurement for pretreatment IMRT verification. These results suggest a reconfigured EPID could be an efficient and accurate dosimeter. Alternatively, optical switching methods could be developed to produce a dual-mode EPID with both dosimetry and imaging capabilities.« less

Purpose:Verification of the mechanical isocenter position is required as part of comprehensive quality assurance programs for stereotactic radiosurgery/radiotherapy (SRS/SRT) treatments. Several techniques have been proposed for this purpose but each of them has certain drawbacks. In this paper, a new efficient and more comprehensive method using cine-EPID images has been introduced for automatic verification of the isocenter with sufficient accuracy for stereotactic applications. Methods: Using a circular collimator fixed to the gantry head to define the field, EPID images of a Winston-Lutz phantom were acquired in cine-imaging mode during 360 deg. gantry rotations. A robust matlab code was developed to analyze the data by finding the center of the field and the center of the ball bearing shadow in each image with sub-pixel accuracy. The distance between these two centers was determined for every image. The method was evaluated by comparison to results of a mechanical pointer and also by detection of a manual shift applied to the phantom position. The repeatability and reproducibility of the method were tested and it was also applied to detect couch and collimator wobble during rotation. Results:The accuracy of the algorithm was 0.03 {+-} 0.02 mm. The repeatability was less than 3 {mu}m andmore » the reproducibility was less than 86 {mu}m. The time elapsed for the analysis of more than 100 cine images of Varian aS1000 and aS500 EPIDs were {approx}65 and 20 s, respectively. Processing of images taken in integrated mode took 0.1 s. The output of the analysis software is printable and shows the isocenter shifts as a function of angle in both in-plane and cross-plane directions. It gives warning messages where the shifts exceed the criteria for SRS/SRT and provides useful data for the necessary adjustments in the system including bearing system and/or room lasers. Conclusions: The comprehensive method introduced in this study uses cine-images, is highly accurate, fast, and independent of the observer. It tests all gantry angles and is suitable for pretreatment QA of the isocenter for stereotactic treatments.« less

Background and purpose Dose planning requires a CT scan which provides the electron density distribution for dose calculation. MR provides superior soft tissue contrast compared to CT and the use of MR-alone for prostate planning would provide further benefits such as lower cost to the patient. This study compares the accuracy of MR-alone based dose calculations with bulk electron density assignment to CT-based dose calculations for prostate radiotherapy. Materials and methods CT and whole pelvis MR images were contoured for 39 prostate patients. Plans with uniform density and plans with bulk density values assigned to bone and tissue were compared to the patient’s gold standard full density CT plan. The optimal bulk density for bone was calculated using effective depth measurements. The plans were evaluated using ICRU point doses, dose volume histograms, and Chi comparisons. Differences in spatial uniformity were investigated for the CT and MR scans. Results The calculated dose for CT bulk bone and tissue density plans was 0.1 ± 0.6% (mean ± 1 SD) higher than the corresponding full density CT plan. MR bulk bone and tissue density plans were 1.3 ± 0.8% lower than the full density CT plan. CT uniform density plans and MR uniform density plans were 1.4 ± 0.9% and 2.6 ± 0.9% lower, respectively. Paired t -tests performed on specific points on the DVH graphs showed that points on DVHs for all bulk electron density plans were equivalent with two exceptions. There was no significant difference between doses calculated on Pinnacle and Eclipse. The dose distributions of six patients produced Chi values outside the acceptable range of values when MR-based plans were compared to the full density plan. Conclusions MR-alone bulk density planning is feasible provided bone is assigned a density, however, manual segmentation of bone on MR images will have to be replaced with automatic methods. The major dose differences for MR bulk density plans are due to differences in patient external contours introduced by the MR couch-top and pelvic coil.

Purpose To assess the accuracy of the initial CT plan dose–volume histograms (DVH's) for prostate, rectum and bladder by comparison to delivered doses determined from cone beam CT (CBCT) scans acquired during image-guided treatment. Materials and methods Twelve prostate patients were treated using daily implanted fiducial guidance and following local protocol for bladder and rectal preparation. CBCT scans were acquired twice weekly and contoured for prostate, rectum and bladder. The planned beams were applied to all CBCT scans to determine the delivered doses. Prostate dose coverage was assessed by the proportion of the CTV fully encompassed by the 95% and 98% isodose lines. Rectal and bladder volumes receiving 40Gy, 60Gy and 70Gy at treatment were compared to the initial plan, with significance determined using the one-sample t -test. Results Four patients showed marginally compromised CTV coverage by the 95% isodose at all CBCT plans. For nine patients the initial plan rectal DVH was significantly outside the range of the treatment DVH's. Conclusions Dose coverage of the prostate was not achieved for all patients. Observed rectal and bladder doses were higher than predicted. The initial treatment plan cannot be assumed to represent accurate normal tissue doses.