Characterizing a deformable registration algorithm for surface-guided breast radiotherapy

Surface-guided radiation therapy (SGRT) is a nonionizing imaging approach for patient setup guidance, intra-fraction monitoring, and automated breath-hold gating of radiation treatments. SGRT employs the premise that the external patient surface correlates to the internal anatomy, to infer the treatment isocenter position at time of treatment delivery. Deformations and posture variations are known to impact the correlation between external and internal anatomy. However, the degree, magnitude, and algorithm dependence of this impact are not intuitive and currently no methods exist to assess this relationship. The primary aim of this work was to develop a framework to investigate and understand how a commercial optical surface imaging system (C-RAD, Uppsala, Sweden), which uses a nonrigid registration algorithm, handles rotations and surface deformations.

A workflow consisting of a female torso phantom and software-introduced transformations to the corresponding digital reference surface was developed. To benchmark and validate the approach, known rigid translations and rotations were first applied. Relevant breast radiotherapy deformations related to breast size, hunching/arching back, distended/deflated abdomen, and an irregular surface to mimic a cover sheet over the lower part of the torso were investigated. The difference between rigid and deformed surfaces was evaluated as a function of isocenter location.

For all introduced rigid body transformations, C-RAD computed isocenter shifts were determined within 1 mm and 1˚. Additional translational shifts to correct for rotations as a function of isocenter location were determined with the same accuracy. For yaw setup errors, the difference in shift corrections between a plan with an isocenter placed in the center of the breast (BrstIso) and one located 12 cm superiorly (SCFIso) was 2.3 mm/1˚ in lateral direction. Pitch setup errors resulted in a difference of 2.1 mm/1˚ in vertical direction. For some of the deformation scenarios, much larger differences up to 16 mm and 7˚ in the calculated shifts between BrstIso and SCFIso were observed that could lead to large unintended gaps or overlap between adjacent matched fields if uncorrected.

The methodology developed lends itself well for quality assurance (QA) of SGRT systems. The deformable C-RAD algorithm determined accurate shifts for rigid transformations, and this was independent of isocenter location. For surface deformations, the position of the isocenter had considerable impact on the registration result. It is recommended to avoid off-axis isocenters during treatment planning to optimally utilize the capabilities of the deformable image registration algorithm, especially when multiple isocenters are used with fields that share a field edge.