Iterative simulations to estimate the elastic properties from a series of MRI images followed by MRI-US validation

The modeling of breast deformations is of interest in medical applications such as image-guided biopsy, or image registration for diagnostic purposes. In order to have such information, it is needed to extract the mechanical properties of the tissues. In this work, we propose an iterative technique based on finite element analysis that estimates the elastic modulus of realistic breast phantoms, starting from MRI images acquired in different positions (prone and supine), when deformed only by the gravity force. We validated the method using both a single-modality evaluation in which we simulated the effect of the gravity force to generate four different configurations (prone, supine, lateral, and vertical) and a multi-modality evaluation in which we simulated a series of changes in orientation (prone to supine). Validation is performed, respectively, on surface points and lesions using as ground-truth data from MRI images, and on target lesions inside the breast phantom compared with the actual target segmented from the US image. The use of pre-operative images is limited at the moment to diagnostic purposes. By using our method we can compute patient-specific mechanical properties that allow compensating deformations

A Proof of Concept of a Non-Invasive Image-Based Material Characterization Method for Enhanced Patient-Specific Computational Modeling

PURPOSE: Computational models of cardiovascular structures rely on their accurate mechanical characterization. A validated method able to infer the material properties of patient-specific large vessels is currently lacking. The aim of the present study is to present a technique starting from the flow-area (QA) method to retrieve basic material properties from magnetic resonance (MR) imaging. METHODS: The proposed method was developed and tested, first, in silico and then in vitro. In silico, fluid-structure interaction (FSI) simulations of flow within a deformable pipe were run with varying elastic modules (E) between 0.5 and 32 MPa. The proposed QA-based formulation was assessed and modified based on the FSI results to retrieve E values. In vitro, a compliant phantom connected to a mock circulatory system was tested within MR scanning. Images of the phantom were acquired and post-processed according to the modified formulation to infer E of the phantom. Results of in vitro imaging assessment were verified against standard tensile test. RESULTS: In silico results from FSI simulations were used to derive the correction factor to the original formulation based on the geometrical and material characteristics. In vitro, the modified QA-based equation estimated an average E = 0.51 MPa, 2% different from the E derived from tensile tests (i.e. E = 0.50 MPa). CONCLUSION: This study presented promising results of an indirect and non-invasive method to establish elastic properties from solely MR images data, suggesting a potential image-based mechanical characterization of large blood vessels

Novel Ultrasound Elastography Imaging System for Breast Cancer Assessment

Abstract Most conventional methods of breast cancer screening such as X-ray, Ultrasound (US) and MRI have some issues ranging from weaknesses associated with tumour detection or classification to high cost or excessive time of image acquisition and reconstruction. Elastography is a non- invasive technique to visualize suspicious areas in soft tissues such as the breast, prostate and myocardium using tissue stiffness as image contrast mechanism. In this study, a breast Elastography system based on US imaging is proposed. This technique is fast, expected to be cost effective and more sensitive and specific compared to conventional US imaging. Unlike current Elastography techniques that image relative elastic modulus, this technique is capable of imaging absolute Young\u27s modulus (YM). In this technique, tissue displacements and surface forces used to mechanically stimulate the tissue are acquired and used as input to reconstruct the tissue YM distribution. For displacements acquisition, two techniques were used in this research: 1) a modified optical flow technique, which estimates the displacement of each node from US pre- and post-compression images and 2) Radio Frequency (RF) signal cross-correlation technique. In the former, displacements are calculated in 2 dimensions whereas in the latter, displacements are calculated in the US axial direction only. For improving the quality of elastography images, surface force data was used to calculate the stress distribution throughout the organ of interest by using an analytical model and a statistical numerical model. For force data acquisition, a system was developed in which load cells are used to measure forces on the surface of the breast. These forces are input into the stress distribution models to estimate the tissue stress distribution. By combining the stress field with the strain field calculated from the acquired displacements using Hooke\u27s law, the YM can be reconstructed efficiently. To validate the proposed technique, numerical and tissue mimicking phantom studies were conducted. For the numerical phantom study, a 3D breast-shape phantom was created with synthetic US pre- and post-compression images where the results showed the feasibility of reconstructing the absolute value of YM of tumour and background. In the tissue mimicking study, a block shape gelatine- agar phantom was constructed with a cylindrical inclusion. Results obtained from this study also indicated reasonably accurate reconstruction of the YM. The quality of the obtained elasticity images shows that image quality is improved by incorporating the adapted stress calculation techniques. Furthermore, the proposed elastography system is reasonably fast and can be potentially used in real-time clinical applications

Position-based simulation of deformations for autonomous robotic ultrasound scanning

Realistic and fast simulation of anatomical deformations due to ultrasound probe pressure is of outstanding importance for testing and validation of autonomous robotic ultrasound systems. We propose a deformation model which relies on the position-based dynamics (PBD) approach to simulate the probetissue interaction and predict the displacement of internal targets during US acquisition. Performances of the patient-specific PBD anatomical model are evaluated in comparison to two different simulations relying on the traditional finite element (FE) method, in the context of breast ultrasound scanning. Localization error obtained when applying the PBD model remains below 11 mm for all the tumors even for input displacements in the order of 30 mm. The proposed method is able to achieve a better trade-off among accuracy, computation time and generalization capabilities with respect to the two FE models. Position-based dynamics approach has proved to be successful in modeling breast tissue deformations during US acquisition. It represents a valid alternative to classical FE methods for simulating the interaction between US probe and tissues

Biomechanical modelling of probe to tissue interaction during ultrasound scanning

Purpose: Biomechanical simulation of anatomical deformations caused by ultrasound probe pressure is of outstanding importance for several applications, from the testing of robotic acquisition systems to multi-modal image fusion and development of ultrasound training platforms. Different approaches can be exploited for modelling the probe-tissue interaction, each achieving different trade-offs among accuracy, computation time and stability. Methods: We assess the performances of different strategies based on the finite element method for modelling the interaction between the rigid probe and soft tissues. Probe\u2013tissue contact is modelled using (i) penalty forces, (ii) constraint forces, and (iii) by prescribing the displacement of the mesh surface nodes. These methods are tested in the challenging context of ultrasound scanning of the breast, an organ undergoing large nonlinear deformations during the procedure. Results: The obtained results are evaluated against those of a non-physically based method. While all methods achieve similar accuracy, performance in terms of stability and speed shows high variability, especially for those methods modelling the contacts explicitly. Overall, prescribing surface displacements is the approach with best performances, but it requires prior knowledge of the contact area and probe trajectory. Conclusions: In this work, we present different strategies for modelling probe\u2013tissue interaction, each able to achieve different compromises among accuracy, speed and stability. The choice of the preferred approach highly depends on the requirements of the specific clinical application. Since the presented methodologies can be applied to describe general tool\u2013tissue interactions, this work can be seen as a reference for researchers seeking the most appropriate strategy to model anatomical deformation induced by the interaction with medical tools

Bridging spatiotemporal scales in biomechanical models for living tissues : from the contracting Esophagus to cardiac growth

Appropriate functioning of our body is determined by the mechanical behavior of our organs. An improved understanding of the biomechanical functioning of the soft tissues making up these organs is therefore crucial for the choice for, and development of, efficient clinical treatment strategies focused on patient-specific pathophysiology. This doctoral dissertation describes the passive and active biomechanical behavior of gastrointestinal and cardiovascular tissue, both in the short and long term, through computer models that bridge the cell, tissue and organ scale. Using histological characterization, mechanical testing and medical imaging techniques, virtual esophagus and heart models are developed that simulate the patient-specific biomechanical organ behavior as accurately as possible. In addition to the diagnostic value of these models, the developed modeling technology also allows us to predict the acute and chronic effect of various treatment techniques, through e.g. drugs, surgery and/or medical equipment. Consequently, this dissertation offers insights that will have an unmistakable impact on the personalized medicine of the future.Het correct functioneren van ons lichaam wordt bepaald door het mechanisch gedrag van onze organen. Een verbeterd inzicht in het biomechanisch functioneren van deze zachte weefsels is daarom van cruciale waarde voor de keuze voor, en ontwikkeling van, efficiënte klinische behandelingsstrategieën gefocust op de patiënt-specifieke pathofysiologie. Deze doctoraatsthesis brengt het passieve en actieve biomechanisch gedrag van gastro-intestinaal en cardiovasculair weefsel, zowel op korte als lange termijn, in kaart via computermodellen die een brug vormen tussen cel-, weefsel- en orgaanniveau. Aan de hand van histologische karakterisering, mechanische testen en medische beeldvormingstechnieken worden virtuele slokdarm- en hartmodellen ontwikkeld die het patiënt-specifieke orgaangedrag zo accuraat mogelijk simuleren. Naast de diagnostische waarde van deze modellen, laat de ontwikkelde modelleringstechnologie ook toe om het effect van verschillende behandelingstechnieken, via medicatie, chirurgie en/of medische apparatuur bijvoorbeeld, acuut en chronisch te voorspellen. Bijgevolg biedt deze doctoraatsthesis inzichten die een onmiskenbare impact zullen hebben op de gepersonaliseerde geneeskunde van de toekomst

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

In vivo morphometric and mechanical characterization of trabecular bone from high resolution magnetic resonance imaging

La osteoporosis es una enfermedad ósea que se manifiesta con una menor densidad ósea y el deterioro de la arquitectura del hueso esponjoso. Ambos factores aumentan la fragilidad ósea y el riesgo de sufrir fracturas óseas, especialmente en mujeres, donde existe una alta prevalencia. El diagnóstico actual de la osteoporosis se basa en la cuantificación de la densidad mineral ósea (DMO) mediante la técnica de absorciometría dual de rayos X (DXA). Sin embargo, la DMO no puede considerarse de manera aislada para la evaluación del riesgo de fractura o los efectos terapéuticos. Existen otros factores, tales como la disposición microestructural de las trabéculas y sus características que es necesario tener en cuenta para determinar la calidad del hueso y evaluar de manera más directa el riesgo de fractura. Los avances técnicos de las modalidades de imagen médica, como la tomografía computarizada multidetector (MDCT), la tomografía computarizada periférica cuantitativa (HR-pQCT) y la resonancia magnética (RM) han permitido la adquisición in vivo con resoluciones espaciales elevadas. La estructura del hueso trabecular puede observarse con un buen detalle empleando estas técnicas. En particular, el uso de los equipos de RM de 3 Teslas (T) ha permitido la adquisición con resoluciones espaciales muy altas. Además, el buen contraste entre hueso y médula que proporcionan las imágenes de RM, así como la utilización de radiaciones no ionizantes sitúan a la RM como una técnica muy adecuada para la caracterización in vivo de hueso trabecular en la enfermedad de la osteoporosis. En la presente tesis se proponen nuevos desarrollos metodológicos para la caracterización morfométrica y mecánica del hueso trabecular en tres dimensiones (3D) y se aplican a adquisiciones de RM de 3T con alta resolución espacial. El análisis morfométrico está compuesto por diferentes algoritmos diseñados para cuantificar la morfología, la complejidad, la topología y los parámetros de anisotropía del tejido trabecular. En cuanto a la caracterización mecánica, se desarrollaron nuevos métodos que permiten la simulación automatizada de la estructura del hueso trabecular en condiciones de compresión y el cálculo del módulo de elasticidad. La metodología desarrollada se ha aplicado a una población de sujetos sanos con el fin de obtener los valores de normalidad del hueso esponjoso. Los algoritmos se han aplicado también a una población de pacientes con osteoporosis con el fin de cuantificar las variaciones de los parámetros en la enfermedad y evaluar las diferencias con los resultados obtenidos en un grupo de sujetos sanos con edad similar.Los desarrollos metodológicos propuestos y las aplicaciones clínicas proporcionan resultados satisfactorios, presentando los parámetros una alta sensibilidad a variaciones de la estructura trabecular principalmente influenciadas por el sexo y el estado de enfermedad. Por otra parte, los métodos presentan elevada reproducibilidad y precisión en la cuantificación de los valores morfométricos y mecánicos. Estos resultados refuerzan el uso de los parámetros presentados como posibles biomarcadores de imagen en la enfermedad de la osteoporosis.Alberich Bayarri, Á. (2010). In vivo morphometric and mechanical characterization of trabecular bone from high resolution magnetic resonance imaging [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/8981Palanci

An Efficient Hemodynamic Workflow in Computational Surgery

For few decades, it has been shown that atherosclerosis is the cause of the majority of clinical cardiovascular diseases including peripheral arterial diseases. The diagnosis and treatment for vascular disease has evolved significantly over the past years considering the rapid advances in imaging technologies. In recent years, computational fluid dynamics has been increasingly used as a simulation tool for blood flows. Numerous researches connect wall shear stress quantities to endovascular diseases such as stenosis, aneurism, and atherosclerosis. A thorough knowledge of vascular anatomy and hemodynamic would be beneficial for understanding the development and progression of the disease, the therapeutic decision process and follow up. The objective of this dissertation is to propose a computational fluid dynamic framework that includes: Understanding how streamline efficiently hemodynamic simulation for main arteries to produce database for clinical study/Providing some confidence estimate on numerical results/Extending the state of the art of clinical study by including motion and particles analysis.Computer Science, Department o