Brevet FR2598249A1 — publié le 19871106

Dispositif de reconstruction d'image par modulation d'un faisceau incident de radiations La présente invention concerne un dispositif de reconstruction d'image par modulation spatiale et/ou temporelle, d'un faisceau incident de radiations, ce faisceau étant composé d'électrons ou de photons visibles ou invisibles. Elle s'applique notamment à la radiologie et particulièrement à la réalisation de filtres spatiaux de rayons X, de conformateurs de radiodiagnostic et de systèmes de visualisation en général, particulièrement simples, dans le cas o le senseur d'image est un tomographe. Les modulateurs de faisceau permettant la visualisation d'une 10 image peuvent être classés en deux catégories, suivant qu'ils modulent ou non, directement chaque point ou pixel de celle-ci. Une première catégorie se rapporte au cas o la modulation porte sur un ensemble de pixels, chaque pixel n'étant pas directement accessible. Ce mcyen de visualisation est adapté au cas o l' Pinformation 15 permettant la modulation du faisceau n'est pas disponible sous forme de niveau de gris par pixel comme c'est le cas par exemple en tomographie appliqué à Pimagerie du corps humain, avant traitement informatique. La deuxième catégorie est relative au domaine de la visualisation point par point des images.Dans ce cas, la modulation de faisceau 20 s'effectue toujours dans un espace à deux dimensions, par des écrans à cristaux liquides, des tubes cathodiques etc... et, elle suppose nécessairement l'existence de modulateurs élémentaires, dont la taille correspond à la taille des pixels ou points de l'image, ainsi qu'à celle de moyens appropriés pour les adresser. Elle englobe tous les dispositifs actuellement 25 connus qui permettent la modulation locale d'un faisceau lumineux, en faisant varier son émission, sa transmission ou sa réflexion. Le besoin d'obtenir des images de bonne qualité augmente les difficultés de réalisation des modulateurs de faisceau, car plus la taille des modulateurs élémentaires répartis sur une surface donnée est réduite, 30 plus leur nombre augmente et plus le dispositif d'adressage et de commande de ceux-ci devient compliqué. On peut constater également que leur emploi est actuellement limité à la modulation de faisceaux de photons d'énergie lumineuse. Le but de ['invention est de pallier les inconvénients précités. A cet effet, rinvention a pour objet, un dispositif de reconstruc5 tion d'image par modulation d'un faisceau incident de radiations caracté' ' risé en ce qu'il comprend au moins un réseau de modulateurs, les modulateurs d'un réseau étant disposés suivant des lignes parallèles de manière a recouvrir la surface d'un plan à peu près perpendiculaire à l'axe du faisceau, chaque réseau pouvant être décalé angulairement autour de laxe du faisceau et des moyens de commande pour faire varier les coefficients de transmission et/ou de réflexion des modulateurs en fonction de limage à reconstruire. Le dispositif selon l'invention a pour principal avantage qu'il peut être mis en oeuvre dans tout type de modulation de faisceau grâce au fait que la modulation employée ne porte pas directement sur des pixels, mais sur des lignes de limage à reconstruire, linformation permettant cette modulation pouvant être la connaissance de l'image elle-même, ou une décomposition de celle-ci en ensemble de lignes tel qu'on peut l'obtenir par un procédé tomographique. Il présente également comme avantage dcêtre simple à réaliser, dès que dans les applications, et notamment pour celle concernant le filtrage spatial d'un faisceau de rayon X o chaque modulateur élémentaire doit être capable d'une absorption variable, le modulateur élémentaire de la taille d'un pixel présente une certaine complexité de mise en oeuvre, du fait de ses dimensions ou de sa disposition au voisinage des autres pixels qui le rendent difficilement accessible. D'autres caractéristiques et avantages de linvention apparaîtront ciaprès a l'aide de la description faite en regard des dessins annexés qui représentent: - la figure 1, un schéma de principe d'un mode particulier de réalisation dun dispositif selon linvention, - la figure 2, une illustration d'un mode de reconstruction d'une image de taille 3 x 3 à laide d'un réseau de trois modulateurs linéaires, - la figure 3, un schéma pour illustrer la fonction de détection 35 d'un modulateur élémentaire, la figure 4, un schéma pour illustrer le principe de reconstruction d'image selon rinvention, - les figures Sa et 5b, une application du dispositif selon l'invention à la réalisation d'un filtre centré, - la figure 6 un exemple d'image obtenue à l'aide dun dispositif selon rinvention dans lequel chaque modulateur est constitué par un conduit rempli d'un liquide de contraste à ses deux extrémités. Selon l'invention, chaque modulateur élémentaire, au lieu de correspondre à un pixel, a les dimensions dune ligne de l'image et le 10 modulateur complet se présente comme un réseau de modulateurs élémentaires disposés suivant des lignes parallèles entre elles, ces lignes n'étant pas nécessairement rectilignes. Pour assurer la reconstruction de l'image, les modulateurs élémentaires doivent pouvoir occuper un certain nombre d'orientations, ce qui peut se faire successivement dans le temps, ou 15 simultanément par plusieurs réseaux juxtaposés décalés angulairement entre eux, toutes les solutions intermédiaires étant évidemment possibles. Un mode de réalisation correspondant à un type de solution intermédiaire est représenté à la figure 1. Le dispositif selon l'invention qui est représenté à la figure 1 20 comprend un ensemble de plusieurs réseaux plans 1, 2 et 3 de modulateurs élémentaires 4. Les plans des réseaux sont orientés parallèlement entre eux et les modulateurs 4 sont disposés suivant des lignes parallèles dans chacun des réseaux. Ces derniers sont orientés de façon que les directions des modulateurs 4 d'un réseau soient décalées angulairement par rapport à 25 celles des autres réseaux. Les plans des réseaux 1, 2 et 3 sont disposés dans un faisceau de radiations 5 perpendiculairement à leur direction de propagation, ce faisceau de radiations pouvant éventuellement être émis par une source de rayons X, ou toute source susceptible d'émettre un rayonnement d'électrons ou de photons visibles ou invisibles. Les modulateurs élémentaires 4 sont constitués de façon connue, soit par des volets en matériau lourd comme le plomb, qui, en position ouverte, vue par la tranche, présentent une bonne transmission et en position fermée, en faible incidence, absorbent les photons, soit encore par des conduits de matière absorbante introduite aux extrémités des 35 modulateurs élémentaires. Suivant encore dautres modes de réalisation possibles de l'invention la commande des volets peut être assurée par des moyens piézoélectriques et dans ce cas les volets sont plaqués de façon connue sur une structure piézoélectrique de type par exemple bimorphe, travaillant avec -. 0 S une extrémité libre ou en flambement, ou encore, par des moyens:.hydropneumatiques, électrostatiques ou électromagnétiques, éventuellement constitués par un ou plusieurs moteurs déportés aux extrémités des -. volets. Sur la figure 1, les réseaux 1, 2, 3 sont maintenus solidairement 10 les uns aux autres par tous moyens de fixation connus, non représentés. Ils sont mis en rotation par pas angulaires successifs, autour crd'un axe de _- rotation (X', 0, X), parallèle à la direction de propagation du faisceau 5, -';:par des moyens de mise en rotation 6, de types connus, constitués éventuellement par un moteur pas à pas ou tout dispositif équivalent. Les moyens de mise en rotation 6 entraînent un plateau ou une couronne 7 transparents aux radiations sur lesquels reposent ou à l'intérieur desquels sont fixés les réseaux plans 1, 2 et 3 de modulateurs 4. Ces modulateurs 4 sont commandés par des moyens de commande 8 formés, de façon connue, par un organe de traitement de l'information, du type ordinateur ou tout dispositif équivalent, comprenant au moins une mémoire 9 pour stocker les valeurs des coefficients de transmission de chaque ligne de modulateur 4 nécessaires à la reconstruction de rimage. Les moyens de commande 8 assurent également la commande des pas de rotation des moyens de mise en rotation 6 en synchronisme avec la commande des lignes de modulateur 4. Un senseur dimage 10, formé par exemple par une caméra, est. 'S S placé suivant Paxe (X',0,X) du faisceau 5. Il permet de visualiser limage reconstruite résultant de l'excitation par les moyens de commande 8 des différents modulateurs 4. De préférence le senseur d'image sera pourvu de moyens d'intégration pour permettre la juxtaposition des différentes images retransmises à chaque pas de rotation du plateau ou de la couronne 7 et à chaque commande des modulateurs 4. Le fonctionnement du dispositif selon rinvention est explicité ciaprès, en prenant à titre d'exemple, le cas d'un unique réseau de modula35 teurs élémentaires linéaires, parallèles entre eux, animé d'un mouvement de rotation pas à pas autour de l'axe de rotation (X',0,X), le réseau étant constitué par K lignes de modulateurs élémentaires parallèles, pouvant 2irn présenter par rapport à une direction de référence des angles 2-n, avec 0 désigne par a (k,n) l'intensité de la ligne k, lorsque l'angle est égal à -2 et par I (X,Y) Pintensité du point de coordonnées X,Y reconstruit sur le senseur d'image 10, et le mot "intensité", désigne ici une grandeur physique liée à l'information en termes de niveau de gris sur l'image, que le senseur d'images 10 est capable d'intégrer. Dans cette hypothèse, la commande du modulateur peut être réalisée, du moins en théorie, de trois façons différentes. Selon une première approche, ou approche algébrique, on peut considérer que le dispositif précédemment décrit forme un système discret, pour lequel la relation I (X,Y) = N1 a[k(n,X,Y),n] (1) est toujours n=0 vérifiée. Dans la relation (1), k est une fonction de l'angle (donc de n) et du point (XY) car, la ligne qui contribue à l'intensité d'un point de coordonnées (X,Y) donné est différente suivant l'angle de rotation du réseau. I (X,Y) a dans ces conditions les dimensions dune matrice M x M, la relation (1) étant représentative d'un système de M équations à k x N 20 inconnues. Elle permet de reconstruire l'image donnée I (X,Y) par des valeurs de modulations a (k,n) qui sont les solutions du système à M2 équations. Ces solutions ne peuvent exister que si K x N> M2, ce qui donne une condition nécessaire sur le nombre de positions que le système doit 25 occuper lorsqu'il dispose dun nombre de lignes donné, pour assurer la reconstruction. La recherche de la solution, si elle existe, passe par l'inversion, si elle est possible, d'une matrice de taille M2, binaire et en général lacunaire. Cette approche est un bon moyen pour tester la méthode de reconstruction, cependant elle met en évidence des défauts notamment sur la forme et le positionnement des pixels reconstruits, du type de ceux qui apparaissent sur l'image reconstruite de la figure 2 qui résulte de l'application de la méthode précédente au cas o M = 3, k = 3 et N = 3. En effet, selon l'image reconstruite de la figure 2, la matrice[I] de pixels A, B, C, D, E, F, G, H est obtenue à partir d'un réseau de k = 3 6 k modulateurs linéaires prenant N = 3 orientations successives: i = (i-l) avec 1I obtenue pour chaque position du réseau est momentanément mémorisée par le senseur dimage 10 d'un temps suffisant pour pouvoir être juxtaposée aux images suivantes fournies au cours des autres positions successives du réseau, la ligne l(resp 2,3) du réseau étant caractérisée par une intensité ai(resp Ai', i) lorsque son inclinaison vaut (i-l)' avec 1\ La répartition des intensités sur chaque modulateur du réseau 10 pour les trois orientations Q. est alors définie par: 1' 1 2 OEZ" 2 ' 1 3 â3 i r,,j" [11]={:i'$i' J;['2j = Y2,i2B2 2 e ['3] 13 P3 PB3 " I" I j L ' il, lie e 21 lsz3 2 '(2 l 3 3 l 3 d'o il résulte Pimage I qui est formée par la juxtaposition de [i, l[12 et 15 [I31; cette juxtaposition vérifiant naturellement la relation matricielle A 111000000 1 B 1 10001000 a2 C 100010001 a3 D 0 1 1 1 00000 EE = 0 0 1 1 10 0 0 x 2 F 000100011 3 G 0 0 1 0 1 0 1 0 0 H 00000 1 1 10 l 000000111 l3 C'est-à-dire les relations A = al +a2 + a3 C = a1 + L2 + 3 C = ai + "2 + '3 1 +2 + 3 qui traduisent le recouvrement des différentes intensités d'éclairement ai ' /i et Xi des modulateurs combinés trois à trois. Cependant comme le montre la figure 2, la surface commune de recouvrement varie d'un pixel à rautre, de sorte que l'image reconstruite obtenue par cette méthode n'a pas une qualité irréprochable. Une deuxième approche ou approche analytique consiste à supposer que l'on dispose d'une méthode permettant, à partir d'une "image" caractérisée par une fonction m (x,y), de calculer les intensités a(k,n) - le problème consiste à déterminer à partir d'une image donnée 10 i(x,y) à reconstruire, de quelle pseudo-image m(x,y) il faut partir pour qu'un réseau de modulateurs linéaires, de valeurs d'intensités c(k,n) calculées par une méthode donnée, reconstruise l'image i(x,y). Cette formulation du problème est en apparence plus compliquée, mais elle a l'avantage de toujours présenter une solution m (x,y). Soit, en effet la fonction de détection an,X,Y d'un modulateur linéaire 4 d'orientation 2n- passant par un point (X,Y) tel que représenté à la figure 3, (an, X,y(X,y) = 1 lorsque le modulateur passe par le point (X,Y) et a n,X,y(X,y) = 0 dans le cas contraire). soit a k(n, X, Y), ni] f an Xy(x,y).m(x,y)dx.dy (2) la relation entre l'intensité du modulateur linéaire k (passant par X, Y), d'orientation 2-n et la fonction m(x,y). N On a alors, si r(x,y) désigne l'intensité de l'image reconstruite. +oo N-1 r(x,y) =? W ae an,Xy(x,y).m(x,y).dx.dy (3) Comme on cherche à obtenir une image donnée i(X, Y), la 25 fonction m(x, y) cherchée est solution de: r(X,Y) = i(X,Y) (4) Un modulateur élémentaire, d'inclinaison 2--n donnée, passe par N de nombreux points X, Y. La surface Sk,n qu'il occupe peut alors être - représentée par la relation Sk,n= [(XY)E Ox,Oy tels que k(n,X,Y) = Cti. Si XknYk,n est un point particulier de Skn pour n fixé, par exemple, son centre, on a V(X,Y)(Sk,n, an,X,Y= anXknYk,n Soit Pn la fonction qui, à tout point (X,Y) du plan, fait correspondre le point (Xkn, Y) à n fixé; pn est une projection sur une famille de kn k,n points, que Pon peut définir comme suit: Yn([4,N-lJ et V(X,Y)(e Support de limage, il existe un unique modulateur k tel que (X,Y) ESk n Alors Pn(XY) = (Xkn, Yk,n), centre de Skn On a: Vn, V (X,Y)E support de l'image, anXY = anp (XY) fd'autre part, la fonction de détection d'un modulateur élémentaire 10 centré en Pn(X,Y) se déduit par translation de celle d'un modulateur centré en O, et de même orientation: l'équation (4) s'écrit alors: n(xy =0 >: wn[Pn(XlY)-(xy0l m(xy).dx.dy i(X,Y)_= j an[:OXY) _Y) avec an(xy) anoo [(xy N-I -15 Soit i = 2 (aen x m)o Pn n=0 ou a désigne Popérateur de convolution et o désigne ropérateur de composition. En considérant la surface Skn d'un modulateur linéaire comme un rectangle de largeur d, et de longueur D, un tel modulateur, centré à 2 l'origine d'un repère (O,x,y) et aligné sur Oy aura comme fonction de détection ao(x,y) = Rect k(x).Rect2(y) Si D est très grand devant les dimensions de l'image efficace (ce qui revient à prolonger i(X,Y) et m(x,y) par la valeur O au delà du support efficace), et si d est très petit (ce qui suppose un réseau avec beaucoup de 25 modulateurs linéaires), on ne fait pas beaucoup d'erreur en considérant que chaque modulateur linéaire passant par (X,Y) est centré en (X,Y). On a alors Pn = Identité pour tout n, et léquation (5) devient: N-I i= (: an)Im (6) n=0 n L'image reconstruite s'écrit alors comme la convolution d'une 30 image source à définir par un filtre pondérateur, illustré aux figures Sa, b. L'équation (6) a pour Transformée de Fourier: N-I 1= (5 A'n).M (7) n=O n ao0(X,Y) a pour Transformée de Fourier Ao(fx,fy) = d.D.sinc (d.fx).sinc (D.fy) et an(X,Y) se déduit de ao(X,Y) par une rotation d'angle 2N sinc étant le symbole de la fonction sinus cardinal définie comme suit sinc(x) = x V x R" X sinc(0) = 1 Donc an(X,Y) a pour Transformée de Fourier 21rn An(fx,fy) = d.D.sinc(d(fx. cos _ sin 2n)sinc (D(f.sin 2n cos 2n n x y x N Y' N x N Y' N d'o on déduit: N-I dfr I(fxfy)=M(xfy) x' d.D.sinc d(fx.COS 2 sin _ sinc (f.sin n n=O xc N Y N IDfx N + fy.cos 2n)] puisque A'n(f fy) = An(-f -f). nx, y n x, y ceci permet, en prenant les précautions nécessaires de remonter à la fonction cherchée: m(*,y)'nf 2Co n - 21rn M(XY) TF (f jf(2 d.D.sinc.d(f.o N f 'sin 21rn) sinc[D(fx. sin -2 + fy.cos N)])} (8) Les valeurs des modulations des différents modulateurs peuvent alors être déduites en appliquant la relation (2). Cependant la méthode précédente donne une solution approchée difficile à calculer et on pourra lui préférer une troisième méthode dictée par des considérations géométriques qui permettront d'approcher plus 25 grossièrement le résultat mais avec ravantage qu'elle envisage le problème d'une façon plus simple. On peut en effet chercher à reconstruire l'image I sous la forme d'une matrice M x M en prenant en considération que le nombre N d'orientations minimum est tel qu'a tout point de 1, il existe une direction permettant d'atteindre n'importe quel autre point de l'image. La projection p. de 'image I 'elon la direction n, est alors définie par: P = m.j dans la direction n passant par le point (iojo) (io, o) Pn peut être vue comme l'image du sytème de modulation figé dans la position angulaire -n. Comme chaque position du système est pondérée par le poids I, et que le senseur d'image 10 peut être considéré comme un dispositif intégrateur, chaque coefficient de la matrice Pn vue comme image du modulateur vaut dans ces conditions: 10 Pn N m. On aalors: N m.. N10 p. p 1 +iN-' c (9) n. pn N-- n=l o I est la matrice image, et C la matrice unité. Cette dernière methode apparaît dans ces conditions comme un 15 moyen extrêment simple de reconstruction puisque pour chaque valeur n, les intensités at(k,n) sont les k coefficients différents de Pn' Elle est entachée cependant d'un effet parasite que l'on peut appeler "fond continu. Cependant on peut pallier cet inconvénient en modifiant les 20 projections Pn' Cette modification peut s'effectuer: -- en multipliant les projections Pn par un coefficient "amplificateur de dynamique". - en supprimant la composante continue de chacune des projections pn (pour se ramener à une situation: telle que Y n iflf pn =0), ce 25 qui diminue la composante continue de la reconstruction exprenée par léquation (9). - en diminuant davantage le fond continu par l'incorporation de connaissances "à priori" sur l'image, par exemple, les points périphériques, lorsqu'ils sont en dehors de la zone à observer, pourront être considérés 30 comme non pertinents, et on pourra de cette façon améliorer la qualité de la zone centrale de l'image. Il La méthode qui vient dêtre exposée est illustrée ci-après à l'aide de l'exemple suivant. Dans cet exemple de reconstruction, l'image à reconstruire est définie par une matrice de coefficient mij ayant les valeurs numériques définies par le tableau: r0101 M= 1 32 2 2 1 et cette matrice est projetée suivant un nombre N = 8 projections. Les matrices de projection obtenues pour des valeurs de 2wn N gaes,ww w w 2w 3w. 5w '3 '2' 3' 46,-sont respectivement: 6 '4 '- '-i ' î I 4 il 2wn pour 2n = 0, pour N ' pour 21mn = pour N = 1 N = 4 ' 2wn pour 15 pour - N -3' 2wn 2w pour =-f, N - 23' 2wn 3w pour 2 N- 4 ' et M 1661 M= 5 0 1 1 M 2 = 133 4 4 2 1 2o 2 S M 3 = 664 3 2 M 4 = -3 2 1 3 6 3 M 6 =,33 2 2d M 7 = 4 2 3 4 212n 5w pour M = N 6' 8 2 2 2 En additionnant chaque matrice élémentaire il vient Mm c 7 21 14 + 12 Lîî i=l I 14 14 1 1 Pour simplifier il est possible d'enlever un fond continu de la matrice M1 correspondant à une fois la matrice unité et un fond continu de la matrice M5 correspondant à 3 fois la matrice unité. En considérant que les points (1,1) et (1,3) de l'image, de valeur nulle, ne sont pas significatifs et correspondent à une zone o il n'y a pas d'objet à observer on peut encore soustraire respectivement les matrices 1 01 22 222II et des matrices M2, M3, M4, M6, M7 et M8 On, obtient alors une représentation matricielle simplifiée telle que: 721 14 + 100 o 1414 0I 0 Soit après pondération " l avec 0 C= 0 0 ' 0 0 (La matrice E contenant les points non pertinents qui représentent Pécart par rapport à la solution). Dans le cas o la matrice image ne contient pas de O, la matrice est identiquement nulle. Le fond continu, lorsque l'on a N projections vaut alors en chaqJe point: 259O249 C = N 1 Min(mij) ij D'autres considérations interviennent, pour clore cette approche du problème: - Le "fond continu" en lui-même peut encore être compensé par 5 un effet de seuil du senseur d'image intégrateur. Dans le cas d'une image "en négatif", on pourra le compenser en augmentant la puissance du faisceau incident. - Enfin, le problème du filtrage n'apparaît pas ici du fait des hypothèses de départ simplificatrices. Cependant, dans le cas d'un filtre 10 spatial de rayons X dont on ne veut pas qu'il ait une structure trop bien définie, le flou résultant de cette convolution peut être un atout. (Chaque pixel a une absorption qui dépend de celle des pixels d'un certain voisinage). Dans le cas particulier o l'invention est appliquée à un conforma15 teur pour le radiodiagnostic, o l'image à reconstruire est à 2 niveaux de gris (Noir et Blanc), et représente une forme assez simple, (exemple le coeur) on pourra opter pour un dispositif statique formé par plusieurs réseaux juxtaposés plut8t que de faire tourner un seul réseau autour de raxe du faisceau lorsque la forme à représenter est convexe (majorité des 20 cas), le nombre de réseaux de modulateurs dépendra alors de la précision souhaitée sur sa définition, mais pourra cependant être restreint. La présente invention n'est pas limitée aux exemples de réalisation qui viennent d'être décrits et il va de soi que ceux-ci pourront subir de nombreuses modifications et adaptations en fonction des applications. 25 Par exemple, pour toutes les formes ayant la propriété suivante: "en chaque point du contour de la forme, il passe au moins une droite dont Pintersection avec la forme est connexe", comme celle représentée sur la figure 6, une solution avantageuse pourra consister à utiliser un "pseudomodulateur" linéaire, constitué crd'un réservoir linéaire autorisant l'admis30 sion de liquide de contraste par ses 2 extrémités. Chaque modulateur pourra alors être séparé en 3 zones de dimensions variables, dont les 2 extrêmes seront opaques, et la zone centrale transparente. Le nombre de réseaux de tels modulateurs dépend de la forme en question, mais peut être estimé à la moitié de ceux requis dans le cas précédent. Suivant d'autres modes de réalisations, on pourra également éviter Pemploi de la modulation temporelle précédemment décrite, en disposant simultanément des'réseaux de modulateurs linéaires en série, décalés en rotation les uns par rapport aux autres. On n'aura plus besoin dans ce cas de senseur d'image intégrateur. De mêmes toutes les solutions intermédiaires seront également possibles: en disposant par exemple N réseaux de modulateurs légèrement décalés les uns par rapport aux autres, l'ensemble étant animé d'un mouvement de rotation, la vitesse de rotation et la fréquence de commutation seront alors divisées par N. 10 Les approches permettant de concevoir des algorithmes de reconstruction de limage qui génèrent les commandes des modulateurs linéaires en fonction de leur angle ont été données à titre indicatif. Rien n'empêche en particulier d'imaginer une suite (p'l) telle que VN' la n relation, N" Norme. ( 1 P' - I) soit minimale n=l ou Norme est une norme de l'espace des matrices à définir et I est la matrice image que lon souhaite obtenir en fin de reconstruction. Ceci permet d'envisager un système de reconstruction dont la durée de 20 reconstruction deviendrait un paramètre indicateur de qualité, donc que l'on n'aurait plus besoin de connaître à priori. Dans le cas o les modulateurs utilisés ne fonctionnent qu'en tout ou rien, plusieurs solutions sont encore possibles: Si le temps de commutation est suffisamment bref, pour une 25 position donnée, on peut faire varier la durée pendant laquelle le modulateur est a "1" (modulation temporelle pendant la durée d'exposition à une position angulaire donnée). À Si lon sait faire des modulateurs linéaires dont la position "1" (absorption) peut varier en absorption d'un type de modulateur à l'autre, 30 un empilement de quatre réseaux dont les modulations nominales sont bien choisies suffit à reconstituer une image à 16 niveaux de gris. On peut aussi composer rimage en 4 fois, ce qui multiplie fréquence et vitesse de rotation par 4, en changeant l'énergie incidente à chaque flash, en agissant directement sur le tube, ou par l'intermédiaire d'un modulateur, ou en modifiant la durée, ou en combinant ces effets, ce qui donne aussi une image à 16 niveaux de gris. L'invention en question permet de reconstruire une image numérique présente dans la mémoire d'un ordinateur, mais aussi une image dont 5 on ne connait que la décomposition en réseaux de lignes, telle qu'un tomographe peut la fournir. Il est possible d'envisager le cas o le couplage tomographe/dispositif de reconstruction serait purement analogique. La reconstruction envisagée est principalement une modulation 10 d'absorption, mais rien ne s'oppose à ce que ron puisse prévoir le cas o l'on disposerait de modulateurs actifs: émetteurs de lumière lampes à décharge...) réflecteurs, modulateurs de phase (biréfringence...). Egalement on a pu remarquer que les modulateurs décrits ont la forme de lignes droites, mais rien ne s'oppose à ce que le même principe de l'invention s'applique aussi à des modulateurs de formes diverses, dans la mesure o ils recouvrent la surface de Pimage, et o il existe une famille de positions pour cette surface réduisant les intersections des modulateurs'à la taille d'un pixel. Dans le même esprit ces modulateurs linéaires décrits ont une 20 forme constante tout au long du mouvement, on-concevra que cette forme n'est pas limitative et Pon pourra envisager les cas o leur forme varie d'une position à l'autre. La taille des pixels sera alors donnée par l'intersection des formes pour toutes les positions. Par exemple, en optique, un modulateur linéaire peut être, dans un matériau homogène de 25 type magnétooptique, un domaine o l'aimantation est uniforme. Ces domaines peuvent voir leurs largeurs varier au cours du mouvement. Encore une autre manière de procéder sera dengendrer un réseau dans un matériau acoustooptique, à l'aide d'un transducteur piezoélectrique large bande. Ce réseau sera alors de pas variable. REVENDICATIONS 1. Dispositif de reconstruction d'image par modulation dun faisceau incident de radiations (5), caractérisé en ce qu'il comprend au moins un réseau (1, 2, 3) de modulateurs (4), les modulateurs d'un réseau étant disposés suivant des lignes parallèles, de manière à recouvrir la surface ( d'un plan à peu près perpendiculaire à laxe du faisceau (5), chaque réseau pouvant &tre décalé angulairement autour de raxe du faisceau, et des moyens de commande (8) pour faire varier les coefficients de transmission et/ou de réflexion des modulateurs en fonction de limage à reconstruire. 2. Dispositif selon la revendication 1, caractérisé en ce que 1 occupation du nombre déterminé d'orientation d'un réseau s'effectue successivement dans le temps par rotations (6, 7) d'angles élémentaires successifs du réseau autour de laxe longitudinal du faisceau. 3. Dispositif selon la revendication 1, caractérisé en ce qu'il comprend plusieurs réseaux plans (1, 2, 3) de modulateurs (4) maintenus 15 juxtaposés les uns aux autres, la direction des modulateurs de chaque réseau étant décalée relativement à la direction des modulateurs des autres réseaux. 4. Dispositif selon l'une quelconque des revendications I à 3, caractérisé en ce que les modulateurs (1) sont disposés suivant des lignes 20 droites.. 5. Dispositif selon l'une quelconque des revendications I à 4, caractérisé en ce que les modulateurs (4) sont constitués par des volets inclinables, opaques aux radiations. 6. Dispositif selon l'une quelconque des revendications I à 4, 25 caractérisé en ce que les modulateurs (4) sont constitués par des volets inclinables réfléchissant une source lumineuse. 7. Dispositif selon l'une quelconque des revendications 1 à 4, caractérisé en ce que les modulateurs (4) sont constitués par des tubes remplis de matière absorbante aux radiations. - 30 -8. Dispositif selon l'une quelconque des revendications I à 6, caractérisé en ce que les modulateurs (4) comprennent des guides de lumière émettant des rayons lumineux en réponse à une excitation des moyens de commande. 9. Dispositif selon la revendication 8, caractérisé en ce que les 35 guides de lumière sont constitués par des lampes à décharges. 10. Dispositif selon l'une quelconque des revendications 1 à 9, caractérisé en ce que les moyens de commande (8) sont constitués par un ordinateur comprenant au moins une mémoire (9), la mémoire (9) étant chargée à partir d'informations fournies par un tomographe. 11. Dispositif selon l'une quelconque des revendications précédentes, caractérisé en ce que le faisceau de radiation est un faisceau de rayons X. 12. Dispositif selon l'une quelconque des revendications I à 11, caractérisé en ce qu'il comprend un senseur d'image (10) placé en regard 10 des modulateurs pour permettre la visualisation de l'image reconstruite.