Abstract:Purpose: Often, the inverse deformation vector field (DVF) is needed together with the corresponding forward DVF in 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 non-translational 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 6 patients at end of expiration and end of inspiration. Results: NTDC-adaptive iterations are shown to attain a larger convergence region at a faster pace compared to previous non-adaptive 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 non-small NTDCs. The adaptive iterations or the spectral measures, or both, may potentially be incorporated into deformable image registration methods.

Discretize-Then-Optimize Modeling for Dynamic Force Inversion Based on Runge–Kutta Explicit Time Integration

Stable rank-adaptive Dynamically Orthogonal Runge-Kutta schemes

Momentum‐preserving inversion alleviation for elastic material simulation

New explicit time integration algorithm for the rotor dynamics problem

Krylov-based Adaptive-Rank Implicit Time Integrators for Stiff Problems with Application to Nonlinear Fokker-Planck Kinetic Models

A Framework for Data-Driven Computational Dynamics Based on Nonlinear Optimization

A reduced-order model for advection-dominated problems based on Radon-Cumulative-Distribution Transform

A direct linear inversion for discontinuous elastic parameters recovery from internal displacement information only

A TPAA-SBFEM based full-scale forward and inverse interval analysis for dynamic viscoelastic problems

An Inverse Model and Mathematical Solution for Inferring Viscoelastic Properties and Dynamic Deformations of Heterogeneous Structures

The Generalization of Diagonally Implicit Runge–Kutta–Nyström Method with Controllable Numerical Dissipation for Structural Dynamics

Data-driven model reduction via non-intrusive optimization of projection operators and reduced-order dynamics

A nonintrusive distributed reduced‐order modeling framework for nonlinear structural mechanics—Application to elastoviscoplastic computations

Iterative Inversion of Deformation Vector Fields with Feedback Control

A gradient based numerical algorithm to solve inverse dynamic viscoelastic problems of multi-variable identification

Knowledge-dominated and data-driven rigid-flexible coupling dynamics for rotating flexible structure

Friction-adaptive descent: a family of dynamics-based optimization methods

Solving viscoelastic problems in a step-inverse Laplace transform approach supplanted with ARX models: a way to upgrade Finite Element or spectral codes

Accelerating force calculation for dislocation dynamics simulations

Lagrangian interpolation for kinetic dynamic relaxation method with the variable load factor

Time discretization in visco-elastodynamics at large displacements and strains in the Eulerian frame