Dataset Open Access

# Virtual cohort of adult healthy four-chamber heart meshes from CT images

Rodero, Cristobal; Strocchi, Marina; Marciniak, Maciej; Longobardi, Stefano; Whitaker, John; O'Neill, Mark D.; Gillette, Karli; Augustin, Christoph; Plank, Gernot; Vigmond, Edward J.; Lamata, Pablo; Niederer, Steven A.

Dataset Description: We present the first database of four-chamber healthy heart models suitable for electro-mechanical (EM) simulations. Our database consists of twenty four-chamber heart models generated from end-diastolic CT acquired from patients who went to the emergency room with acute chest pains. Since no cardiac conditions were detected in follow-up, these patients were taken as representative of "healthy" (or asymptomatic) hearts. These meshes were used for EM simulations and to build a statistical shape model (SSM). The output of the simulations and the weights of the SSM are also provided.

Cardiac meshes: We segmented end-diastolic CT. The segmentation was then upsampled and smoothed. The final multi-label segmentation was used to generate a tetrahedral mesh. The resulting meshes had an average edge length of 1 mm. The elements of all the twenty meshes are labelled as follows:

1. Left ventricle myocardium
2. Right ventricle myocardium
3. Left atrium myocardium
4. Right atrium myocardium
5. Aorta wall
6. Pulmonary artery wall
7. Mitral valve plane
8. Tricuspid valve plane
9. Aortic valve plane
10. Pulmonary valve plane
11. Left atrium appendage "inlet"
12. Left superior pulmonary vein inlet
13. Left inferior pulmonary vein inlet
14. Right inferior pulmonary vein inlet
15. Right superior pulmonary vein inlet
16. Superior vena cava inlet
17. Inferior vena cava inlet
18. Left atrial appendage border
19. Right inferior pulmonary vein border
20. Left inferior pulmonary vein border
21. Left superior pulmonary vein border
22. Right superior pulmonary vein border
23. Superior vena cava border
24. Inferior vena cava border

Ventricular fibres were generated using a rule-based method, with a fibre orientation varying transmurally from endocardium to epicardium from 80˚ to -60˚, respectively.  We defined a system of universal ventricular coordinates on the meshes: an apico-basal coordinate (Z) varying continuously from 0 at the apex to 1 at the base (defined with the mitral and tricuspid valve); a transmural coordinate ($$\rho$$) varying continuously from 0 at the endocardium to 1 at the epicardium; a rotational coordinate ($$\phi$$) varying continuously from – π at the left ventricular free wall, 0 at the septum and then back to + π at the left ventricular free wall; intra-ventricular coordinate (V) defined at -1 at the left ventricle and +1 at the right ventricle. This coordinate system was assigned to the ventricles in the four-chamber meshes and all the other labels were assigned with -10.

We provide a zipped folder for each mesh, A VTK file for each mesh was included (in ASCII) as an UNSTRUCTURED GRID. In all the cases the following fields were included:

• POINTS, with the coordinates of the points in mm.
• CELL_TYPES, having all of the points the value 10 since they are tetrahedra.
• CELLS, with the indices of the vertices of every element.
• CELL_DATA, corresponding to the meshing tags.
• VECTORS, with the directions of the fibres.
• POINT_DATA, with four LOOKUP_TABLE subfields corresponding to the UVC in the order $$\rho$$, $$\phi$$, Z and V.

Cardiac simulations: For the cardiac EM simulations we used CARP (Cardiac Arrhythmia Research Package). We used the reaction-eikonal model for electrophysiology, stimulating as initial condition the bottom third (Z < 0.33) of the endocardium. We simulated the large deformation mechanics in a Lagrangian reference system. The ventricular myocardium was modelled as a hyperelastic transversely isotropic material with Guccione's strain energy function. The remaining tissues were modelled as non-contracting neo-Hookean materials. Simulations of meshes #09 and #10 failed to converge. Details on the specific parametrisation can be found in the supplements of the reference paper.

We provide comma-separated-values files with the output of the simulations used in the reference paper for validation purposes. The simulations of the cases that did not converge were not included. The acronyms used in the names of columns are:

• EDP: End-diastolic pressure
• EDV: End-diastolic volume
• Myo_vol: Myocardial volume of the ventricle (as sum of its elements)
• ESV: End-systolic volume
• SV: Stroke volume
• EF: Ejection fraction
• V1: Volume at time of peak flow
• EF1: First-Phase Ejection Fraction
• ESP: End-systolic pressure
• dPdtmax: Maximum increase of pressure
• dPdtmin: Maximum decrease of pressure
• PeakP: Peak pressure
• tpeak: Time to peak pressure
• ET: Ejection time
• ICT: Isovolumic contraction time
• IRT: Isovolumic relaxation time
• tsys: Duration of systole
• QRS: QRS duration
• AT1090: Time taken to activate from 10% to 90% of the mesh
• AT: Activation time of the left ventricle

Besides the output value name, in each column is specified the ventricle where that output was extracted from with the suffixes "_LV" or "_RV".

Statistical shape model: All the meshes but #20 were used to build a statistical shape model of four-chambers cardiac meshes. In short, we rigidly aligned the meshes and extracted the surfaces, representing them as deRham currents. The registration between meshes and computation of the average shape (also called atlas or template) was done using a Large Deformation Diffeomorphic Metric Mapping method. Each one of the meshes can be approximated as a linear combination of the shape modes, extracted using Principal Component Analysis on the space where the meshes are located. More details on the Statistical Shape Model are provided in the supplement of the reference paper. The average heart and extreme cases are provided in the repository named "Virtual cohort of extreme and average four-chamber heart meshes from statistical shape model". We have added 1000 more meshes from the same statistical shape model, modifying the weights from the PCA randomly within the 2SD range. These meshes are provided in the repository named "Virtual cohort of 1000 synthetic heart meshes from the adult human healthy population".

We provide the weights of the modes for each of the 19 meshes in a comma-separated-values file.

GP received support from the Austrian Science Fund (FWF) (https://fwf.ac.at/en/); grant number PI2760-B30. PL is supported by BHF [PG/16/75/32383]. SAN is supported by NIH R01-HL152256, ERC PREDICT-HF 453 (864055) and BHF (RG/20/4/34803).
Files (1.5 GB)
Name Size
01.tar.gz
md5:12df3f7f99fa8d2784d87f32ea98de71
91.1 MB
02.tar.gz
md5:50772706775d3e2a1ea90a984e799ac2
103.2 MB
03.tar.gz
md5:068522705fbd325c82b8694de02c3c74
85.9 MB
04.tar.gz
md5:02e3f9b60f99d85f3372ce903efc2ab0
102.6 MB
05.tar.gz
md5:1457ccc773e26b42ef2dafe399e6c2e1
54.9 MB
06.tar.gz
md5:a81424b8955b418e70b2909fcdf5a7bd
63.8 MB
07.tar.gz
md5:4853411c3ebea9271ac35b8385137575
61.4 MB
08.tar.gz
md5:4b8f883766e3387b1e8e65ef975e35d0
80.0 MB
09.tar.gz
md5:630c6322918571b9e930aa43430fe017
74.2 MB
10.tar.gz
md5:27651d985f2967042b2cebf993345264
96.1 MB
11.tar.gz
md5:22ac3f1246e7d13af9becb4e8eff0bd3
79.6 MB
12.tar.gz
md5:696769149674e27f9d33fc0979b75aca
59.5 MB
13.tar.gz
md5:8c0ea2f240e8c738960c04e994bb8923
84.9 MB
14.tar.gz
md5:6dcffc6e10e7f015d0ef1fe41561f8fd
94.7 MB
15.tar.gz
md5:d734e23870e0e4cc85285467f6357877
61.9 MB
16.tar.gz
md5:7ac94827d4e0ffa4bd60b2a1b270b2ca
74.5 MB
17.tar.gz
md5:f5eeefca986d8b4849ab263f365ff01e
70.1 MB
18.tar.gz
md5:004856c5ebb3f150a87147370b348d6c
53.7 MB
19.tar.gz
md5:73ebd0b20ac2c53fa5a4430e180693b5
70.4 MB
20.tar.gz
md5:6dbfa0c7ec88d045fd0f38c997fb1667
69.6 MB
modes_weights.csv
md5:51690579387fa6dc3435e172f29e6a43
4.4 kB
simulation_output.csv
4.0 kB
• Rodero et al. (2021), "Linking statistical shape models and simulated function in the healthy adult human heart". DOI: 10.1371/journal.pcbi.1008851

1,339
1,121
views