MR imaging of left-ventricular function : novel image acquisition and analysis techniques.

Many cardiac diseases, such as myocardial ischemia, secondary to coronary artery disease, may be identified and localized through the analysis of cardiac deformations. Early efforts for quantifying ventricular wall motion used surgical implantation and tracking of radiopaque markers with X-ray imaging in canine hearts [1]. Such techniques are invasive and affect the regional motion pattern of the ventricular wall during the marker tracking process and, clearly are not feasible clinically. Noninvasive imaging techniques are vital and have been widely applied to the clinic. MRI is a noninvasive imaging technique with the capability to monitor and assess the progression of cardiovascular diseases (CVD) so that effective procedures for the care and treatment of patients can be developed by physicians and researchers. It is capable of providing 3D analysis of global and regional cardiac function with great accuracy and reproducibility. In the past few years, numerous efforts have been devoted to cardiac motion recovery and deformation analysis from MR imaging sequences. In order to assess cardiac function, there are two categories of indices that are used: global and regional indices. Global indices include ejection fraction, cavity volume, and myocardial mass [2]. They are important indices for cardiac disease diagnosis. However, these global indices are not specific for regional analysis. A quantitative assessment of regional parameters may prove beneficial for the diagnosis of disease and evaluation of severity and the quantification of treatment [3]. Local measures, such as wall deformation and strain in all regions of the heart, can provide objective regional quantification of ventricular wall function and relate to the location and extent of ischemic injury. This dissertation is concerned with the development of novel MR imaging techniques and image postprocessing algorithms to analyze left ventricular deformations. A novel pulse sequence, termed Orthogonal CSPAMM (OCSPAMM), has been proposed which results in the same acquisition time as SPAMM for 2D deformation estimation while keeping the main advantages of CSPAMM [4,5]: i.e., maintaining tag contrast through-out the ECG cycle. Different from CSPAMM, in OCSPAMM the second tagging pulse orientation is rotated 90 degrees relative to the first one so that motion information can be obtained simultaneously in two directions. This reduces the acquisition time by a factor of two as compared to the traditional CSPAMM, in which two separate imaging sequences are applied per acquisition. With the application of OCSPAMM, the effect of tag fading encountered in SPAMM tagging due to Tl relaxation is mitigated and tag deformations can be visualized for the entire cardiac cycle, including diastolic phases. A multilevel B-spline fitting method (MBS) has been proposed which incorporates phase-based displacement information for accurate calculation of 2D motion and strain from tagged MRI [6, 7]. The proposed method combines the advantages of continuity and smoothness of MBS, and makes use of phase information derived from tagged MR images. Compared to previous 2D B-spline-based deformation analysis methods, MBS has the following advantages: 1) It can simultaneously achieve a smooth deformation while accurately approximating the given data set; 2) Computationally, it is very fast; and 3) It can produce more accurate deformation results. Since the tag intersections (intersections between two tag lines) can be extracted accurately and are more or less distributed evenly over the myocardium, MBS has proven effective for 2D cardiac motion tracking. To derive phase-based displacements, 2D HARP and SinMod analysis techniques [8,9] were employed. By producing virtual tags from HARP /SinMod and calculating intersections of virtual tag lines, more data points are obtained. In the reference frame, virtual tag lines are the isoparametric curves of an undeformed 2D B-spline model. In subsequent frames, the locations of intersections of virtual tag lines over the myocardium are updated with phase-based displacement. The advantage of the technique is that in acquiring denser myocardial displacements, it uses both real and virtual tag line intersections. It is fast and more accurate than 2D HARP and SinMod tracking. A novel 3D sine wave modeling (3D SinMod) approach for automatic analysis of 3D cardiac deformations has been proposed [10]. An accelerated 3D complementary spatial modulation of magnetization (CSPAMM) tagging technique [11] was used to acquire complete 3D+t tagged MR data sets of the whole heart (3 dynamic CSPAMM tagged MRI volume with tags in different orientations), in-vivo, in 54 heart beats and within 3 breath-holds. In 3D SinMod, the intensity distribution around each pixel is modeled as a cosine wave front. The principle behind 3D SinMod tracking is that both phase and frequency for each voxel are determined directly from the frequency analysis and the displacement is calculated from the quotient of phase difference and local frequency. The deformation fields clearly demonstrate longitudinal shortening during systole. The contraction of the LV base towards the apex as well as the torsional motion between basal and apical slices is clearly observable from the displacements. 3D SinMod can automatically process the image data to derive measures of motion, deformations, and strains between consecutive pair of tagged volumes in 17 seconds. Therefore, comprehensive 4D imaging and postprocessing for determination of ventricular function is now possible in under 10 minutes. For validation of 3D SinMod, 7 3D+t CSPAMM data sets of healthy subjects have been processed. Comparison of mid-wall contour deformations and circumferential shortening results by 3D SinMod showed good agreement with those by 3D HARP. Tag lines tracked by the proposed technique were also compared with manually delineated ones. The average errors calculated for the systolic phase of the cardiac cycles were in the sub-pixel range

On motion in dynamic magnetic resonance imaging: Applications in cardiac function and abdominal diffusion

La imagen por resonancia magnética (MRI), hoy en día, representa una potente herramienta para el diagnóstico clínico debido a su flexibilidad y sensibilidad a un amplio rango de propiedades del tejido. Sus principales ventajas son su sobresaliente versatilidad y su capacidad para proporcionar alto contraste entre tejidos blandos. Gracias a esa versatilidad, la MRI se puede emplear para observar diferentes fenómenos físicos dentro del cuerpo humano combinando distintos tipos de pulsos dentro de la secuencia. Esto ha permitido crear distintas modalidades con múltiples aplicaciones tanto biológicas como clínicas. La adquisición de MR es, sin embargo, un proceso lento, lo que conlleva una solución de compromiso entre resolución y tiempo de adquisición (Lima da Cruz, 2016; Royuela-del Val, 2017). Debido a esto, la presencia de movimiento fisiológico durante la adquisición puede conllevar una grave degradación de la calidad de imagen, así como un incremento del tiempo de adquisición, aumentando así tambien la incomodidad del paciente. Esta limitación práctica representa un gran obstáculo para la viabilidad clínica de la MRI. En esta Tesis Doctoral se abordan dos problemas de interés en el campo de la MRI en los que el movimiento fisiológico tiene un papel protagonista. Éstos son, por un lado, la estimación robusta de parámetros de rotación y esfuerzo miocárdico a partir de imágenes de MR-Tagging dinámica para el diagnóstico y clasificación de cardiomiopatías y, por otro, la reconstrucción de mapas del coeficiente de difusión aparente (ADC) a alta resolución y con alta relación señal a ruido (SNR) a partir de adquisiciones de imagen ponderada en difusión (DWI) multiparamétrica en el hígado.Departamento de Teoría de la Señal y Comunicaciones e Ingeniería TelemáticaDoctorado en Tecnologías de la Información y las Telecomunicacione

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

Computer-Assisted Electroanatomical Guidance for Cardiac Electrophysiology Procedures

Cardiac arrhythmias are serious life-threatening episodes affecting both the aging population and younger patients with pre-existing heart conditions. One of the most effective therapeutic procedures is the minimally-invasive catheter-driven endovascular electrophysiology study, whereby electrical potentials and activation patterns in the affected cardiac chambers are measured and subsequent ablation of arrhythmogenic tissue is performed. Despite emerging technologies such as electroanatomical mapping and remote intraoperative navigation systems for improved catheter manipulation and stability, successful ablation of arrhythmias is still highly-dependent on the operator’s skills and experience. This thesis proposes a framework towards standardisation in the electroanatomical mapping and ablation planning by merging knowledge transfer from previous cases and patient-specific data. In particular, contributions towards four different procedural aspects were made: optimal electroanatomical mapping, arrhythmia path computation, catheter tip stability analysis, and ablation simulation and optimisation. In order to improve the intraoperative electroanatomical map, anatomical areas of high mapping interest were proposed, as learned from previous electrophysiology studies. Subsequently, the arrhythmic wave propagation on the endocardial surface and potential ablation points were computed. The ablation planning is further enhanced, firstly by the analysis of the catheter tip stability and the probability of slippage at sparse locations on the endocardium and, secondly, by the simulation of the ablation result from the computation of convolutional matrices which model mathematically the ablation process. The methods proposed by this thesis were validated on data from patients with complex congenital heart disease, who present unusual cardiac anatomy and consequently atypical arrhythmias. The proposed methods also build a generic framework for computer guidance of electrophysiology, with results showing complementary information that can be easily integrated into the clinical workflow.Open Acces

Electromagnetic Wave Propagation for Industry and Biomedical Applications

This book highlights original research and high-quality technical briefs on electromagnetic wave propagation, radiation, and scattering, and their applications in industry and biomedical engineering. It also presents recent research achievements in the theoretical, computational, and experimental aspects of electromagnetic wave propagation, radiation, and scattering. The book is divided into three sections. Section 1 consists of chapters with general mathematical methods and approaches to the forward and inverse problems of wave propagation. Section 2 presents the problems of wave propagation in superconducting materials and porous media. Finally, Section 3 discusses various industry and biomedical applications of electromagnetic wave propagation, radiation, and scattering

Analysis of first pass myocardial perfusion imaging with magnetic resonance

Early diagnosis and localisation of myocardial perfusion defects is an important step in the treatment of coronary artery disease. Thus far, coronary angiography is the conventional standard investigation for patients with known or suspected coronary artery disease and it provides information about the presence and location of coronary stenoses. In recent years, the development of myocardial perfusion CMR has extended the role of MR in the evaluation of ischaemic heart disease beyond the situations where there have already been gross myocardial changes such as acute infarction or scarring. The ability to non-invasively evaluate cardiac perfusion abnormalities before pathologic effects occur, or as follow-up to therapy, is important to the management of patients with coronary artery disease. Whilst limited multi-slice 2D CMR perfusion studies are gaining increased clinical usage for quantifying gross ischaemic burden, research is now directed towards complete 3D coverage of the myocardium for accurate localisation of the extent of possible defects. In 3D myocardial perfusion imaging, a complete volumetric data set has to be acquired for each cardiac cycle in order to study the first pass of the contrast bolus. This normally requires a relatively large acquisition window within each cardiac cycle to ensure a comprehensive coverage of the myocardium and reasonably high resolution of the images. With multi-slice imaging, long axis cardiac motion during this large acquisition window can cause the myocardium imaged in different cross- sections to be mis-registered, i.e., some part of the myocardium may be imaged more than twice whereas other parts may be missed out completely. This type of mis-registration is difficult to correct for by using post-processing techniques. The purpose of this thesis is to investigate techniques for tracking through plane motion during 3D myocardial perfusion imaging, and a novel technique for extracting intrinsic relationships between 3D cardiac deformation due to respiration and multiple ID real-time measurable surface intensity traces is developed. Despite the fact that these surface intensity traces can be strongly coupled with each other but poorly correlated with respiratory induced cardiac deformation, we demonstrate how they can be used to accurately predict cardiac motion through the extraction of latent variables of both the input and output of the model. The proposed method allows cross-modality reconstruction of patient specific models for dense motion field prediction, which after initial modelling can be use in real-time prospective motion tracking or correction. In CMR, new imaging sequences have significantly reduced the acquisition window whilst maintaining the desired spatial resolution. Further improvements in perfusion imaging will require the application of parallel imaging techniques or making full use of the information content of the ¿-space data. With this thesis, we have proposed RR-UNFOLD and RR-RIGR for significantly reducing the amount of data that is required to reconstruct the perfusion image series. The methods use prospective diaphragmatic navigator echoes to ensure UNFOLD and RIGR are carried out on a series of images that are spatially registered. An adaptive real-time re-binning algorithm is developed for the creation of static image sub-series related to different levels of respiratory motion. Issues concerning temporal smoothing of tracer kinetic signals and residual motion artefact are discussed, and we have provided a critical comparison of the relative merit and potential pitfalls of the two techniques. In addition to the technical and theoretical descriptions of the new methods developed, we have also provided in this thesis a detailed literature review of the current state-of-the-art in myocardial perfusion imaging and some of the key technical challenges involved. Issues concerning the basic background of myocardial ischaemia and its functional significance are discussed. Practical solutions to motion tracking during imaging, predictive motion modelling, tracer kinetic modelling, RR-UNFOLD and RR-RIGR are discussed, all with validation using patient and normal subject data to demonstrate both the strength and potential clinical value of the proposed techniques.Open acces

Estimation of tissue contractility from cardiac cine-MRI using a biomechanical heart model

International audienceThe objective of this paper is to propose and assess an estimation procedure - based on data assimilation principles - well-suited to obtain some regional values of key biophysical parameters in a beating heart model, using actual Cine-MR images. The motivation is twofold: (1) to provide an automatic tool for personalizing the characteristics of a cardiac model in order to achieve predictivity in patient-specific modeling, and (2) to obtain some useful information for diagnosis purposes in the estimated quantities themselves. In order to assess the global methodology we specifically devised an animal experiment in which a controlled infarct was produced and data acquired before and after infarction, with an estimation of regional tissue contractility - a key parameter directly affected by the pathology - performed for every measured stage. After performing a preliminary assessment of our proposed methodology using synthetic data, we then demonstrate a full-scale application by first estimating contractility values associated with 6 regions based on the AHA subdivision, before running a more detailed estimation using the actual AHA segments. The estimation results are assessed by comparison with the medical knowledge of the specific infarct, and with late enhancement MR images. We discuss their accuracy at the various subdivision levels, in the light of the inherent modeling limitations and of the intrinsic information contents featured in the data

Advancements and Breakthroughs in Ultrasound Imaging

Ultrasonic imaging is a powerful diagnostic tool available to medical practitioners, engineers and researchers today. Due to the relative safety, and the non-invasive nature, ultrasonic imaging has become one of the most rapidly advancing technologies. These rapid advances are directly related to the parallel advancements in electronics, computing, and transducer technology together with sophisticated signal processing techniques. This book focuses on state of the art developments in ultrasonic imaging applications and underlying technologies presented by leading practitioners and researchers from many parts of the world

Modified mass-spring system for physically based deformation modeling

Mass-spring systems are considered the simplest and most intuitive of all deformable models. They are computationally efficient, and can handle large deformations with ease. But they suffer several intrinsic limitations. In this book a modified mass-spring system for physically based deformation modeling that addresses the limitations and solves them elegantly is presented. Several implementations in modeling breast mechanics, heart mechanics and for elastic images registration are presented

Modified mass-spring system for physically based deformation modeling

Mass-spring systems are considered the simplest and most intuitive of all deformable models. They are computationally efficient, and can handle large deformations with ease. But they suffer several intrinsic limitations. In this book a modified mass-spring system for physically based deformation modeling that addresses the limitations and solves them elegantly is presented. Several implementations in modeling breast mechanics, heart mechanics and for elastic images registration are presented