Quantifying Tensor Field Similarity with Global Distributions and Optimal Transport

Strain tensor fields quantify tissue deformation and are important for functional analysis of moving organs such as the heart and the tongue. Strain data can be readily obtained using medical imaging. However, quantification of similarity between different data sets is difficult. Strain patterns vary in space and time, and are inherently multidimensional. Also, the same type of mechanical deformation can be applied to different shapes; hence, automatic quantification of similarity should be unaffected by the geometry of the objects being deformed. This work introduces the application of global distributions used to classify shapes and vector fields in the pattern recognition literature, in the context of tensorial strain data. In particular, the distribution of mechanical properties of a field are approximated using a 3D histogram, and the Wasserstein distance from optimal transport theory is used to measure the similarity between histograms. To measure the method's consistency in matching deformations across different objects, the proposed approach was evaluated by sorting strain fields according to their similarity. Performance was compared to sorting via maximum shear distribution (a 1D histogram) and tensor residual magnitude (in perfectly registered objects). The technique was also applied to correlate muscle activation to muscular contraction observed via tagged MRI. The results show that the proposed approach accurately matches deformation regardless of the shape of the object being deformed. Sorting accuracy surpassed 1D shear distribution and was on par with residual magnitude, but without the need for registration between objects.

[1]  Benjamin J. Ellis,et al.  FEBio: finite elements for biomechanics. , 2012, Journal of biomechanical engineering.

[2]  Aart J. Nederveen,et al.  Validation of continuously tagged MRI for the measurement of dynamic 3D skeletal muscle tissue deformation. , 2012, Medical physics.

[3]  Jerry L. Prince,et al.  Construction of An Unbiased Spatio-Temporal Atlas of the Tongue During Speech , 2015, IPMI.

[4]  K. T. Ramesh,et al.  A Three-Dimensional Computational Human Head Model That Captures Live Human Brain Dynamics. , 2017, Journal of neurotrauma.

[5]  Jerry L Prince,et al.  Measuring tongue motion from tagged cine-MRI using harmonic phase (HARP) processing. , 2007, The Journal of the Acoustical Society of America.

[6]  Benjamin Berkels,et al.  Survey of Non-Rigid Registration Tools in Medicine , 2017, Journal of Digital Imaging.

[7]  Yudong Tian,et al.  Performance Metrics, Error Modeling, and Uncertainty Quantification , 2016 .

[8]  Arnold David Gomez,et al.  Motion Estimation with Finite-Element Biomechanical Models and Tracking Constraints from Tagged MRI. , 2017, Computational biomechanics for medicine. From algorithms to models and applications.

[9]  Wei Zeng,et al.  Optimal Mass Transport for Shape Matching and Comparison , 2015, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[10]  Heath B. Henninger,et al.  Validation of computational models in biomechanics , 2010, Proceedings of the Institution of Mechanical Engineers. Part H, Journal of engineering in medicine.

[11]  Bernard Chazelle,et al.  Shape distributions , 2002, TOGS.

[12]  Ryutarou Ohbuchi,et al.  Shape-similarity search of 3D models by using enhanced shape functions , 2003, Proceedings of Theory and Practice of Computer Graphics, 2003..

[13]  David Saloner,et al.  Regional left ventricular myocardial contractility and stress in a finite element model of posterobasal myocardial infarction. , 2011, Journal of biomechanical engineering.

[14]  Marco Cuturi,et al.  Sinkhorn Distances: Lightspeed Computation of Optimal Transportation , 2013, NIPS 2013.

[15]  E. H. Ibrahim,et al.  Myocardial tagging by Cardiovascular Magnetic Resonance: evolution of techniques--pulse sequences, analysis algorithms, and applications , 2011, Journal of cardiovascular magnetic resonance : official journal of the Society for Cardiovascular Magnetic Resonance.

[16]  H. Quynh Dinh,et al.  Measuring the Similarity of Vector Fields Using Global Distributions , 2008, SSPR/SPR.

[17]  Gallagher Pryor,et al.  3D nonrigid registration via optimal mass transport on the GPU , 2009, Medical Image Anal..