Sensitivity Analysis of Mesh Warping and Subsampling Strategies for Generating Large Scale Electrophysiological Simulation Data

The analysis of large-scale simulation data from virtual populations can be effective to gain computational insight into disease mechanisms and treatment strategies, which can serve for generating hypotheses for and focusing subsequent clinical trials. This can be instrumental in shortening the critical path in medical product development and more cost-effective clinical trials. A previously published pipeline established point correspondence among volumetric meshes to enable meaningful statistics on cardiac electrophysiological simulations on the anatomical distribution of a large-scale virtual population. Thin Plate Splines (TPS), derived from surface deformations, were used to warp a template volumetric mesh, removing the costly operation of repeated volumetric meshing from the pipeline, but potentially at the cost of the volumetric mesh quality. In this work we compare (1) the influence of using TPS versus volumetric meshing of deformed surface meshes, and (2) the influence of surface mesh subsampling prior to the TPS computation. Our results suggest that warping of a template volumetric mesh introduces errors in electrophysiological simulation results of around 4 ms, while having computational times per mesh on the order of seconds, at surface subsampling rates of up to 80%.

[1]  William E. Lorensen,et al.  Decimation of triangle meshes , 1992, SIGGRAPH.

[2]  D. Durrer,et al.  Total Excitation of the Isolated Human Heart , 1970, Circulation.

[3]  N. Simionescu,et al.  The Cardiovascular System , 1983 .

[4]  Dd. Streeter,et al.  Gross morphology and fiber geometry of the heart , 1979 .

[5]  Fred L. Bookstein,et al.  Principal Warps: Thin-Plate Splines and the Decomposition of Deformations , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

[6]  Oscar Camara,et al.  Statistical Atlases and Computational Models of the Heart, First International Workshop, STACOM 2010, and Cardiac Electrophysiological Simulation Challenge, CESC 2010, Held in Conjunction with MICCAI 2010, Beijing, China, September 20, 2010. Proceedings , 2010, STACOM/CESC.

[7]  Alejandro F. Frangi,et al.  Influence of Geometric Variations on LV Activation Times: A Study on an Atlas-Based Virtual Population , 2010, STACOM/CESC.

[8]  Bruce H Smaill,et al.  Three Distinct Directions of Intramural Activation Reveal Nonuniform Side-to-Side Electrical Coupling of Ventricular Myocytes , 2009, Circulation. Arrhythmia and electrophysiology.

[9]  J A Sethian,et al.  A fast marching level set method for monotonically advancing fronts. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[10]  Arthur C. Guyton,et al.  Handbook of Physiology—The Cardiovascular System , 1985 .

[11]  Peter Lindstrom,et al.  Out-of-core simplification of large polygonal models , 2000, SIGGRAPH.

[12]  Cheuk-Man Yu,et al.  Understanding Nonresponders of Cardiac Resynchronization Therapy—Current and Future Perspectives , 2005, Journal of cardiovascular electrophysiology.

[13]  Nassir Navab,et al.  Medical Image Computing and Computer-Assisted Intervention - MICCAI 2010, 13th International Conference, Beijing, China, September 20-24, 2010, Proceedings, Part III , 2010, MICCAI.

[14]  Alejandro F. Frangi,et al.  Efficient computational fluid dynamics mesh generation by image registration , 2007, Medical Image Anal..

[15]  Nick C. Fox,et al.  Phenomenological Model of Diffuse Global and Regional Atrophy Using Finite-Element Methods , 2006, IEEE Transactions on Medical Imaging.

[16]  Yehoshua Y. Zeevi,et al.  The farthest point strategy for progressive image sampling , 1997, IEEE Trans. Image Process..

[17]  Jack Lee,et al.  Personalization of Cubic Hermite Meshes for Efficient Biomechanical Simulations , 2010, MICCAI.