Brevet FR3157630A1 — publié le 20250627

L’invention concerne un procédé (100) de défloutage d’une image médicale expérimentale (Mn) fruit d’une séquence d’acquisition par un appareil d’imagerie médicale(1) par : calcul des transformées de Fourier (TFMn, TFM0) respectives de ladite image expérimentale (Mn) et d’une image de référence (M0) de même résolution ; détermination des modules (YMn, YM0) respectifs des transformées de Fourier de l’image expérimentale et de l’image de référence ; obtention d’un masque de défloutage (D) par le calcul du ratio des modules (YMn, YM0) ; application dudit ratio (D) à la transformée de Fourier de l’image expérimentale (TFMn) par multiplication, en résultant une transformée de Fourier corrigée (TFMnc) ; obtention de l’image expérimentale défloutée (Mnc) par application de la transformée de Fourier inverse sur ladite transformée de Fourier corrigée (TFMnc). Figure pour l’abrégé : Fig 4 Procédé de défloutage d’une image résultant d’une séquence d’acquisition par un appareil d’imagerie médicale La présente invention a trait, de manière générale, au traitement d’image dans le domaine de l’imagerie médicale. L’invention décrit plus particulièrement un procédé de correction d’un flou ou procédé de défloutage, autrement nommé « deblurring » selon une terminologie anglo-saxonne, d’une image médicale résultant d’une séquence d’acquisition via un appareil d’imagerie médicale. En imagerie médicale, une séquence d’acquisition permet de recueillir et stocker des données brutes sur l’anatomie et/ou sur le fonctionnement d’un organe permettant de reconstruire une représentation visuelle tridimensionnelle dudit organe, en d’autres termes une image dudit organe. De telles acquisitions peuvent induire, dans le domaine fréquentiel, une perte ou un filtrage de certaines fréquences, hautes et/ou basses, induisant un flou sur ladite image. Bien que l’invention puisse s’appliquer à toutes techniques ou modalités basées sur l’imagerie et dont une image résultante aurait subi l’application virtuelle d’un filtre atténuant fortement ou supprimant certaines fréquences dans le domaine fréquentiel, nous distinguerons principalement dans la suite de la présente demande, à titre illustratif, l’imagerie par résonance magnétique, plus connue sous le sigle IRM ou l’acronyme anglo-saxon MRI pour « Magnetic Resonance Imaging » , et notamment l’IRM de perfusion et plus précisément la perfusion par marquage de spins artériels, connue sous l’acronyme ASL pour perfusion par « arterial spin labeling » , selon une terminologie anglo-saxonne. 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 par Résonance Magnétique Nucléaire, tel que l’appareil 1 d’un système d’analyse d’imagerie médicale SAIM 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. L’acquisition IRM se fait par échantillonnage d’un volume parallélépipédique dans un plan de coupe donné. La FIG. 3 illustre les plans de coupe les plus couramment utilisés en imagerie médicale : Vue A : projection de face ou en coupe coronale (ou frontale), la droite du patient est présentée à gauche sur l’image : comme si l’observateur regardait simplement le patient de face ; Vue B : sur une vue axiale (ou transverse), la droite du patient est toujours située à gauche sur l’image : comme si l’observateur regardait le patient – allongé sur le dos – depuis ses pieds ; Vue C : sur une vue sagittale ou de profil, l’image est présentée avec la partie antérieure du patient à gauche, comme si l’observateur regardait le patient de profil gauche. Ledit appareil 1 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 1 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 d’intérêt, que l’on nomme communément un voxel, du volume imagé. Comme l’indique la FIG. 1 , un appareil d'imagerie par Résonance Magnétique Nucléaire 1 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, 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 structurelles 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 dudit système d’analyse d’imagerie. On entend par « moyens de mémorisation » toute mémoire informatique volatile ou, avantageusement, non volatile. Une mémoire non volatile est une mémoire informatique dont la technologie permet de conserver ses données en l’absence d’une alimentation en énergie électrique. Elle peut contenir des données résultant de saisies, de calculs, de mesures et/ou des instructions de programmes. Les principales mémoires non volatiles actuellement disponibles sont de type inscriptible électriquement, telles que l’EPROM (« Erasable Programmable Read-Only Memory », selon une terminologie anglo-saxonne) ou encore inscriptibles et effaçables électriquement, telle que l’EEPROM (« Electrically-Erasable Programmable Read-Only Memory »), flash, SSD (« Solid-State Drive », selon une terminologie anglo-saxonne), etc. Les mémoires non volatiles se distinguent des mémoires dites « volatiles » dont les données sont perdues en l’absence d’une alimentation électrique. Les principales mémoires volatiles actuellement disponibles sont de type RAM (« Random Access Memory » selon une terminologie anglo-saxonne ou encore nommées « mémoires vives »), DRAM (mémoire vive dynamique, nécessitant une réactualisation régulière), SRAM (mémoire vive statique nécessitant une telle réactualisation lors d’une sous-alimentation électrique), DPRAM ou VRAM (particulièrement adaptées à la vidéo), etc. Une « mémoire de données », dans la suite du document, peut être volatile ou non volatile selon l’application visée. 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 ou sortir, in fine, un rendu, par exemple graphique et/ou sonore, d’une estimation ou d’une quantification d’un biomarqueur ou d’un paramètre pharmacocinétique QI é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 FIG. 2 , 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. Plus spécifiquement, l’IRM de perfusion, quant à elle, donne accès à des informations sur la microcirculation capillaire des tissus par analyse de leur microcirculation. Les paramètres quantitatifs essentiels évalués sont les volumes sanguins et des données temporelles (telles que le temps de transit, le temps jusqu’au pic de contraste…). L’objectif final de l’IRM de perfusion est de mesurer, ou d’estimer, le débit sanguin qui irrigue l’organe exploré, exprimé en millilitre/minute/100 grammes de tissu, correspondant à la perfusion tissulaire. La différenciation des tissus perfusés et non perfusés repose sur l’utilisation d’un marqueur intravasculaire : exogène dans un premier cas d’application : produit de contraste injecté ; endogène dans un deuxième cas d’application : par marquage des noyaux d’hydrogène de l’eau de la circulation sanguine. À ce titre, l’imagerie de perfusion ASL se place dans ce deuxième cas d’application pour mesurer la perfusion cérébrale, c’est-à-dire l’apport de sang aux tissus cérébraux, sans utilisation d’agents de contraste ni d’éventuels effets secondaires pharmacologiques. Cette technique est ainsi non invasive ce qui rend son acquisition plus rapide et ce qui permet de la répéter. L’ASL permet de mettre en évidence de manière entièrement non invasive des anomalies de la perfusion dans de nombreuses pathologies, en particulier encéphaliques, ce qui en fait une techniques particulièrement intéressante d’un point de vue clinique. Il s’agit d’une technique différentielle au cours de laquelle deux acquisitions sont réalisées : une acquisition avec marquage des protons artériels et une acquisition contrôle. L’imagerie de perfusion ASL quantifie donc la perfusion dans le cerveau en marquant les noyaux d’hydrogène de l’eau intravasculaire artérielle par une ou des impulsions radiofréquences. Ces impulsions servent à inverser l’aimantation longitudinale du sang artériel, lorsque ce dernier circule à travers un seul plan (ou une plaque de tissu) juste en amont de la région d’intérêt dont l’aimantation mesurable et le temps de relaxation T1 des protons (temps de retour à l’équilibre des protons) seront ainsi modifiés. Ce processus est appelé « marquage » et consiste à convertir l’eau présente dans le sang circulant dans les artères et à travers le cou en un traceur endogène. Ainsi, après un délai permettant au sang de circuler du lieu de marquage jusqu’à son lieu final d’échange avec le cerveau, nommé temps de transit artériel, autrement connu sous l’acronyme anglo-saxon ATT pour « Arterial Transit Time », des images sont acquises et contiennent des signaux provenant à la fois de l’eau « marquée » et de l’eau issue des tissus statiques. Il est également associé à ce processus un délai dit de post-marquage, autrement connu sous l’acronyme anglo-saxon PLD pour « Post Labeling Delay », définissant le délai entre le marquage et l’acquisition d’une image. La deuxième acquisition contrôle, dite image « témoin », est réalisée sans marquage artériel préalable : les protons artériels de la région d’intérêt sont alors à l’équilibre, complètement relaxés. La soustraction de l’acquisition de marquage et de l’acquisition contrôle annule le signal des tissus statiques et permet d’obtenir une image pondérée en fonction de la perfusion. Il est à noter que la différence de signal est de l’ordre de quelques pour cent, d’où un rapport signal/bruit faible. Il est donc nécessaire d’acquérir plusieurs dizaines de couples images marquées–images « témoin », puis de les moyenner pour obtenir un rapport signal/bruit satisfaisant. Ainsi, l’acquisition des données par technique ASL nécessite deux principales étapes : le marquage des protons du sang circulant et l’acquisition des images. Concernant le marquage des protons, il en existe trois types principaux : le marquage continu (connu sous l’acronyme anglo-saxon CASL), le marquage est réalisé de façon continue au niveau des vaisseaux du cou, cette technique permet d’obtenir un bon rapport signal/bruit mais génère un dépôt d’énergie élevé dans les tissus ; le marquage pulsé (connu sous l’acronyme anglo-saxon PASL) utilise de brèves impulsions de radiofréquences limitant le dépôt d’énergie dans les tissus ; et le marquage pseudo-continu (connu sous l’acronyme anglo-saxon pCASL) utilise des trains d’impulsions de radiofréquence de très courte durée. Concernant l’acquisition des images, comme précisé précédemment il est nécessaire d’acquérir de nombreux couples d’images marquées–images « témoin ». En raison d’un rapport signal/bruit satisfaisant et d’une rapidité d’acquisition limitant les artefacts de mouvement, la technique « echo-planar imaging » selon une terminologie anglo-saxonne et connue sous l’acronyme EPI, est principalement utilisée. Toutefois, afin d’améliorer la qualité des images, des séquences 3D ont été développées. Il est possible d’utiliser des séquences 3D « single-shot », ultra-rapides, combinant une acquisition en écho de gradient et spin-écho (connue sous l’acronyme anglo-saxon GRASE). L’utilisation de telles séquences facilite la suppression du signal des tissus statiques. Néanmoins, en raison d’un temps d’écho, représentant l’intervalle de temps entre l’excitation des protons et la survenue du signal, trop long en rapport à la diminution du temps de relaxation T2, soit le temps de déphasage des protons, les images ASL peuvent être affectées par un flou important, visible à la fois dans le plan et à travers le plan, phénomène d’autant plus marqué avec la séquence GRASE, phénomène mentionné, par exemple dans l’article scientifique “Reducing blurring artifacts in 3D-GRASE ASL by integrating new acquisition and analysis strategies” publié dans Proceedings of the International Society for Magnetic Resonance in Medicine volume 22 (2014), page 2704. Un tel flou agit comme un filtre passe-bas dans le domaine fréquentiel, laissant passer les valeurs basses fréquences et coupant toutes les valeurs hautes fréquences. Ainsi, la distinction entre les différents tissus cérébraux, tels que la substance blanche, la substance grise, une tumeur et le liquide céphalo-rachidien (LCR) devient difficile. Le flou est particulièrement problématique à la frontière entre la substance grise et la substance blanche qui présentent des niveaux de perfusion différents (un long temps de transit dans la substance blanche affectant de manière significative le marquage) et pour lesquels les signaux de la substance blanche sont affectés par les signaux de la substance grise. Afin de remédier à cette problématique de flou sur une image expérimentale médicale, il est connu d’avoir recours à des méthodes de suppression du flou. De telles méthodes effectuent généralement une déconvolution dans le domaine image après avoir estimé le noyau de déconvolution dans le domaine fréquentiel, nous pouvons notamment citer les travaux de Chappell et. Al avec pour référence l’article scientifique “Reducing blurring artifacts in 3D-GRASE ASL by integrating new acquisition and analysis strategies” publié dans Proceedings of the International Society for Magnetic Resonance in Medicine volume 22 (2014), page 2704. De plus, toute méthode de suppression du flou d'image bidimensionnelle pourrait être potentiellement appliquée à ladite image expérimentale, en appliquant ladite méthode à chaque tranche de l’image indépendamment. À ce titre, les travaux de Wang et Tao présentent une série de techniques de défloutage dans le contexte du défloutage naturel d'images au sein de l’article scientifique « Recent progress in image deblurring » publié dans arXiv : 1409.6838. Toutefois, de telles méthodes de suppression du flou se concentrent sur la correction du flou traversant uniquement le plan. Ainsi, le développement d’un procédé de suppression et/ou correction du flou tridimensionnel (correction simultanément du flou dans le plan et à travers le plan) d’une image résultant d’une séquence d’acquisition via un appareil d’imagerie médicale trouve tout son intérêt. À ce titre, un tel procédé permet de supprimer les artefacts flous des images afin d’en améliorer la lisibilité et d’en faciliter l’interprétation en vue d’exploiter de telles images dans le domaine clinique 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. A cette fin, l’invention prévoit un procédé de défloutage d’une image médicale expérimentale fruit d’une séquence d’acquisition par un appareil d’imagerie décrivant une pluralité de volumes élémentaires d’intérêt d’un organe d’un patient, nommés « voxels », ledit procédé étant mis en œuvre par une unité de traitement d’un système d’analyse d’imagerie médicale. Un tel procédé comporte : une étape de calcul de la transformée de Fourier de ladite image expérimentale ; une étape de calcul de la transformée de Fourier d’une image de référence de même résolution que l’image expérimentale ; une étape de détermination d’un premier module de la transformée de Fourier de l’image expérimentale et d’un deuxième module de la transformée de Fourier de l’image de référence ; une étape d’obtention d’un masque de défloutage consistant à calculer le ratio entre les premier et deuxième modules ; une étape d’application dudit ratio à la transformée de Fourier de l’image expérimentale par multiplication pour obtenir une transformée de Fourier corrigée ; une étape d’obtention d’une image expérimentale défloutée consistant à appliquer la transformée de Fourier inverse sur ladite transformée de Fourier corrigée. Dans un mode particulier de l’invention, ladite image de référence peut émaner d’un atlas, d’un patient témoin ou d’un référentiel anatomique artificiel. Dans un mode privilégié de l’invention, ledit procédé peut comporter une étape préalable aux étapes de calcul de transformée de Fourier consistant en une étape d’acquisition d’une image structurelle dudit patient de même résolution que ladite image expérimentale. Ladite image structurelle correspond alors à ladite image de référence. Préférentiellement, afin d’assurer la correspondance voxel par voxel de l’image de référence par rapport à l’image expérimentale, le procédé selon l’invention peut comporter une étape préalable aux étapes de calcul de transformée de Fourier consistant en une étape de recalage de l’image expérimentale par rapport à l’image de référence. Selon un mode de réalisation avantageux, ledit procédé peut comporter une étape subsidiaire préalable à l’étape d’obtention de l’image expérimentale défloutée et subséquente à l’étape d’obtention d’un masque de défloutage consistant à appliquer un filtre de lissage sur ledit masque de défloutage. Selon le mode de réalisation précédent, le filtre de lissage peut être un filtre gaussien. Préférentiellement, les transformées de Fourier de l’image expérimentale et de l’image de référence sont respectivement calculées avec l’algorithme de transformée de Fourier rapide. Selon un mode de réalisation préféré, la séquence d’acquisition de l’image expérimentale est une séquence d’acquisition de perfusion par un appareil d’imagerie médicale par résonance magnétique durant laquelle un traceur endogène circule au sein de ladite pluralité de volumes élémentaires d’intérêt. Selon un deuxième objet, l’invention concerne un produit programme d’ordinateur comportant une ou plusieurs instructions de programme exécutables par l’unité de traitement d’un ordinateur, lesdites instructions de programme étant chargeables dans une mémoire non volatile dudit ordinateur et dont l’exécution par ladite unité de traitement provoque la mise en œuvre d’un procédé de post-traitement selon l’invention. Selon un troisième objet, l’invention concerne un support de mémorisation lisible par un ordinateur comportant les instructions d’un tel produit programme d’ordinateur. Enfin, l’invention concerne en outre un système d’analyse d’imagerie médicale comportant une unité de traitement agencée pour déflouter une image médicale expérimentale fruit d’une séquence d’acquisition par un appareil d’imagerie médicale en lien avec une pluralité de volumes élémentaires d’intérêt, nommés voxels, ladite unité de traitement étant configurée pour : calculer la transformée de Fourier de ladite image expérimentale ; calculer la transformée de Fourier d’une image de référence de même résolution que l’image expérimentale ; déterminer un premier module de la transformée de Fourier de l’image expérimentale et un deuxième module de la transformée de Fourier de l’image de référence ; obtenir un masque de défloutage en calculant le ratio entre les premier et deuxième modules ; appliquer ledit ratio à la transformée de Fourier de l’image expérimentale par multiplication pour obtenir une transformée de Fourier corrigée ; obtenir l’image expérimentale défloutée en appliquant la transformée de Fourier inverse sur ladite transformée de Fourier corrigée. L’invention sera mieux comprise et d’autres caractéristiques et avantages de celle-ci apparaîtront à la lecture de la description suivante de modes de réalisation particuliers de l’invention, donnés à titre d’exemples illustratifs et non limitatifs, et faisant référence aux dessins annexés, parmi lesquels : FIG. 1 , 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, FIG. 2 , 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, FIG. 3 d’ores et déjà décrite, illustre les plans de coupe les plus couramment utilisés en imagerie médicale, FIG. 4 illustre un exemple d’algorithme fonctionnel d’un procédé de défloutage selon l’invention, FIG. 5 illustre des variantes de réalisation d’un algorithme fonctionnel d’un procédé de défloutage selon l’invention, FIG. 6 illustre un gain découlant de la mise en œuvre d’un procédé de défloutage 100 conforme à l’invention, FIG. 7 illustre un gain découlant de la mise en œuvre d’un procédé de défloutage 100 conforme à l’invention. Afin de simplifier la description, une même référence est utilisée dans différentes figures pour désigner un même objet, élément ou étape. Ainsi, lorsque la description cite un objet, élément ou étape référencé, cet objet, élément ou étape pourra être identifié sur plusieurs figures. En outre, les figures ainsi que la description sont données à titre d’exemples non limitatifs de réalisation. En préambule, au sens de l’invention, une image est une représentation numérique matricielle BGR d’éléments nommés « voxels », d’épaisseur k (selon un axe transversal), BGR(i,j), i et j étant des indices de valeurs entières pour identifier le voxel situé à la ligne i et à la colonne j de la matrice BGR, soit une représentation graphique sous la forme d’un tableau ou matrice, d’éléments BGR(i,j) encodant chacun un triplet de valeurs entières comprises entre zéro et deux cents cinquante-cinq, selon le codage de couleurs RVB, acronyme pour « Rouge Vert Bleu », également connu sous le sigle RGB, acronyme anglo-saxon de « Red Green Blue ». Un tel codage informatique des couleurs est le plus proche des matériels actuellement disponibles pour constituer des interfaces homme-machine de sortie telles que des écrans d’ordinateur. En général, ces derniers reconstituent une couleur par synthèse additive à partir de trois couleurs primaires, rouge, vert et bleu, formant sur l'écran une mosaïque généralement trop petite pour être discriminée par un œil humain. Le codage RVB indique une valeur d’intensité lumineuse pour chacune de ces couleurs primaires. Une telle valeur est généralement codée sur un octet et appartient donc à un intervalle de valeurs entières comprises entre zéro et deux cents cinquante-cinq. D’autres codages de couleurs pourraient en variante être exploités. Dans ce cas, chaque voxel ou élément du tableau BGR décrirait un ensemble de valeurs numériques adapté audit codage en lieu et place du triplet évoqué ci-dessus pour le codage RVB. Un procédé 100 de défloutage d’une image médicale, conforme à l’invention et illustré en FIG. 4 , se traduit avantageusement sous la forme d’un produit programme d’ordinateur dont les instructions de programme sont destinées à être implantées dans la mémoire de programmes d’un élément d’un système d’imagerie médicale, tel que le système S selon les figures 1 et 2, par exemple un ordinateur ou un serveur informatique ou, plus généralement, de tout objet électronique disposant d’une puissance de calcul suffisante. La FIG. 4 illustre ainsi un procédé 100, conforme à l’invention, de défloutage d’une image expérimentale Mn fruit d’une séquence d’acquisition par un appareil d’imagerie, tel que l’appareil 1 du système d’analyse d’imagerie médicale S illustré par les figures 1 et 2, en lien avec une pluralité de n voxels. Une telle image expérimentale Mn peut être considérée comme tout ou partie de l’image « totale » résultant de ladite séquence d’acquisition. L’image Mn peut donc être l’image totale ou une portion de l’image totale représentant une zone d’intérêt. Un tel procédé 100 peut être mis en œuvre par l’unité de traitement 4 d’un tel système d’analyse d’images médicales S. À des fins illustratives mais non limitatives, un tel procédé 100 selon l’invention est avantageusement décrit dans la suite de la présente et mis en œuvre pour le traitement d’une image médicale, représentant le cerveau d’un patient, fruit d’une séquence d’acquisition de perfusion par un appareil d’imagerie médicale par résonance magnétique 1 durant laquelle un traceur endogène passe au sein de ladite pluralité de n voxels. Un tel procédé de défloutage 100 selon l’invention consiste à analyser une image médicale expérimentale Mn dans le domaine fréquentiel et de comparer ladite image Mn à une image de référence M0 de même résolution, à savoir le même nombre de voxels par unité de longueur, que ladite image Mn afin de récupérer les fréquences perdues ou tronquées de l’image expérimentale Mn. Ainsi, ladite image de référence M0 représente une zone d’intérêt similaire à ladite image Mn, pour exemple un cerveau tel qu’illustré sur la FIG. 4 , selon le même plan de coupe (les plans de coupe ont été explicités précédemment et illustrés en FIG. 3 ) pour exemple en coupe coronale telle qu’illustrée en FIG. 4 . Dans un mode de réalisation préféré afin que le procédé de défloutage 100 soit le plus efficient, il est possible que ladite image de référence M0 puisse être une image structurelle dudit patient. Une telle image structurelle M0 représente la même zone d’intérêt et est de même résolution que ladite image expérimentale Mn. Dans ce cas de figure illustré en FIG. 5 , ledit procédé 100 peut comporter au préalable une étape 101. Une telle étape 101 consiste en une étape d’acquisition d’une image structurelle dudit patient fruit d’une séquence d’acquisition par un appareil d’imagerie médicale. L’acquisition de l’image structurelle de référence M0 peut être réalisée en même temps que l’acquisition de l’image expérimentale Mn par un même appareil d’imagerie médicale, tel que par exemple l’appareil d’imagerie médicale par résonance magnétique 1. À titre illustratif mais non limitatif, l’image structurelle M0 du patient peut être : une image dite « pondérée en T1 », souvent utilisée pour l’anatomie, dans la mesure où une telle image résulte d’une séquence d’acquisition, dite « anatomique », favorisant la détection de l’eau peu mobile c’est-à-dire intracellulaire. Pour une telle image, la graisse apparaît hyper-intense, de couleur claire et l’eau hypo-intense, de couleur sombre ; une image dite « pondérée en T2 », souvent utilisées comme image fonctionnelle, dans la mesure où une telle image résulte d’une séquence d’acquisition favorisant la détection de l’eau mobile c’est-à-dire extracellulaire ou intravasculaire. Pour une telle image, l’eau apparaît hyper-intense, de couleur claire et la graisse apparaît un peu plus sombre que l’eau ; une image dite « FLAIR » pour « Fluid Attenuated Inversion Recovery » selon une terminologie anglo-saxonne, très adaptée à l’imagerie cérébrale, dans la mesure où une telle image résulte d’une séquence d’acquisition supprimant le signal provenant du liquide céphalo-rachidien. Pour une telle image, les lésions de la substance blanche, les ramollissements (Accident Vasculaire Cérébrale), démyélinisation inflammatoire (sclérose en plaques) apparaissent hyper-intenses et sont particulièrement bien mises en exergue. L’invention ne saurait être limitée à l’utilisation de ce type d’images structurelles. Tout autre image structurelle pourra être utilisée pour la mise en œuvre du procédé 100 selon l’invention. En variante ou en complément, une telle image de référence M0 peut émaner : d’un atlas, à savoir une image type résultant d’une moyenne d’images acquises d’après diverses populations de patients ou d’un même patient ; d’un sujet sain autrement nommé patient « témoin » ; ou encore d’un référentiel anatomique artificiel. L’utilisation de telles images peut trouver tout son intérêt par exemple lorsque l’image structurelle du patient est perdue ou malencontreusement effacée, ou lorsque des problèmes d’acquisition de l’image structurelle du patient sont rencontrés, dus par exemple à des mouvements du patient lors de l’acquisition, ces mouvements pouvant être volontaires (par exemple dans le cas de claustrophobie) ou involontaires (par exemple dans le cas de certaines maladies telles que la maladie de Parkinson). En complément, afin de garantir un alignement parfait entre l’image de référence M0 et l’image expérimentale Mn, de manière optionnelle, tel qu’illustré en FIG. 5 , ledit procédé 100 peut comporter une étape préalable 102 de recalage de l’image expérimentale Mn par rapport à l’image de référence M0 ou inversement de l’image de référence M0 par rapport à l’image expérimentale Mn afin d’assurer la correspondance des deux images M0 et Mn voxel par voxel. En effet, plus l’image de référence M0 est proche de l’image expérimentale Mn, meilleur le résultat de défloutage sera. A l’issue d’un tel recalage optionnel, les deux images étant sensiblement de mêmes dimensions sont superposables car présentant des orientations communes. Comme évoqué précédemment, l’invention consiste à analyser une image médicale expérimentale Mn dans le domaine fréquentiel. Ainsi, tel qu’illustré en figures 4 et 5, ledit procédé 100 comporte des première 110 et deuxième 120 étapes de calcul des transformées de Fourier TFMn et TFM0 respectivement de l’image expérimentale Mn et l’image de référence M0. À ce titre, la transformée de Fourier est une technique mathématique permettant de déterminer le spectre de fréquences d'un signal résultant d’une séquence d’acquisition. Ledit signal résultant d’une séquence d’acquisition étant un signal discret, la transformée de Fourier discrète pourra être privilégiée. En effet, l’analyse de signaux discrets par transformée de Fourier discrète permet de représenter un signal discret dans le domaine fréquentiel : il est alors possible de traiter directement les fréquences de l’image. À ce titre, pour obtenir les fréquences contenues dans ledit signal, une opération nommée transformation est appliquée et le résultat de cette transformation est appelé transformée. En pratique, pour calculer une transformée de Fourier discrète, il est avantageusement utilisé un algorithme appelé « transformée de Fourier rapide », connu sous l’acronyme anglo-saxon FFT pour « Fast Fourier Transform » , permettant un calcul plus rapide. Toutefois, l’invention ne saurait être limitée à l’utilisation de cet algorithme. Tout autre algorithme permettant de calculer une transformée de Fourier pourra être utilisé pour la mise en œuvre du procédé 100 selon l’invention. Les images obtenues en appliquant une transformée de Fourier donnent une image complexe. Ainsi, en général, le module de la transformée de Fourier est calculé et représenté. Un tel module est calculé comme étant la racine carrée de la somme du carré de la partie imaginaire et du carré de la partie réelle. À ce titre, tel qu’illustré en FIG. 4 , ledit procédé 100 conforme à l’invention comporte une étape 130 de détermination du module YMn de la transformée de Fourier TFMn de l’image expérimentale Mn et du module YM0 de la transformée de Fourier TFM0 de l’image expérimentale M0. Avantageusement, une telle étape 130 peut comporter un traitement de normalisation desdits modules YMn et YM0 afin que lesdits modules partagent une échelle de valeurs commune. Un telle normalisation est réalisée pour chaque module YMn et YM0 et consiste par exemple à diviser les valeurs de chaque module YMn et YM0 respectivement par leur valeur de fréquence centrale, autrement dite « fréquence zéro ». En variante, ladite normalisation peut consister à utiliser tout système de normalisation permettant aux deux modules YMn et YM0 de partager la même échelle de valeurs. En complément, chaque pixel dudit module agissant en tant que diviseur est préférentiellement non-nul pour prévenir toute division par zéro. Pour cela, il peut être envisagé d’ajouter une constante, telle que, par exemple, un dixième ou un centième, afin de ne pas déséquilibrer l’échelle. Une fois que lesdits modules YMn et YM0 sont déterminés, ledit procédé 100 comporte une étape 140 d’obtention d’un masque de défloutage D, tel qu’illustré en FIG. 4 . Une telle étape 140 consiste à calculer le ratio D entre les modules YMn et YM0. Préférentiellement mais non limitativement, le masque de défloutage D ou ratio D est égal à la division du module YM0 de la transformée de Fourier TFM0 de l’image de référence M0 par le module YMn de la transformée de Fourier TFMn de l’image expérimentale Mn. Avantageusement, ladite étape 140 peut comporter un traitement subsidiaire de limitation des valeurs dudit ratio ou dudit masque de défloutage D consistant à écrêter les valeurs de D à une amplitude maximale. Préférentiellement, une telle amplitude maximale est comprise entre 2 et 100, idéalement 5 à 20 pour éviter tout risque d’amplification de bruit. Il est à noter qu’une amplitude de valeur 10 semble être la valeur idéale. Ensuite, ledit procédé 100 comporte une étape 150 d’application dudit masque de défloutage ou dudit ratio D à la transformée de Fourier TFMn de l’image expérimentale Mn par multiplication. À l’issue d’une telle étape 150, on obtient une transformation de Fourier corrigée TFMnc. Ainsi, ledit masque de défloutage D est appliqué à l’entièreté de la transformée de Fourier TFMn soit à l’ensemble des fréquences et donc à l’ensemble de l’image expérimentale Mn. Ladite étape 150 pourrait aussi ne concerner que certaines fréquences. Enfin, ledit procédé 100 comporte une dernière étape 160 d’obtention de l’image expérimentale défloutée Mnc. Une telle étape 160 consiste à appliquer la transformée de Fourier inverse sur ladite transformée de Fourier corrigée TFMnc. En effet, la transformée de Fourier inverse est une opération mathématique permettant de reconstruire l’image d’origine (le signal d’origine) à partir de sa représentation fréquentielle. Dans le cas de l’invention, on obtient donc une image expérimentale défloutée Mnc. En complément, tel qu’illustré en FIG. 5 , ledit procédé 100 peut comporter une étape optionnelle 142, préalable à l’étape 150 d’application dudit ratio D et subséquente à l’étape 140 d’obtention d’un masque de défloutage. Une telle étape 142 consiste en l’application d’un filtre de lissage sur ledit masque de défloutage D (autrement nommé ledit ratio D). À l’issue d’une telle étape 142, on obtient donc un masque de défloutage lissé D L , permettant de diminuer les artefacts. Une telle opération 142 de lissage peut avantageusement consister en l’application d’un filtre gaussien. Un tel filtre gaussien présente, préférentiellement, un écart-type compris entre 1 et 50, au-delà d’un écart-type de 50 le lissage perdrait tout son intérêt, et idéalement un écart-type de valeur 10. L’invention ne saurait être limitée à l’utilisation de ce type de filtre de lissage ni même aux valeurs énoncées comme idéales. Tout autre filtre de lissage pourrait être utilisé pour la mise en œuvre du procédé 100 selon l’invention. En outre, une telle étape 142 peut avantageusement consister à normaliser ledit ratio D. Dans ce cas, une telle normalisation peut consister par exemple à diviser les valeurs dudit ratio D par la valeur de sa fréquence centrale, dite fréquence zéro ou consister en toute autre technique de normalisation. En variante, si plusieurs images expérimentales Mn1, Mn2, …, Mnx ont été acquises et donc disponibles pour traitement, les étapes 110, 130 et 140 (ainsi que les étapes optionnelles 101, 102 et/ou 142 si nécessaire) sont réalisées indépendamment pour chaque image expérimentale afin de calculer pour chacune un masque de défloutage individuel D1, D2, …, Dx. Le masque de défloutage final Df est alors égal à la moyenne, avantageusement arithmétique, voire la médiane, de tous les masques de défloutage individuels calculés D1, D2, …, Dx. Les figures 6 et 7 permettent d’illustrer le bénéfice découlant de la mise en œuvre d’un procédé de défloutage 100 conforme à l’invention. En effet, pour mesurer et illustrer l’impact d’un tel procédé de défloutage 100, la FIG. 6 montre la différence de contraste entre une image expérimentale Mn et l’image expérimentale défloutée Mnc. Il est également possible, tel qu’illustré en FIG. 7 , de calculer les modules respectifs YMn et YMnc des transformées de Fourier de l’image expérimentale Mn et de l’image expérimentale défloutée Mnc. De tels modules YMn et YMnc sont représentés sur la FIG. 7 et les cercles en pointillés blancs mettent bien en évidence que des fréquences perdues dans le cas de l’image expérimentale Mn ont été restaurées dans le cas de l’image expérimentale défloutée Mnc, l’image représentant le module YMnc apparaîssant plus contrastée que celle représentant le module YMn. Il sera apprécié de l’Homme du métier que la présente divulgation n’est pas limitée à ce qui est particulièrement montré et décrit ci-dessus. D’autres modifications peuvent être envisagées sans sortir du cadre de la présente invention définie par les revendications ci-annexées. Notamment, dans l’exemple préféré décrit précédemment, le procédé 100 est mis en œuvre pour déflouter une image, représentant le cerveau d’un patient, fruit d’une séquence d’acquisition de perfusion ASL par un appareil d’imagerie médicale par résonance magnétique. Or, un tel procédé 100 conforme à l’invention pourrait être utilisé dans le cadre de toutes techniques d’imagerie médicale, dès lors qu’une image présente un flou, par exemple pour la technique PET acronyme anglo-saxon pour tomographie par émission de positrons, pour laquelle la qualité de l’image acquise peut être limitée par exemple par la non-colinéarité des paires de photons émises, la diffusion intercristalline, ou encore la pénétration des cristaux. Procédé (100) de défloutage d’une image médicale expérimentale (Mn) fruit d’une séquence d’acquisition par un appareil d’imagerie médicale (1) décrivant une pluralité de volumes élémentaires d’intérêt d’un organe d’un patient, nommés « voxels », ledit procédé (100) étant mis en œuvre par une unité de traitement (4) d’un système d’analyse d’imagerie médicale (S), ledit procédé (100) comportant : une étape (110) de calcul de la transformée de Fourier (TFMn) de ladite image expérimentale (Mn) ; une étape de calcul (120) de la transformée de Fourier (TFM0) d’une image de référence (M0) de même résolution que l’image expérimentale (Mn) ; une étape de détermination (130) d’un premier module (YMn) de la transformée de Fourier (TFMn) de l’image expérimentale (Mn) et d’un deuxième module (YM0) de la transformée de Fourier (TFM0) de l’image de référence (M0) ; une étape d’obtention (140) d’un masque de défloutage (D) consistant à calculer le ratio (D) entre les premier et deuxième modules (YMn, YM0) ; une étape d’application (150) dudit ratio (D) à la transformée de Fourier (TFMn) de l’image expérimentale (Mn) par multiplication pour obtenir une transformée de Fourier corrigée (TFMnc) ; une étape d’obtention (160) d’une image expérimentale défloutée (Mnc) consistant à appliquer la transformée de Fourier inverse sur ladite transformée de Fourier corrigée (TFMnc). Procédé (100) selon la revendication précédente, pour lequel l’image de référence (M0) émane d’un atlas, d’un patient témoin ou d’un référentiel anatomique artificiel. Procédé (100) selon la revendication 1 comportant une étape (101) préalable aux étapes (110, 120) de calcul de transformée de Fourier consistant en une étape (101) d’acquisition d’une image structurelle dudit patient de même résolution que l’image expérimentale (Mn), ladite image structurelle correspondant à ladite image de référence (M0). Procédé (100) selon l’une des revendications précédentes comportant une étape (102) préalable aux étapes (110, 120) de calcul de transformée de Fourier consistant en une étape (102) de recalage de l’image expérimentale (Mn) par rapport à l’image de référence (M0) ou de l’image de référence (M0) par rapport à l’image expérimentale (Mn) de sorte à assurer leur correspondance voxel par voxel. Procédé (100) selon l’une des revendications précédentes comportant une étape (142) préalable à l’étape (150) d’application dudit ratio (D) et subséquente à l’étape (140) d’obtention d’un masque de défloutage (D) consistant à appliquer un filtre de lissage sur ledit masque de défloutage (D) pour obtenir un masque de défloutage lissé (D L ). Procédé (100) selon la revendication précédente, pour lequel le filtre de lissage est un filtre gaussien. Procédé (100) selon l’une des revendications précédentes, pour lequel les transformées de Fourier (TFMn, TFM0) de l’image expérimentale (Mn) et de l’image de référence (M0) sont respectivement calculées avec l’algorithme de transformée de Fourier rapide. Procédé (100) selon l’une des revendications précédentes pour lequel la séquence d’acquisition de l’image expérimentale (Mn) est une séquence d’acquisition de perfusion par un appareil d’imagerie médicale par résonance magnétique (1) durant laquelle un traceur endogène passe au sein de ladite pluralité de volume élémentaire d’intérêt. Produit programme d’ordinateur comportant une ou plusieurs instructions de programme exécutables par l’unité de traitement (4) d’un ordinateur, lesdites instructions de programme étant chargeables dans une mémoire non volatile dudit ordinateur et dont l’exécution par ladite unité de traitement (4) provoque la mise en œuvre d’un procédé (100) selon l’une quelconque des revendications précédentes. Support de mémorisation lisible par un ordinateur comportant les instructions d’un produit programme d’ordinateur selon la revendication précédente. Système d’analyse d’imagerie médicale (S) comportant une unité de traitement (4) agencée pour déflouter une image médicale expérimentale (Mn) fruit d’une séquence d’acquisition par un appareil d’imagerie médicale (1) en lien avec une pluralité de volumes élémentaires d’intérêt, nommés voxels, ladite unité de traitement (4) étant configurée pour : calculer la transformée de Fourier (TFMn) de ladite image expérimentale (Mn) ; calculer la transformée de Fourier (TFM0) d’une image de référence (M0) de même résolution que l’image expérimentale (Mn) ; déterminer un premier module (YMn) de la transformée de Fourier (TFMn) de l’image expérimentale (Mn) et un deuxième module (YM0) de la transformée de Fourier (TFM0) de l’image de référence (M0) ; obtenir un masque de défloutage (D) en calculant le ratio (D) entre les premier et deuxième modules (YMn, YM0) ; appliquer ledit ratio (D) à la transformée de Fourier (TFMn) de l’image expérimentale (Mn) par multiplication pour obtenir une transformée de Fourier corrigée (TFMnc) ; obtenir une image expérimentale défloutée (Mnc) en appliquant la transformée de Fourier inverse sur ladite transformée de Fourier corrigée (TFMnc). Système d’analyse d’imagerie médicale (S) selon la revendication précédente comportant une mémoire de programmes comportant les instructions de programme d’un produit programme d’ordinateur conforme à la revendication 9.