Brevet FR2947658A1 — publié le 20110107

PROCEDE ET DISPOSITIF D'IMAGERIE A RAYONS X A TRAITEMENT TRIDIMENSIONNEL L'invention se rapporte à un procédé d'imagerie à rayons X, du type scanographie, dans lequel, un corps à examiner étant reçu par un support, - on émet, à partir du foyer d'une source, un faisceau de rayons X i et j en direction du corps à examiner, - on détecte, à l'aide d'un détecteur irradié par le faisceau, une intensité atténuée en fonction d'une traversée des rayons X à travers le corps à examiner, sous différentes paires de positions angulaires 0a, 0b du corps à examiner ou de l'ensemble source ù détecteur autour d'un axe de rotation, les deux positions angulaires 0a, 0b des différentes paires présentant un écart angulaire 4 égal de préférence à 90 degrés, on convertit, à l'aide d'un convertisseur analogiqueùnumérique, les intensités détectées en données c(a) , p(b) et, - on traite, à l'aide d'un ordinateur programmé, les données c(a) , p(b) pour aboutir à une image exprimant un coefficient d'atténuation des rayons X en un ensemble de points {P0} du corps à examiner, en calculant, pour une paire de positions angulaires 0a, 0b, une valeur C du coefficient d'atténuation en un point PO de l'ensemble {P0} par application de la formule suivante : p(b) C(a) n n C n + n 2 n2 (k=1 ca> + 11 Pubs ) dans laquelle : C(a) est la donnée provenant de l'intensité détectée au point d'impact T(a) d'un premier rayon i passant par le point considéré PO dans la première position angulaire 0a, p(b) est la donnée provenant de l'intensité détectée au point d'impact R(b) d'un deuxième rayon j passant par le point considéré PO dans la deuxième position angulaire 0b, n k c a) est la charge associée au premier point d'impact T(a), définie par k=1 la somme des données ce provenant des intensités détectées dans la première position angulaire 6a, le long de (n) points T a) d'une première ligne de calcul de charge qui passe au plus près du premier point d'impact T(a) , et Un procédé de ce type est connu de la demande de brevet française FR 10 2 888 374 et de la demande internationale WO 2007/006560 déposées par le demandeur. Le coefficient d'atténuation C au point considéré PO du corps à examiner est calculé par la méthode dite de la génération directe. Le calcul ne nécessite que la connaissance des données c(a) et p(b) provenant des intensités détectées aux deux points d'impact des deux rayons X passant par le point 15 considéré PO sous les première et deuxième positions angulaires du corps à examiner ou de l'ensemble source ù détecteur autour de l'axe de rotation a La méthode de la génération directe est utilisée pour calculer le coefficient d'atténuation d'un ensemble {P0} de points appartenant à un plan de coupe du corps à examiner, perpendiculairement à l'axe de rotation S2. Dans ces conditions, 20 les charges associées aux premier et deuxième points d'impact sont calculées le long de deux lignes de calcul de charge définies par la trace sur le détecteur du plan de coupe, sous les première et deuxième positions angulaires du corps à examiner ou de l'ensemble source ù détecteur autour de l'axe de rotation û. La valeur des charges associées est constante pour tous les points PO appartenant au plan de 25 coupe. Les opérations élémentaires de calcul du coefficient d'atténuation C conduisent à un gain de temps très important. Le terme soustractif qui rend compte des charges associées aux premier et deuxième points d'impact est calculé une fois pour toutes et sa valeur est simplement rappelée dans le calcul du coefficient 30 d'atténuation C de tous les points PO du plan de coupe. De surcroît, le calcul du n E pub) est la charge associée au deuxième point d'impact R(b) , définie u=1 par la somme des données p(t,b) provenant des intensités détectées dans la deuxième position angulaire Ob, le long de (n) points Re d'une deuxième ligne de calcul de charge qui passe au plus près du deuxième point d'impact R(b) coefficient d'atténuation C ne nécessite pas d'étape de rétro- projection de l'image obtenue sous la première position angulaire dans l'image obtenue sous la deuxième position angulaire. Le calcul direct à partir des données c(a) et p(b) provenant des intensités détectées aux deux points d'impact conduit à une réduction importante des erreurs de calcul. Cependant, le calcul du coefficient d'atténuation C dans un plan de coupe perpendiculaire à l'axe de rotation û impose de déplacer, parallèlement à cet axe, le corps à examiner ou l'ensemble source û détecteur pour obtenir des images de toute une série de plans de coupe, autrement dit pour explorer un volume du corps à examiner. L'un de buts de l'invention est de modifier le procédé connu qui vient d'être décrit pour s'affranchir du calcul du coefficient d'atténuation dans un plan de coupe perpendiculaire à l'axe de rotation tout en conservant la rapidité de temps de calcul et la réduction des erreurs de calcul apportées par la méthode de la génération directe. A cet effet, l'invention a pour objet un procédé conforme à celui rappelé en introduction, caractérisé en ce que l'ordinateur est programmé pour effectuer les étapes suivantes : (1) pré- calculer, pour tout point générique d'un volume irradié par le faisceau émis par la source, la position sur le détecteur d'un premier point d'impact T d'un rayon i et d'un deuxième point d'impact R d'un rayon j, les deux rayons i et j passant par le point générique et se déduisant l'un de l'autre par une rotation du point générique ou du de l'ensemble source û détecteur autour de l'axe de rotation S2 dont l'angle est égal à l'écart angulaire A entre les deux positions angulaires (ea,eb) des différentes paires de positions angulaires du corps à examiner ou de l'ensemble source û détecteur autour de l'axe de rotation S2, (2) pré- calculer sur le détecteur la position de premières et de deuxièmes lignes de calcul de charge de points de détection Tpk;k=1û*n;p=1-*m et Rqu;u=1ù> n;q=1û m se prolongeant pour converger respectivement en un premier et un deuxième points géométriques dont les positions sont déterminées par la distance du foyer à l'axe de rotation et la distance du détecteur à l'axe de rotation, chaque ligne étant inclinée par rapport à l'axe de rotation d'un angle de pente yp, yq qui s'échelonne en une série de valeurs discrètes déterminées par une résolution spatiale du détecteur parallèlement à l'axe de rotation, (3) effectuer le calcul de la valeur du coefficient d'atténuation en considérant, en fonction de la position du point PO dans le volume irradié par le faisceau, les positions pré- calculées des premier et deuxième points d'impact T(a) ; R(b) et les première p et deuxième q lignes de calcul de charge pré- construites Tpa;k =1 n ; Re ;U =1-* n qui passent au plus près des premier et deuxième points d'impact (T(a) ; R(b) ) (4) recommencer l'étape (3) pour les différentes paires de positions angulaires (ea,eb) du corps à examiner ou de l'ensemble source - détecteur autour de l'axe de rotation et calculer une valeur moyenne du coefficient d'atténuation C au point considéré PO pour les différentes paires de positions angulaires, et (5) recommencer l'étape (4) et afficher à l'écran de l'ordinateur l'image exprimant les valeurs moyennes de l'ensemble {P0} des points P0. L'étape (1) de programmation de l'ordinateur permet de construire une base géométrique faisant correspondre à tout point du volume irradié par le faisceau les deux points d'impact des deux rayons qui passent par ce point. Les m premières et deuxièmes lignes de calcul de charge définies à l'étape (2) marquent la trace sur le détecteur de m plans inclinés par rapport à l'axe de rotation. L'écart angulaire entre deux plans est fixé en fonction de la résolution spatiale du détecteur parallèlement à l'axe de rotation. Il est choisi de préférence le plus faible possible pour former des lignes de calcul de charge les plus serrées possible. La base de correspondance géométrique pré- calculée à l'étape (1) et les m lignes de calcul de charge pré- construites à l'étape (2) sont utilisées à l'étape (3) de programmation de l'ordinateur pour calculer le coefficient d'atténuation C en n'importe quel point PO du corps à examiner. En fonction de la position de ce point PO dans le volume irradié par le faisceau, l'ordinateur recherche dans la base géométrique les premier et deuxième points d'impact correspondants, appelle la donnée c(a) provenant de l'intensité détectée au premier point d'impact sous la première position angulaire, appelle la donnée p(b) provenant de l'intensité détectée au deuxième point d'impact sous la deuxième position angulaire, recherche la ligne de calcul de charge qui passe le plus près du premier point d'impact et calcule la charge associé à ce point d'impact, recherche la ligne de calcul de charge qui passe le plus près du deuxième point d'impact et calcule la charge associée à ce point, et finalement calcule la valeur du coefficient d'atténuation C au point considéré sous la paire de positions angulaires considérées. Les étapes (4) et (5) permettent de passer d'une paire de positions angulaires à une autre pour calculer une valeur moyenne du coefficient d'atténuation au point considéré, puis de passer à un autre point du corps à examiner. Ainsi, le procédé selon l'invention permet de suivre n'importe quel point du corps à examiner se trouvant dans le volume d'irradiation ou d'illumination du faisceau émis par la source. La programmation de l'ordinateur conduit à affranchir le calcul du coefficient d'atténuation de la condition d'appartenance des points à un même plan de coupe perpendiculaire à l'axe de rotation. La pré- construction des lignes de calcul de charge et de la base de correspondance géométrique permet de conserver la rapidité et la réduction des erreurs de calcul de la méthode de la génération directe. Le temps de calcul est limité à ce qui est strictement nécessaire au nombre de points PO dont les valeurs du coefficient d'atténuation sont recherchées sans qu'il soit nécessaire de les extraire d'un ensemble de calculs de points beaucoup plus large. De surcroît, le volume à irradier est limité à ce qui est strictement nécessaire pour obtenir l'information pertinente à l'exécution du calcul. Par comparaison avec le procédé tomodensitométrique décrit dans les demandes antérieures citées en introduction, l'invention apporte un procédé spatiodensitométrique. D'autres avantages de l'invention apparaîtront à la lumière de la description des modes de réalisation illustrés par les dessins. La figure 1 est une vue dans le plan O,x,y de la géométrie d'un dispositif d'imagerie conçu pour la mise en oeuvre d'un procédé selon l'invention, le dispositif comportant deux ensembles source ù détecteur. La figure 2 est une vue en élévation, parallèlement à l'axe de rotation S2, d'un premier rayon i et de son point d'impact sur l'un des détecteurs du dispositif illustré par la figure 1. La figure 3 est une vue en élévation, parallèlement à l'axe de rotation û, d'un deuxième rayon j et de son point d'impact sur l'autre détecteur du dispositif illustré par la figure 1. La figure 4 est un tableau indiquant les coordonnées de deux points d'impact sous une succession de positions angulaires. La figure 5 est une vue de face de l'un des détecteurs, sur lequel sont reportés le point d'impact du premier rayon i et une ligne de calcul de charge passant au plus près de ce point d'impact. La figure 6 est une vue de face de l'autre détecteur, sur lequel sont reportés le point d'impact du deuxième rayon j et une ligne de calcul de charge passant au plus près de ce point d'impact. La figure 7 est un tableau indiquant des largeurs de marche, des nombres de points de détection par marche, un nombre de marches par ligne de calcul de charge et une hauteur atteinte par deux lignes de calcul de charge. La figure 8 est un tableau indiquant les coordonnées de deux points d'impact, une altitude Z* de deux points de référence correspondant et l'indice p ou q de deux lignes de calcul de charge passant au plus près des deux points de référence. La figure 9 est une vue dans le plan O,x,y de la géométrie d'un dispositif d'imagerie conçu pour la mise en oeuvre d'un procédé selon l'invention, le dispositif ne comportant qu'un ensemble source ù détecteur. La figure 10 est une vue dans le plan O,x,y de la géométrie du dispositif illustré par la figure 9, un point générique G ayant effectué une rotation A autour de l'axe de rotation S2. La figure 11 est une vue de face du détecteur, sur lequel sont reportés les points d'impact des premier et deuxième rayons i et j ainsi que deux lignes de calcul de charge passant au plus près de ces points d'impact. La figure 12 est une vue schématique en coupe d'un premier dispositif conçu pour mettre en oeuvre un procédé selon l'invention. La figure 13 est une vue schématique en coupe d'un deuxième dispositif conçu pour mettre en oeuvre un procédé selon l'invention. Le procédé selon l'invention est décrit tout d'abord par référence à un dispositif d'imagerie à rayons X comprenant deux ensembles source ù détecteur. Cependant, l'invention n'est pas limitée à une telle géométrie et l'exposé sera poursuivi par référence à un dispositif d'imagerie à rayons X ne comprenant qu'un ensemble source ù détecteur. Un dispositif d'imagerie à rayons X, du type scanographie, comprend, figure 1, deux sources 5 émettant un faisceau de rayons X à partir de deux foyers FI et F2. Les deux foyers F1 et F2 sont disposés respectivement sur deux axes y et x, à une même distance D d'un centre O. Le repère O,x,y définit le plan de référence du dispositif. L'axe perpendiculaire au plan de référence et passant par le centre O définit l'axe de rotation O. Le dispositif comprend également deux détecteurs 9, disposés perpendiculairement aux deux axes x et y sur lesquels les foyers FI et F2 deux sources sont disposées. Les deux détecteurs 9 sont disposés à une même distance L du centre O. Ils s'étendent sur une largeur égale à 21 et sur une hauteur égale à 2h. Lorsqu'ils sont irradiés par les faisceaux émis par les sources, ils détectent une intensité atténuée en fonction d'une traversée des rayons X à travers un corps à examiner disposé entre source et détecteur. La figure 1 étant volontairement limitée à la géométrie du dispositif, le corps à examiner n'y est pas représenté, pas plus que le support sur lequel ce corps est disposé. Cependant, le support est monté mobile en rotation autour de l'axe de rotation S2 pour explorer le corps à examiner sous différentes positions angulaires. Alternativement, on monte les deux ensembles source ù détecteur mobiles en rotation autour de l'axe de rotation Q. La figure 1 illustre une géométrie dans laquelle la rotation du support ou des deux ensembles source ù détecteur conduit par construction à des paires de positions angulaires dont les deux positions angulaires de toute paire présentent un écart angulaire A égal à 90 degrés. Non représentés sur la figure 1, le dispositif comprend encore un convertisseur analogiqueùnumérique pour convertir les intensités détectées en données et un ordinateur programmé pour traiter les données et aboutir à une image exprimant un coefficient d'atténuation des rayons X en un ensemble {P0} de points PO du corps à examiner. II convient de rappeler ici que le coefficient d'atténuation en un point PO du corps à examiner est calculé, pour une paire de positions angulaires 0a, 0b par application de la formule suivante : p(b) c(a) n n C n + n 2n2 (k 1ca)+ E1pub)) dans laquelle : c(a) est la donnée provenant de l'intensité détectée au point d'impact T(a) d'un premier rayon i passant par le point considéré PO dans la première position angulaire 0a, p(b) est la donnée provenant de l'intensité détectée au point d'impact R(b) d'un deuxième rayon j passant par le point considéré PO dans la deuxième position angulaire 0b, n k ce est la charge associée au premier point d'impact T(a), définie par k=1 la somme des données c a) provenant des intensités détectées dans la première position angulaire 0a, le long de (n) points Ta) d'une première ligne de calcul de charge qui passe au plus près du premier point d'impact T(a) , et n E p(ub) est la charge associée au deuxième point d'impact R(b) , définie u=1 par la somme des données p(ub) provenant des intensités détectées dans la deuxième position angulaire 0b, le long de (n) points le) d'une deuxième ligne de calcul de charge qui passe au plus près du deuxième point d'impact R(b) Selon l'invention, l'ordinateur est programmé pour pré- calculer, lors d'une première étape, la position sur les détecteurs 9 du point d'impact T du rayon i et du point d'impact R du rayon j, passant par tout point générique G d'un volume irradié par les faisceaux émis par les sources 5. Figure 2, le rayon i issu du foyer F1 et passant par le point générique G de coordonnées (X,Y,Z) frappe l'un des détecteurs 9 au point d'impact T de coordonnées (X1,L,Z1). En ce point d'impact, la donnée provenant de l'intensité détectée par le détecteur est c. L'altitude Z1 se calcule par la formule suivante : Z1 L + D = Z y +D='" De même, figure 3, le rayon j issu du foyer F2 et passant par le point ~o générique G de coordonnées (X,Y,Z) frappe l'autre détecteur 9 au point d'impact R de coordonnées (L,Y2,Z2). En ce point d'impact, la donnée provenant de l'intensité détectée par le détecteur est p. L'altitude Z2 se calcule par la formule suivante : Z2 ù L+D _ k2 Z X +D' 15 Les deux rayons i et j passant par le point générique G se déduisent l'un de l'autre par une rotation de l'un des ensembles source 5 ù détecteur 9 vers l'autre ensemble source 5 ù détecteur 9 autour de l'axe de rotation 0, l'angle de cette rotation étant égal à 90 degrés, c'est à dire à l'écart angulaire 4 entre les deux positions angulaires °a, Ob des différentes paires de positions angulaires du corps à 20 examiner 3 ou des deux ensembles source 5 ù détecteur 9 autour de l'axe de rotation Q. Figure 1, le point générique G (X,Y,Z) et le point d'impact T (X1,L,Z1) du rayon i sur l'un des détecteurs se projettent dans le plan médian respectivement aux points g (X,Y,O) et t (X1,L,0). Le rayon i issu du foyer F1 et passant par le point 25 générique G se projette dans le plan médian pour passer par les points projetés g (X,Y,O) et t (X1,L,0). La coordonnée X1 du point d'impact T se calcule alors par une formule analogue à celle de l'altitude Z1 : L+D =k~ X y +D De même, le point d'impact R (L,Y2,Z2) du rayon j sur l'autre détecteur se projette dans le plan médian au point r (L,Y2,O). Le rayon j issu du foyer F2 et passant par le point générique G se projette dans le plan médian pour passer par les points projetés g (X,Y,O) et r (L,Y2,0). La coordonnée Y2 du point d'impact R se calcule alors par une formule analogue à celle de l'altitude Z2 : Y2 L+D =k2 y x+D La première étape de programmation de l'ordinateur permet ainsi de construire une base de correspondance géométrique entre tout point du volume irradié par les faisceaux émis par les deux sources 5 et les points d'impact sur les deux détecteurs 9, des deux rayons i et j qui passent par ce point en étant issus des foyers FI et F2 des deux sources. Cette base de correspondance est déterminée par la géométrie du dispositif, en l'occurrence par la distance D des foyers FI et F2 à l'axe de rotation S2 et la distance L des détecteurs 9 à cet axe de rotation. Indépendante du corps à examiner et des différentes paires de positions angulaires du corps à examiner ou des deux ensembles source û détecteur autour de l'axe de rotation S2, elle est programmée une fois pour toutes et réutilisée pour tous les points du volume exploré du corps à examiner sous toutes les paires de positions angulaires. A titre d'exemple numérique, on considère la géométrie illustrée par les figures 1 à 3 dans laquelle D = 150 cm et L = 25 cm. Le tableau de la figure 4 indique les rapports entre les coordonnées du point générique G (X,Y,Z) et celles des deux points d'impact T (X1,L,Z1) et R (L,Y2,Z2) sous une succession de positions angulaires e du point générique G allant de 10 degrés en 10 degrés autour de l'axe de rotation S2. L'ordinateur est également programmé pour pré- calculer, lors d'une deuxième étape, la position sur l'un des détecteur 9, de m premières lignes de calcul de charge constituées de n points de détection TP k; k =1-* n; p =1û~ m et sur l'autre détecteur 9, de m deuxièmes lignes de calcul de charge constituées de n points de détection Rq u ; u = û n ; q =1--> m. Le nombre n de points de détection le long d'une ligne de calcul de charge est déterminé par la résolution spatiale des détecteurs 9 dans le sens de la largeur. Le nombre m de lignes de calcul de charges est déterminé par la résolution des détecteurs dans le sens de la hauteur, parallèlement à l'axe de rotation Q. Les premières et deuxièmes lignes de calcul de charge se prolongent pour converger respectivement en un premier S1 et un deuxième S2 points géométriques dont les positions sont déterminées par la distance D des foyers F1 et F2 à l'axe de rotation S2 et la distance L des détecteurs 9 à l'axe de rotation Q. Chaque ligne p ou q de calcul de charge est inclinée par rapport à l'axe de rotation S2 d'un angle de pente yp ou yq qui s'échelonne en une série de valeurs discrètes déterminées par la résolution spatiale du détecteur 9 dans le sens de la hauteur. Figure 5, l'un des détecteurs 9 est représenté de face, sur une demi-hauteur. Il s'agit du détecteur disposé perpendiculairement à l'axe y sur lequel est disposé le foyer F1 de l'une des deux sources. Conventionnellement, on définit le côté gauche, respectivement droit, du détecteur par le côté qui s'étend à gauche, respectivement à droite, de l'axe de rotation û projeté sur le détecteur par le faisceau émis par le foyer FI. Sur la figure 5 sont également représentés, le point d'impact T (X1,L,Z1) du rayon i passant par le point générique G (X,Y,Z) et la ligne de calcul de charge qui passe au plus près de ce point d'impact. Cette ligne indicée p est constituée par les n points de détection TP k; k = 1 - net la charge calculée le long de cette ligne n est égale à E cp k . La ligne p de calcul de charge se prolonge vers la gauche au- k=1 delà du détecteur pour converger vers le point géométrique S1 dont les coordonnées sont (-(D+L),L,O). Comme visible sur la figure 1, ce point géométrique appartient au plan de référence O,x,y et forme l'intersection de la droite reliant les deux foyers F1 et F2 avec le plan du détecteur 9, perpendiculairement à l'axe y. L'inclinaison de la ligne p de calcul de charge est fixée en fonction de la résolution spatiale du détecteur dans le sens de la hauteur. Dans l'exemple illustré par la figure 5, le bord droit du détecteur, soit encore le bord opposé au point géométrique de convergence SI, a été retenu pour fixer l'échelonnement des angles de pente yp des m premières lignes de calcul de charge. Ce choix conduit à des angles de pente les plus faibles possibles et procure le réseau de m lignes de calcul de charge le plus serré. Si l'on retient le bord gauche du détecteur, c'est à dire le bord disposé du même côté que le point géométrique de convergence S1, pour fixer l'échelonnement des angles de pente, on obtient des angles de pente plus forts et un réseau de m lignes de calcul de charge moins serré. 30 De même, figure 6, l'autre des détecteurs 9 est représenté de face, sur une demi-hauteur également. Il s'agit du détecteur disposé perpendiculairement à l'axe x sur lequel est disposé le foyer F2 de l'autre source. Conventionnellement, on définit le côté gauche, respectivement droit, du détecteur par le côté qui s'étend à gauche, respectivement à droite, de l'axe de rotation S2 projeté sur le détecteur par le faisceau émis par le foyer F2. Sur la figure 6 sont également représentés, le point d'impact R (L,Y2,Z2) du rayon j passant par le point générique G (X,Y,Z) et la ligne de calcul de charge qui passe au plus près de ce point d'impact. Cette ligne indicée q est constituée par les n points de détection Rq u; u = 1 ù> n et la charge calculée n le long de cette ligne est égale à pq u . La ligne q de calcul de charge se u=1 prolonge vers la droite au-delà du détecteur pour converger vers le point géométrique S2 dont les coordonnées sont (L,(D+L),O). Comme visible sur la figure 1, ce point géométrique appartient au plan de référence O,x,y et forme l'intersection de la droite reliant les deux foyers F1 et F2 avec le plan du détecteur 9, perpendiculairement à l'axe x. L'inclinaison de la ligne q de calcul de charge est fixée en fonction de la résolution spatiale du détecteur dans le sens de la hauteur. Dans l'exemple illustré par la figure 6, le bord gauche du détecteur, soit encore le bord opposé au point géométrique de convergence S2, a été retenu pour fixer l'échelonnement des angles de pente yq des m deuxièmes lignes de calcul de charge. Ce choix conduit à des angles de pente les plus faibles possibles et procure le réseau de m lignes de calcul de charge le plus serré. Si l'on retient le bord droit du détecteur, c'est à dire le bord disposé du même côté que le point géométrique de convergence S2, pour fixer l'échelonnement des angles de pente, on obtient des angles de pente plus forts et un réseau de m lignes de calcul de charge moins serré. Si l'on appelle Rh la résolution spatiale des détecteurs dans le sens de la hauteur, l'angle de pente yp de la ligne p des m premières lignes de calcul de charge est donné, dans l'hypothèse où le bord droit de l'un des détecteurs a été retenu pour fixer l'échelonnement des angles de pente, par la formule yp=Arctg( pxRh );P=1-*m= h D+L+.e Rh 5 10 De même, l'angle de pente 7q de la ligne q des m deuxièmes lignes de calcul de charge est donné, dans l'hypothèse où le bord gauche de l'autre détecteur a été retenu pour fixer l'échelonnement des angles de pente, par la formule rq= ùArctg( gXRh );p=1ù> m = h D+L+.e Rh Pour pré- construire les premières et deuxièmes lignes de calcul de charge, on approche chaque ligne par une succession de marches d'escalier. La hauteur des marches est fixée à la valeur de la résolution spatiale Rh dans le sens de la hauteur. La largeur des marches est donnée par la formule D+L+e ou D+L+.e .ep= = p eq q Si l'on appelle Rt la résolution spatiale des détecteurs dans le sens de la largeur, le nombre de points de détection, ou pixels, par marche est égal à ou 15 Np_ .ep R� Nq = .eq RI La largeur des marches ainsi que le nombre de points de détection par marche diminuent avec l'augmentation de l'angle de pente des lignes de calcul de charge. Les m premières lignes et les m deuxièmes lignes de calcul de charge 20 définies précédemment marquent la trace sur les deux détecteurs de m plans inclinés par rapport à l'axe de rotation O. Tous ces plans se déduisent les uns des autres par une rotation autour de la droite reliant les deux foyers F1 et F2. L'écart angulaire entre deux plans inclinés est fixé en fonction de la résolution spatiale des détecteurs suivant la hauteur Rh. Un plan incliné coïncide avec le plan porté par 25 les deux rayons i et j passant par le point générique G (X,Y,Z) lorsque les deux points d'impact T (Xl,L,Zl) et R (L,Y2,Z2) des deux rayons i et j tombent sur les deux lignes de calcul de charge indicées p et q qui marquent la trace du plan incliné sur les deux détecteurs. Autrement dit, à la résolution spatiale près des détecteurs dans le sens de la hauteur, les lignes de calcul de charge marquent la trace sur les deux détecteurs des plans portés par tout couple de rayons i et j issus respectivement des foyers F1 et F2 et se coupant en tout point générique G du volume irradié. Pour un tel couple de rayons i et j, la ligne p de calcul de charge qui passe par le point d'impact T (X1,L,Z1) du rayon i et la ligne q de calcul de charge qui passe par le point d'impact R (L,Y2,Z2) du rayon j appartiennent au même plan incliné, ont le même indice p ou q et comme le montre les figures 4 et 5, le même angle de pente en valeur absolue : x h p=q ; yp=yq=Arctg( p+ R~) ;p=1--> m=Rh A titre d'exemple numérique, on considère la géométrie illustrée par les figures 1 à 5 dans laquelle D = 150 cm ; L = = h = 25 cm ; Rh = Rt = 81 �m. Le tableau de la figure 7 indique, pour différents angles de pente yp ou yq, les largeurs de marche ,ep ou eq , les nombres N p ou Nq de points de détection par marche et 15 le nombre E de marches par escalier. Le tableau indique également la hauteur atteinte par les deux lignes de calcul de charge, respectivement sur le bord droit de l'un des détecteurs et sur le bord gauche de l'autre détecteur, et l'indice p ou q de ces deux lignes, obtenu par arrondi à la valeur entière la plus proche de la hauteur atteinte rapportée à la résolution spatiale Rh du détecteur dans le sens de la 20 hauteur. L'indice p ou q est ainsi égal à la hauteur en points de détection ou en pixels, atteinte par les deux lignes de calcul de charge respectivement sur le bord droit de l'un des détecteurs et sur le bord gauche de l'autre détecteur. La base de correspondance géométrique pré- calculée à l'étape (1) et les m premières et m deuxièmes lignes de calcul de charge pré- construites à l'étape (2) 25 sont utilisées à l'étape (3) de programmation de l'ordinateur pour calculer le coefficient d'atténuation C en n'importe quel point PO du corps à examiner. Il faut considérer à présent que le point PO du corps à examiner, par exemple un point du rein d'un individu, occupe sur la figure 1, la position du point générique G dans le volume irradié par les faisceaux. Le corps à examiner se trouve 30 simultanément dans la position angulaire 0a vis-à-vis de l'ensemble foyer F1 ù détecteur et dans la position angulaire 0b vis-à-vis de l'ensemble foyer F2 ù10 détecteur. L'ordinateur recherche dans la base de correspondance géométrique les premier T(a) et deuxième R(b) points d'impact du point P0, appelle la donnée c(a) provenant de l'intensité détectée au premier point d'impact T(a) sous la première position angulaire 6a, appelle la donnée p(b) provenant de l'intensité détectée au deuxième point d'impact R(b) sous la deuxième position angulaire Ob, recherche la ligne p de calcul de charge T pa); k = 1 n qui passe le plus près du n premier point d'impact T(a) et calcule la charge E c(a~ associé à ce point k=1 p' d'impact, recherche la ligne q de calcul de charge Re ;U = 1 -* n qui passe le plus n près du deuxième point d'impact R(b) et calcule la charge E pgbu associée à ce u=1 point d'impact, et finalement calcule la valeur du coefficient d'atténuation C au point PO considéré sous la paire de positions angulaires considérées ea et 6b. De préférence, l'ordinateur 13 est programmé pour effectuer l'étape (3) à l'aide de deux étapes supplémentaires suivantes. D'une part, la position sur l'un des détecteurs 9 d'un premier point de référence T* est pré- calculée. Ce point de référence T* (t ,L, Z*1) est, figure 5, à l'intersection du bord droit du détecteur, pris pour échelonner les valeurs discrètes des angles de pente yp des m premières lignes de calcul de charge TP k; k = 1 n; p = 1 m et de la droite passant par le premier point géométrique de convergence S1 (-(D+L)L,0) et par le premier point d'impact T (X1,L,Z1). Au premier point d'impact T est associée la charge calculée le long de la ligne indicée p de calcul de charge qui passe au plus près du point de référence T* compte tenu de la résolution spatiale Rh du détecteur dans le sens de la hauteur. D'autre part, la position sur l'autre détecteur 9 d'un deuxième point de référence R* est également pré- calculée. Ce point de référence R* (- ,L, Z*2) est, figure 6, à l'intersection du bord gauche de l'autre détecteur 9, pris pour échelonner les valeurs discrètes des angles de pente Yb des m deuxièmes lignes de calcul de charge Rq u; u =1-> n ; q =1- m et de la droite passant par le deuxième point géométrique de convergence S2 (D+L,L,O) et par le deuxième point d'impact R (L,Y2,Z2). Au deuxième point d'impact R est associée la charge calculée le long de la ligne indicée q de calcul de charge qui passe par le point de référence R* compte tenu de la résolution spatiale Rh du détecteur dans le sens de la hauteur. Etant donné que les deux points d'impact T (X1,L,Z1) et R (L,Y2,Z2) appartiennent au même plan incliné, les indices p et q des deux lignes de charge sont égaux de même que l'altitude Z*1 du premier point de référence T et celle Z*2 du deuxième point de référence R* . Cette altitude est donnée par la formule : Ê+L+D Z Z X+Y+D En poursuivant l'exemple numérique précédent, le tableau de la figure 8 indique les rapports entre les coordonnées du point générique G (X,Y,Z), celles des deux points d'impact T (X1,L,Z1) et R (L,Y2,Z2), l'altitude Z* des point de référence correspondant T* et R et l'indice p ou q des deux lignes de calcul de charge passant au plus près des deux points de référence, pour une succession de positions angulaires O du point générique G allant de 10 degrés en 10 degrés autour de l'axe de rotation S2. L'indice p ou q est obtenu par arrondi à la valeur entière de la plus proche de l'altitude Z* rapportée à la résolution spatiale Rh du détecteur dans le sens de la hauteur. Il faut relever que le calcul du coefficient d'atténuation C est effectué pour des points d'impact et des lignes de calcul de charge associées à ces points qui tous appartiennent au plan porté par les deux rayons i et j se coupant au point PO considéré, à la résolution spatiale près des détecteurs suivant le sens de la hauteur. Ce plan prenant selon le point PO considéré dans le volume d'irradiation, n'importe quelle inclinaison vis-à-vis de l'axe de rotation S2, le procédé selon l'invention s'affranchit d'un calcul du coefficient d'atténuation limité à un plan de coupe perpendiculaire à l'axe de rotation et apporte ainsi une imagerie par spatiodensitométrie. A l'étape (4), le calcul du coefficient d'atténuation C au point PO pour la paire des positions angulaires ea et Ob est recommencé pour les différentes paires de positions angulaires obtenues par la rotation simultanée des deux ensembles source ù détecteur autour de l'axe de rotation û. Alternativement, les différentes paires de positions angulaires sont obtenues en laissant fixes les deux ensembles source ù détecteur et en faisant tourner le corps à examiner et son support autour de l'axe de rotation S2. Enfin à l'étape (5), on recommence le calcul de la valeur moyenne du coefficient d'atténuation en n'importe quel point du corps à examiner inscrit dans le volume d'irradiation. Le procédé selon l'invention permet ainsi de suivre un point PO quelconque du corps à examiner inscrit dans le volume d'irradiation pour aboutir à une valeur moyenne du coefficient d'atténuation en ce point. De préférence, l'ordinateur 13 est programmé pour effectuer l'étape (5) à l'aide des étapes supplémentaires suivantes. Premièrement, les termes d'un premier et d'un deuxième vecteur de charge sont pré- calculés pour les deux positions angulaires 0a et 0b de chacune des différentes paires de positions angulaires obtenues par la rotation des deux ensembles source ù détecteur. Pour chaque paire, les m termes du premier vecteur et du deuxième vecteur de charge sont donnés par les charges k cW et u et)j k=1 p' u=1 p le long des m premières et des m deuxièmes lignes de calcul de charge. Deuxièmement, lors du calcul du coefficient d'atténuation de n'importe quel point PO de l'ensemble considéré, les termes des deux vecteurs de charge pré-calculés le long des deux lignes de calcul de charge qui passent au plus près des deux points d'impact T(a) et R(b) du point P0, sous chacune des différentes paires de positions angulaires, sont simplement appelés sans être recalculés d'un point PO à un autre point PO de l'ensemble considéré. Par ces étapes supplémentaires, on simplifie l'exécution de l'étape (5) de deux façons. D'une part, on calcule une fois pour toutes les charges de tous les points PO dont les deux points d'impact T(a) et R(b) appartiennent à un même plan incliné et possèdent ainsi le même terme du premier vecteur de charge et le même terme du deuxième vecteur de charge. D'autre part, on appelle ou rappelle les termes pré- calculés des vecteurs de charge pour tout autre ensemble de points considérés, par exemple lors de l'examen d'une autre partie du corps d'un individu, ou pour le même ensemble de points considérés mais à un autre moment de la vie de l'individu. A noter que les deux détecteurs 9 reçoivent des rayons i et j dont la longueur varie en fonction de la position des points d'impact par rapport à la position des foyers F1 et F2. Une correction est apportée aux données provenant des intensités détectées, en majorant la valeur en un point d'impact par référence au point central M en la multipliant par le rapport du carré des longueurs du foyer FI ou F2 au point d'impact à la longueur minimale, c'est-à-dire celle qui relie le foyer F1 ou F2 au point central M. Il faut également tenir compte de l'allongement du rayon qui crée une augmentation de la dose apparente. On prendra comme coefficient de correction soit la racine carrée de la valeur indiquée ci-dessus, soit plus simplement une valeur égale à 1+e/2 au lieu de 1+e. L'exposé de l'invention est poursuivi en référence à un dispositif d'imagerie à rayons X ne comportant qu'un ensemble source û détecteur. Figure 9, le foyer FI de la source 5 est disposé sur l'axe y, à la distance D du centre O. Le détecteur 9 est disposé perpendiculairement à l'axe y à la distance L du centre O. Il s'étend sur une largeur égale à 21 et sur une hauteur égale à 2h. Bien que non représenté sur la figure 9, le support destiné à recevoir le corps à examiner est monté mobile en rotation autour de l'axe de rotation Q pour explorer le corps à examiner sous différentes positions angulaires O. Alternativement, on monte l'ensemble source û détecteur mobile en rotation autour de l'axe de rotation S2. Le rayon i issu du foyer FI et passant par le point générique G de coordonnées (X;,Y;,Z) frappe le détecteur 9 au point d'impact T de coordonnées (XI,L,Z1). En ce point d'impact, la donnée provenant de l'intensité détectée par le détecteur est c. L'altitude Z1 se calcule par la formule suivante : Z1 _ L+D Z Y; + D=k1 De même, après la rotation d'angle A du point générique G autour de l'axe de rotation ç2, le rayon j issu du foyer F1 et passant par le point générique G de coordonnées (Xi,YY,Z) frappe le détecteur 9 au point d'impact R de coordonnées (X2,L,Z2). En ce point d'impact, la donnée provenant de l'intensité détectée par le détecteur est p. L'altitude Z2 se calcule par la formule suivante : Z2 ù L+D _ k2 Z Yj +D Les deux rayons i et j passant par le point générique G se déduisent l'un de l'autre par la rotation du point générique G autour de l'axe de rotation S, l'angle de cette rotation étant égal à 90 degrés, c'est à dire à l'écart angulaire A entre les deux positions angulaires 0a et 0b des différentes paires de positions angulaires du corps à examiner 3 autour de l'axe de rotation S2. Figure 9, le point générique G (X;,Y;,Z) et le point d'impact T (X1,L,Z1) du rayon i sur le détecteur se projettent dans le plan de référence respectivement aux points g (X;,Y;,O) et t (X1,L,0). Le rayon i issu du foyer F1 et passant par le point générique G se projette dans le plan médian pour passer par les points projetés g (X;,Y;,O) et t (X1,L,0). La coordonnée X1 du point d'impact T se calcule alors par une formule analogue à celle de l'altitude Z1 : X1 _ L+D =1(1 X; y;+D De même, figure 10, le point générique G (X;,Yi,Z) et le point d'impact R (X2,L,Z2) du rayon j sur le détecteur se projettent dans le plan médian respectivement aux points g (Xi,YY,O) et r (X2,L,0). Le rayon j issu du foyer F2 et passant par le point générique G se projette dans le plan de référence pour passer par les points projetés g (X,Y,O) et r (X2,L,0). La coordonnée Y2 du point d'impact R se calcule alors par une formule analogue à celle de l'altitude Z2 : X2 L+D =k2 xi yi+D Figure 11, le détecteur 9 est représenté de face, sur une demi-hauteur. Conventionnellement, on définit le côté gauche, respectivement droit, du détecteur par le côté qui s'étend à gauche, respectivement à droite, de l'axe de rotation S2 projeté sur le détecteur par le faisceau émis par le foyer F1. Le rayon i passe par le point générique G (X;,Y;,Z) et frappe le détecteur au point d'impact T (X1,L,Z1). La ligne de calcul de charge indicée p qui passe au plus près de point d'impact est constituée par les n points de détection T(3' ; k =1 ù> n. La ligne de calcul de charge se prolonge vers la gauche au-delà du détecteur pour converger vers le point géométrique S1 dont les coordonnées sont (-(D+L),L,O). De préférence, la position sur le détecteur 9 d'un premier point de référence T* est pré- calculée. Ce point de référence T* (.e ,L, Z*1) est, figure 11, à l'intersection du bord droit du détecteur, pris pour échelonner les valeurs discrètes des angles de pente yp des m premières lignes de calcul de charge T(4 ; k =1-* n; p =1ù> m et de la droite passant par le premier point géométrique de convergence S1 (-(D+L)L,O) et par le premier point d'impact T (X1,L,Z1). Au premier point n d'impact T est associée la charge E cW calculée le long de la ligne de calcul de k=1 charge indicée p qui passe au plus près du point de référence T* compte tenu de la 10 résolution spatiale Rh du détecteur dans le sens de la hauteur. De même, Le rayon j passe par le point générique G (Xi,Yi,Z) et frappe le détecteur au point d'impact R (X2,L,Z2). La ligne de calcul de charge indicée q qui passe au plus près de point d'impact est constituée par les n points de détection RW; u = -+ n . La ligne de calcul de charge se prolonge vers la droite au-delà du 15 détecteur pour converger vers le point géométrique S'1 dont les coordonnées sont (D+L,L,O). De préférence, la position le détecteur 9 d'un deuxième point de référence R* est pré- calculée. Ce point de référence R* (-.r? ,L, Z*2) est, figure 11, à l'intersection du bord gauche du détecteur, pris pour échelonner les valeurs discrètes des angles de pente yq des m deuxièmes lignes de calcul de charge 20 RW; u =1ù> n ; q =1ù> m , avec la droite passant par le deuxième point géométrique de convergence S'1 (D+L,L,O) et par le deuxième point d'impact R n (X2,L,Z2). Au deuxième point d'impact R est associée la charge E pq'u calculée u=1 le long de la ligne de calcul de charge indicée q qui passe au plus près du point de référence R* compte tenu de la résolution spatiale Rh du détecteur dans le sens 25 de la hauteur. A la résolution spatiale près du détecteur dans le sens de la hauteur, les m premières et m deuxièmes lignes de calcul de charge marquent la trace sur le détecteur de m plans portés par tout couple de rayons i et j issus du foyers F1 passant par le point générique G et se déduisant l'un de l'autre par la rotation 30 d'angle A du point générique G autour de l'axe de rotation O. Pour un tel couple de rayons i et j, la ligne p de calcul de charge qui passe par le point d'impact T (X1,L,Z1) du rayon i et la ligne q de calcul de charge qui passe par le point d'impact R (X2,L,Z2) du rayon j appartiennent au même plan incliné, ont le même indice p ou q, le même angle de pente en valeur absolue, comme le montre la figure 11. De même, les deux points de références T*V ,L, Z*1) et R* (-t ,L, Z*2) ont la même altitude Z*1 ou Z*2. L'exposé de l'invention est poursuivi par la description de deux dispositifs spécialement conçus pour la mise en oeuvre du procédé d'imagerie à rayons X. Le premier dispositif, figure 12, est un scanner dit � intégral � permettant en un seul passage d'obtenir les images initiales d'un patient de la tête aux pieds. Ce ~o dispositif comprend dans une cabine 17 de 220 cm de hauteur : - un support 1 pour recevoir un patient 3 à examiner, monté mobile en rotation autour d'un de rotation n, - deux sources 5 alimentant deux foyers F émettant deux faisceaux 7 de rayons X, 15 - deux détecteurs 9 irradiés par les faisceaux pour détecter une intensité atténuée en fonction d'une traversée des rayons X à travers le patient 3, et - un convertisseur analogiqueùnumérique 11 pour convertir les intensités détectées en données pour aboutir à une image exprimant un coefficient d'atténuation des rayons X en un ensemble de points PO du patient. 20 Les deux sources 5 et les deux détecteurs 9 sont fixés à des supports 19 de la cabine 17 et forment deux ensembles superposés suivant la hauteur du dispositif, c'est-à-dire parallèlement à l'axe de rotation O. Ils réalisent chacun une géométrie telle qu'elle a été décrite en référence aux figures 6 à 8. Ainsi les foyers F sont disposés sur les axes x perpendiculairement aux détecteurs 9. Ces derniers 25 s'étendent sur une hauteur 2h égale à 110 cm. Ils comprennent une ou plusieurs barrettes de détection 21 montées mobiles en translation parallèlement aux axes y pour couvrir une zone de détection égale en largeur 21 à 60 cm. La distance D entre les foyers F et l'axe de rotation O est égale à 150 cm. La distance L entre les détecteurs 9 et l'axe de rotation S2 est égale à 30cm. 30 Le patient 3 se tient debout sur le support 1 et est immobilisé par des tubes transparents verticaux non représentés sur la figure 12. Les deux sources 5 sont par exemple constituées par deux tubes à rayons X à grand angle d'ouverture, par exemple 40° ou 56 °. Si l'on prend par exemple un angle d'ouverture totale de 40°, à 180 cm de distance, une hauteur de 131 cm peut être irradiée. En plaçant les deux tubes à des hauteurs convenables, on peut irradier une hauteur supérieure à 2 m avec un chevauchement. Ce chevauchement pourra être annulé par des procédés informatiques, ou réduit au minimum en plaçant une plaque horizontale en acier et plomb à mi-hauteur du détecteur. On laissera un léger espace entre cette plaque et le détecteur de sorte à limiter le chevauchement à un petit nombre de lignes facilement détectables. Sur chaque tube est fixé un réverbérateur 15 en tronc de pyramide par exemple. Devant ce réverbérateur coulisse une plaque, non représentée sur la figure 12, étanche aux rayons X comprenant une fente verticale dont la largeur sera par exemple égale à 2,5 mm. On déplace simultanément le détecteur 9 et la plaque à fente verticale pour limiter l'irradiation à la barrette de détection du détecteur dont la dimension est de l'ordre de 6,6 mm. Dans ces conditions, la dose de rayons X reçue par le patient 3 est limitée par le temps de passage de la fente verticale. En pratique, aucun point du patient 3 n'est atteint pendant un temps supérieur à 125 millisecondes. Les expériences ont montré que pour une hauteur de 46 cm, la dose reçue par le patient est largement inférieure à 100 micro grays par prise de vue. Le fonctionnement est le suivant. On place le patient 3 dans une première position angulaire 01 et on prend une image intégrale en activant les deux détecteurs et les deux tubes, ou une image partielle en activant un détecteur et un tube. Si le détecteur comprend plusieurs barrettes de détection 21, on déplace une ou plusieurs de ces barrettes en fonction de la zone du patient à explorer. On fait tourner le patient de 10° ou 18° et l'on prend une deuxième image dans un deuxième position angulaire 02 et ainsi de suite jusqu'à obtenir 36 ou 20 images décalées régulièrement sur 360°. Ces 36 images serviront au calcul des valeurs des coefficients d'atténuation en tout point PO du patient choisi à l'avance selon le procédé décrit précédemment. Dans la pratique, on réalisera par exemple un ensemble de 36 images de référence de l'ensemble du patient qui pourront être conservées pour des usages ultérieurs. Ces images auront une définition de base de 81 microns, c'est-à-dire une résolution spatiale Rh en hauteur et une résolution spatiale Re en largeur égales à 81 microns. On pourra aussi réaliser une image de l'ensemble du patient à 162 microns de définition et, pour un volume restreint, par exemple l'orbite, le coeur, le foie, le rein etc..., réaliser une image à 81 microns. Avec un système informatique comprenant 40 PC mis en parallèle, on peut traiter 125 méga pixels par seconde, ce qui autorise de scanner un patient à 1 mm de définition en moins d'une seconde. A 81 microns, le temps de calcul global sera de l'ordre de 1128 secondes. A 162 microns, il sera de l'ordre de 141 secondes. De manière standard, on pourra ainsi scanner un individu de la tête aux pieds en 4 minutes de prise de vue à 81 microns et restituer une image à 162 microns en moins de 2 minutes. L'examen d'une partie du patient ayant un volume de 125 cm3, par exemple une orbite oculaire, le temps de calcul à 81 microns sera de l'ordre de 1,6 seconde. En définitive, la réalisation d'un scanner intégral permettra au radiologue de travailler dans les conditions suivantes : - prise de vue en 4 minutes, ce qui autorise de traiter 6 à 10 patients par heure de manière acceptable mais qui suppose une bonne planification des opérations, - restitution d'une image totale à 162 microns, en moins de 2 minutes, - restitution d'une image partielle, par exemple l'orbite, le coeur, le foie ou l'articulation du genou, avec possibilité d'amplifier les contrastes, dans des temps de quelques dizaines de seconde. La concomitance pratique de la prise de vue et du calcul autorise le radiologue à rédiger son rapport en temps réel, en se fondant sur des images à haute définition globales ou partielles selon les besoins. Par exemple, pour vérifier la présence d'un cancer secondaire dans les os, on réalise une image intégrale à 162 microns, sans amplification de contraste. Le dispositif qui vient d'être décrit comprend un ordinateur 13 programmé pour effectuer les étapes du procédé selon l'invention. Plus particulièrement, l'ordinateur 13 est pourvu d'une mémoire de masse pour stocker les positions des deux points d'impacts T et R de tout point générique G d'un volume irradié pouvant contenir un être humain. Si l'on considère un cylindre ayant une hauteur de 2 m et un rayon de 30 cm, à la définition de 81 microns, le nombre de positions du point générique est égal à 1,06 1012 et le nombre de points d'impact à 2,12 1012. En fonction de la partie du patient dont l'imagerie X est recherchée, l'ordinateur prélève dans la mémoire de masse l'information pertinente, savoir les positions des points d'impact T et R des différents points PO de la partie du patient à examiner, et les transfère à une mémoire à accès direct de l'ordinateur. Le deuxième dispositif, figure 13, est un scanner pour l'imagerie X du sein. Ce dispositif comprend : - une table d'examen 1 pour recevoir une patiente 3 à examiner, munie d'une ouverture 23 pour permettre le logement du sein de la patiente dans une coquille 25 transparente aux rayons X et destinée à comprimer légèrement le sein, ou à l'aspirer légèrement pour que la peau vienne au contact de la coquille et ainsi immobiliser le sein, - une source 5 alimentant un foyer F émettant un faisceau 7 de rayons X, - un détecteur 9 irradié par le faisceau pour détecter une intensité atténuée en fonction d'une traversée des rayons X à travers le sein du patient 3, et - un convertisseur analogiqueûnumérique 11 pour convertir les intensités détectées en données pour aboutir à une image exprimant un coefficient d'atténuation des rayons X en un ensemble de points PO du sein de la patiente. La table d'examen 1 est portée par deux supports 27. La source 5 et le détecteur 9 forment un ensemble mobile autour d'un axe de rotation SI Ils sont fixés à une cage 29 de forme cylindrique à base circulaire 31 actionnée en rotation par l'intermédiaire d'un arbre 33 et d'un moteur 35. La cage comprend une paroi supérieure, une paroi inférieure et une paroi latérale fabriquées en plomb pour être imperméables aux rayons X et empêcher leur propagation à l'extérieur. Ils réalisent une géométrie telle qu'elle a été décrite en référence aux figures 6 à 8, dans laquelle l'ensemble source û détecteur est mobile et le support fixe. La source 5 est par exemple un tube de faible puissance, typiquement de quelques milliampères sous une tension de 75 kV. Le foyer F est disposé sur l'axe x perpendiculairement au détecteur 9. Ce dernier est constitué d'une plaque de détection 37 de 23 cm de hauteur et de 30 cm de largeur et possède une définition de 162 microns. La distance D entre le foyer F et l'axe de rotation S2 est égale à 37,5 cm. La distance L entre le détecteur 9 et l'axe de rotation S2 est aussi égale à 37,5 cm. Dans l'exemple de la figure 13, le foyer F est disposé de préférence dans un plan perpendiculaire à l'axe de rotation S2 et passant par le sommet S du détecteur 9 disposé juste au dessous de la table d'examen 1. Cet agencement en demi-champ permet d'irradier le sein sur une grande hauteur sans trop augmenter l'angle d'ouverture du réverbérateur 15 fixé au tube à rayons X. Mais on peut également disposer le foyer F dans le plan médian, c'est-à-dire dans le plan perpendiculaire à l'axe de rotation S2 et passant à mi-hauteur du détecteur. Le fonctionnement est le suivant. Dans une première position angulaire 01 de l'ensemble source û détecteur, on prend une première image. On fait tourner la cage de 10° ou 18° et l'on prend une deuxième image dans un deuxième position angulaire 02 et ainsi de suite jusqu'à obtenir 36 ou 20 images décalées régulièrement sur 360°. Ces 36 images serviront au calcul des valeurs des coefficients d'atténuation en tout point PO du sein de la patiente. On réalisera par exemple un ensemble de 36 images de référence de l'ensemble du sein qui pourront être conservées pour des usages ultérieurs. A noter que chaque image ayant, dans l'exemple cité précédemment, une largeur de 30 cm, les 36 images se chevauchent. Le dispositif qui vient d'être décrit comprend un ordinateur 13 programmé pour effectuer les étapes du procédé selon l'invention. Plus particulièrement, l'ordinateur 13 est pourvu d'une mémoire à accès rapide pour stocker les positions des deux points d'impacts T et R de tout point générique G d'un volume irradié correspondant à un sein. Si l'on considère une coquille en demi-sphère ayant un rayon de 15 cm, à la définition de 162 microns, le nombre de positions du point générique est égal à 2,8 108 et le nombre de points d'impact à 5,6 108. L'ordinateur prélève dans la mémoire à accès rapide l'information pertinente, savoir les positions des points d'impact T et R des différents points PO du sein de la patiente et les transfère à une mémoire à accès direct de l'ordinateur. Lorsqu'on applique le procédé selon l'invention à un organe de petite dimension comme le foie ou le coeur, des tissus proches peuvent conduire à des données voisines. Dans ces conditions, il est intéressant de réaliser un zoom de contraste basé sur le principe suivant : on considère la valeur du coefficient d'absorption en un point et on l'élève à une puissance comprise entre 1,2 et 1,8 par exemple. On illustre l'amplification de contraste comme suit : le niveau d'absorption le plus bas est égal à 10, le niveau le plus élevé est égal à 100, le contraste est égal à 10 sans amplification, et si l'on emploie un exposant égal à 1,5 le niveau le plus bas devient égal à 31,62 le niveau le plus haut égal à 1000 et l'amplification égale à 31,62. 20 REVENDICATIONS1. Procédé d'imagerie à rayons X, du type scanographie, dans lequel, un corps à examiner (3) étant reçu par un support (1), - on émet, à partir du foyer (F1,F2) d'une source (5), un faisceau (7) de rayons X (i,j) en direction du corps à examiner (3), - on détecte, à l'aide d'un détecteur (9) irradié par le faisceau, une intensité atténuée en fonction d'une traversée des rayons X à travers le corps à examiner (3), sous différentes paires de positions angulaires (6a,Ob) du corps à examiner (3) ou de ~o l'ensemble source (5) û détecteur (9) autour d'un axe de rotation (û), les deux positions angulaires (ea,Ob) des différentes paires présentant un écart angulaire (A) égal de préférence à 90 degrés, - on convertit, à l'aide d'un convertisseur (11) analogiqueûnumérique, les intensités détectées en données (c(a) , p(b)) et, 15 - on traite, à l'aide d'un ordinateur programmé (13), les données (C(a), p(b� pour aboutir à une image exprimant un coefficient d'atténuation des rayons X en un ensemble de points {P0} du corps à examiner, en calculant, pour une paire de positions angulaires (ea,Ob), une valeur (C) du coefficient d'atténuation en un point (P0) de l'ensemble {P0} par application de la formule suivante : P C(a) 1 n n C n + n 2n2 (k=1C r+u=lpub� dans laquelle : c(a) est la donnée provenant de l'intensité détectée au point d'impact T(a) d'un premier rayon i passant par le point considéré (P0) dans la première position angulaire (Oa), 25 p(b) est la donnée provenant de l'intensité détectée au point d'impact R(b) d'un deuxième rayon j passant par le point considéré (P0) dans la deuxième position angulaire (Ob),n E ce est la charge associée au premier point d'impact (T(a)), définie par la k=1 somme des données (ce ) provenant des intensités détectées dans la première position angulaire (ea), le long de (n) points (T a)) d'une première ligne de calcul de charge qui passe au plus près du premier point d'impact (T(a)), et n p(ub) est la charge associée au deuxième point d'impact (R(b) ), définie par la u=1 somme des données (pub)) provenant des intensités détectées dans la deuxième position angulaire (Ob), le long de (n) points (Rh d'une deuxième ligne de calcul de charge qui passe au plus près du deuxième point d'impact (R(b) ) caractérisé en ce que l'ordinateur (13) est programmé pour effectuer les étapes suivantes : (1) pré- calculer, pour tout point générique (G) d'un volume irradié par le faisceau émis par la source (5), la position sur le détecteur (9) d'un premier point d'impact (T) d'un rayon (i) et d'un deuxième point d'impact (R) d'un rayon (j), les deux rayons (i) et (j) passant par le point générique (G) et se déduisant l'un de l'autre par une rotation du point générique (G) ou du de l'ensemble source (5) ù détecteur (9) autour de l'axe de rotation (S2) dont l'angle est égal à l'écart angulaire (4) entre les deux positions angulaires (ea,eb) des différentes paires de positions angulaires du corps à examiner (3) ou de l'ensemble source (5) ù détecteur (9) autour de l'axe de rotation (e), (2) pré- calculer sur le détecteur (9) la position de (m) premières et de (m) deuxièmes lignes de calcul de charge de (n) points de détection (Tp,k;k=lù>n;p=1ù>m; Rq,u;u=l-+n;q=1ù>m) se prolongeant pour converger respectivement en un premier (Si) et un deuxième (S2) points géométriques dont les positions sont déterminées par la distance (D) du foyer (F1,F2) à l'axe de rotation (Q) et la distance (L) du détecteur (9) à l'axe de rotation (Q), chaque ligne (p,q) étant inclinée par rapport à l'axe de rotation (S2) d'un angle de pente (Yp,7q) qui s'échelonne en une série de valeurs discrètes déterminées par une résolution spatiale (Rh) du détecteur (9) parallèlement à l'axe de rotation (S2),(3) effectuer le calcul de la valeur (C) du coefficient d'atténuation en considérant, en fonction de la position du point (P0) dans le volume irradié par le faisceau, les positions pré- calculées des premier et deuxième points d'impact (T(a) ; R(b)) et les première (p) et deuxième (q) lignes de calcul de charge pré- construites (T (pa~; k =1--> n ; RW ;U = n) qui passent au plus près des premier et deuxième points d'impact (T(a) ; R(b) ) (4) recommencer l'étape (3) pour les différentes paires de positions angulaires (8a,Ob) du corps à examiner (3) ou de l'ensemble source (5) - détecteur (9) autour de l'axe de rotation (û) et calculer une valeur moyenne du coefficient d'atténuation (C) au point considéré (P0) pour les différentes paires de positions angulaires, et (5) recommencer l'étape (4) et afficher à l'écran de l'ordinateur l'image exprimant les valeurs moyennes de l'ensemble {P0} des points (P0). 2. Procédé d'imagerie à rayons X selon la revendication 1, caractérisé en ce que l'ordinateur (13) est programmé pour effectuer l'étape (3) à l'aide des étapes suivantes : (6) pré- calculer la position sur le détecteur (9) d'un point de référence (T ) du premier point d'impact (T ), défini par l'intersection de la droite passant par le premier point géométrique de convergence (SI) et par le premier point d'impact (T) et une ligne de points du détecteur, parallèle à l'axe de rotation û, prise pour échelonner les valeurs discrètes des angles de pente (yp) des (m) premières lignes de calcul de charge (T P k; k =1-* n; p =1--> m) et associer au premier point d'impact (T ), la charge calculée le long de la ligne de calcul de charge qui passe au plus près du point de référence (T ), et (7) pré- calculer la position sur le détecteur (9) d'un point de référence (R ) du deuxième point d'impact (R ), défini par l'intersection de la droite passant par le deuxième point géométrique de convergence (S2) et par le deuxième point d'impact (R) et une ligne de points du détecteur, parallèle à l'axe de rotation û, prise pour échelonner les valeurs discrètes des angles de pente (7q) des (m) deuxièmes lignes de calcul de charge (Rq u ; u = n ;q = 1 ù> m) et associer au deuxième point d'impact (R ), lacharge calculée le long de la ligne de calcul de charge qui passe au plus près du point de référence (R ). 3. Procédé d'imagerie à rayons X selon la revendication 1 ou 2, caractérisé en ce que l'ordinateur est programmé pour effectuer l'étape (2) ou les étapes (6) et (7) en considérant les points de détection alignés parallèlement à l'axe de rotation (û) sur un bord du détecteur (9) opposé au point géométrique de convergence (S1,S2) pour échelonner les valeurs discrètes des angles de pente (7p,yq) des (m) premières et (m) deuxièmes lignes de calcul de charge. 4 Procédé d'imagerie à rayons X selon la revendication 1 ou 2, caractérisé en ce que l'ordinateur (13) est programmé pour effectuer l'étape (5) à l'aide des étapes suivantes : les termes d'un premier et d'un deuxième vecteur de charge, donnés par les charges (E cp(a~ , pgbu ) le long des (m) premières et des (m) k=1 p' u=1 deuxièmes lignes de calcul de charge, sont pré- calculés pour les deux positions angulaires (Oa,Ob) de chacune des différentes paires de positions angulaires, et (9) le coefficient d'atténuation (C) de n'importe quel point (P0) de l'ensemble considéré est calculé en appelant les termes des deux vecteurs de charge pré- calculés le long des deux lignes de calcul de charge qui passent au plus près des deux points d'impact (T(a) , R(b)) du point (P0), sous chacune des différentes paires de positions angulaires (0a,eb). 5. Procédé d'imagerie à rayons X selon l'une des revendications 1 à 4, caractérisé en ce que la rotation d'angle (A) définie à l'étape (1) correspond à deux positions angulaires (Oa, Ob) occupées successivement dans le temps du corps à examiner (3) ou de l'ensemble source (5) ù détecteur (9) autour de l'axe de rotation (Q). 6. Procédé d'imagerie à rayons X selon l'une des revendications 1 à 4, caractérisé en ce que la rotation d'angle (A) définie à l'étape (1) correspond à deux positions angulaires (Ba, Ob) occupées simultanément par deux ensembles sources (5) - détecteurs (9) par rapport au corps à examiner (3). (8)7. Ordinateur (13) programmé pour la mise en oeuvre d'un procédé selon l'une des revendications 1 à 6. 8. Programme d'ordinateur pour la mise en oeuvre d'un procédé selon l'une des revendications 1 à 6 lorsqu'il est chargé dans un ordinateur (13).