Evaluation of advanced Lukas-Kanade optical flow on thoracic 4D-CT

Extensive use of high frequency imaging in medical applications permit the estimation of velocity fields which corresponds to motion of landmarks in the imaging field. The focus of this work is on the development of a robust local optical flow algorithm for velocity field estimation in medical applications. Local polynomial fits to the medical image intensity-maps are used to generate convolution operators to estimate the spatial gradients. A novel polynomial window function with a compact support is used to differentially weight the optical flow gradient constraints in the region of interest. Tikhonov regularization is exploited to synthesize a well posed optimization problem and to penalize large displacements. The proposed algorithm is tested and validated on benchmark datasets for deformable image registration. The ten datasets include large and small deformations, and illustrate that the proposed algorithm outperforms or is competitive with other algorithms tested on this dataset, when using mean and variance of the displacement error as performance metrics.