EP4308949A1 - Method for generating multi-ray spin-excitation sequences and application thereof to magnetic resonance imaging - Google Patents

Method for generating multi-ray spin-excitation sequences and application thereof to magnetic resonance imaging

Info

Publication number
EP4308949A1
EP4308949A1 EP22717061.0A EP22717061A EP4308949A1 EP 4308949 A1 EP4308949 A1 EP 4308949A1 EP 22717061 A EP22717061 A EP 22717061A EP 4308949 A1 EP4308949 A1 EP 4308949A1
Authority
EP
European Patent Office
Prior art keywords
pulses
nuclear spins
excitation
sequence
functions
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
EP22717061.0A
Other languages
German (de)
French (fr)
Inventor
Vincent Gras
Franck MAUCONDUIT
Nicolas Boulant
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
Original Assignee
Commissariat a lEnergie Atomique CEA
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Commissariat a lEnergie Atomique CEA, Commissariat a lEnergie Atomique et aux Energies Alternatives CEA filed Critical Commissariat a lEnergie Atomique CEA
Publication of EP4308949A1 publication Critical patent/EP4308949A1/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/5659—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the RF magnetic field, e.g. spatial inhomogeneities of the RF magnetic field
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/482—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory
    • G01R33/4822—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory in three dimensions
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
    • G01R33/4833—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
    • G01R33/4835—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices of multiple slices
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5611—Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
    • G01R33/5612—Parallel RF transmission, i.e. RF pulse transmission using a plurality of independent transmission channels

Abstract

Method for determining a global and regular parameterisation for a family of spin-excitation pulse sequences in magnetic resonance imaging, each pulse sequence of said family being a multi-ray sequence suitable for selectively exciting nuclear spins in a section (CP) of a volume of interest (ROI) of a body (C) immersed in a static magnetic field (B0) and comprising radio-frequency pulses (IRF1, IRF2, IRF3) at a Larmor frequency of said nuclear spins alternated with magnetic-gradient pulses (IG1, IG2, IG3). The global parameterisation minimises a function representative of a mean absolute deviation from a setpoint of the excitation of the nuclear spins, the mean being computed for said volume of interest and for all possible orientations and positions of said sections. It makes it possible to design straightforwardly a selective excitation sequence for a section of arbitrary orientation and position.

Description

DESCRIPTION DESCRIPTION
Titre de l’invention : Procédé de génération de séquences d’excitation de spin de type multi-rayon et son application à l’imagerie par résonance magnétiqueTitle of the invention: Method for generating multi-ray type spin excitation sequences and its application to magnetic resonance imaging
[0001] L’invention relève du domaine de l’imagerie par résonance magnétique (IRM) et concerne plus particulièrement la phase d’excitation des spins nucléaires. Plus particulièrement encore, elle concerne la conception et l’application de séquences d’impulsions d’excitation de spins nucléaires (« séquences d’excitation ») sélectives de type « multi-rayon » (« multi-spokes » ou « fast-kz » dans la littérature en langue anglaise). Elle s’applique notamment à l’IRM à champ magnétique élevé (plusieurs Teslas, par exemple 3 Teslas (T) ou plus voire 5 Teslas ou plus ou même 7 Teslas ou plus). The invention relates to the field of magnetic resonance imaging (MRI) and relates more particularly to the phase of excitation of nuclear spins. More particularly, it relates to the design and application of selective nuclear spin excitation pulse sequences (“excitation sequences”) of the “multi-spoke” (“multi-spokes” or “fast-kz”) type. in English-language literature). It applies in particular to high magnetic field MRI (several Teslas, for example 3 Teslas (T) or more or even 5 Teslas or more or even 7 Teslas or more).
[0002] En IRM, le profil d’excitation des spins est subordonné à la distribution du champ radiofréquence d’excitation des spins dans la partie du corps à imager. A bas champ magnétique (B0<1T), la fréquence de résonance magnétique (fréquence de Larmor) du spin du proton étant faible, une distribution très uniforme et relativement indépendante de l’objet à imager peut être obtenue avec un design adapté d’antenne RF de transmission. A champ magnétique plus élevé, la longueur d’onde d’un rayonnement à la fréquence de Larmor devient comparable aux dimensions caractéristiques du corps à imager ; il en résulte des effets d’interférence qui provoquent des inhomogénéités importantes. Cela se vérifie autour de 7T pour l’imagerie cérébrale et dès 3T pour celle de l’abdomen. Dans ces conditions, il n’est plus possible d’assurer une excitation homogène au moyen d’un design adapté de l’antenne d’émission et il est nécessaire d’avoir recours à des techniques alternatives. [0002] In MRI, the spin excitation profile depends on the distribution of the spin excitation radiofrequency field in the part of the body to be imaged. At low magnetic field (B0<1T), the magnetic resonance frequency (Larmor frequency) of the proton spin being low, a very uniform and relatively independent distribution of the object to be imaged can be obtained with a suitable antenna design transmitting RF. At a higher magnetic field, the wavelength of radiation at the Larmor frequency becomes comparable to the characteristic dimensions of the body to be imaged; this results in interference effects that cause significant inhomogeneities. This is verified around 7T for cerebral imaging and from 3T for that of the abdomen. Under these conditions, it is no longer possible to ensure homogeneous excitation by means of a suitable design of the transmitting antenna and it is necessary to resort to alternative techniques.
[0003] Une première technique consiste à remplacer l’antenne volumique RF de transmission par un réseau d’antennes dont chaque élément peut être excité par un signal qui lui est propre, indépendamment des autres éléments du réseau. On parle de transmission parallèle ou, en anglais, de « RF-shimming ». [0003] A first technique consists in replacing the volumic RF transmission antenna with an array of antennas, each element of which can be excited by a signal which is specific to it, independently of the other elements of the array. This is referred to as parallel transmission or, in English, “RF-shimming”.
[0004] Une deuxième technique, dite du contrôle optimal, touche directement à la dynamique des spins (équation de Bloch) et vise à exposer les spins à champs B (champ radiofréquence émis) et à des gradients de champ magnétique variables dans le temps de manière à créer dans une région d’intérêt un angle de bascule des spins nucléaires conforme à une consigne, à une erreur près que l’on cherche à minimiser. [0004] A second technique, called optimal control, directly affects the dynamics of the spins (Bloch's equation) and aims to expose the spins to B fields (emitted radiofrequency field) and to variable magnetic field gradients over time of so as to create in a region of interest an angle of tilting of the nuclear spins in accordance with a setpoint, except for an error that we seek to minimize.
[0005] Il est possible en IRM de restreindre l’étendue de l’excitation des spins à une coupe dont la position, l’épaisseur et l’orientation peuvent être arbitrairement choisies, ce grâce à l’application simultanée de gradient de champ magnétique et d’impulsions radiofréquence bien choisis. On parle ici d’excitation sélective. Cette technique permet - mais ce n’est pas là le seul intérêt - de partager l’espace en coupes et ainsi de ramener le problème de l’homogénéisation de l’excitation à celui de l’homogénéisation de l’excitation des spins à l’intérieur de chaque coupe, ce indépendamment de ce qui se passe dans les autres coupes. [0005] It is possible in MRI to restrict the extent of the excitation of the spins to a section whose position, thickness and orientation can be chosen arbitrarily, thanks to the simultaneous application of a magnetic field gradient and well-chosen radio frequency pulses. We are talking here about selective excitation. This technique makes it possible - but this is not the only advantage - to divide the space into sections and thus to reduce the problem of the homogenization of the excitation to that of the homogenization of the excitation of the spins at inside each cup, regardless of what happens in the other cups.
[0006] Ces techniques d’homogénéisation peuvent être utilisées séparément ou en combinaison. Cela conduit à l’utilisation d’une séquence globale multi-coupes constituée de plusieurs séquences d’excitation sélectives, chacune optimisée pour exciter de manière aussi homogène que possible une coupe du volume d’intérêt. Ces séquences peuvent exploiter, ou pas, la transmission parallèle. Voir par exemple [Cao 2015] [0006] These homogenization techniques can be used separately or in combination. This leads to the use of a global multi-slice sequence consisting of several selective excitation sequences, each optimized to excite as homogeneously as possible a slice of the volume of interest. These sequences may or may not use parallel transmission. See for example [Cao 2015]
[0007] La définition spatiale de la ou des coupes à imager, c’est-à-dire la position et l’inclinaison de coupe, peut donner lieu à une optimisation spécifique des propriétés des signaux RF et des signaux de gradient de champ magnétique statique. Notant N le nombre de coupes considérées, l’implémentation de la séquence globale multi- coupes requiert donc la définition de N ensembles indépendants de signaux (1 signal complexe par canal de transmission au niveau du réseau d’antenne de transmission RF and 3 signaux réels pour les 3 axes de gradient de champ magnétique). Cette approche permet effectivement d’optimiser l’homogénéité d’excitation des spins dans chaque coupe mais présente néanmoins certains désavantages par rapport à une acquisition multi-coupe non optimisée (même signaux de RF et de gradient magnétique pour toutes les coupes). [0007] The spatial definition of the section(s) to be imaged, that is to say the position and the inclination of the section, can give rise to a specific optimization of the properties of the RF signals and of the magnetic field gradient signals. static. Noting N the number of cuts considered, the implementation of the global multi-cut sequence therefore requires the definition of N independent sets of signals (1 complex signal per transmission channel at the level of the RF transmission antenna network and 3 real signals for the 3 magnetic field gradient axes). This approach effectively optimizes the homogeneity of spin excitation in each slice but nevertheless has certain disadvantages compared to a non-optimized multi-slice acquisition (same RF and magnetic gradient signals for all slices).
[0008] En premier lieu, l’obtention de ces N ensembles de signaux peut donner lieu à des calculs conséquents en termes de ressources computationnelles. Ceci peut entraver la qualité de l’examen, ou limiter de manière assez arbitraire les possibilités offertes par l’IRM. [0009] Par ailleurs, ces N ensembles de signaux étant obtenus indépendamment, l’analyse conjointe des N images peut dans ces conditions révéler des discontinuités de contraste (en amplitude et en phase) entre coupes. Ce comportement n’est pas souhaitable car il rend difficile ou même impossible une recombinaison de ces N images en une image volumique en vue d’une lecture d’image dans des plans non- alignés au plan de coupe. [0008] Firstly, obtaining these N sets of signals can give rise to substantial calculations in terms of computational resources. This can hamper the quality of the examination, or quite arbitrarily limit the possibilities offered by MRI. [0009] Furthermore, since these N sets of signals are obtained independently, the joint analysis of the N images can under these conditions reveal contrast discontinuities (in amplitude and in phase) between slices. This behavior is undesirable because it makes it difficult or even impossible to recombine these N images into a volume image with a view to image reading in planes not aligned with the cutting plane.
[0010] L’invention vise à apporter une solution à ces deux problèmes. The invention aims to provide a solution to these two problems.
[0011] Conformément à l’invention, ce but est atteint en définissant une paramétrisation régulière (c’est-à-dire continue) des signaux RF et de gradient en fonction des paramètres d’inclinaison et de position de la coupe. Cette paramétrisation étant calculée, il est possible d’implémenter tout protocole multi- coupes en évaluant simplement la paramétrisation sur les N inclinaisons/positions de coupe définies par le protocole, ce qui nécessite un calcul extrêmement simple et quasiment instantanée. On peut parler de concept de « méta-impulsion » ou « paramétrisation globale » en ce sens que la paramétrisation couvre une famille entière d’excitations sélectives optimisées pour créer un angle de bascule uniforme. Cette meta-impulsion est convertie en une impulsion à part entière dès lors que les propriétés de coupe (position et inclinaison) sont fournies. L’approche proposée est entièrement compatible avec les deux techniques d’homogénéisation de l’excitation que sont la transmission parallèle et le contrôle optimal. Elle est en outre compatible avec les « impulsions universelles » divulguées par EP 15306569 ou avec le procédé de conception et d’ajustement d’impulsions radiofréquence universelles en IRM en transmission parallèle divulgué par EP3594710. According to the invention, this object is achieved by defining a regular (that is to say continuous) parametrization of the RF and gradient signals as a function of the inclination and position parameters of the cut. This parametrization being calculated, it is possible to implement any multi-cut protocol by simply evaluating the parametrization on the N inclinations/cut positions defined by the protocol, which requires an extremely simple and almost instantaneous calculation. We can speak of the concept of “meta-pulse” or “global parametrization” in the sense that the parametrization covers an entire family of selective excitations optimized to create a uniform flip angle. This meta-pulse is converted into a full-fledged pulse as soon as the cutting properties (position and inclination) are provided. The proposed approach is fully compatible with the two excitation homogenization techniques that are parallel transmission and optimal control. It is also compatible with the “universal pulses” disclosed by EP 15306569 or with the method for designing and adjusting universal radiofrequency pulses in MRI in parallel transmission disclosed by EP3594710.
[0012] L’invention s’appuie sur la technique des impulsions multi-rayons (en anglais « multi-spokes » ou « Fast-kz ») - connue notamment de [Setsompop 2008] - qui consiste à construire des séquences d’excitations sélectives constituées d’impulsions RF séparées par l’application d’impulsions de gradient magnétique qui permettent de faire interférer les contributions des différentes impulsions RF. En optimisant les moments (aire sous la courbe) des impulsions de gradient magnétique, il est possible de créer ainsi des impulsions composites offrant une bonne homogénéité d’excitation. Typiquement, cette approche est utilisée conjointement au « RF shimming » : chaque impulsion RF est en fait constituée de plusieurs sous-impulsions indépendantes émises par des canaux d’émission (ou transmission, les deux termes seront utilisés de manière indifférenciée) distincts. L’utilisation d’un seul canal d’émission est néanmoins possible, mais cela nécessite l’utilisation de séquences d’impulsions beaucoup plus longues. The invention is based on the technique of multi-ray pulses (in English "multi-spokes" or "Fast-kz") - known in particular from [Setsompop 2008] - which consists in constructing sequences of excitations made up of RF pulses separated by the application of magnetic gradient pulses which allow the contributions of the different RF pulses to interfere. By optimizing the moments (area under the curve) of the magnetic gradient pulses, it is thus possible to create composite pulses offering good excitation homogeneity. Typically, this approach is used in conjunction with “RF shimming”: each RF pulse is in fact made up of several independent sub-pulses emitted by emission channels (or transmission, the two terms will be used interchangeably) distinct. The use of a single transmission channel is nevertheless possible, but this requires the use of much longer pulse sequences.
[0013] Il convient de noter que des techniques de conception globale de séquences multi-coupes ont été proposées par le passé, voir par exemple [Poser 2014], [Guérin 2015], [Gras 2017] Cependant ces techniques visent à résoudre un problème tout à fait différent, à savoir la minimisation du débit d’absorption spécifique et ne permettent pas de remédier aux discontinuités de contraste entre coupes. [0013] It should be noted that techniques for the global design of multi-cut sequences have been proposed in the past, see for example [Poser 2014], [Guérin 2015], [Gras 2017] However, these techniques aim to solve a problem completely different, namely the minimization of the specific absorption rate and do not make it possible to remedy the discontinuities of contrast between sections.
[0014] Un objet de l’invention est un procédé de détermination d’une paramétrisation globale d’une famille de séquences d’impulsions d’excitations de spin en imagerie par résonance magnétique, chaque séquence d’impulsions de ladite famille étant une séquence de type multi-rayon adaptée pour exciter sélectivement des spins nucléaires dans une coupe respective d’un volume d’intérêt d’un corps immergé dans un champ magnétique stationnaire et comprenant un même nombre prédéterminé d’impulsions radiofréquence à une fréquence de Larmor desdits spins nucléaires alternées à des impulsions de gradient magnétique, le procédé comprenant les étapes consistant à : An object of the invention is a method for determining an overall parameterization of a family of sequences of spin excitation pulses in magnetic resonance imaging, each sequence of pulses of said family being a sequence of multi-ray type adapted to selectively excite nuclear spins in a respective section of a volume of interest of a body immersed in a stationary magnetic field and comprising the same predetermined number of radiofrequency pulses at a Larmor frequency of said spins pulses alternated with magnetic gradient pulses, the method comprising the steps of:
A) Acquérir au moins une carte d’intensité d’un rayonnement radiofréquence à la fréquence de Larmor dans ledit volume d’intérêt ; B) Définir un premier ensemble de paramètres caractérisant les impulsions radiofréquence et un second ensemble de paramètres caractérisant les impulsions de gradient magnétique, la valeur de chaque paramètre de chaque dit ensemble étant exprimée par une série tronquée de fonctions d’une coordonnée de position et d’une paire de coordonnées angulaires d’orientation d’une dite coupe, chaque série étant déterminée par ses coefficients ; A) Acquire at least one intensity map of radiofrequency radiation at the Larmor frequency in said volume of interest; B) Define a first set of parameters characterizing the radio frequency pulses and a second set of parameters characterizing the magnetic gradient pulses, the value of each parameter of each said set being expressed by a truncated series of functions of a position coordinate and d a pair of angular orientation coordinates of a said section, each series being determined by its coefficients;
C) A partir de ladite ou desdites cartes d’intensité, calculer les valeurs des coefficients desdites séries tronquées de fonctions qui minimisent une fonction représentative d’un écart moyen à une consigne prédéterminée (par exemple, fournie par un utilisateur) de l’excitation des spins nucléaires, l’écart moyen étant calculé sur ledit volume d’intérêt et sur toutes les orientations et positions possibles desdites coupes ; lesdites valeurs constituant ladite paramétrisation globale. [0015] Un autre objet de l’invention est un procédé de conception d’une séquence d’impulsions d’excitations de spin de type multi-rayon en imagerie par résonance magnétique, ladite séquence d’impulsions étant adaptée pour exciter sélectivement des spins nucléaires dans une coupe respective d’un volume d’intérêt à l’intérieur d’un corps immergé dans un champ magnétique stationnaire et comprenant un même nombre prédéterminé d’impulsions radiofréquence à une fréquence de Larmor desdits spins nucléaires alternées à des impulsions de gradient magnétique, la coupe étant identifiée par des valeurs d’une coordonnée de position et d’une paire de coordonnées angulaires d’orientation, le procédé étant caractérisé en ce qu’il comporte les étapes consistant à : a) Procurer une paramétrisation globale d’une famille de telles séquences d’excitation, ladite paramétrisation globale étant constituée par une première et une deuxième pluralité d’ensembles de coefficients de séries tronquées respectives de fonctions desdites coordonnées de position et d’orientation d’une dite coupe, chaque ensemble de la première pluralité d’ensembles étant constitué des coefficients d’une série tronquée de fonctions définissant un paramètre d’une impulsion radiofréquence et chaque ensemble de la deuxième pluralité d’ensembles étant constitué des coefficients d’une série tronquée de fonctions définissant un paramètre d’une impulsion de gradient magnétique ; et b) Calculer les valeurs desdits paramètres d’impulsions de gradient et d’impulsions radiofréquence à partir de ladite paramétrisation globale, desdites séries tronquées de fonctions et des valeurs de ladite coordonnée de position de ladite paire de coordonnées angulaires d’orientation de la coupe. C) From said intensity map or maps, calculate the values of the coefficients of said truncated series of functions which minimize a function representative of an average deviation from a predetermined setpoint (for example, supplied by a user) of the excitation nuclear spins, the mean deviation being calculated on said volume of interest and on all possible orientations and positions of said slices; said values constituting said global parametrization. Another object of the invention is a method for designing a sequence of spin excitation pulses of the multi-ray type in magnetic resonance imaging, said sequence of pulses being adapted to selectively excite spins nuclear spins in a respective section of a volume of interest inside a body immersed in a stationary magnetic field and comprising the same predetermined number of radio frequency pulses at a Larmor frequency of said nuclear spins alternating with gradient pulses magnetic, the section being identified by values of a position coordinate and of a pair of angular coordinates of orientation, the method being characterized in that it comprises the steps consisting in: a) Procuring a global parametrization of a family of such excitation sequences, said global parametrization being constituted by a first and a second plurality of sets of coefficients of respective truncated series of functions of said coordinates of position and orientation of a said cut, each set of the first plurality of sets consisting of the coefficients of a truncated series of functions defining a parameter of a radio frequency pulse and each set of the second plurality sets consisting of the coefficients of a truncated series of functions defining a parameter of a magnetic gradient pulse; and b) Calculating the values of said gradient pulse and radio frequency pulse parameters from said global parameterization, said truncated series of functions and the values of said position coordinate of said pair of angular coordinates of orientation of the cut .
[0016] La paramétrisation globale peut notamment être obtenue par un procédé tel qu’esquissé ci-dessus. [0016] The global parametrization can in particular be obtained by a method as outlined above.
[0017] Encore un autre objet de l’invention est un procédé d’imagerie par résonance magnétique comprenant les étapes suivantes : i) Concevoir une pluralité de séquences d’impulsions d’excitations de spin de type multi-rayon, chaque dite séquence d’impulsions étant adaptée à exciter sélectivement des spins nucléaires dans une coupe respective d’un volume d’intérêt d’un corps à imager immergé dans un champ magnétique stationnaire et comprenant des impulsions radiofréquence à une fréquence de Larmor desdits spins nucléaires alternées à des impulsions de gradient magnétique, chaque coupe étant identifiée par des valeurs respectives d’une coordonnée de position et d’une paire de coordonnées angulaires d’orientation ; ii) Appliquer audit corps, en succession, lesdites séquences d’impulsions d’excitations de spin, un gradient magnétique de sélection de coupe étant appliqué en même temps que chaque impulsion radiofréquence ; iii) Après chaque dite séquence d’impulsions d’excitation de spin, acquérir un signal de résonance magnétique émis par les spins nucléaires ; et iv) Traiter les signaux ainsi acquis pour reconstruire une image d’une portion dudit volume de référence définie par l’union desdites coupes ; caractérisé en ce que l’étape i) est mise en œuvre par un procédé de conception d’une séquence d’impulsions d’excitations de spin tel que mentionné ci-dessus. [0017] Yet another object of the invention is a magnetic resonance imaging method comprising the following steps: i) designing a plurality of sequences of multi-ray type spin excitation pulses, each said sequence of pulses being adapted to selectively excite nuclear spins in a respective section of a volume of interest of a body to be imaged immersed in a stationary magnetic field and comprising radiofrequency pulses at a Larmor frequency of said nuclear spins alternated with pulses of magnetic gradient, each cut being identified by respective values of a position coordinate and a pair of angular orientation coordinates; ii) Applying to said body, in succession, said sequences of spin excitation pulses, a magnetic cut selection gradient being applied at the same time as each radio frequency pulse; iii) After each said sequence of spin excitation pulses, acquiring a magnetic resonance signal emitted by the nuclear spins; and iv) processing the signals thus acquired to reconstruct an image of a portion of said reference volume defined by the union of said sections; characterized in that step i) is implemented by a method of designing a sequence of spin excitation pulses as mentioned above.
[0018] Encore un autre objet de l’invention est un ordinateur programmé pour : a) Recevoir en entrée au moins une carte d’intensité d’un rayonnement radiofréquence dans un volume d’intérêt d’un corps immergé dans un champ magnétique stationnaire, ledit rayonnement radiofréquence étant à une fréquence de Larmor de spins nucléaires dudit corps ; b) Définir une paramétrisation globale d’une famille de séquences d’impulsions d’excitations de spin de type multi-rayon en imagerie par résonance magnétique, chaque séquence d’impulsions de ladite famille étant adaptée pour exciter sélectivement des spins nucléaires dans une coupe dudit volume d’intérêt et comprenant des impulsions radiofréquence à une fréquence de Larmor desdits spins nucléaires alternées à des impulsions de gradient magnétique, ladite paramétrisation globale comprenant un premier ensemble de paramètres caractérisant lesdites impulsions radiofréquence et un second ensemble de paramètres caractérisant lesdites impulsions de gradient magnétique, la valeur de chaque paramètre de chaque dit ensemble étant exprimée par une série tronquée de fonctions d’une coordonnée de position et d’une paire de coordonnées angulaires d’orientation d’une dite coupe, chaque série étant déterminée par ses coefficients ; g) A partir de ladite ou desdites cartes d’intensité, calculer les valeurs des coefficients desdites séries tronquées de fonctions qui minimisent une fonction représentative d’un écart moyen à une consigne de l’excitation des spins nucléaires, la moyenne étant calculée sur ledit volume d’intérêt et sur toutes les orientations et positions possibles desdites coupes. [0019] Encore un autre objet de l’invention est un ordinateur programmé pour la mise en oeuvre d’un procédé de conception d’une séquence d’impulsions d’excitations de spin tel qu’esquissé plus haut. Yet another object of the invention is a computer programmed to: a) receive as input at least one intensity map of radiofrequency radiation in a volume of interest of a body immersed in a stationary magnetic field , said radiofrequency radiation being at a Larmor frequency of nuclear spins of said body; b) Define a global parametrization of a family of sequences of spin excitation pulses of the multi-ray type in magnetic resonance imaging, each sequence of pulses of said family being adapted to selectively excite nuclear spins in a section of said volume of interest and comprising radiofrequency pulses at a Larmor frequency of said nuclear spins alternating with magnetic gradient pulses, said global parameterization comprising a first set of parameters characterizing said radiofrequency pulses and a second set of parameters characterizing said gradient pulses magnetic, the value of each parameter of each said set being expressed by a truncated series of functions of a position coordinate and a pair of angular orientation coordinates of a said section, each series being determined by its coefficients; g) From said intensity map or maps, calculate the values of the coefficients of said truncated series of functions which minimize a function representative of an average deviation from a setpoint of the excitation of the nuclear spins, the average being calculated on said volume of interest and on all possible orientations and positions of said sections. Yet another object of the invention is a computer programmed for implementing a method for designing a sequence of spin excitation pulses as outlined above.
[0020] Un ultérieur objet de l’invention est un appareil d’imagerie par résonance magnétique équipé d’un tel ordinateur. A further object of the invention is a magnetic resonance imaging device equipped with such a computer.
[0021] D’autres caractéristiques, détails et avantages de l’invention ressortiront à la lecture de la description faite en référence aux dessins annexés donnés à titre d’exemple et qui représentent, respectivement : Other characteristics, details and advantages of the invention will become apparent on reading the description made with reference to the appended drawings given by way of example and which represent, respectively:
[0022] [Fig.1] un appareil, ou « scanner » d’IRM convenant à la mise en œuvre de l’invention ; [0022] [Fig.1] an MRI device or “scanner” suitable for implementing the invention;
[0023] [Fig. 2], les coordonnées permettant d’identifier une coupe d’un corps à imager. [0023] [Fig. 2], the coordinates making it possible to identify a section of a body to be imaged.
[0024] [Fig. 3], une portion d’une séquence d’excitation selon un mode de réalisation de l’invention ; [0024] [Fig. 3], a portion of an excitation sequence according to one embodiment of the invention;
[0025] [Fig. 4], [Fig. 5a], [Fig. 5b], [Fig. 6], [Fig. 7a], [Fig. 7b] et [Fig. 7c], des graphiques illustrant l’effet technique de l’invention. [0025] [Fig. 4], [Fig. 5a], [Fig. 5b], [Fig. 6], [Fig. 7a], [Fig. 7b] and [Fig. 7c], graphics illustrating the technical effect of the invention.
[0026] [Fig. 8], un organigramme d’un procédé selon un mode de réalisation de l’invention. [0026] [Fig. 8], a flowchart of a method according to one embodiment of the invention.
[0027] La [Fig. 1] illustre très schématiquement un appareil (ou « scanner ») d’IRM SI pouvant convenir à la mise en œuvre de l’invention. Cet appareil comprend : un aimant AL pour générer un champ magnétique stationnaire B0 d’aimantation, dit champ magnétique longitudinale, orienté selon une direction « z » ; des bobines BGx (non représentée), BGy et BGz pour générer des gradients de champ magnétique selon les directions mutuellement orthogonales « x », « y » et « z » ; une pluralité d’antennes d’émission et réception à radiofréquence CTx/CRx entourant une région ou volume d’intérêt ROI pouvant contenir un corps C à imager (on entend par « corps » tout objet matériel, biologique ou pas, pouvant être imagé par IRM ; il s’agira typiquement d’un corps humain ou d’une partie de corps humain, par exemple une tête) ; en particulier, la région d’intérêt contient l’isocentre IC (référence visible sur la [Fig. 2]), qui est le seul point où le champ magnétique reste constant quel que soit le gradient appliqué par les bobines BGx, BGy, BGz ; et un ordinateur de contrôle OC qui pilote les bobines de gradient BGx, BGy et BGz et les antennes RF CTx/CRx. [0027] The [Fig. 1] very schematically illustrates an IS MRI device (or “scanner”) that may be suitable for implementing the invention. This device comprises: a magnet AL for generating a stationary magnetic field B 0 of magnetization, called longitudinal magnetic field, oriented in a “z” direction; coils BGx (not shown), BGy and BGz for generating magnetic field gradients along mutually orthogonal “x”, “y” and “z” directions; a plurality of radiofrequency transmission and reception antennas CTx/CRx surrounding a region or volume of interest ROI which can contain a body C to be imaged (the term “body” is understood to mean any material object, biological or not, which can be imaged by MRI; this will typically be a human body or part of a human body, for example a head); in particular, the region of interest contains the IC isocenter (reference visible in [Fig. 2]), which is the only point where the magnetic field remains constant regardless of the gradient applied by the coils BGx, BGy, BGz; and a control computer OC which drives the gradient coils BGx, BGy and BGz and the RF antennas CTx/CRx.
[0028] Par l’intermédiaire de circuits électroniques non représentés, l’ordinateur OC pilote les bobines de gradient BGx, BGy et BGz pour qu’elles génèrent des impulsions de gradient et les antennes RF CTx/CRx pour qu’elles génèrent des impulsions radiofréquence de manière à appliquer au corps C une séquence d’excitation obtenue par le procédé de l’invention. Par ailleurs, l’ordinateur OC reçoit des antennes CTx/CRx des signaux de résonance magnétique émis par les spins nucléaires de l’objet suite à leur excitation, et traite ces signaux pour reconstituer - de manière conventionnelle - une image d’IRM du corps C. [0028] Via electronic circuits not shown, the computer OC drives the gradient coils BGx, BGy and BGz so that they generate gradient pulses and the RF antennas CTx/CRx so that they generate pulses radiofrequency so as to apply to the body C an excitation sequence obtained by the method of the invention. Furthermore, the OC computer receives from the CTx/CRx antennas magnetic resonance signals emitted by the nuclear spins of the object following their excitation, and processes these signals to reconstitute - in a conventional manner - an MRI image of the body vs.
[0029] Des séquences d’excitation appropriées permettent par ailleurs à l’ordinateur OC d’obtenir des cartes de la distribution spatiale du champ RF B émis par chaque antenne. Ces cartes sont nécessaires à la mise en œuvre de l’invention, comme cela sera expliqué en détail plus loin. Le document [Amadon 2012] décrit un procédé d’acquisition de cartes de B . [0029] Appropriate excitation sequences also allow the computer OC to obtain maps of the spatial distribution of the RF field B emitted by each antenna. These cards are necessary for the implementation of the invention, as will be explained in detail later. The document [Amadon 2012] describes a process for acquiring maps of B .
[0030] L’appareil de la figure 1 comprend une pluralité d’antennes RF permettant d’avoir recours à des techniques de transmission parallèle, mais l’utilisation d’une seule antenne est également possible, bien que moins avantageuse. Par ailleurs, les antennes émettrices et réceptrices peuvent être distinctes. The apparatus of FIG. 1 comprises a plurality of RF antennas making it possible to use parallel transmission techniques, but the use of a single antenna is also possible, although less advantageous. Furthermore, the transmitting and receiving antennas may be distinct.
[0031] La [Fig. 2] permet de définir les coordonnées qui identifient une coupe CP du corps C. [0031] The [Fig. 2] allows to define the coordinates that identify a CP section of the body C.
[0032] Soit n le vecteur unitaire perpendiculaire à la coupe. L’orientation de ce vecteur est définie par les coordonnées angulaires sphériques Q (azimut ou longitude) et f (colatitude) ayant pour origine l’isocentre IC. Soit r le vecteur position d’un point quelconque de la coupe par rapport à l’isocentre ; le produit scalaire n r donne z, la distance de la coupe à l’origine. Le scalaire z peut être considéré une coordonnée de position de la coupe. Let n be the unit vector perpendicular to the cut. The orientation of this vector is defined by the spherical angular coordinates Q (azimuth or longitude) and f (colatitude) having for origin the isocenter IC. Let r be the position vector of any point of the section with respect to the isocenter; the scalar product n r gives z, the distance from the cut to the origin. The scalar z can be considered a position coordinate of the cut.
[0033] La [Fig. 3] est un schéma d’une simple séquence d’excitation selon un mode de réalisation de l’invention, comprenant trois impulsions radiofréquence IRF1 , IRF2 et IRF3 et trois impulsions de gradient IG1 , IG2 et IG3 alternés avec les impulsions radiofréquence. Des gradients de sélection de coupe GS sont appliqués en même temps que les impulsions radiofréquence. [0033] The [Fig. 3] is a diagram of a simple excitation sequence according to one embodiment of the invention, comprising three radiofrequency pulses IRF1, IRF2 and IRF3 and three gradient pulses IG1, IG2 and IG3 alternated with the pulses radio frequency. GS cut selection gradients are applied at the same time as the radio frequency pulses.
[0034] Dans un souci de simplicité, la [Fig. 3] montre uniquement les impulsions radiofréquence associées à un seul canal d’émission, identifié par un indice « i ».[0034] For the sake of simplicity, [Fig. 3] shows only the radiofrequency pulses associated with a single emission channel, identified by an index “i”.
Les impulsions radiofréquence IRF1, IRF2 et IRF3 présentent une même forme temporelle (un sinus cardinal) et des amplitudes complexes différentes Les parties réelles et imaginaires des impulsions - d’amplitude Re(aj,j) et lm(aj,j), - sont représentées en trait continu et haché, respectivement. Pour exciter une coupe décalée par rapport à l’isocentre ( z ¹ 0) on multiplie les formes d’onde des impulsions radiofréquence par une porteuse exp (ίD/ 1) avec D/ = étant la résultant des gradients de sélection appliqués par les différentes bobines, g le rapport gyromagnétique de noyaux et i l’unité imaginaire. The IRF1, IRF2 and IRF3 radiofrequency pulses have the same temporal shape (a cardinal sine) and different complex amplitudes The real and imaginary parts of the pulses - of amplitude Re(aj, j ) and lm(aj, j ), - are represented in continuous and dashed lines, respectively. For excite a section shifted with respect to the isocenter ( z ¹ 0) the waveforms of the radio frequency pulses are multiplied by a carrier exp (ίD/ 1) with D/ = being the resultant of the selection gradients applied by the different coils, g the gyromagnetic ratio of nuclei and i the imaginary unit.
[0035] Les gradients de sélection GS sont appliqués en même temps que les impulsions radiofréquence. Seule une composante, identifiée par un indice i' e {x,y,z}, et d’amplitude maximale Gmaxj’ est représentée. L’amplitude de la résultante des trois composantes détermine l’épaisseur de la coupe. On remarque que les impulsions de gradient présentent une polarité alternée ; d’une manière connue en soi, cela permet de se passer d’impulsions de gradient de « rembobinage » de la phase qui seraient nécessaire si toutes les impulsions de gradient avaient la même polarité ([Gras 2017 B]). En principe, lorsqu’on change la polarité des impulsions de gradient, il est nécessaire de retourner temporellement les impulsions radiofréquence, mais dans le cas considéré ici ces dernières sont symétriques. The selection gradients GS are applied at the same time as the radiofrequency pulses. Only one component, identified by an index i' e {x,y,z}, and of maximum amplitude G maxj' is represented. The amplitude of the resultant of the three components determines the thickness of the cut. Note that the gradient pulses have an alternating polarity; in a known way, this makes it possible to dispense with phase "rewind" gradient pulses which would be necessary if all the gradient pulses had the same polarity ([Gras 2017 B]). In principle, when the polarity of the gradient pulses is changed, it is necessary to time-reverse the radiofrequency pulses, but in the case considered here the latter are symmetrical.
[0036] Les impulsions de gradient IG1 , IG2 et IG3, qui ne se superposent pas aux impulsions radiofréquence, déterminent des positions respectives dans l’espace k à partir desquelles les gradients de sélection « tracent » des rayons (« spokes »). La position dans l’espace k déterminée par une impulsion de gradient dépend uniquement de l’aire Br de chacune de ses trois composantes. Il convient de noter que les impulsions de gradient sont beaucoup moins intenses que les gradients de sélection ; dans la [Fig. 3], dans un souci de visibilité, elles ont été amplifiées d’un facteur 100. The gradient pulses IG1, IG2 and IG3, which do not overlap the radio frequency pulses, determine respective positions in k space from which the selection gradients “draw” rays (“spokes”). The position in k-space determined by a gradient pulse depends only on the area Br of each of its three components. It should be noted that gradient pulses are much less intense than selection gradients; in [Fig. 3], for the sake of visibility, they have been amplified by a factor of 100.
[0037] Un procédé de détermination d’une paramétrisation globale d’une famille de séquences d’impulsions d’excitation de spin selon un mode de réalisation de l’invention sera décrit en détail dans la suite. La description repose dans l’approximation des petits angles de bascule des spins, qui permet de linéariser les équations de Bloch, mais une généralisation aux grands angles est tout à fait possible. A method for determining an overall parameterization of a family of spin excitation pulse sequences according to one embodiment of the invention will be described in detail below. The description rests in the approximation of the small angles of flipping of the spins, which makes it possible to linearize the Bloch equations, but a generalization to the large angles is quite possible.
[0038] Soit Nc le nombre de canaux et Ns le nombre d’impulsions radiofréquence utilisées. L’information sur le moment de chaque impulsion RF et de chaque impulsion de gradient) peut être représentée sous forme matricielle, à savoir une matrice complexe A e (CNcXNs et une matrice B e ®>3 x(Ns-i) respectivement. Dans cette représentation, chaque colonne indique les moments par canal de transmission (par axe - x, y, z) de l’impulsion RF (de l’impulsion de gradient). Let N c be the number of channels and N s the number of radio frequency pulses used. The information on the timing of each RF pulse and each gradient pulse) can be represented in matrix form, namely a complex matrix A e (C NcXNs and a matrix B e ® > 3 x ( N si) respectively. In this representation, each column indicates the moments per transmission channel (by axis - x, y, z) of the RF pulse (of the gradient pulse).
[0039] Pour satisfaire la condition selon laquelle le moment total du gradient magnétique le long de chaque axe est nul, la dernière sous-impulsion RF est suivie d’une impulsion de gradient dont les composantes suivant x, y, et z est égale à : To satisfy the condition that the total moment of the magnetic gradient along each axis is zero, the last RF sub-pulse is followed by a gradient pulse whose components along x, y, and z is equal to :
[0043] Le fait que cette condition soit satisfaite n’est pas essentiel (l’angle de bascule n’en dépend pas) mais peut s’avérer nécessaire suivant la séquence IRM (séquences de type « steady-state free precession »). [0043] The fact that this condition is satisfied is not essential (the tilt angle does not depend on it) but may prove necessary depending on the MRI sequence (“steady-state free precession” type sequences).
[0044] On note B le vecteur des coefficients de gradient magnétiques complété de l’impulsion finale ; chaque élément du vecteur est égal à la somme des aires des impulsions de gradient le long d’un axe. We note B the vector of the gradient coefficients magnets completed with the final impulse; each element of the vector is equal to the sum of the areas of the gradient pulses along an axis.
[0045] On note B^ (r) e (ClxNc les valeurs du champ magnétique résonant en r sur chaque canal de transmission. Ces valeurs sont normalisées à 1 Volt transmis sur chaque canal ; elles sont donc exprimées en Tesla par Volt (T/V). We note B^ (r) e (C lxNc the values of the resonant magnetic field in r on each transmission channel. These values are normalized to 1 Volt transmitted on each channel; they are therefore expressed in Tesla per Volt (T /V).
[0046] On définit la position des impulsions RF dans l’espace K de transmission par : la matrice triangulaire inférieure : The position of the RF pulses in the transmission space K is defined by: the lower triangular matrix:
[0049] T [0050] En l’absence de gradient de sélection de coupe, l’angle de bascule des spins nucléaires a(r) créé en r par la succession des Ns impulsions RF s’écrit, dans l’approximation des petits angles de bascule : [0049] T In the absence of a cut selection gradient, the flip angle of the nuclear spins a(r) created at r by the succession of N s RF pulses is written, in the small flip angle approximation :
[0051 ] [0051 ]
[0052] où A est la matrice de terme général: [0052] where A is the general term matrix:
[0053] A [0053] A
[0054] où dénote le coefficent de la ligne i’ et de la colonne j dans la matrice K, et [0054] where denotes the coefficient of row i' and column j in the matrix K, and
[0055] où [0055] where
[0056] On note FAr l’opérateur ( qui dépend de la distribution du champ radiofréquence dans le corps à imager, c’est-à-dire d’une ou plusieurs (en cas de transmission parallèle) cartes B We note FA r the operator ( which depends on the distribution of the field radio frequency in the body to be imaged, i.e. one or more (in case of parallel transmission) B cards
[0057] En présence d’un gradient de sélection de coupe, en tout point de l’axe de la coupe C l’angle de bascule est encore donné par In the presence of a cut selection gradient, at any point on the axis of the cut C the tilt angle is again given by
[0058] Pour être entièrement défini, un protocole multi-coupes à N coupes requiert la définition de N paramétrisations distinctes P To be completely defined, a multi-cut protocol with N cuts requires the definition of N distinct parametrizations P
[0059] Pour chaque coupe C chaque coefficient complexe de la matrice A(c) et chaque coefficient réel b de la matrice B(c) s’écrit comme une fonction des paramètres géométriques de la coupe [0059] For each cut C each complex coefficient of the matrix A (c) and each real coefficient b of the matrix B (c) is written as a function of the geometric parameters of the cut
[0060] [0060]
[0061] [0061]
[0062] où et sont des fonctions continues de n et de respectivement à valeurs complexes et réelles (les impulsions de gradient étant à valeur réelle). [0062] where and are continuous functions of n and of respectively at complex and real values (the gradient pulses being at real value).
[0063] Les coefficients complexes et les coefficients réels constituent deux ensembles de paramètres qui, conjointement, déterminent une séquence d’excitation. [0064] En remarquant que les fonctions doivent de plus satisfaire la condition de symétrie: [0063] The complex coefficients and the real coefficients constitute two sets of parameters that together determine an excitation sequence. [0064] By noting that the functions must no longer satisfy the symmetry condition:
[0065] [0065]
[0067] Bien que d’autre paramétrisations, soient possible, il peut être judicieux de prendre des fonctions exprimées par des séries tronquées de la forme : [0067] Although other parametrizations are possible, it may be wise to take functions expressed by truncated series of the form:
[0070] Où est une longueur suffisamment grande pour que la boule centrée sur l’isocentre et de rayon z0 contienne la région d’intérêt, [0070] Where is a length large enough for the ball to be centered on the isocenter and of radius z 0 contains the region of interest,
[0072] et sont des fonctions continues sur la sphère unité satisfaisant la condition: [0072] and are continuous functions on the sphere unit satisfying the condition:
[0074] et [0074] and
[0075] k [0075] k
[0076] Les harmoniques sphériques formant une base de fonctions continues sur la sphère unité, il convient de choisir pour fi j k des sommes finies (jusqu’à l’ordre q) d’harmoniques sphériques complexes en ne retenant pour k pair que les harmoniques dont le degré est pair et pour k impair, que les harmoniques dont le degré est impair: The spherical harmonics forming a basis of continuous functions on the unit sphere, it is advisable to choose for f ijk finite sums (up to the order q) of complex spherical harmonics by retaining for k even only the harmonics whose degree is even and for k odd, than the harmonics whose degree is odd:
[0078] Avec la contrainte: [0078] With the constraint:
[0079] [0080] et : [0079] [0080] and:
[0082] où P™ est le polynôme de Legendre d’ordre (1, m). Avec la contrainte ci- dessous, on peut vérifier aisément que la contrainte de symétrie des fonction est bien satisaite pout tout i et pour tout j. [0082] where P™ is the Legendre polynomial of order (1, m). With the constraint below, we can easily verify that the symmetry constraint of the functions is well satisfied for all i and for all j.
[0083] Pour gk, la même règle s’applique à différence que la décomposition se fait sur les harmoniques sphériques réelles, jusqu’à un certain ordre q', et que les coefficients gi'j',kM,m sont réels : For g k , the same rule applies except that the decomposition is done on the real spherical harmonics, up to a certain order q', and that the coefficients gi' j ' ,k M ,m are real :
[0091] En conclusion, chaque terme de chaque série tronquée exprimant fy et g^y est donné par le produit d’un monôme de la coordonnée de position z et d’une somme finie d’harmoniques sphériques complexes ou réelles (plus généralement, de fonctions appartenant à une base de fonctions continues des coordonnées angulaires d’orientation). In conclusion, each term of each truncated series expressing fy and g^y is given by the product of a monomial of the position coordinate z and a finite sum of complex or real spherical harmonics (more generally, of functions belonging to a basis of continuous functions of the angular coordinates of orientation).
[0092] Ainsi, dans cette approche, tout protocole multi-coupe est entièrement caractérisé dès lors que les matrices : Thus, in this approach, any multi-cut protocol is fully characterized as soon as the matrices:
[0093] F = (f4Um) [0093] F = (f 4Um )
[0094] et : [0096] sont définies. Soit X le vecteur réel composé des coefficients de F suivis de ceux de G. Notons par et par n^( ) les matrices A et B obtenues par évaluation des fonctions Notons la paramétrisation complète ainsi obtenue, suffisante pour caractériser l’action de la séquence d’impulsions dans l’axe de la coupe [0094] and: [0096] are defined. Let X be the real vector composed of the coefficients of F followed by those of G. Let us denote by and by n ^( ) the matrices A and B obtained by job evaluation Note it complete parametrization thus obtained, sufficient to characterize the action of the pulse sequence in the axis of the cut
[0097] Se pose ensuite le problème de l’optimisation des coefficients de ces matrices de manière à obtenir la meilleure conformité possible de l’angle de bascule à une consigne donnée d’angle de bascule aT, quelle que soit la coupe considérée. The problem then arises of optimizing the coefficients of these matrices so as to obtain the best possible compliance of the rocking angle with a given rocking angle set point a T , whatever the cut considered.
[0098] Pour une position donnée et pour un vecteur unitaire n donné, soit Cn(r) la coupe normale à n passant par r. Sa position étant donnée par : For a given position and for a given unit vector n, let C n (r) the cut normal to n through r. Its position being given by:
[0099] [0099]
[0100 n’est autre que C(n,n · r). Etant donné une paramétrisation X, on peut définir l’angle de bascule an r(X) obtenu en r en « jouant » la séquence c’est-à-dire en appliquant cette séquence au corps à imager. [0100 is nothing but C(n,n · r). Given a parametrization X, we can define the tilt angle a nr (X) obtained in r by "playing" the sequence that is, by applying this sequence to the body to be imaged.
[0101] Le point r étant un point dans l’axe de la coupe Cn(r) on obtient : The point r being a point in the axis of the cut C n (r) we obtain:
[0102] ) [0102] )
[0103] On peut alors définir en chaque point l’écart quadratique moyen à la consigne d’angle de bascule lorsque l’axe de sélection de coupe parcourt la sphère unitaire : [0103] It is then possible to define at each point the mean square deviation from the rocking angle setpoint when the cutting selection axis traverses the unit sphere:
[0104] [0104]
[0105] Par sommation dans la région d’intérêt ROI (désignée dans les équations) on peut définir comme suit l’écart quadratique moyen à la consigne : [0105] By summation in the ROI region of interest (designated in the equations) the root mean square deviation from the setpoint can be defined as follows:
[0106] [0106]
[0107] L’obtention d’une paramétrisation optimale X peut être finalement obtenue par la résolution du problème d’optimisation : [0107] Obtaining an optimal parametrization X can finally be obtained by solving the optimization problem:
[0108] Gh1hc [0108] Gh1h c
[0109] Avantageusement, des contraintes sont adjointes à ce problème pour obtenir des solutions admissibles. Ces contraintes portent notamment sur l’énergie des impulsions RF. [0110] On note s: [0,1] ® IR la forme temporelle d’une générique impulsion RF et T sa durée. Soit M1(s) = les moments d’ordre 1 et 2 de l’impulsion s(t). Pour un choix de paramétrisation on peut définir (on suppose ici une adaptation d’impédance à 50 Ohm de tous la chaîne de transmission RF): 2 [0109] Advantageously, constraints are added to this problem to obtain acceptable solutions. These constraints relate in particular to the energy of the RF pulses. We denote s: [0,1] ® IR the temporal form of a generic RF pulse and T its duration. Let M 1 (s) = the moments of order 1 and 2 of the pulse s(t). For a choice of parametrization we can define (we here assumes 50 Ohm impedance matching of the entire RF transmission chain): 2
[0111] L’énergie par canal et par sous-impulsion [0111] Energy per channel and per sub-pulse
[0112] L’énergie par canal [0112] Energy per channel
[0113] L’énergie totale du puise composite [0113] The total energy of the composite pulse
[0114] Etant donné une paramétrisation X, on peut définir pout tout i et pout tout j: [0114] Given a parametrization X, we can define for all i and for all j:
[0116] Il convient de noter que le fait que z soit borné permet d’assurer que P et donc , soit borné (alors même que les fonctions polynomiales à la base de la construction de ne le sont pas). Cet ensemble de fonctions permet de définir un ensemble de contraintes restreignant l’espace des solutions à un domaine donnant lieu à des paramétrisations admissibles en termes de bilan d’énergie radiofréquence. On peut ainsi définir : [0116] It should be noted that the fact that z is bounded ensures that P and therefore , is bounded (even though the polynomial functions at the base of the construction of are not). This set of functions makes it possible to define a set of constraints restricting the space of solutions to a domain giving rise to admissible parametrizations in terms of radiofrequency energy balance. We can thus define:
[0118] L’ajout de contraintes explicites du type : [0118] The addition of explicit constraints of the type:
[0120] permet par exemple limiter l’énergie de chaque sous-impulsion pour l’orientation de coupe et position de coupe la plus défavorable parmi toutes celles envisageables ayant une intersection avec la région d’intérêt [0120] makes it possible, for example, to limit the energy of each sub-pulse for the most unfavorable cutting orientation and cutting position among all those that can be envisaged having an intersection with the region of interest
[0121] L’effet technique de l’invention a été démontré dans un cas simple en utilisant des données (carte de champ radiofréquence d’émission) acquise sur le système Siemens Magnetom 7T équipé d’un système de transmission parallèle à 8 canaux d’un réseau d’antenne tête à 8 canaux de transmission et 32 canaux de réception de Nova Medical. Cette preuve de concept est limitée à un axe de coupe dans le plan YZ, autrement dit à f = 0. Ceci a pour conséquence de réduire de 1 la dimension des fonction f et g. Ainsi, dans cette démonstration, les fonctions apparaissent comme des fonctions de z et de Q et les fonctions fi j k comme des fonctions continues sur le cercle unité. Elles admettent la décomposition: [0121] The technical effect of the invention has been demonstrated in a simple case using data (emission radio frequency field map) acquired on the Siemens Magnetom 7T system equipped with an 8-channel parallel transmission system of an 8 transmit and 32 receive channel head coil array from Nova Medical. This proof of concept is limited to a cutting axis in the YZ plane, in other words at f = 0. This has the consequence of reducing the dimension by 1 functions f and g. Thus, in this demonstration, the functions appear as functions of z and Q and the functions f ijk as continuous functions on the unit circle. They admit the decomposition:
[0123] Le même raisonnement s’applique sur g. Cette fonction étant réelle, on a : The same reasoning applies to g. This function being real, we have:
[0125] Avec réels pout tout 1. On se limite au cas Ns = 2 (séquence composées de deux impulsions) et on prend 8 cm. Pour la sous-impulsion, T = 1 ms, [0125] With reals for all 1. We limit ourselves to the case N s = 2 (sequence composed of two pulses) and we take 8 cm. For the under-pulse, T = 1 ms,
[0126] [0126]
[0127] Où représente la fonction d’apodisation de Planning. On a pour ce choix de sous-impulsion : 2.2894. Where represents the apodization function of Planning. We have for this choice of subpulse: 2.2894.
[0129] La démonstration s’appuie sur la carte de B (5 mm de résolution isotrope) mesurée avec la séquence XFL sur la tête d’un sujet humain ; le problème posé consiste à calculer la paramétrisation X minimisant l’écart quadratique moyen à un angle de bascule cible de 30° sur la région d’intérêt composée de l’ensemble des voxels du cerveau (métrique e). On ajoute à l’objective la contrainte : The demonstration is based on the map of B (5 mm of isotropic resolution) measured with the XFL sequence on the head of a human subject; the problem posed consists in calculating the parametrization X minimizing the root mean square deviation at a target flip angle of 30° on the region of interest composed of all the voxels of the brain (metric e). We add to the objective the constraint:
[0131] Avec : [0131] With:
[0132] Esp = 33 mJ [0132] E sp = 33 mJ
[0133] Le calcul de la métrique e(C), impliquant une intégration sur cercle unité, cette intégration est ramenée à une intégration sur le demi-cercle en vertu du fait que et est discrétisée de manière assez grossière en posant comme nouvelle fonction de coût à minimiser : ). [0136] Pour le calcul de est évalué ensemble de valeurs discrètes de z régulièrement distribué dans l’intervalle [-<o, <o] avec un pas de 5 mm. The calculation of the metric e(C), implying an integration on a unit circle, this integration is reduced to an integration on the semi-circle by virtue of the fact that and is discretized in a rather crude way by setting as a new cost function to be minimized : ). [0136] For the calculation of is evaluated set of discrete values of z regularly distributed in the interval [-<o, <o] with a step of 5 mm.
[0137] Pour la minimisation sous contrainte de on utilise l’implémentation Matlab (marque déposée) (R2019a Update 7) de l’algorithme SQP (Sequential Quadratic Programming). For the constrained minimization of we use the implementation Matlab (registered trademark) (R2019a Update 7) of the SQP (Sequential Quadratic Programming) algorithm.
[0138] Pour l’initialisation de l’optimisation, on adopte l’approche en deux étapes détaillée ci-dessous : For the initialization of the optimization, we adopt the two-step approach detailed below:
[0139] i) (une seule sous-impulsion et par conséquent pas de « saut » dans l’espace k de transmission puisque ce dernier suppose [0139] i) (a single sub-pulse and therefore no "jump" in the transmission k-space since the latter supposes
[0140] ii) pour tout est indépendant de Q et de i et est défini par l’équation: [0140] ii) for all is independent of Q and i and is defined by the equation:
[0141 ] [0141 ]
[0142] Où : [0142] Where:
[0144] Où est le mode de transmission standard de l’antenne. Where is the standard transmission mode of the antenna.
[0145] Partant de cette initialisation, une première optimisation SQP avec et est appliquée en limitant le nombre d’itérations à 50. Cette optimisation recherche une solution dans un espace de dimension [0145] Starting from this initialization, a first SQP optimization with and is applied by limiting the number of iterations to 50. This optimization searches for a solution in a space of dimension
[0146] La solution obtenue ) après convergence de l’algorithme est utilisée pour définir une nouvelle initialisation à deux impulsions, toutes deux en k = 0 (g = 0) et toutes deux d’amplitude f sur tous les canaux. On vérifie aisément que le profil d’angle de bascule n’est pas modifié par cette opération. Cette deuxième itération de l’algorithme SQP est appliquée avec jusqu’à la convergence de l’algorithme. Cette optimisation relève d’un problème de dimension 396 (X de dimension 396). [0146] The solution obtained after convergence of the algorithm is used to define a new initialization with two pulses, both in k = 0 (g = 0) and both of amplitude f on all channels. It is easily verified that the tilt angle profile is not modified by this operation. This second iteration of the SQP algorithm is applied with up to the convergence of the algorithm. This optimization is a 396 dimension problem (X dimension 396).
[0147] Afin de montrer un intérêt pratique de cette approche paramétrique pour définir toute excitation sélective optimisée vis-à-vis de l’homogénéité de l’angle de bascule dans la coupe, la méthode avec un design classique par coupe. On considère le cas de l’excitation de 52 coupes axiales (Q = 0) d’épaisseur 1 mm et respectant une distance de 1 ,2 mm (centre à centre) entre deux coupes adjacentes. Pour chaque position de coupe, on calcule la meilleure paramétrisation n s no suivant le critère : [0149] Où In order to show the practical interest of this parametric approach for defining any selective excitation optimized with respect to the homogeneity of the tilt angle in the slice, the method with a classic design per slice. We considers the case of the excitation of 52 axial slices (Q = 0) with a thickness of 1 mm and respecting a distance of 1.2 mm (center to center) between two adjacent slices. For each cutting position, the best parametrization n s no is calculated according to the criterion: [0149] Where
[0150]^2d = 1 P C(n0, z) [0150]^ 2d = 1 PC(n 0 , z)
[0151] Est résolu de manière numérique. [0151] Is solved numerically.
[0152] Comme pour l’optimisation précédente, l’algorithme SQP est utilisé pour le calcul des nz s no optimales (52 optimisations indépendantes). Pour son initialisation, on part de la paramétrisation P0 = (A0, 0) où A0 est une matrice constante et qui satisfait : As for the previous optimization, the SQP algorithm is used for calculating the non- optimal n z s (52 independent optimizations). For its initialization, we start from the parametrization P 0 = (A 0 , 0) where A 0 is a constant matrix and which satisfies:
[0154] En partant d’une valeur initiale du critère d’optimisation de 0,26, l’algorithme SQP converge à une valeur de 0.04 en 180 itérations environ. Starting from an initial value of the optimization criterion of 0.26, the SQP algorithm converges to a value of 0.04 in approximately 180 iterations.
[0155] Etant donné la solution trouvée, Xopt, on peut calculer pour tout couple ( z , Q) le profil d’angle de bascule dans la coupe réalisé par On peut alors calculer l’écart quadratique moyen à la consigne de l’angle de bascule, normalisé à la consigne d’angle de bascule (FA-NRMSE). Cette carte de FA-NRMSE est représentée à la [Fig.4] On vérifier que, sauf pour des valeurs très extrêmes de z, qui sont peut significative car l’intersection entre la région d’intérêt (cerveau) et la coupe est très petite, l’écart quadratique moyen à la consigne n’excède pas 5%. [0155] Given the solution found, X opt , we can calculate for any couple ( z , Q) the tilt angle profile in the cut performed by We can then calculate the mean square deviation from the setpoint of the angle tilt, normalized to the tilt angle setpoint (FA-NRMSE). This map of FA-NRMSE is represented in [Fig.4] We check that, except for very extreme values of z, which are not significant because the intersection between the region of interest (brain) and the cut is very small, the root mean square deviation from the setpoint does not exceed 5%.
[0156] Sur le même principe, il est possible de calculer pour tout , sur chaque canal et chaque sous-impulsion, l’énergie des sous-impulsions de nz n0(Xopt). On trouve que le maximum des énergies des impulsions en fonction de ( z , Q) n’excède pas 0,16 Joules. [0156] On the same principle, it is possible to calculate for all , on each channel and each sub-pulse, the energy of the sub-pulses of n z n0 (X opt ). It is found that the maximum of the pulse energies as a function of (z, Q) does not exceed 0.16 Joules.
[0157] Les [Fig. 5a] et [Fig. 5b] montrent une représentation graphique partielle de pour ([Fig. 5a]) et pour = 0 ([Fig. 5b]). On peut vérifier en particulier le caractère constant de lorsque Par ailleurs, on vérifie que est p -périodique en Q. [0157] The [Figs. 5a] and [Fig. 5b] show a partial graphical representation of for ([Fig. 5a]) and for = 0 ([Fig. 5b]). You can check in particular the constant character of when Moreover, we check that is p -periodic in Q.
[0158] La [Fig. 6] permet de comparer c’est-à-dire deux séquences définissant des protocoles multi-coupes (52 coupes) obtenus conformément à l’invention (paramétrisation globale) et à l’art antérieur (optimisation coupe par coupe), respectivement. Le panneau de gauche représente les coefficients An (partie réelle) de la matrice A et le panneau de gauche les coefficients Bu de la matrice B. On vérifie que l’approche paramétrique retourne bien une évolution continue et régulière des matrices A et B lorsque l’on parcourt les coupes, et que ce n’est pas le cas avec l’approche non-paramétrique où est optimisée sans prendre en compte les coupes voisines. [0158] The [Fig. 6] allows to compare that is to say two sequences defining multi-cut protocols (52 cuts) obtained in accordance with the invention (global parameterization) and the prior art (optimization cut by cut), respectively. The left panel represents the coefficients An (real part) of the matrix A and the left panel the coefficients Bu of the matrix B. We check that the parametric approach indeed returns a continuous and regular evolution of the matrices A and B when the 'the cuts are traversed, and this is not the case with the non-parametric approach where is optimized without taking into account neighboring cuts.
[0159] Les [Fig. 7a], [Fig. 7b] et [Fig. 7c] sont des cartes de la distribution de l’angle de bascule dans un plan sagittal , un plan coronal y=0 ([Fig. 7b]) et un plan axial obtenues conformément à l’invention (images de gauche) et par l’approche non-paramétrique définie par n s n0 (images de droite). Ces figures montrent l’effet du caractère discontinu de n s n0 sur le profil d’angle de bascule. L’angle de bascule a bien un profil régulier dans les plans de coupe mais montre en revanche des discontinuités lorsqu’il est observé dans un plan orthogonal, par exemple sagittal ou coronal. Ce problème est résolu avec la paramétrisation globale selon l’invention. [0159] The [Figs. 7a], [Fig. 7b] and [Fig. 7c] are maps of the tilt angle distribution in a sagittal plane , a coronal plane y=0 ([Fig. 7b]) and an axial plane obtained in accordance with the invention (left images) and by the non-parametric approach defined by n s n0 (right images). These figures show the effect of the discontinuous nature of n s n0 on the tilt angle profile. The tilt angle does indeed have a regular profile in the section planes but, on the other hand, shows discontinuities when it is observed in an orthogonal plane, for example sagittal or coronal. This problem is solved with the global parametrization according to the invention.
[0160] La [Fig. 8] est un organigramme d’un procédé d’imagerie par résonance magnétique mettant en œuvre différents aspects de l’invention. [0160] The [Fig. 8] is a flowchart of a magnetic resonance imaging method implementing various aspects of the invention.
[0161] Une première étape i du procédé consiste à concevoir une pluralité de séquences d’impulsions d’excitations de spin de type multi-rayon, chaque dite séquence d’impulsions correspondant à une coupe respective. Conformément à l’invention, cette première étape se décompose en deux sous-étapes : une première sous-étape a) consistant à déterminer une paramétrisation globale ou « méta- impulsion » par la méthode décrite plus haut ; et une deuxième étape b) consistant à obtenir les séquences d’excitations correspondant à chaque coupe à partir de ce paramétrage global et des coordonnées de position et d’orientation de la coupe. [0162] La sous-étape a) comporte trois phases : A first step i of the method consists in designing a plurality of sequences of spin excitation pulses of the multi-ray type, each said sequence of pulses corresponding to a respective slice. In accordance with the invention, this first step is broken down into two sub-steps: a first sub-step a) consisting in determining a global parametrization or “meta-pulse” by the method described above; and a second step b) consisting in obtaining the sequences of excitations corresponding to each cut from this global parameter setting and the position and orientation coordinates of the chopped off. [0162] Sub-step a) comprises three phases:
[0163] A) : l’acquisition d’une carte de pour chaque canal d’émission. [0163] A): the acquisition of a map for each transmission channel.
[0164] B) : La définition d’un premier ensemble de paramètres caractérisant les impulsions radiofréquence et d’un second ensemble de paramètres caractérisant les impulsions de gradient magnétique. Comme expliqué plus haut la valeur de chaque paramètre de chaque dit ensemble étant exprimée par une série tronquée de fonctions d’une coordonnée de position et d’une paire de coordonnées angulaires d’orientation d’une dite coupe, chaque série étant déterminée par ses coefficients C) Le calcul des valeurs des coefficients de manière à minimiser l’écart quadratique moyen (ou une autre fonction de coût adaptée) à une consigne de l’excitation des spins nucléaires, la moyenne étant calculée sur ledit volume d’intérêt et sur toutes les orientations et positions possibles desdites coupes. [0164] B): The definition of a first set of parameters characterizing the radiofrequency pulses and a second set of parameters characterizing magnetic gradient pulses. As explained above, the value of each parameter of each said set being expressed by a truncated series of functions of a position coordinate and a pair of angular orientation coordinates of a said section, each series being determined by its coefficients C) The calculation of the values of the coefficients so as to minimizing the root mean square deviation (or another suitable cost function) at a setpoint of the excitation of the nuclear spins, the average being calculated over said volume of interest and over all the possible orientations and positions of said slices.
[0165] La deuxième étape consiste à appliquer au corps, en succession, les séquences d’impulsions d’excitations de spin ainsi calculée, avec les gradients de sélection correspondants. [0165] The second stage consists in applying to the body, in succession, the sequences of spin excitation pulses thus calculated, with the corresponding selection gradients.
[0166] La troisième étape iii) consiste à acquérir les signaux de résonance magnétique après chaque séquence d’excitation. The third step iii) consists in acquiring the magnetic resonance signals after each excitation sequence.
[0167] Et la quatrième étape iv) consiste à traiter les signaux ainsi acquis pour reconstruire une image d’une portion dudit volume de référence définie par l’union des coupes. And the fourth step iv) consists in processing the signals thus acquired to reconstruct an image of a portion of said reference volume defined by the union of the slices.
[0168] Les étapes ii) - iv) sont conventionnelles en soi, la contribution technique de l’invention résidant principalement dans la détermination de la paramétrisation globale et dans son application à la conception des séquences d’excitation effectivement utilisées. The steps ii) - iv) are conventional in themselves, the technical contribution of the invention residing mainly in the determination of the global parametrization and in its application to the design of the excitation sequences actually used.
[0169] L’invention a été décrite en référence à un mode de réalisation particulier, mais des variantes sont possibles. Par exemple : The invention has been described with reference to a particular embodiment, but variants are possible. For instance :
[0170] Un seul canal d’émission peut être utilisé, au lieu d’une pluralité de canaux réalisant une transmission parallèle. [0171] Il est possible de réaliser une excitation 2d multi bande (adressage simultané de plusieurs coupe parallèles) [Moller2010]. Dans ce cas il est nécessaire d’imposer pour tout ^ e la forme c’est-à-dire indépendant de Cette condition assure en effet que toutes les impulsions des séquences appliquées de manière simultanée (une séquence par coupe) aient les mêmes positions dans l’espace K de transmission. A single transmission channel can be used, instead of a plurality of channels performing parallel transmission. [0171] It is possible to perform a 2d multi-band excitation (simultaneous addressing of several parallel cuts) [Moller2010]. In this case it is necessary to impose for all ^ e form i.e. independent of This condition in fact ensures that all the pulses of the sequences applied simultaneously (one sequence per cut) have the same positions in the transmission space K.
[0172] La démonstration de principe a été faite dans le cas d’une optimisation de l’excitation spécifique à une carte de B mesurée. Cependant, la méthode se généralise au cas du design simultané, donnant lieu à des solutions dites universelles. Pour ce faire, suivant l’approche décrite dans EP 15306569 et EP3594710, il est nécessaire d’acquérir des cartes de B pour une pluralité de corps (sujets d’IRM) ayant des propriétés électromagnétiques différentes et optimiser conjointement la paramétrisation globale pour l’ensemble des sujets. Cela permet d’obtenir une paramétrisation globale « universelle », adaptée à une catégorie de corps à imager, ou pouvant être « personnalisée » en fonction de propriétés aisément mesurables d’un corps déterminé. De cette façon, il n’est pas nécessaire de procéder à une acquisition de cartes de B pour chaque corps à imager. The demonstration of principle was made in the case of an optimization of the excitation specific to a map of B measured. However, the method generalizes to the case of simultaneous design, giving rise to so-called universal solutions. To do this, following the approach described in EP 15306569 and EP3594710, it is necessary to acquire maps of B for a plurality of bodies (MRI subjects) with different electromagnetic properties and jointly optimize the global parameterization for the set of subjects. This makes it possible to obtain a “universal” global parametrization, adapted to a category of body to be imaged, or that can be “customized” according to easily measurable properties of a given body. In this way, it is not necessary to acquire B maps for each body to be imaged.
[0173] La démonstration de principe a été faite dans le cas où l’approximation des petits angles de bascule est valable, mais une généralisation aux grands angles de bascule ne pose pas de difficulté fondamentale. The demonstration of principle has been made in the case where the approximation of small tilt angles is valid, but a generalization to large tilt angles does not pose any fundamental difficulty.
[0174] La forme de la paramétrisation choisie (polynomiale en z et harmoniques sphériques en q, f) est particulièrement adapté au problème mais n’est pas la seule qui soit éligible. La troncature en z est un artifice utile pour pallier au fait que les fonctions polynomiales divergent à l’infini. Là encore, d’autres approches sont potentiellement envisageables. The form of the parametrization chosen (polynomial in z and spherical harmonics in q, f) is particularly suited to the problem but is not the only one that is eligible. The truncation in z is a useful artifice to compensate for the fact that the polynomial functions diverge at infinity. Here again, other approaches are potentially possible.
[0175] L’optimisation de l’écart quadratique moyen de l’excitation à la cible est un critère usuel et avantageux sur le plan mathématique, mais d’autres fonctions de l’écart à la cible peuvent être utilisées. Par ailleurs, bien que ce soit le cas le plus usuel, il n’est pas essentiel que l’excitation-cible soit constante sur toute la coupe. The optimization of the mean square deviation of the excitation from the target is a usual criterion and advantageous from a mathematical point of view, but other functions of the deviation from the target can be used. Moreover, although this is the most common case, it is not essential that the target excitation be constant over the entire slice.
[0176] L’invention ne présuppose pas un type particulier d’impulsions. Le mode de réalisation décrit en détail est le plus simple, impliquant un gradient magnétique constant combiné à une forme d’onde radiofréquence de type sinus cardinal apodisé. Néanmoins, l’invention peut être mis en oeuvre avec d’autres impulsions, par exemple de type SLR [Shinnar1994, Sharma2015] et/ou VERSE [Conolly1988] The invention does not presuppose a particular type of pulse. The embodiment described in detail is the simplest, involving a constant magnetic gradient combined with an apodized cardinal sine type radio frequency waveform. Nevertheless, the invention can be implemented with other pulses, for example of the SLR type [Shinnar1994, Sharma2015] and/or VERSE [Conolly1988]
[0177] Références [0177] References
[0178] [Conollyl 988] S. Conolly, D. Nishimura, A. Macovski, et G. Glover, « Variable- rate sélective excitation », Journal of Magnetic Résonance (1969), vol. 78, no 3, p. 440-458, 1988. [0178] [Conollyl 988] S. Conolly, D. Nishimura, A. Macovski, and G. Glover, “Variable-rate selective excitation”, Journal of Magnetic Resonance (1969), vol. 78, No. 3, p. 440-458, 1988.
[0179] [Shinnarl 994] M. Shinnar, « Reduced power sélective excitation radio frequency puises », Magnetic résonance in medicine, vol. 32, no 5, Art. no 5, 1994. [0179] [Shinnarl 994] M. Shinnar, “Reduced power selective excitation radio frequency pulses”, Magnetic resonance in medicine, vol. 32, No. 5, Art. No. 5, 1994.
[0180] [Setsompop 2008] Setsompop, Kawin, Vijayanand Alagappan, Borjan Gagoski, Thomas Witzel, Jonathan Polimeni, Andréas Potthast, Franz Hebrank, et al. « Slice-selective RF puises for in vivo B 1+ inhomogeneity mitigation at 7 tesla using parallel RF excitation with a 16-element coil ». Magnetic Résonance in Medicine 60, no 6 (2008): 1422-1432. [0180] [Setsompop 2008] Setsompop, Kawin, Vijayanand Alagappan, Borjan Gagoski, Thomas Witzel, Jonathan Polimeni, Andréas Potthast, Franz Hebrank, et al. “Slice-selective RF pulses for in vivo B 1+ inhomogeneity mitigation at 7 tesla using parallel RF excitation with a 16-element coil”. Magnetic Resonance in Medicine 60, no 6 (2008): 1422-1432.
[0181 ] [Moller 2010] S. Moeller et al., « Multiband multislice GE-EPI at 7 tesla, with 16-fold accélération using partial parallel imaging with application to high spatial and temporal whole-brain fMRI », Magn. Reson. Med., vol. 63, no 5, Art. no 5, mai 2010[0181 ] [Moller 2010] S. Moeller et al., “Multiband multislice GE-EPI at 7 tesla, with 16-fold acceleration using partial parallel imaging with application to high spatial and temporal whole-brain fMRI”, Magn. Resound. Med., vol. 63, No. 5, Art. No. 5, May 2010
[0182] [Amadon2012] [1]A. Amadon, M. A. Cloos, N. Boulant, M.-F. Hang, C. J. Wiggins, et H. -P. Fautz, « Validation of a very fast B1-mapping sequence for parallel transmission on a human brain at 7T », In Proceedings of the 20th Annual Meeting of ISMRM, p. 3358, 2012. [0182] [Amadon2012] [1]A. Amadon, M.A. Cloos, N. Boulant, M.-F. Hang, C.J. Wiggins, and H. -P. Fautz, “Validation of a very fast B1-mapping sequence for parallel transmission on a human brain at 7T”, In Proceedings of the 20th Annual Meeting of ISMRM, p. 3358, 2012.
[0183] [Poser 2014] Poser, Benedikt A., Robert James Anderson, Bastien Guérin, Kawin Setsompop, Weiran Deng, Azma Mareyam, Peter Serano, Lawrence L. Wald, et V. Andrew Stenger. « Simultaneous multislice excitation by parallel transmission ». Magnetic Résonance in Medicine 71, no 4 (2014): 1416-1427. [0183] [Poser 2014] Poser, Benedikt A., Robert James Anderson, Bastien Guérin, Kawin Setsompop, Weiran Deng, Azma Mareyam, Peter Serano, Lawrence L. Wald, and V. Andrew Stenger. "Simultaneous multislice excitation by parallel transmission". Magnetic Resonance in Medicine 71, no 4 (2014): 1416-1427.
[0184] [Guérin 2015] Guérin, Bastien, Kawin Setsompop, Huihui Ye, Benedikt A. Poser, Andrew V. Stenger, et Lawrence L. Wald. « Design of parallel transmission puises for simultaneous multislice with explicit control for peak power and local spécifie absorption rate ». Magnetic Résonance in Medicine 73, no 5 (2015): 1946- 1953. [0184] [Guérin 2015] Guérin, Bastien, Kawin Setsompop, Huihui Ye, Benedikt A. Poser, Andrew V. Stenger, and Lawrence L. Wald. “Design of parallel transmission pulses for simultaneous multislice with explicit control for peak power and local specifies absorption rate”. Magnetic Resonance in Medicine 73, no 5 (2015): 1946-1953.
[0185] [Sharma2015] A. Sharma, R. Bammer, V. A. Stenger, et W. A. Grissom, « Low peak power multiband spokes puises for B inhomogeneity-compensated simultaneous multislice excitation in high field MRI: Low Peak Power Multiband Spokes Puises », Magnetic Résonance in Medicine, vol. 74, n° 3, Art. n° 3, sept.[0185] [Sharma2015] A. Sharma, R. Bammer, VA Stenger, and WA Grissom, “Low peak power multiband spokes pulses for B inhomogeneity-compensated simultaneous multislice excitation in high field MRI: Low Peak Power Multiband Spokes Pulses”, Magnetic Resonance in Medicine, vol. 74, No. 3, Art. No. 3, Sept.
2015 2015
[0186] [Gras 2017] Gras, Vincent, Alexandre Vignaud, Alexis Amadon, Franck Mauconduit, Denis Le Bihan, et Nicolas Boulant. « In Vivo Démonstration of Whole- Brain Multislice Multispoke Parallel Transmit Radiofrequency Puise Design in the Small and Large Flip Angle Régimes at 7 Tesla ». Magnetic Résonance in Medicine 78, n° 3 (2017): 1009-19. [0186] [Gras 2017] Gras, Vincent, Alexandre Vignaud, Alexis Amadon, Franck Mauconduit, Denis Le Bihan, and Nicolas Boulant. “In Vivo Demonstration of Whole-Brain Multislice Multispoke Parallel Transmit Radiofrequency Pulse Design in the Small and Large Flip Angle Regimes at 7 Tesla”. Magnetic Resonance in Medicine 78, No. 3 (2017): 1009-19.
[0187] [Gras2017 B] V. Gras, A. Vignaud, A. Amadon, F. Mauconduit, D. Le Bihan, et N. Boulant, « New method to characterize and correct with sub-ps précision gradient delays in bipolar multispoke RF puises », Magn. Reson. Med, vol. 78, p. 2194-2202, janv. 2017. [0187] [Gras2017 B] V. Gras, A. Vignaud, A. Amadon, F. Mauconduit, D. Le Bihan, and N. Boulant, “New method to characterize and correct with sub-ps precision gradient delays in bipolar multispoke RF pulses”, Magn. Resound. Med, vol. 78, p. 2194-2202, Jan. 2017.
[0188][Cao 2015] : Z. Cao et al. « Joint design of large-tip-angle parallel RF puises and blipped gradient trajectories », Magnetic Résonance in Medicine, vol. 75, no. 3, 27 avril 2015. [0188][Cao 2015]: Z. Cao et al. "Joint design of large-tip-angle parallel RF pulses and blipped gradient trajectories", Magnetic Resonance in Medicine, vol. 75, no. 3, April 27, 2015.

Claims

REVENDICATIONS
1. Procédé de détermination d’une paramétrisation globale d’une famille de séquences d’impulsions d’excitations de spin en imagerie par résonance magnétique, chaque séquence d’impulsions de ladite famille étant une séquence de type multi-rayon adaptée pour exciter sélectivement des spins nucléaires dans une coupe (CP) respective d’un volume d’intérêt (ROI) d’un corps (C) immergé dans un champ magnétique stationnaire (B0) et comprenant un même nombre prédéterminé d’impulsions radiofréquence (IRF1, IRF2, IRF3) à une fréquence de Larmor desdits spins nucléaires alternées à des impulsions de gradient magnétique (IG1, IG2, IG3), le procédé comprenant les étapes consistant à : 1. Method for determining an overall parametrization of a family of sequences of spin excitation pulses in magnetic resonance imaging, each sequence of pulses of said family being a sequence of the multi-ray type adapted to selectively excite nuclear spins in a respective section (CP) of a volume of interest (ROI) of a body (C) immersed in a stationary magnetic field (B 0 ) and comprising the same predetermined number of radiofrequency pulses (IRF1, IRF2, IRF3) at a Larmor frequency of said nuclear spins alternating with magnetic gradient pulses (IG1, IG2, IG3), the method comprising the steps of:
A) Acquérir au moins une carte d’intensité d’un rayonnement radiofréquence à la fréquence de Larmor dans ledit volume d’intérêt ; A) Acquire at least one intensity map of radiofrequency radiation at the Larmor frequency in said volume of interest;
B) Définir un premier ensemble de paramètres caractérisant les impulsions radiofréquence et un second ensemble de paramètres caractérisant les impulsions de gradient magnétique, la valeur de chaque paramètre de chaque dit ensemble étant exprimée par une série tronquée de fonctions d’une coordonnée de position (z) et d’une paire de coordonnées angulaires d’orientation (f, Q) d’une dite coupe, chaque série étant déterminée par ses coefficients ; B) Define a first set of parameters characterizing the radiofrequency pulses and a second set of parameters characterizing the magnetic gradient pulses, the value of each parameter of each said set being expressed by a truncated series of functions of a position coordinate (z ) and a pair of angular orientation coordinates (f, Q) of a said section, each series being determined by its coefficients;
C) A partir de ladite ou desdites cartes d’intensité, calculer les valeurs des coefficients desdites séries tronquées de fonctions qui minimisent une fonction représentative d’un écart moyen à une consigne prédéterminée de l’excitation des spins nucléaires, ledit écart moyen étant calculé sur ledit volume d’intérêt et sur toutes les orientations et positions possibles desdites coupes ; lesdites valeurs constituant ladite paramétrisation globale. C) From said intensity map(s), calculate the values of the coefficients of said truncated series of functions which minimize a function representative of an average deviation from a predetermined setpoint of the excitation of the nuclear spins, said average deviation being calculated on said volume of interest and on all possible orientations and positions of said sections; said values constituting said global parametrization.
2. Procédé selon la revendication 1 dans lequel chaque terme de chaque dite série tronquée est donné par le produit d’un monôme de ladite coordonnée de position et d’une somme finie de fonctions appartenant à une base de fonctions continues desdites coordonnées angulaires d’orientation. 2. Method according to claim 1, in which each term of each said truncated series is given by the product of a monomial of the said position coordinate and of a finite sum of functions belonging to a basis of continuous functions of the said angular coordinates of orientation.
3. Procédé selon la revendication 2 dans lequel chaque terme de chaque dite série tronquée exprimant chaque paramètre du premier ensemble est donné par le produit d’un monôme de ladite coordonnée de position et d’une somme finie d’harmoniques sphériques complexes de symétrie paire et chaque terme de chaque dite somme tronquée exprimant chaque paramètre du second ensemble est donné par le produit d’un monôme de ladite coordonnée de position et d’une somme finie d’harmoniques sphériques réelles de symétrie paire. 3. Method according to claim 2 wherein each term of each said truncated series expressing each parameter of the first set is given by the product of a monomial of said position coordinate and a finite sum of complex spherical harmonics of even symmetry and each term of each said truncated sum expressing each parameter of the second set is given by the product of a monomial of said position coordinate and a finite sum of real spherical harmonics of even symmetry.
4. Procédé selon l’une des revendications précédentes dans lequel l’étape C) comporte la minimisation d’un écart quadratique moyen de l’excitation des spins nucléaires à une consigne, la moyenne étant calculée sur ledit volume d’intérêt et sur toutes les orientations possibles desdites coupes. 4. Method according to one of the preceding claims, in which step C) comprises the minimization of a mean square deviation of the excitation of the nuclear spins at a setpoint, the average being calculated on said volume of interest and on all the possible orientations of said slices.
5. Procédé selon l’une des revendications précédentes dans lequel l’étape C) comprend une minimisation sous contrainte de ladite fonction représentative d’un écart moyen de l’excitation des spins nucléaires à une consigne. 5. Method according to one of the preceding claims, in which step C) comprises a minimization under constraint of said function representing an average deviation of the excitation of the nuclear spins from a setpoint.
6. Procédé selon la revendication 5 dans lequel la minimisation est effectuée sous une ou plusieurs contraintes d’énergie desdites impulsions radiofréquence. 6. Method according to claim 5, in which the minimization is carried out under one or more energy constraints of said radio frequency pulses.
7. Procédé selon l’une des revendications précédentes dans lequel l’étape A) comprend l’acquisition d’une pluralité de cartes d’intensité du rayonnement radiofréquence associées à des canaux d’émission respectifs de rayonnement radiofréquence, chaque impulsion radiofréquence de la séquence étant constituée de sous-impulsions destinées à être émises par des canaux d’émission respectifs. 7. Method according to one of the preceding claims, in which step A) comprises the acquisition of a plurality of intensity maps of the radiofrequency radiation associated with respective emission channels of radiofrequency radiation, each radiofrequency pulse of the sequence consisting of sub-pulses intended to be transmitted by respective transmission channels.
8. Procédé selon l’une des revendications précédentes dans lequel l’étape8. Method according to one of the preceding claims, in which the step
A) comprend l’acquisition d’une pluralité de cartes d’intensité du rayonnement radiofréquence associées à des corps respectifs présentant des propriétés électromagnétiques différentes, l’écart moyen à une consigne de l’excitation des spins nucléaires étant également calculée sur l’ensemble desdits corps. A) comprises the acquisition of a plurality of intensity maps of the radiofrequency radiation associated with respective bodies having different electromagnetic properties, the mean deviation from a setpoint of the excitation of the nuclear spins also being calculated over the whole of said bodies.
9. Procédé de conception d’une séquence d’impulsions d’excitations de spin de type multi-rayon en imagerie par résonance magnétique, ladite séquence d’impulsions étant adaptée pour exciter sélectivement des spins nucléaires dans une coupe (CP) respective d’un volume d’intérêt (ROI) d’un corps (C) immergé dans un champ magnétique stationnaire (B0) et comprenant un même nombre prédéterminé d’impulsions radiofréquence (IRF1, IRF2, IRF3) à une fréquence de Larmor desdits spins nucléaires alternées à des impulsions de gradient magnétique (IG1, IG2, IG3), la coupe étant identifiée par des valeurs d’une coordonnée de position (z) et d’une paire de coordonnées angulaires d’orientation (f, Q), le procédé étant caractérisé en ce qu’il comporte les étapes consistant à : a) Procurer une paramétrisation globale d’une famille de telles séquences d’excitation, ladite paramétrisation globale étant constituée par une première et une deuxième pluralité d’ensembles de coefficients de séries tronquées respectives de fonctions desdites coordonnées de position et d’orientation d’une dite coupe, chaque ensemble de la première pluralité d’ensembles étant constitué des coefficients d’une série tronquée de fonctions définissant un paramètre d’une impulsion radiofréquence et chaque ensemble de la deuxième pluralité d’ensembles étant constitué des coefficients d’une série tronquée de fonctions définissant un paramètre d’une impulsion de gradient magnétique ; et b) Calculer les valeurs desdits paramètres d’impulsions de gradient et d’impulsions radiofréquence à partir de ladite paramétrisation globale, desdites séries tronquées de fonctions et des valeurs de ladite coordonnée de position de ladite paire de coordonnées angulaires d’orientation de la coupe. 9. Method for designing a sequence of spin excitation pulses of the multi-ray type in magnetic resonance imaging, said sequence of pulses being adapted to selectively excite nuclear spins in a respective slice (CP) of a volume of interest (ROI) of a body (C) immersed in a stationary magnetic field (B 0 ) and comprising the same predetermined number of radio frequency pulses (IRF1, IRF2, IRF3) at a Larmor frequency of said nuclear spins alternated with magnetic gradient pulses (IG1, IG2, IG3), the cut being identified by values of a position coordinate (z) and a pair of angular orientation coordinates (f, Q), the method being characterized in that it includes the steps of: a) providing a global parametrization of a family of such excitation sequences, said global parametrization being constituted by a first and a second plurality of sets of truncated series coefficients r respective functions of said coordinates of position and orientation of a said cut, each set of the first plurality of sets consisting of the coefficients of a truncated series of functions defining a parameter of a radio frequency pulse and each set of the second plurality of sets consisting of the coefficients of a truncated series of functions defining a parameter of a magnetic gradient pulse; and b) Calculating the values of said gradient pulse and radio frequency pulse parameters from said global parameterization, said truncated series of functions and the values of said position coordinate of said pair of angular coordinates of orientation of the cut.
10. Procédé selon la revendication 9, dans lequel la paramétrisation globale est obtenue par un procédé selon l’une des revendications 1 à 8. 10. Method according to claim 9, in which the global parametrization is obtained by a method according to one of claims 1 to 8.
11. Procédé d’imagerie par résonance magnétique comprenant les étapes suivantes : i) Concevoir une pluralité de séquences d’impulsions d’excitations de spin de type multi-rayon, chaque dite séquence d’impulsions étant adaptée à exciter sélectivement des spins nucléaires dans une coupe respective (CP) d’un volume d’intérêt (ROI) d’un corps (C) à imager immergé dans un champ magnétique stationnaire (B0) et comprenant des impulsions radiofréquence (IRF1 , IRF2, IRF3) à une fréquence de Larmor desdits spins nucléaires alternées à des impulsions de gradient magnétique (IG1, IG2, IG3), chaque coupe étant identifiée par des valeurs respectives d’une coordonnée de position (z) et d’une paire de coordonnées angulaires d’orientation (f, Q) ; ii) Appliquer audit corps, en succession, lesdites séquences d’impulsions d’excitations de spin, un gradient magnétique de sélection de coupe étant appliqué en même temps que chaque impulsion radiofréquence ; iii) Après chaque dite séquence d’impulsions d’excitation de spin, acquérir un signal de résonance magnétique émis par les spins nucléaires ; et iv) Traiter les signaux ainsi acquis pour reconstruire une image d’une portion dudit volume de référence définie par l’union desdites coupes ; caractérisé en ce que l’étape i) est mise en œuvre selon une des revendications 9 ou 10. 11. Magnetic resonance imaging method comprising the following steps: i) Designing a plurality of sequences of spin excitation pulses of the multi-ray type, each said sequence of pulses being adapted to selectively excite nuclear spins in a respective section (CP) of a volume of interest (ROI) of a body (C) to be imaged immersed in a stationary magnetic field (B 0 ) and comprising radiofrequency pulses (IRF1, IRF2, IRF3) at a frequency of Larmor of said alternating nuclear spins at magnetic gradient pulses (IG1, IG2, IG3), each cut being identified by respective values of a position coordinate (z) and a pair of angular coordinates of orientation (f , Q); ii) Applying to said body, in succession, said sequences of spin excitation pulses, a magnetic cut selection gradient being applied at the same time as each radiofrequency pulse; iii) After each said sequence of spin excitation pulses, acquiring a magnetic resonance signal emitted by the nuclear spins; and iv) processing the signals thus acquired to reconstruct an image of a portion of said reference volume defined by the union of said sections; characterized in that step i) is implemented according to one of claims 9 or 10.
12. Ordinateur (OC) programmé pour : a) Recevoir en entrée au moins une carte d’intensité d’un rayonnement radiofréquence dans un volume d’intérêt (ROI) d’un corps (C) immergé dans un champ magnétique stationnaire (B0), ledit rayonnement radiofréquence étant à une fréquence de Larmor de spins nucléaires dudit corps ; b) Définir une paramétrisation globale d’une famille de séquences d’impulsions d’excitations de spin de type multi-rayon en imagerie par résonance magnétique, chaque séquence d’impulsions de ladite famille étant adaptée pour exciter sélectivement des spins nucléaires dans une coupe dudit volume d’intérêt et comprenant des impulsions radiofréquence (IRF1, IRF2, IRF3) à une fréquence de Larmor desdits spins nucléaires alternées à des impulsions de gradient magnétique (IG1, IG2, IG3), ladite paramétrisation globale comprenant un premier ensemble de paramètres caractérisant lesdites impulsions radiofréquence et un second ensemble de paramètres caractérisant lesdites impulsions de gradient magnétique, la valeur de chaque paramètre de chaque dit ensemble étant exprimée par une série tronquée de fonctions d’une coordonnée de position (z) et d’une paire de coordonnées angulaires d’orientation (f, Q) d’une dite coupe, chaque série étant déterminée par ses coefficients ; g) A partir de ladite ou desdites cartes d’intensité, calculer les valeurs des coefficients desdites séries tronquées de fonctions qui minimisent une fonction représentative d’un écart moyen à une consigne de l’excitation des spins nucléaires, la moyenne étant calculée sur ledit volume d’intérêt et sur toutes les orientations et positions possibles desdites coupes. 12. Computer (OC) programmed to: a) Receive as input at least one intensity map of radiofrequency radiation in a volume of interest (ROI) of a body (C) immersed in a stationary magnetic field (B 0 ), said radiofrequency radiation being at a Larmor frequency of nuclear spins of said body; b) Define a global parametrization of a family of sequences of spin excitation pulses of the multi-ray type in magnetic resonance imaging, each sequence of pulses of said family being adapted to selectively excite nuclear spins in a section said volume of interest and comprising radiofrequency pulses (IRF1, IRF2, IRF3) at a Larmor frequency of said nuclear spins alternating with magnetic gradient pulses (IG1, IG2, IG3), said global parameterization comprising a first set of parameters characterizing said radio frequency pulses and a second set of parameters characterizing said magnetic gradient pulses, the value of each parameter of each said set being expressed by a truncated series of functions of a position coordinate (z) and a pair of orientation angular coordinates (f, Q) of a so-called cut, each series being determined by its coefficients; g) From said intensity map or maps, calculate the values of the coefficients of said truncated series of functions which minimize a function representative of an average deviation from a setpoint of the excitation of the nuclear spins, the average being calculated on said volume of interest and on all possible orientations and positions of said sections.
13. Ordinateur selon la revendication 12, également programmé pour : d) Recevoir en entrée une pluralité de valeurs de ladite coordonnée de position et de paires respectives desdites coordonnées angulaires définissant une pluralité de coupes ; e) Pour chaque dite coupe, concevoir une séquence d’impulsions d’excitation de spin en calculant des valeurs desdits paramètres d’impulsions radiofréquence et d’impulsions de gradient magnétique à partir de la paramétrisation globale calculée à l’étape g, desdites séries tronquées de fonctions et des valeurs de ladite coordonnée de position et de ladite paire de coordonnées angulaires d’orientation de la coupe ; z) Appliquer audit corps, en succession, lesdites séquences d’impulsions d’excitations de spin, un gradient magnétique de sélection de coupe étant appliqué en même temps que chaque impulsion radiofréquence ; h) Après chaque impulsion, acquérir un signal de résonance magnétique émis par les spins nucléaires ; 13. A computer according to claim 12, also programmed to: d) receive as input a plurality of values of said position coordinate and respective pairs of said angular coordinates defining a plurality of cuts; e) For each said slice, design a sequence of spin excitation pulses by calculating values of said parameters of radio frequency pulses and magnetic gradient pulses from the global parametrization calculated in step g, of said series truncated functions and values of said position coordinate and of said pair of angular coordinates of orientation of the cut; z) Applying to said body, in succession, said sequences of spin excitation pulses, a magnetic cut selection gradient being applied at the same time as each radio frequency pulse; h) After each pulse, acquire a magnetic resonance signal emitted by the nuclear spins;
Q) Traiter les signaux ainsi acquis pour reconstruire une image d’une portion dudit volume de référence définie par l’union desdites coupes. Q) Process the signals thus acquired to reconstruct an image of a portion of said reference volume defined by the union of said slices.
14. Ordinateur (OC) programmé pour la mise en oeuvre d’un procédé selon l’une des revendications 9 ou 10. 14. Computer (OC) programmed for the implementation of a method according to one of claims 9 or 10.
15. Appareil d’imagerie par résonance magnétique (SI) équipé d’un ordinateur selon l’une des revendications 13 et 14. 15. Magnetic resonance imaging device (SI) equipped with a computer according to one of claims 13 and 14.
EP22717061.0A 2021-03-17 2022-03-11 Method for generating multi-ray spin-excitation sequences and application thereof to magnetic resonance imaging Pending EP4308949A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR2102654A FR3120949B1 (en) 2021-03-17 2021-03-17 Method for generating multi-ray type spin excitation sequences and its application to magnetic resonance imaging
PCT/EP2022/056381 WO2022194711A1 (en) 2021-03-17 2022-03-11 Method for generating multi-ray spin-excitation sequences and application thereof to magnetic resonance imaging

Publications (1)

Publication Number Publication Date
EP4308949A1 true EP4308949A1 (en) 2024-01-24

Family

ID=75850333

Family Applications (1)

Application Number Title Priority Date Filing Date
EP22717061.0A Pending EP4308949A1 (en) 2021-03-17 2022-03-11 Method for generating multi-ray spin-excitation sequences and application thereof to magnetic resonance imaging

Country Status (3)

Country Link
EP (1) EP4308949A1 (en)
FR (1) FR3120949B1 (en)
WO (1) WO2022194711A1 (en)

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3594710B1 (en) 2018-07-09 2024-04-17 Commissariat à l'Energie Atomique et aux Energies Alternatives A computer-implemented method of building a database of pulse sequences for magnetic resonance imaging, and a method of performing magnetic resonance imaging using such a database

Also Published As

Publication number Publication date
FR3120949B1 (en) 2023-03-17
WO2022194711A1 (en) 2022-09-22
FR3120949A1 (en) 2022-09-23

Similar Documents

Publication Publication Date Title
US7375523B1 (en) System and method for fast MR coil sensitivity mapping
US8698495B2 (en) Flip angle imaging with improved B1 mapping for multi-RF transmit systems
US20130144156A1 (en) Method and Apparatus for Correcting B1-Inhomogeneity in Slice-Selective Nuclear Magnetic Resonance Imaging
CH634412A5 (en) APPARATUS FOR DETERMINING THE DENSITY OF CORES BY SPIN ECHO FOR TOMOGRAM ELABORATION.
CN106990375B (en) For speculating without dictionary MR parameter for fingerprint reconstruct
JP6356809B2 (en) Zero echo time MR imaging with water / fat separation
EP2426563A1 (en) Method for calibrating a device with nuclear operation
Wang et al. Fast rotary nonlinear spatial acquisition (FRONSAC) imaging
Malik et al. Spatially resolved extended phase graphs: modeling and design of multipulse sequences with parallel transmission
JP2003052667A (en) Magnetic resonance imaging system
US7038453B2 (en) Method and magnetic resonance tomography apparatus for spatially resolved measurement of the B1 field distribution
Layton et al. Single shot trajectory design for region‐specific imaging using linear and nonlinear magnetic encoding fields
Malik et al. Direct signal control of the steady‐state response of 3D‐FSE sequences
FR3036189A1 (en) METHOD AND DEVICE FOR MAGNETIC RESONANCE IMAGING
Malik et al. Phase relaxed localized excitation pulses for inner volume fast spin echo imaging
JP6762284B2 (en) Magnetic resonance imaging device and noise removal method
EP2399142B1 (en) Use of strongly modulating pulses in mri for providing chemical shift selective flip angles
EP0316312B1 (en) radio-frequency energization method in an nmr experiment
EP4308949A1 (en) Method for generating multi-ray spin-excitation sequences and application thereof to magnetic resonance imaging
FR2579753A1 (en) METHOD AND DEVICE FOR NUCLEAR MAGNETIC RESONANCE IMAGING
EP0208601B1 (en) Nuclear magnetic resonance imaging method
JP2023134495A (en) Zero echo time MR imaging with water/fat separation
De Santis et al. Shannon information of the geomagnetic field for the past 7000 years
EP2539727B1 (en) Method for creating a module for controlling a high-frequency antenna for a nuclear magnetic resonance imaging apparatus
JP5342155B2 (en) Magnetic resonance imaging system

Legal Events

Date Code Title Description
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: UNKNOWN

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE

PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20231006

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR