 Research
 Open access
 Published:
Mathematical model of the anatomy and fibre orientation field of the left ventricle of the heart
BioMedical Engineering OnLine volume 12, Article number: 54 (2013)
Abstract
Background
One of the main factors affecting propagation of electrical waves and contraction in ventricles of the heart is anisotropy of cardiac tissue. Anisotropy is determined by orientation of myocardial fibres. Determining fibre orientation field and shape of the heart is important for anatomically accurate modelling of electrical and mechanical function of the heart. The aim of this paper is to introduce a theoretical rulebased model for anatomy and fibre orientation of the left ventricle (LV) of the heart and to compare it with experimental data. We suggest explicit analytical formulae that allow us to obtain the left ventricle form and its fibre direction field. The ventricle band concept of cardiac architecture given by TorrentGuasp is chosen as the model postulate.
Methods
In our approach, anisotropy of the heart is derived from some general principles. The LV is considered as a set of identical spiral surfaces, each of which can be produced from the other by rotation around one vertical axis. Each spiral surface is filled with nonintersecting curves which represent myocardial fibres.
For model verification, we use experimental data on fibre orientation in human and canine hearts.
Results
LV shape and anisotropy are represented by explicit analytical expressions in a curvilinear 3D coordinate system. The derived fibre orientation field shows good qualitative agreement with experimental data. The model reveals the most thorough quantitative simulation of fibre angles at the LV middle zone.
Conclusions
Our analysis shows that the band concept can generate realistic anisotropy of the LV. Our model shows good qualitative agreement between the simulated fibre orientation field and the experimental data on LV anisotropy, and the model can be used for various numerical simulations to study the effects of anisotropy on cardiac excitation and mechanical function.
Background
Modern models of complex physiological systems, such as the heart, integrate description from the molecular to the whole organ level and allow researchers to study mechanisms of both mechanical and electrical cardiac activity in normal and pathological hearts.
Over the last several years, a number of models describing electrical and/or mechanical function of the whole heart or its chambers have been proposed [1–12]. The most recent of them consider detailed description of cardiac anatomy and fibre orientation fields as crucial factors for correct representation of the physiological features that are central to heart function.
All approaches to representing cardiac anatomy and anisotropy can be subdivided roughly into two large groups: the individual map approaches, in which fibre orientation is directly measured in the heart using various experimental techniques; and theoretical approaches, in which fibre orientation is generated by algorithms.
In this article, we suggest a theoretical model for anatomy and fibre orientation of the LV. The model is based on the ventricle band concept of cardiac architecture given by TorrentGuasp [13]. In 1972, TorrentGuasp proposed an anatomic concept in which both right and left heart ventricles were considered segments of a single myofibre band twisted and wrapped into a double helical coil [13]. Since that time, this concept has been a subject of intense discussion. Many cardiac anatomists [14, 15] consider the TorrentGuasp hypothesis a gross simplification, and a number of imaging scientists propose a more complex organization of the LV microarchitecture [16]. Another group of researchers has a favorable view on the ventricle band concept [17–19]. For example, an article signed by more than 20 prominent scientists [20] concludes that ‘models such as that of TorrentGuasp et al., which proposes conduction along fibre orientation in a single muscular band and defies conventional concepts of activation, should be investigated’. In spite of that interest, the TorrentGuasp model was never formalized and compared to data on measured cardiac anatomy. Note that in our view, features of the model, such as the possibility of representing realistic fibre orientations by a single warped band, can be proved or disproved only by means of mathematical modelling. Regardless of the outcome, such a formulation will be useful.
In this article, we follow TorrentGuasp’s approach to build cardiac anatomical models of increasing complexity that also use later measurements by Streeter [21]. We show that this description allows one to represent such properties of heart anisotropy as fibre rotation and its dependence on the latitude, spiralling of fibres at the apex and fibre’s maximal angle of torsion about the LV axis. We also perform quantitative comparisons with data from Streeter [21] and Hunter [22], showing good correspondence of the measured fibre orientation fields with that given by our model.
In our model, both the anatomy and fibre orientation field are precisely formulated mathematically. This allows a researcher to apply analytical methods to investigate cardiac electrophysiology and mechanics. In addition, any variations in the shape of the LV and the anisotropy pattern can be reproduced easily by this approach.
The construction of the LV model
The description of the model consists of several steps, starting with simple shapes and approaching the final LV model. Initial steps follow the representation of Pettigrew’s idea (one can see a copy of his figure in Streeter’s work [21] (see Figure three)). We then modify its description and obtain the LV model.
A semicircle with chords on it
Pettigrew began his construction from a semicircle with a set of curves on it [23], as one can see in Streeter’s paper ([21], Figure three a). We describe the figure analytically as follows.
Consider a semicircle with radius K given in the polar coordinate system (P,Φ):
Following Pettigrew [23], let us construct a sequence of horizontal chords (Figure 1):
Here, Φ_{0} and Φ_{1} are polar angles of the right and the left ends of a chord.
A global idea of the model is to wrap this semicircle onto a surface (e.g. conical). The curves will give fibre orientation on that surface, and then rotation of such a surface will give a 3D structure of the heart. The first surface that we will construct is a simple cone. We use a semicircle and a cone because the possibility to wrap a sector of circle to a simple cone is a proven mathematical fact and the LV form closely resembles a cone. We consider the semicircle as the part of the TorrentGuasp “unique myocardial band” which corresponds to the LV. After this wrapping, we will convert the cone to a more complex surface by a nonlinear transformation which will allow us to obtain a more feasible LV model.
Wrapping the semicircle to a cone
In ([21] Figure three (bc)), one can see a wrapping of the semicircle to a surface. We propose the following analytical description of the procedure.
Let us imagine that the semicircle is made of paper. We can wrap this semicircle to a right circular cone (maybe partial or with an overlap) so that the cone vertex corresponds to the semicircle centre (Figure 2). Let us denote an angle along the cone arc as cone twist angle ϕ_{max}. We get a cone that becomes closed, if ϕ_{max} = 2 π. We are going to consider only the case where ϕ_{max} > π. The cone is specified in a cylindrical coordinate system (ρ, ϕ, z) as follows.
First note that curves described by Φ = const do not bend during wrapping. Therefore, these curves Φ = const are the generatrices of the cone. The parametric equations of the cone are
An explicit equation for the object is
Let us note that the semicircumference, limiting the given semicircle, transforms to the cone planar arc, and the centre of the circle becomes the cone’s apex.
In order to model not only muscular layers, but also myofibres, one has to look for the position of the transformed chords after the wrapping.
Let us find the chord Y = const images on the conical surface, by substituting (4) into (5), (6), (7):
The results can be seen in Figure 3. It is not difficult to see that rotation of such a simple conical surface around the vertical axis does not give a good representation of the heart anatomy, as it will produce only a conical surface, i.e. a body of zero thickness. To improve that, we will generalize the procedure by introducing dependency of the generatrices on the rotation angle.
Construction of spiral surfaces
In Figure 2, the cone’s generatrix was z = k ρ, where k was a constant (see (8)), so that it did not depend on the angle ϕ. Let us consider the more general situation when z = _{ ϕ }(ρ). Such a generatrix will generate a spiral surface, which will finally give us a proper representation of the heart’s geometry.
Let us consider a few examples.
We can first assume that generatrix changes its slope as shown in Figure 4. It can be formally represented as
where H,R,r are positive constants, and
with 0 ≤ γ ≤ 1. The domain of _{ ϕ } is taken
such that the codomain is
The conical shape produced by such generatrices is shown in Figure 4, on the right. We will call it a pseudoconical surface.
If we rotate a pseudoconical surface around the vertical axis, we get a conical body that has some resemblance to the LV (it has an LV cavity), but its thickness at the apex will be zero. To improve, we modify the generatrices as follows.
Let us move the end A (see Figure 4) down by a value h γ, i.e. proportional to the angle ϕ = γ ϕ_{max} (h is a positive constant), as shown in Figure 5:
As a result, the thickness of the LV at the apex will become h > 0, which improves our representation. However, a real LV surface is not conical, so we have to use curves as generatrices instead of straight line segments.
Let us connect the same points A and B, as in the previous example, but by an arc (Figure 6):
where the function _{ ϕ } represents a curved generatrix.
The following properties are imposed on the function _{ ϕ }(ρ)

1.
D(_{ ϕ }(ρ)) = E(_{ ϕ }(ρ)) = [0,1]

2.
_{ ϕ }(ρ) is continuous and differentiable on [0,1]

3.
_{ ϕ }(ρ) increases monotonically on [0,1].
We refer to the functions with these properties as generating functions (GF). In the first example {\mathcal{F}}_{\varphi}\left(t\right)=t,\phantom{\rule{1em}{0ex}}{\mathcal{Z}}_{\varphi}\left(\rho \right)={\mathcal{F}}_{\varphi}\left(\frac{\rho}{r+\gamma (Rr)}\right)\xb7H; in the second example _{ ϕ }(t) = t, and is represented through in the same way as in the third example (see (15)). Now let us choose a proper GF to construct a more realistic LV.
A model of the LV
To represent the shape of the epicardial LV surface, Streeter ([21], pp. 91 – 92) used the following functions:
If ϵ = 0, the curve AB is a line segment, and if ϵ = 1, it is a quarter of an ellipse. For intermediate ϵ values, we get intermediate curves. Let us use the following equation of the endocardium in analogy with Streeter’s description:
where “latitude” ψ takes values 0 ° ≤ ψ ≤ 90° (Figure 7). Parameters determining the shape are: an outer radius R_{ b } near the equator; a thickness L near the equator; a height Z_{ b }; a thickness h at the apex.
The form of any intermediate layer between epi and endocardium is described by the equations (see also Additional file 1):
Elimination of ψ delivers an explicit equation for the surface:
where _{ϵ sp}: [0,1] → [0,1] is the GF of the epicardium, the endocardium and every intermediate layer of the LV wall; ρ_{epi}/R_{ b } ∈ [0,1], ρ_{endo}/(R_{ b }  L) ∈ [0,1].
As a result, we get the following definition of the ϵ  spiral surface (ESS, see Additional files 2 and 3):
The spiral surface’s border with an angle ϕ = 0 is located at the endocardial side, and the border with an angle ϕ = ϕ_{max} lies at the epicardial side. The LV model is made as a body of revolution of an ESS. See Figure 8 for an example of a single and multiple nested ESS. We form the LV by using shifted layers (or the rotated spiral surfaces, which is the same) because a thick muscular layer cannot be wrapped to a body of revolution with overlapping, but without any shift of its sheets.
Our description reproduces the form of epicardium from Streeter’s work (see Figure forty in [21]), but for endocardium, Streeter uses more complex equations with an additional parameter that he called the ‘angle of taper’. However, we found that even without this parameter our function reasonably reproduces the form of the endocardium, and, thus, we decided to use equation (23) without further modifications.
To represent fibre orientation and to compare it with anatomical data, Streeter [21] used the following angles: true fibre angle, α; the helix angle, α_{1}; and the longitudinal angle, α_{2} (see Figure 9 for their definitions). We follow the same approach in this work.
Together with (22), (23) forms the basis for the LV model. Finally, we need to set proper fibre angles at the epicardial and endocardial surfaces.
In the model given by (22) and (23), angle α depends on point position in the LV wall and changes from 90° on the endocardium to approximately 0° in the middle of the wall and then to 90° on the epicardium. In real hearts, the rotation of fibre is less, and its values at the endocardium and epicardium are about 60° and 70° ([21], Figure thirtythree). To account for that, we use only part of the interval 0≤γ_{0} ≤ γ ≤ γ_{1} ≤ 1. For example, if γ_{0} = 0.1 ≤ γ ≤ 0.75 = γ_{1}, then the extreme angle values in the equatorial area are 55° at the endocardium and 75° at the epicardium (see Figure 10 and Results for more details).
As a change in the range of γ changes the anatomy of the heart, we need to rescale it to the normal heart using:
After choosing of γ_{0} and γ_{1}, we completely specify parameters of our model and can use it for generation of anisotropy.
Comparison of the theoretical model with experimental data
In this section, we indicate the parameter values used in our study and compare theoretically obtained results with experimental data from [21] and [22].
Verification of the model: comparison with Streeter’s data
We used the following parameter values reported in ([21], Table two): external radius of LV at the equator {R}_{b}^{e}=33 mm, thickness of LV wall at the equator L^{e} = 10 mm, height of underequatorial LV part {Z}_{b}^{e}=60 mm, thickness of LV wall at the apex h^{e} = 7 mm, ϵ^{e} = 0.9; and we set angle of spiral surface torsion {\varphi}_{\text{max}}^{e}=3\pi according to ([21], Figure three c).
See Figure 11 and 3 additional movie files for the spiral surfaces we made using these parameter values.
We used the parameters of the subepicardial and subendocardial boundaries γ_{0} = 0.1, and γ_{1} = 0.75, which gives α = 55° at the endocardium and θ = 75° at the epicardium that fits with ([21], Figure fortyfour) taking into account the experimental measurement error and data variability.
We compared three angular characteristics of our myofibres field with experimental data from [21]. The comparison was made in two areas of the LV: the equatorial (58 mm ≤z ≤ 60 mm) and the bottom (18 mm ≤ z ≤ 21 mm) parts of the LV. In both regions, angles were compared along a line orthogonal to the epicardium, which is common in anatomical studies.
The results of comparison are shown in Figures 10, 12 and 13.
LV top zone
Note that in the LV top zone, the true fibre angle α varies nonmonotonically, as in Streeter’s data [21]. In particular, α on the endocardium was close to 55°, on the epicardium it was close to 75° and decreased to 5° approximately at the middle of the LV wall, between 35% and 40% of the wall depth (Figure 10, on the left). We see some differences in the slope of curves in our model and in measured data; however, a qualitative correspondence is observed.
Figure 12, on the left, shows similar results for the helix angle. We see that going from endocardium to epicardium, the helix angle α_{1} decreased monotonically from +55° to 75° and was equal to 0° approximately at the middle of the LV wall, between 35% and 40% of the wall depth. We also see a good agreement of our model with the experimental results from [21].
As we see in Figure 13, on the left, the longitudinal angle, α_{2}, going from endocardium to the middle of the LV wall, decreased from 0° to 90°. At the middle of the LV wall, between 35% and 40% of the wall depth, the experimental and modelled angle abruptly changed from 90 to +90°, which is an artefact. We found the jump at 43%, while it was 38% in Streeter’s data. At the exterior half of the LV wall, the α_{2} angle decreased at a decelerating rate to 0°. Here too, we observe good qualitative correspondence of our model with experimental data, although some quantitative differences in the slopes of the dependencies are present.
LV middle zone
The true fibre angle α also reaches a minimum in the midwall region, both in our model and in Streeter’s data. On the endocardium, it was close to 35°, on the epicardium, it was close to 50° and decreased to 5° approximately at the middle of the LV wall, at 0.4 of the wall depth (Figure 10, on the right).
Figure 12, on the right, displays similar results for the helix angle. We see that going from endocardium to epicardium, the helix angle, α_{1}, decreased monotonically from +25° to 50° and was equal to 0° approximately at the middle of the LV wall, at 0.4 of the wall depth. We also see a good agreement of our model with the experimental results from [21].
In Figure 13, on the right, the longitudinal angle, α_{2}, decreased from 5° to 90° between endocardium and midwall. At the middle of the LV wall, between 0.35 and 0.4 of the wall depth, the angle again abruptly changed to +90°. On the exterior half of the LV wall, the α_{2} angle decreased at a decelerating rate to 0°. A good qualitative correspondence with anatomical data is found in the outer twothirds of the LV wall.
Our model also successfully reproduces the distribution of fibre directions in the thick, radially placed layer of the LV wall, as shown in Figure 14. Streeter studied the pattern and named it ‘Japanese fan’ ([21], Figure fortytwo c).
Overall we can claim that our model adequately reproduces the direction of myocardial fibres in the human LV.
Comparison with Auckland canine dataset
We also compared our model to the data used in [22]. The comparison was conducted in the following way. We aligned our LV model with the Auckland model along the vertical axis and found that reasonable fit of our anatomy to the Auckland heart model occurs for: {R}_{b}^{e}=45 mm, {Z}_{b}^{e}=80 mm, h^{e} = 12 mm, L^{e} = 15 mm, ϵ^{e} = 0.85, {\varphi}_{\text{max}}^{e}=3\pi, γ_{0} = 0.15, γ_{1} = 0.9. Subsequently we constructed five meridional halfplanes so that they divided the LV free wall region and the dihedral angle to four equal dihedral angles. We examined only points lying near the five halfplanes. We used the same comparison procedure that we described in Section “Verification of the model: comparison with Streeter’s data”. The fibre orientation was compared at three different “latitudes” (close to apex, in the middle and close to base) along lines orthogonal to the heart surface. We computed the x values (they show the position of a point in the ventricular wall) and the two angles, α and α_{1}, and plotted their values from the experimental dataset and from the theoretical model. The results obtained are shown in Figure 15.
We see that the plots of the α angle have characteristic Vshaped forms, both in the model and in four of five meridional sections at all three “latitudes”. The model angle is within the limits of the section angles in almost all the positions, except the external onethird of the wall on the left graph. The plots of the α_{1} angle show the same pattern: they descend from some positive values to some negative ones, speed of the decrease is higher in the endo and epicardial area and slightly less in the midmyocardium.
Note that we show only one line for our model. This is because our model is rotationally symmetric, and all five lines for different rotation angles are the same. The model is an idealization and it represents only some averaged characteristics of LV and neglects possible individual aspects of the LV form.
Our simulation programme is available on request.
Discussion
In this section, we discuss limitations of our model, experimental methods and models that can be used to verify it, its possible applications and ways of improving it.
Limitations
Our model adequately reproduces fibre angles at the LV middle zone, but agreement for the LV top and apex zones is mainly qualitative.
Inaccurate reproduction of fibre direction in the basal and apical zones can be caused by different reasons, one of which can be the idealized axisymmetrical LV form in the model. Real ventricles have essential deflections from the axis symmetry; they are somewhat individual, that is specific for every LV. In this work, we tried to construct an axisymmetrical model that could maximally suit the middle (by height) LV area. The apex and base of a real LV could considerably move away from the axis. We are currently working on extensions of the model which will adapt the idealized model to individual peculiarities of real ventricles. We hope to achieve more accurate imitation of the transmural run of the fibres at the apical and basal zones.
Experimental methods that can be used to verify our model
Currently, there are several experimental methods that can be used to measure fibre orientation in the heart. One of them is the diffusion tensor imaging (DTI) technique. In this approach, a researcher finds the diffusion matrix of water molecules in the heart. The main directions of diffusion are determined by the structure of the tissue [24–27], and by calculating the eigenvectors of the matrix corresponding to the largest eigenvalue, the fibre direction can be found. Diffusion tensor magnetic resonance imaging (DTMRI) measurements can be done with spatial resolutions up to 200 μ m [28]. Another advanced technique is the micro computerized tomography (MicroCT) imaging. MicroCT measurements can be done with spatial resolutions up to 36–70 μ m [29], and both these methods produce highquality data that can be used in computer models.
Direct measurement of anisotropy also can be conducted via tedious histological studies of fibre direction in 3D [21, 22]. In this method, the researcher makes a series of parallel sections of the heart. In each section, angles of fibre slope are measured, which gives a full picture of the fibre directions in the heart. Recently, Smaill et al. developed a combined highresolution serial imaging microscopy technique [30]. In this method, after heart fixation, they perform crosssection and make successive, highresolution images of the heart. Then, using computer processing, the data are collected to form an overall 3D dataset.
In [31], the myofibres in the foetal human heart are investigated using quantitative polarized light microscopy. The hearts are embedded in a transparent resin, polymerized and then sectioned. Afterward, the elevation and azimuth angles are measured by means of polarized light (see [31] for details).
Comparison with other models
Experimental data on the fibre orientation obtained as described above may be in different ways used for construction of anatomical computational models:

either as a discrete dataset in finite element models [32–34];

or for the verification of rulebased models, i.e. the models formed on the base of some constitutive rule [35–38].
One of the most recent rulebased methods is a Laplace–Dirichlet algorithm [37], which takes a noisy DTIderived fibre orientation field as input data and yields, firstly, the transmural and apicobasal directions for the entire myocardium and, secondly, a smooth and continuous fibre orientation field. Another approach was used by Peskin, who derived fibre orientation field from the principles of mechanical equilibrium [39]. One more anatomic model based on a principle of the mechanical activity of the heart was a model by Chadwick, who considered a cylindrical LV and specified the helix angle linearly depending on point position in the LV wall [40]. Beyar et al. shaped LV into a spheroid and also used linear dependence of the helix angle on the distance between any point in the ventricle wall and the endocardium [2]. An interesting example of the theoretical approach was developed by Arts et al. in 1992 [41]. They constructed a model of an ellipsoidal LV, complicated the law of helix angle change to a piecewise quadratic one and quantified orientation of the muscle fibres via the helix fibre angle distribution, which was found upon application of the mechanical adaptative principle suggested by Arts et al. in 1982 [42].
In this work we present our rulebased model focused on the LV morphology including simulation of the ventricle shape and fibre orientation in its wall. The developed formalism is substantially associated with both ventricle band concept of cardiac architecture given by TorrentGuasp [13] and anatomic observations presented by Streeter in his classical work [21]. In our approach to the modelling of the LV architecture, anisotropy of the heart was derived from some general principles. In our model, the LV is considered a set of identical spiral surfaces combined with each other by rotation about the vertical axis. Every spiral surface is defined analytically and represents a mapping of a half disc. The first step of the transformation is the mapping of the semicircle to a conical surface. In the second step, the conical surface is transformed to a curved spiral surface, representing the quasielliptical shape of the LV boundary surface. Finally, every spiral surface is filled with myocardial fibres, represented by the transformed images of the chords that were parallel to the diameter in the initial semicircle (see Figure 1).
Our model is not the only wrappingbased myocardium model. Sinha et al. proposed in [43] a model of one myocardial layer which had a rectangular form and was wrapped around a (truncated) cone. They used this very simple model to study termination of reentrant waves rotating around obstacles in cases of isotropy and anisotropy but without any linkage to the real fibre pattern in the heart ventricles.
We used experimental data from the above cited work by Streeter [21], as well as from other more recent works [22] for the model verification. In particular, the model proved to reproduce adequately both the looping arrangement of the muscle fibres and the specific 3D pattern of the relative positions of the fibres in the transmural direction through the ventricle wall.
These accurate reproductions allow us to consider the model a touchstone in validating the ventricle band concept of cardiac architecture originated by TorrentGuasp, because the model, based on this concept, yields an adequate fibre field as a consequence of the postulates.
It seems reasonable to compare our model with another rulebased model that assigns fibre orientation locally, particularly with the very interesting and promising model by Bayer et al., mentioned earlier in this section [37]. For this comparison, only the reproduced fibre orientation in various parts of the LV of the two models can be used. The Bayer model is based on DTMRI data of anisotropy in two ventricles of a canine heart. The average angle divergence between the model and the DTMRI fibre directions is 23°; that is there is not a complete quantitative matching of the real experimental and reproduced data, but there is reasonable concordance. Specifically, Bayer’s model quantitatively reproduces fibre anisotropy in the basal and apical LV zones better than our model (see [37], Figure three). In the section describing the limitations of our model, we point to this quantitative inaccuracy in our model and propose some ways to eliminate it. At the same time, our model better reproduces the experimental data in the middle LV zone; namely, we obtain the specific slike plot of the angle α_{1} in the transmural direction (see Figure 15, the bottom right panel). In Bayer’s model, this dependency is linear by definition. Moreover, if we follow the cited paper by Bayer et al. ([37], Eq. (1), (2)), all plots for angle α_{1} presented in Figure 12 reveal independence of the angle from both the latitude and the longitude of the intramural position within the wall. Bayer and coauthors suggest that their model can be easily improved to take the nonlinearity of angle α_{1} into account. But, it also is necessary to make the anisotropy latitudinally and longitudinally independent, and it is not easy to do so. Our model reproduces such a dependence (see Figures 10, 12, 13, 14 and 15), which is proper for real hearts, and does it quite fairly for the middle zone of the LV.
One more simplification of Bayer’s model concerns transmural rotation of the fibres’ directions, named ‘Japanese fan’ by Streeter ([21], Figure fortytwo c). In that model, the rotation is defined in one plane only, that is, around only one axis, settled transmurally. This plane lies tangentially to the surface determined by the condition d = const, where d is a term specified in the cited paper by Bayer et al. [37] and presents there the occurrence depth of particular locus within the wall; for example, d = 0 on the endocardium and d = 1 on the epicardium. Moreover, if we assess results obtained in Bayer’s model by means of the angle α_{3} defined by Streeter [21] and determine transmural direction of the fibre orientation, it will prove to be constantly 0, which is a substantial simplification. This feature does not allow mapping the 3D pattern of the relative positions of the fibres in the transmural direction through the ventricular wall.
In contrast, in our model the 3D pattern is taken into consideration (see Figures 10, 12, 13, 14 and 15, and note that in the middle LV zone these angles are reproduced quite well).
Thus, we believe that both models have their own virtues and their own limitations, and further development of the models would be useful to overcome the limitations.
Development and uses of our model
We suggest that analytical representation of the geometry presented here can be used for development of new numerical methods to study electrical and mechanical activity of the heart. As our model provides analytical mapping of a rectangle in (γ,ψ,ϕ) space into the curvilinear heart shape, one can formulate a rectangular numerical scheme in (γ,ψ,ϕ) space (where representation of boundary conditions is simplest) and account for anisotropy by explicit analytical formulae. The model can also be used to generate various anisotropic properties of the heart and modulations of the LV shape (via model parameter variations) and to study their effects on electrical and mechanical heart functions.
Of course, our model is an idealization; it represents some averaged characteristics of the LV and neglects possible individual aspects of the LV form. Model adaptation to individual characteristics is the subject of a particular line of research. We are developing methods to further modifying the LV model to customize it to individual hearts.
Conclusions
We have constructed one of the simplest analytical descriptions of cardiac anatomy based on the TorrentGuasp’s ventricular band concept. The model can be used for verification of the band concept as well as for various numerical simulations to study the effects of anisotropy on cardiac excitation and mechanical function. Our model shows good qualitative agreement between the simulated fibre orientation field and the experimental data on LV anisotropy.
Appendix
Mapping a point on the semicircle to the spiral surface
In this section, we give all needed formulae to map a point on the semicircle to the cone and then to the spiral surface. The following input parameters are used

external radius of LV at the equator, R_{ b };

thickness of LV wall at the equator, L;

height of LV, Z_{ b };

thickness of LV wall at the apex, h;

conicityellipticity parameter, ϵ;

angle of spiral surface torsion, ϕ_{max}.
Let us consider a point (P,Φ) on a semicircle with a radius K. In order to get outer radius R_{ b } of the LV on the equator, we need to use
The image of the point on the ESS has the following coordinates (see also (23)):
The mapping (x,y,z) → (γ,ψ,ϕ) is given by
The fibre direction \overrightarrow{v}=(\mathit{\text{vx}},\mathit{\text{vy}},\mathit{\text{vz}}) at a point (γ,ψ,ϕ):
Abbreviations
 LV:

Left ventricle
 DT(I):

Diffusion tensor (imaging)
 MRI:

Magnetic resonance imaging
 CT:

Computer tomography
 GF:

Generating function.
References
Trayanova N, Constantino J, Gurev V: Electromechanical models of the ventricles. Am J Physiol Heart Circ Physiol 2011, 301: H279H286. 10.1152/ajpheart.00324.2011
Beyar R, Sideman S: A computer study of the left ventricular performance based on fiber structure, sarcomere dynamics, and transmural electrical propagation velocity. Circ Res 1984, 55: 358–375. 10.1161/01.RES.55.3.358
Greenstein J, Winslow R: Integrative systems models of cardiac excitationcontraction coupling. Circ Res 2011, 108: 70–84. 10.1161/CIRCRESAHA.110.223578
Hunter P, McCulloch A, Keurs HT: Modelling the mechanical properties of cardiac muscle. Prog Biophys Mol Biol 1998, 69: 289–331. 10.1016/S00796107(98)000133
Tusscher KT, Noble D, Noble P, Panfilov A: A model for human ventricular tissue. Am J Physiol Heart Circ Physiol 2004, 286: H1573H1589.
Tusscher KT, Panfilov A: Alternans and spiral breakup in a human ventricular tissue model. Am J Physiol Heart Circ Physiol 2006, 291: H1088H1100. 10.1152/ajpheart.00109.2006
Grandi E, Pasqualini F, Bers D: A novel computational model of the human ventricular action potential and Ca transient. J Mol Cell Cardiol 2010, 48: 112–121. [Special Issue: Ion Channels] 10.1016/j.yjmcc.2009.09.019
O’Hara T, Virag L, Varro A, Rudy Y: Simulation of the undiseased human cardiac ventricular action potential: model formulation and experimental validation. PLoS Comput Biol 2011, 7(5):e1002061. 10.1371/journal.pcbi.1002061
Sulman T, Katsnelson L, Solovyova O, Markhasin V: Mathematical modeling of mechanically modulated rhythm disturbances in homogeneous and heterogeneous myocardium with attenuated activity of Na+  K+ pump. Bull Math Biol 2008, 70(3):910–949. 10.1007/s115380079285y
Niederer S, Hunter P, Smith N: A quantitative analysis of cardiac myocyte relaxation: a simulation study. Biophys J 2006, 90(5):1697–1722. 10.1529/biophysj.105.069534
Keurs HT, Shinozaki T, Zhang Y, Zhang M, Wakayama Y, Sugai Y, Kagaya Y, Miura M, Boyden P, Stuyvers B, Landesberg A: Sarcomere mechanics in uniform and nonuniform cardiac muscle: a link between pump function and arrhythmias. Prog Biophys Mol Biol 2008, 97(2–3):312–331. 10.1016/j.pbiomolbio.2008.02.013
Rice J, de Tombe P: Approaches to modeling crossbridges and calciumdependent activation in cardiac muscle. Prog Biophys Mol Biol 2004, 85(2–3):179–195.
TorrentGuasp F: The Cardiac Muscle. Madrid: Fundacion Juan March; 1973.
Anderson R, SanchezQuintana D, Redmann K, Lunkenheimer P: How are the myocytes aggregated so as to make up the ventricular mass? Semin Thorac Cardiovasc Surg Pediatr Card Surg Annu 2007, 10: 76–86. 10.1053/j.pcsu.2007.01.016
Anderson RH, Smerup M, SanchezQuintana D, Loukas M, Lunkenheimer PP: The threedimensional arrangement of the myocytes in the ventricular walls. Clin Anatomy 2009, 22: 64–76. 10.1002/ca.20645
LeGrice I, Smaill B, Chai L, Edgar S, Gavin J, Hunter P: Laminar structure of the heart: ventricular myocyte arrangement and connective tissue architecture in the dog. Am J Physiol 1995, 269(2):H571H582.
Kocica MJ, Corno AF, Lackovic V, Kanjuh VI: The helical ventricular myocardial band of TorrentGuasp. Semin Thorac Cardiovasc Surg: Pediatr Card Surg Annu 2007, 10: 52–60. 10.1053/j.pcsu.2007.01.006
Corno A, Kocica M, TorrentGuasp F: The helical ventricular myocardial band of TorrentGuasp: potential implications in congenital heart defects. Eur J Cardiothorac Surg 2006, 29(Supplement 1):S61S68.
Buckberg GD: Architecture must document functional evidence to explain the living rhythm. Eur J CardioThorac Surg 2005, 27(2):202–209. 10.1016/j.ejcts.2004.10.052
Buckberg G, Myron L, Ballester M, Beyar R, Burkhoff D, Coghlan H, Doyle M, Epstein N, Gharib M, Ideker R, Ingels N, LeWinter M, McCulloch A, Pohost G, Reinlib L, Sahn D, Sopko G, Spinale F, Spotnitz H, TorrentGuasp F, Shapiro E: Left ventricular form and function : Scientific priorities and strategic planning for development of new views of disease. Circulation 2004, 110(14):e333e336. 10.1161/01.CIR.0000143625.56882.5C
Streeter D: Handbook of Physiology. Sec. 2. Vol. I. The Heart.. Bethesda, Maryland: Am Physiol Soc. Chap. Gross morphology and fiber geometry of the heart; 1979:. 61–112
Nielsen P, LeGrice I, Smaill B, Hunter P: Mathematical model of the geometry and fibrous structure of the heart. Am J Physiol 1991, 260: H1365H1378.
Pettigrew J: On the arrangement of the muscular fibers of the ventricular portion of the heart of the mammal. Proc Roy Soc, London 1860, 10: 433–440.
Helm P, Tseng H, Younes L, McVeigh E, Winslow R: Ex Vivo 3D diffusion tensor imaging and quantification of cardiac laminar structure. Magn Reson Med 2005, 54: 850–859. 10.1002/mrm.20622
Lunkenheimer PP, Redmann K, Kling N, Jiang X, Rothaus K, Cryer CW, Wübbeling F, Niederer P, Heitz PU, Yen Ho S, Anderson RH: Threedimensional architecture of the left ventricular myocardium. Anat Rec Part A 2006, 288: 565–578.
Helm P, Beg M, Miller M, Winslow R: Measuring and mapping cardiac fiber and laminar architecture using diffusion tensor MR imaging. Ann N Y Acad Sci 2005, 1047: 296–307. 10.1196/annals.1341.026
Vadakkumpadan F, Arevalo H, Prassl AJ, Chen J, Kickinger F, Kohl P, Plank G, Trayanova N: Imagebased models of cardiac structure in health and disease. Wiley Interdiscip Rev Syst Biol Med 2010, 2(4):489–506. 10.1002/wsbm.76
Gilbert S, Sands G, LeGrice I, Smaill B, Bernus O, Trew M: A framework for myoarchitecture analysis of high resolution cardiac MRI and comparison with Diffusion Tensor MRI. In Engineering in Medicine and Biology Society (EMBC), 2012 Annual International Conference of the IEEE.. : ; 2012:4063–4066.
Aslanidi O, Nikolaidou T, Zhao J, Smaill B, Gilbert S, Holden A, Lowe T, Withers P, Stephenson R, Jarvis J, Hancox J, Boyett M, Zhang H: Application of microcomputed tomography with iodine staining to cardiac imaging, segmentation, and computational model development. Med Imaging, IEEE Trans 2013, 32: 8–17.
Trew M, Caldwell B, Sands G, LeGrice I, Smaill B: Threedimensional cardiac tissue image registration for analysis of in vivo electrical mapping. Ann Biomed Eng 2011, 39: 235–248. 10.1007/s1043901001637
Jouk P, Usson Y, Michalowicz G, Grossi L: Threedimensional cartography of the pattern of the myofibres in the second trimester fetal human heart. Anat Embryol (Berl) 2000, 202(2):103–118. 10.1007/s004290000103
Vicky YW: Modelling In Vivo Cardiac Mechanics using MRI and FEM. PhD thesis, Auckland Bioengineering Institute, The University of Auckland, New Zealand; 2012
Gurev V, Lee T, Constantino J, Arevalo H, Trayanova N: Models of cardiac electromechanics based on individual hearts imaging data: Imagebased electromechanical models of the heart. Biomech Model Mechanobiol 2011, 10: 295–306. 10.1007/s1023701002355
Zhang Y, Liang X, Ma J, Jing Y, Gonzales MJ, Villongco C, Krishnamurthy A, Frank LR, Nigam V, Stark P, Narayan SM, McCulloch AD: An atlasbased geometry pipeline for cardiac Hermite model construction and diffusion tensor reorientation. Med Image Anal 2012, 16(6):1130–1141. 10.1016/j.media.2012.06.005
Seemann G: Modeling of electrophysiology and tension development in the human heart. PhD thesis, Universitat Karlsruhe; 2005
Bishop MJ, Plank G, Burton RAB, Schneider JE, Gavaghan DJ, Grau V, Kohl P: Development of an anatomically detailed MRIderived rabbit ventricular model and assessment of its impact on simulations of electrophysiological function. Am J Physiol Heart Circ Physiol 2010, 298(2):H699H718. 10.1152/ajpheart.00606.2009
Bayer J, Blake R, Plank G, Trayanova N: A novel rulebased algorithm for assigning myocardial fiber orientation to computational heart models. Ann Biomed Eng 2012, 40(10):2243–2254. 10.1007/s1043901205935
Hren R: A realistic model of the human ventricular myocardium: application to the study of ectopic activation. PhD thesis, Halifax, Nova Scotia, Canada: Dalhousie University; 1996
Peskin C: Fiber architecture of the left ventricular wall: An asymptotic analysis. Commun Pure Appl Math 1989, 42: 79–113. 10.1002/cpa.3160420106
Chadwick R: Mechanics of the left ventricle. Biophys J 1982, 39: 279–288. 10.1016/S00063495(82)845189
Bovendeerd P, Arts T, Huyghe J, Campen DV, Reneman R: Dependence of local left ventricular wall mechanics on myocardial fiber orientation: a model study. J Biomechanics 1992, 25(10):1129–1140. 10.1016/00219290(92)90069D
Arts T, Veenstra P, Reneman R: Epicardial deformation and left ventricular wall mechanics during ejection in the dog. Am J Physiol Heart Circ Physiol 1982, 243(12):H379H390.
Sinha S, Stein KM, Christini DJ: Critical role of inhomogeneities in pacing termination of cardiac reentry. Chaos: An Interdiscip J Nonlinear Sci 2002, 12(3):893–902. 10.1063/1.1501176
Acknowledgements
This work was supported by the Presidium of Urals Branch of the Russian Academy of Sciences (project 12M142009), the Flemish Community of Belgium (grant 1F2B8M/JDW/20102011/10BTLRUS01), Ghent University (grant 01SF1511), Russian Foundation for Basic Research (grant 130196048) and the Government of Sverdlovsk Region. The authors thank Hans Dierckx for his help in the preparation of the article.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
SFP did the mathematical and programming work and wrote a draft of the main part of the article. VIB took part in the mathematical conception formulation and discussion of the results. AVP participated at a later stage of the research on model verification, proposed the idea on epi and endocardial anisotropy fit and took part in writing the paper. LBK formulated the basic ideas of the mathematical construction, took part in writing Abstract, Introduction, Discussion, and Conclusions, editing the paper and discussion of the results. OES took part in discussion of the mathematical construction, proposed the ideas for numerical algorithms, experiments and visualisation, took part in discussion of the results and editing the paper. VSM took part in the physiological conception formulation, comparison of the results with the real morphological data and discussion of the model validity. All authors read and approved the final manuscript.
Electronic supplementary material
12938_2013_675_MOESM1_ESM.avi
Additional file 1: Rotating plane curve in space. On the left, changing curve on plane; on the right, the same curve moves and rotates in space. (AVI 2 MB)
12938_2013_675_MOESM2_ESM.avi
Additional file 2: Traced rotating plane curve in space. The same movement of the same curve, but with its trace, which is a spiral surface. Colour is linked with height z, red codes 0, purple codes 60 mm. (AVI 1 MB)
12938_2013_675_MOESM3_ESM.avi
Additional file 3: Lower part of the spiral surface. It is the lower part (0 ≤ z ≤ 15 mm) of the forming spiral surface from Additional file 2. Here we see the apical zone clearer. Colour is linked with height z, red codes 0, purple codes 15 mm. (AVI 1 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Pravdin, S.F., Berdyshev, V.I., Panfilov, A.V. et al. Mathematical model of the anatomy and fibre orientation field of the left ventricle of the heart. BioMed Eng OnLine 12, 54 (2013). https://doi.org/10.1186/1475925X1254
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1475925X1254