A Simple Penalty that Encourages Local Invertibility and Considers Sliding Effects for Respiratory Motion

Nonrigid image registration is a key tool in medical imaging. Because of high degrees of freedom in nonrigid transforms, there have been many efforts to regularize the deformation based on some reasonable assumptions. Especially, motion invertibility and local tissue rigidity have been investigated as reasonable priors in image registration. There have been several papers on exploiting each constraint separately. These constraints are reasonable in respiratory motion estimation because breathing motion is invertible and there are some rigid structures such as bones. Using both constraints seems very attractive in respiratory motion registration since using invertibility prior alone usually causes bone warping in ribs. Using rigidity prior seems natural and straightforward. However, the “sliding effect” near the interface between rib cage and diaphragm makes problem harder because it is not locally invertible. In this area, invertibility and rigidity priors have opposite forces. Recently, we proposed a simple piecewise quadratic penalty that encourages the local invertibility of motions. In this work we relax this penalty function by using a Geman-type function that allows the deformation to be piecewise smooth instead of globally smooth. This allows the deformation to be discontinuous in the area of the interface between rib cage and diaphragm. With some small sacrifice of regularity, we could achieve more realistic discontinuous motion near diaphragm, better data fitting error as well as less bone warping. We applied this Geman-type function penalty only to the x- and y-direction partial derivatives of the z-direction deformation to address the sliding effect. 192 × 128 × 128 3D CT inhale and exhale images of a real patient were used to show the benefits of this new penalty method.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85922/1/Fessler238.pd

A Simple Regularizer for B-spline Nonrigid Image Registration That Encourages Local Invertibility

Nonrigid image registration is an important task for many medical imaging applications. In particular, for radiation oncology it is desirable to track respiratory motion for thoracic cancer treatment. B-splines are convenient for modeling nonrigid deformations, but ensuring invertibility can be a challenge. This paper describes sufficient conditions for local invertibility of deformations based on B-spline bases. These sufficient conditions can be used with constrained optimization to enforce local invertibility. We also incorporate these conditions into nonrigid image registration methods based on a simple penalty approach that encourages diffeomorphic deformations. Traditional Jacobian penalty methods penalize negative Jacobian determinant values only at grid points. In contrast, our new method enforces a sufficient condition for invertibility directly on the deformation coefficients to encourage invertibility globally over a 3-D continuous domain. The proposed penalty approach requires substantially less compute time than Jacobian penalties per iteration.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85951/1/Fessler21.pd

Regularized Methods for Topology-Preserving Smooth Nonrigid Image Registration Using B-Spline Basis

B-splines are a convenient tool for nonrigid registration, but ensuring invertibility can be challenge. This paper describes a new penalty method that is devised to enforce a sufficient condition for local invertibility and smoothness of nth order B-spline based deformations. Traditional direct Jacobian penalty methods penalize negative Jacobian determinant values only at grid points. In contrast, our new penalty method enforces the sufficient condition for invertibility directly on the B-spline coefficients by using a modified quadratic penalty function so that it enforces invertibility globally over a 3D continuous domain. This approach also saves computation time and memory compared to using Jacobian determinant values. We apply this method to 3D CT images of a thorax at inhale and exhale.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85874/1/Fessler232.pd

Estimation of vector fields in unconstrained and inequality constrained variational problems for segmentation and registration

Vector fields arise in many problems of computer vision, particularly in non-rigid registration. In this paper, we develop coupled partial differential equations (PDEs) to estimate vector fields that define the deformation between objects, and the contour or surface that defines the segmentation of the objects as well.We also explore the utility of inequality constraints applied to variational problems in vision such as estimation of deformation fields in non-rigid registration and tracking. To solve inequality constrained vector field estimation problems, we apply tools from the Kuhn-Tucker theorem in optimization theory. Our technique differs from recently popular joint segmentation and registration algorithms, particularly in its coupled set of PDEs derived from the same set of energy terms for registration and segmentation. We present both the theory and results that demonstrate our approach