A publicly available virtual cohort of four-chamber heart meshes for cardiac electro-mechanics simulations

Motivation: Computational models of the heart are increasingly being used in the development of devices, patient diagnosis and therapy guidance. While software techniques have been developed for simulating single hearts, there remain significant challenges in simulating cohorts of virtual hearts from multiple patients.

Dataset Description: We present the first database of four-chamber heart models suitable for electro-mechanical simulations. Our database consists of twenty-four four-chamber heart models generated from end-diastolic CT acquired from heart failure patients recruited for cardiac resynchronization therapy upgrade. We also provide a higher resolution version for each of the twenty-four 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.1mm. The elements of all the twenty-four 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) Left atrium appendage ring 8) Left superior pulmonary vein ring 9) Left inferior pulmonary vein ring 10) Right inferior pulmonary vein ring 11) Right superior pulmonary vein ring 12) Superior vena cava ring 13) Inferior vena cava ring 14) Mitral valve plane 15) Tricuspid valve plane 16) Aortic valve plane 17) Pulmonary valve plane 18) Left atrial appendage valve plane 19) Left superior pulmonary vein valve plane 20) Left inferior pulmonary vein valve plane 21) Right inferior pulmonary vein valve plane 22) Right superior pulmonary vein valve plane 23) Superior vena cava valve plane 24) Inferior vena cava valve plane

Meshes 12 and 14 had only one left pulmonary vein and therefore had only 22 labels. 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 varying continuously from 0 at the apex to 1 at the base; a transmural coordinate varying continuously from 0 at the endocardium to 1 at the epicardium; a rotational coordinate 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 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 -100.

We also refined each mesh from 1.1mm resolution down to 0.39mm resolution. Each refined mesh has tags defined on its elements (same numbering as described above) and ventricular fibres. Suitability of 1mm-resolution meshes for electro-mechanics simulations were tested using test simulations for ventricular electrical activation and free mechanical contraction. We built the first cohort of twenty-four four-chamber heart meshes from heart failure patients CT data. These geometries can be used for large cohort computational studies and as a data set for verification of numerical methods for cardiac electro-mechanics.