L’élastographie est une technique d’imagerie qui vise à cartographier in vivo les propriétés mécaniques des tissus biologiques dans le but de fournir des informations diagnostiques additionnelles. Depuis son introduction en imagerie ultrasonore dans les années 1990, l’élastographie a trouvé de nombreuses applications. Cette modalité a notamment été utilisée pour l’étude du sein, du foie, de la prostate et des artères par imagerie ultrasonore, par résonance magnétique ou en tomographie par cohérence optique. Dans le contexte des maladies cardiovasculaires, cette modalité a un fort potentiel diagnostique puisque l’athérosclérose modifie la structure des tissus biologiques et leurs propriétés mécaniques bien avant l’apparition de tout symptôme. Quelle que soit la modalité d’imagerie utilisée, l’élastographie repose sur : l’excitation mécanique du tissu (statique ou dynamique), la mesure de déplacements et de déformations induites, et l’inversion qui permet de recouvrir les propriétés mécaniques des tissus sous-jacents. Cette thèse présente un ensemble de travaux d’élastographie dédiés à l’évaluation des tissus de l’appareil cardiovasculaire. Elle est scindée en deux parties. La première partie intitulée « Élastographie vasculaire » s’intéresse aux pathologies affectant les artères périphériques. La seconde, intitulée « Élastographie cardiaque », s’adresse aux pathologies du muscle cardiaque. Dans le contexte vasculaire, l’athérosclérose modifie la physiologie de la paroi artérielle et, de ce fait, ses propriétés biomécaniques. La première partie de cette thèse a pour objectif principal le développement d’un outil de segmentation et de caractérisation mécanique des composantes tissulaires (coeur lipidique, tissus fibreux et inclusions calciques) de la paroi artérielle, en imagerie ultrasonore non invasive, afin de prédire la vulnérabilité des plaques. Dans une première étude (Chapitre 5), nous présentons un nouvel estimateur de déformations, associé à de l’imagerie ultrarapide par ondes planes. Cette nouvelle méthode d’imagerie permet d’augmenter les performances de l’élastographie non invasive. Dans la continuité de cette étude, on propose une nouvelle méthode d’inversion mécanique dédiée à l’identification et à la quantification des propriétés mécaniques des tissus de la paroi (Chapitre 6). Ces deux méthodes sont validées in silico et in vitro sur des fantômes d’artères en polymère. Dans le contexte cardiaque, les ischémies et les infarctus causés par l’athérosclérose altèrent la contractilité du myocarde et, de ce fait, sa capacité à pomper le sang dans le corps (fonction myocardique). En échocardiographie conventionnelle, on évalue généralement la fonction myocardique en analysant la dynamique des mouvements ventriculaires (vitesses et déformations du myocarde). L’abscence de contraintes physiologiques agissant sur le myocarde (contrairement à la pression sanguine qui contraint la paroi vasculaire) ne permet pas de résoudre le problème inverse et de retrouver les propriétés mécaniques du tissu. Le terme d’élastographie fait donc ici référence à l’évaluation de la dynamique des mouvements et des déformations et non à l’évaluation des propriétés mécanique du tissu. La seconde partie de cette thèse a pour principal objectif le développement de nouveaux outils d’imagerie ultrarapide permettant une meilleure évaluation de la dynamique du myocarde. Dans une première étude (Chapitre 7), nous proposons une nouvelle approche d’échocardiographie ultrarapide et de haute résolution, par ondes divergentes, couplée à de l'imagerie Doppler tissulaire. Cette combinaison, validée in vitro et in vivo, permet d’optimiser le contraste des images mode B ainsi que l’estimation des vitesses Doppler tissulaires. Dans la continuité de cette première étude, nous proposons une nouvelle méthode d’imagerie des vecteurs de vitesses tissulaires (Chapitre 8). Cette approche, validée in vitro et in vivo, associe les informations de vitesses Doppler tissulaires et le mode B ultrarapide de l’étude précédente pour estimer l’ensemble du champ des vitesses 2D à l’intérieur du myocarde. / Elastography is an imaging technique that aims to map the in vivo mechanical properties of biological tissues in order to provide additional diagnostic information. Since its introduction in ultrasound imaging in the 1990s, elastography has found many applications. This method has been used for the study of the breast, liver, prostate and arteries by ultrasound imaging, magnetic resonance imaging (MRI) or optical coherence tomography (OCT). In the context of cardiovascular diseases (CVD), this modality has a high diagnostic potential as atherosclerosis, a common pathology causing cardiovascular diseases, changes the structure of biological tissues and their mechanical properties well before any symptoms appear. Whatever the imaging modality, elastography is based on: the mechanical excitation of the tissue (static or dynamic), the measurement of induced displacements and strains, and the inverse problem allowing the quantification of the mechanical properties of underlying tissues. This thesis presents a series of works in elastography for the evaluation of cardiovascular tissues. It is divided into two parts. The first part, entitled « Vascular elastography » focuses on diseases affecting peripheral arteries. The second, entitled « Cardiac elastography » targets heart muscle pathologies. In the vascular context, atherosclerosis changes the physiology of the arterial wall and thereby its biomechanical properties. The main objective of the first part of this thesis is to develop a tool that enables the segmentation and the mechanical characterization of tissues (necrotic core, fibrous tissues and calcium inclusions) in the vascular wall of the peripheral arteries, to predict the vulnerability of plaques. In a first study (Chapter 5), we propose a new strain estimator, associated with ultrafast plane wave imaging. This new imaging technique can increase the performance of the noninvasive elastography. Building on this first study, we propose a new inverse problem method dedicated to the identification and quantification of the mechanical properties of the vascular wall tissues (Chapter 6). These two methods are validated in silico and in vitro on polymer phantom mimicking arteries. In the cardiac context, myocardial infarctions and ischemia caused by atherosclerosis alter myocardial contractility. In conventional echocardiography, the myocardial function is generally evaluated by analyzing the dynamics of ventricular motions (myocardial velocities and deformations). The abscence of physiological stress acting on the myocardium (as opposed to the blood pressure which acts the vascular wall) do not allow the solving the inverse problem and to find the mechanical properties of the fabric. Elastography thus here refers to the assessment of motion dynamics and deformations and not to the evaluation of mechanical properties of the tissue. The main objective of the second part of this thesis is to develop new ultrafast imaging tools for a better evaluation of the myocardial dynamics. In a first study (Chapter 7), we propose a new approach for ultrafast and high-resolution echocardiography using diverging waves and tissue Doppler. This combination, validated in vitro and in vivo, optimize the contrast in B-mode images and the estimation of myocardial velocities with tissue Doppler. Building on this study, we propose a new velocity vector imaging method (Chapter 8). This approach combines tissue Doppler and ultrafast B-mode of the previous study to estimate 2D velocity fields within the myocardium. This original method was validated in vitro and in vivo on six healthy volunteers.

L'imagerie médicale ultrasonore est une modalité en perpétuelle évolution et notamment en post-traitement où il s'agit d'améliorer la résolution et le contraste des images. Ces améliorations devraient alors aider le médecin à mieux distinguer les tissus examinés améliorant ainsi le diagnostic médical. Il existe déjà une large palette de techniques "hardware" et "software". Dans ce travail nous nous sommes focalisés sur la mise en oeuvre de techniques dites de "déconvolution aveugle", ces techniques temporelles utilisant l'enveloppe du signal comme information de base. Elles sont capables de reconstruire des images parcimonieuses, c'est-à-dire des images de diffuseurs dépourvues de bruit spéculaire. Les principales étapes de ce type de méthodes consistent en i) l'estimation aveugle de la fonction d'étalement du point (PSF), ii) l'estimation des diffuseurs en supposant l'environnement exploré parcimonieux et iii) la reconstruction d'images par reconvolution avec une PSF "idéale". La méthode proposée a été comparée avec des techniques faisant référence dans le domaine de l'imagerie médicale en utilisant des signaux synthétiques, des séquences ultrasonores réelles (1D) et images ultrasonores (2D) ayant des statistiques différentes. La méthode, qui offre un temps d'exécution très réduit par rapport aux techniques concurrentes, est adaptée pour les images présentant une quantité réduite ou moyenne des diffuseurs. / The ultrasonic imaging knows a continuous advance in the aspect of increasing the resolution for helping physicians to better observe and distinguish the examined tissues. There is already a large range of techniques to get the best results. It can be found also hardware or signal processing techniques. This work was focused on the post-processing techniques of blind deconvolution in ultrasound imaging and it was implemented an algorithm that works in the time domain and uses the envelope signal as input information for it. It is a blind deconvolution technique that is able to reconstruct reflectors and eliminate the diffusive speckle noise. The main steps are: the estimation of the point spread function (PSF) in a blind way, the estimation of reflectors using the assumption of sparsity for the examined environment and the reconstruction of the image by reconvolving the sparse tissue with an ideal PSF. The proposed method was tested in comparison with some classical techniques in medical imaging reconstruction using synthetic signals, real ultrasound sequences (1D) and ultrasound images (2D) and also using two types of statistically different images. The method is suitable for images that represent tissue with a reduced amount or average scatters. Also, the technique offers a lower execution time than direct competitors.

Este trabalho trata da modelagem de sinais ultrassônicos gerados por transdutores circulares (plano e côncavo) no modo de pulso-eco, inspecionando uma região de interesse predeterminada. Para essa análise, dois modelos da resposta impulsiva de um transdutor piezelétrico foram implementados: um do transdutor circular plano e outro do transdutor circular côncavo focalizado. Este último sendo o modelo proposto pelo autor com uma geometria baseada em anéis concêntricos como elementos de área do emissor. A adição de diversos anéis concêntricos, deslocados ao longo de seu eixo e de raios sucessivamente menores, permitiu calcular o campo acústico, gerado por uma abertura côncava, bem como seu eco refletido. A resposta impulsiva de cada anel resulta da diferença entre as respostas impulsivas calculadas para um emissor circular grande e um pequeno. O modelo implementado para o cálculo dos sinais de eco foi utilizado na varredura de uma região contendo um conjunto de pontos com refletividade acústica igual a 1. A reconstrução da imagem foi feita com esses sinais aplicando-lhes a técnica de regularização de Tikhonov. A qualidade das reconstruções das imagens obtidas foi avaliada e comparada a das imagens convencionais. Dentre as métricas de avaliação das imagens reconstruídas estão a influência na variação da velocidade de propagação da onda acústica no meio, a extensão e a discretização do grid e o parâmetro de regularização alfa. Todas as reconstruções foram analisadas segundo o Erro Médio Quadrático (MSE). Finalmente, ensaios experimentais foram conduzidos para a obtenção de A-scans (imagens em modo de amplitude) os quais foram inseridas no modelo teórico para a reconstrução de imagens e analisadas pelo MSE. / This work deals with the modeling of ultrasonic signals generated by circular transducers (planar and concave pistons) in pulse-echo mode, inspecting a predetermined region of interest (ROI). For this analysis, two models of the impulse response of a transducer were implemented: one using a plane piston transducer and the other, a model proposed by the author of this work, using a concave transducer with concentric rings as elements of the emitting area. The addition of several concentric rings moved along its axis allowed us to calculate the acoustic field generated by a concave opening and the echo reflected from each point in space. The impulse response of each ring represents the difference between the impulse responses calculated for a large circular transmitter and a small one. The model implemented for calculating the echo signals is used to scan a region, within a ROI, containing a set of points with acoustic reflectivity of 1. Simulations of the regions are made with these signals by applying the Tikhonov regularization method. To evaluate the quality of image reconstruction, the images are compared with the conventional images. Among the metrics to evaluate the reconstructed images are the influence of the variation of the acoustic wave propagation in the media, the grid range and discretization and the parameter of regularization alpha. All of the image reconstructions were analyzed through the Mean Square Error (MSE) criterion. Finally, experiments were conducted in order to obtain A-scans which were then re-inserted in the theoretical model to reconstruct and analyze the images.

Depuis le milieu du XXème siècle, les techniques d’imagerie fonctionnelles ont un rôle grandissant dans notre compréhension sur les fonctions du cerveau en conditions physiologique et pathologique. Bien que l’IRMf fasse partie des techniques les plus communément utilisées pour l’imagerie du cerveau complet lors d’études préclinique et clinique, cette modalité souffre de sa résolution spatiotemporelle et sa sensibilité pour enregistrer finement les fonctions et activités cérébrales. Récemment l’imagerie fonctionnelle par ultrason (ifUS) a subi des développements permettant d’être complémentaires à l’IRMf ainsi qu’aux autres techniques d’imagerie cérébrales classiquement employées. Contrairement aux ultrasons focalisés conventionnels, l’imagerie hémodynamique proposé par l’ifUS repose sur une illumination ultrasonore plane permettant la détection des globules rouges en mouvement et la mesure de leur vitesse dans les micro-vaisseaux cérébraux. De ce fait, l’ifUS est indirectement lié à l’activité cérébrale d’où l’importance d’une meilleure compréhension des mécanismes du couplage neuro-vasculaire liant l’activité neuronale et les variations cérébrales d’apport en sang. De plus, cette technique a le potentiel pour fournir des informations précises sur les processus de certaines pathologies à la fois sur des modèles précliniques et chez l’homme. Dans un premier temps, j’exposerais mes travaux sur les récents développements techniques permettant l’ifUS in vivo (i) en condition chronique, (ii) sur l’animal éveillé, libre de mouvement et effectuant une tache comportementale et (iii) des vaisseaux capillaires chez le rongeur et l’homme. Dans un second temps, je démontrerais que l’ifUS in vivo peut fournir des informations nouvelles sur des pathologies telles que l’accident vasculaire cérébrale. / Since the middle of the 20th century, functional imaging technologies are making an increasing impact on our understanding on brain functions in both physiological and pathological conditions. Even if fMRI is nowadays one of the most used tool for whole brain imaging in pre-clinical and clinical studies, it lacks sufficient spatiotemporal resolution and sensitivity to assess fine brain function and activity. Functional ultrasound imaging (fUSi) has been recently developed and presents a potential to complement fMRI and other existing brain imaging modalities. Contrary to conventional ultrasound using focus beams, fUSi relies on hemodynamic imaging based on ultrasound plane-wave illumination to detect red blood cells movement and velocity in brain micro-vessels. Consequently, the fUSi signal is indirectly related to brain activity and it is therefore important to better understand the mechanisms of the neurovascular coupling linking neural activity and cerebral blood changes. Here again, fUSi may provide relevant information about disease processes in preclinical models but also in humans. First, I will present recent technical developments allowing in vivo fUSi (i) in chronic condition, (ii) in freely moving and behaving rats and (iii) in rodents and human brain capillaries. Second, I will demonstrate how fUSi could provide new insights in brain pathologies such as stroke.

L’imagerie échographique en trois dimensions (3D) est une modalité d’imagerie médicale en plein développement. En plus de ses nombreux avantages (faible cout, absence de rayonnement ionisant, portabilité) elle permet de représenter les structures anatomiques dansleur forme réelle qui est toujours 3D. Les sondes à balayage mécaniques, relativement lentes, tendent à être remplacées par des sondes bidimensionnelles ou matricielles qui sont unprolongement dans les deux directions, latérale et azimutale, de la sonde classique 1D. Cetagencement 2D permet un dépointage du faisceau ultrasonore et donc un balayage 3D del’espace. Habituellement, les éléments piézoélectriques d’une sonde 2D sont alignés sur unegrille et régulièrement espacés d’une distance (en anglais le « pitch ») soumise à la loi del’échantillonnage spatial (distance inter-élément inférieure à la demi-longueur d’onde) pour limiter l’impact des lobes de réseau. Cette contrainte physique conduit à une multitude d’éléments de petite taille. L’équivalent en 2D d’une sonde 1D de 128 éléments contient128x128=16 384 éléments. La connexion d’un nombre d’éléments aussi élevé constitue unvéritable défi technique puisque le nombre de canaux dans un échographe actuel n’excède querarement les 256. Les solutions proposées pour contrôler ce type de sonde mettent en oeuvredu multiplexage ou des techniques de réduction du nombre d’éléments, généralement baséessur une sélection aléatoire de ces éléments (« sparse array »). Ces méthodes souffrent dufaible rapport signal à bruit du à la perte d’énergie qui leur est inhérente. Pour limiter cespertes de performances, l’optimisation reste la solution la plus adaptée. La première contribution de cette thèse est une extension du « sparse array » combinéeavec une méthode d’optimisation basée sur l’algorithme de recuit simulé. Cette optimisation permet de réduire le nombre nécessaire d’éléments à connecter en fonction des caractéristiques attendues du faisceau ultrasonore et de limiter la perte d’énergie comparée à la sonde complète de base. La deuxième contribution est une approche complètement nouvelle consistant à adopter un positionnement hors grille des éléments de la sonde matricielle permettant de supprimer les lobes de réseau et de s’affranchir de la condition d’échantillonnage spatial. Cette nouvelles tratégie permet d’utiliser des éléments de taille plus grande conduisant ainsi à un nombre d’éléments nécessaires beaucoup plus faible pour une même surface de sonde. La surface active de la sonde est maximisée, ce qui se traduit par une énergie plus importante et donc unemeilleure sensibilité. Elle permet également de balayer un angle de vue plus important, leslobes de réseau étant très faibles par rapport au lobe principal. Le choix aléatoire de la position des éléments et de leur apodization (ou pondération) reste optimisé par le recuit simulé.Les méthodes proposées sont systématiquement comparées avec la sonde complète dansle cadre de simulations numériques dans des conditions réalistes. Ces simulations démontrent un réel potentiel pour l’imagerie 3D des techniques développées. Une sonde 2D de 8x24=192 éléments a été construite par Vermon (Vermon SA, ToursFrance) pour tester les méthodes de sélection des éléments développées dans un cadreexpérimental. La comparaison entre les simulations et les résultats expérimentaux permettentde valider les méthodes proposées et de prouver leur faisabilité. / 3D Ultrasound imaging is a fast-growing medical imaging modality. In addition to its numerous advantages (low cost, non-ionizing beam, portability) it allows to represent the anatomical structures in their natural form that is always three-dimensional. The relativelyslow mechanical scanning probes tend to be replaced by two-dimensional matrix arrays that are an extension in both lateral and elevation directions of the conventional 1D probe. This2D positioning of the elements allows the ultrasonic beam steering in the whole space. Usually, the piezoelectric elements of a 2D array probe are aligned on a regular grid and spaced out of a distance (the pitch) subject to the space sampling law (inter-element distancemust be shorter than a mid-wavelength) to limit the impact of grating lobes. This physical constraint leads to a multitude of small elements. The equivalent in 2D of a 1D probe of 128elements contains 128x128 = 16,384 elements. Connecting such a high number of elements is a real technical challenge as the number of channels in current ultrasound scanners rarely exceeds 256. The proposed solutions to control this type of probe implement multiplexing or elements number reduction techniques, generally using random selection approaches (« spars earray »). These methods suffer from low signal to noise ratio due to the energy loss linked to the small number of active elements. In order to limit the loss of performance, optimization remains the best solution. The first contribution of this thesis is an extension of the « sparse array » technique combined with an optimization method based on the simulated annealing algorithm. The proposed optimization reduces the required active element number according to the expected characteristics of the ultrasound beam and permits limiting the energy loss compared to the initial dense array probe.The second contribution is a completely new approach adopting a non-grid positioningof the elements to remove the grating lobes and to overstep the spatial sampling constraint. This new strategy allows the use of larger elements leading to a small number of necessaryelements for the same probe surface. The active surface of the array is maximized, whichresults in a greater output energy and thus a higher sensitivity. It also allows a greater scansector as the grating lobes are very small relative to the main lobe. The random choice of the position of the elements and their apodization (or weighting coefficient) is optimized by the simulated annealing.The proposed methods are systematically compared to the dense array by performing simulations under realistic conditions. These simulations show a real potential of the developed techniques for 3D imaging.A 2D probe of 8x24 = 192 elements was manufactured by Vermon (Vermon SA, Tours,France) to test the proposed methods in an experimental setting. The comparison between simulation and experimental results validate the proposed methods and prove their feasibility. / L'ecografia 3D è una modalità di imaging medicale in rapida crescita. Oltre ai vantaggiin termini di prezzo basso, fascio non ionizzante, portabilità, essa permette di rappresentare le strutture anatomiche nella loro forma naturale, che è sempre tridimensionale. Le sonde ascansione meccanica, relativamente lente, tendono ad essere sostituite da quelle bidimensionali che sono una estensione in entrambe le direzioni laterale ed azimutale dellasonda convenzionale 1D. Questo posizionamento 2D degli elementi permette l'orientamentodel fascio ultrasonico in tutto lo spazio. Solitamente, gli elementi piezoelettrici di una sondamatriciale 2D sono allineati su una griglia regolare e separati da una distanza (detta “pitch”) sottoposta alla legge del campionamento spaziale (la distanza inter-elemento deve esseremeno della metà della lunghezza d'onda) per limitare l'impatto dei lobi di rete. Questo vincolo fisico porta ad una moltitudine di piccoli elementi. L'equivalente di una sonda 1D di128 elementi contiene 128x128 = 16.384 elementi in 2D. Il collegamento di un così grandenumero di elementi è una vera sfida tecnica, considerando che il numero di canali negliecografi attuali supera raramente 256. Le soluzioni proposte per controllare questo tipo disonda implementano le tecniche di multiplazione o la riduzione del numero di elementi, utilizzando un metodo di selezione casuale (« sparse array »). Questi metodi soffrono di unbasso rapporto segnale-rumore dovuto alla perdita di energia. Per limitare la perdita di prestazioni, l’ottimizzazione rimane la soluzione migliore. Il primo contributo di questa tesi è un’estensione del metodo dello « sparse array » combinato con un metodo di ottimizzazione basato sull'algoritmo del simulated annealing. Questa ottimizzazione riduce il numero degli elementi attivi richiesto secondo le caratteristiche attese del fascio di ultrasuoni e permette di limitare la perdita di energia.Il secondo contributo è un approccio completamente nuovo, che propone di adottare un posizionamento fuori-griglia degli elementi per rimuovere i lobi secondari e per scavalcare il vincolo del campionamento spaziale. Questa nuova strategia permette l'uso di elementi piùgrandi, riducendo così il numero di elementi necessari per la stessa superficie della sonda. La superficie attiva della sonda è massimizzata, questo si traduce in una maggiore energia equindi una maggiore sensibilità. Questo permette inoltre la scansione di un più grande settore,in quanto i lobi secondari sono molto piccoli rispetto al lobo principale. La scelta casualedella posizione degli elementi e la loro apodizzazione viene ottimizzata dal simulate dannealing. I metodi proposti sono stati sistematicamente confrontati con la sonda completaeseguendo simulazioni in condizioni realistiche. Le simulazioni mostrano un reale potenzialedelle tecniche sviluppate per l'imaging 3D.Una sonda 2D di 8x24 = 192 elementi è stata fabbricata da Vermon (Vermon SA, ToursFrance) per testare i metodi proposti in un ambiente sperimentale. Il confronto tra lesimulazioni e i risultati sperimentali ha permesso di convalidare i metodi proposti edimostrare la loro fattibilità.

L'échographie est aujourd'hui une technique d'imagerie de diagnostic répandue. Si l'imagerie dite ‘classique' basée sur la réponse linéaire des tissus est couramment utilisée, l'imagerie harmonique, basée sur la réponse non linéaire des tissus, est maintenant aussi utilisée en routine clinique. L'estimation du paramètre de non linéarité d'un milieu par une technique ultrasonore amène de nouvelles perspectives en termes d'imagerie et de diagnostic. Cependant, la méthode de mesure du paramètre de non linéarité est limitée par deux facteurs, la présence du speckle et la concentration de l'énergie à une profondeur donnée (la zone focale). Cette thèse a pour objectifs de répondre aux deux limitations mentionnées précédemment en proposant de nouvelles méthodes de lissage de l'image pour réduire le speckle et d'améliorer l'estimation du paramètre de non linéarité en mode écho par de nouvelles méthodes d'émission. Dans un premier temps, il a été proposé d'utiliser une méthode de filtrage spatiale basée sur des filtres orthogonaux (filtres de Thomson) lors de la formation de voie en réception pour lisser le speckle. Ce filtrage spatiale intervient après la transmission d'ondes planes sous différents angles pour améliorer la résolution spatiale et le contraste tout en accélérant la cadence d'imagerie. Dans un deuxième temps, l'estimation du paramètre de non linéarité est faite avec une méthode comparative. Le champ de pression du second harmonique d'une zone de référence est comparé avec le champ de pression d'une zone dont le paramètre de non linéarité est inconnu. Cependant, dans le cas des images échographiques, le champ de pression du second harmonique n'est pas accessible. Nous faisons l'hypothèse que la pression acoustique locale est liée à l'intensité de l'image échographique si le speckle est réduit et lissé. La transmission d'ondes planes et l'application de filtres orthogonaux permet de mieux délimiter le paramètre de non linéarité par rapport à une transmission focalisée / Nowadays, ultrasound imaging is a common diagnostic tool thanks to its non-invasive behavior and relatively cheap equipment. Classic medical echographic imaging is based on the linear response of the biological tissue. However harmonic imaging, based on the harmonic frequencies generated by the nonlinear properties of the tissue, is more and more used for clinical application. The quantification of nonlinearity is based on the evaluation of the nonlinearity parameter which strongly influences the harmonics generation. The nonlinearity parameter estimation using an echographic approach would bring new modalities for imaging and diagnosis. However the echographic method for nonlinearity estimation is limited by two factors: the presence of speckles in the image and the focalization used during transmission, which concentrates the energy at one particular depth. The objectives of this thesis work are developing novel approaches to reduce the speckle noise using original smoothing techniques and improving the nonlinearity parameter estimation in echo mode using new transmission-reception strategies. Firstly, new speckle noise reduction approaches were investigated. The Thomson’s multitaper approach was proposed, consisting in using several different orthogonal apodizations during beamforming. This method was combined to a coherent plane-wave compounding transmission-reception strategy improving the spatial resolution and the contrast while improving the frame rate. In a second time, the nonlinearity parameter was estimated using a comparative method. The second-harmonic pressure field of a reference area was compared to the pressure field of an area where the nonlinearity parameter is unknown. However in echo-mode, the pressure field of the medium is unknown. It is assumed in this thesis work that the local pressure can be derived from envelope image local amplitude if the speckle noise is smoothed. The nonlinearity parameter estimation has been improved using plane-wave transmission and orthogonal apodizations compared to the use of a single focalization transmission / Oggigiorno, le tecniche di imaging ad ultrasuoni sono un comune strumento di diagnosi, grazie alla loro non invasività e alla relativa economicità dei sistemi. La risposta lineare dei tessuti biologici è la base per le tecniche di imaging ecografico tradizionali. La generazione di frequenze ad armoniche superiori da parte dei tessuti può essere sfruttata per sviluppare tecniche di imaging innovative (i.e., imaging armonico), che sono sempre più utilizzate per applicazioni cliniche. Tali tecniche sono basate sul metodo di valutazione del parametro di non linearità che influenza fortemente la generazione delle armoniche all’interno dei tessuti. I metodi per la stima dei suddetti parametri sfruttano solitamente un approccio ecografico tradizionale. Di conseguenza, gli effetti legati alla focalizzazione impiegata durante la trasmissione, che concentra l’energia ad una particolare profondità, e la presenza di speckle nell’immagine finale, rendono più incerta la stima del parametro di non linearità. In questa tesi sono proposti metodi innovativi finalizzati a due scopi: ridurre, nelle immagini, il rumore dovuto a speckle, tramite l’adozione di nuove tecniche di smoothing; migliorare la stima dei parametri di non linearità, tramite l’impiego di nuove strategie di beamforming in trasmissione e ricezione. Per ridurre il rumore dovuto a speckle, è stato proposto un approccio di filtraggio spaziale basato sull’impiego dei filtri di Thomson. Tale tecnica consiste nell’impiego di numerose apodizzazioni ortogonali fra di loro in fase di beamforming. Il metodo è stato in particolare combinato con la tecnica di imaging coherent plane-wave compounding, con lo scopo di migliorare la risoluzione spaziale e il contrasto e, al contempo, incrementare il frame rate. Il parametro di non linearità è stato misurato tramite un approccio comparativo. Il campo di pressione della seconda armonica in un’area di riferimento dell’immagine è stato confrontato con quello di un’area in cui il parametro di non linearità è ignoto. In questa tesi, grazie alla riduzione del rumore speckle, è stato possibile assumere che il campo di pressione

L'objectif de cette thèse était d'explorer de nouvelles applications cardiaques pour l'histotripsie et de développer les outils permettant leur mise en place non-invasive. La thérapie ultrasonore cardiaque est en effet encore assez peu développée aujourd’hui, à cause de la difficulté à traiter un organe en mouvement permanent, et très bien protégé derrière la cage thoracique.Nous avons d'abord montré in vivo, sur un modèle ovin, que l’on pouvait sectionner les cordages mitraux de manière non-invasive ainsi que traiter la sténose aortique calcifiée. Engendrer de la cavitation sur les feuillets valvulaires permet effectivement d’agir à distance sur les calcifications, et de globalement assouplir la valve.Simultanément, nous avons développé un dispositif pour la thérapie cardiaque non invasive, fondé sur le concept de cavité à retournement temporel. Ce dispositif permet l'émission d'impulsions ultrasonores de haute intensité dans un très grand volume d’intérêt. L’on peut ainsi déplacer le point de thérapie en 3 dimensions de manière entièrement électronique, et sans déplacer mécaniquement l’appareil. Après optimisation, ce dispositif a permis de créer des lésions mécaniques bien contrôlées dans une région d'intérêt de 2 000 cm3.Pour faire face au défi que représente la cage thoracique, nous avons développé une méthode de focalisation adaptative et l'avons mise en œuvre dans un prototype 2D du dispositif. Avec cette méthode, nous pouvons non seulement construire un front d'onde ultrasonore adaptatif qui se propage de manière préférentielle à travers les espaces intercostaux, mais grâce aux propriétés des cavités à retournement temporel, nous pouvons également augmenter la pression focale obtenue sur la cible de thérapie.Enfin, pour approfondir ce travail sur la focalisation adaptative, et nous avons considéré le cas de l'imagerie transcrânienne. Pour cette application, nous avons choisi d’utiliser la focalisation par retournement temporel dans le bruit de speckle, pour corriger les aberrations induites par le crâne. En simulations numériques, nous avons pu calculer les modulations de phase et d'amplitude induites par les os et améliorer le contraste et la résolution d'une image B-mode. / The objective of this thesis was to explore new applications for cardiac histotripsy, and to develop the tools making it possible non-invasively. Cardiac ultrasound therapy indeed still remains limited due to the tremendous challenge of treating a constantly and rapidly moving organ, well protected behind the ribcage.We first showed in vivo, on a large animal model, that histotripsy could be used non-invasively to cut mitral chordae, and to treat calcified aortic stenosis in a beating heart. Cavitation on the valve leaflets can indeed locally and remotely act on the calcifications, and globally soften the valve. Simultaneously, we developed a therapeutic device allowing completely non-invasive cardiac shock-wave therapy based on the time reversal cavity concept. In particular, this device allows the emission of high intensity ultrasound pulses, and provides 3D electronical steering of the therapy focal spot in a large volume. After a thorough optimisation process, this device was capable of creating well controlled mechanical lesions over a 2 000 cm3 region of interest. To tackle the challenge of ultrasound propagation through the rib cage, we developed an adaptive focusing method (DORT method through a time reversal cavity), and implemented it in a 2D prototype of the device. With this method, we not only could build an adaptive ultrasonic wavefront propagating preferentially through the intercostal spaces, but due to time reversal cavities properties, we could also increase the peak pressure obtained on target.Finally, we pushed our work on adaptive focusing further, and considered the case of transcranial imaging. For this application, we chose to use the time reversal of speckle noise technique, to correct the aberrations induced by the skull. In numerical simulations, we were able to derive the phase and amplitude modulations induced by the bones, and could improve the contrast and resolution of a B-mode image.

Vascular diseases are among the leading causes of death in Canada and around the globe. A major underlying cause of most such medical conditions is atherosclerosis, a gradual accumulation of plaque on the walls of blood vessels. Particularly vulnerable to atherosclerosis is the carotid artery, which carries blood to the brain. Dangerous narrowing of the carotid artery can lead to embolism, a dislodgement of plaque fragments which travel to the brain and are the cause of most strokes. If this pathology can be detected early, such a deadly scenario can be potentially prevented through treatment or surgery. This not only improves the patient's prognosis, but also dramatically lowers the overall cost of their treatment. Medical imaging is an indispensable tool for early detection of atherosclerosis, in particular since the exact location and shape of the plaque need to be known for accurate diagnosis. This can be achieved by locating the plaque inside the artery and measuring its volume or texture, a process which is greatly aided by image segmentation. In particular, the use of ultrasound imaging is desirable because it is a cost-effective and safe modality. However, ultrasonic images depict sound-reflecting properties of tissue, and thus suffer from a number of unique artifacts not present in other medical images, such as acoustic shadowing, speckle noise and discontinuous tissue boundaries. A robust ultrasound image segmentation technique must take these properties into account. Prior to segmentation, an important pre-processing step is the extraction of a series of features from the image via application of various transforms and non-linear filters. A number of such features are explored and evaluated, many of them resulting in piecewise smooth images. It is also proposed to decompose the ultrasound image into several statistically distinct components. These components can be then used as features directly, or other features can be obtained from them instead of the original image. The decomposition scheme is derived using Maximum-a-Posteriori estimation framework and is efficiently computable. Furthermore, this work presents and evaluates an algorithm for segmenting the carotid artery in 3D ultrasound images from other tissues. The algorithm incorporates information from different sources using an energy minimization framework. Using the ultrasound image itself, statistical differences between the region of interest and its background are exploited, and maximal overlap with strong image edges encouraged. In order to aid the convergence to anatomically accurate shapes, as well as to deal with the above-mentioned artifacts, prior knowledge is incorporated into the algorithm by using weak geometric priors. The performance of the algorithm is tested on a number of available 3D images, and encouraging results are obtained and discussed.

Vascular diseases are among the leading causes of death in Canada and around the globe. A major underlying cause of most such medical conditions is atherosclerosis, a gradual accumulation of plaque on the walls of blood vessels. Particularly vulnerable to atherosclerosis is the carotid artery, which carries blood to the brain. Dangerous narrowing of the carotid artery can lead to embolism, a dislodgement of plaque fragments which travel to the brain and are the cause of most strokes. If this pathology can be detected early, such a deadly scenario can be potentially prevented through treatment or surgery. This not only improves the patient's prognosis, but also dramatically lowers the overall cost of their treatment. Medical imaging is an indispensable tool for early detection of atherosclerosis, in particular since the exact location and shape of the plaque need to be known for accurate diagnosis. This can be achieved by locating the plaque inside the artery and measuring its volume or texture, a process which is greatly aided by image segmentation. In particular, the use of ultrasound imaging is desirable because it is a cost-effective and safe modality. However, ultrasonic images depict sound-reflecting properties of tissue, and thus suffer from a number of unique artifacts not present in other medical images, such as acoustic shadowing, speckle noise and discontinuous tissue boundaries. A robust ultrasound image segmentation technique must take these properties into account. Prior to segmentation, an important pre-processing step is the extraction of a series of features from the image via application of various transforms and non-linear filters. A number of such features are explored and evaluated, many of them resulting in piecewise smooth images. It is also proposed to decompose the ultrasound image into several statistically distinct components. These components can be then used as features directly, or other features can be obtained from them instead of the original image. The decomposition scheme is derived using Maximum-a-Posteriori estimation framework and is efficiently computable. Furthermore, this work presents and evaluates an algorithm for segmenting the carotid artery in 3D ultrasound images from other tissues. The algorithm incorporates information from different sources using an energy minimization framework. Using the ultrasound image itself, statistical differences between the region of interest and its background are exploited, and maximal overlap with strong image edges encouraged. In order to aid the convergence to anatomically accurate shapes, as well as to deal with the above-mentioned artifacts, prior knowledge is incorporated into the algorithm by using weak geometric priors. The performance of the algorithm is tested on a number of available 3D images, and encouraging results are obtained and discussed.

Este trabalho trata da modelagem de sinais ultrassônicos gerados por transdutores circulares (plano e côncavo) no modo de pulso-eco, inspecionando uma região de interesse predeterminada. Para essa análise, dois modelos da resposta impulsiva de um transdutor piezelétrico foram implementados: um do transdutor circular plano e outro do transdutor circular côncavo focalizado. Este último sendo o modelo proposto pelo autor com uma geometria baseada em anéis concêntricos como elementos de área do emissor. A adição de diversos anéis concêntricos, deslocados ao longo de seu eixo e de raios sucessivamente menores, permitiu calcular o campo acústico, gerado por uma abertura côncava, bem como seu eco refletido. A resposta impulsiva de cada anel resulta da diferença entre as respostas impulsivas calculadas para um emissor circular grande e um pequeno. O modelo implementado para o cálculo dos sinais de eco foi utilizado na varredura de uma região contendo um conjunto de pontos com refletividade acústica igual a 1. A reconstrução da imagem foi feita com esses sinais aplicando-lhes a técnica de regularização de Tikhonov. A qualidade das reconstruções das imagens obtidas foi avaliada e comparada a das imagens convencionais. Dentre as métricas de avaliação das imagens reconstruídas estão a influência na variação da velocidade de propagação da onda acústica no meio, a extensão e a discretização do grid e o parâmetro de regularização alfa. Todas as reconstruções foram analisadas segundo o Erro Médio Quadrático (MSE). Finalmente, ensaios experimentais foram conduzidos para a obtenção de A-scans (imagens em modo de amplitude) os quais foram inseridas no modelo teórico para a reconstrução de imagens e analisadas pelo MSE. / This work deals with the modeling of ultrasonic signals generated by circular transducers (planar and concave pistons) in pulse-echo mode, inspecting a predetermined region of interest (ROI). For this analysis, two models of the impulse response of a transducer were implemented: one using a plane piston transducer and the other, a model proposed by the author of this work, using a concave transducer with concentric rings as elements of the emitting area. The addition of several concentric rings moved along its axis allowed us to calculate the acoustic field generated by a concave opening and the echo reflected from each point in space. The impulse response of each ring represents the difference between the impulse responses calculated for a large circular transmitter and a small one. The model implemented for calculating the echo signals is used to scan a region, within a ROI, containing a set of points with acoustic reflectivity of 1. Simulations of the regions are made with these signals by applying the Tikhonov regularization method. To evaluate the quality of image reconstruction, the images are compared with the conventional images. Among the metrics to evaluate the reconstructed images are the influence of the variation of the acoustic wave propagation in the media, the grid range and discretization and the parameter of regularization alpha. All of the image reconstructions were analyzed through the Mean Square Error (MSE) criterion. Finally, experiments were conducted in order to obtain A-scans which were then re-inserted in the theoretical model to reconstruct and analyze the images.