Les techniques ultrasonores offrent de nombreux avantages, à la fois par leur utilisation facile et la sécurité du patient. De plus, la recherche, visant à étendre les possibles champs d’applications, est particulièrement active. Cependant, l’accès à des équipements adaptés et supportant des logiciels est conditio sine qua non pour l’expérimentation de nouvelles techniques. Ce projet de thèse traite des problématiques de traitement du signal et d'image dans un contexte d'imagerie médicale et vise à répondre à deux objectifs scientifiques: le premier consiste à contribuer au développement d'une puissante plateforme de recherche (ULA-OP) alors que le second a pour objectif d'introduire et de valider, grâce à cette plateforme, des méthodes de traitement non-standard qui ne pourraient pas être évaluées avec des équipements médicaux commerciaux. ULA-OP, un équipement recherche qui donne accès aux développeurs à une grande liberté de contrôle et de configuration l’ensemble des parties actives du système, de la transmission aux traitements des signaux échographiques. Il offre aussi la possibilité d’accéder aux signaux bruts à n’importe quel niveau de la chaîne de réception. Durant cette thèse, les capacités du système ont été améliorées en implémentant des outils logiciels comme des simulateurs de champ acoustique (propagation linéaire et non-linéaire), et en développant des programmes de génération de signaux post-échographique. L’ULA-OP a été crucial pour développer et tester différentes techniques non-standard telles qu’un schéma adaptatif de formation de voie et une méthode d’imagerie Doppler couleur/vecteur, qui seront détaillés dans le manuscrit. En particulier, une nouvelle méthode a été développée pour des applications d’élastographie quasi-statique. Cette méthode, basée sur un algorithme d’estimation du mouvement dans le domaine fréquentiel et combinée à une méthode d’imagerie haute fréquence, a permis d’améliorer la qualité des élastogrammes obtenus. Cette nouvelle méthode a d’abord été testée in-vitro par des traitements hors ligne des signaux reçus et pour ensuite être implémentée en temps réel sur le ULA-OP. Les résultats obtenus montrent que cette technique est performante et que les élastogrammes présentent une qualité supérieure comparée à ceux obtenues avec les méthodes connues de la littérature / Ultrasound techniques offer many advantages, in terms of both ease of realization and patients’ safety. The research aimed at expanding the fields of application, is nowadays particularly active. The availability of suitable hardware and supporting software tools is condicio sine qua non for the experimentation of new techniques. This Ph.D project addresses signal/image processing issues in medical ultrasound and seeks to achieve two major scientific goals: the first is to contribute to the development of a powerful ultrasound research platform (ULA¬OP), while the second is introducing and validating, through this platform, non-standard methods which could not be tested with commercial equipment. ULA-OP is a research system, which gives developers great freedom in terms of management and control of every section, from signal transmission to echo-signal processing; it also offers the possibility to access raw data at any point in the receive chain. During the thesis, the capabilities of the system were improved by creating advanced software tools, such as acoustic field simulators (for linear and nonlinear propagation), and by developing echo-signals post-elaboration programs. ULA-OP was crucial to develop and test various non-standard techniques such as an adaptive beamforming scheme and a color/vector Doppler imaging method, which will be detailed in this thesis. In particular, a novel technique was developed for quasi-static elastography applications. This technique, based on a frequency domain displacement estimation algorithm, combined with a high-frame-rate averaging method, aims at improving the quality of the elastograms. The new method was first tested in-vitro by offline processing the received signals, and then it was implemented in real-time on ULA-OP. The results show that this technique is effective and that the obtained elastograms present higher quality compared with those obtained with standard algorithms

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'imagerie photoacoustique est une modalité d'imagerie fonctionnelle basée sur la génération d'ondes acoustiques par des tissus soumis à une illumination optique (impulsion laser). L'utilisation de différentes longueurs d'ondes optiques permet la discrimination des milieux imagés. Cette modalité est prometteuse pour de nombreuses applications médicales liées, par exemple, à la croissance, au vieillissement et à l'évolution de la vascularisation des tissus. En effet, l'accès à l'oxygénation du sang dans les tissus est rendu possible par l'imagerie photoacoustique. Cela permet, entre autres applications, la discrimination de tumeurs bénignes ou malignes et la datation de la mort tissulaire (nécrose). Ce travail de thèse a pour objectif principal la construction d'une chaîne de traitement des données photoacoustiques multispectrales pour le calcul de l'oxygénation du sang dans les tissus. Les principales étapes sont, d'une part, la discrimination des données (clustering), pour extraire les zones d'intérêt, et d'autre part, la quantification des différents constituants présents dans celles-ci (unmixing). Plusieurs méthodes non supervisées de discrimination et de quantification ont été développées et leurs performances comparées sur des données photoacoustiques multispectrales expérimentales. Celles-ci ont été acquises sur la plateforme photoacoustique du laboratoire, lors de collaborations avec d'autres laboratoires et également sur un système commercial. Pour la validation des méthodes développées, de nombreux fantômes contenant différents absorbeurs optiques ont été conçus. Lors du séjour de cotutelle de thèse en Italie, des modes d'imagerie spécifiques pour l'imagerie photoacoustique 2D et 3D temps-réel ont été développés sur un échographe de recherche. Enfin, des acquisitions in vivo sur modèle animal (souris) au moyen d'un système commercial ont été réalisées pour valider ces développements. / Photoacoustic imaging is a functional technique based on the creation of acoustic waves from tissues excited by an optical source (laser pulses). The illumination of a region of interest, with a range of optical wavelengths, allows the discrimination of the imaged media. This modality is promising for various medical applications in which growth, aging and evolution of tissue vascularization have to be studied. Thereby, photoacoustic imaging provides access to blood oxygenation in biological tissues and also allows the discrimination of benign or malignant tumors and the dating of tissue death (necrosis). The present thesis aims at developing a multispectral photoacoustic image processing chain for the calculation of blood oxygenation in biological tissues. The main steps are, first, the data discrimination (clustering), to extract the regions of interest, and second, the quantification of the different media in these regions (unmixing). Several unsupervised clustering and unmixing methods have been developed and their performance compared on experimental multispectral photoacoustic data. They were acquired on the experimental photoacoustic platform of the laboratory, during collaborations with other laboratories and also on a commercial system. For the validation of the developed methods, many phantoms containing different optical absorbers have been produced. During the co-supervision stay in Italy, specific imaging modes for 2D and 3D real-time photoacoustic imaging were developed on a research scanner. Finally, in vivo acquisitions using a commercial system were conducted on animal model (mouse) to validate these developments.

L'objectif de cette thèse était d'explorer le potentiel de l’imagerie du cerveau humain par ultrasons. L'anatomie, le flux sanguin et la rigidité des tissus mous ont déjà été étudiés avec l'imagerie ultrasonore ultrarapide chez l'homme et validés sur plusieurs organes, tels que le sein et le foie, mais pas encore sur le cerveau adulte. La principale limitation de l'imagerie échographique transcrânienne est aujourd'hui le très fort artefact d'aberration induit par le crâne. En effet, l’os, de par sa composition ne permet pas la propagation des ultrasons comme ailleurs dans le corps humain. Dans cette thèse, nous avons utilisé l'imagerie ultrasonore ultrarapide pour l'évaluation de la rigidité des tissus mous et l'imagerie neurofonctionnelle dans le cerveau humain adulte, lors de chirurgies du cerveau afin de contourner dans un premier temps le problème des aberrations induites par le crâne. La dernière partie de cette thèse était axée sur la correction d’aberration pour l’échographie quantitative et l’imagerie ultrasonore transcrânienne. Nous avons tout d’abord fourni plusieurs preuves de l'intérêt d'utiliser l’élastographie par onde de cisaillement pendant la chirurgie du cerveau. Nous avons également présenté notre nouvelle technique d’élastographie par onde de cisaillement en 3D à l'aide d'une sonde matricielle dans le but de pouvoir dépasser les limitations du 2D et notamment être moins dépendant de l’opérateur.Dans un second volet, nous avons démontré la capacité des ultrasons ultrarapides à identifier, cartographier et différencier en profondeur les régions d'activation corticales en réponse à un stimulus, à la fois chez les patients éveillés et chez les patients anesthésiés. Nous avons démontré que l'imagerie neurofonctionnelle par ultrasons a le potentiel de devenir une modalité complète de neuroimagerie avec des avantages majeurs pour une utilisation peropératoire.Dans un troisième volet, nous avons utilisé une technique en trois étapes pour calculer précisément la vitesse du son (SSE) dans un milieu. Cette technique a été testée dans des fantômes ultrasonores et in vivo dans les foies de patients. Dans les deux cas, notre méthode a été capable de trouver la vitesse du son correspondant au milieu. Nous avons démontré que la SSE était liée à la fraction de graisse. Cette analyse a permis de conclure que la SSE était en mesure de distinguer un foie sain et d’un foie malade aussi bien avec la biopsie qu’avec l’IRM comme méthode de référence. Combiné à l'utilisation de la formule de Wood, nous avons même pu avoir accès à une fraction de graisse mesurée par ultrasons de manière non invasive. Puis nous avons combiné la correction d’aberration de phase, d'amplitude et de vitesse du son pour faire de l’imagerie transcrânienne en simulation numérique. Nous avons atteint notre objectif en obtenant des images représentant fidèlement le milieu (position latérale et profondeur) et caractérisées par une résolution et un contraste similaires à ceux obtenus avec une source ponctuelle dans le milieu / The objective of this thesis was to explore the potential of human brain ultrasound imaging. Anatomy, blood flow and soft tissue stiffness have already been studied with ultrafast ultrasound imaging in humans and validated in several organs, such as, the breast and liver but not yet on the adult brain. The main limitation of transcranial ultrasound imaging is today the very strong skull-induced aberration artefact. Indeed, the bone, due to its composition, does not allow for ultrasound propagation as elsewhere in the human body. Therefore, this thesis was focused on the development of ultrafast ultrasound imaging for the evaluation of soft tissue stiffness and neurofunctional imaging in the adult human brain, during brain surgery to bypass the problem of skull aberration, and on an aberration correction technique for transcranial ultrasound imaging.We first provided several evidence of the benefit of using shear wave elastography during brain surgery. We also presented our new technique for 3D shear wave elastography using a matrix array in order to be able to overcome the limitations of 2D imaging and in particular to reduce the operator dependence.In a second phase, we demonstrated the capability of ultrasound to identify, map and differentiate in depth cortical regions of activation in response to a stimulus, both in awake patients and in anaesthetized patients. We have demonstrated that ultrasound neurofunctional imaging has the potential to become a comprehensive modality of neuroimaging with major benefits for intraoperative use. In a third part, we developed a new sound speed estimation (SSE) technique, based on a three-step technique that estimates the sound speed accurately corresponding to the illuminated medium. This technique was tested in ultrasound phantoms and in vivo in patient’s liver. In both cases, our method was able to find the sound speed corresponding to the medium. We demonstrated that SSE was related to the fat fraction. This analysis led to the conclusion that SSE was able to distinguish a healthy liver from a diseased liver with both biopsy and MRI as gold standard. Combined with the use of the Wood’s formula, we were even able to access a fat fraction measured by non-invasive ultrasound. Finally, by combining the phase, the amplitude and the sound speed estimation, we have developed a new aberration correction algorithm to perform transcranial ultrasound imaging. By performing numerical simulations, we obtained images that faithfully represented the medium (lateral position and depth) and characterized by one resolution and one contrast similar to those obtained with a punctual source in the medium

Cette thèse est consacrée au développement d'outils de traitement pour l'Imagerie par Résonance Magnétique du Tenseur de Diffusion (IRM-TD). Cette technique d'IRM récente est d'une grande importance pour comprendre le fonctionnement du cerveau ou pour améliorer le diagnostic de pathologies neurologiques. Nous proposons des méthodes de traitement basées sur la géométrie Riemannienne, les équations aux dérivées partielles et les techniques de propagation de front. La première partie de ce travail est théorique. Après des rappels sur le système nerveux humain, l'IRM et la géométrie différentielle, nous étudions l'espace des lois normales multivariées. L'introduction d'une structure Riemannienne sur cet espace nous permet de définir des statistiques et des schémas numériques intrinsèques qui sont à la base des algorithmes proposés dans la seconde partie. Les propriétés de cet espace sont importantes pour l'IRM-TD car les tenseurs de diffusion sont les matrices de covariance de lois normales modélisant la diffusion des molécules d'eau en chaque voxel du milieu imagé. La seconde partie est méthodologique. Nous y introduisons des approches originales pour l'estimation et la régularisation d'IRM-TD. Puis nous montrons comment évaluer le degré de connectivité entre aires corticales et introduisons un modèle statistique d'évolution de surface permettant de segmenter ces images. Finalement, nous proposons une méthode de recalage non-rigide. La dernière partie de cette thèse est consacrée à l'analyse des connexions entre le cortex cérébral et les noyaux gris centraux, impliquées dans des tâches motrices, et à l'étude du réseau anatomo-fonctionnel du cortex visuel humain..

Le travail présenté dans ce mémoire de thèse s'intéresse au problème de l'imagerie microonde bidimensionnelle polarimétrique. Il est scindé en deux parties, la première est dédiée aux études expérimentales liées à la détection d'objets enfouis tandis que la seconde présente le développement d'un algorithme d'imagerie quantitative. Ainsi le travail effectué dans la première partie du manuscrit porte sur la réalisation d'antennes ultra-large bande et la détermination de leur comportement lors de mesures polarimétriques sur sites. Dans ce dernier cas, plusieurs tests ont été réalisés et permettent de conclure sur les qualités des capteurs développés. Dans une deuxième phase, un code d'imagerie bidimensionnelle prenant en compte à la fois la polarisation 2D-TM et 2D-TE est présenté. Le développement d'un algorithme de diffraction électromagnétique est tout d'abord exposé. Son association avec un algorithme de problème inverse existant déjà au laboratoire permet alors de générer le code d'imagerie quantitative. Plusieurs configurations de reconstruction d'objet homogène, isotrope ou non, permettent ensuite de mettre en évidence ses capacités.

Depuis 1995, environ 350 exoplanètes ont été détectées mais seule une dizaine l'a été directement pour plusieurs raisons. D'abord, les projections de l'étoile hôte et de sa planète sur la sphère céleste sont très proches -- quelques fractions de secondes d'arc. Ceci impose un diamètre minimum pour le télescope et un système compensant les perturbations atmosphériques. Ensuite, le flux lumineux de l'étoile hôte est entre un million et dix milliards de fois plus fort que celui de la planète. Pour réduire ce flux stellaire sans affecter le flux planétaire, nous utilisons un coronographe dont les performances sont limitées par les défauts optiques qui doivent être compensés ou estimés. Pendant ma thèse, j'ai étudié sur des plans théoriques et expérimentaux deux techniques: le coronographe à quatre quadrants à étages multiples (MFQPM) pour atténuer le flux stellaire et la self coherent camera (SCC) qui minimise l'impact des aberrations optiques en utilisant l'incohérence entre lumières stellaires et planétaires. J'ai montré en laboratoire qu'un prototype non optimisé du MFQPM fournissait une extinction achromatique de l'étoile centrale de mille à dix milles entre 550 et 750nm. Puis, j'ai montré par simulations que la SCC s'associe aisément avec un coronographe de type Lyot et qu'en utilisant ses deux modes de fonctionnement -- analyseur de surface d'onde en plan focal et imagerie différentielle --, des planètes de type Terre pouvaient être détectées de l'espace sous des conditions réalistes. J'ai également mesuré un défaut réel de phase via la SCC sur le banc d'Imagerie Très Haute Dynamique que j'ai développé à l'Observatoire de Paris.

Les imageurs embarqués sur les appareils de traitement par radiothérapie sont des dispositifs très efficaces pour améliorer la précision géométrique des irradiations. Cependant, ils ne sont pertinents que s'ils font l'objet de contrôles de performances précis et réguliers. Notre propos a donc été de développer des solutions logicielles permettant d'automatiser l'analyse des images de contrôle de qualité des principaux modes utilisés en imagerie de repositionnement : le mode 2D-MV qui mesure l'image radiante bidimensionnelle haute énergie transmise par le patient pendant le traitement, les modes 2D-kV (images bidimensionnelles basse énergie) et mode 3D (images tridimensionnelles MV ou kV) qui permettent de réajuster la position du patient avant de traiter. Nous avons également automatisé l'analyse des images du test de Winston & Lutz utilisé pour évaluer la précision mécanique des appareils de traitement sur lesquels s'exercent des contraintes supplémentaires dues au poids supplémentaire de l'imageur embarqué. Enfin nous avons développé et mis en oeuvre une méthodologie originale et performante utilisant des objetstest numériques pour évaluer les performances des solutions logicielles mises au point. Au bilan, l'automatisation des contrôles de qualité des imageurs embarqués avec les solutions logicielles développées dans le cadre de ces travaux de thèse permet de réduire d'un facteur de 100 le temps consacré par l'équipe de physique médicale à l'analyse des résultats des contrôles tout en améliorant leur précision grâce à l'utilisation d'analyses objectives et reproductibles et leur traçabilité grâce à l'édition automatique de rapports de contrôle et d'études statistiques.

Ce travail présente les résultats obtenus dans le cadre de développements méthodologiques en imagerie par résonance magnétique. Tout d'abord, la mise en place de la technique d'imagerie rapide utilisant un échantillonnage de l'espace k le long d'une trajectoire spirale à densité variable est exposée. L'optimisation des paramètres d'acquisition et la comparaison de différentes techniques de reconstruction ont été menées via des simulations numériques. Une approche originale de la calibration de la trajectoire dans l'espace k est proposée. L'imagerie spirale a été utilisée pour mettre en œuvre une méthode de mesure de la perméabilité de la BHE au Gd-DOTA. Ce protocole a été intégré aux mesures de volume sanguin et d'index de taille des vaisseaux utilisant le Sinerem®. Les résultats obtenus montrent des différences au niveau des paramètres de la microvascularisation sur les modèles de tumeurs C6 et RG2. La présence du Sinerem® conduit à une diminution des valeurs de constante d'échange à travers la paroi vasculaire (Ktrans) dans la tumeur de 24% en moyenne. Cette étude a également mis en évidence l'extravasation du Sinerem®, pendant les deux premières heures après injection du produit, uniquement dans les tumeurs RG2.

Cette these traite de la modelisation statistique de formes en imagerie<br />cerebrale.<br /><br />Dans une premiere partie, nous proposons un modele statistique de la forme des sillons corticaux. Le modele est bati par apprentissage a partir de sillons extraits d'images IRM et dotes d'une representation parametrique. La definition d'un repere intrinseque a la forme sillon permet d'aligner l'ensemble des formes extraites et de construire une population d'apprentissage coherente sur laquelle appliquer une analyse en composantes principales afin de deriver le modele. Ce modele statistique est ensuite etendu a un graphe de sillons afin de decrire non plus seulement les caracteristiques morphologiques d'un sillon, mais aussi les relations de position et d'orientation entre sillons principaux.<br />L'analyse presentee ici porte sur un sous-graphe defini par un couple de sillons. <br /><br />Dans une seconde partie, trois applications de la modelisation proposee sont envisagees. D'une part, nous l'utilisons dans un cadre d'evaluation de methodes de recalage global inter-sujets. Pratiquee sur des amers locaux, l'analyse statistique fournit un indicateur de la similarite des formes au sein des populations recalees, et produit un critere de comparaison entre les methodes. D'autre part, nous exploitons la connaissance statistique apportee par le modele sur les sillons dans le contexte de la construction d'atlas anatomiques et fonctionnels. Nous proposons une methode locale et non-lineaire de recalage inter-sujets de donnees fonctionnelles, exprimees sous forme de dipoles MEG (localisations d'activations fonctionnelles), base sur la modelisation des amers anatomiques que sont les sillons corticaux. Exprimentee sur une population de 18 sujets, cette methode s'est averee apte a reduire la variabilite fonctionnelle inter-individuelle observee. Enfin, nous appliquons la methodologie proposee dans le cas des sillons a la modelisation statistique de la forme de frontieres fonctionnelles delimitant des aires visuelles de bas-niveau.