Ocean Eddy Identification and Tracking using Neural Networks

Global climate change plays an essential role in our daily life. Mesoscale ocean eddies have a significant impact on global warming, since they affect the ocean dynamics, the energy as well as the mass transports of ocean circulation. From satellite altimetry we can derive high-resolution, global maps containing ocean signals with dominating coherent eddy structures. The aim of this study is the development and evaluation of a deep-learning based approach for the analysis of eddies. In detail, we develop an eddy identification and tracking framework with two different approaches that are mainly based on feature learning with convolutional neural networks. Furthermore, state-of-the-art image processing tools and object tracking methods are used to support the eddy tracking. In contrast to previous methods, our framework is able to learn a representation of the data in which eddies can be detected and tracked in more objective and robust way. We show the detection and tracking results on sea level anomalies (SLA) data from the area of Australia and the East Australia current, and compare our two eddy detection and tracking approaches to identify the most robust and objective method.Comment: accepted for International Geoscience and Remote Sensing Symposium 201

Flow pattern analysis for magnetic resonance velocity imaging

Blood flow in the heart is highly complex. Although blood flow patterns have been investigated by both computational modelling and invasive/non-invasive imaging techniques, their evolution and intrinsic connection with cardiovascular disease has yet to be explored. Magnetic resonance (MR) velocity imaging provides a comprehensive distribution of multi-directional in vivo flow distribution so that detailed quantitative analysis of flow patterns is now possible. However, direct visualisation or quantification of vector fields is of little clinical use, especially for inter-subject or serial comparison of changes in flow patterns due to the progression of the disease or in response to therapeutic measures. In order to achieve a comprehensive and integrated description of flow in health and disease, it is necessary to characterise and model both normal and abnormal flows and their effects. To accommodate the diversity of flow patterns in relation to morphological and functional changes, we have described in this thesis an approach of detecting salient topological features prior to analytical assessment of dynamical indices of the flow patterns. To improve the accuracy of quantitative analysis of the evolution of topological flow features, it is essential to restore the original flow fields so that critical points associated with salient flow features can be more reliably detected. We propose a novel framework for the restoration, abstraction, extraction and tracking of flow features such that their dynamic indices can be accurately tracked and quantified. The restoration method is formulated as a constrained optimisation problem to remove the effects of noise and to improve the consistency of the MR velocity data. A computational scheme is derived from the First Order Lagrangian Method for solving the optimisation problem. After restoration, flow abstraction is applied to partition the entire flow field into clusters, each of which is represented by a local linear expansion of its velocity components. This process not only greatly reduces the amount of data required to encode the velocity distribution but also permits an analytical representation of the flow field from which critical points associated with salient flow features can be accurately extracted. After the critical points are extracted, phase portrait theory can be applied to separate them into attracting/repelling focuses, attracting/repelling nodes, planar vortex, or saddle. In this thesis, we have focused on vortical flow features formed in diastole. To track the movement of the vortices within a cardiac cycle, a tracking algorithm based on relaxation labelling is employed. The constraints and parameters used in the tracking algorithm are designed using the characteristics of the vortices. The proposed framework is validated with both simulated and in vivo data acquired from patients with sequential MR examination following myocardial infarction. The main contribution of the thesis is in the new vector field restoration and flow feature abstraction method proposed. They allow the accurate tracking and quantification of dynamic indices associated with salient features so that inter- and intra-subject comparisons can be more easily made. This provides further insight into the evolution of blood flow patterns and permits the establishment of links between blood flow patterns and localised genesis and progression of cardiovascular disease.Open acces

Topology Based Flow Analysis and Superposition Effects

Using topology for feature analysis in flow fields faces several problems. First of all, not all features can be detected using topology based methods. Second, while in flow feature analysis the user is interested in a quantification of feature parameters like position, size, shape, radial velocity and other parameters of feature models, many of these parameters can not be determined using topology based methods alone. Additionally, in some applications it is advantageous to regard the vector field as a superposition of several, possibly simple, features. As topology based methods are quite sensitive to superposition effects, their precision and usability is limited in these cases. In this paper, topology based analysis and visualization of flow fields is estimated and compared to other feature based approaches demonstrating these problems

Extract and Characterize Hairpin Vortices in Turbulent Flows

Hairpin vortices are one of the most important vortical structures in turbulent flows. Extracting and characterizing hairpin vortices provides useful insight into many behaviors in turbulent flows. However, hairpin vortices have complex configurations and might be entangled with other vortices, making their extraction difficult. In this work, we introduce a framework to extract and separate hairpin vortices in shear driven turbulent flows for their study. Our method first extracts general vortical regions with a region-growing strategy based on certain vortex criteria (e.g., $\lambda_2$) and then separates those vortices with the help of progressive extraction of ($\lambda_2$) iso-surfaces in a top-down fashion. This leads to a hierarchical tree representing the spatial proximity and merging relation of vortices. After separating individual vortices, their shape and orientation information is extracted. Candidate hairpin vortices are identified based on their shape and orientation information as well as their physical characteristics. An interactive visualization system is developed to aid the exploration, classification, and analysis of hairpin vortices based on their geometric and physical attributes. We also present additional use cases of the proposed system for the analysis and study of general vortices in other types of flows.Comment: Accepted for presentation at IEEE VIS 2023. The paper will appear in IEEE Transactions on Visualization and Computer Graphic

Feature-based Visualization of Dense Integral Line Data

Feature-based visualization of flow fields has proven as an effective tool for flow analysis. While most flow visualization techniques operate on vector field data, our visualization techniques make use of a different simulation output: Particle Tracers. Our approach solely relies on integral lines that can be easily obtained from most simulation software. The task is the visualization of dense integral line data. We combine existing methods for streamline visualization, i.e. illumination, transparency, and halos, and add ambient occlusion for lines. But, this only solves one part of the problem: because of the high density of lines, visualization has to fight with occlusion, high frequency noise, and overlaps. As a solution we propose non-automated choices of transfer functions on curve properties that help highlighting important flow features like vortices or turbulent areas. These curve properties resemble some of the original flow properties. With the new combination of existing line drawing methods and the addition of ambient occlusion we improve the visualization of lines by adding better shape and depth cues. The intelligent use of transfer functions on curve properties reduces visual clutter and helps focusing on important features while still retaining context, as demonstrated in the examples given in this work

Comprehensive 4D velocity mapping of the heart and great vessels by cardiovascular magnetic resonance

<p>Abstract</p> <p>Background</p> <p>Phase contrast cardiovascular magnetic resonance (CMR) is able to measure all three directional components of the velocities of blood flow relative to the three spatial dimensions and the time course of the heart cycle. In this article, methods used for the acquisition, visualization, and quantification of such datasets are reviewed and illustrated.</p> <p>Methods</p> <p>Currently, the acquisition of 3D cine (4D) phase contrast velocity data, synchronized relative to both cardiac and respiratory movements takes about ten minutes or more, even when using parallel imaging and optimized pulse sequence design. The large resulting datasets need appropriate post processing for the visualization of multidirectional flow, for example as vector fields, pathlines or streamlines, or for retrospective volumetric quantification.</p> <p>Applications</p> <p>Multidirectional velocity acquisitions have provided 3D visualization of large scale flow features of the healthy heart and great vessels, and have shown altered patterns of flow in abnormal chambers and vessels. Clinically relevant examples include retrograde streams in atheromatous descending aortas as potential thrombo-embolic pathways in patients with cryptogenic stroke and marked variations of flow visualized in common aortic pathologies. Compared to standard clinical tools, 4D velocity mapping offers the potential for retrospective quantification of flow and other hemodynamic parameters.</p> <p>Conclusions</p> <p>Multidirectional, 3D cine velocity acquisitions are contributing to the understanding of normal and pathologically altered blood flow features. Although more rapid and user-friendly strategies for acquisition and analysis may be needed before 4D velocity acquisitions come to be adopted in routine clinical CMR, their capacity to measure multidirectional flows throughout a study volume has contributed novel insights into cardiovascular fluid dynamics in health and disease.</p

Assessing the performance of ultrafast vector flow imaging in the neonatal heart via multiphysics modeling and In vitro experiments

Ultrafast vector flow imaging would benefit newborn patients with congenital heart disorders, but still requires thorough validation before translation to clinical practice. This paper investigates 2-D speckle tracking (ST) of intraventricular blood flow in neonates when transmitting diverging waves at ultrafast frame rate. Computational and in vitro studies enabled us to quantify the performance and identify artifacts related to the flow and the imaging sequence. First, synthetic ultrasound images of a neonate's left ventricular flow pattern were obtained with the ultrasound simulator Field II by propagating point scatterers according to 3-D intraventricular flow fields obtained with computational fluid dynamics (CFD). Noncompounded diverging waves (opening angle of 60 degrees) were transmitted at a pulse repetition frequency of 9 kHz. ST of the B-mode data provided 2-D flow estimates at 180 Hz, which were compared with the CFD flow field. We demonstrated that the diastolic inflow jet showed a strong bias in the lateral velocity estimates at the edges of the jet, as confirmed by additional in vitro tests on a jet flow phantom. Furthermore, ST performance was highly dependent on the cardiac phase with low flows (< 5 cm/s), high spatial flow gradients, and out-of-plane flow as deteriorating factors. Despite the observed artifacts, a good overall performance of 2-D ST was obtained with a median magnitude underestimation and angular deviation of, respectively, 28% and 13.5 degrees during systole and 16% and 10.5 degrees during diastole