Iterative inversion of deformation vector fields with feedback control

Purpose Often, the inverse deformation vector field (DVF) is needed together with the corresponding forward DVF in four‐dimesional (4D) reconstruction and dose calculation, adaptive radiation therapy, and simultaneous deformable registration. This study aims at improving both accuracy and efficiency of iterative algorithms for DVF inversion, and advancing our understanding of divergence and latency conditions. Method We introduce a framework of fixed‐point iteration algorithms with active feedback control for DVF inversion. Based on rigorous convergence analysis, we design control mechanisms for modulating the inverse consistency (IC) residual of the current iterate, to be used as feedback into the next iterate. The control is designed adaptively to the input DVF with the objective to enlarge the convergence area and expedite convergence. Three particular settings of feedback control are introduced: constant value over the domain throughout the iteration; alternating values between iteration steps; and spatially variant values. We also introduce three spectral measures of the displacement Jacobian for characterizing a DVF. These measures reveal the critical role of what we term the nontranslational displacement component (NTDC) of the DVF. We carry out inversion experiments with an analytical DVF pair, and with DVFs associated with thoracic CT images of six patients at end of expiration and end of inspiration. Results The NTDC‐adaptive iterations are shown to attain a larger convergence region at a faster pace compared to previous nonadaptive DVF inversion iteration algorithms. By our numerical experiments, alternating control yields smaller IC residuals and inversion errors than constant control. Spatially variant control renders smaller residuals and errors by at least an order of magnitude, compared to other schemes, in no more than 10 steps. Inversion results also show remarkable quantitative agreement with analysis‐based predictions. Conclusion Our analysis captures properties of DVF data associated with clinical CT images, and provides new understanding of iterative DVF inversion algorithms with a simple residual feedback control. Adaptive control is necessary and highly effective in the presence of nonsmall NTDCs. The adaptive iterations or the spectral measures, or both, may potentially be incorporated into deformable image registration methods.

[1]  Jean-Philippe Thirion,et al.  Image matching as a diffusion process: an analogy with Maxwell's demons , 1998, Medical Image Anal..

[2]  Dietmar Georg,et al.  Adaptive radiation therapy. , 2018, Zeitschrift fur medizinische Physik.

[3]  Michael Brady,et al.  Deformable image registration by combining uncertainty estimates from supervoxel belief propagation , 2016, Medical Image Anal..

[4]  Fréderic Duprez,et al.  Deformation field validation and inversion applied to adaptive radiation therapy , 2013, Physics in medicine and biology.

[5]  W P Segars,et al.  Realistic CT simulation using the 4D XCAT phantom. , 2008, Medical physics.

[6]  Jurgen Fripp,et al.  Robust inverse-consistent affine CT-MR registration in MRI-assisted and MRI-alone prostate radiation therapy , 2015, Medical Image Anal..

[7]  Gary E. Christensen,et al.  Consistent image registration , 2001, IEEE Transactions on Medical Imaging.

[9]  S. Jolly,et al.  Assessment of dosimetric impact of system specific geometric distortion in an MRI only based radiotherapy workflow for prostate , 2017, Physics in medicine and biology.

[10]  Weiguo Lu,et al.  Deformable registration of the planning image (kVCT) and the daily images (MVCT) for adaptive radiation therapy , 2006, Physics in medicine and biology.

[11]  Fang-Fang Yin,et al.  A limited-angle intrafraction verification (LIVE) system for radiation therapy. , 2014, Medical physics.

[12]  Nikos Paragios,et al.  Deformable Medical Image Registration: A Survey , 2013, IEEE Transactions on Medical Imaging.

[13]  Mert R. Sabuncu,et al.  Asymmetric Image-Template Registration , 2009, MICCAI.

[14]  R. Castillo,et al.  Four-dimensional deformable image registration using trajectory modeling , 2010, Physics in medicine and biology.

[15]  M. V. van Herk,et al.  Precise and real-time measurement of 3D tumor motion in lung due to breathing and heartbeat, measured during radiotherapy. , 2002, International journal of radiation oncology, biology, physics.

[16]  Weiguo Lu,et al.  A simple fixed-point approach to invert a deformation field. , 2007, Medical physics.

[17]  Gary E. Christensen,et al.  Consistent Linear-Elastic Transformations for Image Matching , 1999, IPMI.

[18]  E Weiss,et al.  A pseudoinverse deformation vector field generator and its applications. , 2010, Medical physics.

[19]  Thomas Guerrero,et al.  A reference dataset for deformable image registration spatial accuracy evaluation using the COPDgene study archive , 2013, Physics in medicine and biology.

[20]  Paul M. Thompson,et al.  Inverse Consistent Mapping in 3D Deformable Image Registration: Its Construction and Statistical Properties , 2005, IPMI.

[21]  Brian B. Avants,et al.  Symmetric diffeomorphic image registration with cross-correlation: Evaluating automated labeling of elderly and neurodegenerative brain , 2008, Medical Image Anal..