Brevet FR3099621A1 — publié le 20210205

Description Titre de l'invention : PROCEDE DE GENERATION D'UN BIO- MARQUEUR, SYSTEME. [0001] Le domaine de l'invention concerne le domaine des procédés pour générer des bio- marqueurs de la région pulmonaire. Plus particulièrement, l'invention se rapporte aux méthodes de traitement d'une image afin de quantifier la présence de caractéristiques susceptibles d'être associées à une pathologie donnée. Le domaine de l'invention s'applique plus particulièrement au traitement des images acquises par IRM. [0002] Actuellement, la quantification des phénomènes inflammatoires et de remodelage des voies aériennes est principalement réalisée à partir d'un système d'imagerie par scanner. Cette technique plus connue sous le nom de tomodensitométrie, TDM, ou CTscan permet d'obtenir des images tridimensionnelles offrant une résolution et des niveaux de contrastes suffisants pour élaborer des stratégies de quantification. D'une manière plus générale, l'imagerie CT-scan permet d'obtenir des indicateurs de l'extension 3D des anomalies des voies pulmonaires. [0003] Toutefois, cette technique présente l'inconvénient d'être irradiante. D'autre part, cette technique ne permet pas de discriminer l'inflammation du remodelage des voies aériennes. Il existe un besoin d'élaborer des biomarqueurs utilisant une imagerie non irradiante, et permettant de caractériser le signal des phénomènes inflammatoires actifs versus le signal des lésions cicatricielles définitives, telle que le permet l'imagerie par IRM. [0004] Or, à ce jour, il est difficile de définir des modalités permettant de traiter les images des régions pulmonaires par IRM du fait que les tissus et les structures cellulaires de ces régions sont variés et difficilement segmentables sans perdre de données. D'autre part, le signal en IRM n'est pas calibré comme pour un scanner, il n'est donc pas possible d'utiliser les données brutes des valeurs de signal sur une image donnée, tandis que les cartographies paramétriques des valeurs réelles du Tl et du T2, sur un poumon, sont de natures à être « artéfactés », c'est-à-dire ayant des erreurs d'acquisition, par la susceptibilité magnétique du poumon normal. Il n'existe pas à cc jour de modalité d'acquisition qui permet de définir une segmentation adéquate des voies aériennes offrant des possibilités de générer des biomarqueurs suffisamment efficaces pour suppléer une acquisition CT-scan. [0005] Il existe donc un besoin de définir une solution qui puisse répondre à cette pro- blématique. La présente invention vise à résoudre les inconvénients précités. [0006] Selon un premier aspect, l'invention concerne un procédé de génération d'un bio- marqueur comportant : 2 [0007] - Acquisition d'une image par un système IRM traitement de ladite image IRM pour générer une image tridimensionnelle du poumon ; génération d'une première fonction correspondant à la distribution des différentes valeurs de l'intensité du signal de chaque voxel d'une portion de l'image tridimensionnelle acquise ; calcul automatique d'au moins un seuil de filtrage de ladite première fonction à partir d'au moins une seconde fonction de distribution de différentes valeurs d'intensité du signal de chaque voxel d'une portion de l'image tridimensionnelle acquise ; segmentation d'un volume comportant : - Un volume principal, correspondant au volume pulmonaire ; - Un volume filtré d'un volume de voxels quantifiés par la première fonction et filtré au moyen d'au moins du seuil de filtrage calculé, normalisation des valeurs de l'image tridimensionnelle du volume pulmonaire à partir des valeurs absolues des valeurs d'intensité du signal des voxels de l'image et d'au moins le seuil de filtrage calculé; génération d'un biomarqueur indiquant un ratio de volume segmenté normalisé. [0008] Le volume pulmonaire correspond au volume délimité par les parois du ou des poumon(s). [0009] Avantageusement, l'invention permet une automatisation de la génération d'un bio- marqueur tout en définissant des seuils adaptifs selon les patients. Un autre avantage est de permettre de configurer une acquisition pour générer différents biomarqueurs selon le cas d'usage que l'on souhaite obtenir. [0010] L'étape de génération d'une première fonction correspondant à la distribution des différentes valeurs de l'intensité du signal de chaque voxel d'une portion de l'image tridimensionnelle acquise correspond, autrement dit, à l'application d'une première fonction entrainant la génération d'une distribution des différentes valeurs de l'intensité du signal de chaque voxel d'une portion de l'image tridimensionnelle acquise. [0011] Selon un exemple, le procédé est mis en oeuvre par une fonction au sein d'une console d'exploitation des images IRM acquises ou au sein d'un équipement comportant une interface physique ou sans fil avec la console d'exploitation du système IRM. Afin de mettre en oeuvre le procédé de l'invention, une interface pour recevoir les images IRM et une mémoire pour les stocker peuvent être utilisées. [0012] L'équipement permettant de mettre en oeuvre le procédé de l'invention comprend avantageusement un afficheur et une interface graphique afin de configurer un pa- 3 ramétrage selon le biomarqueur que l'on souhaite générer. [0013] Selon un mode de réalisation, le procédé comprend [0014] - l'acquisition d'une image tridimensionnelle par un système IRM configure° par : - une pondération en T, ; - un temps d'écho TE supérieur à seuil prédéfini ; le calcul automatique d'au moins un seuil comprenant : - acquisition d'un volume de référence ; - génération d'une fonction de référence correspondant à une distribution des différentes valeurs de l'intensité du signal de chaque voxel du volume de référence, ladite fonction de référence étant la seconde fonction de distribution ; - calcul de l'écart-type de la fonction de référence ; - détermination d'un seuil de référence à partir de l'écart-type calculé de la fonction de référence. [0015] Un avantage est de permettre de générer un biomarqueur permettant de quantifier un volume de charge en mucus et/ou d'oedème sans mettre en oeuvre une technique irradiante. [0016] Selon un exemple de ce mode de réalisation, l'acquisition est pararnetrée et comprend : [0017] - Une acquisition d'une image composée d'une pluralité d'images acquises sur un nombre de cycle d'au moins 4 temps d'écho, les temps d'échos étant configurés selon des durées croissantes ; Une séquence Echo de spin ; Un paramétrage visant à émettre un signal pour pré-saturer ou une saturer du signal acquis. [0018] Un avantage est notamment d'améliorer le traitement des volumes circulants afin de les supprimer du traitement d'images. [0019] Selon un exemple de ce même mode de réalisation, une seconde acquisition de l'image est réalisée avec une première configuration pour détourer le volume pulmonaire, une étape de traitement d'image étant réalisée pour combiner l'image acquise en pondération T2 avec l'image acquise par la seconde acquisition. La première configuration définit un paramétrage d'une acquisition en pondération TI ou en densité de protons. Dit autrement, la première acquisition permet de configurer une pondération de type TI ou une pondération en densité de proton. [0020] Un avantage est de construire un masque permettant de détourer le volume pulmonaire et des organes comme le coeur afin de traiter l'image acquise en pondération T,. En effet, généralement l'image acquise en pondération T, ne permet pas d'obtenir une segmentation du volume pulmonaire qui permet à elle seule de détourer les contours des organes. [0021] Selon un exemple de ce mode de réalisation, une opération de fusion entre au moins une image acquise en pondération T1 avec un temps d'écho ultra court UTE et au moins une image acquise par une pondération T2 est réalisé pour générer une image dont les données provenant de chacune des images acquises ont été combinées pour segmenter le volume pulmonaire. [0022] Un avantage est de générer une image à partir de la fusion de deux images, l'image produite offrant une capacité d'isolation optimale des zones d'intérêt au sein d'un volume pulmonaire. [0023] Selon un exemple de ce mode de réalisation, le seuil de référence est déterminé à partir d'une combinaison entre une valeur de répartition de référence de la fonction de référence Ft, et une valeur comprise entre 10 et 20 fois la valeur de l'écart-type 6A de la fonction de référence FA. [0024] Un avantage est de permettre d'éliminer un niveau de bruit. Un autre avantage est de définir un seuil adaptatif afin de générer un seuil compatible avec un grand nombre de patients. [0025] Selon un exemple de ce mode de réalisation, la valeur de répartition de référence de la fonction de référence FA est le mode principal MF de la répartition des valeurs d'intensité du signal de l'image acquise en pondération T2 au sein du volume pulmonaire V,. Un avantage est d'utiliser une valeur de référence pour réaliser un filtre qui puisse être appliqué à différents types de patients ayant des profils différents. Ainsi, le filtre est adaptatif et permet de s'affranchir de la configuration d'acquisition ou des profils de patients. La sélection du mode principal permet de définir un seuil adaptif pour chaque patient. [0026] Selon un exemple de ce mode de réalisation, le procédé comprend une étape de nor- malisation, la normalisation comprenant le calcul d'un indice volumique du signal VIP à partir des valeurs absolues du signal du volume filtré VF, du volume résultant du volume filtré VF et du volume pulmonaire V. Un avantage est de quantifier des niveaux d'inflammation, par exemple, d'une zone du poumon. Un autre avantage est qu'il est possible de quantifier la propagation d'une inflammation dans le temps. [0027] Selon un autre mode de réalisation, le procédé de génération d'un biomarqueur est caractérisé en ce que : [0028] - l'acquisition d'une image tridimensionnelle par un système IRM, est configurée par : - une pondération en densité de proton ou en TI.le temps d'écho ; - un temps d'écho inférieur à seuil prédéfini ; le calcul automatique d'au moins un seuil comprend : - modélisation d'au moins deux fonctions gaussiennes par ajustement de la première fonction ; - détermination du premier seuil de filtrage par un calcul de l'intersection de la première fonction gaussienne et de la seconde fonction gaussienne ; - détermination d'un second seuil correspondant à une valeur minimale de la première fonction gaussienne et une valeur minimale d'une valeur de l'intensité du signal d'un voxel ; [0029] - Le volume filtré est un premier volume correspondant aux voxels quantifiés par la première fonction gaussienne comprise entre le premier seuil et le second seuil, lesdits voxels correspondant à un milieu aérien, [0030] - normalisation des valeurs de l'image tridimensionnelle du volume pulmonaire à partir du premier seuil et du second seuil calculés, génération d'un premier biomarqueur indiquant un ratio d'un volume caractéristique du volume segmenté normalisé, ledit ratio étant calculé entre le volume caractéristique et le volume pulmonaire. [0031] Un avantage de cette configuration est d'obtenir une bonne segmentation du poumon tout en isolant des volumes caractéristiques grâce à un biomarqueur qui peut être représentatif d'un phénomène lié à l'emphysème. [0032] Selon un exemple de cc mode de réalisation, le temps d'écho est inférieur à lms. Un avantage d'une configuration d'acquisition avec un TE relativement court est d'être insensible aux courants de Foucault. Cette configuration permet de maximaliser le signal dans le parenchyme pulmonaire [0033] Selon un exemple de ce mode de réalisation, la segmentation comprend la dé- finition d'un second volume correspondant aux voxels quantifiés par la seconde fonction gaussienne supérieure au premier seuil, lesdits voxels correspondant à un milieu intermédiaire ou graisseux. Un avantage est d'améliorer la segmentation des différents volumes dans le poumon. [0034] Selon un exemple de ce mode de réalisation, la segmentation comprend une étape d'extraction d'un volume caractéristique comportant des voxels du volume filtré dont la valeur d'intensité du signal est inférieure à un troisième seuil prédéfini, ledit seuil prédéfini étant déterminé sur une échelle normalisée de [0 ;1]. Un avantage est d'obtenir un biomarqueur représentatif du phénomène d'emphysème. [0035] Selon un exemple de ce mode de réalisation, la segmentation comprend une étape d'exclusion/suppression des voxels déconnectés à leur voisinage de voxels d'une même quantification. Un avantage est de de réduire les artefacts de l'acquisition et de sélectionner les volumes correspondant réellement à un volume caractéristique d'une pathologie. 6 [0036] Selon un exemple de ce mode de réalisation, la modélisation des fonctions gaussiennes comprend : [0037] - Un lissage gaussien appliqué à l'image acquise avec réduction du temps de codage ponctuel et avec acquisition radiale pour le débruitage ; Détourage des contours par application d'un filtre local ; Utilisation de la méthode d'ajustement de courbe représentant la fréquence de distribution de voxels. [0038] Un avantage est d'obtenir une modélisation réaliste et fidèle d'une répartition de l'intensité des valeurs du signal dans l'image. La fonctionnalisation du traitement à partir des fonctions gaussiennes permet de générer un traitement identique pour chaque patient qui soit représentatif et en même temps adaptatif. [0039] Selon un mode de réalisation, le procédé comprend l'acquisition des valeurs d'intensité du signal de chaque voxel quantifiées par la première fonction, les valeurs d'intensité correspondant à une donnée de contraste d'image. [0040] Selon un mode de réalisation, l'acquisition est réalisée de manière synchronisée avec un respirateur. Selon un exemple, le respirateur est un navigateur ou une ceinture respiratoire. Un avantage est de s'affranchir des mouvements du poumon lors de l'acquisition de l'image. [0041] Selon un mode de réalisation, une étape d'extraction d'une image volumétrique est réalisée à partir des images acquises par l'IRM, ladite image extraite étant réalisée à un instant déterminé de la séquence. [0042] Selon un autre mode de réalisation, l'acquisition de l'image tridimensionnelle est réalisée au moyen d'un empilement d'images 2D acquises, l'épaisseur de la coupe étant égale au moins à la largeur d'un voxel. [0043] Selon un autre aspect, l'invention concerne un système comportant au moins un cal- culateur et une mémoire et une interface pour recevoir une image acquise par un système IRM, ledit système étant configuré pour mettre en oeuvre les étapes du procédé de l'invention. Avantageusement, un afficheur permet de visualiser les images IRM acquises et les images traitées. Les biomarqueurs peuvent être représentés graphiquement au sein d'une image générée. Selon un exemple, une interface graphique permet de régler les paramètres d'acquisition et de définir les zones d'intérêt ou une pathologie recherchée. [0044] D'autres caractéristiques et avantages de l'invention ressortiront à la lecture de la description détaillée qui suit, en référence aux figures annexées, qui illustrent : [0045] Fig. 1 : des étapes d'un mode de réalisation d'un procédé de l'invention ; [0046] Fig. 2 : un premier mode de réalisation du procédé de l'invention appliqué pour quantifier un volume pulmonaire aérien afin de générer un premier biomarqueur ; [0047] Fig. 3 : un second mode de réalisation du procédé de l'invention appliqué pour 7 quantifier un volume d'encombrement des voies aériennes afin de générer un second biomarqueur [0048] Fig. 4 : une coupe 2D d'une image tridimensionnelle du poumon, la coupe ou l'image tridimensionnelle étant acquise au moyen d'un système IRM ; [0049] Fig. 5 : un exemple d'utilisation de deux fonctions pour déterminer un premier seuil afin de générer un premier hiomarqueur normalisé selon le premier mode de réalisation [0050] Fig. 6A : un exemple d'une courbe quantifiant les différents contrastes de l'image 2D obtenue par un scanner ainsi qu'un seuil permettant de générer un indicateur représentant une quantification de l'emphysème d'un patient atteint d'une bronchopneumopathie chronique obstructive ; [0051] Fig. 6B : un exemple d'une courbe quantifiant les différents contrastes de l'image 2D obtenue par un scanner ainsi qu'un seuil permettant de générer un indicateur représentant une quantification de l'emphysème d'un patient sain ; [0052] Fig. 7A : un exemple d'une courbe quantifiant les différents contrastes de l'image 2D obtenue par une IRM et traitée selon le procédé de l'invention pour calculer automatiquement un seuil permettant de générer un indicateur représentant une quantification de l'emphysème d'un patient atteint d'une bronchopneumopathic chronique obstructive ; [0053] Fig. 7B : un exemple d'une courbe quantifiant les différents contrastes de l'image 2D obtenue par une IRM et traitée selon le procédé de l'invention pour calculer automatiquement un seuil permettant de générer un indicateur représentant une quantification de l'emphysème d'un patient sain ; [0054] Fig. 8 : une représentation d'une quantité d'un volume encombré (ex : sécrétion de mucus) d'un volume pulmonaire chez un patient dans un plan de coupe d'une IRM, [0055] Fig. 9 : une représentation d'un seuil de référence permettant de filtrer les voxels de hautes intensités afin de générer un masque. Définitions [0056] Dans la suite de la description, une image tridimensionnelle comporte une pluralité de pixels volumiques, appelés également voxel. On nommera indifféremment un pixel ou un voxel de l'image. [0057] Dans la suite de la description, on nommera un « hyposignal » un signal de faible intensité, apparaissant sur l'image IRM comme une plage noire ou gris foncé. C'est le cas en l'absence de protons (air), des substances à T, court (tissus calcifiés (os cortical, émail, calcifications), des tissus riches en collagène (tendons, ligaments, fascias...), des liquides riches en protéines, des liquides stagnants (urine, LCS, épanchements, oedème, inflammation...), des substances à effet paramagnétique (gadolinium concentré), du fer en concentration élevée (hémochromatose), des certains hématomes, 8 etc. [0058] Dans la suite de la description, on nommera un « hypersignal », un signal de forte intensité, apparaissant sur l'image IRM comme une plage blanche ou gris pâle. En TI, c'est le signal des substances à TI court : lipides (graisse, tumeurs graisseuses, moelle osseuse grasse...), liquides riches en protéines (mucocèles, kystes tumoraux...), substances à effet paramagnétiques (gadolinium...), fer (hématome subaigu), radicaux libres (mélanine), posthypophyse.... En T2, c'est le cas de l'eau libre (urine, LCS, épanchements, liquide synovial...) de l'eau interstitielle (oedème, inflammation...), du sang stagnant, des liquides pauvres en protéines, etc. [0059] La figure 1 représente les principales étapes d'un exemple de réalisation du procédé de l'invention. Les étapes comprennent une acquisition d'une image ou d'une séquence d'images par un système MM. L'acquisition est, par exemple, réalisée par acquisition de coupes 2D qui sont ensuite empilées pour former une image tridimensionnelle. Selon un autre exemple, l'acquisition est une acquisition directement tridimensionnelle qui permet d'obtenir une image 3D. L'invention est compatible également d'une acquisition dite « 4D » dans laquelle l'échelle temporelle est prise en compte pour générer et reconstruire une image par un algorithme de traitement d'images. [0060] Selon un mode de réalisation, un dispositif de mesure du rythme respiratoire est utilisé pour acquérir le signal avec une synchronisation respiratoire. Un respirateur peut être utilisé ou tout autre dispositif permettant de synchroniser le rythme respiratoire avec l'acquisition d'image. Cette caractéristique permet d'obtenir des acquisitions moins sensibles aux mouvements. [0061] Le procédé de l'invention comporte une étape visant à générer une première fonction de distribution F, d'une répartition de voxels de l'image en fonction de l'intensité du signal. La résolution du paramètre est configurée selon la modalité d'acquisition de l'image IRM. Selon d'autres modes de réalisation, l'intensité du signal prise en compte est représentée à partir d'un niveau de teinte, de saturation ou de luminosité ou une combinaison de ces trois facteurs. Selon un exemple, l'intensité du signal est un paramètre codé selon une échelle absolue selon l'image acquise. Selon un autre exemple, le paramètre est codé selon une échelle relative en fonction de l'amplitude de variation du paramètre sur l'ensemble des voxels de l'image acquise. [0062] La première fonction F, permet donc de rendre compte de la distribution de cette intensité du signal dans l'image. [0063] Le procédé de l'invention comprend une étape de génération d'une seconde fonction F., permettant de calculer automatiquement un seuil appliqué à la fonction F1 afin de normaliser les valeurs de distribution obtenues. [0064] La seconde fonction F2 est avantageusement une fonction décrivant une répartition de l'intensité du signal de l'image. Selon différentes réalisations, cette fonction peut être 9 mise en oeuvre et définie de manière à optimiser le calcul du seuil, à générer un seuil permettant une normalisation optimisée ou encore en fonction d'un objectif d'identification d'un motif dans l'image ou d'une singularité liée à la répartition de la fonction F1. [0065] Le procédé de l'invention comprend une étape, DET S, permettant de déterminer au moins un seuil automatiquement. Le au moins un seuil déterminé est notamment utilisé pour filtrer la distribution de l'intensité du signal représentée par la première fonction F1 ou par la seconde fonction F2. [0066] Selon différentes réalisations, une pluralité de seuils peut être déterminée. Les seuils peuvent être déterminés par une opération ou une combinaison d'opérations sur au moins la seconde fonction F2. Les opérations peuvent comprendre des calculs d'intersection de courbes, des calculs d'écart-type ou encore la détermination d'une valeur caractéristique, ou encore tout autre type de calculs pouvant être réalisés sur une fonction numérique. [0067] Le procédé comprend une étape SEG visant à segmenter un volume de l'image acquise, en sélectionnant les pixels relatifs à la fonction F,, filtrée avec au moins un seuil précédemment déterminé. [0068] La segmentation comprend la segmentation du volume pulmonaire V, et d'au moins un volume filtré V, Le volume filtré V, peut correspondre à un volume obtenu après filtrage de la fonction F, ou F2. Cette segmentation peut être réalisée en plusieurs étapes. Une première étape peut correspondre à la segmentation du volume pulmonaire V, c'est-à-dire à l'enveloppe du poumon. [0069] Selon un mode de réalisation, le procédé de l'invention comprend une étape de nor- malisation qui permet de normaliser notamment la proportion d'un volume filtré VF segmentée sur le volume pulmonaire V,. [0070] Selon un mode de réalisation, le procédé de l'invention comprend une étape visant à extraire un volume caractéristique du volume filtré qui soit représentatif d'un signe physiologique. Il peut s'agir d'un biomarqueur permettant de quantifier ou qualifier la présence d'une pathologie. [0071] Le procédé permet donc de générer un biomarqueur permettant d'évaluer la proportion d'un volume caractéristique V. [0072] Selon un premier mode de réalisation de l'invention, une mise en oeuvre du procédé permet de générer un premier biomarqueur BIND . Selon une autre mise en oeuvre, soit un second mode de réalisation, le procédé de l'invention permet de générer un second biomarqueur BIND/ - Premier Biomarqueur [0073] La figure 2 représente les principales étapes du premier mode de réalisation de l'invention. Dans cette réalisation, une première configuration CONFI du système d'imagerie par 1RM est réalisée. La première configuration CONE, permet de configurer une pondération de type T1 ou une pondération en densité de proton. Dans les deux cas de pondération, le temps d'écho, noté TE, est défini inférieur à un seuil prédéfini. Le temps de répétition, noté TR, est lui configuré selon une plage de valeurs considérée comme courte ou longue selon la pondération utilisée. L'homme de l'art saura établir l'association entre la configuration TE, TR et la pondération choisie. [0074] Acquisition [0075] Cette acquisition 1RM peut être volumique, c'est-à-dire que l'image obtenue est di- rectement une image tridimensionnelle. Cette configuration CONE, permet de définir un temps d'écho ultracourt, dit "UTE", ou bien une séquence Zéro TE, également appelee ZTE. [0076] La séquence UTE peut être utilisée pour acquérir une image 3D isotrope. Les images 3D peuvent être obtenues, par exemple, au moyen d'une acquisition cylindrique avec des plans radiaux empilés ou une acquisition sphérique en utilisant des projections radiales dans les trois directions. Un choix possible est de prendre les extrémités des rayons sur la surface de la sphère répartis sur des spirales. De tels modes d'acquisition sont décrits par exemple dans la publication « /RN/pulmonaire 3D à tempsdécho ultracourt par acquisition spiralée ou sphérique de l'espace-k », Joseph Yazbek. Le remplissage du volume se fait par le gradient de lecture. On aboutit alors à une surface discrétisée de manière sensiblement uniforme en coordonnées sphériques. [0077] La séquence ZTE est également une technique d'acquisition 3D. cette technique comporte l'application d'un gradient d'amplitude stable et l'envoie d'une impulsion RF ayant une durée très courte et un angle de bascule faible. L'acquisition commence après une durée donnée qui correspond à la bascule entre l'émission et la réception. Entre deux TRs, l'amplitude du gradient de lecture augmente graduellement. Cette configuration présente l'avantage d'être insensible aux courants de Foucault du fait qu'il n'y a pas de contrainte sur la vitesse de variation de l'intensité du gradient pendant la durée d'acquisition. [0078] Un intérêt d'une séquence UTE ou ZTE est de maximiser le signal disponible dans le parenchyme pulmonaire [0079] L'image acquise peut être traitée par des algorithmes de traitements d'image afin de réduire le bruit, supprimer des artefacts ou encore adapter une échelle de contraste. [0080] Selon un mode de réalisation, lorsque l'acquisition est réalisée en 2D et non en 3D, une étape de reconstruction d'une image 3D peut être réalisée par empilement de couches 2D. L'image tridimensionnelle obtenue, et éventuellement traitée, comporte alors une pluralité de voxels. [0081] Traitement d'image [0082] La figure 4 représente une coupe d'une image tridimensionnelle de chaque poumon 11 ayant un volume pulmonaire V. Le volume de la cage thoracique est par exemple ici représenté par le volume V,. Dans cette image est illustré le contour de l'enveloppe de chaque poumon notamment grâce à une discontinuité de l'intensité du signal des voxels délimitant le volume de chaque poumon avec le volume de la cage thoracique. [0083] Selon un mode de réalisation, le système d'imagerie IRM est configuré avec une séquence PETRA. Selon un exemple de réalisation, un filtre rchausseur de contours peut être appliqué à l'image acquise. Un intérêt est d'obtenir une meilleure définition des contrastes de l'image. Le détourage des contours peut être obtenu à partir de l'application d'un filtre local. [0084] Selon un mode de réalisation, le procédé de l'invention comporte une étape de lissage gaussien appliqué à l'image acquise avec réduction du temps de codage ponctuel et avec acquisition radiale pour le débruitage. Un intérêt de ce traitement est de réduire le bruit obtenu lors de l'acquisition par cette méthode de traitement de la régularité supposé de la courbe obtenue. [0085] Extraction du volume pulmonaire [0086] Le procédé de l'invention permet de générer une courbe CI représentant la répartition ou la distribution des voxels selon la valeur de l'intensité du signal ou d'un paramètre colorimétrique. Cette répartition peut être limitée au volume pulmonaire V, en procédant préalablement à l'extraction de l'enveloppe 3D du poumon. Cette extraction peut être réalisée automatiquement par un algorithme de traitement de contours basés sur l'analyse de la valeur de l'intensité du signal ou d'un paramètre colorimétrique des voxels et de leur voisinage. [0087] La valeur de l'intensité du signal ou le paramètre colorimétrique permettant de réaliser une distribution des voxels est avantageusement une donnée de contraste. Tout autre type de paramètre pouvait être utilisé en fonction de l'image acquise. Par exemple, une donnée de luminosité, de teinte ou encore de saturation peut être utilisée en combinaison ou en variante de la donnée de contraste. [0088] Ainsi, une première fonction FI représente la courbe de distribution de voxels dans le volume pulmonaire V, selon une résolution donnée. La résolution donnée con-espond la résolution de l'acquisition de l'image. La résolution est donc dépendante de la configuration d'acquisition de l'image. [0089] La fonction FI peut être une fonction continue ou une fonction discrète. Dans ce dernier cas, les voxels peuvent être rassemblés selon une échelle définissant une unité de contraste. Cette échelle peut être prédéfinie ou automatiquement calculée en fonction, par exemple, d'une valeur minimale et d'une valeur maximale et du volume pulmonaire et/ou en fonction du nombre de voxels considérés. La courbe CI qui est la représentation graphique de la fonction F1 est illustrée à la figure 5. [0090] Selon un mode de réalisation, l'extraction des volumes à l'intérieur du volume 12 pulmonaire V,' tels que le volume filtré V,', ou le volume caractéristique Vc, comprend un traitement d'images permettant d'exclure, voire de supprimer, des voxels d'une même gamme de contraste qui sont déconnectés des voxcls de leur voisinage. Un intérêt est de supprimer les faux positifs, erreurs d'acquisition ou encore les voxcls en dehors de zones d'intérêt représentatives d'un signe physiologique ou pathologique. [0091] Le procédé de l'invention comporte une étape de calcul automatique d'au moins un seuil SI permettant de filtrer les voxcls de la distribution représentée par la fonction FI ou d'une fonction modélisant la fonction FI. [0092] Dans ce premier mode de réalisation, deux seuils SI et S, sont calculés. Le seuil SI définit une valeur maximale et le second seuil définit une valeur minimale de la valeur de l'intensité du signal ou du paramètre colorimétrique considéré pour générer la distribution des voxels. [0093] Définition de la seconde fonction [0094] Selon ce premier mode de réalisation, afin de déterminer des seuils optimaux, le procédé comporte une étape visant à modéliser au moins deux fonctions gaussiennes Fg,, Fg2 par ajustement de la première fonction F,. [0095] Ces deux fonctions gaussiennes Fg,, Fg2 sont obtenues de sorte qu'elles modélisent chacune une répartition du signal du volume pulmonaire V, c'est-à-dire de la distribution des voxels dans le volume pulmonaire V,. [0096] La première fonction gaussienne Fg, modélise essentiellement les signaux du volume pulmonaire V, correspondant à un milieu aérien. La seconde fonction gaussienne Fg2 modélise essentiellement les signaux du volume pulmonaire V, correspondant à un milieu intermédiaire, tel que le coeur, les muscles, les vaisseaux ou la moelle osseuse. [0097] La modélisation par ajustement de la courbe C, permet d'obtenir deux fonctions Fg, et Fg2 représentant deux répartitions gaussiennes. La figure 5 représente les courbes associées aux fonctions Fg,, Fg2. Une troisième courbe Fg% représentée également sur la figure 5 peut être obtenue par ajustement de la courbe Cl. La fonction Fg, modélise essentiellement les signaux du volume pulmonaire V, correspondant à un milieu graisseux. [0098] - Génération des seuils [0099] Le premier seuil S, est généré automatiquement par un calcul de l'intersection des courbes Fg, et Fg2. Le premier seuil SI définit une valeur maximale d'une quantification du signal dans le volume pulmonaire V. Dans le cas de la figure 5, il s'agit d'une valeur maximale d'un paramètre de contraste. La valeur maximale SI permet de filtrer les voxels de la distribution de CI en dessous de la valeur du premier seuil SI. [0100] Un intérêt de la définition du premier seuil SI est d'obtenir la meilleure segmentation entre les valeurs du tissu parenchymateux pulmonaire et les parties molles. [0101] Le second seuil S2 est généré automatiquement en considérant une valeur minimale 13 de contraste définissant les premiers voxcls ayant un premier niveau de contraste. La figure 5 représente le second seuil So. [0102] Un intérêt est de s'affranchir de la variabilité des valeurs absolues du signal de l'image qui peuvent être dépendantes d'un patient à l'autre ou d'une configuration donnée à une autre. Le procédé de l'invention permet de considérer la valeur minimale du signal de l'image comme une valeur normalisée à O. [0103] Une étape de segmentation SEG du volume pulmonaire est alors réalisée dans l'objectif d'extraire le volume filtré V, définissant les voxels de la courbe C1 dont la valeur de l'intensité du signal ou dont le paramètre colorimétrique est compris entre le second seuil Sn (valeur minimale) et le premier seuil SI (valeur maximale). [0104] Le volume du poumon V.' formant l'enveloppe de chaque poumon est également extrait à cette étape s'il ne l'a pas été précédemment. [0105] Normalisation [0106] Le procédé comporte une étape de normalisation NORM permettant de normaliser les valeurs des signaux comprises entre les deux seuils Sn et Si. La valeur minimale correspondant au seuil Sn est définie à la valeur 0 et la valeur maximale correspondant au seuil SI est définie à la valeur 1. Les valeurs à 0 sont proches des valeurs des voxels correspondant à un milieu aérien intra-pulmonaire. Les valeurs à 1 sont proches des valeurs de voxels correspondant à des parties molles thoraciques. [0107] Le procédé de l'invention permet donc de sélectionner un premier seuil S, adaptatif qui permet de quantifier le volume de signaux dont la valeur est inférieure à ce premier seuil SI. [0108] Avantageusement, le procédé permet de prendre en compte une variation possible du niveau du seuil selon les patients. [0109] Acquisition 4D pour le suivi du biomarqueur [0110] Le procédé de l'invention peut être appliqué soit à partir d'une image IRM tridimen- sionnelle statique soit à partir d'une image IRM tridimensionnelle dynamique avec une composante temporelle supplémentaire. On parle alors d'acquisition 4D. Un intérêt d'une acquisition 4D est de générer un marqueur dynamique qui permet d'évaluer la variation temporelle du premier biomarqueur BIO,ND, pendant une durée donnée. Cette acquisition est particulièrement intéressante lorsque la variation de l'indicateur est mesurée entre l'inspiration et l'expiration d'un patient. Le premier biomarqueur BIO [ND permet alors de quantifier une atteinte fonctionnelle des petites voies aériennes. [0111] Définition d'un troisième seuil pour générer le biomarqueur [0112] Les figures 6A et 6B représentent des courbes de distribution de signal d'un volume segmenté à partir d'une image obtenue par un scanner. Dans cet exemple, un seuil S3 prédéfini sur une échelle donnée permet de discriminer un poumon comportant un volume caractéristique Vc représenté à la figure 6A plus important que le volume ca- 14 ractéristique Vc obtenu dans un cas représenté à la figure 6B. Le volume Vc peut correspondre à un milieu intermédiaire corroborant un signe physiologique donné. Les courbes obtenues aux figures 6A et 6B peuvent correspondre à celles de deux patients. La figure 6A est caractéristique de la présence d'emphysème. La figure 6B est caractéristique d'un sujet sain. [0113] Les figures 7A et 7B représentent des courbes de distribution de signaux compris entre les deux seuils S, et S2 d'un volume segmenté caractéristique Vc obtenu par la mise en oeuvre du procédé de l'invention. La figure 7A correspond, par exemple, au même patient que la figure 6A et la figure 6B au même patient que la figure 7B. Les figures 7A et 7B sont obtenues à partir d'un système d'imagerie IRM conformément au procédé de l'invention. [0114] La figure 7A représente un volume caractéristique Vc obtenu pour des voxcls ayant une valeur inférieure à un seuil prédéfini sur une échelle normalisée de 0 à 1. Dans le cas de la figure 7A et 7B, le seuil prédéfini correspond à une valeur de 0,2. Ce seuil permet d'obtenir une bonne représentativité de l'emphysème des patients lorsque les hypo-signaux sont significativement représentés en dessous de cc troisième seuil 5.3. [0115] Lorsque le seuil 8.3, compris entre les seuils S2 et Si, est bas, c'est-à-dire proche de 0, par exemple entre 0 et 0,30, le premier biomarqueur BIO,,,, peut être utilisé dans la détection de signes propres à un phénomène d'emphysème. [0116] Lorsque le seuil 8.3, compris entre les seuils S2 et Si, est un seuil ayant une valeur in- termédiaire, c'est-à-dire une gamme de valeurs définies autour de la valeur médiane ou moyenne de l'échelle [0 ; 1], par exemple entre 0,35 et 0,65, le biomarqueur BIO,,,,, peut être utilisé comme un marqueur du piégeage des voies aériennes. [0117] Ce niveau de premier seuil permet donc de générer un biomarqueur BIO,N,,, permettant de classifier certains patients selon une donnée physiologique. [0118] Un intérêt de la figure 7A vis-à-vis de la figure 6A est que la portion du volume ca- ractéristique Vc obtenu est plus importante avec le procédé de l'invention qu'avec une méthode de l'état de l'art. Par ailleurs la différence entre un sujet saint représenté à la figure 7B et un sujet atteint d'emphysème est plus significative entre les figures 7A et 7B qu'entre les figures 6A et 6B. [0119] Le seuil S:3 défini à 0,2 sur l'échelle normalisée permet d'obtenir des résultats satis- faisants pour le praticien permettant de discriminer la présence d'un volume caractéristique d'un milieu intermédiaire dans le poumon. Second Biomarqueur [0120] La figure 3 représente les principales étapes du second mode de réalisation de l'invention. Dans cette réalisation, une seconde configuration CONF, du système d'imagerie par IRM est réalisée. La seconde configuration CONF2 permet de configurer une pondération de type T2. [0121] L'objectif de ce mode de réalisation est de générer un biomarqueur représentatif d'un volume de charge en mucus et/ou d'oedème au sein du volume pulmonaire. Un avantage de la pondération de type T2 est de représenter le signal de l'image propre aux liquides stagnants, tels que l'eau, le mucus, le sang non circulant. A cette fin, le procédé de l'invention comporte une étape de filtrage de certains volumes pour isoler le volume représentatif de la charge de mucus. [0122] La seconde configuration CONE, de l'acquisition ACQ, peut être avantageusement complétée par une acquisition complémentaire ACQ, d'une image en pondération T, ou une pondération en densité de protons. Cette acquisition complémentaire ACQ, permet notamment d'améliorer la segmentation du volume pulmonaire V,' elle est préférentiellement configurée comme l'acquisition du premier mode de réalisation représenté à la figure 2. [0123] Acquisition [0124] Selon le second mode de réalisation, la pondération T2 est avantageusement réalisée avec un temps d'écho TE long et un temps de répétition TR considéré comme long. On rappelle que cette terminologie TE long, TE court, TR long, TR court est connu de l'homme de l'art qui saurait déterminer une configuration appropriée du système IRM en pondération T2. Une configuration possible correspond à un temps d'écho : TE > 80 ms et un temps de répétition : TR > 2000 ms. [0125] Cette acquisition ACQ2 est, par exemple, réalisée en 2D ou 3D. Dans le premier cas, une étape de reconstruction d'une image 3D peut être réalisée par empilement de couches 2D. L'image tridimensionnelle obtenue, et éventuellement traitée, comporte alors une pluralité de voxels. Cette étape est représentée à la figure 3 par l'étape GEN 3D. Cette étape vise à générer une image tridimensionnelle à partir d'une acquisition 2D. [0126] Acquisition RTSE [0127] Selon un mode de réalisation, l'acquisition en pondération Ti est une séquence d'acquisition radiale dite « turbo spin écho » dans la terminologie anglo-saxonne. Elle est plus communément notée T2-RTSE. [0128] Cette séquence T2-RTSE est une séquence 2D basée sur l'écho de spin avec une tra- jectoire radiale dans l'espace k. Lors de l'acquisition de données radiales, le centre de l'espace k est obtenu naturellement par un effet de moyennage résultant de la configuration de l'acquisition. Cette caractéristique permet d'améliorer les artefacts liés aux mouvements. L'immunité aux mouvements radiaux est renforcée par l'utilisation d'une technique d'acquisition, dite PACE, et désignant dans la terminologie anglo- saxonne : « Prospective Acquisition and CorrEction ». [0129] Cette technique d'acquisition est compatible avec diverses séquences, dont la séquence RTSE. Elle permet notamment de supprimer certains des mouvements ci- 16 nétiques, d'augmenter le ratio signal sur bruit. [0130] L'acquisition T2-RTSE comprend une densité d'échantillonnage selon des anneaux de rayons croissants, les projections radiales des temps d'échos TEs voisins sont alors obtenus par une interpolation linéaire. Ce processus est réitéré jusqu'à ce que le rayon maximal soit atteint, par exemple, en utilisant un algorithme avec un pseudo-rapport d'angle d'or, par exemple décrit dans la publication : « Multi-Echo Pseudo-Golden Angle Stack of Stars Thermometry with High Spatial and Temporal Résolution Using k-Space Weighted Image Contrast ». [0131] Plusieurs échos sont acquis et sont reconstruits à l'aide d'une méthode de hiérar- chisation dans l'espace k qui préserve les contrastes des images acquises à chaque TE. Cette technique comprend ensuite la génération d'une image composite construite à partir de toutes les images acquises sur une pluralité de TE successifs, ainsi que la reconstruction de l'image en pondération T2 à partir d'acquisitions d'écho multiples. [0132] Selon un exemple, les temps d'échos TE successifs se succèdent avec des durées croissantes. Selon un exemple, le signal est acquis sur 5 TE, par exemple ayant des valeurs respectives de : 20 ms, 50 ms, 85 ms, 100 ms et 150 ms. [0133] L'utilisation de longs temps d'écho TE allant jusqu'à 150 ms permet d'améliorer la suppression des signaux provenant des volumes sanguins circulant dans les vaisseaux. Cette fonction de suppression des signaux provenant des volumes sanguins circulant peut être améliorée par une acquisition avec une saturation spatiale et/ou une double inversion-récupération, et/ou un module de diffusions. La double inversion récupération est une séquence IRM, notée DIR et connue de l'homme de l'art. Inversement, les liquides non circulant ayant des TE de durées courtes à des durées intermédiaires pourraient être effacés sur des TE long. La disponibilité de TE courts dans l'image composite reconstruite sur l'ensemble des TE est donc de nature à permettre de ne pas perdre l'information des TE courts. Ainsi, une acquisition multi-TE permet un compromis optimal entre la visualisation efficace des liquides non circulant et l'effacement des parties molles ou les vaisseaux circulants. [0134] Avantageusement, la seconde configuration CONF2 du système d'imagerie IRM comprend un paramétrage des temps TR et TE. Un paramétrage adapté permet de contribuer à une acquisition visant à supprimer les volumes circulants, par exemple les volumes des vaisseaux sanguins. En effet, le paramètre TR permet de contrôler le phénomène d'entrée de coupe des protons et le TE permet de contrôler le phénomène de sortie de coupe des protons. Dans une configuration d'acquisition en écho de spin, le signal acquis baisse rapidement lorsque les vitesses circulatoires augmentent. En effet, les protons circulants ne sont pas rephasés. [0135] Le phénomène est encore plus marqué lorsque l'on paramètre des temps TE longs, car la quasi-totalité des protons excités sont sortis de la coupe lorsque le signal est 17 recueilli. Du fait de l'utilisation d'une séquence d'écho de spin, l'intensité du signal des vaisseaux est rendue faible. Notamment, l'utilisation de TE longs allant jusqu'à 150 ms, ainsi que la saturation spatiale en amont de chaque train d'écho, permettent de supprimer davantage l'intensité du signal des vaisseaux en circulation. [0136] La seconde configuration CONE, comporte, en outre, avantageusement une pré- saturation ou une saturation. Cette technique est fondée, par exemple, sur le choix d'une fréquence émise qui peut être mise en oeuvre par l'émission d'une impulsion radiofréquence. Cette impulsion permet d'annuler l'aimantation longitudinale des protons. L'impulsion est avantageusement émise sous la forme d'une pré-impulsion avant chaque cycle de TE. L'impulsion de radiofréquence disperse les phases des spins sans être suivie par un gradient de re-phasagc mais avec un spoiler qui élimine l'aimantation résiduelle. Par exemple, une bande de pré-saturation placée en amont d'une pile de coupes permet de brouiller l'aimantation des protons susceptible de rentrer dans la coupe. Le signal de ces protons sera donc réduit, voire nul, sans quoi ils seraient susceptibles d'avoir un signal détectable en l'absence de pré-saturation. Cette configuration permet de diminuer ou de supprimer les flux circulants. En conséquence, le procédé de l'invention comporte une étape visant à annuler le signal des protons qui bougent et à conserver le signal des protons qui ne bougent pas. Cette opération peut donc être réalisée grâce à une configuration de la saturation du signal. [0137] Selon un exemple de réalisation, les paramètres suivants peuvent être définis pour établir la seconde configuration CONF2: TE reconstitués = 20, 50, 85, 100 et 150 ms, TR = 2350 ms, angle de bascule = 145 °, taille du pixel = 1,6 x 1,6 mm2, épaisseur de la coupe = 1,6 mm, aucun entrelacement, temps d'acquisition moyen = 12 minutes. [0138] Selon cet exemple, l'image 3D peut être construite dans un second temps par em- pilement d'images 2D. Un volume composé de voxels isotropes, par exemple de dimension : {1,6 mm}8 peut alors être obtenu. [0139] Acquisition complémentaire en pondération T ou densité de protons [0140] Lorsqu'une seconde acquisition IRM est réalisée en pondération T, ou en pon- dération de protons afin de segmenter l'enveloppe du volume pulmonaire Vp, la configuration en pondération T, peut comprendre, par exemple, un temps d'écho : TE = 10 à 20 ms et un temps de répétition : TR = 400 à 600 ms. Cette seconde acquisition IRM peut être volumique, c'est-à-dire que l'image obtenue est directement une image tridimensionnelle. Cette configuration est alors sensiblement proche de la configuration CONE' du premier mode de réalisation. Cette dernière configuration peut être réalisée en définissant un temps d'écho ultracourt UTE ou bien une séquence Zéro TE. [0141] L'acquisition de l'image en pondération TI ou en pondération de protons est alors recalée avec l'image acquise en pondération T2. Les contours du volume pulmonaire V p extraits de l'acquisition en pondération T1 peuvent être recalés sur l'image de 18 l'acquisition en pondération T2 pour segmenter le volume pulmonaire V' de l'image de l'acquisition en pondération T2. [0142] Un avantage de ce rccalage est de permettre de récupérer le contour du volume pulmonaire dont le contraste a été diminué du fait du traitement du volume sanguin circulant. [0143] Selon un autre exemple, une séquence avec TR = 4,1 ms, TE = 0,05 ms et un angle d'inclinaison = 5 peut être définie. [0144] Identiquement, l'image acquise par cette seconde acquisition IRM en pondération T, peut être traitée par des algorithmes de traitement d'image afin de réduire le bruit, supprimer des artefacts ou encore adapter une échelle de contraste. [0145] Extraction du volume pulmonaire [0146] Les images tridimensionnelles obtenues à partir de la première acquisition en pon- dération T, et à partir de la seconde acquisition en pondération T, peuvent être superposées de manière à filtrer les voxels de l'image à l'intérieur du volume pulmonaire V,. Enfin, cette acquisition en pondération TI permet également de détourer le volume du coeur et de le soustraire au volume d'intérêt. [0147] L'objectif de l'acquisition en pondération TI ou en densité de protons est de réaliser un masque de l'enveloppe pulmonaire. Pour réaliser cc masque 3D, l'acquisition 3D UTE est interpolée pour correspondre à la résolution spatiale de l'acquisition T2-RTSE. Selon un exemple, les images ont été enregistrées en utilisant un algorithme optimisé multimodal basé sur l'intensité. Un algorithme itératif pour sélectionner le seuil optimal d'intensité peut être mis en oeuvre. Cet algorithme calcule la moyenne des deux classes séparant l'image en anière-plan et avant-plan. En outre, il calcule de façon itérative le meilleur seuil de séparation et s'arrête une fois que ce calcul ne change pas. Cet algorithme peut être basé sur la division de l'histogramme en deux images selon deux classes d'intensité de voxels afin de créer un masque pulmonaire. Ensuite, le masque pulmonaire est, par exemple, appliqué aux images enregistrées T2-RTSE afin d'isoler les voxels compris dans le volume pulmonaire V. D'autres techniques permettent d'acquérir un masque du volume pulmonaire Par exemple, le masque peut être obtenu à partir d'un autre système d'imagerie ou à partir d'une image précédemment acquise et enregistrée dans une mémoire. [0148] Première fonction F [0149] Le procédé de l'invention comporte une étape visant à générer une première fonction de distribution F, d'une répartition de voxels de l'image acquise par l'acquisition en pondération T, en fonction de l'intensité du signal. Afin d'obtenir une fonction F, qui soit représentative d'une charge de mucus à l'intérieur du volume pulmonaire, une étape de détourage du volume pulmonaire Vp est réalisée, par exemple, à partir d'une seconde acquisition en pondération Tl, comme énoncé précédemment. Par ailleurs, une 19 étape de filtrage des voxcls correspondant au liquide sanguin circulant dans les vaisseaux est également réalisée du fait du mode d'acquisition de l'image en pondération T, comme détaillé précédemment. [0150] La première fonction F, permet donc de rendre compte de la distribution d'un paramètre colorimétrique, tel que l'intensité du signal, dans l'image représentative des volumes de liquide, tels que l'eau. La fonction F, permet donc de discriminer les volumes de mucus, qui sont formés sensiblement de 95% d'eau, ou les volumes stagnants d'oedème, qui sont formés sensiblement de 90% d'eau, dans le volume pulmonaire. L'intensité du signal est préférentiellement un paramètre de contraste. [0151] Une fois la fonction F, générée, une courbe C, représentative de la distribution des voxcls dans le volume pulmonaire est obtenue. La distribution des voxcls est obtenue pour une échelle absolue de voxels dont la résolution dépend de la configuration de l'acquisition de l'image. [0152] Seconde fonction F 2 [0153] Une fonction F2 permettant de définir un seuil automatique permet d'extraire un volume filtré V, et volume caractéristique Vc. On rappelle que le volume filtré est le volume extrait de l'image segmentée correspondant à une zone d'intérêt. Le volume caractéristique Vc est un volume obtenu à partir du volume filtré de manière à générer le hiomarqueur. [0154] A cette fin, un volume de référence VA est identifié dans l'image acquise. Le volume de référence VA est préférentiellement choisi au niveau d'une zone aérienne, c'est-à-dire un volume correspondant à un volume d'air. Le volume de référence VA correspond à une zone de dimension minimale permettant d'être suffisamment représentative d'une zone correspondant à un volume d'air et non à un artefact de l'image ou une zone singulière. En d'autres termes, on choisit une zone suffisamment grande pour obtenir une valeur représentative du niveau de contraste de l'air. Ce volume de référence VA peut être automatiquement sélectionné à une position donnée à l'intérieur de l'enveloppe pulmonaire V, ou à une position donnée à l'extérieure de l'enveloppe pulmonaire V,. Par exemple, une position à une distance prédéfinie de la limite de la surface de la cage thoracique permet d'obtenir une probabilité importante d'obtenir une zone aérienne. La distance peut être par exemple comprise entre 1 cm et 4 cm, préférentiellement entre 1,5 cm et 3 cm. Une distance de 2 cm permet d'être compatible d'un grand nombre de sujets. Le volume de référence VA permet de sélectionner une distribution d'un paramètre colorimétrique, c'est-à-dire par exemple une distribution des valeurs de l'intensité du signal, du volume de référence VA, telles que les valeurs d'une distribution d'un niveau de contraste. La fonction représentant la distribution des voxels de ce volume de référence VA est dans ce mode de réalisation la seconde fonction F,, également notée dans ce second mode de réalisation FA. [0155] Le procédé de l'invention comporte une étape visant à mesurer l'écart-type GA de la seconde fonction FA. L'écart-type nA est calculé sur le volume de référence VA. Le procédé de l'invention comprend ensuite le calcul d'un seuil de référence SA calculé à partir de l'écart-type calculé 0A. Selon un mode de réalisation, le seuil de référence SA est un multiple de l'écart-type 0A. Selon un exemple, la valeur du seuil de référence SA est comprise entre 8 et 16 fois la valeur de l'écart-type ciA à laquelle est ajoutée la valeur du mode principal de la distribution des voxels du parenchyme pulmonaire A cet effet, un volume de hase VB de voxels est automatiquement sélectionné dans la zone du poumon. Ce volume de hase V. peut être prédéfini, par exemple à une distance minimale de l'enveloppe du poumon, ou encore être déterminé autour d'une zone d'intérêt. La zone d'intérêt peut être définie à proximité ou au voisinage dcs bronchioles respiratoires, des conduits alvéolaires et des alvéoles. Selon un cas d'exemple, l'ensemble du volume de chaque poumon est considéré. Selon un autre exemple, différents volumes de hase V. sont calculés en différents endroits et une valeur moyenne est ensuite générée. Selon un autre exemple, le procédé est répété pour différents volumes de hase V.. Le mode principal du volume de hase V. est alors noté ME. [0156] Le mode principal représente la valeur de la plus grande distribution du paramètre co- lori métrique ou des valeurs de l'intensité du signal, il est noté ME. Cc mode principal permet de générer un seuil adaptatif qui dépend du patient et de ses données physiologiques propres. Ceci constitue un réel avantage de l'invention afin de mesurer et quantifier le volume de mucus relativement à un individu. Afin d'obtenir un volume filtré discriminant un volume de charge de mucus, la valeur du seuil de référence SA est comprise entre 10 et 14 fois la valeur de l'écart-type 0A à laquelle on ajoute le mode principal ME. Selon un exemple préféré, la valeur du seuil de référence S' est définie entre 11 et 13 fois la valeur de l'écart-type a' à laquelle on ajoute le mode principal M .. En particulier, la valeur du seuil de référence SA = 12 0A+ ME. [0157] La figure 9 représente un exemple de seuil SA calculé à partir d'un mode principal obtenu à partir de la répartition la plus représentative du signal du volume pulmonaire V,. En ordonnée CI figure l'intensité du signal acquis en Ti et en abscisse figure une métrique de distance propre à l'image en coupe représentée. La distance « Distance (VP) » est exprimée en millimètre et correspond à la distance depuis l'enveloppe du poumon jusqu'à une autre extrémité du poumon. [0158] Selon un premier exemple, la fonction Distance (VP) est définie au moyen d'un segment de droite défini à partir d'un premier point situé sur la paroi du volume pulmonaire, d'une normale à la paroi et d'un second point défini à l'intersection de la paroi pulmonaire opposée au premier point. [0159] Selon un autre exemple, la fonction distance (VP) est définie à partir d'un segment de 21 droite défini à partir d'un premier point situé sur la paroi du volume pulmonaire et qui passe par un centre du volume pulmonaire. Le centre du volume pulmonaire est un centre fictif pouvant être défini à partir d'une fonction géométrique, telle qu'une fonction déterminant un barycentre, appliquée à l'enveloppe du volume pulmonaire Vp. [0160] Selon un autre exemple, la fonction Distance (VP) est définie en considérant un segment compris entre deux points considérés sur la paroi du volume pulmonaire. [0161] Le mode principal MB peut également être calculé au sein d'un volume partiel du volume pulmonaire Vp. [0162] Le signal SN correspond à l'intensité du signal dans une image de coupe (2D). Il est donc mesuré à partir d'une distance de l'enveloppe du poumon. Le seuil SA est ohtenu par un calcul du mode principal MB, ici représenté à titre indicatif sur la répartition des intensités du signal d'une coupe. Afin d'obtenir le seuil SA, il a été ajouté au mode principal un multiple de l'écart-type du signal mesuré sur un volume de référence de l'air pris en dehors du volume pulmonaire, typiquement entre 12 et 16 fois. [0163] Segmentation [0164] Le procédé comporte une étape de segmentation du volume des voxels supérieur au seuil SA. Le volume ohtenu est le volume filtré V,. Ce volume filtré permet de conserver le voxel d'intensité élevé et permet d'éliminer une partie du signal correspondant à du bruit. [0165] Selon un mode de réalisation, la fonction FA, c'est-à-dire la seconde fonction F2, définit un filtre de la fonction FI. Le filtre est appliqué de manière à extraire les voxels filtrés par la fonction FA.c'est à dire ceux pour lesquels les voxels ont un niveau de contraste supérieur au seuil SA. Dans ce cas, la fonction FA définit un masque appliqué à la fonction FI. [0166] Le volume filtré V, est défini par le volume composé des voxels de la distribution de la fonction F, filtrée par le seuil de référence SA. Les voxels du volume filtré ont donc une intensité de contraste supérieure au seuil de référence SA. Un intérêt est d'obtenir la valeur d'intensité en pondération T2 des zones de mucus présentent dans le volume pulmonaire. La fonction F2 agit comme un masque volumétrique filtrant les hautes intensités de signal acquises. [0167] La figure 8 représente un exemple d'une coupe des deux lobes pulmonaires PA, PR dans lesquels sont représentées des zones 10 en 2D correspondant à des parties de volume de zone d'intensité supérieur au seuil SA. Ainsi, la figure 8 représente un exemple d'un volume filtré pal- la fonction R. [0168] Le procédé comprend une étape de normalisation des valeurs de voxels du volume filtré Vp. Selon une première variante de réalisation, les valeurs d'intensité de contraste sont normalisées sur une échelle comprise, par exemple, entre 0 et 1. [0169] Selon une autre variante de réalisation, qui peut être combinée ou réalisée en al- ternative à la première variante, la normalisation comprend le calcul d'un indice d'intensité volumétrique noté VIP. Cet indice VIP est obtenu par un produit entre le volume et l'intensité T, du signal. Il permet de quantifier l'intensité volumique des hypersignaux et donc de quantifier des niveaux d'inflammation. Il permet, en outre, de quantifier l'évolution d'une zone par la corrélation entre l'augmentation ou la diminution d'un volume de zone inflammée et l'augmentation ou la diminution d'une intensité de signal. [0170] La normalisation comprend donc une normalisation des valeurs absolues des pa- ramètres colorimétriques ou de l'intensité du signal des voxels, ici exprimés par la quantification du signal T2, par exemple exprimé en ms et d'un seuil de filtrage correspondant au masque du volume filtré appliqué, grâce au seuil SA. [0171] Selon un exemple, l'indice VIP peut s'exprimer en millisccondcs et résulter de l'équation suivante : [0172] VIP E max T2ixHSVi VP [0173] - i correspond à une valeur de T2, à cette fin une distribution de valeurs du signal en T2 est réalisée ; max correspond à la valeur maximale du signal en T2. T2i représente le signal de niveau i de l'image tridimensionnelle ; HSV; correspond au volume associé au niveau de valeur de T2 de la distribution, ledit volume étant compris dans le masque de l'image tridimensionnelle correspondant au volume filtré VF des hautes intensités de signaux en T2; Vp correspond au volume pulmonaire. 101741 Un avantage est (l'obtenir une valeur normalisée par un indice reflétant une intensité par unité de volume. Compte tenu du seuil adaptatif obtenu grâce au seuil S, le procédé de l'invention permet la génération d'un biomarqueur qui soit exploitable en limitant la variabilité de l'exploitation et de l'interprétation des images d'un individu à un autre. Par ailleurs, ce biomarqueur permet d'obtenir des points de références temporels pour mesurer l'évolution clans le temps du développement du mucus. 101751 Selon un mode de réalisation, une étape comportant la fusion des images UTE et RTSE est réalisée. Une fusion des images morphologiques à haute résolution obtenues à l'aide d'une séquence à temps d'écho ultra court UTE, et le masque du volume des hypersignaux T, segmentés, est dans ce cas réalisée afin de combiner les deux informations. Cette technique peut être appliquée par l'homme de l'art en PET-TOM. 23 REVENDICATIONS[Revendication 1] Procédé de génération d'un biomarqueur comportant : acquisition (ACQ) d'une image par un système IRM ; traitement de ladite image IRM pour générer (GEN 3D) une image tridimensionnelle du poumon ; génération (GEN FI) d'une première fonction (F,) correspondant à la distribution des différentes valeurs dc l'intensité du signal dc chaque voxcl d'une portion de l'image tridimensionnelle acquise ; calcul automatique (DET S) d'au moins un seuil dc filtrage (S A, SI, S2) de ladite première fonction (F1) à partir d'au moins une seconde fonction (F2, En, Fo, Fg2) de distribution de différentes valeurs d'intensité du signal de chaque voxel d'une portion de l'image tridimensionnelle acquise ; segmentation (SEG) d'un volume comportant : - Un volume principal (V,), correspondant au volume pulmonaire ; - Un volume filtré (V,.) d'un volume de voxels quantifiés par la première fonction (F1) et filtré au moyen d'au moins du seuil de filtrage (SA, SI) calculé normalisation des valeurs de l'image tridimensionnelle du volume pulmonaire (V,i) à partir des valeurs absolues des valeurs d'intensité du signal des voxels de l'image et d'au moins le seuil de filtrage calculé (SA, S,) ; génération d'un biomai-queur (BLND) indiquant un ratio de volume segmenté normalisé. [Revendication 2] Procédé de génération d'un bien-largueur selon la revendication 1, ca- ractérisé en ce que : l'acquisition (ACQ) d'une image tridimensionnelle par un système IRM est configurée par : une pondération en T, un temps d'écho (TE) supérieur à seuil prédéfini le calcul automatique d'au moins un seuil de filtrage (SA) comprend : - acquisition d'un volume de référence ; - - 24 génération d'une fonction de référence (FA) correspondant à une distribution des différentes valeurs de l'intensité du signal de chaque voxcl du volume de référence (VAIR), ladite fonction de référence (FA) étant la seconde fonction de distribution (FA) ; calcul de l'écart-type (cr,,) de la fonction de référence (FA) ; détermination d'un seuil de filtrage, dit seuil de référence (SA) à partir de l'écart-type calculé de la fonction de référence (FA). [Revendication 3] Procédé de génération d'un biomarqueur selon la revendication 2, caractérisé en ce que l'acquisition est paramétrée et comprend : une acquisition d'une image recomposée d'une pluralité d'images acquises sur un nombre de cycle d'au moins 4 temps d'écho (TE), les temps d'échos étant configurés selon des durées croissantes ; une séquence Echo de spin ; un paramétrage visant à émettre un signal pour pré-saturer ou une saturer du signal acquis. - - [Revendication 4] [Revendication 5] [Revendication 6] Procédé de génération d'un biomarqueur selon l'une quelconque des revendications 2 à 3, caractérisé en ce qu'une seconde acquisition (ACQ2) de l'image est réalisée avec une première configuration (Conf.!) pour détourer le volume pulmonaire (Vp), la première configuration (CONF1) définissant un paramétrage d'une acquisition en pondération TI ou en densité de protons, une étape de traitement d'image étant réalisée pour combiner l'image acquise en pondération T, avec l'image acquise par la seconde acquisition (ACQ2). Procédé de génération d'un biomarqueur selon la revendication 4, caractérisé en ce qu'une opération de fusion entre au moins une image acquise en pondération T1 avec un temps d'écho ultra court UTE et au moins une image acquise par une pondération T2 est réalisé pour générer une image dont les données provenant de chacune des images acquises ont été combinées pour segmenter le volume pulmonaire Procédé de génération d'un biomarqueur selon l'une quelconque des re- 25 vendications 2 à 4, caractérisé en ce que le seuil de référence (SA) est établi à partir d'une combinaison entre une valeur de répartition de référence de la fonction de référence (FA) et une valeur comprise entre 10 et 20 fois la valeur de l'écart-type (6,) de la fonction de référence (FA )- Procédé de génération d'un biomarqueur selon la revendication 6, caractérisé en cc que la valeur de répartition de référence de la fonction de référence (FA) est le mode principal (MB) de la répartition des valeurs d'intensité du signal de l'image acquise en pondération T, au sein du volume pulmonaire (V'). Procédé de génération d'un biomarqueur selon l'une quelconque des revendications 2 à 7, caractérisé en ce que le procédé comprend une étape de normalisation, la normalisation comprenant le calcul d'un indice volumique du signal (VIP) à partir des valeurs absolues du signal du volume filtré (VE), du volume résultant du volume filtré (y,) et du volume pulmonaire (Vp). Procédé de génération d'un biomarqueur selon la revendication 1 caractérisé en cc que : l'acquisition (ACQ) d'une image tridimensionnelle par un système IRM, est configurée par : une pondération en densité de proton (TR long ou TE long) ou en Ti(TR court, TE court) ; un temps d'écho inférieur à seuil prédéfini ; le calcul automatique d'au moins un seuil de filtrage (SI, S,) comprend : modélisation (MODEL) d'au moins deux fonctions gaussiennes (Fgl, Fg2) par ajustement de la première fonction (F1) ; détermination du seuil de filtrage, dit premier seuil, par un calcul de l'intersection de la première fonction gaussienne (Fg1) et de la seconde fonction gaussienne (Fg2) ; détermination d'un second seuil (S2) correspondant à une valeur minimale de la première fonction gaussienne (Fg1) et une valeur minimale d'une valeur de l'intensité du signal d'un voxel ; [Revendication 7] [Revendication 8] [Revendication 9] - - - 26 le volume filtré (VF) est un premier volume (V1) correspondant aux voxcls quantifiés (Q) par la première fonction gaussicnne (F81) comprise entre le premier seuil (S1) et le second seuil (S2 ), lesdits voxcls correspondant à un milieu aérien, normalisation des valeurs de l'image tridimensionnelle du volume pulmonaire (Vo) à partir du premier seuil (Si) et du second seuil (S2) calculés (SI, S2) ; génération d'un premier biomarqueur (B,,,i) indiquant un ratio d'un volume caractéristique (Va du volume segmenté normalisé, ledit ratio étant calculé entre le volume caractéristique (Va) et le volume pulmonaire (VF). [Revendication 10] [Revendication 11] [Revendication 12] [Revendication 13] [Revendication 14] Procédé de génération d'un biomarqueur selon la revendication 9, caractérisé en ce que le temps d'écho est inférieur à lms. Procédé de génération d'un biomarqueur selon l'une quelconque des revendications 9 à 10, caractérisé en ce que la segmentation (SEG) comprend la définition d'un second volume (V2) correspondant aux voxels quantifiés (Q) par la seconde fonction gaussienne (Fg2) supérieure au premier seuil (S1), lesdits voxels correspondant à un milieu intermédiaire ou graisseux. Procédé de génération d'un biomarqueur selon l'une quelconque des revendications 9 à 11, caractérisé en ce que la segmentation (SEG) comprend une étape d'extraction d'un volume caractéristique (Va) comportant des voxels du volume filtré (VF) dont la valeur d'intensité du signal est inférieure à un troisième seuil (53) prédéfini, ledit troisième seuil (S3) prédéfini étant déterminé sur une échelle normalisée de [0 ; Procédé de génération d'un biomarqueur selon la revendication 11, caractérisé en ce que la segmentation (SEG) comprend une étape d'exclusion/suppression des voxels déconnectés à leur voisinage de voxels d'une même quantification (QI, Q. Q). Procédé de génération d'un biomarqueur selon la revendication 9, caractérisé en ce que la modélisation des fonctions gaussiennes (Fgi, Fg,) comprend : un lissage gaussien appliqué à l'image acquise avec réduction du temps de codage ponctuel et avec acquisition radiale pour 27 le débruitage ; détourage des contours par application d'un filtre local ; utilisation de la méthode d'ajustement de courbe représentant b fréquence de distribution de voxels. [Revendication 15] Procédé de génération d'un biomarqueur selon l'une quelconque des re- vendications 1 à 14, caractérisé en cc que le procédé comprend l'acquisition (ACQ) des valeurs de l'intensité du signal de chaque voxel quantifiées par la première fonction (F1), les valeurs d'intensité correspondant à une donnée de contraste d'image. [Revendication 16] Procédé de génération d'un biomarqueur selon l'une quelconque des re- vendications 1 à 14, caractérisé en ce que l'acquisition (ACQ) est réalisée de manière synchronisée avec un respirateur, le respirateur étant un navigateur ou une ceinture respiratoire. [Revendication 17] Procédé de génération d'un biomarqueur (13,,,,,) selon l'une quelconque des revendications 1 à 16, caractérisé en ce qu'une étape d'extraction d'une image volumétrique est réalisée à partir des images acquises par l'IRM, ladite image extraite étant réalisée à un instant déterminé de la séquence. [Revendication 18] Procédé de génération d'un biomarqueur (13,,,,,) selon l'une quelconque des revendications 1 à 17, caractérisé en ce que l'acquisition (ACQ) de l'image tridimensionnelle est réalisée au moyen d'un empilement d'images 2D acquises, l'épaisseur de la coupe étant égale au moins à la largeur d'un voxel. [Revendication 19] Système comportant au moins un calculateur et une mémoire et une interface pour recevoir une image acquise par un système IRM, ledit système étant configuré pour mettre en oeuvre les étapes du procédé de l'une quelconque des revendications 1 à 18.