A motion model-guided 4D dose reconstruction for pencil beam scanned proton therapy

Abstract

Objective. 4D dose reconstruction in proton therapy with pencil beam scanning (PBS) typically relies on a single pre-treatment 4DCT (p4DCT). However, breathing motion during fractionated treatment can vary considerably in both amplitude and frequency. We present a novel 4D dose reconstruction method combining delivery log-files with patient-specific motion models, to account for the dosimetric effects of intra- and inter-fractional breathing variability.

Approach. Correlation between an external breathing surrogate and anatomical deformations of the p4DCT is established using principal component analysis. Using motion trajectories of a surface marker acquired during dose delivery by an optical tracking system, deformable motion fields are retrospectively reconstructed and used to generate time-resolved synthetic 4DCTs (’5DCTs’) by warping a reference CT. For three abdominal/thoracic patients, treated with respiratory gating and rescanning, example fraction doses were reconstructed using the resulting 5DCTs and delivery log files. The motion model was validated beforehand using leave-one-out cross-validation (LOOCV) with subsequent 4D dose evaluations. Moreover, besides fractional motion,
fractional anatomical changes were incorporated as proof of concept. 

Main results. For motion model validation, the comparison of 4D dose distributions for the original 4DCT and predicted LOOCV resulted in 3%/3 mm gamma pass rates above 96.2%. Prospective gating simulations on the p4DCT can overestimate the target dose coverage V95% by up to 2.1% compared to 4D dose reconstruction based on observed surrogate trajectories. Nevertheless, for the studied clinical cases treated with respiratory-gating and rescanning, an acceptable target coverage was maintained with V95% remaining above 98.8% for all studied fractions. For these gated treatments, larger dosimetric differences occurred due to CT changes than due to breathing variations.

Significance. To gain a better estimate of the delivered dose, a retrospective 4D dose reconstruction workflow based on motion data acquired during PBS proton treatments was implemented and validated, thus considering both intra- and inter-fractional motion and anatomy changes.

A review of the clinical introduction of 4D particle therapy research concepts

“…The 4D dose reconstruction is essential to evaluate the interplay effect and resulting dose deterioration [27] , [28] , [29] , [30] , [31] , [32] , [33] , [98] . CT-based analytical dose calculation algorithms might pose a challenge due to heterogeneities causing Bragg peak degradation [99] .…”

“…In terms of PT with ultra-high FLASH dose rates even more attention should be paid to the timing of motion and beam delivery potentially changing the complexity of intrafractional motion compensation [24] . In cases when motion compensation is not feasible during initial or adaptive treatment planning, the adequate classification of breathing motion and evaluation of the interplay effect becomes crucial in PT [25] , [26] , [27] , [28] , [29] , [30] , [31] , [32] , [33] . While organ movements can be predicted theoretically, the timing of data extraction, processing, synchronisation and correlation modelling needs consideration for a time-efficient treatment planning process.…”

Surveying the clinical practice of treatment adaptation and motion management in particle therapy

“…In addition to reporting on the RRMM practice patterns, Zhang et al included a DELPHI consensus analysis performed by the authors. It named 4D dose calculation including uncertainty evaluation as the most required software feature in the next two years following many recent publications [20] , [21] , [22] , [23] , [24] and recommendations from the TG290 report on respiratory motion management for PT [1] .…”