Brevet FR3000556A1 — publié le 20140704

DOMAINE TECHNIQUE DE L'INVENTION La présente invention concerne en général les techniques de résonance magnétique nucléaire (RMN) utilisées notamment pour l'analyse d'un échantillon. La présente invention s'applique avantageusement au domaine de la médecine pour l'analyse d'échantillons biologiques du corps humain par exemple ou au domaine de la chimie pour l'analyse de compositions chimiques. Elle s'avère particulièrement avantageuse pour corriger la phase et/ou la fréquence d'un signal de décroissance libre d'induction (FID). Dans une application avantageuse mais non limitative, l'invention sera utilisée pour recombiner efficacement des signaux RMN détectés par des réseaux d'antennes. ÉTAT DE LA TECHNIQUE Découverte dans les années quarante, la résonance magnétique nucléaire (RMN) a pris une importance considérable dans de nombreux domaines scientifiques pour l'analyse de matériaux et particulièrement en médecine pour l'observation in vivo du corps humain. De très nombreuses applications de cette découverte ont été mises en oeuvre. Dans le domaine médical, il s'agit notamment de l'imagerie par résonance magnétique ou IRM et de la spectroscopie RMN qui permet d'analyser la composition des tissus du corps humain. La résonance magnétique nucléaire est un phénomène complexe qui n'est décrit complètement que par les lois de la physique quantique. Pour ce qui est de la compréhension de l'invention et de l'état de la technique il est cependant suffisant de connaître les phénomènes principaux mis en oeuvre par la RMN et très brièvement rappelés ci-après. On notera ici qu'il existe une très abondante littérature technique sur tous les aspects de la RMN. Placés dans un fort champ magnétique, traditionnellement noté Bo, généralement de l'ordre du tesla (T) et plus, certains atomes acquièrent des propriétés particulières telles que s'ils sont perturbés par une onde radio de fréquence adaptée ils renvoient alors un signal que l'on peut analyser. Les atomes concernés sont ceux qui possèdent un noyau atomique ayant un moment magnétique ou spin non nul. Dans le domaine médical, l'hydrogène (H) de l'eau (H2O) qui constitue la majeure partie du corps humain possède cette propriété. La molécule d'eau possède deux noyaux d'hydrogène ou protons ayant un spin non nul animés d'un mouvement de précession de leur axe à une fréquence particulière dite fréquence de Larmor. Placés dans un fort champ magnétique statique extérieur les noyaux voient alors leur spin s'orienter soit dans la direction du champ magnétique extérieur, soit dans la direction opposée tout en subissant un mouvement de précession à la fréquence de Larmor dépendant de l'intensité du champ extérieur Bo. Pour les noyaux d'hydrogène, la fréquence de Larmor s'accroît de 42 MHz (mégahertz = 106 hertz) par tesla. Elle est par exemple proche de 4 MHz à 0,1 T ; 64 MHz à 1,5 T et 600 MHz à 14 T. Les noyaux peuvent alors être perturbés par une onde radio ou radiofréquence (RF) d'autant plus facilement que la fréquence de celle-ci est proche de la fréquence de Larmor et donc de préférence en résonance avec cette dernière. La perturbation a pour effet de changer l'orientation du spin jusqu'à possiblement le faire basculer dans l'autre direction permise par le champ magnétique extérieur Bo. C'est-à-dire que l'orientation est parallèle ou antiparallèle à celui-ci et passe ainsi de l'une à l'autre. Traditionnellement noté Bi, l'induction magnétique de l'onde radio perturbatrice est appliquée perpendiculairement à celle de Bo pour faire basculer le spin des noyaux. C'est ce basculement et le retour à l'équilibre ou relaxation qui génère le signal électrique détecté dans tous les équipements RMN. Ce signal est habituellement qualifié de FID de l'anglais « free induction decay » c'est-à-dire « décroissance libre de l'induction ». La perturbation électromagnétique Bi est appliquée à l'aide de bobines électriques le plus souvent qualifiées d'antennes. Les bobines ou antennes servent en général aussi bien à appliquer la perturbation qu'à détecter les signaux électriques engendrés lors du retour à l'équilibre des noyaux. Pour ce faire, elles sont associées d'une part à des circuits électroniques qui vont permettre d'y faire circuler pendant un temps contrôlé un courant électrique à la fréquence de Larmor, ce qui va générer dans la zone du corps à analyser l'onde électromagnétique et le champ perturbateur B1, et d'autre part à des circuits électroniques qui vont permettre de détecter les signaux électriques engendrés par le retour à l'équilibre des noyaux après que l'on a interrompu la perturbation électromagnétique. On notera ici que les signaux électriques qu'il faut pouvoir détecter avec ces antennes sont généralement faibles voire très faibles et que l'enjeu est de pouvoir les distinguer du niveau de bruit intrinsèque à l'équipement électrique utilisé. Le rapport signal sur bruit (RSB) est le critère clé de la partie détection électrique de tout équipement RMN. Il doit être suffisamment élevé pour que les signaux détectés le soient à un niveau au dessus du bruit qui les rend exploitables par l'application considérée. Comme cela sera décrit en détail par la suite en référence à la figure 1, les techniques connues consistent à effectuer une transformée de Fourier du signal complexe de FID pour obtenir une répartition spectrale de l'échantillon. Cette répartition spectrale permet de comparer la fréquence et l'amplitude d'un signal de référence par rapport à celles des signaux des métabolites à caractériser. Le signal de référence correspond le plus souvent à une espèce majoritaire dans l'échantillon. Pour l'analyse d'un échantillon biologique, le signal de référence est typiquement le signal de l'eau. La comparaison de la fréquence et de l'amplitude du signal de l'eau par rapport à celles des signaux des métabolites à caractériser permet d'identifier la nature et la teneur des métabolites dans l'échantillon. Cette technique est globalement satisfaisante. Elle requiert néanmoins que le signal de référence acquis soit correctement positionné, en phase et en fréquence par rapport à une position théorique. On dit alors que le signal de référence doit être « en résonance ». Pour cela, il faut régler la fréquence de conversion du récepteur, typiquement une bobine, exactement sur la fréquence de résonance du signal de référence. Or, de nombreux paramètres tendent à faire apparaître un glissement en fréquence qui éloigne le signal de référence acquis de sa position théorique. Comme cela sera détaillé par la suite, ce glissement est au moins en partie dû à la différence qui existe en pratique entre la fréquence de démodulation utilisée par le récepteur et la fréquence effective du signal de référence. Ainsi, les techniques connues nécessitent de corriger manuellement la phase et la fréquence du signal complexe de FID avant son exploitation. Ces étapes de correction sont relativement longues et fastidieuses. Par ailleurs, une autre opération de correction manuelle consiste à rephaser le signal de référence. En effet, en pratique il y a forcément un décalage aléatoire entre la phase du signal utilisé par le récepteur et la phase du signal de référence. Cela se traduit, sur une représentation d'un signal complexe de FID, par le fait que le premier point de la FID n'est en général pas réel. En outre, dans les cas où l'on acquiert des signaux provenant de plusieurs entités volumétriques de l'échantillon, habituellement désignées voxels, l'inhomogénéité du champ BO tend également à entraîner des déphasages et des décalages en fréquence. Par ailleurs, des déphasages et décalages en fréquence sont également générés par l'utilisation de multiples antennes. En effet, au cours de décennies de développement les structures d'antennes ont beaucoup évoluées et la recherche et le développement des équipements RMN tendent maintenant à favoriser l'emploi, pour l'analyse d'une zone donnée du corps humain, d'antennes multiples. Ces antennes multiples sont généralement identiques et de petites tailles et on combine les signaux détectés séparément dans le but d'obtenir un signal composite qui a notamment un meilleur RSB que s'il avait été obtenu avec une seule antenne couvrant la même zone d'analyse. Par ailleurs, l'utilisation d'une pluralité d'antennes permet généralement d'augmenter la rapidité de l'acquisition, tout en conservant un bon RSB, ce qui s'avère primordiale lorsque l'échantillon fait partie d'un corps vivant et donc bien souvent non parfaitement immobile. Ce type d'antennes dites « antennes en réseau phasé » traduction de l'anglais « phased array antenna » se caractérise par une complexité notablement accrue de la partie électronique et du traitement informatique des signaux produits séparément par chacune des antennes du réseau. Il faut en effet pouvoir les recombiner efficacement en dépit des inévitables différences géométriques, physiques et électriques qui existent entre antennes et qu'il faut pouvoir compenser pour effectivement obtenir les améliorations attendues du réseau d'antennes. Le document US 2004/0095139 décrit une solution pour recombiner des signaux de FID provenant de plusieurs bobines. Un objectif de la présente invention est donc de proposer une méthode qui permette de simplifier et d'améliorer la correction de phase et/ou de 5 fréquence des signaux de RMN détectés par des antennes tels que des bobines. Il serait en outre particulièrement avantageux, de proposer une solution pour améliorer la recombinaison des signaux RMN détectés séparément par une ou plusieurs antennes sans pour autant complexifier lourdement les 10 équipements connus. Un autre avantage consiste à proposer une solution pour simplifier la complexité des équipements de RMN, notamment ceux utilisant plusieurs bobines réceptrices. Les autres objets, caractéristiques et avantages de la présente 15 invention apparaîtront à l'examen de la description suivante et des dessins d'accompagnement. Il est entendu que d'autres avantages peuvent être incorporés. RÉSUMÉ DE L'INVENTION 20 Selon un mode de réalisation, l'invention porte sur un procédé de correction de phase et/ou de fréquence d'au moins un signal de décroissance libre d'induction (FID) se rapportant à au moins un échantillon à analyser par résonance magnétique nucléaire (RMN), l'échantillon comprenant au moins une espèce à caractériser ainsi qu'une espèce de référence prise parmi l'eau ou un 25 solvant, la teneur de l'espèce de référence dans le solvant étant au moins supérieure à deux fois la teneur de l'au moins une espèce à caractériser, le procédé comprenant les étape suivantes : a. application d'au moins un champ Bo constant sur l'au moins un échantillon, 30 b. acquisition par une ou plusieurs antennes d'un ou plusieurs signaux complexes de décroissance libre d'induction (FID), chaque signal complexe de FID comprenant une partie réelle et une partie imaginaire, l'étape d'acquisition étant effectuée de manière à ce que dans chaque signal complexe de FID l'amplitude du signal de l'espèce de référence soit de préférence au moins deux fois supérieur à l'amplitude de celui de l'au moins une espèce à caractériser. Le procédé comprend également au moins l'étape suivante : c. pour chaque signal complexe de FID, génération d'un module de FID obtenu par extraction du module de chaque signal complexe de FID, la phase de chaque signal complexe de FID étant ainsi remise à 0. Selon un mode de réalisation privilégié, le procédé comprend une étape ultérieure d'application d'une transformée de Fourier à un signal prenant en compte le module de FID. Cette étape est optionnelle. Elle permet de générer une représentation spectrale de l'échantillon. Selon un autre mode de réalisation on conserve le signal dans le domaine temporel. Cet autre mode de réalisation sera par exemple retenu pour générer une carte métabolique par CSI. Les signaux complexes de décroissance libre d'induction (FID) comprennent chacun des données relatives à l'espèce de référence. Ainsi, avant génération du module de FID, le procédé selon l'invention conserve le signal de l'espèce de référence. De préférence, l'invention ne prévoit donc pas de suppression de l'espèce de référence avant extraction du module de la FID. Une suppression partielle est possible, mais il faut néanmoins que l'amplitude du signal de l'espèce de référence soit au moins supérieure à deux fois l'amplitude du signal de l'espèce à caractériser avant extraction du module de la FID. Le signal de l'espèce de référence fait ainsi office de porteuse pour les signaux des espèces à caractériser. Ainsi, avec l'invention, à partir du moment où l'on acquiert un signal de l'espèce de référence, typiquement le signal de l'eau, on peut se passer des parties réelles et imaginaires. Ensuite, après avoir corrigé la phase du signal par prise en compte du module, on peut supprimer ou atténuer fortement le signal de l'espèce de référence. En appliquant la méthode selon l'invention, on obtient un spectre symétrique de part et d'autre de la fréquence correspondant à celle de l'espèce majoritaire, typiquement de l'eau lorsque l'échantillon est un tissu biologique. De manière surprenante, bien qu'une partie des informations captées par les antennes (l'information de phase) soit éliminée lors de la prise en compte du module de la FID, le ratio signal sur bruit du signal (RSB) final demeure parfaitement satisfaisant. Dans le cadre de développement de la présente invention, on s'attendait en effet à ce que cette perte d'information conduise à une baisse trop pénalisante du RSB. L'obtention d'un RSB très satisfaisant provient vraisemblablement du fait que la perte d'information est largement compensée par l'élimination automatique des décalages en phase et en fréquence que permet la séquence d'étapes selon l'invention. L'invention permet en effet de corriger automatiquement les décalages en phase et en fréquence dus notamment aux mouvements de l'échantillon durant la durée de la FID, aux inhomogénéités du champ BO et aux inévitables variations de comportement des différentes antennes dans le cadre d'une acquisition par plusieurs antennes. Ainsi, le spectre obtenu permet d'identifier et de caractériser très aisément les espèces de l'échantillon. De manière particulièrement avantageuses, l'invention permet également de se passer d'étapes supplémentaires spécialement dédiées aux corrections de phase et de fréquence, ces étapes nécessitant habituellement un temps de traitement significatif, voire bien souvent l'intervention manuelle d'un opérateur. En outre, la quantité de données à traiter est moitié de celles qu'il faut traiter de façon conventionnelle où l'on a besoin du signal complexe pour calculer le déphasage du signal. De plus le calcul du module peut être effectué dans le domaine réel puisqu'il suffit de prendre la racine carré de la somme de deux nombres, chacun élevé au carré, alors que le calcul de la phase nécessite de travailler avec des nombres complexes, ce qui est plus difficile à mettre en oeuvre. Le temps et le coût du traitement sont donc considérablement réduits, et le fait que l'intervention manuelle d'un opérateur ne soit plus requise permet d'effectuer un traitement automatique en temps réel sans impacter les utilisateurs. L'invention offre ainsi des avantages significatifs en termes de précision, de fiabilité, de reproductibilité, de temps de traitement et de réduction de coûts. Selon un autre mode de réalisation, l'invention porte sur un procédé d'analyse par résonance magnétique nucléaire (RMN) d'au moins un échantillon comprenant au moins une espèce à caractériser ainsi qu'une espèce de référence prise parmi l'eau ou un solvant, la teneur de l'espèce de référence dans le solvant étant au moins supérieure à deux fois la teneur de l'au moins une espèce à caractériser. Le procédé comprend les étape suivantes : a. application d'au moins un champ Bo constant sur l'au moins un échantillon, b. acquisition par une ou plusieurs antennes d'un ou plusieurs signaux complexes de décroissance libre d'induction (FID), chaque signal complexe de FID comprenant une partie réelle et une partie imaginaire, l'étape d'acquisition étant effectuée de manière à ce que dans chaque signal complexe de FID l'amplitude du signal de l'espèce de référence soit au moins deux fois supérieur à l'amplitude de celui de l'au moins une espèce à caractériser, c. pour chaque signal complexe de FID, génération d'un module de FID obtenu par extraction du module de chaque signal complexe de FID, d. application d'une transformée de Fourier à un signal prenant en compte le module de FID, cette étape étant optionnelle. De manière facultative mais avantageuse, le procédé selon l'invention peut en outre présenter au moins l'une quelconque des caractéristiques optionnelles suivantes : Selon un mode de réalisation particulier, l'échantillon est un échantillon de matière biologique, l'espèce majoritaire est de l'eau et les espèces à caractériser sont des métabolites. Selon un autre mode de réalisation particulier, l'échantillon est une composition chimique, l'espèce majoritaire est de l'eau ou un solvant et les espèces à caractériser sont des composés chimiques. Le module de FID d'un échantillon comprenant l'espèce à caractériser est égal à IIS(t)II qui est défini par l'équation suivante: dans laquelle : Aw = co-wH20 et Acp = cp - WE-120 correspondent respectivement aux décalages de fréquence et de phase entre l'au moins une espèce à caractériser et l'espèce de référence, AH2O(t) est l'amplitude en fonction du temps du signal de FID de l'espèce de référence, Ao(t) est l'amplitude en fonction du temps du signal de FID de l'espèce à caractériser Selon un mode de réalisation avantageux, au cours de l'étape d'acquisition on acquiert une pluralité de signaux complexes de décroissance libre d'induction (FID). Avantageusement, après l'étape de génération d'un module de FID pour chaque signal complexe de décroissance libre d'induction (FID) de la pluralité de signaux complexes de décroissance libre d'induction (FID), on effectue une sommation des modules de FID pour obtenir un signal FID combiné. Avantageusement, on applique ladite transformée de Fourier audit signal FID combiné. De préférence, l'étape d'acquisition est effectuée de manière à ce que dans chaque représentation spectrale du signal complexe de FID l'amplitude du signal de l'espèce de référence soit de préférence au moins deux fois supérieur à l'amplitude de celui de l'au moins une espèce à caractériser. Selon un mode de réalisation avantageux, les antennes sont des bobines et les signaux complexes de FID acquis proviennent de la même bobine. Dans ce mode de réalisation, l'invention permet avantageusement d'apporter une correction efficace et simple au déphasage et au décalage en fréquence dus notamment aux mouvements de l'échantillon pendant l'acquisition. Selon un mode de réalisation avantageux, les antennes sont des bobines et les signaux complexes de FID acquis proviennent de bobines différentes. Dans ce mode de réalisation, l'invention permet avantageusement d'apporter une correction efficace et simple au déphasage de chaque bobine les unes par rapport aux autres, ainsi qu'au déphasage et au décalage en fréquence dus notamment à l'inhomogénéité du champ BO. Selon un mode de réalisation avantageux, les antennes sont des bobines et les signaux complexes de FID acquis proviennent de bobines différentes et plusieurs signaux complexes de FID sont acquis pour au moins certaines bobines. Dans ce mode de réalisation, l'invention permet avantageusement de corriger de manière efficace et simple le déphasage et le décalage en fréquence dus notamment aux mouvements de l'échantillon et à l'inhomogénéité de BO Avantageusement, après génération d'un module de FID pour chaque signal complexe de FID, et avant sommation des modules de FID pour obtenir un signal FID combiné, on effectue une étape de calcul d'un facteur de pondération pour chaque bobine et on pondère chaque module de FID par le facteur de pondération de la bobine dont il provient. Avantageusement, l'étape de calcul d'un facteur de pondération est effectuée par la méthode de la somme des carrés des amplitudes au début du module de FID. Selon un mode de réalisation, au cours de l'étape d'acquisition on acquiert un unique signal complexe de décroissance libre d'induction (FID) et on applique la transformée de Fourier à un module de FID obtenu par extraction du module de cet unique signal complexe de FID. Selon un mode de réalisation, au cours de l'étape d'acquisition on acquiert un signal FID encodé spatialement, on répète cette acquisition plusieurs fois, de préférence au moins deux fois, afin d'obtenir plusieurs signaux encodés, chacun desdits signaux encodés étant ensuite décodé de manière à ce qu'il soit associé à un voxel de l'échantillon et on génère un module de FID pour le signal de FD associé à chaque voxel. De préférence, au cours de l'étape d'acquisition on acquiert plusieurs signaux complexes de décroissance libre d'induction (FID). Ces signaux complexes de FID comprennent un encodage spatial fonction de la position d'une pluralité de voxels de l'échantillon, dans lequel on décode ces signaux complexes de manière à obtenir un signal de FID associé à chaque voxel et dans lequel on génère un module de FID pour le signal de FID associé à chaque voxel. Selon un mode de réalisation, l'échantillon est un échantillon de matière biologique, l'espèce de référence est de l'eau et les espèces à caractériser sont des métabolites. Selon un mode de réalisation alternatif, l'échantillon est une composition chimique, l'espèce de référence est de l'eau ou un solvant et les espèces à caractériser sont des composés chimiques. Selon un mode de réalisation, le champ BO constant est appliqué sur plusieurs voxels d'un échantillon et les spectres des fréquences des signaux FID combinés des différents voxels sont utilisés pour générer une ou plusieurs images spectroscopiques. De préférence, la teneur de l'espèce de référence dans le solvant est au moins supérieure à 103 fois et de préférence 105 fois la teneur de chacune des espèces à caractériser. Dans la présente invention, la teneur de l'au moins une espèce à caractériser est la teneur massique dans l'échantillon. Préalablement à l'acquisition des signaux de FID, on effectue une étape conventionnelle de réglages nécessaires à la réalisation d'une expérience de spectroscopie de RMN. Selon un autre aspect, l'invention porte également sur un produit programme d'ordinateur comprenant des instructions, qui lorsqu'elles sont exécutées par au moins un processeur, effectuent les étapes suivantes : - réception d'un ou plusieurs signaux complexes de décroissance libre d'induction (FID), chaque signal complexe de FID comprenant une partie réelle et une partie imaginaire, tel que dans chaque signal complexe de FID l'amplitude du signal de l'espèce de référence est au moins deux fois supérieur à l'amplitude de celui de l'au moins une espèce à caractériser, - pour chaque signal complexe de FID, génération d'un module de FID obtenu par extraction du module de chaque signal complexe de FID, - application d'une transformée de Fourier à un signal prenant en compte le module de FID. Système de spectroscopie par résonance magnétique nucléaire (RMN) d'au moins un échantillon comprenant au moins une espèce à caractériser ainsi qu'une espèce de référence prise parmi l'eau ou un solvant, la teneur de l'espèce de référence dans le solvant étant au moins supérieure à deux fois la teneur de l'au moins une espèce, le système comprenant au moins une antenne configurée pour acquérir un ou plusieurs signaux complexes de décroissance libre d'induction (FID) dans le domaine temporelle, chaque signal complexe de FID étant généré par application d'au moins un champ BO sur l'au moins un échantillon et comprenant une partie réelle et une partie imaginaire. Le système comprend des moyens de traitement configurés de manière à, pour chaque signal complexe de FID, calculer le module de chaque signal complexe de FID pour générer un module de FID. De manière avantageuse mais optionnelle, les moyens de traitement sont configurés pour ensuite appliquer une transformée de Fourier à un signal prenant en compte le module de FID. De manière préférée, l'étape d'acquisition est effectuée de manière à ce que le signal de l'espèce majoritaire soit au moins deux fois supérieur à celui des espèces à caractériser De manière préférée, le système comprend également des moyens d'application du champ BO constant sur l'échantillon et des moyens d'application d'un champ générant une excitation électromagnétique dans le champ BO constant. De manière préférée, le système est configuré pour appliquer une transformée de Fourier à un signal prenant en compte le module de FID. Le système comprend un ou plusieurs antennes configurées pour acquérir une pluralité de signaux complexes de décroissance libre d'induction (FID). Il est configuré de manière à, après l'étape de génération d'un module de FID pour chaque signal complexe de décroissance libre d'induction (FID) de la pluralité de signaux complexes de décroissance libre d'induction (FID), effectuer une sommation des modules de FID pour obtenir un signal FID combiné et appliquer ladite transformée de Fourier audit signal FID combiné. Dans la présente invention, le terme antenne désigne tout type de récepteur d'ondes électromagnétiques. BRÈVE DESCRIPTION DES FIGURES Les buts, objets, ainsi que les caractéristiques et avantages de l'invention ressortiront mieux de la description détaillée d'un mode de réalisation de cette dernière qui est illustré par les dessins d'accompagnement suivants dans lesquels : La FIGURE 1 décrit brièvement le type de signaux FID détectés par toute antenne des équipements RMN connus. La FIGURE 2 compare, à partir du cas d'une seule antenne et d'une seule 30 acquisition d'un signal FID, la méthode de l'invention et le traitement conventionnel du signal FID. La FIGURE 3 illustre le cas où l'on procède à plusieurs acquisitions successives pour obtenir un signal FID ayant un RSB suffisant pour pouvoir être exploité. La FIGURE 4 compare des résultats expérimentaux qui illustrent l'avantage d'utiliser le module pour le traitement des signaux FID. La FIGURE 5 illustre le cas où l'on procède à l'acquisition d'une image spectroscopique d'un volume de tissus incluant plusieurs voxels. La FIGURE 6 compare des résultats expérimentaux obtenus pour deux voxels situés à des positions différentes dans un cerveau humain. La FIGURE 7 illustre le cas où l'on procède à l'acquisition des signaux FID à partir d'un réseau d'antennes. La FIGURE 8 illustre les étapes du procédé de l'invention pour le traitement des signaux FID capturés à partir d'un réseau d'antennes et les compare à celles de leur traitement conventionnel. La FIGURE 9 compare des résultats expérimentaux obtenus à partir d'un réseau d'antennes. Les dessins joints sont donnés à titre d'exemples et ne sont pas limitatifs de l'invention. DESCRIPTION DÉTAILLÉE DE L'INVENTION La figure 1 décrit brièvement le type de signaux détectés par toute antenne d'un équipement RMN. Le signal peut être observé après chaque émission d'une impulsion électromagnétique B1 utilisée pour perturber les noyaux d'hydrogène ou protons à l'équilibre dans le champ magnétique constant BO. Comme on l'a vu, la fréquence RF utilisée doit être de préférence en résonance avec la fréquence de Larmor des protons pour l'intensité du champ magnétique BO développé par l'équipement RMN considéré et qui est d'autant plus élevée que le champ BO est grand. Comme mentionné elle croît de 42 MHz par tesla. Le signal RMN détecté est qualifié de FID de l'anglais « free induction decay » c'est-à-dire « décroissance libre de l'induction ». Il est représentatif du retour à l'équilibre des noyaux après que la perturbation électromagnétique à laquelle ils sont soumis a cessé. Par des moyens classiques qui ne sont pas décrits et qui ne sont pas nécessaires à la compréhension de l'invention, le signal FID capturé provient d'une entité volumique élémentaire de la zone couverte par l'antenne 101. Cette entité élémentaire est qualifiée de voxel 110. La taille physique d'un voxel dépend de la résolution volumique de l'équipement RMN utilisé. Le signal capturé par l'antenne est classiquement traité par un récepteur électronique 120 comprenant principalement un amplificateur analogique 121 et deux convertisseurs analogique numérique (CAN) 124. De manière conventionnelle, après amplification, on procède en effet à une conversion de fréquence et à une détection en quadrature du signal analogique FID capturé par l'antenne à partir de signaux 122 déphasés de 90°. La conversion de fréquence a pour but de transposer le signal FID dans une gamme de fréquences compatible avec le fonctionnement des CAN qui suivent. Après filtrage 123 on obtient ainsi, pour chaque antenne, un signal complexe ayant une partie dite « réelle » et une partie dite « imaginaire ». Les signaux électriques analogiques ainsi obtenus, tous deux bien réels, sont alors échantillonnés chacun par un CAN 124 de façon à disposer d'un signal FID complexe sous forme numérique 130 ayant une composante « réelle » 132 et une composante « imaginaire » 131, c'est-à-dire en quadrature avec la partie réelle. Typiquement, pour chaque signal FID capturé, le nombre d'échantillons numériques dont on dispose pour chacune des voies est de l'ordre de 210 c'est-à-dire 1024 échantillons. Le signal numérique 130 est ainsi susceptible d'être traité à l'aide de tous les moyens logiciels et matériels qui ont été développés depuis des décennies pour le traitement numérique des signaux analogiques. Ces moyens comprennent notamment des processeurs numériques spécialisés dits « processeurs de signaux » et d'algorithmes spécifiques notamment ceux pour mettre en oeuvre les transformées de Fourier qui permettent de faire une analyse spectrale des signaux temporels reçus. Les diagrammes 140 et 150 représentent, sous forme analogique, les deux composantes d'un exemple de signal FID. Il s'agit, respectivement, de la partie réelle et de la partie imaginaire d'un signal FID capturé par l'antenne après amplification, conversion et détection en quadrature du signal. En première approximation, le signal FID est une simple exponentielle décroissante. En réalité, comme il peut être perturbé par de nombreuses sources extérieures, le fait d'introduire une amplitude qui dépende du temps permet de ne faire aucune supposition sur la forme du signal FID. D'un point de 5 vue calcul mathématique chaque signal d'antenne est ainsi une expression complexe qui est fonction du temps (t) de la forme : A(t) = Ao(t)el(wt+9) ou w représente la fréquence de résonance de l'eau et ço la distorsion de phase. L'exemple de la figure 1 représente le cas simple d'une acquisition unique du signal FID à partir d'une seule antenne. On suppose en effet dans 10 cet exemple simple que le signal capturé est d'amplitude suffisante pour être traité après une seule acquisition et qu'il n'est pas nécessaire de procéder à plusieurs acquisitions successives pour améliorer le RSB comme on le verra dans certains des cas décrits ci-après. On pourra noter que si le patient bouge pendant l'acquisition d'une FID 15 alors on a une modification de la fréquence et de la phase entre le début de la FID et la fin de la FID ce qui peut entrainer une modification de la forme des raies après transformée de Fourrier et pourrait diminuer la précision des résultats obtenus avec le procédé selon l'invention. Ceci est en pratique peu probable car le temps d'acquisition d'une FID est de quelques millisecondes. 20 Par contre si le patient bouge entre le moment où l'on a réglé la fréquence est le moment où l'on commence l'acquisition de la FID (de 1 à plusieurs secondes), on se retrouve dans le cas décrit par la suite. On remarquera aussi que le signal de l'eau est normalement « en résonance » c'est-à-dire que l'on règle la fréquence de conversion du récepteur 25 exactement sur la fréquence de résonance du signal de l'eau. Lorsque ce réglage est exact, le signal de l'eau apparait effectivement comme une exponentielle décroissante. Lorsque ce réglage n'est pas exact, la fréquence du signal de l'eau est légèrement différente de la fréquence de conversion. Cela se traduit par une oscillation du signal FID. L'exponentielle décroissante est alors 30 modulée par un sinus dont la fréquence est égale à la différence entre la fréquence du signal de l'eau et la fréquence de conversion du récepteur. L'oscillation du signal FID de la figure 1 traduit ce décalage en fréquence. Un signal acquis exactement « en résonance » serait semblable à celui du diagramme 222 de la figure 2. En ce qui concerne la phase du signal, elle est arbitraire car la phase du signal FID reçu est complètement indépendante de la phase du signal de conversion du récepteur. La conversion en fréquence au niveau du récepteur se fait en soustrayant du signal FID un signal de fréquence égale à la fréquence de Larmor. Mais le récepteur ne connait pas la phase du signal de FID. Il y a donc forcément un décalage aléatoire entre la phase du signal utilisé par le récepteur et la phase du signal de l'eau. Cela se traduit sur la figure par le fait que le premier point de la FID n'est en général pas réel. Les coordonnées de ce premier point dans le plan complexe sont dans cet exemple, si l'on se reporte aux diagrammes 140 et 150, approximativement (-60,-100). Si le signal était bien en phase les coordonnées de ce point serait (116,0) c'est à dire 116 pour la partie réelle, et 0 pour la partie imaginaire. La figure 2 compare, à partir du cas de la figure 1, c'est-à-dire celui d'une seule antenne et d'une seule acquisition du signal, la méthode de l'invention et le traitement conventionnel du signal FID. Comme on l'a vu dans la figure précédente on dispose pour le signal capturé par une antenne d'une partie réelle 140 et d'une partie imaginaire 150. Le traitement conventionnel du signal FID est donc effectué dans le plan complexe. La transformée de Fourier correspondante 210, qui permet de passer du domaine temporel au domaine fréquentiel et d'obtenir une analyse spectrale du signal reçu, reflète les glissements de fréquence 213 et de phase possiblement présents dans le signal complexe 130 obtenant ainsi le spectre fréquentiel 212 auquel il faut appliquer successivement une correction de phase 214 et de fréquence 216 pour retrouver le signal de l'eau 211 en résonance, c'est-à-dire calé sur la fréquence 0 de la transformée de Fourier Comme mentionné dans la description de la figure 1, le signal complexe de RMN est obtenu en déphasant de 90° une partie du signal reçu. Cependant, le signal recueilli est bien réel, et une fois que toutes les opérations ont été menées dans le plan complexe, pour plus de facilité mathématique, on s'arrange pour que le spectre obtenu soit lui aussi réel. À cette fin, la partie imaginaire introduite pour des raisons de facilité de calcul est alors retirée. Si l'on suppose que le signal FID est une exponentielle décroissante on obtient alors des raies étroites dites « Lorentziennes » telles que celles représentées dans le diagramme 214. Les raies du diagramme 212 sont un mélange de lorentziennes et de raies en dispersion. L'opération de mise en phase ou phasage consiste à récupérer les lorentziennes dans la partie réelle et à mettre les raies en dispersion dans la partie imaginaire. Une fois cette opération réalisée, la partie imaginaire peut être mise de coté et on ne continue à travailler que sur la partie réelle. Le glissement de fréquence est le même que celui précédemment décrit en figure 1. Dans la figure 1 il se traduit par une modulation du signal FID alors que dans la figure 2 c'est sa représentation spectrale qui apparaît. Le glissement en fréquence est dans ce cas, comme on l'a vu, dû à la différence entre la fréquence de démodulation utilisée par le récepteur et la fréquence effective du signal de l'eau. On remarquera ici qu'il faut pouvoir différencier ce qui est phase et fréquence dites « globales », telles que décrites ci-dessus, des phase et fréquence dites « relatives » qui correspondent à une différence de phase ou de fréquence par rapport à un autre signal FID, ou à une autre localisation. Si on acquiert successivement deux signaux FID, ils auront tous deux la même phase et la même fréquence globale puisque le signal de démodulation du récepteur n'a pas changé d'un signal FID à l'autre. Par contre si le patient a bougé pendant l'acquisition, la phase et la fréquence du second signal FID seront légèrement différentes du premier. C'est le cas également pour des signaux FID provenant de voxels différents comme dans le cas d'une analyse dite CSI décrite en figure 5 (une analyse CSI, acronyme de l'anglais « chemical shift imaging » signifie « imagerie de déplacement chimique » de la fréquence de résonance des noyaux comme déjà discuté brièvement plus haut). La phase et la fréquence provenant de voxels différents seront légèrement différentes en raison des variations de champ BO, mais fréquence et phase « globales » seront les mêmes pour tous les voxels. Il y a donc en fait deux corrections à faire : une générale qui consiste à corriger les variations de phase et de fréquence générées par le convertisseur, puis FID par FID ou voxel par voxel, on pratique une correction des petites variations induites par le mouvement du patient ou les variations de champ BO. L'invention permet de réaliser ces deux opérations en une seule fois. Le traitement selon l'invention du signal capturé par l'antenne commence par le calcul du module 220 du signal complexe 130 tel que défini ci-dessus. Dans l'expression complexe mentionné ci-dessus et qui est fonction du temps (t) A(t) = Ao(t)eeut+9), le terme Ao(t) est le module ou amplitude du signal complexe. Ce module est indépendant de la fréquence et de la phase. D'un point de vue pratique, le module à chaque instant t est donc la racine carrée de la somme des carrés des parties réelle et imaginaire du signal FID à cet instant t tel que capturé et converti par le module électronique 120 associé à chaque antenne. Ainsi, à chaque instant t, le module du signal FID est défini par Module de FID(t) =I A(t) I = (Ao(t)2 [cos(wt-F(p)2+sin(wt-F(p)2D1/2 Le terme [cos(wt-F(p)2+sin(wt-F(p)2] étant égal à 1, on obtient le module du signal FID(t) = Ao(t). S'étant affranchi des termes de phase et de fréquence parasites, on peut maintenant se consacrer à la quantification de Ao(t) ce qui est en réalité le but de l'expérience. Le signal obtenu correspond donc à chaque instant au module de la FID. Il est qualifié dans la présente de demande de brevet de « module de FID » ou plus simplement de module. Un exemple de ce type de signal est illustré par le diagramme 222. On applique ensuite sur le module de FID 222 la transformée de Fourier 230 afin d'obtenir son spectre fréquentiel 232. Le module du signal capturé par l'antenne ne portant aucune information de fréquence ni de phase ne nécessite donc pas qu'on y applique les corrections 214 et 216 du mode de traitement conventionnel brièvement décrit ci-dessus. En particulier le spectre de l'eau 231 est alors toujours bien centré sur la fréquence 0. Le spectre du module, qui est un signal ayant une partie réelle uniquement, est caractérisé par le fait qu'il est symétrique 233. Il est directement exploitable par les processus suivants qui sont particuliers à chaque application RMN. L'invention permet ainsi de se passer d'étapes dédiées aux corrections de phase et de fréquence, ces étapes nécessitant habituellement un temps de traitement significatif, voire bien souvent l'intervention manuelle d'un opérateur. De préférence, on choisira un échantillon qui ne présente pas, avec une solution classique basée sur une transformée de Fourier dans le domaine temporel, de signal de part et d'autre de la fréquence correspondant à celle de l'espèce majoritaire. Le spectre du module étant symétrique on n'en retient en général qu'une des moitiés. On remarquera aussi que la transformée de Fourier du module de la FID est plus simple notamment parce que la quantité de données à traiter est moitié de celles qu'il faut traiter de façon conventionnelle où l'on opère simultanément sur les parties réelles et sur les parties imaginaires des signaux FID capturés. Le temps de traitement est donc considérablement réduit ce qui permet d'effectuer un traitement en temps réel sans impacter les utilisateurs. En particulier, on notera que le fait de ne travailler que sur un signal réel divise par deux l'occupation de la mémoire vive de l'ordinateur. Ceci permet de travailler sur un plus gros volume de donnée sans risquer un ralentissement de la machine. De plus, le fait de n'avoir à garder que la moitié du spectre après traitement divise par quatre le volume total de données à sauvegarder. Le gain en temps de calcul n'est effectivement obtenu que si la taille des données est supérieure à la taille de la mémoire allouée par l'ordinateur. Ceci est cependant souvent le cas notamment pour les analyses dites CSI, plus particulièrement décrites en figure 5, qui se font en haute résolution et mettent en jeu des tailles de données pouvant atteindre quatre gigaoctets, c'est-à-dire 4x109 octets. De manière surprenante, bien qu'une partie des informations captées par les antennes (l'information de phase) soit éliminée lors de la prise en compte du module de la FID, le RSB demeure parfaitement satisfaisant. On notera ici que l'invention s'adresse plus particulièrement aux applications RMN qui s'intéressent à l'analyse des pics secondaires 234 qui apparaissent dans la décomposition spectrale et qui sont engendrés par des métabolites possiblement présents dans les tissus examinés. L'eau constituant la majeure partie des tissus génère un signal de très forte amplitude 231 qui agit comme une porteuse et comme une référence pour les pics secondaires 234 que l'on désire détecter et qui sont caractéristiques dans le voxel examiné de la présence de métabolites. Ceux-ci créent un « chemical shift », terme utilisé dans la littérature technique en anglais sur ces sujets pour qualifier un « glissement ou décalage chimique » de la fréquence de précession des noyaux mis en résonance. Le glissement ou décalage fréquentiel apparaît dans la décomposition spectrale sous la forme de pics secondaires décalés qui révèlent la présence de ces métabolites. Le décalage en fréquence est caractéristique de la nature du métabolite. La surface du pic, relativement à celle de l'eau, permet de quantifier leur concentration dans le tissu examiné. La figure 3 illustre le cas où il faut procéder à plusieurs acquisitions 310 successives pour obtenir un signal FID ayant un RSB suffisant pour pouvoir être exploité par l'application. Les traitements conventionnel et selon l'invention du signal FID sont pour l'essentiel identiques à ce qui est décrit dans la figure 2. Les seules différences concernent les étapes de sommation. Dans le traitement conventionnel on doit sommer séparément les parties réelles 320 et imaginaires 330 des signaux FID successifs capturés Au contraire, dans le traitement selon l'invention on procède d'abord au calcul du module 340 de chaque signal de FID. Seuls les modules de FID ont alors besoin d'être sommés 350 avant d'appliquer la transformée de Fourier 230 sur le module résultant, ce qui simplifie sensiblement le processus de calcul dans le cas de l'invention. La figure 4 compare des résultats expérimentaux qui illustrent l'avantage d'utiliser le module 420 pour le traitement des signaux FID au lieu d'utiliser le traitement conventionnel séparé des parties réelles et imaginaires 410. Les figures 410 et 420 mettent l'accent sur la partie intéressante de chacun des spectres c'est-à-dire qu'elles ne montrent que les parties concernant les métabolites. Le signal de l'eau et le signal symétrique du volume ont été éliminés afin de rendre plus visible la région 422 du spectre. Le diagramme 420 montre que le traitement du module améliore sensiblement le RSB du signal FID capturé et permet de diminuer la largeur des raies du spectre. Notamment, comme le montre cet exemple, la diminution de la largeur des raies révèle quelques petites résonances secondaires 422 qui n'étaient pas visibles dans le spectre conventionnel. L'invention permet ainsi de détecter des espèces difficilement détectables voire non détectables avec les solutions de l'art antérieur. La forme des résonances est également améliorée 424. On sait que la correction de phase et de décalage ou glissement en fréquence de chaque FID, illustrée dans la figure 2, avant sommation, améliorerait la largeur des raies du spectre et le RSB. Cette correction n'est cependant en pratique jamais effectuée en raison de la difficulté qu'il y aurait à trouver la phase et la fréquence de chaque FID en temps réel avant sommation. Utiliser le module des signaux FID comme dans l'invention est donc une manière rapide, simple et robuste pour corriger l'effet des mouvements de la cible qui sont principalement à l'origine des glissements de phase et de fréquence observés lors de l'acquisition des signaux FID. La figure 5 illustre le cas où l'on procède à l'acquisition d'une image spectroscopique d'un volume de tissus incluant plusieurs voxels. Ce type d'analyse par CSI se pratique de façon classique en surimposant des gradients de magnétisation à la magnétisation fixe Bo de façon à pouvoir pratiquer une discrimination spatiale des signaux en provenance de chacun des voxels. Ce cas est similaire à celui de la figure 1 mais s'adresse à un volume de voxels 510 pour lesquels on peut constater entre eux des décalages de phase et de fréquence dus aux inhomogénéités du champ magnétique BO. On remarquera ici que tout ce qui est fait dans le domaine spatial est strictement identique à ce qui se fait de manière conventionnel en imagerie RMN. De la même façon que l'on réalise une transformation de Fourier pour passer d'un signal FID à son spectre, on réalise une transformation de Fourier pour passer des signaux codés par les gradients de champ magnétique aux signaux spatialement résolus. La transformation de Fourier qui convertit le signal FID en un spectre travaille sur un signal temporel et le transforme en fréquentiel d'ou le nom de TF temporelle (ou fréquentielle). La transformation de Fourier qui fait passer des signaux codés par les gradients en une image RMN travaille dans la dimension spatiale. Ce traitement est spécifique à une analyse CSI et est effectué en amont, avant même que l'invention ne soit appliquée. La figure 5 montre que la technique de l'invention est compatible avec celle des analyses CSI. Comme on l'a vu précédemment les corrections globales sont celles qui concernent tous les signaux FID de manière identique. Il s'agit principalement des décalages en fréquence et en phase entre le signal de démodulation du récepteur et le signal effectivement acquis. Les corrections locales ou relatives sont celles qui concernent les différences de phase et de fréquence observées entre deux signaux FID acquis, soit à partir de deux bobines différentes, soit successivement par une même bobine ou soit encore si les signaux FID proviennent de régions spatialement différentes. Les différences observées sont, suivant les cas, provoquées par : un mouvement du patient, l'électronique des bobines, les variations spatiales de champs BO ou une combinaison de ces phénomènes. La figure 6 compare des résultats expérimentaux obtenus dans les conditions décrites dans la figure précédente pour deux voxels situés à des positions différentes dans un cerveau humain. Les diagrammes 610 et 620 d'une part et 630 et 640 d'autre part décrivent, respectivement, les signaux obtenus pour deux voxels à des positions différentes (1 et 2) avec un traitement conventionnel des signaux FID capturés et un traitement selon l'invention de leurs modules. On remarque que les décalages de fréquence 650 et de phase 660 qui apparaissent dans le traitement conventionnel des signaux capturés n'apparaissent plus 670 dans le traitement du module ce qui améliore très sensiblement leur quantification. On notera ici que dans le cas décrit dans les figures 5 et 6 ci-dessus, c'est-à-dire pour l'acquisition de signaux CSI, on pourra avoir aussi à procéder à des acquisitions multiples pour améliorer le RSB. L'invention s'appliquera également comme décrit précédemment et en particulier comme décrit dans la figure 3 avec sommation des modules des FID. La figure 7 illustre le cas où l'on procède à l'acquisition des signaux FID à partir d'un réseau d'antennes comme discuté dans le chapitre sur l'état de la technique. L'invention est particulièrement avantageuse pour ce cas complexe où il faut aussi pouvoir corriger facilement les différences importantes, notamment de phase, qui peuvent exister entre les antennes individuelles. Dans ce cas, le signal FID émis par chaque voxel 110 est capturé par plusieurs antennes. Deux antennes, 101 et 102, sont représentées dans cet exemple. Il est courant d'en associer un plus grand nombre, par exemple des réseaux de huit antennes sont fréquemment utilisés. Chaque antenne individuelle est munie de son propre système de réception, 1201 et 1202 dans cet exemple, délivrant chacun un signal FID complexe, 1301 et 1302. Les diagrammes 710 et 730 d'une part et 720 et 740 d'autre part sont, respectivement, des exemples des parties réelles et des parties imaginaires des signaux FID complexes délivrés par chacune des antennes. Comme attendu, on peut constater de fortes différences 750 dans la phase et l'amplitude des signaux capturés par les antennes individuelles du réseau d'antennes. On notera ici qu'il existe un décalage global de phase et de fréquence qui est commun à toutes les antennes individuelles du réseau cependant les différences individuelles de phase et d'amplitude entre antennes doivent toujours être compensées comme on va le voir ci-après. La figure 8 illustre les étapes du procédé de l'invention pour le traitement des signaux FID capturés à partir d'un réseau d'antennes et les compare à celles de leur traitement conventionnel. Le processus de l'invention reste extrêmement simple et n'est pas sensiblement différent de ce qui a été décrit dans la figure 3 pour l'acquisition successive de signaux FID à partir d'une seule antenne. Dans le cas de la figure 8 les signaux sont acquis, de préférence simultanément, et on procède au calcul du module 340 sur chacun des signaux acquis à partir des antennes individuelles 101 et 102 dans cet exemple. Il n'est pas exclu bien sûr de procéder aussi, comme dans la figure 3, à des acquisitions multiples dans le cas d'un réseau d'antennes. Une différence notable dans le cas d'un réseau d'antenne est qu'il faut à l'étape suivante 810 procéder à une correction de l'amplitude des modules provenant d'antennes différentes afin de tenir compte de toutes les différences existant inévitablement entre les antennes individuelles formant le réseau. Une méthode pour appliquer cette correction de l'amplitude consiste à 30 effectuer une pondération de chaque signal de FID. On peut par exemple appliquer à chaque module de FID un facteur de pondération de la bobine dont il provient, ce facteur de pondération étant égal au carré du module de la FID à l'instant t=0 divisé par la somme des carrés des module des chaque FID à l'instant t=0. On procède ensuite, comme précédemment, à la sommation 350 des modules dont l'amplitude vient d'être pondérée. On peut alors calculer la 5 composition spectrale du module résultant de la combinaison des signaux FID à l'aide d'une transformée de Fourier 230. Il convient de noter que la pondération des signaux de FID, bien que préférable, n'est qu'optionnelle. Par comparaison le processus conventionnel, beaucoup plus complexe, 10 figure sur la partie gauche de la figure 8. Cette partie est aussi à comparer avec la figure 2 qui ne concerne qu'un signal FID. On remarquera ici que dans le cas de la figure 2 la correction de phase peut être indifféremment réalisée avant ou après la transformation de Fourier puisqu'il n'y a qu'un seul signal FID. Elle est traditionnellement effectuée manuellement après la transformée de Fourier car 15 il est plus simple de visualiser les variations de phase sur le spectre que sur le signal FID. Dans le cas de la figure 8 cette correction est réalisée avant car on recombine d'abord les signaux des différentes bobines avant d'effectuer le traitement. Toutefois, on pourrait aussi envisager de réaliser la transformée de Fourier de chaque signal de chaque bobine, puis de les phaser et de les 20 sommer. Cependant, cela demande alors de réaliser autant de transformations de Fourier qu'il y a de bobines, par exemple 64 si l'on dispose de 64 bobines. La figure 9 compare des résultats expérimentaux obtenus à partir d'un réseau d'antennes dans les conditions décrites dans la figure précédente. Les diagrammes 910 et 920 montrent les spectres obtenus après 25 sommation des signaux FID, transformation de Fourrier et suppression de la raie de l'eau dans ces exemples pour mieux visualiser les métabolites. Le diagramme 920 montre la partie réelle du spectre obtenue à l'issue du traitement conventionnel des signaux FID. Comme mentionné précédemment, la partie imaginaire est introduite artificiellement par le récepteur et ne contient 30 aucune information supplémentaire. Elle n'est donc utilisée que lors du traitement des signaux, jamais dans la présentation des résultats. Le diagramme 910 montre le résultat du traitement selon l'invention du module des signaux FID. Comme déjà mentionné, le traitement du module conduit à une symétrisation du spectre dont on ne garde qu'une moitié. Les pics présents sur le diagramme 920 obtenu avec les systèmes connus se retrouvent bien dans le diagramme 910 obtenu avec la solution de l'invention. Ainsi, en dépit de la perte d'information induite par la prise en compte du module avant sommation des différentes FID, les métabolites peuvent toujours être caractérisés. On remarque en outre que ses pics sont très étroits, ce qui permet une caractérisation précise et montre que le rapport signal/bruit n'est pas dégradé de manière visible par le traitement selon l'invention du module du signal. Comme le montre la figure 9 il peut être avantageux de supprimer l'eau pour améliorer l'acquisition des signaux provenant des métabolites. En effet, dans le cerveau, le signal de l'eau est très supérieur à celui des métabolites, de l'ordre de 105 fois. Afin de bénéficier de toute la dynamique des convertisseurs analogique numérique (CAN) le signal de l'eau est souvent supprimé partiellement ou complètement avant acquisition du signal FID Le signal de l'eau est typiquement supprimé en deux étapes. Une première étape est effectuée lors de l'acquisition en utilisant une séquence appropriée (il s'agira d'une atténuation dans ce cas). Comme cette suppression n'est jamais parfaite, une seconde suppression est réalisée après sommation pour diminuer le nombre de signaux à traiter, en général juste avant la quantification des spectres. Comme mentionné ci-dessus la première suppression a pour but de diminuer l'amplitude du signal reçu par le convertisseur analogique numérique du récepteur, et donc de garder la pleine plage de codage, qui est typiquement de 16 bits, pour numériser au mieux les signaux les plus petits. La seconde étape est réalisée pour faciliter la quantification des métabolites. Pour cela on modélise le signal de l'eau en réalisant une décomposition en valeur singulière du signal, et en ne retenant dans les composantes principales que celles qui correspondent au signal de l'eau. On reconstruit un signal à partir de ces valeurs que l'on soustrait au spectre. Toutefois, comme on l'a vu, le signal de l'eau peut être aussi utilisé comme signal de référence de phase et de fréquence. On pratique alors une suppression partielle de l'eau pour garder un signal d'amplitude suffisante pour être exploité. Comme déjà aussi mentionné précédemment le signal de l'eau peut également être utilisé comme référence de concentration pour estimer la concentration absolue des métabolites. Il est dans ce cas nécessaire d'obtenir la totalité du signal de l'eau. Souvent, pour concilier les impératifs liés à une meilleure visualisation des métabolites d'une part et à leur quantification d'autre part on peut être amené à pratiquer deux expérimentations successives, l'une avec suppression partielle ou totale de l'eau et l'autre sans. Les aspects théoriques de l'invention sont supportés par les développements mathématiques qui suivent. Si A(t) est le signal acquis depuis un voxel fait d'un seul composé, c'est- à-dire une seule espèce, en l'absence de bruit, il s'exprime comme indiqué dans l'équation [1] ci-après : dans laquelle Ao(t) est la forme d'onde du signal FID, w la fréquence de résonance du composé unique et cp représente toutes les distorsions de phase. Une façon de corriger les distorsions de phase consiste à prendre le module du signal comme indiqué dans l'équation [2] ci-après : [2] Le module enlève également l'information de fréquence et après transformation de Fourier le résultat est un signal « en résonance ». Si on suppose maintenant que le composé unique correspond au pic de l'eau, on obtient, en utilisant la même notation que dans l'équation [1], l'expression suivante : On notera que dans ces équations l'espèce de référence est considérée comme étant l'eau. L'invention n'est cependant pas limitée à l'eau comme espèce de référence. Ainsi, ces équations seront pourront être utilisées en remplaçant l'eau par un solvant faisant office d'espèce de référence. En ajoutant un autre composé au signal de l'eau, le module du signal acquis devient alors : Il I 1-1201) + ) (H2ow+ (f). [4] Dans cette expression le caractère * indique la conjuguée complexe. En utilisant les équations [1] et [3], après développement, le produit sous la racine carrée peut être écrit comme suit : Si on introduit maintenant les termes Aco = co-coH20 et Aço = ço - 9H20 qui 5 correspondent respectivement aux glissements de fréquence et de phase entre le composé ajouté et celui de l'eau on obtient : icos(Aox - A i) + C I ) cos( Awf j{ + 1 I ) C'(:)';( A(0)! (JO j1 + I 2(1- COH A(r)( - = ) +Acof + -',( }sin( A(..i1)) et finalement l'équation [5] II.s 01= (0+ oc.i)cos Accx--t n sin( - A0)) (/ )sicciC Aue + Arp) 2 [5] + )cos Aox cf))) 11+ (^+.',Ît)cos(i91 + .41) - - En faisant l'hypothèse que AH2O(t) est très grand devant A2(t), ce qui est le cas dans toutes les expérimentations où l'on ne supprime pas l'eau, le 10 second terme de l'expression peut être négligé. On obtient alors l'expression approchée suivante : } cos (01 - .A(;!)1 [6] Bien que la description détaillée se réfère à titre d'exemple à la caractérisation de métabolites d'un échantillon biologique, la présente invention 15 s'avère également avantageuse pour l'analyse d'une composition chimique non biologique. En effet, l'invention permet par exemple d'identifier des composés chimiques dans une solution comprenant une espèce de référence telle que de l'eau ou un solvant et de mesurer la teneur de ces composés dans cette solution. L'invention n'est pas limitée aux modes de réalisation décrits mais s'étend à tout mode de réalisation couvert par les revendications. REVENDICATIONS1. Procédé de correction de phase et/ou de fréquence d'au moins un signal de décroissance libre d'induction (FID) se rapportant à au moins un échantillon à analyser par résonance magnétique nucléaire (RMN), l'échantillon comprenant au moins une espèce à caractériser ainsi qu'une espèce de référence prise parmi l'eau ou un solvant, la teneur de l'espèce de référence dans le solvant étant au moins supérieure à deux fois la teneur de l'au moins une espèce à caractériser, le procédé comprenant les étape suivantes : a. application d'au moins un champ Bo constant sur l'au moins un échantillon, b. acquisition (310) par une ou plusieurs antennes d'un ou plusieurs signaux complexes de décroissance libre d'induction (FID), chaque signal complexe de FID comprenant une partie réelle (140) et une partie imaginaire (150), l'étape d'acquisition étant effectuée de manière à ce que dans chaque signal complexe de FID l'amplitude du signal de l'espèce de référence soit au moins deux fois supérieur à l'amplitude de celui de l'au moins une espèce à caractériser, caractérisé en ce que le procédé comprend également au moins l'étape suivante : c. pour chaque signal complexe de FID, génération (340) d'un module (222) de FID obtenu par extraction du module de chaque signal complexe de FID, la phase de chaque signal complexe de FID étant ainsi remise à 0. 2. Procédé selon la revendication précédente, comprenant une étape (230) ultérieure d'application d'une transformée de Fourier à un signal prenant en compte ledit module (222) de FID. 3. Procédé selon l'une quelconque des revendications précédentes, dans lequel le module (222) de FID d'un échantillon comprenant l'espèce à caractériser est égal à IIS(t)II qui est défini par l'équation suivante: 31 ) dans laquelle : Ac° = co-coH2oet Aço = cA - 9H20 correspondent respectivement aux décalages de fréquence et de phase entre l'au moins une espèce à caractériser et l'espèce 5 de référence, AH20(t) est l'amplitude en fonction du temps du signal de FID de l'espèce de référence, Ao(t) est l'amplitude en fonction du temps du signal de FID de l'espèce à caractériser. 10 4. Procédé selon l'une quelconque des revendications précédentes, dans lequel au cours de l'étape d'acquisition (310) on acquiert une pluralité de signaux complexes de décroissance libre d'induction (FID). 15 5. Procédé selon la revendication précédente, dans lequel après l'étape de génération (340) d'un module (222) de FID pour chaque signal complexe de décroissance libre d'induction (FID) de la pluralité de signaux complexes de décroissance libre d'induction (FID), on effectue une sommation (350) des modules de FID pour obtenir un signal FID combiné. 20 6. Procédé selon la revendication précédente, dans lequel on applique ladite transformée de Fourier (230) audit signal FID combiné. 7. Procédé selon l'une quelconque des revendications précédentes, 25 dans lequel les antennes sont des bobines et dans lequel les signaux complexes de FID acquis proviennent de la même bobine. 8. Procédé selon l'une quelconque des revendications 1 à 6, dans lequel les antennes sont des bobines et dans lequel les signaux complexes 30 de FID acquis proviennent de bobines différentes. 9. Procédé selon l'une quelconque des revendications 1 à 6, dans lequel les antennes sont des bobines et dans lequel les signaux complexes de FID acquis proviennent de bobines différentes et plusieurs signaux complexes de FID sont acquis pour au moins certaines bobines. 10. Procédé selon l'une quelconque des deux revendications précédentes, dans lequel après génération d'un module (222) de FID pour chaque signal complexe de FID, et avant sommation (350) des modules de FID pour obtenir un signal FID combiné, on effectue une étape de calcul d'un facteur de pondération (810) pour chaque bobine et on pondère chaque module (222) de FID par le facteur de pondération de la bobine dont il provient. 11. Procédé selon la revendication précédente, dans lequel l'étape de calcul d'un facteur de pondération est effectuée par la méthode de la somme des carrés des amplitudes au début du module (222) de FID. 12. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel au cours de l'étape d'acquisition on acquiert un unique signal complexe de décroissance libre d'induction (FID) et dans lequel on applique la transformée de Fourier à un module (222) de FID obtenu par extraction du module de cet unique signal complexe de FID. 13. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel au cours de l'étape d'acquisition (310) on acquiert un signal FID encodé spatialement, on répète cette acquisition plusieurs fois afin d'obtenir plusieurs signaux encodés, chacun desdits signaux encodés étant ensuite décodé de manière à ce qu'il soit associé à un voxel de l'échantillon et on génère un module (222) de FID pour le signal de FD associé à chaque voxel. 14. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'échantillon est un échantillon de matière biologique, l'espèce de référence est de l'eau et les espèces à caractériser sont des métabolites. 15. Procédé selon l'une quelconque des revendications 1 à 13, dans lequel l'échantillon est une composition chimique, l'espèce de référence est de l'eau ou un solvant et les espèces à caractériser sont des composés chimiques. 16. Procédé selon l'une quelconque des revendications précédentes, dans lequel le champ BO constant est appliqué sur plusieurs voxels d'un échantillon et dans lequel les spectres des fréquences des signaux FID combinés des différents voxels sont utilisés pour générer une ou plusieurs images spectroscopiques. 17. Procédé selon l'une quelconque des revendications précédentes, dans lequel la teneur de l'espèce de référence dans le solvant est au moins supérieure à 103 fois la teneur de chacune des espèces à caractériser. 18. Système de spectroscopie par résonance magnétique nucléaire (RMN) d'au moins un échantillon comprenant au moins une espèce à caractériser ainsi qu'une espèce de référence prise parmi l'eau ou un solvant, la teneur de l'espèce de référence dans le solvant étant au moins supérieure à deux fois la teneur de l'au moins une espèce, le système comprenant au moins une antenne configurée pour acquérir (310) un ou plusieurs signaux complexes de décroissance libre d'induction (FID) dans le domaine temporelle, chaque signal complexe de FID étant généré par application d'au moins un champ BO sur l'au moins un échantillon et comprenant une partie réelle et une partie imaginaire, caractérisé en ce que le système comprend des moyens de traitement configurés de manière à, pour chaque signal complexe de FID, calculer le module (340) de chaque signal complexe de FID pour générer un module (222) de FID. 19. Système selon la revendication précédente, comprenant également des moyens d'application du champ BO constant sur l'échantillonet des moyens d'application d'un champ générant une excitation électromagnétique dans le champ BO constant. 20. Système selon l'une quelconque des deux revendications précédentes, - configuré pour appliquer une transformée de Fourier (230) à un signal prenant en compte le module (222) de FID, - comprenant un ou plusieurs antennes configurées pour acquérir une pluralité de signaux complexes de décroissance libre d'induction (FID), - configuré de manière à, après l'étape de génération d'un module (222) de FID pour chaque signal complexe de décroissance libre d'induction (FID) de la pluralité de signaux complexes de décroissance libre d'induction (FID), effectuer une sommation (350) des modules de FID pour obtenir un signal FID combiné et appliquer ladite transformée de Fourier audit signal FID combiné.