Physics and Imaging in Radiation Oncology最新文献

The use of and research in automation and artificial intelligence (AI) in radiotherapy is moving with incredible pace. Many innovations do, however, not make it into the clinic. One technical reason for this may be the lack of a platform to deploy such software into clinical practice. We suggest RadDeploy as a framework for integrating containerized software in clinical workflows outside of treatment planning systems. RadDeploy supports multiple DICOM as input for model containers and can run model containers asynchronously across GPUs and computers. This technical note summarizes the inner workings of RadDeploy and demonstrates three use-cases with varying complexity.

Background and purpose

Brachytherapy is treatment of choice for early stage nasal vestibule cancer. Over the years improvements were achieved by means of image guided target definition, interstitial implant techniques and also individual mold techniques. The aim of this study was to improve the technique of the implant so that the need for interstitial catheters can be limited by making use of patient individualized 3D-printed applicators.

Materials and Methods

In 19 patients 3D-printed applicators were used to deliver pulse dose rate (PDR) brachytherapy. All patients underwent computed tomography (CT) and magnetic resonance imaging (MRI). A pre-plan with tumor delineation and manually optimized catheter positions to achieve tumor coverage was made. Based on the pre-plan a 3D-printed applicator was manufactured. Dose was evaluated by several indices: Conformity Index, Healthy Tissues Conformity Index, Dose Homogeneity Index, Dose non-uniformity ratio, Conformal index and high dose (HD) index.

Results

A high target coverage was achieved, with a median V100%_CTV of 99.1 % (range, 81.8–100 %) and median CI of 0.99 (range, 0.82–1.00), as well as a median V0.7Gy_GTV of 100 % (range, 93.0–100 %). The median HD was 0.39 (range, 0.20–0.83). Interstitial catheters were needed in 12 patients. None of the patients developed grade ≥ II toxicity within the median follow up of 18 months.

Conclusions

This study shows that using 3D-printed applicators limits the need for interstitial catheters and also limits the high doses in normal tissue.

Background and purpose

Intensity modulated proton therapy (IMPT) enables generation of conformal dose plans with organ at risk (OAR) sparing potential. However, pelvic IMPT robustness is challenged by inter-fraction motion caused by constant anatomical variations. In this study, the dosimetric impact of inter-fraction motion on target coverage and dose to OAR was quantified in the prospective phase II study ReRad-II on dose-escalated proton reirradiation for locally recurrent rectal cancer (LRRC).

Materials and methods

The inter-fraction motion robustness was assessed for the initial twelve patients enrolled in the ReRad-II study. Patients with resectable LRRC were assessed for neoadjuvant IMPT (55 Gy(RBE)/44Fx) and unresectable recurrences for definitive IMPT (57.5–65 Gy(RBE)/ 46-52Fx). Target coverage and dose to OAR were assessed for robustly optimised three-field IMPT, on 12 plan computerized tomography (CT) scans (pCT) − and 47 repetitive control CT scans (cCTs) during the treatment. The target coverage and doses to OAR were re-calculated on each cCT and the mean dose ratio (pCT/cCT-ratio) and target coverage (V_95%) was evaluated.

Results

The target coverage was robust with a mean dose pCT/cCT-ratio of 1.00 (+/-1%). The V_95% target coverage for every cCT were above the accepted worst-case scenario in the robust evaluation. Considerable variation in bladder-, bowel bag-, and bowel loop volume was observed. The OAR with the largest variation in ratio was the bladder (pCT/cCT-ratio: 1.3 (range: 0.5–4.7).

Conclusions

IMPT for dose-escalated reirradiation of LRRC provided anatomically robust target coverage despite OAR changes. Inter-fraction motion resulted in OAR doses varying within clinically acceptable range.

Background and purpose

Delivery of high precision radiotherapy lymph node boosts requires detailed information on the interfraction positional variation of individual lymph nodes. In this study we characterized interfraction positional shifts of suspected malignant lymph nodes for rectal cancer patients receiving long course radiotherapy. Furthermore, we investigated parameters which could affect the magnitude of the position variation.

Materials and Methods

Fourteen patients from a prospective clinical imaging study with a total of 61 suspected malignant lymph nodes in the mesorectum, presacral, and lateral regions, were included. The primary gross tumor volume (GTV_p) and all suspected malignant lymph nodes were delineated on six magnetic resonance imaging scans per patient. Positional variation was calculated as systematic and random errors, based on shifts of center-of-mass, and estimated relative to either bony structures or the GTV_p using a hierarchical linear mixed model.

Results

Depending on location and direction, systematic and random variations (relative to bony structures) were within 0.6–2.8 mm and 0.6–2.9 mm, respectively. Systematic and random variations increased when evaluating position relative to GTV_p (median increase of 0.6 mm and 0.5 mm, respectively). Correlations with scan time-point and relative bladder volume were found in some directions.

Conclusions

Using linear mixed modeling, we estimated systematic and random positional variation for suspected malignant lymph nodes in rectal cancer patients treated with long course radiotherapy. Statistically significant correlations of the magnitude of the lymph node shifts were found related to scan time-point and relative bladder volume.

Background & purpose

Magnetic resonance imaging (MRI) is increasingly used in treatment preparation of ocular proton therapy, but its spatial accuracy might be limited by geometric distortions due to susceptibility artefacts. A correct geometry of the MR images is paramount since it defines where the dose will be delivered. In this study, we assessed the geometrical accuracy of ocular MRI.

Materials & methods

A dedicated ocular 3 T MRI protocol, with localized shimming and increased gradients, was compared to computed tomography (CT) and X-ray images in a phantom and in 15 uveal melanoma patients. The MRI protocol contained three-dimensional T2-weighted and T1-weighted sequences with an isotropic reconstruction resolution of 0.3–0.4 mm. Tantalum clips were identified by three observers and clip-clip distances were compared between T2-weighted and T1-weighted MRI, CT and X-ray images for the phantom and between MRI and X-ray images for the patients.

Results

Interobserver variability was below 0.35 mm for the phantom and 0.30(T1)/0.61(T2) mm in patients. Mean absolute differences between MRI and reference were below 0.27 ± 0.16 mm and 0.32 ± 0.23 mm for the phantom and in patients, respectively. In patients, clip-clip distances were slightly larger on MRI than on X-ray images (mean difference T1: 0.11 ± 0.38 mm, T2: 0.10 ± 0.44 mm). Differences did not increase at larger distances and did not correlate to interobserver variability.

Conclusions

A dedicated ocular MRI protocol can produce images of the eye with a geometrical accuracy below half the MRI acquisition voxel (<0.4 mm). Therefore, these images can be used for ocular proton therapy planning, both in the current model-based workflow and in proposed three-dimensional MR-based workflows.

Current online adaptive radiotherapy (oART) workflows require dedicated equipment. Our aim was to develop and implement an oART workflow for a C-arm linac which can be performed using standard clinically available tools. A workflow was successfully developed and implemented. Three patients receiving palliative radiotherapy for bladder cancer were treated, with 33 of 35 total fractions being delivered with the cone-beam computed tomography (CBCT)-guided oART workflow. Average oART fraction duration was 24 min from start of CBCT acquisition to end of beam on. This work shows how oART could be performed without dedicated equipment, broadening oART availability for application at existing treatment machines.

Background and Purpose

Standard imaging protocols can guarantee the spatial integrity of magnetic resonance (MR) images utilized in radiotherapy. However, the presence of metallic implants can significantly compromise this integrity. Our proposed method aims at characterizing the geometric distortions induced by both passive and active implants commonly encountered in planning images obtained from a low-field 0.35 T MR-linear accelerator (LINAC).

Materials and Methods

We designed a spatial integrity phantom defining 1276 control points and covering a field of view of 20x20x20 cm³. This phantom was scanned in a water tank with and without different implants used in hip and shoulder arthroplasty procedures as well as with active cardiac stimulators. The images were acquired with the clinical planning sequence (balanced steady-state free-precession, resolution 1.5x1.5x1.5 mm³). Spatial integrity was assessed by the Euclidian distance between the control point detected on the image and their theoretical locations. A first plane free of artefact (FPFA) was defined to evaluate the spatial integrity beyond the larger banding artefact.

Results

In the region extending up to 20 mm from the largest banding artefacts, the tested passive and active implants could cause distortions up to 2 mm and 3 mm, respectively. Beyond this region the spatial integrity was recovered and the image could be considered as unaffected by the implants.

Conclusions

We characterized the impact of common implants on a low field MR-LINAC planning sequence. These measurements could support the creation of extra margin while contouring organs at risk and target volumes in the vicinity of implants.

Background and purpose

MRI-guided online adaptive treatments can account for interfractional variations, however intrafraction motion reduces treatment accuracy. Intrafraction plan adaptation methods, such as the Intrafraction Drift Correction (IDC) or sub-fractionation, are needed. IDC uses real-time automatic monitoring of the tumor position to initiate plan adaptations by repositioning segments. IDC is a fast adaptation method that occurs only when necessary and this method could enable margin reduction. This research provides a treatment planning evaluation and experimental validation of the IDC.

Materials and methods

An in silico treatment planning evaluation was performed for 13 prostate patients mid-treatment without and with intrafraction plan adaptation (IDC and sub-fractionation). The adaptation methods were evaluated using dose volume histogram (DVH) metrics. To experimentally verify IDC a treatment was mimicked whereby a motion phantom containing an EBT3 film moved mid-treatment, followed by repositioning of segments. In addition, the delivered treatment was irradiated on a diode array phantom for plan quality assurance purposes.

Results

The planning study showed benefits for using intrafraction adaptation methods relative to no adaptation, where the IDC and sub-fractionation showed consistently improved target coverage with median target coverages of 100.0%. The experimental results verified the IDC with high minimum gamma passing rates of 99.1% and small mean dose deviations of maximum 0.3%.

Conclusion

The straightforward and fast IDC technique showed DVH metrics consistent with the sub-fractionation method using segment weight re-optimization for prostate patients. The dosimetric and geometric accuracy was shown for a full IDC workflow using film and diode array dosimetry.

Background and purpose

Dynamic trajectory radiotherapy (DTRT) has been shown to improve healthy tissue sparing compared to volumetric arc therapy (VMAT). This study aimed to assess and compare the robustness of DTRT and VMAT treatment-plans for head and neck (H&N) cancer to patient-setup (PS) and machine-positioning uncertainties.

Materials and methods

The robustness of DTRT and VMAT plans previously created for 46 H&N cases, prescribed 50–70 Gy to 95 % of the planning-target-volume, was assessed. For this purpose, dose distributions were recalculated using Monte Carlo, including uncertainties in PS (translation and rotation) and machine-positioning (gantry-, table-, collimator-rotation and multi-leaf collimator (MLC)). Plan robustness was evaluated by the uncertainties’ impact on normal tissue complication probabilities (NTCP) for xerostomia and dysphagia and on dose-volume endpoints. Differences between DTRT and VMAT plan robustness were compared using Wilcoxon matched-pair signed-rank test ($\alpha$ = 5 %).

Results

Average NTCP for moderate-to-severe xerostomia and grade ≥ II dysphagia was lower for DTRT than VMAT in the nominal scenario (0.5 %, p = 0.01; 2.1 %, p < 0.01) and for all investigated uncertainties, except MLC positioning, where the difference was not significant. Average differences compared to the nominal scenario were $\le$ 3.5 Gy for rotational PS ($\le$ 3°) and machine-positioning ($\le$ 2°) uncertainties, <7 Gy for translational PS uncertainties ($\le$ 5 mm) and < 20 Gy for MLC-positioning uncertainties ($\le$ 5 mm).

Conclusions

DTRT and VMAT plan robustness to the investigated uncertainties depended on uncertainty direction and location of the structure-of-interest to the target. NTCP remained on average lower for DTRT than VMAT even when considering uncertainties.

Background and purpose

Despite hardware acceleration, state-of-the-art Monte Carlo (MC) dose engines require considerable computation time to reduce stochastic noise. We developed a deep learning (DL) based dose engine reaching high accuracy at strongly reduced computation times.

Materials and methods

Radiotherapy treatment plans and computed tomography scans were collected for 350 treatments in a variety of tumor sites. Dose distributions were computed using a MC dose engine for $\sim$30,000 separate segments at 6 MV and 10 MV beam energies, both flattened and flattening filter free. For dynamic arcs these explicitly incorporated the leaf, jaw and gantry motions during dose delivery. A neural network was developed, combining two-dimensional convolution and recurrence using 64 hidden channels. Parameters were trained to minimize the mean squared log error loss between the MC computed dose and the model output. Full dose distributions were reconstructed for 100 additional treatment plans. Gamma analyses were performed to assess accuracy.

Results

DL dose evaluation was on average 82 times faster than MC computation at a 1 % accuracy setting. In voxels receiving at least 10 % of the maximum dose the overall global gamma pass rate using a 2 % and 2 mm criterion was 99.6 %, while mean local gamma values were accurate within 2 %. In the high dose region over 50 % of maximum the mean local gamma approached a 1 % accuracy.

Conclusions

A DL based dose engine was implemented, able to accurately reproduce MC computed dynamic arc radiotherapy dose distributions at high speed.