FR 2956744 A1 20110826 FR 1051301 A 20100224 PROCEDE DE REALISATION D'UN MODULE DE PILOTAGE D'UNE ANTENNE HAUTE FREQUENCE POUR APPAREIL D'IMAGERIE PAR RESONNANCE MAGNETIQUE NUCLEAIRE. La présente invention concerne le pilotage d'une antenne haute fréquence à voies multiples en réception et/ou en transmission, utilisées dans les appareils de Résonnance Magnétique Nucléaire (RMN), ainsi que leurs applications : Imagerie par Résonance Magnétique (IRM) pour l'homme ou l'animal, Spectroscopie par Résonance Magnétique (SRM), ou encore imagerie du tenseur de diffusion. Le domaine de l'invention est celui de la Résonance Magnétique Nucléaire (RMN) et plus particulièrement celui des antennes haute-fréquence du type réseau phasé utilisées dans les appareils d'Imagerie par Résonance Magnétique (IRM) à très haut champ magnétique, c'est-à-dire supérieur ou égal à 7 Tesla. L'invention s'applique particulièrement au pilotage d'une antenne haute fréquence de type réseau phasé utilisée pour examiner une partie du corps humain comme par exemple la tête d'un patient. D'une façon connue, les antennes des appareils d'IRM sont formées par des éléments rayonnants de forme variable, en cuivre ou de bande de cuivre, jouant soit le rôle de bobinages émetteurs, en résonance avec le circuit qui l'alimente, soit le rôle de récepteurs de signaux de relaxation, soit plus généralement ces deux fonctions alternativement. Elles entourent le patient ou seulement la partie du corps à analyser. Placées dans le champ magnétique longitudinal permanent Bo, d'un dispositif de RMN, elles reçoivent une excitation électrique leur permettant de produire le champ magnétique radiofréquence (RF) B orthogonal à Bo, à la fréquence de résonance de précession (également appelée fréquence de Larmor) des noyaux des atomes étudiés qui se trouvent dans le champ magnétique Bo, et/ou captent l'écho B1, également onde radiofréquence à la fréquence de Larmor, produit par les atomes qui ont été momentanément soumis au champ magnétique B; . On rappelle que sous l'effet du champ magnétique statique Bo, les moments magnétiques de spin des noyaux d'atomes d'hydrogène vont progressivement s'aligner dans une direction initialement parallèle au champ magnétique Bo et donner lieu à une aimantation globale dans la direction du champ Bo, dite direction longitudinale z. Lorsqu'on applique une excitation sous forme d'impulsion, c'est-à-dire un champ magnétique oscillant à la fréquence de Larmor, noté B, les moments magnétiques de spin vont alors résonner et s'écarter progressivement de cet axe longitudinal z pour aller se placer 1 o perpendiculairement à leur axe de départ en décrivant un mouvement, appelé mouvement de précession. Le champ magnétique B; permet donc de faire « basculer » les moments magnétiques de spin pour les placer dans un plan perpendiculaire à la direction du champ Bo, selon un angle de bascule FA. 15 Lorsque l'excitation est interrompue, les moments magnétiques de spin qui se sont écartés de leur axe initial reviennent vers leur position d'équilibre, c'est-à-dire l'axe z, sans cesser de tourner. Ce retour à l'équilibre est appelé relaxation. Il est alors possible de mesurer ce mouvement de rotation des spins sous la forme d'un signal oscillant de radiofréquence qui a 20 la même fréquence que le champ magnétique excitateur, c'est-à-dire la fréquence de résonance de Larmor. Ce signal est ensuite capté par l'antenne. L'un des problèmes rencontrés avec de telles antennes de tête réside dans le fait qu'elles fonctionnent en champ électromagnétique proche, 25 contrairement aux autres antennes, c'est-à-dire près des éléments physiques rayonnant. D'une façon classique, les intensités de champ magnétique Bo utilisées en imagerie médicale sont comprises entre 0,1 et 3 Tesla. A ces faibles intensités de champ magnétique, la fréquence de 30 fonctionnement correspondant à la fréquence de Larmor reste faible. A titre d'exemple, pour un champ magnétique Bo de 1,5 Tesla, la fréquence de Larmor correspondante est de 64 MHz et pour un champ magnétique Bo de 3 Tesla, la fréquence de Larmor correspondante est de 128 MHz. Généralement, les antennes utilisées pour l'examen de la tête à faible champ (entre 0,1 et 3 Tesla) sont des antennes, dites antenne volumique, du type cage à oiseau (en anglais bird cage) ou TEM (pour Traverse Electric and Magnetic en langue anglaise) telles que décrit dans les documents US 4746866, US4751464, US 5557247. Ces antennes du type cage à oiseau sont optimisées afin de produire le plus de champ magnétique avec le moins de puissance. 1 o Ces antennes, du type cage à oiseau, comportent au moins deux éléments rayonnants et sont alimentées par une seule voie, c'est-à-dire une seule amenée de puissance, avec un déphasage du signal de 90° au niveau de deux éléments rayonnants créant un champ magnétique tournant, l'amenée de puissance étant elle-même alimentée par un seul élément de 15 pilotage. Ce type d'antenne permet d'obtenir des résultats satisfaisant et notamment un champ magnétique B; relativement homogène sur toute la zone à étudier à des intensités de champ magnétique dites à bas et à moyen champ magnétique. En effet, avec des champs magnétiques appliqués 20 faibles, c'est-à-dire jusqu'à 3 Tesla, la fréquence de Larmor correspondante reste faible et les artefacts sont tolérables. Cependant, pour des intensités de champ magnétique Bo importantes, c'est-à-dire supérieures ou égales à 7 Tesla, ces antennes deviennent critiques. En effet, la fréquence oscillante du champ magnétique B, devant 25 être ajustée de façon à correspondre à la fréquence de Larmor qui est proportionnelle à l'intensité du champ magnétique statique longitudinal Bo devient alors supérieure à 298 Mhz pour le noyau de proton et la longueur d'onde effective du champ magnétique B, est relativement courte, de l'ordre de quelques dizaines de centimètres, c'est-à-dire de l'ordre de la taille de 30 l'échantillon à étudier, notamment la tête du patient. De plus, dans cette gamme de fréquences, une grande partie de la puissance électromagnétique est perdue dans l'échantillon, et notamment à la périphérie de l'échantillon ce qui implique que les zones éloignées de la périphérie de l'échantillon reçoivent moins de puissance. Dans cette gamme de fréquence, les ondes électromagnétiques ont également une forme irrégulière créant des concentrations locales du champ magnétique au centre de l'échantillon. Il en résulte que ces antennes, du type cage à oiseau, à ces valeurs de champs magnétiques élevées produisent des artéfacts importants et une inhomogénéité dans la répartition du champ d'excitation dans l'échantillon rendant ces antennes pratiquement inutilisables. En effet, dans les zones de l'échantillon où le champ d'excitation sera faible, le signal oscillant radiofréquence de relaxation sera faible, ce qui risque de poser des problèmes si cette zone est pertinente pour l'analyse de l'échantillon. Il a été développé des antennes de type réseau phasé formées par la combinaison d'une pluralité d'antennes de surface, chaque antenne de surface comportant une voie propre, ou canal, d'émission et de réception du signal. Chaque antenne de surface permet ainsi de produire une image de la région anatomique en regard de laquelle elle se trouve, les différentes images sont ensuite combinées par des algorithmes pour former l'image finale. Afin d'exploiter convenablement ces antennes de type réseau phasé, il convient que chaque canal de l'antenne soit relié à un amplificateur de puissance et que chaque amplificateur de puissance soit piloté avec une amplitude et un déphasage particulier par une interface homme/machine. Pour cela, des procédés de pilotage ont donc été développés de façon à piloter ces antennes réseau phasé consistant notamment à utiliser des ajustements d'excitation statiques et dynamiques lors des séquences d'imagerie. L'inconvénient de cette technologie est le coût élevé de réalisation qui est dû notamment à la multiplication des amplificateurs de puissance et surtout à celui des électroniques de pilotage nécessaires pour le pilotage de chaque canal de l'antenne. Ainsi, pour chaque canal de l'antenne, une chaîne de pilotage est attribuée. Afin de réduire ces coûts de réalisation d'une antenne réseau, il a été développé un procédé de pilotage permettant de réduire le nombre de chaînes de pilotage générant des signaux au moyen des matrices de Butler. Ce procédé de pilotage est mentionné notamment dans l'article « Slice-Selective RF pulses for in vivo B1 inhomogeneity mitigation at 7 Tesla using parallel RF excitation with a 16-Element coil », 2008, Kawin Setsompop. Ce procédé de pilotage se base sur l'utilisation de la matrice de Butler qui est une matrice passive dissipant peu d'énergie ce qui permet de l'implanter entre les amplificateurs et les canaux d'antenne. Cela permet de réduire le nombre d'éléments de pilotage mais également le nombre d'amplificateurs de puissance de coût élevés. D'une façon connue, cette matrice consiste en l'excitation de certains modes de propagation électromagnétiques dans un guide d'onde cylindrique. La matrice de Butler permet de convertir des signaux, de façon à ce qu'à une entrée de puissance soit associée plusieurs sorties de même amplitude, et dont la phase varie selon le mode de propagation à exciter. Cependant, ce système est un système fixe, dont les réglages des composants passifs déphaseurs sont longs à effectuer. Ils ne permettent donc aucune modularité en fonction du nombre de canaux de l'antenne, même si différentes matrices de Butler spécifiques ont été développées pour des utilisations à 2, 4, 8, ou 16 canaux et ne sont véritablement optimales que pour des systèmes cylindriques ou des objets assimilables à des guides d'onde cylindriques. Ainsi, la matrice de Butler ne permet pas non plus de s'adapter de façon optimale aux excitations s'établissant dans une tête humaine, la tête humaine n'étant pas un guide d'onde de forme cylindrique et n'étant pas non plus homogène. La figure 1 illustre schématiquement le principe du pilotage par matrice de Butler. Dans l'exemple illustré, quatre signaux issus d'éléments de pilotage la, lb, l c, Id sont amplifiés par des amplificateurs 2a, 2b, 2c, 2d puis convertis dans la matrice de Butler 3 de façon à associer une entrée à deux sorties de même amplitude mais de phase différente. Ainsi, chaque canal 4a, 4b, 4c, 4d, 4e, 4f, 4g, 4h de l'antenne produit un signal d'émission de même amplitude mais avec une phase propre. L'utilisation de la matrice de Butler pour le pilotage d'une antenne réseau phasé d'appareils d'IRM ne permet pas non plus d'obtenir une répartition du champ magnétique d'excitation B; suffisamment homogène à haut champ magnétique dans l'échantillon. Afin de garantir que les zones recevant le moins de puissance reçoivent tout de même assez de puissance pour générer un signal de retour, il est nécessaire d'augmenter fortement la puissance d'émission globale. 1 o En effet, cette technique permet de sélectionner des modes d'excitation parmi les modes proposés par la matrice de Butler (et seulement parmi ceux-là) afin de réduire le nombre de degrés de liberté du système, et donc le nombre d'éléments de puissance. Ainsi, les résultats obtenus avec ce dispositif ne permettent pas non 15 plus d'exploiter ce type d'antenne réseau de façon optimale à très haut champ magnétique, dans la mesure où les modes d'excitation de Butler ne sont pas optimaux pour un cerveau. Dans ce contexte, la présente invention vise à résoudre les problèmes mentionnés précédemment en proposant un procédé de réalisation d'un 20 module de pilotage d'une antenne haute fréquence de type réseau phasé permettant d'optimiser le signal d'émission en fonction de la région d'intérêt (ROI pour Region Of Interest en langue anglaise) à analyser tout en proposant une solution économique par la réduction du nombre d'élément de pilotage nécessaires au pilotage de l'antenne. 25 A cette fin, l'invention propose un procédé de réalisation d'un module de pilotage d'une antenne haute fréquence pour appareil d'imagerie par résonnance magnétique nucléaire comportant une pluralité n d'éléments de pilotage et une pluralité m d'éléments rayonnants, ladite pluralité m d'éléments rayonnants étant supérieure à ladite pluralité n d'éléments de 30 pilotage, ledit procédé comportant successivement : une étape de détermination d'une cartographie de distribution d'un champ magnétique B; ; - une étape de construction d'une matrice B à partir de la sélection d'une pluralité de x points parmi ladite cartographie de B déterminée lors de l'étape précédente, ladite construction de la matrice B étant réalisée à partir de xxm points Bl(i J) avec i représentant l'indice du point sélectionné variant de 1 à x , et j représentant l'indice de l'élément rayonnant ERS variant de 1 à m , de sorte que É(i,i) = B1(i,.i) ; - une étape de détermination des vecteurs singuliers de ladite matrice B par décomposition en valeurs singulières, lesdits vecteurs singuliers de ladite matrice B étant représentés par une matrice (v)k 1 o avec k représentant l'indice d'un élément de pilotage EPk variant de 1 à n ; - une étape d'agencement de moyens de couplage et/ou de moyens de déphasage formant ledit module de pilotage (10) de sorte qu'un signal So; d'un élément rayonnant ERS en sortie dudit module de pilotage est n 15 fonction de la relation Soi _ 1Vk J •Sik ; avec Sik représentant le signal k=1 d'entrée d'un élément de pilotage EPk. On désigne par le terme « élément de pilotage » un sous-ensemble électronique comprenant généralement un modulateur d'amplitude et un modulateur de phase rapides qui sont contrôlés par le système informatique 20 de l'appareil d'IRM et accessibles à l'opérateur à travers une interface homme-machine. L'invention consiste à proposer un procédé de réalisation d'un module de pilotage permettant de piloter m éléments rayonnants ER d'une antenne à voies multiples à l'aide de n éléments de pilotage EP, le nombre m 25 d'éléments rayonnants ER étant supérieur au nombre n d'éléments de pilotage EP. Le passage de n éléments de pilotage et donc de n signaux de pilotage aux m signaux à l'entrée des éléments rayonnants est réalisé au moyen du module de pilotage selon l'invention. Le module de pilotage de l'antenne réseau phasé est réalisé à partir de l'utilisation d'une méthode de décompositions en valeurs singulières permettant de définir les vecteurs singuliers de la matrice B représentant la cartographie du champ magnétique B. La méthode de décomposition en valeurs singulières (SVD pour Singular Value Decomposition en langue anglaise) d'une matrice est connue pour être un outil important de diagonalisation des matrices rectangulaires réelles et complexes. La méthode SVD est utilisée dans le domaine du traitement du signal, aux statistiques, à la météorologie et dans le domaine de la compression d'image. Mais l'utilisation de méthode de décomposition en valeurs singulière appliquée aux champs magnétiques pour l'optimisation d'un module de pilotage pilotant les signaux d'émission dans le domaine particulier de la RMN, qui est très éloigné des domaines d'utilisation classiques, est tout à fait originale. L'utilisation de cette méthode permet en effet d'obtenir un module de pilotage avec lequel on puisse connaître exactement la puissance effectivement transmise à chaque élément rayonnant. Ainsi, lors de la conception du module, il est possible de choisir, parmi différents vecteurs singuliers, les vecteurs singuliers de plus fortes excitations, de façon à obtenir un module de pilotage qui fournit moins de puissance aux éléments rayonnants lors du pilotage de l'antenne et qui améliore le rendement d'une telle antenne réseau phasé. De plus, ce module de pilotage obtenu par le procédé selon l'invention permet de piloter plus librement une antenne quel que soit le nombre de canaux de l'antenne, contrairement à l'utilisation d'un procédé de pilotage utilisant le principe des matrices de Butler. Le procédé selon l'invention peut également présenter une ou plusieurs des caractéristiques ci-dessous, considérées individuellement ou selon toutes les combinaisons techniquement possibles : - ladite étape d'agencement de moyens de couplage et/ou de moyens de déphasage formant ledit module de pilotage est réalisée au moyens de coupleurs, dits coupleur branchline, agissant en tant que moyens de couplage et/ou moyens de déphasage ; ledit procédé comporte une étape de décomposition de ladite matrice V en une pluralité N+1 de matrice selon la relation : V= VNXVN_1X...XVIXVO - ladite décomposition en valeurs singulières de l'étape de détermination des vecteurs singuliers de ladite matrice B est réalisée au moyen de la relation BV = UE ; avec : U matrice unitaire de même dimension que la matrice B , et dont les m colonnes sont unitaires ; E matrice diagonale de dimension mxm et dont les coefficients diagonaux sont des réels positifs ; - ladite étape de détermination des vecteurs singuliers de ladite matrice B est réalisée au moyen de la relation B*B = VEZV* où * désigne l'opérateur hermitien ; - ladite étape de détermination des vecteurs singuliers de ladite matrice B comporte une sous-étape de sélection de n vecteurs singuliers de ladite matrice B ; - ladite sous-étape de sélection de n vecteurs singuliers de ladite matrice B est réalisée par la sélection des n valeurs singulières de plus forte valeur représentées par les coefficients diagonaux de la matrice E - ladite étape de détermination d'une cartographie de distribution d'un champ magnétique B est réalisée par d'un champ électromagnétique B1 via des moyens logiciels ; - ladite étape de détermination d'une cartographie de distribution d'un champ magnétique B est réalisée par mesure. D'autres caractéristiques et avantages de l'invention ressortiront plus clairement de la description qui en est donnée ci-dessous, à titre indicatif et nullement limitatif, en référence aux figures annexées, parmi lesquelles : - la figure 1 illustre un schéma d'un dispositif de pilotage d'une antenne haute fréquence, de type réseau phasé, utilisant le principe des matrices de Butler pour le pilotage de l'antenne ; - la figure 2 illustre un schéma d'une chaîne de pilotage d'une antenne haute fréquence, de type réseau phasé ; - la figure 3 illustre un schéma synoptique présentant les principales étapes du procédé de pilotage selon l'invention ; - la figure 4 représente schématiquement un coupleur de type branchline utilisé lors de la quatrième étape du procédé de réalisation 1 o illustré à la figure 3; - la figure 5 illustre un exemple de valeurs pour une matrice dite matrice supérieureVtu ; - la figure 6 illustre un exemple d'agencement de coupleurs branchline illustré à la figure 4 obtenu par le procédé de réalisation d'un module 15 de pilotage selon l'invention. Dans toutes les figures, les éléments communs portent les mêmes numéros de référence. La figure 1 a déjà été décrite précédemment en référence à la présentation de l'état de la technique. 20 La figure 2 illustre un schéma d'une chaîne de pilotage 100 permettant le pilotage et l'alimentation d'une antenne réseau phasé 300. La chaine de pilotage 100 est formée par une pluralité n d'éléments de pilotage EPk, k étant un indice variant de 1 à n , un module de pilotage 10 obtenu par le procédé de réalisation selon l'invention, relié à la pluralité m de canaux ou 25 d'éléments rayonnants ERS de l'antenne 300, j étant un indice variant de 1 à m. La chaîne de pilotage comporte également des amplificateurs de puissance Ak positionnés entre les éléments de pilotage EPk et le module de pilotage 10. Selon l'invention, le module de pilotage 10 permet de piloter une 3o antenne réseau phasé avec moins d'éléments de pilotage EP que d'éléments rayonnants ER de l'antenne (i.e. m> n ). A titre d'exemple, le procédé de réalisation et le module de pilotage issu du procédé seront décrits pour une antenne comportant quinze canaux indépendants d'émission et quinze éléments rayonnants (m = 15) pilotés par huit éléments de pilotage (n = 8). On rappelle que par la désignation de canal, on entend chaque entrée de puissance de l'antenne communiquant avec un élément rayonnant. Selon l'invention, chaque élément de pilotage EPk (i.e. EP1, EP2, ... EPä) est formé par un sous-ensemble électronique (non représenté) comprenant généralement un modulateur d'amplitude et un modulateur de 1 o phase rapides qui sont contrôlés par le système informatique de l'appareil d'IRM et accessibles à l'opérateur à travers une interface homme/machine. Ainsi dans notre exemple de réalisation illustré à la figure 2, le module de pilotage 10 comporte à son entrée les éléments de pilotage EPk émettant des signaux Sik(i.e Si1, Si2,..., Si8). 15 Le module de pilotage 10 est relié à sa sortie à chaque élément rayonnant ERS de l'antenne qui reçoit un signal Soi (i.e Sol, So2,..., So15). Le module de pilotage a donc pour effet de transformer linéairement une pluralité n de signaux d'entrée Si provenant de n élément de pilotage en une pluralité m de signaux de sortie So pour alimenter m éléments 20 rayonnants ER. La transformation est réalisée au moyen d'une matrice V de sorte que : So =VSi où : 25 - l'ensemble des signaux Si est représenté par une matrice de dimension nxl, c'est-à-dire que le kième élément de la matrice correspond au signal envoyé, dans le module de pilotage 10, par le kième élément de pilotage EPk, avec k variant de 1 à n ; - l'ensemble des signaux So est représenté par une matrice de 30 dimension mxl, c'est-à-dire que le jième élément de la matrice correspond à la puissance envoyée dans le jième élément rayonnant, avec j variant de 1 à m ; - la matrice V est une matrice de dimension mxn unitaire en n. Ainsi dans l'exemple illustré à la figure 2, le module de pilotage 10 comporte des moyens de stockage comportant une matrice de transformation V de dimension 15x8. La figure 3 illustre un schéma synoptique représentant les principales étapes du procédé de réalisation du module de pilotage 10 selon l'invention. Le procédé selon l'invention consiste dans un premier temps à déterminer les coefficients de la matrice V de transformation obtenu au moyen d'une méthode dite de décomposition en valeurs singulières (SVD 1 o pour Singular Value Decomposition en langue anglaise) puis dans un second temps à construire physiquement cette matrice par un ensemble de coupleurs. La méthode de décomposition en valeurs singulières (SVD pour Singular Value Decomposition en langue anglaise) d'une matrice est généralement utilisée dans le domaine du traitement du signal, aux 15 statistiques, à la météorologie et dans le domaine de la compression d'image. Mais l'utilisation de méthode de décomposition en valeurs singulière appliquée aux champs magnétiques pour l'optimisation d'un module de pilotage pilotant les signaux d'émission dans le domaine particulier de la RMN, qui est très éloigné des domaines d'utilisation classiques, est tout à fait 20 originale. Ainsi, le procédé selon l'invention comporte une première étape 210 d'établissement d'une cartographie de la distribution du champ magnétique d'excitation B; dans l'antenne lorsqu'une « charge » est positionnée à l'intérieur de l'antenne. 25 Cette cartographie de la répartition du champ magnétique est réalisée indifféremment par simulation ou par mesure sur un prototype dans un scanner IRM. A titre d'exemple, la simulation est effectuée au moyen du logiciel de simulation numérique HFSS de la société Ansoft. Le principe de la simulation 30 consiste à étudier systématiquement la distribution de la composante magnétique B, telle que : B1 = B1 B1Y les indices « x » et « y » se réfèrent à un repère cartésien fixe, lié au laboratoire, selon les axes X et Y (non représentés) perpendiculaire à l'axe longitudinal Z de l'antenne (non représenté), l'axe longitudinal Z étant la direction de l'aimantation globale du champ Bo. La simulation du champ magnétique B1 permet de déduire B; au moyen d'un calcul intermédiaire, B; correspondant au phaseur du champ magnétique B1 décomposé dans un repère tournant dans le même sens des spins selon la relation : Bl = Bi .e LHCP 1 OU eLHc. = La détermination de la cartographie représentant la distribution du champ magnétique B; peut également être réalisée expérimentalement par mesure à partir d'un prototype. Pour cela, il suffit de déterminer directement 15 la cartographie de la composante magnétique B, notamment au moyen d'une méthode d'acquisition, dite méthode AFI (pour Actual Flip-angle Imaging en langue anglaise) [Vasily L. Yarnykh, Actual Flip-angle Imaging in the Pulsed Steady State. Magn Reson Led 2007 ; 57 : 192-200] connue du domaine de l'imagerie médicale. 20 Pour réaliser une cartographie représentant la distribution du champ magnétique par mesure, il est nécessaire que le prototype ou l'antenne de laboratoire comporte autant d'éléments de pilotage EP que d'élément rayonnant ER (i.e. m =n ). Une deuxième étape 220 du procédé consiste en une étape de 25 construction d'une matrice B à partir de la sélection d'une pluralité de points Bl(t J) parmi la cartographie du champ magnétique complexe B déterminée lors de l'étape précédente 210, la pluralité de points Bl(t J) étant sélectionnée 1 J par un maillage prédéterminé de ladite cartographie par un opérateur au moyen d'une interface homme/machine. Le maillage de la cartographie B; est choisi par l'opérateur en fonction de la résolution désirée ainsi qu'en fonction de la fréquence de 5 résonnance dépendante du champ magnétique Bo appliqué. Avantageusement, un maillage de 5mm, c'est-à-dire un maillage comportant une pluralité de points dont la distance entre deux points est de 5mm, permet d'obtenir une résolution satisfaisante pour une antenne de type tête. 10 Au-delà d'un maillage de 5mm, les résultats risquent de perdre en précision et en deçà d'un maillage de 5mm, les éléments de calculs risquent de devenir conséquents et trop volumineux à traiter par rapport au gain de résolution apporté. Pour cela, chaque voxel extrait de la cartographie du champ 15 établie lors de l'étape 210 précédente, est indexé par la variable ide sorte que chaque voxel correspond à un point Bl(i,.i) On rappelle que l'indice j correspond à l'indice des éléments rayonnant ERS de l'antenne avec j variant de 1 à m et que l'indice k correspond à l'indice des éléments de pilotage EPk avec k variant de 1 à n. 20 On rappelle que le voxel est un pixel en trois dimensions utilisé par la représentation d'espace en trois dimensions notamment en imagerie médicale par traitement numérique de coupes en deux dimensions (2D) issues des appareils d'IRM. Pour chaque voxel i correspond une position xi,yi,zi à laquelle est 25 associé un champ magnétique ponctuel Bl(t j) émis par un élément rayonnant ERS, l'ordre d'indexation des voxels n'ayant aucune influence sur le résultat final. Ainsi, l'élaboration de la matrice est réalisée en fonction de la cartographie réalisée dans l'étape 210 précédente de sorte que : 30 B(i,i) - BI(i,i) La matrice B possède autant de lignes que de voxels sélectionnés dans la région d'intérêt (ROI) dans la cartographie du champ magnétique oscillant B lors de l'étape 210 d'établissement d'une cartographie de la distribution du champ magnétique d'excitation et autant de colonnes (i.e. m 5 colonnes) que d'éléments rayonnants ER de l'antenne 300. Typiquement, dans notre exemple illustré à la figure 2, quinze éléments rayonnants ER sont utilisés et le champ magnétique oscillant B est cartographié avec une résolution en trois dimensions de 32x32x32, ce qui donne une matrice B possédant 32768 lignes et 15 colonnes. 1 o La troisième étape 230 du procédé est une étape de détermination des modes, c'est-à-dire des vecteurs singuliers, de la matrice B par l'utilisation d'une méthode de décomposition en valeurs singulières, dite méthode SVD. On rappelle qu'on entend par mode SVD les vecteurs singuliers de la 15 décomposition en valeurs singulières. Ces modes n'ont pas de lien avec les modes électromagnétiques d'une antenne volumique ou d'une cavité résonnante. Grâce au principe de décomposition en valeurs singulières, on associe à la matrice B les matrices U , E , V telles que : 20 BV =LIE avec : U matrice unitaire de même dimension que la matrice B , dont les m colonnes sont unitaires (i.e. LIU* =Im) ; E matrice diagonale de dimension mxm et dont les coefficients 25 diagonaux sont des réels positifs. V matrice carré de dimension mxm. Les coefficients diagonaux de la matrice E représentent donc les valeurs singulières de la matrice B , et les colonnes de la matrice U et de la matrice V sont respectivement les vecteurs singuliers à gauche et à droite 30 pour les valeurs singulières correspondantes, les vecteurs singuliers à droite de la matrice V étant les modes à déterminer. Ainsi dans notre exemple, la matrice U est une matrice de dimensions 32768 x 15, et les matrices V et E sont des matrices de dimensions 15 x 15. Cette décomposition est unique, à une permutation près des colonnes 5 de V, deU et de E. Les propriétés de la matrice V et de la matrice U permettent à partir de la relation : (BV )*(BV) = (uE)*(uE) , d'obtenir la nouvelle relation suivante : 10 B*B = VE2v* On en déduit alors que la matrice V est également un vecteur propre de la matrice B*B, et les valeurs propres de cette matrice B*B sont les coefficients de E2 . La troisième étape 130 comporte une sous-étape de sélection de n 15 valeurs singulières de la matrice B parmi les coefficients diagonaux de plus forte valeur de la matrice E . Ainsi, on associe pour chaque mode k de la matrice B une valeur singulière ak correspondant à un coefficient diagonal de la matrice E. L'efficacité en puissance du k1ème mode est alors proportionnelle à ak2, c'est 20 pourquoi la sélection des n plus fortes valeurs garantie la plus grande efficacité des modes. Ainsi pour huit éléments de pilotage EP1, EP2, EPB, les huit valeurs singulières al,...,a8 de plus forte valeur parmi les quinze sont sélectionnées, a, étant supérieure à a2 et ainsi de suite jusqu'à a8. 25 La troisième étape 130 peut comporter également une sous-étape de classement des valeurs singulières réelles par ordre croissant par permutation des vecteurs singuliers. Chaque valeur singulière al,...,a8 est associée aux vecteurs singuliers de la matrice réduite U8 à n colonnes (n = 8), U8 étant une matrice de 30 dimension 32768 x 8. Les colonnes de la matrice réduite U8 représentent la distribution du champ magnétique B pour chacun des modes. A la matrice réduite U8 est associée la matrice réduite V8 de dimension 15 x 8 (i.e. mxn) dont les colonnes représentent les modes recherchés. Ainsi, on utilise les vecteurs singuliers, ou modes, de la matrice B , correspondant aux colonnes de la matrice V , pour la transformation des n signaux Si des éléments de pilotage EP en m signaux So alimentant les m éléments rayonnants ER de l'antenne 300, de sorte qu'un élément rayonnant ERS parmi les m élément rayonnant reçoit un signal Soi équivalent à : Soi = IV, ..Sik , k=1 1 o avec j variant de 1 à met k variant de 1 à n. Ainsi, les signaux Soi sont caractérisés par une amplitude proportionnelle à la somme des Vk I et dont la phase est égale à la somme des arg(VkJ). L'intérêt de cette décomposition est l'indépendance de chacun des 15 modes déterminés. En effet, la construction de la matrice V assure que les modes sont orthogonaux au champ magnétique excitateur B; La matrice V étant une matrice unitaire, telle que VV* =I15, elle permet de garantir que la puissance de sortie du module de pilotage 10 est équivalente à la puissance d'entrée selon la démonstration suivante : 20 So z = So:So = (vs, )V Si = Si*V *V Si = Sj Si = si Ainsi, la transformation est réalisée sans perte d'énergie. Cette décomposition est également avantageuse dans la mesure où le choix des modes selon la troisième étape 230 se fait en fonction de leurs rendements énergétiques. En effet, comme décrit précédemment, les valeurs 25 singulières de la matrice B , c'est-à-dire les coefficients diagonaux de la matrice E sont proportionnelles à V f f f 1131 2 dV pour chaque mode étudié. tete On en déduit alors pour une puissance d'entrée arbitraire le rendement énergétique de la charge introduite dans l'antenne, dans notre cas la tête du patient, selon la relation : Em = fff -,uI8112 dV Les champs B; étant orthogonaux pour chaque mode, l'énergie magnétique excitatrice totale reçue, par le patient dans la zone ROI, est égale à la somme des énergies excitatrices de chaque mode, si plusieurs modes sont excités. Ce calcul de rendement énergétique peut être réalisé aussi bien à partir d'une cartographie de distribution de B; obtenue par simulation ou par mesure, puisque seuls les éléments nécessaires à ce calcul sont les valeurs de B; et la différence de phase d'un élément rayonnant ER à un autre. La quatrième étape 240 consiste en une étape de d'agencement de moyens de couplage et de moyens de déphasage de façon à reproduire la transformation de signaux de la matrice V déterminée lors de l'étape précédente 230. On rappelle que la matrice V obtenu est une matrice unitaire en n, c'est-à-dire que la matrice V vérifie la relation : V*V = In , où In est la matrice identité de dimension nxn. Ainsi, en introduisant des sorties intermédiaires Soi ,l de rang l avec 1 variant de 0 à N et des matrices unitaires intermédiaires V1 (1>0) de dimension mxm et Vo (1=0) de dimension mxn telles que : Sol o =Si pour 1=0 Sol 1 =V1So1 1_1 pour N 1 > O SOI N = Soi Si les matrices intermédiaires V1 sont unitaires, on obtient alors VN xVN_1 x...xV1 xVo qui est unitaire de dimension mxn. Ce principe de décomposition en une pluralité de matrice est utilisé pour la construction de la matrice V de sorte que : So = VSi = VN XVN_1 X...XVI XVOSi Ensuite chaque matrice Vo,V1,...,VN est réalisée au moyen de coupleurs de type coupleurs branchline tels qu'illustré à la figure 4. Les coupleurs dits branchline sont des coupleurs directionnels, généralement utilisés par une répartition à 3 décibels de l'énergie. D'une façon générale, le coupleur branchline comporte deux entrées et deux sorties et permet donc de transformer les entrées I en sortie z selon la relation : j cos e sine sine j cos e 1 2 I2 où 8 peut être choisi arbitrairement. Il est également possible d'associer au coupleur branchline un déphasage à partir d'une ligne de transmission provoquant un retard de l'onde. Ce déphasage peut être représenté par la matrice unitaire : e 0 0 e302 Le produit est donc une matrice unitaire : j cos e e'q sine sine jcosO La première matrice Vo doit être de dimension mxn, soit dans notre exemple de réalisation illustré à la figure 2 de dimension 15 x 8. Pour obtenir 15 sorties à partir de 8 entrées avec des coupleurs branchline tels que décrit précédemment, il suffit d'utiliser 7 coupleurs branchline dont chaque coupleur comporte une entrée chargée avec une résistance, par exemple de 50 O. Dans ce cas, la matrice de transformation de chaque coupleur branchline devient pour une entrée I : -01-0 _ Z c(0,01,02)= jcosO sine ùjcosO.1 sine .1 Ainsi, en connectant un coupleur branchline ainsi conçu à 7 des 8 entrées, on obtient 15 sorties, la huitième entrée étant connectée directement à une sortie dans l'utilisation d'un coupleur. En connectant ainsi, la première entrée à la première et à la neuvième sortie, la deuxième entrée à la deuxième et dixième sortie ainsi de suite et la huitième entrée directement à la huitième sortie, on obtient la première matrice Vo : jcos91 0 0 0 0 0 0 0 0 j cos 92 0 0 0 0 0 0 0 0 jcos93 0 0 0 0 0 O O O jcos94 O O O O O O O O j cos 95 O O O O O O O O jcos96 O O O O O O O O j cos 97 O Vo ( 91,...,97 ) = O O O O O 0 0 1 sin 91 0 0 0 0 0 0 0 0 sin 92 0 0 0 0 0 0 0 0 sin 93 0 0 0 0 0 0 0 0 sin 94 0 0 0 0 0 0 0 0 sin 95 0 0 0 0 0 0 0 0 sin 96 0 0 0 0 0 0 0 0 sin 9, 0 Les matrices V1,...,VN_1,VN sont réalisées également au moyen de 1 o coupleurs dits coupleurs branchline. Pour la réalisation de ces matrices, les coupleurs ne comportent pas de charge à leur entrée mais créent un déphasage 8 en sortie. Chacune des matrices V1,...,VN_1,VN peut être réalisée par l'alternance de deux matrices que l'on nommera matrice supérieure Vl et matrice 15 inférieure V1d , 1 variant de 1 à N. La matrice supérieure Vl est construite de sorte qu'un coupleur connecte la première entrée et la deuxième entrée, la troisième entrée et la quatrième entrée et ainsi de suite jusqu'à qu'à la treizième entrée et la quatorzième entrée et laisse la quinzième entrée inchangée. La matrice inférieure V1d est construite de sorte qu'un coupleur connecte la deuxième entrée et la troisième entrée, la quatrième entrée et la cinquième entrée, etc. jusqu'à qu'à la quatorzième entrée et la quinzième entrée et laisse la première entrée inchangée. Cette configuration permet d'éviter de croiser les coupleurs et facilite la conception de la matrice V du module de pilotage 10. Une telle conception est illustrée à titre d'exemple à la figure 6. Ainsi, dans le cas où le nombre m d'éléments rayonnants ER est pair alors la matrice supérieure 1/tu applique une transformation sur toutes les entrées couplant les paires 1-2, 3-4, etc. ; et la matrice inférieure Vtd n'appliquera pas de transformation sur la première et la dernière entrée, les coupleurs reliant les paires 2-3, 4-5, etc. Au final, le placement des coupleurs permet d'obtenir une matrice V1,...,VN_1,VN dont chacune dépend de 7 paramètres 6 de couplage et quinze paramètres c de déphasage. Un exemple de valeurs pour la réalisation de la matrice supérieure Vtu est illustré en référence à la figure 5. Ainsi, les matrices Vo,V1,...,VN sont définies et permettent la fabrication de la matrice V en vérifiant les hypothèses précédentes. Pour chacune des matrices Vo,V1,...,VN , il est toutefois nécessaire de déterminer les valeurs de e et de c permettant d'obtenir la matrice V désirée, c'est-à-dire la matrice V déterminée lors des étapes précédentes. Selon un mode avantageux de l'invention, la matrice V est approximée par une matrice approchée Vapp permettant ainsi de limiter les coûts de réalisation de cette matrice. Ainsi, la matrice approchée Vapp possède les mêmes propriétés de matrice unitaire que la matrice V mais les valeurs sont approchées à une erreurs près de sorte que : Vapp -V2 Ainsi selon notre mode de réalisation, la construction de la matrice V nécessite alors 91 coupleurs et 195 éléments de déphasage. Les 286 paramètres réels devront être optimisés par des algorithmes d'optimisation, comme par exemple par les algorithmes d'optimisation de logiciel de type MatLab. L'invention a été particulièrement décrite pour le pilotage d'une antenne haute fréquence d'IRM pour l'examen de la tête ; toutefois, le procédé de pilotage et notamment la méthode SVD décrite est également applicable pour n'importe quelle partie du corps et n'importe quel patient, enfant, adulte, ou encore une personne de forte corpulence, etc. Les champs magnétiques B, mesurés ou simulés dépendent évidemment de la partie du corps, ou de l'objet, introduit dans l'antenne. Les résultats du calcul varient en fonction de la partie du corps introduite mais également en fonction du patient. Le procédé selon l'invention utilisant une méthode de décomposition en valeurs singulières permet de faciliter l'optimisation du résultat de l'ajustement de l'excitation statique et dynamique et notamment la qualité de l'ajustement des signaux par des algorithmes d'optimisation du type Gauss Newton ou encore de tout type d'algorithme requérant une valeur initiale pour l'optimisation. Dans ce cas, l'optimisation peut être réalisée dans un premier temps pour le premier mode, puis progressivement pour les modes suivants. REVENDICATIONS1. Procédé de réalisation d'un module de pilotage (10) d'une antenne haute fréquence (300) pour appareil d'imagerie par résonnance magnétique comportant une pluralité n d'éléments de pilotage (EP1, EP2,..., EPä) et une pluralité m d'éléments rayonnants (ER1, ER2,..., ERm), ladite pluralité m d'éléments rayonnants étant supérieure à ladite pluralité n d'éléments de pilotage (EP1, EP2, ...EPä), ledit procédé comportant successivement : - une étape (210) de détermination d'une cartographie de distribution d'un champ magnétique Bi ; - une étape (220) de construction d'une matrice B à partir de la sélection d'une pluralité de x points parmi ladite cartographie de B déterminée lors de l'étape précédente, ladite construction de la matrice B étant réalisée à partir de xxm points 131+(i ,i) avec i représentant l'indice du point sélectionné variant de 1 à x , et j représentant l'indice de l'élément rayonnant ERS variant de 1 à m , de sorte que B(~,i) = B1(i,.i) ; - une étape (230) de détermination des vecteurs singuliers de ladite matrice 13- par décomposition en valeurs singulières, lesdits vecteurs singuliers de ladite matrice B étant représentés par une matrice (v)k J avec k représentant l'indice d'un élément de pilotage EPk variant de 1 à n ; - une étape (240) d'agencement de moyens de couplage et/ou de moyens de déphasage formant ledit module de pilotage (10) de sorte qu'un signal Soi d'un élément rayonnant ERS en sortie dudit module de pilotage est fonction de la relation : n Soi = 1 Vkj .Sik k=1avec Sik représentant le signal d'entrée d'un élément de pilotage EPk. 2. Procédé de réalisation d'un module de pilotage (10) d'une antenne haute fréquence (300) pour appareil d'imagerie par résonnance magnétique selon la revendication 1 caractérisé en ce que ladite étape (240) d'agencement de moyens de couplage et/ou de moyens de déphasage formant ledit module de pilotage (10) est réalisée au moyens de coupleurs, dits coupleur branchline, agissant en tant que 1 o moyens de couplage et/ou moyens de déphasage. 3. Procédé de réalisation d'un module de pilotage (10) d'une antenne haute fréquence (300) pour appareil d'imagerie par résonnance magnétique selon l'une des revendications 1 à 2 caractérisé en ce qu'il 15 comporte une étape de décomposition de ladite matrice V en une pluralité N+1 de matrice selon la relation : V= VNXVN_1X...XVIXVO 4. Procédé de réalisation d'un module de pilotage (10) d'une antenne 20 haute fréquence (300) pour appareil d'imagerie par résonnance magnétique selon l'une des revendications 1 à 3 caractérisé en ce que ladite décomposition en valeurs singulières de l'étape (230) de détermination des vecteurs singuliers de ladite matrice B est réalisée au moyen de la relation : 25 BV = UE avec : U matrice unitaire de même dimension que la matrice B , et dont les m colonnes sont unitaires ; E matrice diagonale de dimension mxm et dont les coefficients 30 diagonaux sont des réels positifs. 5. Procédé de réalisation d'un module de pilotage (10) d'une antenne haute fréquence (300) pour appareil d'imagerie par résonnance magnétique selon la revendication 4 caractérisé en ce que ladite l'étape (230) de détermination des vecteurs singuliers de ladite matrice B est réalisée au moyen de la relation : B*B = VEZV* . 6. Procédé de réalisation d'un module de pilotage (10) d'une antenne haute fréquence (300) pour appareil d'imagerie par résonnance 1 o magnétique selon l'une des revendications 1 à 5 caractérisé en ce que ladite étape (230) de détermination des vecteurs singuliers de ladite matrice B comporte une sous-étape de sélection de n vecteurs singuliers de ladite matrice B . 15 7. Procédé de réalisation d'un module de pilotage (10) d'une antenne haute fréquence (300) pour appareil d'imagerie par résonnance magnétique selon la revendication 6 caractérisé en ce que ladite sous- étape de sélection de n vecteurs singuliers de ladite matrice 13- est réalisée par la sélection des n valeurs singulières de plus forte valeur 20 représentées par les coefficients diagonaux de la matrice E . 8. Procédé de réalisation d'un module de pilotage (10) d'une antenne haute fréquence (300) pour appareil d'imagerie par résonnance magnétique selon l'une des revendications 1 à 7 caractérisé en ce que 25 ladite étape (210) de détermination d'une cartographie de distribution d'un champ magnétique B est réalisée par d'un champ électromagnétique B, via des moyens logiciels. 9. Procédé de réalisation d'un module de pilotage (10) d'une antenne 30 haute fréquence (300) pour appareil d'imagerie par résonnance magnétique selon l'une des revendications 1 à 8 caractérisé en ce queladite étape (210) de détermination d'une cartographie de distribution d'un champ magnétique B; est réalisée par mesure. 10. Module de pilotage (10) obtenu par le procédé de réalisation selon l'une des revendications 1 à 9 caractérisé en ce qu'il comporte un agencement de moyens de couplage et de moyens de déphasage de sorte qu'un signal Soi de sortie dudit module de pilotage (10), avec j représentant l'indice d'un signal de sortie parmi une pluralité m de signaux de sortie, est fonction de la relation : n SOi = Vk,j .Sik k=1 avec : Sik représentant le signal d'entrée dudit module de pilotage (10).15