Brevet FR2925707A1 — publié le 20090626

L'invention concerne un procédé d'extrapolation d'un masque bidimensionnel ainsi qu'un procédé d'estimation d'un masque digital tridimensionnel à partir d'une série de masques bidimensionnels extrapolés, procédés destinés à être mis en oeuvre dans un dispositif de radiographie, notamment du type tomosynthèse. Dans n'importe quelle application de la tomosynthèse, il peut être intéressant de pouvoir extrapoler un masque bidimensionnel d'un objet projeté au niveau de vues tronquées de l'objet, puis de pouvoir estimer un masque tridimensionnel à partir d'une série de masques bidimensionnels ainsi extrapolés. En effet, cette extrapolation des masques bidimensionnels est utile lorsqu'il est nécessaire de connaitre une estimation du signal en dehors des limites physiques du détecteur durant la reconstruction. Le masque tridimensionnel peut être utilisé comme information a priori dans les méthodes itératives de reconstruction, afin d'accélérer la reconstruction en se focalisant seulement sur les voxels appartenant à l'objet considéré, ou afin d'éliminer les artefacts de reconstruction situés au delà de l'objet considéré. Les principales difficultés d'une reconstruction fiable du masque tridimensionnel proviennent de l'intervalle angulaire limité parcouru par la source et la troncation de l'objet considéré dans certaines projections du fait des dimensions finies du détecteur. Par exemple, dans le cas de la tomosynthèse numérique du sein, cette tomosynthèse est une nouvelle technique d'imagerie du sein par tomographie tridimensionnelle limitée angulairement. Elle permet de contourner substantiellement le problème de la superposition dans la détection de lésions. Plusieurs vues projetées selon différents angles d'acquisition réduisent potentiellement le nombre de faux positifs dus à l'addition d'artefacts ainsi que le nombre de faux négatifs dus à l'effet de masquage des tissus recouvrants. Durant l'examen de tomosynthèse, une série d'images est acquise pour différents angles d'acquisition dans le but de reconstruire une représentation tridimensionnelle du sein. Pendant la procédure de reconstruction, savoir si, en un voxel v donné de la représentation tridimensionnelle, il y a ou il n'y a pas une information propre au sein peut aider à réduire les artefacts ainsi que le temps de calcul. Cette connaissance est représentable à l'aide d'un masque tridimensionnel du sein M : M[v]=1 si v appartient au sein et M[v]=O sinon. Ce masque M étant reconstruit à partir des masques bidimensionnels issus de la série d'images précédemment acquises, cela suppose que tous les voxels de la région d'intérêt se projettent dans les limites du détecteur. Mais en réalité il y a des voxels de la région d'intérêt qui se projettent en dehors des limites du détecteur lors de l'acquisition de certaines images de la série d'images acquises. De ce fait, il est alors nécessaire d'estimer le masque bidimensionnel du sein au-delà des limites du détecteur. Sans information a priori, l'approche naturelle est de considérer que tous les pixels en dehors des limites du détecteur appartiennent au sein. Ainsi, les masques bidimensionnels et, par conséquence, tridimensionnel ne sous-estiment jamais la forme du sein, car ils sont toujours plus grands que cette forme. Dès lors, lorsque le masque tridimensionnel est construit, des artefacts apparaissent : marches d'escalier. Le contour de masque résultant est alors une mauvaise estimation de ce que cela devrait être et la forme du masque n'est pas continue. Par conséquent, l'utilisation du masque tridimensionnel ainsi obtenu sur la représentation tridimensionnelle produit une ligne de peau non naturelle dans la représentation tridimensionnelle finale, ce qui complique le diagnostique que le praticien effectue sur cette base. Un but de l'invention est de fournir un procédé d'extrapolation de masques bidimensionnels d'un objet dans les vues projetées tronquées (afin que le masque ainsi extrapolé selon une projection particulière est régulière et consistante avec les masques bidimensionnels selon les autres projections) et de reconstruction d'un masque tridimensionnel de l'objet acquis qui ne sous-estime pas l'objet et qui est régulier, au sens où la troncation des vues n'introduit aucune discontinuité dans le masque. A cet effet, il est prévu, selon l'invention, un procédé d'extrapolation d'un masque bidimensionnel Me2 dans un dispositif de radiographie du type comportant une source de rayons X mobile prenant au moins deux positions Si eet Se2 d'un espace associées à des masques bidimensionnels Met et Me2 respectivement, des moyens d'enregistrement sensiblement plans disposés en regard de la source et comportant une limite (d, d'), et un volume d'intérêt, comprenant un objet apte à être radiographié, situé entre la source et les moyens d'enregistrement, le procédé comprenant des étapes de : a) estimation au moins partielle d'une série de masques bidimensionnels My associés à une série de positions Sy de la source situées entre les positions Sel et Set de la source, à partir des masques bidimensionnels Met et M(u, b) évaluation au moins partielle d'un masque tridimensionnel intermédiaire de l'objet à partir de la série de masques bidimensionnels My et des masques bidimensionnels Met et Me2r et c) extrapolation du masque bidimensionnel M92 au-delà de la limite d ou d' des moyens d'enregistrement selon une position relative entre Sel et Se2 à partir du masque tridimensionnel intermédiaire. Avantageusement, mais facultativement, le procédé selon l'invention comporte en outre au moins l'une des caractéristiques suivantes : préalablement à l'étape a), si le masque bidimensionnel Me1 associé à la position de la source Sel n'est pas disponible, le masque bidimensionnel Me1 est alors extrapolé à partir de masques bidimensionnels Me disponibles, - l'extrapolation de l'étape c) comporte une projection du masque tridimensionnel intermédiaire à partir de la position Se2 sur un plan passant par les moyens d'enregistrement, l'estimation au moins partielle de l'étape a) s'effectue sur la limite d ou d' des moyens d'enregistrement selon la position relative entre Sel et S92, - l'estimation au moins partielle de l'étape a) comporte, pour chacun des masques bidimensionnels My, une étape de détermination d'un point Ty situé sur un bord de l'objet projeté sur la limite d ou d' des moyens d'enregistrement, la source occupant la position Sy, - les points Ty sont estimés par une interpolation linéaires entre des points Tel et Tee situés sur un bord de l'objet projeté sur la limite d ou d' des moyens d'enregistrement (10), la source (S) occupant les positions Sel et Set respectivement, - si le masque bidimensionnel Me1 associé a la position Sel de la source n'est pas disponible, le point Tel est alors extrapolé à partir de points Te disponibles, - le procédé comporte une étape supplémentaire de : d) application d'une fermeture sur le masque bidimensionnel Met extrapolé, - préalablement à l'étape a), le procédé comporte une étape de détermination d'une limite d'extrapolation lfin sensiblement parallèle à la limite d ou d' des moyens d'enregistrement et située en dehors des limites des moyens d'enregistrement, - l'étape b) comporte, pour chaque ligne 1 sensiblement parallèle à la limite d'extrapolation lfin et située entre la limite d ou d' des moyens d'enregistrement et la limite d'extrapolation lfin, des sous étapes de : b1) estimation d'un plan P1e2 passant par la position Set et la ligne 1, b2) estimation du masque tridimensionnel intermédiaire pour chaque voxel (v) situé à 30 l'intersection du plan P1e2 et du volume d'intérêt, et - l'étape c) comporte une étape de projection sur chaque ligne 1 sensiblement parallèle à la limite d'extrapolation lfin et située entre la limite d ou d' des moyens d'enregistrement et la limite d'extrapolation lfin du masque tridimensionnel intermédiaire. Il est aussi prévu, selon l'invention, un procédé d'extrapolation d'une série de masques bidimensionnels Me dans un dispositif de radiographie, caractérisé en ce que chaque masque bidimensionnel Me de la série de masques bidimensionnel est extrapolé par un procédé présentant au moins l'une des caractéristiques précédentes. Avantageusement, mais facultativement, le procédé selon l'invention comporte en outre au moins l'une des caractéristiques suivantes : le procédé présentant au moins l'une des caractéristiques précédentes est appliqué de manière 20 itérative, et - à chaque itération, les moyens d'enregistrement sont étendus en leur adjoignant des moyens d'enregistrement virtuels correspondant à une partie commune à tous les masques bidimensionnels Medéjà extrapolés. 25 Il est prévu aussi, selon l'invention, un procédé d'estimation d'un masque digital tridimensionnel à partir d'une série de masques bidimensionnels Me dans un dispositif de radiographie du type comportant une 30 source de rayons X mobile prenant au moins deux positions Si eet Set de l'espace associées à des masques bidimensionnels Met et M(u respectivement, des moyens d'enregistrement sensiblement plans disposés en regard de la source et comportant une limite d ou d', et un volume d'intérêt, comprenant un objet apte à être radiographié, situé entre la source et les moyens d'enregistrement, le procédé comprenant des étapes de : a) extrapolation de la série de masques bidimensionnels Me par un procédé présentant au moins l'une des caractéristiques précédentes ; et b) détermination du masque digital tridimensionnel à partir de la série de masques bidimensionnels Me extrapolés. Avantageusement, mais facultativement, le procédé selon 15 l'invention comporte en outre au moins l'une des caractéristiques suivantes : -l'étape b) comporte des sous-étapes de : b1) pour chaque masque bidimensionnel Me, application d'une dilatation suivie d'un filtre passe-bas 20 pour obtenir une fonction d'appartenance pMe ; b2) évaluation d'une fonction d'appartenance }1M3d à partir de la fonction d'appartenance }2Me en utilisant un opérateur t-norme ; b3) détermination du masque digital tridimensionnel à 25 partir de la fonction d'appartenance }1M3d, et - l'opérateur t-norme est de type probabiliste. Il est aussi prévu, selon l'invention, un dispositif de radiographie du type comportant : 30 • une source de rayons X mobile le long d'un arc de cercle de centre C, • des moyens d'enregistrement sensiblement plans disposés en regard de la source, • un volume d'intérêt, comprenant un objet apte à être radiographié, situé entre la source et les moyens d'enregistrement, • des moyens de mise en oeuvre d'un procédé présentant au moins l'une des caractéristiques précédentes. D'autres caractéristiques et avantages de l'invention apparaîtront lors de la description ci-après d'un mode de réalisation de l'invention et de variantes. Aux dessins annexés . la figure 1 est une vue schématique tridimensionnelle d'un dispositif et du procédé selon l'invention ; la figure 2 est une vue schématique de coté du dispositif et du procédé de la figure 1 permettant l'extrapolation sur un segment 1 d'un masque bidimensionnel Met ; la figure 3 est une vue schématique de coté similaire à la figure 2 mais ne permettant pas l'extrapolation du masque bidimensionnel M92 sur le segment 1 ; la figure 4 est une vue schématique de coté similaire aux figures 2 et 3 illustrant le procédé selon l'invention pour l'extrapolation d'un masque bidimensionnel d'extrémité de la série de masques bidimensionnels à extrapoler par le procédé selon l'invention ; et, la figure 5 est une schématique d'une étape de détermination d'un masque digital tridimensionnel par le procédé selon l'invention. En référence à la figure 1, un dispositif de radiographie apte à réaliser de l'imagerie par tomographie tridimensionnelle comporte des moyens d'enregistrement 10 qui se présentent sous la forme d'un détecteur de type numérique sensiblement plan et définissant un plan PD. Un tel détecteur numérique comprend une matrice de capteurs formant chacun un pixel p uniformément répartis en lignes et colonnes. Le dispositif de radiographie comporte en outre une source de rayons X en regard des moyens d'enregistrement et mobile par rapport à ces mêmes moyens d'enregistrement. La source de rayon X, qui se présente généralement sous la forme d'un tube à rayons X, est apte à se déplacer de manière discrète selon une trajectoire S sensiblement en arc de cercle de rayon r et de centre C. La trajectoire S est limité angulaire dans un intervalle de fonctionnement [emini emax] , où 0 inférieur PB et un plan supérieur PT. Les plans PB et PT sont sensiblement parallèles au plan PD du détecteur formant les moyens d'enregistrement 10. Le volume d'intérêt est d'autre part délimitée par quatre plans (non représentés) passant chacun par des bords d, d', e, e' du détecteur et sensiblement perpendiculaires aux plans PB et PT. Ainsi la région d'intérêt 20 a la forme d'un parallépipède rectangle, comme cela est illustré en figure 1. Dans le cas d'un dispositif pour la tomosynthèse numérique du sein, la région d'intérêt 20 est délimitée par une pelote de support du sein dont une face supérieure est dans le plan PB et une face inférieure dans le plan PD, un plateau de compression dont une face inférieure est dans le plan PT, ainsi que les plans sensiblement parallèles au tronc de la patiente et passant par les première e et dernière e' lignes de pixels du détecteur et les plans sensiblement perpendiculaires aux précédents et passant par les première d' et dernière d colonnes de pixels du détecteur. D'une manière classique, un masque digital tridimensionnel Mid de l'objet situé dans le volume d'intérêt 20 est décrit de la manière suivante : Mid [v] = J 1 si v est dans l'objet 0 sinon pour tout voxel v du volume d'intérêt 20. De là, un masque bidimensionnel Me est une projection de ce masque Mid correspondant à une position Se de la source sur le plan PD du détecteur 10, et est défini par : 1 si p appartient à l' aire où l' objet est projeté Me[p] = 0 sinon Pour tout pixel p du plan PD du détecteur. Soit proj:[0;2r]xR3-PD l'application qui associe à un voxel v de l'espace R3, comprenant le volume d'intérêt 20, et à un angle e la projection sur le plan PD correspondant à la position Se de la source. Nous avons alors . VpEPD,Vvtelquep=proj(B,v) M3d[v]=1~Me[p]=1 Me[p] = 0 M3d[v] = 0 ce qui signifie que, pour une position de la source S donnée, tous les voxels se projetant en dehors d'une région de l'objet sur le plan PD du détecteur n'appartiennent pas au masque digital tridimensionnel de l'objet en question. Notations : un référentiel 0 utilisé dans la suite de la description est décrit par le plan PD du détecteur et le plan qui lui est perpendiculaire, qui passe par la première ligne e de pixels du détecteur et qui contient la trajectoire de la source. Ainsi, un axe x est dirigé selon les colonnes de pixels du détecteur et orienté de la première ligne e vers la dernière ligne e'. Un axe y est situé sur la première ligne e et orienté de la dernière colonne d vers la première colonne d'. L'origine du référentiel est située dans le plan PD de sorte que le centre C de la trajectoire de la source a pour coordonnées (0,0,z0). dans ce référentiel, chaque voxel v de l'espace R3 a pour coordonnées (xv, yV, zV) . VOE [O;TC], Se est la position de la source sur la trajectoire d'angle e par rapport à l'axe y. Ainsi, VGE [O;,r] xse = 0, yse = rcos(9), zse = rsin(9)+zc yd est l'ordonnée de la dernière colonne d, yd, est celle de la première colonne d'. VOE [O;2r], Te est le point situé sur la dernière colonne d (resp. la première colonne d') du détecteur, point correspondant à la projection associée à l'angle e d'un bord de l'objet situé dans le volume d'intérêt (20). VGE [O;TC] les points Ee et Je sont les points d'intersection respectifs du segment de droite SeTe et des plans PT et PB. Pour une colonne 1 de pixel du détecteur, A1= (xe=O, yl, O) est le point correspondant au premier pixel du détecteur dans la colonne 1 et B1= (xe,, yl, O) est le point correspondant au dernier pixel du détecteur dans cette colonne 1. Soit Ple le plan passant par les trois points Ser Al et B. Description du procédé selon l'invention : 25 Nous allons décrire le procédé selon l'invention. Comme les points Te sont sur le bord de l'objet projeté au niveau de la dernière colonne d (resp. la première colonne d'), les segments de droite SeTe sont tangents au bord de l'objet du masque digital tridimensionnel 30 Mid. Comme le volume d'intérêt 20 est limité par les plans PT et PB, cela veut dire que les segments de 15 20 droite E9J9 sont tangents audit bord de l'objet. L'idée du procédé selon l'invention est d'utiliser ces segments pour définir un masque tridimensionnel intermédiaire M enveloppant l'objet et de projeter ce masque tridimensionnel intermédiaire au-delà des limites du détecteur représentés par la dernière colonne d (resp. la première colonne d') de manière à extrapoler les masques bidimensionnels M9, colonne de pixels par colonne de pixels. Supposons que nous souhaitons extrapoler un masque bidimensionnel M92, associé à un angle e2 et une position discrète S92 de la source, au niveau d'une colonne 1 de pixels en dehors des limites du détecteur. Nous avons alors yl yd,>0) dans le référentiel 0 de la figure 1. Supposons un angle e1>e2 (resp. e1 L'extrapolation du masque bidimensionnel M92 sur la colonne 1 de pixels est effectuée en projetant le masque tridimensionnel intermédiaire M sur le plan PD du 10 détecteur 10 en considérant la source située à la position S. Dans le cas d'espèce, nous estimons, dans une première étape, la dernière colonne d (resp. la première colonne 15 d') dans un masque bidimensionnel My associé à la position Sy de la source d'angle y compris entre e1 et e2 . Vp E [AdBd], My[p] = Si 01 xp il existe un angle y compris entre e1 et e2 tel que le voxel v se projette sur la dernière colonne d (resp. la première colonne d') lorsque la source est à une position S. Le point SY est l'intersection entre une droite passant par Ad (resp. Ad') et le point (0, yV, zV) et la trajectoire S de la source. Toutefois, il y a une condition nécessaire pour que l'angle y existe tel que tous les voxels v de l'intersection du volume d'intérêt 20 et du plan p1e2 puissent se projeter sur la dernière colonne d (resp. la première colonne d') . yc > max(yproj(ez,EB, ), Y prof ez J,)) (resp. Yl Si cette condition n'est pas remplie, nous sommes dans une situation similaire à celle illustrée en figure 3. Une partie seulement des voxels v de l'intersection considérée appartiennent à la partie 21 du volume d'intérêt 20. De ce fait, le masque tridimensionnel intermédiaire fi ne peut pas être estimé sur l'ensemble du volume d'intérêt 20 car il manque des informations du fait des voxels v de l'intersection considérée qui sont en dehors de la partie 21 du volume d'intérêt 20 qui ne sont pas � vus � par les moyens d'enregistrement 10 que la source soit à la position Set ou à la position Sel. Donc il existe une colonne lfin au- delà de laquelle il n'est pas possible d'estimer le masque bidimensionnel Met. Cette colonne lfin est selon une droite d'équation Y = MaX (Y prof (82 )' Pm] 9z J~1) (resp. y =min(yproj(92,Ee,)'Yproj(92,Je,)/) De là, pour toute colonne 1 de pixels comprise entre la dernière colonne d (resp. la première colonne d') et la colonne lfin, quelque soit le voxel v appartenant à l'intersection du plan Ple2 et du volume d'intérêt 20, le masque tridimensionnel intermédiaire Al est déterminé par . M[v] _ 1 siM7[proj(y,v)]=1 0 sinon De là, dans une troisième étape, nous projetons le masque tridimensionnel Al sur la colonne 1 de pixel en utilisant la position S92 de la source S : J1 si v e volume d'intérêt (20) tel que proj (62, v) = pet M (v) = 1 0 sinon Ceci permet d'obtenir l'extrapolation du masque bidimensionnel M92 à la colonne 1 de pixels. Extrapolation pour les masques bidimensionnels Mez 15 associés aux projections d'extrémité de la trajectoire S de la source Dans cette situation, l'angle e2 est égal soit à emax soit à emin• Il n'existe donc pas d'angle e1>e2=emax (resp. e1 e2=emax (resp. e1 e3 (resp. e2 position du point T'e1 sur la dernière colonne d (resp. la première colonne d') en utilisant une interpolation linéaire entre Tee et Te3r comme précédemment pour estimer le point Ty à partir des points Tel et Tee. Il suffit alors d'appliquer le procédé selon l'invention précédemment décrit. Formes réelles d'objets aptes à être présent dans le volume d'intérêt 20 : Pour décrire le procédé selon l'invention, nous avons supposé l'existence d'un seul point Te sur la dernière colonne d (resp. la première colonne d') . Toutefois, la plupart du temps, ce n'est pas le cas : il y a plusieurs points Te sur cette colonne d (resp. d'). Afin de palier à ces cas, le procédé selon l'invention précédemment décrit est appliqué en considérant un point Te à la fois donnant un résultat intermédiaire associé. S'il y a, par exemple, trois points Te, le procédé selon l'invention est appliqué trois fois et il y a trois résultats intermédiaires associés. Une fois tous les points Te considérés les uns après les autres, le résultat final est l'application d'un � ET � logique sur l'ensemble des résultats intermédiaires associés. Dans une variante de réalisation, le � ET � logique est 25 appliqué lors de l'évaluation de l'intersection entre le volume d'intérêt 20 et le plan Piea Une fois tous les masques bidimensionnels extrapolés, la limite lfin est en mesure d'être repoussée en considérant 30 un détecteur virtuel couvrant une zone commune entre tous les masques bidimensionnels extrapolés et en itérant le procédé d'extrapolation selon l'invention qui vient d'être décrit précédemment. Une fois la série de masques bidimensionnels M92 ainsi extrapolés, une étape finale peut être appliquée sur ces masques afin d'être sûr qu'ils soient aussi naturel que possible au regard de l'objet radiographié dans le volume d'intérêt 20. Cette étape finale est l'application d'une fermeture morphologique sur la base d'un disque ayant un diamètre en pixels prédéterminé. Une telle fermeture morphologique consiste en une dilatation basée sur ledit disque puis une érosion basée sur le même dit disque. Cette méthode est connue en soit en imagerie. Détermination du masque digital tridimensionnel final : Une reconstruction du masque digital tridimensionnel final peut alors être réalisée à partir de la série de masques bidimensionnels M92 ainsi extrapolés. Une telle reconstruction est connue en soit. Toutefois, afin d'éviter des transitions objet/arrière-plan brutales, une pondération (appelé aussi selon l'anglicisme � fuzzification �) est appliquée est utilisant la théorie des ensembles flous. Pour cela, des fonctions d'appartenance sont associées aux propriétés du masque digital tridimensionnel M3d (M3d[v]=1 si le voxel v appartient à l'objet considéré) et du masque bidimensionnel Me (Me[p]=1 si le pixel p appartient à la surface ou l'objet considéré est projeté sur le plan PD). Soient respectivement PM3d et pMe ces fonctions membres. Ces fonctions d'appartenance fournissent des valeurs qui indiquent combien un élément vérifie la propriété associée à la fonction d'appartenance: ces valeurs sont comprises entre 0 (la propriété associée n'est pas vérifiée) et 1 (la propriété associée est complètement vérifiée). Afin de déterminer la fonction d'appartenance pMe, la transition 1/0 (objet/arrière-plan) présente dans les masques bidimensionnels Me extrapolés comme précédemment décrit est utilisée. Dans l'état actuel de ces masques bidimensionnels, un défaut de précision peut survenir sur les transitions objet/arrière-plan des masques bidimensionnels Me. Ainsi, il est sensé de pondérer ou de rendre flou la transition objet/arrière-plan. Pour cela, les masques bidimensionnels Me sont considérés comme des images dont les pixels prennent des valeurs entre 0 (arrière-plan) et 1 (objet considéré), sur lesquelles ont été appliqués une dilatation sur la base d'un disque puis un filtre passe-bas. La dilatation assure que le noyau de la fonction d'appartenance pMe (quand celle-ci est égale à 1) couvre l'ensemble de l'objet considéré, alors que l'application du filtre passe-bas permet d'obtenir des transitions douces objet/arrière-plan. Une taille du disque utilisé pour la dilatation doit correspondre à une taille d'un noyau utilisé durant l'utilisation du filtre passe-bas (par exemple un filtre moyenneur) si celui-ci est fini. Cette condition assure qu'un résultat de l'application du filtre passe-bas sera égal au résultat dans les zones du masque bidimensionnel Me � marquées � comme appartenant à l'objet considéré avant la dilatation. En variante, si le filtre passe- bas est un filtre à réponse impulsionnelle infinie (comme un filtre gaussien), la taille du disque de dilatation doit être d'une valeur à partir de laquelle les coefficients du noyau du filtre passe-bas peuvent être négligés (dans un filtre gaussien, le noyau peut être petit si l'écart type est petit et si les coefficients sont représentés seulement sur peu de bits informatiquement parlant). Il est à noter que le défaut de précision peut être modélisé par la taille du noyau du filtre passe-bas : des filtres avec de grand noyaux contenant des valeurs significatives gardent seulement les basses fréquences et ainsi étalent (floutent, pondèrent) plus les transitions objet/arrière-plan. Afin de déterminer la fonction d'appartenance }1M3d correspondant à la propriété associée au masque digital tridimensionnel M3dr les informations fournies par les propriétés associés aux masques bidimensionnels Me (pour chaque angle e associé à chaque masque bidimensionnel de la série de masques bidimensionnels extrapolés) doivent être rassemblées en utilisant un opérateur d'agrégation (cf. figure 5). Un tel opérateur d'agrégation peut être un opérateur t-norme T. Cela permet de calculer la fonction membre }1M3d comme une contrepartie floue de l'application d'un � ET � logique. De nombreux opérateurs t-normes existent : probabiliste, drastique, de Zadeh, de Lukasiexicz, etc. Le procédé a été mise en oeuvre en utilisant l'opérateur t-norme probabiliste. Ainsi, la fonction membre pM3d est déterminée de la manière suivante : Vv e volume d'intérêt (20), ,uM3d (v) =10' ,Ume (proj(9,v)) Ainsi, nous obtenons un masque digital tridimensionnel Mid final pondéré. Bien entendu de nombreuses modifications peuvent être 5 apportées à l'invention telle que précédemment décrite sans sortie du cadre de celle-ci. REVENDICATIONS 1. Procédé d'extrapolation d'un masque bidimensionnel MO2 dans un dispositif de radiographie du type comportant une source (S) de rayons X mobile prenant au moins deux positions Si eet S02 d'un espace associées à des masques bidimensionnels Met et MO2 respectivement, des moyens d'enregistrement (10) sensiblement plans disposés en regard de la source (S) et comportant une limite (d, d'), et un volume d'intérêt (20), comprenant un objet apte à être radiographié, situé entre la source (S) et les moyens d'enregistrement (10), le procédé comprenant des étapes de : a) estimation au moins partielle d'une série de masques bidimensionnels My associés à une série de positions Sy de la source (S) situées entre les positions Sel et S02 de la source (S), à partir des masques bidimensionnels Mol et Mot, b) évaluation au moins partielle d'un masque tridimensionnel intermédiaire (21) de l'objet à partir de la série de masques bidimensionnels My et des masques bidimensionnels Met et Me2r et c) extrapolation du masque bidimensionnel M92 au-delà de la limite d ou d' des moyens d'enregistrement (10) selon une position relative entre Sel et S02 à partir du masque tridimensionnel intermédiaire (21). 2. Procédé selon la revendication 1, caractérisé en ce 30 que, préalablement à l'étape a), si le masque bidimensionnel Me1 associé à la position de la source Sel n'est pas disponible, le masque bidimensionnel Met est alors extrapolé à partir de masques bidimensionnels Me disponibles. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que l'extrapolation de l'étape c) comporte une projection du masque tridimensionnel intermédiaire (21) à partir de la position S92 sur un plan passant par les moyens d'enregistrement (10). 4. Procédé selon l'une des revendications 1 à 3, caractérisé en ce que l'estimation au moins partielle de l'étape a) s'effectue sur la limite d ou d' des moyens d'enregistrement (10) selon la position relative entre Soi et So2. 5. Procédé selon la revendication 4, caractérisé en ce que l'estimation au moins partielle de l'étape a) comporte, pour chacun des masques bidimensionnels My, une étape de détermination d'un point Ty situé sur un bord de l'objet projeté sur la limite d ou d' des moyens d'enregistrement (10), la source (S) occupant la position S. 6. Procédé selon la revendication 5, caractérisé en ce que les points Ty sont estimés par une interpolation linéaires entre des points Tel et Tee situés sur un bord de l'objet projeté sur la limite d ou d' des moyens d'enregistrement (10), la source (S) occupant les positions Sel et S92 respectivement. 7. Procédé selon la revendication 6, caractérisé en ce que, si le masque bidimensionnel Me1 associé a la position Sel de la source (S) n'est pas disponible, le point Tel est alors extrapolé à partir de points 5 Te disponibles. 8. Procédé selon l'une des revendications 1 à 7, caractérisé en ce qu'il comporte une étape supplémentaire de : 10 d) application d'une fermeture sur le masque bidimensionnel Met extrapolé. 9. Procédé selon l'une des revendications précédentes, caractérisé en ce que, préalablement à l'étape a), le 15 procédé comporte une étape de détermination d'une limite d'extrapolation lfin sensiblement parallèle à la limite d ou d' des moyens d'enregistrement (10) et située en dehors des limites des moyens d'enregistrement (10). 20 10. Procédé selon la revendication précédente, caractérisé en ce que l'étape b) comporte, pour chaque ligne 1 sensiblement parallèle à la limite d'extrapolation lfin et située entre la limite d ou d' des moyens d'enregistrement (10) et la limite 25 d'extrapolation lfin, des sous étapes de : b1) estimation d'un plan P1e2 passant par la position Set et la ligne 1, b2) estimation du masque tridimensionnel intermédiaire (21) pour chaque voxel (v) situé 30 à l'intersection du plan P1e2 et du volume d'intérêt (20). 11. Procédé selon la revendication 9 ou 10, caractérisé en ce que l'étape c) comporte une étape de projection sur chaque ligne 1 sensiblement parallèle à la limite d'extrapolation lfin et située entre la limite d ou d' des moyens d'enregistrement (10) et la limite d'extrapolation lfin du masque tridimensionnel intermédiaire (21). 12. Procédé d'extrapolation d'une série de masques bidimensionnels Me dans un dispositif de radiographie, caractérisé en ce que chaque masque bidimensionnel Me de la série de masques bidimensionnel est extrapolé par un procédé selon l'une des revendications 1 à 11. 13. Procédé selon la revendication 12, caractérisé en ce que le procédé selon l'une des revendications 1 à 11 est appliqué de manière itérative. 14. Procédé selon la revendication 13, caractérisé en ce que, à chaque itération, les moyens d'enregistrement (10) sont étendus en leur adjoignant des moyens d'enregistrement virtuels correspondant à une partie commune à tous les masques bidimensionnels Me déjà extrapolés. 15. Procédé d'estimation d'un masque digital tridimensionnel à partir d'une série de masques bidimensionnels Me dans un dispositif de radiographie du type comportant une source (S) de rayons X mobile prenant au moins deux positions Sel et Set de l'espace associées à des masques bidimensionnels Mol et M92 respectivement, des moyens d'enregistrement (10) sensiblement plans disposés en regard de la source (S) et comportant une limite d ou d', et un volume d'intérêt (20), comprenant un objet apte à être radiographié, situé entre la source (S) et les moyens d'enregistrement (10), le procédé comprenant des étapes de . a) extrapolation de la série de masques bidimensionnels Me par un procédé selon l'une des revendications 1 à 14 ; et b) détermination du masque digital tridimensionnel à partir de la série de masques bidimensionnels Me extrapolés. 16.Procédé selon la revendication 15, caractérisé en ce que l'étape b) comporte des sous-étapes de : b1)pour chaque masque bidimensionnel Me, application d'une dilatation suivie d'un filtre passe-bas pour obtenir une fonction d'appartenance }1Me ; b2)évaluation d'une fonction d'appartenance PM3d a partir de la fonction d'appartenance pMe en utilisant un opérateur t-norme ; b3) détermination du masque digital tridimensionnel à partir de la fonction d'appartenance }1M3d• 17.Procédé selon la revendication 16, caractérisé en ce que l'opérateur t-norme est de type probabiliste. 18. Dispositif de radiographie du type comportant : 30 • une source (S) de rayons X mobile prenant au moins deux positions Sel et Set de l'espace,25• des moyens d'enregistrement (10) sensiblement plans disposés en regard de la source (S), • un volume d'intérêt (20), comprenant un objet apte à être radiographié, situé entre la source (S) et 5 les moyens d'enregistrement (10), caractérisé en ce qu'il comporte des moyens de mise en oeuvre d'un procédé selon l'une des revendications 1 à 17. 10 15