This is a MATLAB implementation of the Laplace–Dirichlet Rule-Based (LDRB) algorithm for generating ventricular fibers. The original algorithm [1] was adapted to eliminate a discontinuity of fibers in the free walls and to yield a fiber rotation that is directly proportional to the transmural Laplace solution (approximately linear across the wall).
Run dependencies/install.sh to install the dependencies [2,3].
Run examples/example.m to generate fibers for an exemplary biventricular geometry.
All source code is subject to the terms of the GPL-3.0 License.
Copyright 2021 Steffen Schuler, Karlsruhe Institute of Technology.
Steffen Schuler
Institute of Biomedical Engineering
Karlsruhe Institute of Technology
www.ibt.kit.edu
[1] Bayer, J. et al., 2012. A novel rule-based algorithm for assigning myocardial fiber orientation to computational heart models. Ann Biomed Eng.
[2] Jacobson, A., 2018. gptoolbox: Geometry processing toolbox.
[3] Schuler, S., 2020. vtkToolbox: A MEX interface to the VTK library.