Generación de modelos físicos personalizados y segmentación de la pared arterial y el trombo intraluminal en imágenes de aneurismas aórticos abdominales

Esta tesis doctoral está encuadrada dentro del marco general de la ingeniería biomédica aplicada al tratamiento de las enfermedades cardiovasculares, enfermedades que provocan alrededor de 1.9 millones (40%) de muertes al año en la Unión Europea. En este contexto surge el proyecto europeo SCATh-Smart Catheterization, cuyo objetivo principal es mejorar los procedimientos de cateterismo aórtico introduciendo nuevas tecnologías de planificación y navegación quirúrgica y minimizando el uso de fluoroscopía. En particular, esta tesis aborda el modelado y diagnóstico de aneurismas aórticos abdominales (AAA) y del trombo intraluminal (TIL), allí donde esté presente, así como la segmentación de estas estructuras en imágenes preoperatorias de RM. Los modelos físicos específicos del paciente, construidos a partir de imágenes médicas preoperatorias, tienen múltiples usos, que van desde la evaluación preoperatoria de estructuras anatómicas a la planificación quirúrgica para el guiado de catéteres. En el diagnóstico y tratamiento de AAA, los modelos físicos son útiles a la hora de evaluar diversas variables biomecánicas y fisiológicas de las estructuras vasculares. Existen múltiples técnicas que requieren de la generación de modelos físicos que representen la anatomía vascular. Una de las principales aplicaciones de los modelos físicos es el análisis de elementos finitos (FE). Las simulaciones de FE para AAA pueden ser específicas para el paciente y permiten modelar estados de estrés complejos, incluyendo los efectos provocados por el TIL. La aplicación de métodos numéricos de análisis tiene como requisito previo la generación de una malla computacional que representa la geometría de interés mediante un conjunto de elementos poliédricos, siendo los hexaédricos los que presentan mejores resultados. En las estructuras vasculares, generar mallas hexaédricas es un proceso especialmente exigente debido a la compleja anatomía 3D ramificada. La mayoría de los AAA se encuentran situados en la bifurcación de la arteria aorta en las arterias iliacas y es necesario modelar de manera fiel dicha bifurcación. En el caso de que la sangre se estanque en el aneurisma provocando un TIL, éste forma una estructura adyacente a la pared aórtica. De este modo, el contorno externo del TIL es el mismo que el contorno interno de la pared, por lo que las mallas resultantes deben reflejar esta particularidad, lo que se denomina como "mallas conformadas". El fin último de este trabajo es modelar las estructuras vasculares de modo que proporcionen nuevas herramientas para un mejor diagnóstico clínico, facilitando medidas de riesgo de rotura de la arteria, presión sistólica o diastólica, etc. Por tanto, el primer objetivo de esta tesis es diseñar un método novedoso y robusto para generar mallas hexaédricas tanto de la pared aórtica como del trombo. Para la identificación de estas estructuras se utilizan imágenes de resonancia magnética (RM). Deben mantenerse sus propiedades de adyacencia utilizando elementos de alta calidad, prestando especial atención al modelado de la bifurcación y a que sean adecuadas para el análisis de FE. El método tiene en cuenta la evolución de la línea central del vaso en el espacio tridimensional y genera la malla directamente a partir de las imágenes segmentadas, sin necesidad de reconstruir superficies triangulares. Con el fin de reducir la intervención del usuario en el proceso de generación de las mallas, es también objetivo de esta tesis desarrollar un método de segmentación semiautomática de las distintas estructuras de interés. Las principales contribuciones de esta tesis doctoral son: 1. El diseño, implementación y evaluación de un algoritmo de generación de mallas hexaédricas conformadas de la pared y el TIL a partir de los contornos segmentados en imágenes de RM. Se ha llevado a cabo una evaluación de calidad que determine su aplicabilidad a métodos de FE. Los resultados demuestran que el algoritmo desarrollado genera mallas conformadas de alta calidad incluso en la región de la bifurcación, que son adecuadas para su uso en métodos de análisis de FE. 2. El diseño, implementación y evaluación de un método de segmentación automático de las estructuras de interés. La luz arterial se segmenta de manera semiautomática utilizando un software disponible a partir de imágenes de RM con contraste. Los resultados de este proceso sirven de inicialización para la segmentación automática de las caras interna y externa de la pared aórtica utilizando métodos basado en modelos de textura y forma a partir de imágenes de RM sin contraste. Los resultados demuestran que el algoritmo desarrollado proporciona segmentaciones fieles de las distintas estructuras de interés. En conclusión, el trabajo realizado en esta tesis doctoral corrobora las hipótesis de investigación postuladas, y pretende servir como aportación para futuros avances en la generación de modelos físicos de geometrías biológicas. ABSTRACT The frame of this PhD Thesis is the biomedical engineering applied to the treatment of cardiovascular diseases, which cause around 1.9 million deaths per year in the European Union and suppose about 40% of deaths per year. In this context appears the European project SCATh-Smart Catheterization. The main objective of this project is creating a platform which improves the navigation of catheters in aortic catheterization minimizing the use of fluoroscopy. In the framework of this project, the specific field of this PhD Thesis is the diagnosis and modeling of abdominal aortic aneurysm (AAAs) and the intraluminal thrombus (ILT) whenever it is present. Patient-specific physical models built from preoperative imaging are becoming increasingly important in the area of minimally invasive surgery. These models can be employed for different purposes, such as the preoperatory evaluation of anatomic structures or the surgical planning for catheter guidance. In the specific case of AAA diagnosis and treatment, physical models are especially useful for evaluating pressures over vascular structures. There are multiple techniques that require the generation of physical models which represent the target anatomy. Finite element (FE) analysis is one the principal applications for physical models. FE simulations for AAA may be patient-specific and allow modeling biomechanical and physiological variables including those produced by ILT, and also the segmentation of those anatomical structures in preoperative MR images. Applying numeric methods requires the generation of a proper computational mesh. These meshes represent the patient anatomy using a set of polyhedral elements, with hexahedral elements providing better results. In the specific case of vascular structures, generating hexahedral meshes is a challenging task due to the complex 3D branching anatomy. Each patient’s aneurysm is unique, characterized by its location and shape, and must be accurately represented for subsequent analyses to be meaningful. Most AAAs are located in the region where the aorta bifurcates into the iliac arteries and it is necessary to model this bifurcation precisely and reliably. If blood stagnates in the aneurysm and forms an ILT, it exists as a conforming structure with the aortic wall, i.e. the ILT’s outer contour is the same as the wall’s inner contour. Therefore, resulting meshes must also be conforming. The main objective of this PhD Thesis is designing a novel and robust method for generating conforming hexahedral meshes for the aortic wall and the thrombus. These meshes are built using largely high-quality elements, especially at the bifurcation, that are suitable for FE analysis of tissue stresses. The method accounts for the evolution of the vessel’s centerline which may develop outside a single plane, and generates the mesh directly from segmented images without the requirement to reconstruct triangular surfaces. In order to reduce the user intervention in the mesh generation process is also a goal of this PhD. Thesis to develop a semiautomatic segmentation method for the structures of interest. The segmentation is performed from magnetic resonance image (MRI) sequences that have tuned to provide high contrast for the arterial tissue against the surrounding soft tissue, so that we determine the required information reliably. The main contributions of this PhD Thesis are: 1. The design, implementation and evaluation of an algorithm for generating hexahedral conforming meshes of the arterial wall and the ILT from the segmented contours. A quality inspection has been applied to the meshes in order to determine their suitability for FE methods. Results show that the developed algorithm generates high quality conforming hexahedral meshes even at the bifurcation region. Thus, these meshes are suitable for FE analysis. 2. The design, implementation and evaluation of a semiautomatic segmentation method for the structures of interest. The lumen is segmented in a semiautomatic way from contrast filled MRI using an available software. The results obtained from this process are used to initialize the automatic segmentation of the internal and external faces of the aortic wall. These segmentations are performed by methods based on texture and shape models from MRI with no contrast. The results show that the algorithm provides faithful segmentations of the structures of interest requiring minimal user intervention. In conclusion, the work undertaken in this PhD. Thesis verifies the investigation hypotheses. It intends to serve as basis for future physical model generation of proper biological anatomies used by numerical methods.

A new microfluidic methodology to assess the haemocompatibility of magnetic nanoparticles designed for theranostic applications

The remarkable physicochemical properties of magnetic nanoparticles (MNPs) at the nanoscale have boosted the development of new and promising strategies for the simultaneous diagnosis and treatment of diseases, particularly in cancer therapy Ð the so-called theranostic applications (1). In these strategies, the intrinsic superparamagnetic properties of MNPs have been exploited to gain access into multifunctional systems able to simultaneously perform as enhanced magnetic resonance imaging (MRI) contrast agents, efficient nanocarriers for drug delivery and nanoheaters in magnetic hyperthermia based therapy (2), among others.

Repercussões cardiovasculares no adulto : pós reparação de tetralogia de Fallot

Trabalho Final do Curso de Mestrado Integrado em Medicina, Faculdade de Medicina, Universidade de Lisboa, 2014

Applications of spin resonance techniques to the study of radiation effects,

Includes bibliography.

Effects of revascularisation and contractile reserve on left ventricular remodelling in patients with impaired left ventricular function

Objective: We sought to define the influence of revascularisation and contractile reserve on left ventricular (LV) remodelling in patients with LV dysfunction after myocardial infarction. Revascularisation of viable myocardium is associated with improved regional function, but the effect on remodelling is undefined. Methods: We studied 70 patients with coronary artery disease and LV dysfunction, 31 of whom underwent revascularisation. A standard dobutamine stress echocardiogram (DbE) was carried out. All patients underwent standard medical treatment; the decision to revascularise was made clinically, independent of this study. LV volumes and ejection fraction were measured by 3D echocardiography at baseline and after an average of 40 weeks. Results: There was no significant difference in baseline ejection fraction or volumes between patients who underwent revascularisation and the remainder. Compared to medically treated patients, revascularised patients had significant improvements in ejection fraction and end-systolic volume in follow-up. The impact of baseline variables on remodelling was assessed by dividing patients into tertiles of LV ejection fraction and volumes. Revascularised patients in the lowest tertile of ejection fraction at baseline (<38%) had a significant improvement in end-systolic volume and ejection fraction, larger than obtained in medically treated patients with low ejection fraction. Revascularised patients with an ejection fraction >38% did not show significant improvement in volumes compared to baseline. Revascularised patients in the largest tertiles of end-systolic (>88 ml) or end-diastolic volume (>149 ml) at baseline had a significant improvement in end-systolic volume. Conclusion: Remodeling appears to occur independent of the presence of regional contractile reserve but does correlate with the volume response to low-dose dobutamine. (C) 2003 Elsevier Ireland Ltd. All rights reserved.

Reproducibility and accuracy of echocardiographic measurements of left ventricular parameters using real-time three-dimensional echocardiography

OBJECTIVES We sought to determine whether assessment of left ventricular (LV) function with real-time (RT) three-dimensional echocardiography (3DE) could reduce the variation of sequential LV measurements and provide greater accuracy than two-dimensional echocardiography (2DE). BACKGROUND Real-time 3DE has become feasible as a standard clinical tool, but its accuracy for LV assessment has not been validated. METHODS Unselected patients (n = 50; 41 men; age, 64 +/- 8 years) presenting for evaluation of LV function were studied with 2DE and RT-3DE. Test-retest variation was performed by a complete restudy by a separate sonographer within 1 h without alteration of hemodynamics or therapy. Magnetic resonance imaging (MRI) images were obtained during a breath-hold, and measurements were made off-line. RESULTS The test-retest variation showed similar measurements for volumes but wider scatter of LV mass measurements with M-mode and 2DE than 3DE. The average MRI end-diastolic volume was 172 +/- 53 ml; LV volumes were underestimated by 2DE (mean difference, -54 +/- 33; p < 0.01) but only slightly by RT-3DE (-4 +/- 29; p = 0.31). Similarly, end-systolic volume by MRI (91 +/- 53 ml) was underestimated by 2DE (mean difference, -28 +/- 28; p < 0.01) and by RT-3DE (mean difference, -3 +/- 18; p = 0.23). Ejection fraction by MRI was similar by 2DE (p = 0.76) and RT-3DE (p = 0.74). Left ventricular mass (183 +/- 50 g) was overestimated by M-mode (mean difference, 68 +/- 86 g; p < 0.01) and 2DE (16 +/- 57; p = 0.04) but not RT-3DE (0 +/- 38 g; p = 0.94). There was good inter- and intra-observer correlation between RT-3DE by two sonographers for volumes, ejection fraction, and mass. CONCLUSIONS Real-time 3DE is a feasible approach to reduce test-retest variation of LV volume, ejection fraction, and mass measurements in follow-up LV assessment in daily practice. (C) 2004 by the American College of Cardiology Foundation.

Clinical decision-making and myocardial viability: current perspectives

Not all myocardium involved in a myocardial infarction is dead or irreversibly damaged. The balance between the amount of scar and live tissue, and the nature of the live tissue, determine the likelihood that contractile function will improve after revascularisation. This improvement (which defines viability) may be predicted with about 80% accuracy using several techniques. This review examines the determinants of functional recovery and how they may be integrated in making decisions regarding revascularisation. (Intern Med J 2005; 35: 118–125)

Relationship of extent and nature of dysfunctional myocardium to brain natriuretic peptide in patients with ischemic left ventricular dysfunction

We studied the relationship between brain natriuretic peptide (BNP) levels and viable myocardium and ischemic myocardium, regional scar and regional contractile function. Fifty-nine patients underwent dobutamine echocardiography and magnetic resonance imaging and resting BNP levels were determined. By magnetic resonance imaging, total extent of dysfunctional myocardium correlated strongest with BNP (r = 0.60, p < 0.0001). The extent of scar, viability and ischemia also correlated. At dobutamine echocardiography, a composite of dysfunctional and ischemic myocardium was the strongest correlate of BNP (r = 0.48, p < 0.0001), with less strong correlations by global parameters. The extent of dysfunctional myocardium, rather than its nature determines BNP levels.

Use of real-time three-dimensional echocardiography to measure left atrial volume: Comparison with other echocardiographic techniques

Objectives: Left atrial (LA) volume (LAV) is a prognostically important biomarker for diastolic dysfunction, but its reproducibility on repeated testing is not well defined. LA assessment with 3-dimensional. (3D) echocardiography (3DE) has been validated against magnetic resonance imaging, and we sought to assess whether this was superior to existing measurements for sequential echocardiographic follow-up. Methods: Patients (n = 100; 81 men; age 56 +/- 14 years) presenting for LA evaluation were studied with M-mode (MM) echocardiography, 2-dimensional (2D) echocardiography, and 3DE. Test-retest variation was performed by a complete restudy by a separate sonographer within 1 hour without alteration of hemodynamics or therapy. In all, 20 patients were studied for interobserver and intraobserver variation. LAVs were calculated by using M-mode diameter and planimetered atrial area in the apical. 4-chamber view to calculate an assumed sphere, as were prolate ellipsoid, Simpson's biplane, and biplane area-length methods. All were compared with 3DE. Results: The average LAV was 72 +/- 27 mL by 3DE. There was significant underestimation of LAV by M-mode (35 +/- 20 mL, r = 0.66, P < .01). The 3DE and various 2D echocardiographic techniques were well correlated: LA planimetry (85 +/- 38 mL, r = 0.77, P < .01), prolate ellipsoid (73 +/- 36 mL, r = 0.73, P = .04), area-length (64 +/- 30 mL, r = 0.74, P < .01), and Simpson's biplane (69 +/- 31 mL, r = 0.78, P = .06). Test-retest variation for 3DE was most favorable (r = 0.98, P < .01), with the prolate ellipsoid method showing most variation. Interobserver agreement between measurements was best for 3DE (r = 0.99, P < .01), with M-mode the worst (r = 0.89, P < .01). Intraobserver results were similar to interobserver, the best correlation for 3DE (r = 0.99, P < .01), with LA planimetry the worst (r = 0.91, P < .01). Conclusions. The 2D measurements correlate closely with 3DE. Follow-up assessment in daily practice appears feasible and reliable with both 2D and 3D approaches.

Analysis of lipoproteins using 2D diffusion-edited NMR spectroscopy and multi-way chemometrics

This study represents the first application of multi-way calibration by N-PLS and multi-way curve resolution by PARAFAC to 2D diffusion-edited H-1 NMR spectra. The aim of the analysis was to evaluate the potential for quantification of lipoprotein main- and subtractions in human plasma samples. Multi-way N-PLS calibrations relating the methyl and methylene peaks of lipoprotein lipids to concentrations of the four main lipoprotein fractions as well as 11 subfractions were developed with high correlations (R = 0.75-0.98). Furthermore, a PARAFAC model with four chemically meaningful components was calculated from the 2D diffusion-edited spectra of the methylene peak of lipids. Although the four extracted PARAFAC components represent molecules of sizes that correspond to the four main fractions of lipoproteins, the corresponding concentrations of the four PARAFAC components proved not to be correlated to the reference concentrations of these four fractions in the plasma samples as determined by ultracentrifugation. These results indicate that NMR provides complementary information on the classification of lipoprotein fractions compared to ultracentrifugation. (C) 2004 Elsevier B.V. All rights reserved.

Comparison between myocardial integrated backscatter and strain rate imaging in the quantitative assessment of subepicardial function post-myocardial infarction

Transmural extent of infarction (TME) may be an important determinant of functional recovery and remodeling. Recent animal data suggest that strain rate imaging (SRI) maybe able to identify subendocardial ischemia.We compared SRI and cyclic variation of integrated backscatter (CVIB) for predicting TME in the quantitative assessment of regional subepicardial function. Forty-nine (n = 49) postmyocardial infarct patients (61±10 years, EF 41±10%) underwent tissue Doppler echocardiography (TDE) and contrast enhanced magnetic resonance imaging (CMR). A15 mm×2mm sampling volume (tracked to wall motion) was placed over the long axis subepicardial region of each segment during TDE offline analysis to measure peak longitudinal systolic strain rate (SR), peak longitudinal systolic strain (PS), and CVIB. Findingswere compared with TME classified into two categories of scar thickness by CMR: Non-transmural (TME≤50%), and transmural (TME > 50%). Of 213 segments identified with resting wall motion abnormalities, 145 segments showed delayed hyperenhancement on CMR. SR, PS and CVIB were similar with no significant differences between transmural and non-transmural infarcts regardless of the echo modality.

Coronary revascularization can improve regional left ventricular volumes in non-transmural scars without contractile reserve

Revascularization (RVS) of scar segts does not lead to recovery of left ventricular (LV) function, but its effect on post-infarct remodeling is unclear. We examined the impact of RVS on regional remodeling in different transmural extents of scar (TME). Dobutamine echo (DbE) and contrast enhanced magnetic resonance imaging (ce- MRI) were performed in 72 pts post MI (age 63±10, EF 49±12%). Pts were selected for RVS (n = 31) or medical treatment (n = 41). Segts were classified as scar if there were no contractile reserve during lowdose DbE.TMEwas measured by ce-MRI; a cutoff of 75% was used to differentiate transmural (TM) from non-transmural (NT) scars. Regional end systolic (ESV) and end diastolic volumes (EDV) were measured at baseline and 12 months follow up.Of 218 segts identified as scar on DbE, 164wereNTand 54 were TM on ce-MRI. Revascularization was performed to 62 NT and 11 TM segts. In the RVS group, there was reverse remodeling with significant reduction in LV volumes in NT (ESV, 6.8±3.2 ml versus 5.8±3.7 ml, p = 0.002; EDV, 10.9±4.9 ml versus 9.8±5.6 ml, p = 0.02), but no significant change in volumes in TM (ESV, 6.9±3.7 ml versus 5.4±2.1 ml, p = 0.09; EDV, 10.2±4.4 ml versus 9.4±4.3 ml, p = 0.5). In the medically treated group, there were no changes in LV volumes in both NT (ESV, 12.0±11.9 ml versus 12.7±13.8 ml, p = 0.3; EDV, 12.5±7.8 ml versus 12.6±9.7 ml, p = 0.8) and TM (ESV, 8.0±3.8 ml versus 7.9±4.6 ml, p = 0.8; EDV, 10.3±4.8 ml versus 10.4±5.4 ml, p = 0.9). Despite absence of contractile reserve on DbE, NT benefit from coronary revascularization with regional reverse LV remodeling.

Local geometric changes in left ventricular remodeling post-infarction can lead to apparent reduction in scar thickness

Serial reduction in scar thickness has been shown in animal models. We sought whether this reduction in scar thickness may be a result of dilatation of the left ventricle (LV) with stretching and thinning of the wall. Contrast enhanced magnetic resonance imaging (CMRI) was performed to delineate radial scar thickness in 25 patients (age 63±10, 21 men) after myocardial infarction. The LV was divided into 16 segts and the absolute radial scar thickness (ST) and percentage scar to total wall thickness (%ST) were measured. Regional end diastolic (EDV) and end systolic volumes (ESV) of corresponding segments were measured on CMRI. All patients underwent revascularization and serial changes in ST, %ST, and regional volumes were assessed with a mean follow up of 15±5 months. CMRI identified a total of 93 scar segments. An increase in EDV or ESV was associated with a serial reduction inST(versusEDV, r =−0.3, p = 0.01; versusESV, r =−0.3, p = 0.005) and%ST(versusEDV, r =−0.2, p = 0.04; versus ESV, r =−0.3, p = 0.001). For segts associated with a positive increase in EDV (group I) or ESV (group II) there was a significant decrease in ST and %ST, but in those segts with stable EDV (group III) or ESV (group IV) there were no significant changes in ST and %ST (Table).

Tailoring magnetic field gradient design to magnet cryostat geometry

Eddy currents induced within a magnetic resonance imaging (MRI) cryostat bore during pulsing of gradient coils can be applied constructively together with the gradient currents that generate them, to obtain good quality gradient uniformities within a specified imaging volume over time. This can be achieved by simultaneously optimizing the spatial distribution and temporal pre-emphasis of the gradient coil current, to account for the spatial and temporal variation of the secondary magnetic fields due to the induced eddy currents. This method allows the tailored design of gradient coil/magnet configurations and consequent engineering trade-offs. To compute the transient eddy currents within a realistic cryostat vessel, a low-frequency finite-difference time-domain (FDTD) method using total-field scattered-field (TFSF) scheme has been performed and validated