Medical imaging studies often rely on a single sample per subject, assuming it is representative of their physiological traits. However, variations in how input descriptors are defined or computed (e.g. due to a lack of consensus in the scientific field) may have a crucial impact on the analysis, and are hardly considered in practice. In this paper, we propose an original strategy based on representation learning to estimate a parametric map reflecting the impact of such definitional differences on a given physiological descriptor, previously extracted from medical images. We consider the different definitions or computations of such physiological descriptors as different high-dimensional data, potentially of heterogeneous types. We specifically focus on myocardial deformation (strain), for which there is limited agreement on its definition. We first use manifold alignment to match the latent representations associated with the different definitions of this descriptor. Then, we formulate plausible distributions in the latent space to represent definitional divergence across descriptors, from which we reconstruct a high-dimensional parametric map to visualize such definitional divergence. Due to the lack of proper ground truth for this specific clinical application, we first demonstrate this methodology on toy experiments and then expand the evaluation on right ventricular strain data from subjects obtained from 3D echocardiographic image sequences, for which different types of strain are available at each point of the right ventricle endocardial surface mesh. Beyond this illustrative application, our methodology has the potential to be generalised to many other population analyses considering heterogeneous high-dimensional descriptors.

Coarctation of the aorta (CoA) is a common congenital heart defect characterized by aortic narrowing. Prenatally, it has mild hemodynamic effects as right ventricular disproportion and ductus arteriosus (DA) dilation occur as adaptive mechanisms, but their impact on CoA hemodynamics remains poorly understood. To investigate this, we built a closed 0D computational model of fetal circulation and simulated different CoA cardiovascular remodeling patterns, including aortic isthmus (AoI) narrowing, ventricular disproportion, and DA dilation. Our results showed mild AoI narrowing (80% of reference diameter) required up to 1.7 right/left ventricular end-diastolic volume ratio and 115% DA dilation to maintain physiological pressures, wall shear stresses, and organ perfusion. In contrast, severe narrowing (20% of reference AoI diameter) required up to 5 right/left ventricular end-diastolic volume ratio and 125% DA dilation, highlighting the necessity of co-occurrence of prenatal ventricular disproportion and DA dilation to compensate for AoI narrowing. These physiological regions were validated with ultrasonographic measurements from 7 controls and 9 CoA patients. We compared blood pressures, velocities, and volumetric flow rates across different fetoplacental anatomical sites. AoI velocity showed a delayed retrograde flow peak and increased antegrade diastolic velocity with greater AoI narrowing, which may aid in diagnosing CoA. Minimal differences were observed in other velocities and pressures. Volumetric flow rates across varying degrees of AoI narrowing decreased in the AoI and mitral and aortic valves, remained stable in the middle cerebral and umbilical arteries, and increased in the DA and tricuspid and pulmonary valves. Therefore, we corroborated that in fetal CoA a redistribution of blood flow occurs to ensure perfusion of the brain and placenta, without a significant alteration in fetal hemodynamics (blood pressure and velocities) except for increased diastolic velocities in the AoI.

Cardiac imaging protocols usually result in several types of acquisitions and descriptors extracted from the images. The statistical analysis of such data across a population may be challenging, and can be addressed by fusion techniques within a dimensionality reduction framework. However, directly combining different data types may lead to unfair comparisons (for heterogeneous descriptors) or over-exploitation of information (for strongly correlated modalities). In contrast, physicians progressively consider each type of data based on hierarchies derived from their experience or evidence-based recommendations, an inspiring approach for data fusion strategies. In this paper, we propose a novel methodology for hierarchical data fusion and unsupervised representation learning. It mimics the physicians' approach by progressively integrating different high-dimensional data descriptors according to a known hierarchy. We model this hierarchy with a Hierarchical Gaussian Process Latent Variable Model (GP-LVM), which links the estimated low-dimensional latent representation and high-dimensional observations at each level in the hierarchy, with additional links between consecutive levels of the hierarchy. We demonstrate the relevance of this approach on a dataset of 1726 magnetic resonance image slices from 123 patients revascularized after acute myocardial infarction (MI) (first level in the hierarchy), some of them undergoing reperfusion injury (microvascular obstruction (MVO), second level in the hierarchy). Our experiments demonstrate that our hierarchical model provides consistent data organization across levels of the hierarchy and according to physiological characteristics of the lesions. This allows more relevant statistical analysis of myocardial lesion patterns, and in particular subtle lesions such as MVO.

3D echocardiography is an increasingly popular tool for assessing cardiac remodelling in the right ventricle (RV). It allows quantification of the cardiac chambers without any geometric assumptions, which is the main weakness of 2D echocardiography. However, regional quantification of geometry and function is limited by the lower spatial and temporal resolution and the scarcity of identifiable anatomical landmarks, especially within the ventricular cavity. We developed a technique for regionally assessing the volume of 3 relevant RV volumetric regions: apical, inlet and outflow. The proposed parcellation method is based on the geodesic distances to anatomical landmarks that are easily identifiable in the images: the apex and the tricuspid and pulmonary valves, each associated to a region. Based on these distances, we define a partition in the endocardium at end-diastole (ED). This partition is then interpolated to the blood cavity using the Laplace equation, which allows to compute regional volumes. For obtaining an end-systole (ES) partition, the endocardial partition is transported from ED to ES using a commercial image-based tracking software, and then the interpolation process is repeated. We assessed the intra- and inter-observer reproducibility using a 10-subjects dataset containing repeated quantifications of the same images, obtaining intra- and inter- observer errors (7-12% and 10-23% respectively). Finally, we propose a novel synthetic mesh generation algorithm that deforms a template mesh imposing a user-defined strain to a template mesh. We used this method to create a new dataset for involving distinct types of remodelling that were used to assess the sensitivity of the parcellation method to identify volume changes affecting different parts. We show that the parcellation method is adequate for capturing local circumferential and global circumferential and longitudinal RV remodelling, which are the most clinically relevant cases.

Two-dimensional representation of 3D anatomical structures is a simple and intuitive way for analysing patient information across populations and image modalities. It also allows convenient visualizations that can be included in clinical reports for a fast overview of the whole structure. While cardiac ventricles, especially the left ventricle, have an established standard representation (e.g. bull's eye plot), the 2D depiction of the left atrium (LA) is challenging due to its sub-structural complexity including the pulmonary veins (PV) and the left atrial appendage (LAA). Quasi-conformal flattening techniques, successfully applied to cardiac ventricles, require additional constraints in the case of the LA to place the PV and LAA in the same geometrical 2D location for different cases. Some registration-based methods have been proposed but 3D (or 2D) surface registration is time-consuming and prone to errors if the geometries are very different. We propose a novel atrial flattening methodology where a quasi-conformal 2D map of the LA is obtained quickly and without errors related to registration. In our approach, the LA is divided into 5 regions which are then mapped to their analogue two-dimensional regions. A dataset of 67 human left atria from magnetic resonance images (MRI) was studied to derive a population-based 2D LA template representing the averaged relative locations of the PVs and LAA. The clinical application of the proposed methodology is illustrated on different use cases including the integration of MRI and electroanatomical data.