Brevet FR3115884A1 — publié le 20220506

L’invention concerne un procédé (100) de post-traitement d’un premier ensemble d’échantillons mesurant la magnitude d’un signal WASAB1 délivré par un appareil d’imagerie médicale par résonance magnétique (1). Pour permettre notamment à la technologie CEST d’accéder au domaine clinique, un tel procédé comporte une étape (110) pour détecter les échantillons du premier ensemble Z pour lesquels les polarités respectives des valeurs du signal mesuré sont connues et construire un deuxième ensemble Y d’échantillons « polarisés ». Un tel procédé (100) comporte dès lors une étape (120) pour ajuster un modèle déterminé audit deuxième ensemble Y, les deux paramètres dudit modèle déterminé décrivant respectivement les champs magnétiques statique B0 et d’excitation B1 de l’appareil d’imagerie médicale par résonance magnétique (1) et produire une estimation desdits paramètres B0 et B1 dudit modèle. Un tel procédé (100) concerne toute application d’imagerie par résonance magnétique pour laquelle une correction d’inhomogénéités desdits champs B0 et B1 est requise. Figure à publier avec l’abrégé : Fig. 9 Procédé de post-traitement d’un échantillonnage d’une séquence d’acquisition WASAB1 L’invention concerne un procédé de post-traitement d’une séquence d’acquisition WASAB1, telle que divulguée dans le document « Simultaneous mapping of water shift and B1 (WASAB1) — application to field‐inhomogeneity correction of CEST MRI data. Magnetic resonance in medicine. 2017 Feb - Schuenke P, Windschuh J, Roeloffs V, Ladd ME, Bachert P, Zaiss M. », ci-après désigné par « Publication ». Une telle séquence WASAB1 permet de cartographier simultanément le champ magnétique statique B0 et le champ magnétique d’excitation B1 d’un appareil d’imagerie par résonance magnétique, plus connu sous le sigle IRM ou MRI, pour Magnetic Resonance Imaging , selon une terminologie anglo-saxonne. Un tel procédé de post-traitement trouve naturellement son application dans le cadre de l’imagerie CEST (Transfert de Saturation d’Échanges Chimiques ou Chemical Exchange Saturation Transfer , selon une terminologie anglo-saxonne), pour laquelle la séquence WASAB1 a été inventée mais, plus généralement, dans toute application d’imagerie par résonance magnétique pour laquelle une correction d’inhomogénéités desdits champs B0 et B1 est requise. L’application dudit procédé de post-traitement selon l’invention, lorsqu’il est appliqué à l’imagerie CEST, permet d’exploiter cette technique dans le domaine clinique et ainsi de quantifier un biomarqueur en lien avec une espèce chimique d’intérêt d’un organe humain ou animal pour, in fine , caractériser différentes lésions aux propriétés métaboliques altérées telles que des tumeurs, des tissus ischémiques ou une sclérose en plaques. L'imagerie par résonance magnétique est basée sur une analyse de la réponse du proton d'une molécule d'eau lorsqu'il est excité dans un champ magnétique. Cette réponse dépend de l'environnement d’un tel proton et permet ainsi de différencier différents types de tissus. Un appareil d'imagerie 1 par Résonance Magnétique Nucléaire, tel qu’illustré à titre d’exemple non limitatif par les figures 1 et 2, est généralement utilisé. Celui-ci délivre une pluralité de séquences d’images numériques 12 d'une ou plusieurs parties du corps d’un patient, à titre d’exemples non limitatifs, le cerveau, le cœur, les poumons. Ledit appareil applique pour cela une combinaison d'ondes électromagnétiques à haute fréquence sur la partie du corps considérée et mesure le signal réémis par certains atomes, tels qu’à titre d’exemple non limitatif, l’hydrogène pour l’imagerie par Résonance Magnétique Nucléaire. L’appareil permet ainsi de déterminer les propriétés magnétiques et, par voie de conséquence, la composition chimique des tissus biologiques et donc leur nature, en chaque volume élémentaire, que l’on nomme communément un voxel, du volume imagé. L’appareil 1 d'imagerie par Résonance Magnétique Nucléaire est commandé à l’aide d’une console 2. Un utilisateur 6, par exemple un opérateur, praticien ou chercheur, peut ainsi choisir des commandes 11 pour piloter l’appareil 1, à partir de paramètres ou de consignes 16 saisis via une interface homme-machine d’entrée 8 du système d’analyse. Une telle interface homme-machine 8 peut consister par exemple en un clavier informatique, un dispositif de pointage, d’un écran tactile, un microphone ou, plus généralement, toute interface agencée pour traduire une gestuelle ou une consigne émise par un humain 6 en données de commande ou de paramétrage. A partir d’informations 10 produites par ledit appareil 1, on obtient une pluralité de séquences d’images numériques 12 d'une partie d’un corps d’un humain ou d’un animal. Nous nommerons également « données expérimentales » de telles informations 10 ou images 12. Les séquences d’images 12 peuvent optionnellement être stockées au sein d’un serveur 3, c’est-à-dire un ordinateur doté de moyens de mémorisation propres et constituer un dossier médical 13 d’un patient. Un tel dossier 13 peut comprendre des images de différents types, telles que des images fonctionnelles mettant en évidence l’activité des tissus ou des images anatomiques reflétant les propriétés des tissus. Les séquences d’images 12 ou, plus généralement les données expérimentales, sont analysées par une unité de traitement 4 agencée à cette fin. Une telle unité de traitement 4 peut, par exemple, consister en un ou plusieurs microprocesseurs ou microcontrôleurs mettant en œuvre des instructions de programmes d’application idoines chargées dans des moyens de mémorisation, tels qu’une mémoire non volatile, dudit système d’analyse d’imagerie. Ladite unité de traitement 4 comporte des moyens pour communiquer avec le monde extérieur pour recueillir les images. Lesdits moyens pour communiquer permettent en outre à l’unité de traitement 4 de délivrer, in fine, un rendu, par exemple graphique et/ou sonore, d’une estimation ou d’une quantification d’un biomarqueur élaborée par ladite unité de traitement 4 à partir des données expérimentales 10 et/ou 12 obtenues par Imagerie par Résonance Magnétique, à un utilisateur 6 du système d’analyse d’imagerie par l’intermédiaire d’une interface homme-machine de sortie 5. Dans tout le document, on entend par « interface homme-machine de sortie », tout dispositif, employé seul ou en combinaison, permettant de sortir ou de délivrer une représentation graphique, haptique, sonore ou, plus généralement, perceptible par l’humain, d’un signal physiologique reconstruit, en l’espèce un biomarqueur, à un utilisateur 6 d’un système d’analyse d’Imagerie par Résonance Magnétique. Une telle interface homme-machine de sortie 5 peut consister de manière non exhaustive en un ou plusieurs écrans, haut-parleurs ou autres moyens alternatifs adaptés. Ledit utilisateur 6 du système d’analyse peut ainsi confirmer ou infirmer un diagnostic, décider d’une action thérapeutique qu’il jugera adéquate, approfondir des travaux de recherche, peaufiner des paramètres de réglage d’un équipement de mesure, etc. Optionnellement, cet utilisateur 6 peut également paramétrer le fonctionnement de l’unité de traitement 4 ou de l’interface homme-machine de sortie 5, au moyen de paramètres 16 de fonctionnement et/ou d’acquisition. Par exemple, il peut ainsi définir des seuils d’affichage ou choisir les biomarqueurs, indicateurs ou paramètres estimés ou quantifiés pour lesquels il souhaite disposer d’une représentation. L’utilisateur exploite pour cela l’interface homme-machine d’entrée 8 précédemment évoquée ou une deuxième interface d’entrée prévue pour cela. Avantageusement, les interfaces homme-machine d’entrée 8 et de sortie 5 peuvent ne constituer qu’une seule et même entité physique. Lesdites interfaces homme-machine d’entrée 8 et de sortie 5 du système d’analyse d’imagerie peuvent également être intégrées à la console d’acquisition 2. Il existe une variante, décrite en liaison avec la , pour laquelle un système d’imagerie, tel que décrit précédemment, comporte en outre une unité de prétraitement 7 pour analyser les séquences d’images 12, en déduire des signaux expérimentaux 15 et délivrer ces derniers à l’unité de traitement 4 qui est ainsi déchargée de cette tâche. Parmi les techniques ou modalités basées sur l’imagerie par résonance magnétique, on distingue l’imagerie par Transfert de Saturation d’Echanges Chimiques ou CEST. Une telle technique consiste à appliquer une impulsion radiofréquence selon différentes fréquences de résonance déterminées ω, de sorte que des espèces chimiques d’intérêt atteignent un état de saturation. Lesdites fréquences de résonance déterminées ω sont en effet respectivement associées à différentes espèces chimiques d'intérêt, telles que, de manière non limitative, le groupe amide (-NH), le groupe amine (-NH2) ou le groupe hydroxyle (-OH). Leurs protons d’hydrogènes labiles ainsi excités sont échangés avec des protons d’hydrogène non excités de l'eau. Une telle application d’une impulsion radiofréquence, à une fréquence de résonance déterminées ω autre que celle ω0 associée à l’eau, continuellement répétée pendant une durée déterminée, par exemple quelques secondes, conduit à une accumulation de saturation en eau d’une espèce chimique d’intérêt. Une concentration de cette dernière peut être indirectement mesurée par la diminution du signal de l'eau, phénomène appelé l’« effet CEST » qui est analysé sous la forme d’un spectre Z (« Z-spectrum » selon une terminologie anglo-saxonne). Une telle diminution peut être facilement détectée par des séquences d’acquisition d'imagerie par résonance magnétique rapides telles que l'imagerie écho-planaire à coup unique également connue sous la terminologie anglo-saxonne « single shot echo-planar imaging ». La présente un exemple d’un spectre Z sous la forme d’un ensemble d’échantillons Z(Δω) (ou « data set » selon une terminologie anglo-saxonne) selon des décalages fréquentiels relatifs Δω par rapport à la fréquence de l'eau, de tels décalages Δω étant exprimés en ppm. Lorsque les échantillons Z(Δω) sont ordonnés selon des Δω décroissants, comme l’indique la , ledit spectre Z peut être décrit sous la forme d’un signal discret décrivant des mesures, avantageusement normalisées par un signal mesuré sans saturation par radiofréquence et exprimées en pourcentages, de la magnitude d’un signal expérimental délivré par un appareil de mesure 1 pour un volume élémentaire d’un organe. Selon la , cet ensemble d’échantillons, lorsque ces derniers sont ordonnés selon des décalages fréquentiels Δω décroissants, constitue un signal discret Z décrivant sensiblement un ‘V’ dont le minimum Z(Δω)m est associé au décalage fréquentiel Δω=1 ppm. En toute logique, ledit minimum Z(Δω)m devrait être associé à un décalage fréquentiel Δω nul si ledit appareil était parfaitement réglé et/ou le champ magnétique statique était homogène. La technique CEST améliore la détection de certains métabolites du corps humain dont la concentration est insuffisante pour être détectée par des séquences d'imagerie par résonance magnétique traditionnelles. Cette technique CEST permet ainsi d’obtenir des informations précieuses pour un praticien cherchant à établir un diagnostic et à prendre une décision thérapeutique dans le traitement de pathologies. Dans le domaine clinique, de tels biomarqueurs doivent cependant être produits de manière rapide et précise. La mise en œuvre de la technique CEST nécessite toutefois selon l’état de la technique une puissance et/ou des temps de calculs incompatibles avec les exigences imposées par le domaine clinique notamment. De nombreux effets interagissent avec le transfert de saturation durant une application continue d'impulsions radiofréquences, tels que la saturation directe en eau (DS pour « Direct water Saturation » selon une terminologie anglo-saxonne) ou encore le contraste de transfert d'aimantation (MT pour « Magnetization Transfer » selon une terminologie anglo-saxonne). D’autres phénomènes se combinent auxdits effets et sont inhérents à la mise en œuvre d’une technique d’imagerie par résonance magnétique, tels que des inhomogénéités de champs magnétiques. En effet, le contraste CEST est basé sur l'échange chimique entre les protons labiles et l'eau. Le contraste MT provient quant à lui de l'échange de magnétisation entre les protons de l'eau et les protons dans un environnement solide ou semi-solide. Des corrections sont donc nécessaires pour un obtenir un signal CEST fiable. Par ailleurs, bien qu'en théorie, pour extraire un signal caractérisant la concentration d’une molécule d'intérêt, une saturation en eau ne nécessite que quelques localisations spectrales, par exemple, principalement autour de la fréquence de résonance de l’eau, en pratique une telle extraction est complexe en raison de la présence d'inhomogénéités du champ magnétique statique B0 des appareils d’imagerie par résonance magnétique. Ces inhomogénéités sont partiellement causées par les propriétés tissulaires du patient et entraînent un décalage des échantillons spectraux acquis, comme l’indique la . Pour surmonter ce problème, un échantillonnage dense est effectué autour de la fréquence de résonance de l'eau. Les mesures acquises sont ensuite normalisées par un signal mesuré sans saturation par radiofréquence. Avec une interpolation mathématique appropriée comme l’indique la , le spectre Z peut être obtenu, voxel par voxel, pour une plage de décalages fréquentiels d’intérêt. Les échantillons de ce dernier, après ordonnancement selon des décalages fréquentiels Δω décroissants, décrivent un signal discret de forme ou courbe sensiblement en « V ». Ledit spectre Z est alors « recentré » pour que son minimum Z(Δω)m corresponde à la fréquence de l'eau comme l’indique la . Le spectre Z peut ainsi être décrit sous une forme relative, la fréquence de l’eau étant égale à un décalage fréquentiel Δω nul, les fréquences ω étant exprimées en ppm sous la forme d’un décalage fréquentiel relatif Δω, positif ou négatif au regard de la fréquence de l’eau. Un tel décalage fréquentiel relatif est généralement exprimé en ppm. Une correction du paramètre B0, telle que l’illustre la , peut ainsi être appliquée pour une pluralité de voxels ou volumes élémentaires. Des cartes appelées cartes MB0 intrinsèques, c’est-à-dire des images matricielles, peuvent être calculées puis affichées, en utilisant la position de décalage du minimum Z(Δω)m du spectre Z pour un ensemble de volumes élémentaires. La illustre la détermination d’un premier exemple non limitatif d’un biomarqueur produit à partir d’un spectre Z tel que décrit précédemment en lien avec les figures 3 à 5. Un tel biomarqueur peut correspondre à la différence entre une première aire ZA1, que l’on nomme généralement « Label » selon une terminologie anglo-saxonne, sous la courbe interpolée du spectre Z à partir des échantillons Z(Δω) et une deuxième aire ZA2, que l’on nomme généralement « Reference » selon une terminologie anglo-saxonne, sous ladite courbe, pour des plages de décalages fréquentiels Δω respectives et symétriques au regard du décalage fréquentiel nul. En l’espèce, sur l’exemple de la , l’aire ZA1 est associée aux décalages fréquentiels compris entre +3 ppm et +4 ppm et l’aire ZA2 est associée aux décalages fréquentiels compris entre -3 ppm et -4 ppm. Un tel choix de plages est usuellement associé aux amides. En effet, la fréquence de résonnance des amides correspond généralement en moyenne au décalage fréquentiel Δω=3,5 ppm. En variante, si l’on s’intéresse aux amines dont la fréquence de résonnance correspond généralement en moyenne au décalage fréquentiel Δω=2 ppm, lesdites plages pourront être choisies comme étant comprises respectivement de 1,5 ppm à 2,5 ppm pour l’aire ZA1 et de -1,5 ppm à -2,5 ppm pour l’aire ZA2. Des cartes paramétriques associées à ce biomarqueur peuvent être constituées et affichées éventuellement pour décrire ledit biomarqueur pour une pluralité de volumes élémentaires d’intérêt. La illustre la détermination d’un deuxième exemple non limitatif d’un biomarqueur produit à partir d’un spectre Z tel que décrit précédemment en lien avec les figures 3 à 5. Contrairement à l’exemple précédent selon lequel le biomarqueur est obtenu par soustraction de deux aires ZA1 et ZA2, ledit biomarqueur peut être produit par soustraction de deux valeurs interpolées d’échantillons Z1 et Z2 du spectre Z, lesdits échantillons étant également choisis de sorte que les décalages fréquentiels Δω1 et Δω2 qui leurs sont respectivement associés soient symétriques au regard du décalage fréquentiel Δω=0. A l’instar de l’exemple précédent, si l’on s’intéresse aux amides, comme l’indique la , les décalages fréquentiels Δω1=3,5ppm et Δω2=-3,5ppm seront privilégiés. En variante ou en complément, si l’on s’intéresse aux amides, on pourra choisir des décalages fréquentiels Δω1=2ppm et Δω2=-2ppm. Un tel biomarqueur est nommé « Magnetic Transfer Ratio asymmetry (MTRasym) » selon une terminologie anglo-saxonne. De tels choix peuvent faire l’objet de paramétrages prédéterminés ou d’options renseignées par un utilisateur au moyen d’une interface homme-machine d’entrée appropriée. L'augmentation de la densité des échantillons du spectre Z conduit à une résolution plus élevée des cartes MB0 calculées mais l'acquisition de plusieurs échantillonnages du spectre Z induit un temps de calcul important et n’est pas toujours synonyme de résultats de cartographie du paramètre B0 optimaux. Outre les difficultés induites par les inhomogénéités du champ magnétique statique B0 des appareils d’imagerie par résonance magnétique, un spectre Z est également affecté par les inhomogénéités du champ d'excitation B1. Pour surmonter ce problème, une correction du paramètre B1 est également généralement mise en œuvre pour parvenir à extraire un signal CEST fiable. Pour pallier ces inconvénients, la séquence WASAB1 a été proposée pour permettre une détermination plus précise de la fréquence de résonance des protons d'eau « en vrac ». Comme évoqué précédemment, les inhomogénéités du champ magnétique statique B0 et/ou du champ d’excitation B1 sont inévitables en imagerie par résonance magnétique. Elles provoquent des distorsions dans l'intensité des images acquises. Il est donc nécessaire de mettre en œuvre un procédé de post-traitement d’une séquence WASAB1 pour augmenter la pertinence desdites intensités en déterminer précisément des informations spatiales en lien avec de la distribution des champs magnétiques B0 et B1 des appareils d’imagerie par résonnance magnétique. L'acquisition WASAB1 est une séquence d’acquisition en imagerie par résonance magnétique qui permet de cartographier simultanément le champ magnétique statique B0 et le champ d'excitation B1 sur la base d’oscillations de Rabi induites par une irradiation hors résonance qui peut être utilisée pour corriger les inhomogénéités précédemment évoquées. Un spectre Z, résultant d’une séquence WASAB1, découle d’un échantillonnage selon plusieurs décalages fréquentiels (offsets selon une terminologie anglo-saxonne) autour de la fréquence de résonance de l'eau. Celui-ci est ensuite normalisé en utilisant une image M0 acquise pour des fréquences éloignées de la fréquence de l'eau. Les intensités du spectre Z(Δω) WASAB1 d'un seul voxel peuvent être modélisées sur la base des équations de Bloch comme, par exemple et de manière non limitative comme le suggère la Publication : où Δω sont les décalages fréquentiels relatifs autour de la fréquence de l'eau, c, d, B1 et B0 sont quatre paramètres libres. Les paramètres c et d décrivent la modulation d'amplitude indépendamment de tout décalage fréquentiel. Le paramètre B1 change la périodicité et B0 l'axe de symétrie de la fonction, γ est le rapport gyromagnétique et tp est la durée de l'impulsion. La fonction ou modèle précédent peut être écrit sous une forme simplifiée telle que : Comme l’indique la , en ajustant ces paramètres pour chaque volume élémentaire d’intérêt ayant fait l’objet d’une séquence WASAB1, tels que les volumes élémentaires Va et Vb référencés à titre d’exemples sur la , des cartes d'inhomogénéité MB0 et MrB1 peuvent être produites pour être utilisées à des fins de contrôle qualité pour corriger la séquence d'imagerie acquise. La carte MrB1 est obtenue en divisant la carte MB1 de WASAB1 par l'amplitude nominale du champ d’excitation. Généralement le champ magnétique d’excitation B1 nominal utilisé pour acquérir des séquences WASAB1 est de 3,7 μT. Malgré des temps d’acquisition rapides de séquences WASAB1, les temps de calcul pour mettre en œuvre un post-traitement sont très lents car basés sur un ajustement d’un modèle à quatre paramètres c, d, B0 et B1 dérivés des équations de Bloch. A l’heure actuelle, une seule méthode de post-traitement a été proposée dans la littérature pour déterminer les paramètres B0 et B1 d’une séquence WASAB1, à savoir la méthode exposée dans la Publication. Une telle méthode consiste à minimiser la somme quadratique des erreurs résiduelles entre mesures et modèle à l’aide d’un procédé d’optimisation non linéaire. Un tel post-traitement connu comprend deux étapes. Selon une première étape, l'espace des paramètres c, d, B0 et B1 est recherché de manière exhaustive au moyen d’une table ou grille de consultation à quatre dimensions pour trouver les meilleurs paramètres c, d, B0 et B1, pour chaque voxel ou volume élémentaire d’intérêt. En effet, pour éviter les minima locaux pendant la mise en œuvre de l’étape d'ajustement, compte-tenu de la non-linéarité du modèle du signal, les auteurs ont créé cette table ou grille de recherche avec des spectres Z WASAB1 pour toutes les combinaisons possibles des quatre paramètres d'ajustement avant d'appliquer l'algorithme d'ajustement en tant que tel. La meilleure combinaison desdits paramètres est celle qui produit la plus petite erreur quadratique moyenne (« Means Square Error » ou MSE selon une terminologie anglo-saxonne) entre le spectre Z WASAB1 acquis pour un voxel et un spectre simulé pour un ensemble donné de paramètres sélectionnés dans ladite table ou grille de consultation. La deuxième étape d’un tel procédé de post-traitement selon l’état de l’art est une étape d'ajustement du modèle aux données selon laquelle la meilleure combinaison des paramètres c, d, B0 et B1 précédemment trouvée est utilisée comme condition initiale pour un procédé de minimisation bornée utilisant l'algorithme de Levenberg-Marquardt. Les nombreuses itérations qui doivent être mises en œuvre, voxel par voxel, cumulées à un temps de calcul élevé pour mettre en œuvre l’ajustement du modèle aux données, conduisent à une durée de mise en œuvre du procédé de post-traitement supérieure à une heure. Une telle durée est prohibitive pour une application dans le domaine clinique. En outre, le modèle analytique qui décrit le signal expérimental résultant d’une séquence WASAB1 prédit que pour certains ensembles de paramètres, ledit signal - que l’on pourrait qualifier de « théorique » et dont on ne mesure et connait que la magnitude pour différents décalages fréquentiels constituant ainsi un signal expérimental non polarisé - peut devenir négatif. Comme l’indique la , un tel signal expérimental théorique peut décrire une polarité variée pour un premier volume élémentaire Va, contrairement à ce que traduit les mesures de la magnitude d’un tel signal expérimental théorique en lien avec un deuxième volume élémentaire Vb. Un appareil d’imagerie par résonance magnétique mesurant la magnitude du signal expérimental théorique, le signe ou polarité dudit signal expérimental théorique est donc « perdu » pendant l'acquisition. Le spectre Z étant un signal ou ensemble d’échantillons que l’on peut qualifier de « non polarisés ». Ainsi, comme l’indique le spectre Za, associé à un volume élémentaire Va, des échantillons entre des points de rebroussement C1 et C2 (« cusps » selon une terminologie anglo-saxonne) décrivent des magnitudes d’un signal expérimental théorique de polarité négative. La mesure réalisée par l’appareil applique donc un filtre en valeurs absolues sur ledit signal expérimental théorique pour constituer un signal expérimental non polarisé, d’où l’apparition desdits points de rebroussement C1 et C2 sur ce dernier. Par conséquent, le modèle fait intervenir une valeur absolue pour tenir compte de cet effet de mesure. L'utilisation d'une valeur absolue rend le modèle complexe à estimer, car il introduit des discontinuités (points de rebroussement) dans la fonction que l'on cherche à minimiser. De plus, cela compromet toute mise en œuvre d’optimisations pertinentes qui seraient possibles si le signe dudit signal n’était pas « perdu ». La illustre ainsi les échantillons de mesure de la magnitude de deux signaux expérimentaux théoriques WASAB1 respectivement associés à deux volumes élémentaires distincts Va et Vb. Le spectre Za, associé à un premier volume élémentaire Va, décrit la magnitude d’un signal expérimental théorique à la polarité changeante ou variée, c’est-à-dire tantôt positive, tantôt négative, contrairement au spectre Zb, associé à un deuxième volume d’intérêt Vb, qui décrit des mesures de la magnitude d’un signal expérimental théorique de polarité systématiquement positive. Le spectre Za présente ainsi deux minima, sous la forme des points de rebroussement C1 et C2, entre lesquels la magnitude du signal mesurée correspond à une amplitude négative du signal expérimental théorique associé au premier volume élémentaire Va, contrairement aux échantillons situés de part et d’autre desdits deux points de rebroussement C1 et C2 qui décrivent des magnitudes d’un signal expérimental théorique dont l’amplitude, pour les décalages fréquentiels considérés, est de polarité positive. La présente en outre les cartes paramétriques MB0 et MrB1 précédemment évoquées pour l’ensemble des volumes élémentaires d’une tranche d’un cerveau humain. Nous pouvons constater que le paramètre B0, décrit un champ magnétique statique inhomogène compris entre -0,1 ppm et +0,3 ppm. L’invention permet de répondre à tout ou partie des inconvénients soulevés par les solutions connues ou précédemment mentionnées. Parmi les nombreux avantages apportés par l’invention, nous pouvons mentionner que : - un procédé de post-traitement des échantillons d’un signal expérimental, fruit d’une séquence WASAB1, est agencé pour produire estimation robuste et très rapide, au regard des temps et difficultés de mise en œuvre selon l’état de l’art, des paramètres B0 et B1 décrivant respectivement le champ magnétique statistique et le champ magnétique d’excitation de l’appareil d’imagerie médicale réalisant la mesure échantillonnée d’un signal expérimental théorique ; - la mise en œuvre d’un tel procédé de post-traitement ouvre au domaine clinique notamment la technologie CEST particulièrement prometteuse pour délivrer une aide à la décision et au diagnostic en lien avec de nombreuses pathologies ; - l’estimation desdits paramètres B0 et B1, selon l’invention, s’appuie sur l’exploitation d’un modèle plus simple puisqu’il ne comporte que deux paramètres au lieu de quatre et, également plus stable que celui décrit selon l’état de l’art, réduisant de manière significative les aberrations dans l’estimation desdits paramètres B0 et B1 et, par voie de conséquence, permettant la production de cartes paramétriques mises à la disposition du personnel médical, de grandes qualités et précises ; - la mise en œuvre d’un procédé de post-traitement selon l’invention trouve une application directe et performante dans toute application d’imagerie par résonance magnétique pour laquelle une correction d’inhomogénéités desdits champs B0 et B1 est requise. Pour cela, l’invention concerne tout d’abord un procédé (100) de post-traitement d'un premier ensemble Z d'échantillons Z(Δω) d'un signal expérimental, fruit d'une séquence d'acquisition WASAB1 par un appareil d'imagerie médicale par résonance magnétique, en lien avec un volume élémentaire d'un organe, les valeurs respectives des échantillons Z(Δω) dudit premier ensemble Z décrivant la magnitude d’un signal expérimental théorique selon des décalages fréquentiels Δω par rapport à la fréquence de résonance de l'eau, ledit procédé étant mis en œuvre par une unité de traitement d'un système d'analyse d'imagerie. Pour une mise en œuvre particulièrement efficace et rapide d’un tel procédé et, in fine , produire des estimations précises et robustes, un tel procédé comporte : - une étape pour détecter les échantillons Z(Δω) du premier ensemble Z pour lesquels les polarités des valeurs du signal expérimental théorique dont la magnitude est mesurée pour les décalages fréquentiels Δω considérés sont certaines et construire un deuxième ensemble Y d'échantillons Y(Δω) tels que Y(Δω)=Z(Δω) lorsque la polarité du signal expérimental théorique est positive et Y(Δω)=-Z(Δω) si la polarité du signal expérimental théorique est négative pour lesdits décalages fréquentiels Δω considérés ; - une étape pour ajuster un modèle déterminé au deuxième ensemble d’échantillons Y, les deux paramètres dudit modèle déterminé décrivant respectivement les champs magnétiques statique B0 et d'excitation B1 de l'appareil d'imagerie médicale par résonance magnétique et produire une première estimation desdits paramètres B0 et B1 dudit modèle. Le cas échéant, pour améliorer la pertinence des estimations produites, un procédé selon l’invention peut comporter : - une étape pour estimer la polarité des valeurs du signal expérimental théorique mesuré pour tout ou partie des échantillons Z(Δω) à partir du modèle ainsi ajusté et construire le deuxième ensemble Y d'échantillons comprenant un ensemble d'échantillons Y(Δω) tels que Y(Δω)=Z(Δω) si la polarité du signal expérimental théorique est estimée positive et Y(Δω)=-Z(Δω) si la polarité du signal expérimental théorique est estimée négative ; - une étape pour ajuster le modèle déterminé au deuxième ensemble Y d'échantillons et produire une deuxième estimation desdits paramètres B0 et B1. Selon un premier mode de réalisation, l'étape pour détecter les échantillons Z(Δω) du premier ensemble Z, pour lesquels les polarités des valeurs du signal expérimental théorique dont la magnitude est mesurée sont certaines, peut comporter : - une première sous-étape, pour tout échantillon Z(Δω), depuis le plus grand décalage fréquentiel négatif en valeur absolue vers un tel décalage fréquentiel nul, pour extrapoler la valeur dudit échantillon à partir de la ou des valeurs d'un ou plusieurs échantillons respectivement associés à différents décalages fréquentiels ; - une deuxième sous-étape, pour tout échantillon Z(Δω), depuis le plus grand décalage fréquentiel vers un tel décalage fréquentiel nul, pour extrapoler la valeur dudit échantillon à partir de la ou des valeurs d'un ou plusieurs échantillons respectivement associés à différents décalages fréquentiels. Selon ce premier mode de réalisation, la polarité du signal expérimental théorique pour un décalage fréquentiel Δω peut être déterminée par le signe de valeur extrapolée Z(Δω)' de l'échantillon Z(Δω) pour ledit décalage fréquentiel. Selon ce premier mode de réalisation de l'étape pour détecter les échantillons Z(Δω) du premier ensemble Z pour lesquels les polarités des valeurs du signal expérimental théorique dont la magnitude est mesurée sont certaines : - l'extrapolation de la valeur d'un échantillon Z(Δω) peut consister, pour les première et deuxième sous-étapes, en une évaluation préalable du plus grand écart entre les valeurs respectives de deux échantillons Z(Δω) associés à des décalages fréquentiels successifs et en la soustraction dudit plus grand écart à la valeur de l'échantillon associé au décalage fréquentiel supérieur en valeur absolue à celui associé à l'échantillon Z(Δω) dont la valeur est extrapolée ; - la polarité positive du signal expérimental théorique pour un décalage fréquentiel donné Δω est considérée certaine tant qu'une valeur extrapolée Z(Δω)' d'un échantillon Z(Δω) demeure positive. En variante, l’invention prévoit un deuxième mode de réalisation de l'étape pour détecter les échantillons Z(Δω) du premier ensemble Z pour lesquels les polarités des valeurs du signal expérimental théorique mesuré sont certaines. Selon ce deuxième mode, ladite étape peut comporter : - une première sous-étape pour calculer la dérivée d'un premier signal discret résultant d'un ordonnancement des échantillons dudit premier ensemble Z selon des décalages fréquentiels croissants ; - une deuxième sous-étape pour détecter toute discontinuité de la polarité de ladite dérivée ; - la polarité de la mesure du signal expérimental théorique pour un ou plusieurs décalages fréquentiels successifs Δω étant inversée lors d'une discontinuité de la polarité de ladite dérivée. Quel que soit le mode de réalisation de l’étape pour détecter les échantillons Z(Δω) du premier ensemble Z traduisant des mesures de la magnitude du signal expérimental théorique dont les polarités des valeurs dudit signal expérimental théorique sont certaines, les étapes nécessaires à l'estimation des paramètres B0 et B1 en lien avec un volume élémentaire de l'organe peuvent être itérées pour produire une pluralité d'estimations desdits paramètres B0 et B1 en lien respectivement avec une pluralité de volumes élémentaires dudit organe. Dans ce cas, afin de délivrer une aide précieuse pour calibrer l’appareil d’imagerie ou pour diagnostiquer certaines pathologies humaines ou animales, un procédé conforme à l’invention peut comporter une étape pour produire une image sous la forme d'une carte paramétrique dont les pixels encodent respectivement les valeurs du paramètre B0 ou B1 estimé pour des volumes élémentaires d'intérêts. L’invention prévoit alors qu’un tel procédé puisse comporter une étape pour provoquer un affichage de ladite image produite par une interface homme-machine de sortie du système d'analyse d'imagerie. Selon un deuxième objet, l’invention concerne un système d'analyse d'imagerie comportant une unité de traitement, des moyens pour communiquer avec le monde extérieur et des moyens de mémorisation. Les moyens pour communiquer d’un tel système sont agencés pour recevoir dudit monde extérieur des données expérimentales résultant d'une séquence d'acquisition WASAB1 de données réalisées par échantillonnage spectral d'un volume élémentaire d'un organe et les moyens de mémorisation comportent des instructions de programme dont l'interprétation ou l'exécution par ladite unité de traitement provoque la mise en œuvre d'un procédé de post-traitement d'un premier ensemble Z d'échantillons Z(Δω) d'un signal expérimental fruit d'une séquence d'acquisition WASAB1 conforme à l’invention. Selon un troisième objet, l’invention concerne un produit programme d'ordinateur comportant une ou plusieurs instructions de programme interprétables par l'unité de traitement d'un système d'analyse d'imagerie, ledit programme étant chargeable dans des moyens de mémorisation dudit système, caractérisé en ce que l'interprétation desdites instructions par ladite unité de traitement provoque la mise en œuvre d'un procédé de post-traitement d'un premier ensemble Z d'échantillons Z(Δω) d'un signal expérimental théorique continu fruit d'une séquence d'acquisition WASAB1 conforme à l’invention. Selon un quatrième objet, l’invention concerne un support de mémorisation lisible par un ordinateur comportant les instructions de programme d'un tel produit programme d'ordinateur. D’autres caractéristiques et avantages apparaîtront plus clairement à la lecture de la description qui suit et à l'examen des figures qui l'accompagnent parmi lesquelles : , d’ores-et-déjà décrite, illustre une description simplifiée d’un système d’analyse d'images obtenues par Résonance Magnétique Nucléaire ; , d’ores-et déjà-décrite, illustre une description simplifiée d’une variante d’un système d’analyse d'images obtenues par Résonance Magnétique Nucléaire ; , d’ores-et-déjà décrite, présente un premier exemple d’un spectre Z sous la forme d’un ensemble d’échantillons Z(Δω) ordonnés en lien avec un volume élémentaire d’intérêt d’un organe ; , d’ores-et-déjà décrite, présente le même exemple d’un spectre Z sous la forme d’un ensemble d’échantillons Z(Δω) ordonnés que celui de la , superposé à un signal expérimental continu estimé après ajustement d’un modèle déterminé aux échantillons dudit spectre Z conformément à l’état de l’art ; , d’ores-et-déjà décrite, présente le même exemple d’un spectre Z sous la forme d’un ensemble d’échantillons Z(Δω) ordonnés que celui de la , superposé à un signal expérimental continu estimé après ajustement d’un modèle déterminé aux échantillons dudit spectre Z conformément à l’état de l’art, après correction et réglage de l’appareil d’imagerie médicale ; , d’ores-et-déjà décrite, présente une estimation d’un premier biomarqueur à partir d’échantillons d’un spectre Z conformément à l’état de l’art après correction et réglage de l’appareil d’imagerie ; , d’ores-et-déjà décrite, présente une estimation d’un deuxième biomarqueur à partir d’échantillons d’un spectre Z conformément à l’état de l’art après correction et réglage de l’appareil d’imagerie ; , d’ores-et-déjà décrite, présente plusieurs exemples d’un spectre Z sous la forme d’un ensemble d’échantillons Za(Δω), Zb(Δω) ordonnés pour une pluralité de volumes élémentaires Va, Vb d’un cerveau humain, ainsi que des cartes paramétriques associées aux champs statique B0 et d’excitation B1 d’un appareil d’imagerie médicale ; décrit un premier mode de réalisation d’une étape d’un procédé de post-traitement d’un spectre WASAB1, selon l’invention, visant à détecter les échantillons dudit spectre décrivant des magnitudes d’un signal expérimental théorique dont la polarité est considérée positive aux décalages fréquentiels considérés ; décrit un deuxième mode de réalisation d’une étape d’un procédé de post-traitement d’un spectre WASAB1, selon l’invention, visant à détecter les échantillons dudit spectre décrivant des magnitudes d’un signal expérimental théorique dont la polarité est considérée connue aux décalages fréquentiels considérés ; illustre un exemple d’algorithme fonctionnel d’un procédé de post-traitement d’un spectre WASAB1, selon l’invention. Décrivons à présent, en lien avec la , un exemple d’un procédé de post-traitement 100, selon l’invention, d’un premier ensemble Z d’échantillons Z(Δω) d’un signal expérimental 10, 12, fruit d’une séquence d’acquisition WASAB1 par un appareil d’imagerie médicale par résonance magnétique 1, tel que celui décrit précédemment en liaison avec la ou 2. Un tel procédé 100 comporte un traitement 100a qui peut être mis en œuvre itérativement pour une pluralité de volumes élémentaires d’un organe d’intérêt. Un tel traitement 100a, mis en œuvre pour un voxel ou volume élémentaire, comporte deux d’étapes 110 et 120, éventuellement complétées par deux étapes 130 et 140 comme nous le verrons ultérieurement en lien avec deux modes de réalisation avantageux. Un premier objectif de ce traitement 100a consiste à construire, en une étape 110, un deuxième ensemble d’échantillons Y à partir de l’ensemble des échantillons d’un premier spectre Z découlant une séquence WASAB1. Contrairement audit spectre Z pour lequel les échantillons Z(Δω) peuvent décrire une amplitude d’un signal expérimental théorique 10, 12 dont la polarité est incertaine ou variable, le deuxième ensemble Y d’échantillons comporte des échantillons Y(Δω) décrivant des amplitudes dudit signal expérimental théorique 10, 12 dont la polarité est connue. Nous étudierons successivement une telle construction du deuxième ensemble Y en liaison avec les figures 8 et 9 qui décrivent deux modes de réalisation distincts. Selon un premier mode de réalisation illustré conjointement par les figures 8 et 10, le deuxième ensemble d’échantillons Y peut être constitué sous la forme d’un sous-ensemble du spectre Z pour lequel les échantillons Y(Δω) sont respectivement égaux aux échantillons Z(Δω) lorsque la polarité du signal expérimental théorique dont la magnitude est mesurée est certaine et positive. Pour cela, l’étape 110 comporte des sous-étapes 111 et 112, de sorte que certains échantillons du premier ensemble Z soient écartés et donc non reproduits dans le deuxième ensemble Y si la polarité du signal expérimental théorique dont la magnitude est mesurée, aux décalages fréquentiels considérés Δω, est incertaine, voire vraisemblablement négative. En variante, en lien avec l’exemple de la , lorsque l’étape 110 comporte des sous-étapes 113 et 114, les échantillons Y(Δω) du deuxième ensemble d’échantillons Y sont respectivement égaux aux échantillons Z(Δω) du premier ensemble Z lorsque la polarité du signal expérimental théorique dont la magnitude est mesurée est certainement positive. Lesdits échantillons Y(Δω) sont respectivement égaux aux -Z(Δω) lorsque la polarité dudit signal expérimental théorique dont la magnitude est mesurée est certainement négative. Le deuxième ensemble Y d’échantillons est donc « polarisé » contrairement au spectre Z. Quel que soit le mode de réalisation retenu, le modèle pour estimer les paramètres B0 et B1, lorsque celui-ci est ajusté en une étape 120 voire 140, au deuxième ensemble Y d’échantillons dont on connaît la polarité du signal expérimental théorique dont la magnitude est mesurée aux décalages fréquentiels considérés, n’intègre plus la notion de valeurs absolues, contrairement au modèle selon l’état de l’art. Grâce à l’apport de l’invention, ledit modèle est donc plus simple à ajuster et robuste. Nous verrons également que ce modèle peut être « ramené » de quatre paramètres à deux paramètres, décuplant d’un facteur cent les performances de mise en œuvre d’un procédé de post-traitement 100 conforme à l’invention au regard de l’état de l’art. Etudions en liaison avec les figures 8 et 10, un premier mode de réalisation de la première étape 110 d’un procédé de post-traitement 100 conforme à l’invention. Ladite étape 110 consiste à détecter les échantillons Z(Δω) du premier ensemble Z pour lesquels la polarité des valeurs du signal expérimental théorique dont la magnitude est mesurée est certaine et positive. Ladite étape 110 consiste en outre à construire un deuxième ensemble Y d’échantillons Y(Δω) tels que Y(Δω)=Z(Δω) à partir des seuls échantillons Z(Δω) précédemment détectés comme décrivant des mesures d’un signal expérimental théorique de polarité positive aux décalages fréquentiels Δω considérés. La illustre, dans la partie gauche, un spectre Z dont les échantillons Z(Δω) sont représentés par des points et sont ordonnés selon des décalages fréquentiels croissants, en l’espèce depuis un décalage fréquentiel égal -2 ppm vers un décalage fréquentiel égal à +2 ppm. Ce premier ensemble d’échantillons Z ne comporte que des valeurs positives traduisant des mesures de la magnitude d’un signal expérimental théorique WASAB1 dont on ne connaît pas la polarité a priori. Sur la partie droite de ladite est représenté, selon un même ordonnancement de décalages fréquentiels Δω, un deuxième ensemble Y d’échantillons construit à partir du premier ensemble Z. L’objectif de cette étape 110 est de construire l’ensemble Y de telle manière qu’il ne comporte que des échantillons Y(Δω) identiques aux Z(Δω) du premier ensemble Z, si et seulement si lesdits échantillons Z(Δω) traduisent respectivement des mesures de la magnitude d’un signal expérimental théorique dont la polarité est certaine et positive aux décalages fréquentiels Δω considérés. Les échantillons Z(Δω) du premier ensemble Z, traduisant des mesures d’un signal expérimental théorique dont la polarité est incertaine, voire vraisemblablement négative, aux décalages fréquentiels considérés, sont écartés. Le deuxième ensemble d’échantillons Y est donc, comme l’indique la , un sous-ensemble du premier ensemble Z. Nous pouvons qualifier l’ensemble Y d’ensemble d’échantillons polarisés puisque la polarité du signal expérimental théorique dont la magnitude est mesurée est connue pour les décalages fréquentiels considérés. Pour détecter les échantillons Z(Δω) décrivant des mesures d’un signal expérimental théorique de polarité certainement positive aux décalages fréquentiels Δω considérés, l’étape 110 comporte une première sous-étape 111 consistant, pour tout échantillon Z(Δω), depuis le plus grand décalage fréquentiel négatif Δω, en valeur absolue, vers un tel décalage fréquentiel nul, à extrapoler la valeur dudit échantillon Z(Δω) à partir de la ou des valeurs d’un ou plusieurs échantillons respectivement associés à des décalages fréquentiels négatifs supérieurs en valeurs absolues. A titre d’exemple préféré mais non limitatif, selon la , la valeur extrapolée Z(Δω)’, représentée par un ‘x’ sur la , d’un échantillon Z(Δω), peut consister en une évaluation préalable du plus grand écart, référencé Gm sur la , entre les valeurs respectives de deux échantillons Z(Δω) associés à des décalages fréquentiels successifs. Ladite étape 111 consiste alors à calculer la valeur extrapolée Z(Δω)’ en soustrayant ledit plus grand écart Gm à la valeur de l’échantillon associé au décalage fréquentiel supérieur en valeur absolue à celui associé à l’échantillon Z(Δω) dont la valeur est extrapolée. Les barres verticales, connectées aux valeurs extrapolées et symbolisées par des ‘x’ sur la , illustrent ainsi cette opération d’extrapolation, réalisée de la gauche vers la droite, des échantillons Z(Δω) depuis le décalage fréquentiel Δω proche de -2 ppm vers le décalage fréquentiel proche de 0 ppm. La polarité du signal expérimental théorique dont la magnitude mesurée correspond aux échantillon Z(Δω) est considérée comme certaine et positive tant que les valeurs extrapolée Z(Δω)’ demeurent positives ou tant que le décalage fréquentiel considéré Δω n’est pas nul. Dès qu’une valeur extrapolée Z(Δω)’ est négative, cela signifie que le signal expérimental discret décrit par les échantillons du premier ensemble Z présente un point de rebroussement. L’étape 111 s’interrompt, les échantillons subséquents, c’est-à-dire associés aux décalages fréquentiels Δω supérieurs, sont considérés comme décrivant des mesures d’un signal expérimental théorique dont la polarité est incertaine ou négative. Seuls les échantillons précédents sont détectés comme décrivant des mesures de la magnitude d’un signal expérimental théorique de polarité positive aux décalages fréquentiels considérés. Ces échantillons détectés sont ainsi reproduits dans le deuxième ensemble Y tel que Y(Δω)=Z(Δω). Ainsi, dans l’exemple de la , à l’issue de la mise en œuvre de la sous-étape 111, seuls les seize premiers échantillons sont considérés favorablement pour constituer l’ensemble Y. Le premier échantillon Z(Δω), c’est-à-dire celui qui est associé au plus grand décalage fréquentiel en valeur absolue, est, par hypothèse, détecté favorablement pour initialiser l’étape 111. L’étape 110 comporte en outre une sous-étape 112 décrivant un comportement miroir à celui de la sous-étape 111. Elle consiste à parcourir les échantillons Z(Δω), de la droite vers la gauche sur la , c’est-à-dire depuis le décalage fréquentiel Δω le plus grand vers un décalage fréquentiel nul, pour détecter les échantillons Z(Δω) du premier ensemble Z traduisant des mesures de la magnitude d’un signal expérimental théorique dont la polarité est certaine et positive aux décalages fréquentiels considérés. Elle consiste en outre à enrichir le deuxième ensemble d’échantillons Y à partir de ces seuls échantillons détectés, tels que Y(Δω)=Z(Δω) et à écarter les échantillons Z(Δω) décrivant des mesures de la magnitude d’un signal expérimental théorique de polarité incertaine ou vraisemblablement négative pour les décalages fréquentiels considérés. Pour cela, la sous-étape 112 consiste, pour tout échantillon Z(Δω), depuis le plus grand décalage fréquentiel positif Δω vers un tel décalage fréquentiel nul, à extrapoler la valeur dudit échantillon Z(Δω) à partir de la ou des valeurs d’un ou plusieurs échantillons Z(Δω) respectivement associés à des décalages fréquentiels Δω positifs supérieurs. A titre d’exemple préféré mais non limitatif, selon la , la valeur extrapolée Z(Δω)’, représentée par un ‘x’ sur la , d’un échantillon Z(Δω) peut consister en le calcul de la valeur extrapolée Z(Δω)’ d’un échantillon Z(Δω) en soustrayant le plus grand écart Gm entre les valeurs respectives de deux échantillons associés à deux décalages fréquentiels consécutifs, écart Gm d’ores et déjà évoqué dans le cadre de la sous-étape 111, à la valeur de l’échantillon associé au décalage fréquentiel inférieur à celui associé à l’échantillon Z(Δω) dont la valeur est extrapolée. Les barres verticales, connectées aux valeurs extrapolées et symbolisées par des ‘x’ sur la , illustrent ainsi cette opération d’extrapolation, réalisée de la droite vers la gauche, des échantillons Z(Δω) depuis le décalage fréquentiel Δω proche de +2 ppm vers le décalage fréquentiel Δω proche de 0 ppm. La polarité du signal expérimental théorique dont la magnitude Z(Δω) est mesurée pour les décalages fréquentiels Δω est considéré comme certaine et positive tant que les valeurs extrapolées Z(Δω)’ demeurent positives et tant que le décalage fréquentiel considéré Δω n’est pas nul. Dès qu’une valeur extrapolée Z(Δω)’ est négative, cela signifie que le signal discret, décrit par les échantillons du premier ensemble Z, présente un point de rebroussement. La sous-étape 112 s’interrompt, les échantillons subséquents, c’est-à-dire associés aux décalages fréquentiels Δω inférieurs, sont considérés comme décrivant des mesures d’un signal expérimental théorique dont la polarité est incertaine ou négative. Seuls les échantillons précédents sont détectés comme décrivant des mesures d’un signal expérimental théorique dont la polarité est positive aux décalages fréquentiels considérés. Ces échantillons sont ainsi reproduits dans le deuxième ensemble Y tel que Y(Δω)=Z(Δω). Ainsi, dans l’exemple de la , à l’issue de la mise en œuvre de la sous-étape 112, seuls les quinze premiers échantillons sont considérés favorablement pour enrichir l’ensemble Y, l’échantillon Z(Δω) associé au plus grand décalage fréquentiel étant, par hypothèse, détecté favorablement pour initialiser la sous-étape 112. Les sous-étapes 111 et 112 ont été décrites au travers d’un exemple de mise en œuvre préféré selon lequel l’élaboration de la valeur extrapolée Z(Δω)’, représentée par un ‘x’ sur la , d’un échantillon Z(Δω) peut consister en une évaluation préalable du plus grand écart, référencé Gm sur la , entre les valeurs respectives de deux échantillons Z(Δω) associés à des décalages fréquentiels successifs puis en la soustraction dudit plus grand écart Gm à la valeur de l’échantillon associé au décalage fréquentiel supérieur en valeur absolue à celui associé à l’échantillon Z(Δω) dont la valeur est extrapolée. D’autres techniques d’extrapolation pourraient, en variante ou en complément, être mises en œuvre. Nous pouvons de manière plus générale considérer que la valeur extrapolée Z(Δωi)’ d’un échantillon Z(Δωi), c’est-à-dire le ième échantillon associé au décalage fréquentiel Δωi parmi les n échantillons disponibles, peut s’écrire selon la relation : Z(Δωi)’ = ƒ (Z(Δω1), Z(Δω2, …, Z(Δωn), P) où n est un entier positif, ƒ une fonction quelconque et P un paramètre quelconque, éventuellement pluriel. Ainsi, selon l’exemple précédent, la fonction ƒ vise à calculer l’écart Gm et à réaliser l’opération de soustraction et P=1 détermine le rang de l’échantillon précédent de celui dont la valeur Z(Δωi) est extrapolée. En variante, selon un deuxième exemple de mise en œuvre non limitatif, la fonction ƒ pourrait consister en le calcul de la valeur moyenne de valeurs de P échantillons subséquents à l’échantillon dont on cherche à extrapoler la valeur. D’autres combinaisons du couple (ƒ, P) pourraient en variante être exploitées sans sortir du cadre de la présente invention. La illustre un deuxième mode de réalisation de l’étape 110 d’un procédé de post-traitement 100 conforme à l’invention. Ladite étape 110 consiste à détecter les échantillons Z(Δω) du premier ensemble Z pour lesquels la polarité des valeurs du signal expérimental théorique dont la magnitude est mesurée, est certaine. Ladite étape 110 consiste en outre à construire un deuxième ensemble Y d’échantillons Y(Δω) tels que Y(Δω)=Z(Δω) lorsque les échantillons Z(Δω) décrivent des mesures d’un signal expérimental théorique de polarité positive et tels Y(Δω)=-Z(Δω) à partir des seuls échantillons Z(Δω) décrivant des mesures d’un signal expérimental théorique de polarité négative. La , à l’instar de la précédemment décrite, illustre dans la partie gauche un spectre Z dont les échantillons Z(Δω) sont représentés par des points et sont ordonnés selon des décalages fréquentiels croissants, en l’espèce depuis -2 ppm vers +2 ppm. Ce premier ensemble d’échantillons Z comporte des valeurs positives traduisant des mesures de la magnitude d’un signal expérimental théorique WASAB1 dont on ne connaît pas la polarité a priori. Sur la partie droite de ladite figure est représenté, selon un même ordonnancement des décalages fréquentiels Δω, un deuxième ensemble d’échantillons Y construit à partir du premier ensemble Z. L’objectif de cette étape 110 est de construire l’ensemble Y de telle manière qu’il comporte uniquement des échantillons Y(Δω) identiques à ceux Z(Δω) du premier ensemble Z si et seulement si lesdits échantillons Z(Δω) traduisent respectivement des mesures de la magnitude d’un signal expérimental théorique dont la polarité est certaine et positive aux décalages fréquentiels Δω considérés. Les échantillons Z(Δω) du premier ensemble Z, traduisant des mesures de la magnitude d’un signal expérimental théorique dont la polarité est certainement négative aux décalages fréquentiels considérés, sont tels que Y(Δω) = -Z(Δω). Le deuxième ensemble d’échantillons Y est donc, comme l’indique la , un ensemble d’échantillons que l’on peut qualifier de « polarisés » puisque la polarité du signal expérimental théorique est connue aux décalages fréquentiels considérés. Pour détecter les échantillons Z(Δω) décrivant des mesures de la magnitude d’un signal expérimental théorique de polarité certainement positive ou négative pour les décalages fréquentiels Δω considérés, l’étape 110 comporte une première sous-étape 113 consistant, pour tout échantillon Z(Δω) depuis le plus grand décalage fréquentiel négatif Δω, en le calcul de la dérivée du signal discret Z résultant d’un ordonnancement des échantillons Z(Δω) selon des décalages fréquentiels Δω croissants. Cette dérivée apparaît sous la forme d’une courbe discontinue sur la . Compte tenu de la présence de points de rebroussement C1 et C2 dudit signal discret Z, la dérivée présente deux discontinuités franches dd1 et dd2 de son signe ou polarité. Il en serait tout autre en l’absence de tels points de rebroussements, c’est-à-dire si la polarité du signal expérimental théorique échantillonné était constante. L’étape 110 comporte une deuxième sous-étape 114 pour détecter de telles discontinuités de la polarité de ladite dérivée, la polarité de la mesure du signal expérimental théorique pour un ou plusieurs décalages fréquentiels successifs Δω étant inversée lors d’une discontinuité de la polarité de ladite dérivée. Ainsi, selon l’exemple illustré par la , l’étape 110 permet de considérer respectivement, en amont et en aval, desdits points de rebroussements C1 et C2, c’est-à-dire respectivement en amont et en aval desdites discontinuités dd1 et dd2, que la polarité du signal expérimental théorique est certaine et positive. Réciproquement, ladite polarité du signal expérimental théorique est considérée certaine et négative pour la plage, apparaissant grisée, des décalages fréquentiels comprise entre lesdits points de rebroussement C1 et C2. La forme du signal mesuré étant sensiblement en ‘V’, deux discontinuités du signe de la dérivée sont, au plus, attendues. C’est ainsi que les échantillons Z(Δω) détectés comme décrivant des mesures de la magnitude d’un signal expérimental théorique de polarité positive aux décalages fréquentiels considérés sont reproduits dans le deuxième ensemble Y tel que Y(Δω)=Z(Δω) et que les échantillons Z(Δω) détectés comme décrivant des mesures de la magnitude d’un signal expérimental théorique de polarité négative aux décalages fréquentiels considérés sont reproduits dans le deuxième ensemble Y tel que Y(Δω)=-Z(Δω). Ainsi, dans l’exemple de la , à l’issue de la mise en œuvre de la sous-étape 114, l’ensemble des échantillons du premier ensemble Z est considéré favorablement pour constituer le deuxième ensemble Y d’échantillons polarisés. Par la mise en œuvre de l’étape 110, la polarité du signal expérimental théorique dont découle le deuxième ensemble Y d’échantillons Y(Δω) est connue à chaque décalage fréquentiel considéré. Contrairement à l’état de l’art pour lequel le modèle des intensités d’un spectre Z WASAB1, sur la base des équations de Bloch, s’exprime, sous une forme simplifiée, de telle sorte qu’un échantillon dudit spectre Z(Δω) est fonction, en valeur absolue, de quatre paramètres libres c, d, B1 et B0, l’exploitation du deuxième ensemble d’échantillons « polarisé » Y, en lieu et place du spectre Z permet d’évacuer le problème de complexité lié à la valeur absolue puisque la polarité du signal expérimental théorique est connue pour chaque échantillon Y(Δω). Ainsi, le modèle des intensités du spectre Y peut être écrit sous une forme simplifiée telle que Y(Δω)=c-d.ƒ(B0, B1, Δω). Les paramètres de ce modèle peuvent avantageusement s’obtenir en minimisant la somme de l’erreur quadratique : pour laquelle : est la somme des erreurs quadratiques. En supposant les paramètres B0 et B1 constants, le modèle devient linéaire en les paramètres c et d. Ces derniers peuvent être remplacés par leur estimation au minimum des moindres carrés qui s’exprime analytiquement en fonction de B0 et B1. En remplaçant ces paramètres c et d par leurs estimations en fonction de B0 et B1, le modèle n’a plus que deux paramètres explicites : B0 et B1. Les deux paramètres linéaires ont été rendus implicites. Selon l’invention, le modèle n’est plus un modèle à quatre paramètres mais à deux paramètres, ce qui rend l’estimation de B0 et B1 plus rapide de deux ordres de grandeurs, plus robuste et moins sensible aux minima locaux. En effet, le modèle peut être écrit, pour T volumes élémentaires, sous une forme matricielle telle que Y = M(B0, B1)β, où Y=[Y(Δω1), …,Y(Δωn)]^T décrit un vecteur contenant les n échantillons Y(Δω), par volume élémentaire, pour les décalages fréquentiels Δω1 à Δωn, M est une matrice de taille 2n contenant les régresseurs ou variables du modèle soit : En outre, β=[c,d]^T est un vecteur de taille deux, pour chaque volume élémentaire, contenant les paramètres linéaires c et d. La matrice M dépend explicitement des paramètres non linéaires du modèle, en l’espèce B0 et B1. Ainsi, la somme des erreurs quadratiques, précédemment exprimée, peut s’écrire à présent sous la forme : Pour B0 et B1 constants, l’estimation optimale de β qui minimise la somme des erreurs quadratiques peut être calculée analytiquement telle que : Si l’on remplace cette estimation de β dans l’expression précédente on obtient : Les paramètres c et d peuvent être estimés avec : où : De cette manière, un procédé de post-traitement 100 selon l’invention comporte une étape 120 pour ajuster un modèle déterminé au deuxième ensemble Y d’échantillons polarisés, ensemble préalablement produit à l’étape 110, les deux paramètres dudit modèle déterminé décrivant respectivement les champs magnétiques statique B0 et d’excitation B1 de l’appareil d’imagerie médicale par résonance magnétique 1 tel que décrit à titre d’exemple par la . Ladite étape 120 produit ainsi une première estimation desdits paramètres B0 et B1 dudit modèle. A titre d’exemple préféré, bien que ce mode de réalisation avantageux soit optionnel, l’invention prévoit qu’un procédé 100 comporte deux étapes 130 et 140 visant, conjointement à produire une deuxième estimation desdits paramètres B0 et B1, après enrichissement éventuel du deuxième ensemble Y d’échantillons polarisés. En effet, puisque l’étape 120 a délivré une première estimation des paramètres du modèle, il est possible de déterminer avec précision les échantillons du spectre Z correspondant à des mesures de la magnitude du signal expérimental théorique alors que la polarité de celui est négative pour les décalages fréquentiels considérés. Cela peut être particulièrement avantageux d’enrichir le deuxième ensemble Y d’échantillons, notamment si ce dernier résulte d’une mise en œuvre du premier exemple de réalisation de l’étape 110, c’est-à-dire lorsque ladite étape 110 comporte les sous-étapes 111 et 112. Dans ce cas, ledit deuxième ensemble d’échantillons Y est un sous-ensemble du premier ensemble Z dont les échantillons Y(Δω) sont égaux respectivement aux échantillon Z(Δω) lorsque la polarité du signal expérimental théorique est positive, les autres échantillons de Z étant écartés pour produire l’ensemble Y. Dans ce cas, l’invention prévoit que le procédé 100 puisse comporter une étape optionnelle 130 pour déterminer la polarité des valeurs du signal expérimental théorique pour tout ou partie des échantillons Z(Δω) à partir du modèle ainsi ajusté à l’étape 120 et construire de nouveau le deuxième ensemble Y d’échantillons tels que Y(Δω)=Z(Δω) si la polarité la polarité du signal expérimental théorique est estimée positive et Y(Δω)=-Z(Δω) si la polarité la polarité (Δω) du signal expérimental théorique est estimée négative. Grâce à la mise en œuvre successive des étapes 111, 112, 120 et 130, le deuxième ensemble d’échantillons Y est plus riche que celui obtenu après la seule mise en œuvre de l’étape 110, à l’image de celui obtenu lorsque ladite étape 110 comporte les sous-étapes 113 et 114 précédemment décrites. Un procédé 100 peut alors comporter une nouvelle itération de l’étape 120, représentée par une étape 140 sur la , subséquente à l’étape 130, pour ajuster le modèle à l’ensemble Y d’échantillons polarisés ainsi enrichi. L’étape 140 produit ainsi une nouvelle estimation des paramètres du modèle en utilisant davantage d’échantillons sans que soit pénalisée de manière significative la mise en œuvre du procédé 100 compte tenu du gain découlant de l’exploitation d’un modèle simplifié au regard de celui selon l’état de l’art. Les paramètres B0 et B1 étant notamment estimés de manière précise et rapide à l’issue d’une instance du traitement 100a pour un volume élémentaire d’intérêt, l’invention prévoit qu’un procédé 100 conforme à l’invention comporte autant d’itérations dudit traitement 100a pour produire des estimations notamment des paramètres B0 et B1 pour une pluralité de volumes d’intérêts. Un procédé 100 peut alors comporter une étape 150 pour produire une image sous la forme d'une carte paramétrique, telle que la carte MB0 ou MrB1 précédemment décrite en lien avec la , dont les pixels encodent respectivement les valeurs du paramètre B0 ou B1 estimés pour des volumes élémentaires d’intérêts. Cette carte ou ces cartes paramétriques peuvent en outre être affichées par une interface homme-machine de sortie adaptée, telle que l’interface 5 décrite en lien avec les figures 1 et 2, pour qu’un utilisateur 6 d’un système d’imagerie S puisse prendre connaissance des estimations desdits paramètres B0 ou B1, voire encore de tout autre biomarqueur que l’on peut produire à partir desdits paramètres. La mise en œuvre d’un procédé de post-traitement d’un premier ensemble Z d’échantillons Z(Δω) décrivant la magnitude d’un signal expérimental théorique fruit d’une séquence d’acquisition WASAB1 par un appareil d’imagerie médicale par résonance magnétique 1 devient parfaitement exploitable dans le domaine clinique, ouvrant ainsi la porte audit domaine clinique à la technique CEST notamment, pour quantifier un ou plusieurs biomarqueurs en lien avec une espèce chimique d’intérêt d’un organe humain ou animal pour, in fine, caractériser différentes lésions aux propriétés métaboliques altérées. L’invention ne saurait être réduite à cette seule exploitation. Elle trouve ainsi une place de choix dans toute application d’imagerie par résonance magnétique pour laquelle une correction d’inhomogénéités de champs statique B0 et d’excitation B1 est requise. Procédé (100) de post-traitement d’un premier ensemble Z d’échantillons Z(Δω) d’un signal expérimental (10, 12) fruit d’une séquence d’acquisition WASAB1 par un appareil d’imagerie médicale par résonance magnétique (1), en lien avec un volume élémentaire d’un organe, les valeurs respectives des échantillons Z(Δω) dudit premier ensemble Z décrivant la magnitude d’un signal expérimental théorique selon des décalages fréquentiels Δω par rapport à la fréquence de résonance de l’eau, ledit procédé étant mis en œuvre par une unité de traitement (4) d’un système d’analyse d’imagerie, caractérisé en ce qu’il comporte : - une étape (110) pour détecter les échantillons Z(Δω) du premier ensemble Z pour lesquels les polarités des valeurs du signal expérimental théorique (10, 12) sont certaines, pour les décalages fréquentiels Δω considérés, et construire un deuxième ensemble Y d’échantillons Y(Δω) tels que Y(Δω)=Z(Δω) si la polarité du signal expérimental théorique est positive et Y(Δω)=-Z(Δω) si la polarité dudit signal expérimental théorique (10, 12) est négative pour lesdits décalages fréquentiels Δω considérés ; - une étape (120) pour ajuster un modèle déterminé au deuxième ensemble Y, les deux paramètres dudit modèle déterminé décrivant respectivement les champs magnétiques statique B0 et d’excitation B1 de l’appareil d’imagerie médicale par résonance magnétique (1) et produire une première estimation desdits paramètres B0 et B1 dudit modèle. Procédé selon la revendication précédente comportant : - une étape (130) pour estimer la polarité des valeurs du signal expérimental théorique pour tout ou partie des échantillons Z(Δω) à partir du modèle ainsi ajusté et construire le deuxième ensemble d’échantillons Y comprenant un ensemble d’échantillons Y(Δω) tels que Y(Δω)=Z(Δω) si la polarité du signal expérimental théorique est estimée positive et Y(Δω)=-Z(Δω) si la polarité dudit signal expérimental théorique est estimée négative ; - une étape (140) pour ajuster le modèle déterminé au deuxième ensemble Y d’échantillons et produire une deuxième estimation desdits paramètres B0 et B1. Procédé selon l’une quelconque des revendications précédentes, pour lequel l’étape (110), pour détecter les échantillons Z(Δω) du premier ensemble Z pour lesquels les polarités des valeurs du signal expérimental théorique sont certaines, comporte : - une première sous-étape (111), pour tout échantillon Z(Δω), depuis le plus grand décalage fréquentiel négatif (Δω) en valeur absolue vers un tel décalage fréquentiel nul, pour extrapoler la valeur dudit échantillon à partir de la ou des valeurs d’un ou plusieurs échantillons respectivement associés à différents décalages fréquentiels ; - une deuxième sous-étape (112), pour tout échantillon Z(Δω), depuis le plus grand décalage fréquentiel (Δω) vers un tel décalage fréquentiel nul, pour extrapoler la valeur dudit échantillon à partir de la ou des valeurs d’un ou plusieurs échantillons respectivement associés à différents décalages fréquentiels ; la polarité du signal expérimental théorique pour tout décalage fréquentiel Δω est déterminée par le signe de valeur extrapolée Z(Δω)’ de l’échantillon Z(Δω) pour ledit décalage fréquentiel. Procédé selon la revendication précédente, pour lequel : - l’extrapolation de la valeur d’un échantillon Z(Δω) consiste, pour les première (111) et deuxième (112) sous-étapes, en une évaluation préalable du plus grand écart (Gm) entre les valeurs respectives de deux échantillons Z(Δω) associés à des décalages fréquentiels successifs et en la soustraction dudit plus grand écart (Gm) à la valeur de l’échantillon associé au décalage fréquentiel supérieur en valeur absolue à celui associé à l’échantillon Z(Δω) dont la valeur est extrapolée ; - la polarité positive du signal expérimental théorique pour un décalage fréquentiel donné Δω est considérée certaine tant qu’une valeur extrapolée Z(Δω)’ d’un échantillon Z(Δω) demeure positive. Procédé selon l’une quelconque des revendications 1 à 2, pour lequel l’étape (110), pour détecter les échantillons Z(Δω) du premier ensemble Z pour lesquels les polarités des valeurs du signal expérimental théorique mesurées sont certaines, comporte : - une première sous-étape (113) pour calculer la dérivée d’un premier signal discret résultant d’un ordonnancement des échantillons dudit premier ensemble Z selon des décalages fréquentiels croissants ; - une deuxième sous-étape (114) pour détecter tout discontinuité de la polarité de ladite dérivée ; la polarité du signal expérimental théorique pour un ou plusieurs décalages fréquentiels successifs Δω étant inversée lors d’une discontinuité de la polarité de ladite dérivée. Procédé selon l’une quelconque des revendications précédentes, pour lequel les étapes (100a, 110, 120, 130, 140) nécessaires à l’estimation des paramètres B0 et B1 en lien avec un volume élémentaire de l’organe sont itérées pour produire une pluralité d’estimations desdits paramètres B0 et B1 en lien respectivement avec une pluralité de volumes élémentaires dudit organe. Procédé selon la revendication précédente, comportant une étape (160) pour produire une image sous la forme d’une carte paramétrique (MB0, MrB1) dont les pixels encodent respectivement les valeurs du paramètre B0 ou B1 estimé pour des volumes élémentaires d’intérêts. Procédé selon la revendication précédente comportant une étape (160) pour provoquer un affichage de ladite image (MB0, MrB1) produite par une interface homme-machine de sortie (5) du système d’analyse d’imagerie. Système d'analyse d'imagerie (S) comportant une unité de traitement (4), des moyens pour communiquer avec le monde extérieur et des moyens de mémorisation, caractérisé en ce que : - les moyens pour communiquer sont agencés pour recevoir dudit monde extérieur des données expérimentales résultant d’une séquence d’acquisition WASAB1 de données réalisées par échantillonnage spectral d'un volume élémentaire d'un organe ; - les moyens de mémorisation comportent des instructions de programme dont l'interprétation par ladite unité de traitement provoque la mise en œuvre d'un procédé (100) selon l'une quelconque des revendications précédentes. Produit programme d’ordinateur comportant une ou plusieurs instructions interprétables par l’unité de traitement (4) d’un système d'analyse d'imagerie (S) conforme à la revendication précédente, ledit programme d’ordinateur étant chargeable dans des moyens de mémorisation dudit système, caractérisé en ce que l’interprétation desdites instructions de programme par ladite unité de traitement provoque la mise en œuvre d’un procédé (100) selon l’une quelconque des revendications 1 à 8. Support de mémorisation lisible par un ordinateur comportant les instructions de programme d’un produit programme d’ordinateur selon la revendication précédente.