Deep learning cardiac motion analysis for human survival prediction

Motion analysis is used in computer vision to understand the behaviour of moving objects in sequences of images. Optimising the interpretation of dynamic biological systems requires accurate and precise motion tracking as well as efficient representations of high-dimensional motion trajectories so that these can be used for prediction tasks. Here we use image sequences of the heart, acquired using cardiac magnetic resonance imaging, to create time-resolved three-dimensional segmentations using a fully convolutional network trained on anatomical shape priors. This dense motion model formed the input to a supervised denoising autoencoder (4Dsurvival), which is a hybrid network consisting of an autoencoder that learns a task-specific latent code representation trained on observed outcome data, yielding a latent representation optimised for survival prediction. To handle right-censored survival outcomes, our network used a Cox partial likelihood loss function. In a study of 302 patients the predictive accuracy (quantified by Harrell's C-index) was significantly higher (p < .0001) for our model C=0.73 (95$\%$ CI: 0.68 - 0.78) than the human benchmark of C=0.59 (95$\%$ CI: 0.53 - 0.65). This work demonstrates how a complex computer vision task using high-dimensional medical image data can efficiently predict human survival

Image-Based Cardiac Diagnosis With Machine Learning: A Review

Cardiac imaging plays an important role in the diagnosis of cardiovascular disease (CVD). Until now, its role has been limited to visual and quantitative assessment of cardiac structure and function. However, with the advent of big data and machine learning, new opportunities are emerging to build artificial intelligence tools that will directly assist the clinician in the diagnosis of CVDs. This paper presents a thorough review of recent works in this field and provide the reader with a detailed presentation of the machine learning methods that can be further exploited to enable more automated, precise and early diagnosis of most CVDs

Texture analysis of cardiac cine magnetic resonance imaging to detect nonviable segments in patients with chronic myocardial infarction

[EN] Purpose: To investigate the ability of texture analysis to differentiate between infarcted nonviable, viable, and remote segments on cardiac cine magnetic resonance imaging (MRI). Methods: This retrospective study included 50 patients suffering chronic myocardial infarction. The data were randomly split into training (30 patients) and testing (20 patients) sets. The left ventricular myocardium was segmented according to the 17-segment model in both cine and late gadolinium enhancement (LGE) MRI. Infarcted myocardium regions were identified on LGE in short-axis views. Nonviable segments were identified as those showing LGE 50%, and viable segments those showing 0 < LGE < 50% transmural extension. Features derived from five texture analysis methods were extracted from the segments on cine images. A support vector machine (SVM) classifier was trained with different combination of texture features to obtain a model that provided optimal classification performance. Results: The best classification on testing set was achieved with local binary patterns features using a 2D + t approach, in which the features are computed by including information of the time dimension available in cine sequences. The best overall area under the receiver operating characteristic curve (AUC) were: 0.849, sensitivity of 92% to detect nonviable segments, 72% to detect viable segments, and 85% to detect remote segments. Conclusion: Nonviable segments can be detected on cine MRI using texture analysis and this may be used as hypothesis for future research aiming to detect the infarcted myocardium by means of a gadolinium-free approach.This work was supported in part by the Spanish Ministerio de Economia y Competitividad (MINECO) and FEDER funds under grant BFU2015-64380-C2-2-R, by Instituto de Salud Carlos III and FEDER funds under grants FIS PI14/00271 and PIE15/00013 and by the Generalitat Valenciana under grant PROMETEO/2013/007. The first author, Andres Larroza, was supported by grant FPU12/01140 from the Spanish Ministerio de Educacion, Cultura y Deporte (MECD).Larroza, A.; López-Lereu, M.; Monmeneu, J.; Gavara-Doñate, J.; Chorro, F.; Bodi, V.; Moratal, D. (2018). Texture analysis of cardiac cine magnetic resonance imaging to detect nonviable segments in patients with chronic myocardial infarction. Medical Physics. 45(4):1471-1480. https://doi.org/10.1002/mp.12783S14711480454Castellano, G., Bonilha, L., Li, L. M., & Cendes, F. (2004). Texture analysis of medical images. Clinical Radiology, 59(12), 1061-1069. doi:10.1016/j.crad.2004.07.008Hodgdon, T., McInnes, M. D. F., Schieda, N., Flood, T. A., Lamb, L., & Thornhill, R. E. (2015). Can Quantitative CT Texture Analysis be Used to Differentiate Fat-poor Renal Angiomyolipoma from Renal Cell Carcinoma on Unenhanced CT Images? Radiology, 276(3), 787-796. doi:10.1148/radiol.2015142215Larroza, A., Moratal, D., Paredes-Sánchez, A., Soria-Olivas, E., Chust, M. L., Arribas, L. A., & Arana, E. (2015). Support vector machine classification of brain metastasis and radiation necrosis based on texture analysis in MRI. Journal of Magnetic Resonance Imaging, 42(5), 1362-1368. doi:10.1002/jmri.24913Thevenot, J., Hirvasniemi, J., Pulkkinen, P., Määttä, M., Korpelainen, R., Saarakkala, S., & Jämsä, T. (2014). Assessment of Risk of Femoral Neck Fracture with Radiographic Texture Parameters: A Retrospective Study. Radiology, 272(1), 184-191. doi:10.1148/radiol.14131390Kassner, A., & Thornhill, R. E. (2010). Texture Analysis: A Review of Neurologic MR Imaging Applications. American Journal of Neuroradiology, 31(5), 809-816. doi:10.3174/ajnr.a2061Pfeiffer, M. P., & Biederman, R. W. W. (2015). Cardiac MRI. Medical Clinics of North America, 99(4), 849-861. doi:10.1016/j.mcna.2015.02.011Flett, A. S., Hasleton, J., Cook, C., Hausenloy, D., Quarta, G., Ariti, C., … Moon, J. C. (2011). Evaluation of Techniques for the Quantification of Myocardial Scar of Differing Etiology Using Cardiac Magnetic Resonance. JACC: Cardiovascular Imaging, 4(2), 150-156. doi:10.1016/j.jcmg.2010.11.015Engan K Eftestøl T Ørn S Kvaloy JT Woie L Exploratory data analysis of image texture and statistical features on myocardium and infarction areas in cardiac magnetic resonance images 2010Kotu LP Engan K Eftestøl T Ørn S Woie L Segmentation of scarred and non-scarred myocardium in LG enhanced CMR images using intensity-based textural analysis 2011Kotu, L., Engan, K., Skretting, K., Måløy, F., Ørn, S., Woie, L., & Eftestøl, T. (2013). Probability mapping of scarred myocardium using texture and intensity features in CMR images. BioMedical Engineering OnLine, 12(1), 91. doi:10.1186/1475-925x-12-91Schofield, R., Ganeshan, B., Kozor, R., Nasis, A., Endozo, R., Groves, A., … Moon, J. C. (2016). CMR myocardial texture analysis tracks different etiologies of left ventricular hypertrophy. Journal of Cardiovascular Magnetic Resonance, 18(S1). doi:10.1186/1532-429x-18-s1-o82Larroza, A., Materka, A., López-Lereu, M. P., Monmeneu, J. V., Bodí, V., & Moratal, D. (2017). Differentiation between acute and chronic myocardial infarction by means of texture analysis of late gadolinium enhancement and cine cardiac magnetic resonance imaging. European Journal of Radiology, 92, 78-83. doi:10.1016/j.ejrad.2017.04.024Baessler, B., Mannil, M., Oebel, S., Maintz, D., Alkadhi, H., & Manka, R. (2018). Subacute and Chronic Left Ventricular Myocardial Scar: Accuracy of Texture Analysis on Nonenhanced Cine MR Images. Radiology, 286(1), 103-112. doi:10.1148/radiol.2017170213Hervas, A., Ruiz-Sauri, A., de Dios, E., Forteza, M. J., Minana, G., Nunez, J., … Bodi, V. (2015). Inhomogeneity of collagen organization within the fibrotic scar after myocardial infarction: results in a swine model and in human samples. Journal of Anatomy, 228(1), 47-58. doi:10.1111/joa.12395Heiberg, E., Sjögren, J., Ugander, M., Carlsson, M., Engblom, H., & Arheden, H. (2010). Design and validation of Segment - freely available software for cardiovascular image analysis. BMC Medical Imaging, 10(1). doi:10.1186/1471-2342-10-1Bodí, V., Sanchis, J., López-Lereu, M. P., Losada, A., Núñez, J., Pellicer, M., … Llácer, À. (2005). Usefulness of a Comprehensive Cardiovascular Magnetic Resonance Imaging Assessment for Predicting Recovery of Left Ventricular Wall Motion in the Setting of Myocardial Stunning. Journal of the American College of Cardiology, 46(9), 1747-1752. doi:10.1016/j.jacc.2005.07.039Rangayyan, R. M., Nguyen, T. M., Ayres, F. J., & Nandi, A. K. (2009). Effect of Pixel Resolution on Texture Features of Breast Masses in Mammograms. Journal of Digital Imaging, 23(5), 547-553. doi:10.1007/s10278-009-9238-0Materka A Strzelecki M On the importance of MRI nonuniformity correction for texture analysis 2013Collewet, G., Strzelecki, M., & Mariette, F. (2004). Influence of MRI acquisition protocols and image intensity normalization methods on texture classification. Magnetic Resonance Imaging, 22(1), 81-91. doi:10.1016/j.mri.2003.09.001Vallières M MATLAB programming tools for radiomics analysis https://github.com/mvallieres/radiomicsZhao G Pietikainen M Center for machine vision and signal analysis http://www.cse.oulu.fi/CMV/Downloads/LBPMatlabZwanenburg A Leger S Vallières M Löck S Image biomarker standardisation initiative 2017 http://arxiv.org/abs/1612.07003Vallières, M., Freeman, C. R., Skamene, S. R., & El Naqa, I. (2015). A radiomics model from joint FDG-PET and MRI texture features for the prediction of lung metastases in soft-tissue sarcomas of the extremities. Physics in Medicine and Biology, 60(14), 5471-5496. doi:10.1088/0031-9155/60/14/5471Zhao, G., & Pietikainen, M. (2007). Dynamic Texture Recognition Using Local Binary Patterns with an Application to Facial Expressions. IEEE Transactions on Pattern Analysis and Machine Intelligence, 29(6), 915-928. doi:10.1109/tpami.2007.1110Ojala T Pietikäinen M Mäenpää T A generalized local binary pattern operator for multiresolution gray scale and rotation invariant texture classificationDuan, K.-B., Rajapakse, J. C., Wang, H., & Azuaje, F. (2005). Multiple SVM-RFE for Gene Selection in Cancer Classification With Expression Data. IEEE Transactions on Nanobioscience, 4(3), 228-234. doi:10.1109/tnb.2005.853657Guyon, I., Weston, J., Barnhill, S., & Vapnik, V. (2002). Machine Learning, 46(1/3), 389-422. doi:10.1023/a:1012487302797Wang, S., & Summers, R. M. (2012). Machine learning and radiology. Medical Image Analysis, 16(5), 933-951. doi:10.1016/j.media.2012.02.005Kuhn, M. (2008). Building Predictive Models inRUsing thecaretPackage. Journal of Statistical Software, 28(5). doi:10.18637/jss.v028.i05Colby J (multiple) Support Vector Machine Recursive Feature Elimination - mSVM-RFE http://www.colbyimaging.com/wiki/statistics/msvm-rfeSalzberg, S. L. (1997). Data Mining and Knowledge Discovery, 1(3), 317-328. doi:10.1023/a:1009752403260Bodí, V., Husser, O., Sanchis, J., Núñez, J., López-Lereu, M. P., Monmeneu, J. V., … Llácer, A. (2010). Contractile Reserve and Extent of Transmural Necrosis in the Setting of Myocardial Stunning: Comparison at Cardiac MR Imaging. Radiology, 255(3), 755-763. doi:10.1148/radiol.10091191Bodi, V., Monmeneu, J. V., Ortiz-Perez, J. T., Lopez-Lereu, M. P., Bonanad, C., Husser, O., … Chorro, F. J. (2016). Prediction of Reverse Remodeling at Cardiac MR Imaging Soon after First ST-Segment–Elevation Myocardial Infarction: Results of a Large Prospective Registry. Radiology, 278(1), 54-63. doi:10.1148/radiol.2015142674Shriki, J. E., Surti, K. S., Farvid, A. F., Lee, C. C., Samadi, S., Hirschbeinv, J., & Colletti, P. M. (2011). Chemical Shift Artifact on Steady-State Free Precession Cardiac Magnetic Resonance Sequences as a Result of Lipomatous Metaplasia: A Novel Finding in Chronic Myocardial Infarctions. Canadian Journal of Cardiology, 27(5), 664.e17-664.e23. doi:10.1016/j.cjca.2010.12.074Goldfarb, J. W., McLaughlin, J., Gray, C. A., & Han, J. (2011). Cyclic CINE-balanced steady-state free precession image intensity variations: Implications for the detection of myocardial edema. Journal of Magnetic Resonance Imaging, 33(3), 573-581. doi:10.1002/jmri.22368Gillies, R. J., Kinahan, P. E., & Hricak, H. (2016). Radiomics: Images Are More than Pictures, They Are Data. Radiology, 278(2), 563-577. doi:10.1148/radiol.201515116

Advances in machine learning applications for cardiovascular 4D flow MRI

Four-dimensional flow magnetic resonance imaging (MRI) has evolved as a non-invasive imaging technique to visualize and quantify blood flow in the heart and vessels. Hemodynamic parameters derived from 4D flow MRI, such as net flow and peak velocities, but also kinetic energy, turbulent kinetic energy, viscous energy loss, and wall shear stress have shown to be of diagnostic relevance for cardiovascular diseases. 4D flow MRI, however, has several limitations. Its long acquisition times and its limited spatio-temporal resolutions lead to inaccuracies in velocity measurements in small and low-flow vessels and near the vessel wall. Additionally, 4D flow MRI requires long post-processing times, since inaccuracies due to the measurement process need to be corrected for and parameter quantification requires 2D and 3D contour drawing. Several machine learning (ML) techniques have been proposed to overcome these limitations. Existing scan acceleration methods have been extended using ML for image reconstruction and ML based super-resolution methods have been used to assimilate high-resolution computational fluid dynamic simulations and 4D flow MRI, which leads to more realistic velocity results. ML efforts have also focused on the automation of other post-processing steps, by learning phase corrections and anti-aliasing. To automate contour drawing and 3D segmentation, networks such as the U-Net have been widely applied. This review summarizes the latest ML advances in 4D flow MRI with a focus on technical aspects and applications. It is divided into the current status of fast and accurate 4D flow MRI data generation, ML based post-processing tools for phase correction and vessel delineation and the statistical evaluation of blood flow

Medical image segmentation and analysis using statistical shape modelling and inter-landmark relationships

The study of anatomical morphology is of great importance to medical imaging, with applications varying from clinical diagnosis to computer-aided surgery. To this end, automated tools are required for accurate extraction of the anatomical boundaries from the image data and detailed interpretation of morphological information. This thesis introduces a novel approach to shape-based analysis of medical images based on Inter- Landmark Descriptors (ILDs). Unlike point coordinates that describe absolute position, these shape variables represent relative configuration of landmarks in the shape. The proposed work is motivated by the inherent difficulties of methods based on landmark coordinates in challenging applications. Through explicit invariance to pose parameters and decomposition of the global shape constraints, this work permits anatomical shape analysis that is resistant to image inhomogeneities and geometrical inconsistencies. Several algorithms are presented to tackle specific image segmentation and analysis problems, including automatic initialisation, optimal feature point search, outlier handling and dynamic abnormality localisation. Detailed validation results are provided based on various cardiovascular magnetic resonance datasets, showing increased robustness and accuracy.Open acces

Automated Segmentation of Left and Right Ventricles in MRI and Classification of the Myocarfium Abnormalities

A fundamental step in diagnosis of cardiovascular diseases, automated left and right ventricle (LV and RV) segmentation in cardiac magnetic resonance images (MRI) is still acknowledged to be a difficult problem. Although algorithms for LV segmentation do exist, they require either extensive training or intensive user inputs. RV segmentation in MRI has yet to be solved and is still acknowledged a completely unsolved problem because its shape is not symmetric and circular, its deformations are complex and varies extensively over the cardiac phases, and it includes papillary muscles. In this thesis, I investigate fast detection of the LV endo- and epi-cardium surfaces (3D) and contours (2D) in cardiac MRI via convex relaxation and distribution matching. A rapid 3D segmentation of the RV in cardiac MRI via distribution matching constraints on segment shape and appearance is also investigated. These algorithms only require a single subject for training and a very simple user input, which amounts to one click. The solution is sought following the optimization of functionals containing probability product kernel constraints on the distributions of intensity and geometric features. The formulations lead to challenging optimization problems, which are not directly amenable to convex-optimization techniques. For each functional, the problem is split into a sequence of sub-problems, each of which can be solved exactly and globally via a convex relaxation and the augmented Lagrangian method. Finally, an information-theoretic based artificial neural network (ANN) is proposed for normal/abnormal LV myocardium motion classification. Using the LV segmentation results, the LV cavity points is estimated via a Kalman filter and a recursive dynamic Bayesian filter. However, due to the similarities between the statistical information of normal and abnormal points, differentiating between distributions of abnormal and normal points is a challenging problem. The problem was investigated with a global measure based on the Shannon\u27s differential entropy (SDE) and further examined with two other information-theoretic criteria, one based on Renyi entropy and the other on Fisher information. Unlike the existing information-theoretic studies, the approach addresses explicitly the overlap between the distributions of normal and abnormal cases, thereby yielding a competitive performance. I further propose an algorithm based on a supervised 3-layer ANN to differentiate between the distributions farther. The ANN is trained and tested by five different information measures of radial distance and velocity for points on endocardial boundary

Developing advanced mathematical models for detecting abnormalities in 2D/3D medical structures.

Detecting abnormalities in two-dimensional (2D) and three-dimensional (3D) medical structures is among the most interesting and challenging research areas in the medical imaging field. Obtaining the desired accurate automated quantification of abnormalities in medical structures is still very challenging. This is due to a large and constantly growing number of different objects of interest and associated abnormalities, large variations of their appearances and shapes in images, different medical imaging modalities, and associated changes of signal homogeneity and noise for each object. The main objective of this dissertation is to address these problems and to provide proper mathematical models and techniques that are capable of analyzing low and high resolution medical data and providing an accurate, automated analysis of the abnormalities in medical structures in terms of their area/volume, shape, and associated abnormal functionality. This dissertation presents different preliminary mathematical models and techniques that are applied in three case studies: (i) detecting abnormal tissue in the left ventricle (LV) wall of the heart from delayed contrast-enhanced cardiac magnetic resonance images (MRI), (ii) detecting local cardiac diseases based on estimating the functional strain metric from cardiac cine MRI, and (iii) identifying the abnormalities in the corpus callosum (CC) brain structure—the largest fiber bundle that connects the two hemispheres in the brain—for subjects that suffer from developmental brain disorders. For detecting the abnormal tissue in the heart, a graph-cut mathematical optimization model with a cost function that accounts for the object’s visual appearance and shape is used to segment the the inner cavity. The model is further integrated with a geometric model (i.e., a fast marching level set model) to segment the outer border of the myocardial wall (the LV). Then the abnormal tissue in the myocardium wall (also called dead tissue, pathological tissue, or infarct area) is identified based on a joint Markov-Gibbs random field (MGRF) model of the image and its region (segmentation) map that accounts for the pixel intensities and the spatial interactions between the pixels. Experiments with real in-vivo data and comparative results with ground truth (identified by a radiologist) and other approaches showed that the proposed framework can accurately detect the pathological tissue and can provide useful metrics for radiologists and clinicians. To estimate the strain from cardiac cine MRI, a novel method based on tracking the LV wall geometry is proposed. To achieve this goal, a partial differential equation (PDE) method is applied to track the LV wall points by solving the Laplace equation between the LV contours of each two successive image frames over the cardiac cycle. The main advantage of the proposed tracking method over traditional texture-based methods is its ability to track the movement and rotation of the LV wall based on tracking the geometric features of the inner, mid-, and outer walls of the LV. This overcomes noise sources that come from scanner and heart motion. To identify the abnormalities in the CC from brain MRI, the CCs are aligned using a rigid registration model and are segmented using a shape-appearance model. Then, they are mapped to a simple unified space for analysis. This work introduces a novel cylindrical mapping model, which is conformal (i.e., one to one transformation and bijective), that enables accurate 3D shape analysis of the CC in the cylindrical domain. The framework can detect abnormalities in all divisions of the CC (i.e., splenium, rostrum, genu and body). In addition, it offers a whole 3D analysis of the CC abnormalities instead of only area-based analysis as done by previous groups. The initial classification results based on the centerline length and CC thickness suggest that the proposed CC shape analysis is a promising supplement to the current techniques for diagnosing dyslexia. The proposed techniques in this dissertation have been successfully tested on complex synthetic and MR images and can be used to advantage in many of today’s clinical applications of computer-assisted medical diagnostics and intervention

Automatic initialization and quality control of large-scale cardiac MRI segmentations

Continuous advances in imaging technologies enable ever more comprehensive phenotyping of human anatomy and physiology. Concomitant reduction of imaging costs has resulted in widespread use of imaging in large clinical trials and population imaging studies. Magnetic Resonance Imaging (MRI), in particular, offers one-stop-shop multidimensional biomarkers of cardiovascular physiology and pathology. A wide range of analysis methods offer sophisticated cardiac image assessment and quantification for clinical and research studies. However, most methods have only been evaluated on relatively small databases often not accessible for open and fair benchmarking. Consequently, published performance indices are not directly comparable across studies and their translation and scalability to large clinical trials or population imaging cohorts is uncertain. Most existing techniques still rely on considerable manual intervention for the initialization and quality control of the segmentation process, becoming prohibitive when dealing with thousands of images. The contributions of this paper are three-fold. First, we propose a fully automatic method for initializing cardiac MRI segmentation, by using image features and random forests regression to predict an initial position of the heart and key anatomical landmarks in an MRI volume. In processing a full imaging database, the technique predicts the optimal corrective displacements and positions in relation to the initial rough intersections of the long and short axis images. Second, we introduce for the first time a quality control measure capable of identifying incorrect cardiac segmentations with no visual assessment. The method uses statistical, pattern and fractal descriptors in a random forest classifier to detect failures to be corrected or removed from subsequent statistical analysis. Finally, we validate these new techniques within a full pipeline for cardiac segmentation applicable to large-scale cardiac MRI databases. The results obtained based on over 1200 cases from the Cardiac Atlas Project show the promise of fully automatic initialization and quality control for population studies

Automated deep phenotyping of the cardiovascular system using magnetic resonance imaging

Across a lifetime, the cardiovascular system must adapt to a great range of demands from the body. The individual changes in the cardiovascular system that occur in response to loading conditions are influenced by genetic susceptibility, and the pattern and extent of these changes have prognostic value. Brachial blood pressure (BP) and left ventricular ejection fraction (LVEF) are important biomarkers that capture this response, and their measurements are made at high resolution. Relatively, clinical analysis is crude, and may result in lost information and the introduction of noise. Digital information storage enables efficient extraction of information from a dataset, and this strategy may provide more precise and deeper measures to breakdown current phenotypes into their component parts. The aim of this thesis was to develop automated analysis of cardiovascular magnetic resonance (CMR) imaging for more detailed phenotyping, and apply these techniques for new biological insights into the cardiovascular response to different loading conditions. I therefore tested the feasibility and clinical utility of computational approaches for image and waveform analysis, recruiting and acquiring additional patient cohorts where necessary, and then applied these approaches prospectively to participants before and after six-months of exercise training for a first-time marathon. First, a multi-centre, multi-vendor, multi-field strength, multi-disease CMR resource of 110 patients undergoing repeat imaging in a short time-frame was assembled. The resource was used to assess whether automated analysis of LV structure and function is feasible on real-world data, and if it can improve upon human precision. This showed that clinicians can be confident in detecting a 9% change in EF or a 20g change in LV mass. This will be difficult to improve by clinicians because the greatest source of human error was attributable to the observer rather than modifiable factors. Having understood these errors, a convolutional neural network was trained on separate multi-centre data for automated analysis and was successfully generalizable to the real-world CMR data. Precision was similar to human analysis, and performance was 186 times faster. This real-world benchmarking resource has been made freely available (thevolumesresource.com). Precise automated segmentations were then used as a platform to delve further into the LV phenotype. Global LVEFs measured from CMR imaging in 116 patients with severe aortic stenosis were broken down into ~10 million regional measurements of structure and function, represented by computational three-dimensional LV models for each individual. A cardiac atlas approach was used to compile, label, segment and represent these data. Models were compared with healthy matched controls, and co-registered with follow-up one year after aortic valve replacement (AVR). This showed that there is a tendency to asymmetric septal hypertrophy in all patients with severe aortic stenosis (AS), rather than a characteristic specific to predisposed patients. This response to AS was more unfavourable in males than females (associated with higher NT-proBNP, and lower blood pressure), but was more modifiable with AVR. This was not detected using conventional analysis. Because cardiac function is coupled with the vasculature, a novel integrated assessment of the cardiovascular system was developed. Wave intensity theory was used to combine central blood pressure and CMR aortic blood flow-velocity waveforms to represent the interaction of the heart with the vessels in terms of traveling energy waves. This was performed and then validated in 206 individuals (the largest cohort to date), demonstrating inefficient ventriculo-arterial coupling in female sex and healthy ageing. CMR imaging was performed in 236 individuals before training for a first-time marathon and 138 individuals were followed-up after marathon completion. After training, systolic/diastolic blood pressure reduced by 4/3mmHg, descending aortic stiffness decreased by 16%, and ventriculo-arterial coupling improved by 14%. LV mass increased slightly, with a tendency to more symmetrical hypertrophy. The reduction in aortic stiffness was equivalent to a 4-year reduction in estimated biological aortic age, and the benefit was greater in older, male, and slower individuals. In conclusion, this thesis demonstrates that automating analysis of clinical cardiovascular phenotypes is precise with significant time-saving. Complex data that is usually discarded can be used efficiently to identify new biology. Deeper phenotypes developed in this work inform risk reduction behaviour in healthy individuals, and demonstrably deliver a more sensitive marker of LV remodelling, potentially enhancing risk prediction in severe aortic stenosis