Higher order dynamic mode decomposition: From fluid dynamics to heart disease analysis

Deriving phenotype-representative left ventricular flow patterns by reduced-order modeling and classification

BACKGROUND

Extracting phenotype-representative flow patterns and their associated numerical metrics is a bottleneck in the clinical translation of advanced cardiac flow imaging modalities. We hypothesized that reduced-order models (ROMs) are a suitable strategy for deriving simple and interpretable clinical metrics of intraventricular flow suitable for further assessments. Combined with machine learning (ML) flow-based ROMs could provide new insight to help diagnose and risk-stratify patients.

METHODS

We analyzed 2D color-Doppler echocardiograms of 81 non-ischemic dilated cardiomyopathy (DCM) patients, 51 hypertrophic cardiomyopathy (HCM) patients, and 77 normal volunteers (Control). We applied proper orthogonal decomposition (POD) to build patient-specific and cohort-specific ROMs of LV flow. Each ROM aggregates a low number of components representing a spatially dependent velocity map modulated along the cardiac cycle by a time-dependent coefficient. We tested three classifiers using deliberately simple ML analyses of these ROMs with varying supervision levels. In supervised models, hyperparameter grid search was used to derive the ROMs that maximize classification power. The classifiers were blinded to LV chamber geometry and function. We ran vector flow mapping on the color-Doppler sequences to help visualize flow patterns and interpret the ML results.

RESULTS

POD-based ROMs stably represented each cohort through 10-fold cross-validation. The principal POD mode captured >80 % of the flow kinetic energy (KE) in all cohorts and represented the LV filling/emptying jets. Mode 2 represented the diastolic vortex and its KE contribution ranged from <1 % (HCM) to 13 % (DCM). Semi-unsupervised classification using patient-specific ROMs revealed that the KE ratio of these two principal modes, the vortex-to-jet (V2J) energy ratio, is a simple, interpretable metric that discriminates DCM, HCM, and Control patients. Receiver operating characteristic curves using V2J as classifier had areas under the curve of 0.81, 0.91, and 0.95 for distinguishing HCM vs. Control, DCM vs. Control, and DCM vs. HCM, respectively.

CONCLUSIONS

Modal decomposition of cardiac flow can be used to create ROMs of normal and pathological flow patterns, uncovering simple interpretable flow metrics with power to discriminate disease states, and particularly suitable for further processing using ML.

Camera-Based Dynamic Vibration Analysis Using Transformer-Based Model CoTracker and Dynamic Mode Decomposition

Accelerometers are commonly used to measure vibrations for condition monitoring in mechanical and civil structures; however, their high cost and point-based measurement approach present practical limitations. With rapid advancements in computer vision and deep learning, research into tracking the motion of individual pixels with vision cameras has increased. The recently developed CoTracker, a transformer-based model, has demonstrated excellence in motion tracking, yet its performance in measuring structural vibrations has not been fully explored. This paper investigates the efficacy of the CoTracker model in extracting full-field structural vibrations using cameras. It is initially applied to capture the dense point movements in video sequences of a cantilever beam recorded using a high-speed camera. Subsequently, modal analysis using delay-embedding dynamic mode decomposition (DMD) is conducted to extract modal parameters including natural frequencies, damping ratios, and mode shapes. The results, benchmarked against those from a reference accelerometer and the Finite Element Method (FEM) result, demonstrate CoTracker's high potential for general applicability in structural vibration measurements.

Deriving Explainable Metrics of Left Ventricular Flow by Reduced-Order Modeling and Classification

Background

Extracting explainable flow metrics is a bottleneck to the clinical translation of advanced cardiac flow imaging modalities. We hypothesized that reduced-order models (ROMs) of intraventricular flow are a suitable strategy for deriving simple and interpretable clinical metrics suitable for further assessments. Combined with machine learning (ML) flow-based ROMs could provide new insight to help diagnose and risk-stratify patients.

Methods

We analyzed 2D color-Doppler echocardiograms of 81 non-ischemic dilated cardiomyopathy (DCM) patients, 51 hypertrophic cardiomyopathy (HCM) patients, and 77 normal volunteers (Control). We applied proper orthogonal decomposition (POD) to build patient-specific and cohort-specific ROMs of LV flow. Each ROM aggregates a low number of components representing a spatially dependent velocity map modulated along the cardiac cycle by a time-dependent coefficient. We tested three classifiers using deliberately simple ML analyses of these ROMs with varying supervision levels. In supervised models, hyperparameter gridsearch was used to derive the ROMs that maximize classification power. The classifiers were blinded to LV chamber geometry and function. We ran vector flow mapping on the color-Doppler sequences to help visualize flow patterns and interpret the ML results.

Results

POD-based ROMs stably represented each cohort through 10-fold cross-validation. The principal POD mode captured >80% of the flow kinetic energy (KE) in all cohorts and represented the LV filling/emptying jets. Mode 2 represented the diastolic vortex and its KE contribution ranged from <1% (HCM) to 13% (DCM). Semi-unsupervised classification using patient-specific ROMs revealed that the KE ratio of these two principal modes, the vortex-to-jet (V2J) energy ratio, is a simple, interpretable metric that discriminates DCM, HCM, and Control patients. Receiver operating characteristic curves using V2J as classifier had areas under the curve of 0.81, 0.91, and 0.95 for distinguishing HCM vs. Control, DCM vs. Control, and DCM vs. HCM, respectively.

Conclusions

Modal decomposition of cardiac flow can be used to create ROMs of normal and pathological flow patterns, uncovering simple interpretable flow metrics with power to discriminate disease states, and particularly suitable for further processing using ML.

Dynamic mode decomposition of dynamic MRI for assessment of pulmonary ventilation and perfusion

PURPOSE

To introduce dynamic mode decomposition (DMD) as a robust alternative for the assessment of pulmonary functional information from dynamic non-contrast-enhanced acquisitions.

METHODS

Pulmonary fractional ventilation and normalized perfusion maps were obtained using DMD from simulated phantoms as well as in vivo dynamic acquisitions of healthy volunteers at 1.5T. The performance of DMD was compared with conventional Fourier decomposition (FD) and matrix pencil (MP) methods in estimating functional map values. The proposed method was evaluated based on estimated signal amplitude in functional maps across varying number of measurements.

RESULTS

Quantitative assessments performed on phantoms and in vivo measurements indicate that DMD is capable of successfully obtaining pulmonary functional maps. Specifically, compared to FD and MP methods, DMD is able to reduce variations in estimated amplitudes across different number of measurements. This improvement is evident in the fractional ventilation and normalized perfusion maps obtain from phantom simulations with frequency variations and noise, as well as in the maps obtained from in vivo measurements.

CONCLUSIONS

A robust method for accurately estimating pulmonary ventilation and perfusion related signal changes in dynamic acquisitions is presented. The proposed method uses DMD to obtain functional maps reliably, while reducing amplitude variations caused by differences in number of measurements.

A novel data-driven method for the analysis and reconstruction of cardiac cine MRI

Cardiac cine magnetic resonance imaging (MRI) can be considered the optimal criterion for measuring cardiac function. This imaging technique can provide us with detailed information about cardiac structure, tissue composition and even blood flow, which makes it highly used in medical science. But due to the image time acquisition and several other factors the MRI sequences can easily get corrupted, causing radiologists to misdiagnose 40 million people worldwide each and every single year. Hence, the urge to decrease these numbers, researchers from different fields have been introducing novel tools and methods in the medical field. Aiming to the same target, we consider in this work the application of the higher order dynamic mode decomposition (HODMD) technique. The HODMD algorithm is a linear method, which was originally introduced in the fluid dynamics domain, for the analysis of complex systems. Nevertheless, the proposed method has extended its applicability to numerous domains, including medicine. In this work, HODMD in used to analyze sets of MR images of a heart, with the ultimate goal of identifying the main patterns and frequencies driving the heart dynamics. Furthermore, a novel interpolation algorithm based on singular value decomposition combined with HODMD is introduced, providing a three-dimensional reconstruction of the heart. This algorithm is applied (i) to reconstruct corrupted or missing images, and (ii) to build a reduced order model of the heart dynamics.