Abstract
Extracting spatial heterogeneities from patient-specific data is challenging. In most cases, it is unfeasible to achieve an arbitrary level of detail and accuracy. This lack of perfect knowledge can be treated as an uncertainty associated with the estimated parameters and thus be modeled as a spatially-correlated random field superimposed to them. In order to quantify the effect of this uncertainty on the simulation outputs, it is necessary to generate several realizations of these random fields. This task is far from trivial, particularly in the case of complex geometries. Here, we present two different approaches to achieve this. In the first method, we use a stochastic partial differential equation, yielding a method which is general and fast, but whose underlying correlation function is not readily available. In the second method, we propose a geodesic-based modification of correlation kernels used in the truncated Karhunen-Loève expansion with pivoted Cholesky factorization, which renders the method efficient even for complex geometries, provided that the correlation length is not too small. Both methods are tested on a few examples and cardiac applications.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Notes
- 1.
Python code is available at the address https://github.com/pezzus/fimh2019.
References
Benito, E.M., et al.: Left atrial fibrosis quantification by late gadolinium-enhanced magnetic resonance: a new method to standardize the thresholds for reproducibility. Europace 19(8), 1272–1279 (2017). https://doi.org/10.1093/europace/euw219
Bolin, D., Kirchner, K., Kovács, M.: Numerical solution of fractional elliptic stochastic PDEs with spatial white noise. IMA J. Numer. Anal. 1–23 (2018). https://doi.org/10.1093/imanum/dry091
Canuto, C., Hussaini, M.Y., Quarteroni, A., Zang, T.A.: Spectral Methods. Springer, Heidelberg (2006). https://doi.org/10.1007/978-3-540-30726-6
Clayton, R.H.: Dispersion of recovery and vulnerability to re-entry in a model of human atrial tissue with simulated diffuse and focal patterns of fibrosis. Front. Physiol. 9(AUG), 1–16 (2018). https://doi.org/10.3389/fphys.2018.01052
Croci, M., Giles, M.B., Rognes, M.E., Farrell, P.E.: Efficient white noise sampling and coupling for multilevel monte carlo with nonnested meshes. SIAM/ASA J. Uncertainty Quant. 6(4), 1630–1655 (2018). https://doi.org/10.1137/18M1175239
Dietrich, C.R., Newsam, G.N.: Fast and exact simulation of stationary gaussian processes through circulant embedding of the covariance matrix. SIAM J. Sci. Comput. 18(4), 1088–1107 (1997). https://doi.org/10.1137/S1064827592240555
Feischl, M., Kuo, F.Y., Sloan, I.H.: Fast random field generation with H-matrices. Numer. Math. 140(3), 639–676 (2018). https://doi.org/10.1007/s00211-018-0974-2
Gneiting, T.: Strictly and non-strictly positive definite functions on spheres. Bernoulli 19(4), 1327–1349 (2013). https://doi.org/10.3150/12-BEJSP06
Haissaguerre, M., et al.: Intermittent drivers anchoring to structural heterogeneities as a major pathophysiological mechanism of human persistent atrial fibrillation. J. Physiol. 594(9), 2387–2398 (2016). https://doi.org/10.1113/JP270617
Harbrecht, H., Peters, M., Schneider, R.: On the low-rank approximation by the pivoted Cholesky decomposition. Appl. Numer. Math. 62(4), 428–440 (2012). https://doi.org/10.1016/j.apnum.2011.10.001
Harbrecht, H., Peters, M., Siebenmorgen, M.: Efficient approximation of random fields for numerical applications. Numer. Linear Algebra Appl. 22(4), 596–617 (2015). https://doi.org/10.1002/nla.1976
Khristenko, U., Scarabosio, L., Swierczynski, P., Ullmann, E., Wohlmuth, B.: Analysis of boundary effects on PDE-based sampling of Whittle-Matérn random fields. arXiv e-prints arXiv:1809.07570 (2018)
Lindgren, F., Rue, H., Lindström, J.: An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach. J. Roy. Stat. Soc. B (Stat. Methodol.) 73(4), 423–498 (2011). https://doi.org/10.1111/j.1467-9868.2011.00777.x
Pathmanathan, P., Gray, R.A.: Validation and trustworthiness of multiscale models of cardiac electrophysiology. Front. Physiol. 9, 1–19 (2018). https://doi.org/10.3389/fphys.2018.00106
Schotten, U., Verheule, S., Kirchhof, P., Goette, A.: Pathophysiological mechanisms of atrial fibrillation: a translational appraisal. Physiol. Rev. 91(1), 265–325 (2011). https://doi.org/10.1152/physrev.00031.2009
Vigmond, E., Pashaei, A., Amraoui, S., Cochet, H., Haïssaguerre, M.: Percolation as a mechanism to explain atrial fractionated electrograms and reentry in a fibrosis model based on imaging data. Heart Rhythm 13, 1536–1543 (2016)
Acknowledgement
The authors acknowledge financial support by the Theo Rossi di Montelera Foundation,the Metis Foundation Sergio Mantegazza, the Fidinam Foundation, and the Horten Foundation to the Center for Computational Medicine in Cardiology. This work was also supported by grants from the Swiss National Supercomputing Centre (CSCS) under project ID s778. Finally, we would like to thank Prof. Michael Multerer for the fruitful discussion.
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2019 Springer Nature Switzerland AG
About this paper
Cite this paper
Pezzuto, S., Quaglino, A., Potse, M. (2019). On Sampling Spatially-Correlated Random Fields for Complex Geometries. In: Coudière, Y., Ozenne, V., Vigmond, E., Zemzemi, N. (eds) Functional Imaging and Modeling of the Heart. FIMH 2019. Lecture Notes in Computer Science(), vol 11504. Springer, Cham. https://doi.org/10.1007/978-3-030-21949-9_12
Download citation
DOI: https://doi.org/10.1007/978-3-030-21949-9_12
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-030-21948-2
Online ISBN: 978-3-030-21949-9
eBook Packages: Computer ScienceComputer Science (R0)