FR 2918204 A1 20090102 FR 0704576 A 20070626 PROCEDE ET DISPOSITIF D'IMAGERIE X OU INFRAROUGE A ACQUISITION ET TRAITEMENT DE DONNEES PAR PAIRES D'ARCS CHEVAUCHES DEUX A DEUX L'invention se rapporte à un procédé et à un dispositif d'imagerie X ou infrarouge conforme à la demande française FR 0606238 déposée par le demandeur. Rappelons que la scannographie (ou tomodensitométrie), a été découverte par un ingénieur de la firme EMI, G. N. Hounsfield, en 1968 (confère les brevets US 3 924 131 ou US 3 919 552). Le brevet de 1972 s'intitule : "A method and apparatus for examination of a body by radiation such as X or gamma-radiation". Cette invention valut en 1979 le prix NOBEL à son inventeur. Le principe est le suivant : Un faisceau de rayons X, balaye un plan défini, il traverse un organe de manière linéaire et frappe une plaque ou un détecteur radiographique. La traversée de l'organe provoque une atténuation du faisceau, dont la mesure peut être effectuée grâce au détecteur. Le balayage de manière croisée, dans le plan de coupe, produit une série d'informations traitée par des logiciels appropriés sur un ordinateur associé. En effet suivant chaque axe de balayage dans un milieu hétérogène l'atténuation peut s'exprimer par une loi exponentielle, tenant compte de l'absorption photoélectrique et de la diffusion par effet Compton. Soit 10, la valeur de référence x la valeur en un point X, On peut écrire la relation suivante : f A (x) dx = F I = Io e-F In=IoE 1 t Ln = b = i A (x) dx Dont on déduit par discrétisation : ft Ln ù = JO A (x) dx = A,X1 + A2X2 + +AnXn In Les valeurs successives, A,, A2, An correspondent aux valeurs de chaque segment définis par X,, X2, >C,,. On peut alors exprimer par une série d'équations, les profils de chaque balayage associé à un angle (ou à une position) précis. On peut définir une échelle particulière par la valeur, relative à une valeur de référence du coefficient d'atténuation, celle de l'eau par exemple, ou toute autre molécule convenablement choisie. L'échelle, utilisée le plus souvent est celle relative à une molécule, abondante dans tous les organismes vivants, l'eau. Si l'on appelle A (H20), le coefficient d'atténuation de l'eau, on peut utiliser une échelle relative telle que : Bn = [An ù A(H20)] * 1000 / A(H20) La valeur du coefficient de l'eau peut être définie comme égale à 1 ou à 0, créant ainsi un système de notation facile à utiliser puisque l'eau est un composant essentiel du corps humain. D'autres systèmes peuvent d'ailleurs être utilisés, en fonction de la manière dont l'information obtenue sera exprimée (de manière visuelle). On choisit souvent la valeur de 1000 pour l'os et de -1000 pour l'air. Le traitement informatique d'un nombre suffisant de balayages croisés, définissant en fait des petites cellules ou zones élémentaires, permet la résolution d'un ensemble d'équations linéaires à condition que le nombre de balayages soit égal aux nombre de cellules. 2 L'édition et l'utilisation de l'information sont effectuées par un ordinateur associé. L'ordinateur collecte l'ensemble des données et calcule ainsi la valeur du coefficient d'atténuation de chaque zone élémentaire. L'information résultant des calculs est traduite par une carte du plan de coupe tomographique. L'ensemble des cartes constitue l'image scanner en trois dimensions de l'analyse, ce qui autorise des coupes longitudinales ou transversales. L'interprétation médicale s'appuie ainsi sur une véritable image intérieure des tissus. De telles images permettent de vérifier l'état de certains os, aussi bien que l'état du cerveau, pour détecter tumeur ou autre anomalie. Les explorations sont précédées ou complétées par d'autres explorations par exemple échographie ultrasonore ou imagerie par résonance magnétique. Le scanner et les méthodes qu'il a initiées restent un outil essentiel de l'exploration médicale. A l'origine on procédait à une série de déplacements angulaires de l'ordre de 3 degrés, répétés une centaine de fois. Les perfectionnements apportés depuis permettent d'associer plusieurs faisceaux à des barrettes de détection d'une longueur suffisante pour multiplier le nombre de mesures effectuées au même instant grâce à des détecteurs multiples. Dans les scanners de cinquième génération on utilise des barrettes de détecteur perpendiculaires aux plans de coupe pour éviter tout déplacement. L'image obtenue résulte d'un processus par étapes: Obtention des valeurs des atténuations pour chaque projection, Calculs des valeurs d'un profil, Représentation matricielle de chaque plan de coupe, Traduction de chaque représentation par une carte particulière, Etablissement d'un système cartographique spatialisé. La demande française FR 0606238 divulgue plus particulièrement un procédé d'imagerie X ou infrarouge d'un corps dans lequel, un corps à examiner étant reçu par un support, - on irradie ou illumine le corps à examiner à l'aide d'une source émettant un faisceau de rayons X ou lumineux suivant une direction de propagation, - on détecte une intensité atténuée en fonction d'une traversée des rayons X ou lumineux à travers le corps à examiner, à l'aide d'un détecteur irradié ou illuminé par le faisceau, 3 -on fait tourner d'un angle de rotation le support monté mobile autour d'un axe de rotation par rapport à la source et au détecteur montés sur un bâti ou tourner d'un angle de rotation la source et le détecteur montés sur un bâti mobile autour d'un axe de rotation par rapport au support, on convertit les intensités détectées en données pour déterminer une atténuation par le corps à examiner des rayons X ou lumineux dans un plan de coupe perpendiculaire à l'axe de rotation, à l'aide d'un convertisseur analogiqueûnumérique, et - à l'aide d'un ordinateur dûment programmé, on effectue les étapes 10 suivantes : (1) construire un vecteur générateur colonne et un vecteur générateur ligne dont les termes sont constitués par les données acquises sous deux angles de rotation, de préférence différents de 90 degrés, du support ou du bâti mobile autour de l'axe de rotation, 15 (2) construire une matrice initiale (n,m) avec les termes des deux vecteurs générateurs, en affectant à chaque zone élémentaire un terme de ligne et de colonne (Bij) représentant un coefficient d'atténuation et défini par la demi-somme, du terme homologue (i) du vecteur générateur colonne divisé par le nombre (m) de termes du vecteur générateur ligne et du terme homologue (j) du vecteur 20 générateur ligne divisé par le nombre (n) de termes du vecteur générateur colonne, et (3) ajuster le coefficient d'atténuation en chaque zone élémentaire par une méthode des moindres carrés en tenant compte, des valeurs de bordure de ligne définies par la somme des termes (Bij) sur chacune des lignes 25 de la matrice initiale et des valeurs de bordure de colonne définies par la somme des termes (Bij) sur chacune des colonnes de la matrice initiale, des termes des vecteurs générateurs comme contraintes de ligne ou de colonne et en utilisant la formule suivante : * (pj ù I,Bij)+ 1 *(ci ù lm Bij) - m J-' nm) (\ n~ Cij = Bij + 30 (~m 1 p j û ) 4 où, dans cette formule, Cij = la valeur recherchée Bij = la valeur estimée initialement (n) = le nombre de lignes de la matrice initiale (m) = le nombre de colonnes de la matrice initiale n Cij = pj pour toutes les valeurs de i, la contrainte de la colonne j a=1 m Cij = ci pour toutes les valeurs de j, la contrainte de la ligne i, j=l. pour aboutir à une matrice ajustée pour laquelle les valeurs de bordure de ligne et de colonne calculées à l'aide des valeurs ajustées (Cij) sont égales, pour 10 chaque ligne et pour chaque colonne, respectivement aux termes des vecteurs générateurs colonne et ligne. Les terrnes des vecteurs générateurs colonne ou ligne sont constitués par les données provenant de la conversion des intensités détectées dans une bande du détecteur, respectivement pour le premier ou le deuxième angle de 15 rotation. Le plan de coupe de l'objet examiné est défini par le plan qui est à la fois perpendiculaire à l'axe de rotation et qui contient la bande du détecteur considérée. Les étapes (1) à (3) sont répétées pour des données acquises sous différentes paires d'angles de rotation. Les matrices ajustées, obtenues pour les différentes paires d'angles, sont ensuite traitées par un opérateur de rotation pour 20 toutes être superposées sur une même paire d'angles. On procède enfin à une moyenne terme à terme des matrices ajustées et superposées pour aboutir à une matrice de synthèse exprimant une image des coefficients d'atténuation du corps examiné dans le plan de coupe. Pour obtenir une image du corps à examiner dans le plan de coupe 25 comportant (m x n) points, on construit, aux étapes (1) et (2), une matrice initiale de (m x n) termes, par une estimation de chacun des termes des vecteurs générateurs colonne et ligne à partir des données obtenues à la suite seulement de deux irradiations du corps à examiner sous deux angles de rotation, par exemple décalés 5 de 90 degrés. A l'étape (3), on procède à un ajustement des termes estimés de la matrice initiale par rapport aux n et m valeurs moyennes d'intensité détectées en tenant compte des valeurs de bordure. On aboutit ainsi à une première matrice ajustée correspondant à une première image ajustée. La description qui va suivre expose la méthode de calcul d'ajustement proprement dite. Cette méthode joue un rôle important dans le traitement des signaux résultant de la mesure, par les détecteurs radiographiques, de l'intensité ou de la valeur résiduelle du faisceau élémentaire produit par le dispositif à rayons X, après le parcours dans l'organisme à étudier. Si l'on veut traiter une matrice de dimensions n lignes et m colonnes si on appelle BU la valeur estimée à la ligne i et à la colonne j, si on appelle la valeur la plus probable du terme correspondant de la matrice, si on appelle pj la somme des termes de la colonne j, si on appelle ci la somme des termes de la ligne i. L'estimation de Bij résulte d'une méthode d'ajustement linéaire ou polynomial, ou plus simplement d'une moyenne arithmétique ou géométrique. Dans le cas présent, on recherchera la solution des valeurs de Ci/ en tenant compte des contraintes de lignes et de colonnes, c'est à dire le minimum de la fonction : I (Cij ù BU)2 Pour toutes valeurs de i et de j sous les contraintes : 1Cij = pj pour toutes les valeurs de j 1 Cij = ci pour toutes les valeurs de i La recherche d'un minimum de la fonction sous contraintes sera faite en utilisant la méthode des multiplicateurs de Lagrange, le Lagrangien s'écrira : 6 L = - l ij Cj ù Bij)2 + 1 Xj (~i (Cil - pj))+ ~n 1 i (Cij Cette fonction est composée de deux parties, une première n'a pas un caractère gauche et la deuxième est un ensemble de relations linéaires. Le Lagrangien est donc dérivable pour les variables Cij et Âj et ,ui, multiplicateurs de Lagrange associés aux contraintes de lignes et colonnes (nous disposons en effet de deux groupes de contraintes, les contraintes de ligne et les contraintes de colonne). Dans ces conditions nous sommes en mesure d'obtenir un ensemble de relations linéaires concernant les Cij par dérivation du Lagrangien et un ensemble de valeur de relations relatives aux valeurs de contraintes, ce qui s'écrit : En précisant que le dLldCij désigne une dérivée partielle, de la fonction L pour la variable Cij. dL=-2(Cij ùBij)++ i=0 1 dCij et les contraintes Cij = pj , pour tout j Cij = ci , pour tout i 1 a Cij = Bij+ (Xj + 0/2 L'ensemble des n*m relations correspondant aux dérivées partielles plus les n+m relations de contraintes est linéaire et n'admet qu'une solution correspondant aux nm+n+m variables. Si par exemple on veut traiter une matrice où n, le nombre de lignes, est 25 égal à 25 et m, le nombre de colonnes, est égal à 30, la solution par l'algèbre linéaire consiste à traiter : n i=1 12 Em j=1 7 0 variables Cij 25 variables correspondant aux multiplicateurs de lignes, les p i 30 variables correspondant aux multiplicateurs de colonnes, les j Un prernier objectif est déjà atteint puisque seulement 55 profils doivent 5 être établis au lieu de 750. Nous disposons au total de 750 relations correspondant aux dérivés partielles et de 55 relations correspondant aux contraintes, pour 805 variables. La résolution de ce problème en faisant appel au calcul matriciel est la solution la plus évidente mais elle implique des calculs très lourds, légèrement plus lourds que ceux 10 impliqués par les méthodes classiques. L'auteur a d'abord tablé sur l'amélioration rapide des processus de calculs, mais au delà de ce qui était son objectif essentiel, la limitation de taux d'irradiation en cours d'examen, il a poursuivi sa recherche pour tenter d'améliorer le temps de calcul. De la combinaison des relations 1 et 2 on déduit : 15 Bij + n *X,j+(El ,u i/2)=Pj Em Bij + m *pi+(Em 2j/2)=cI i=l 2 i=l On peut déduire de ces relations : *(2*(Pj- E ,Ba/)- En,,ui) X11 J * (2 (CI - E~ ~ BZJ) - E~ j) 20 Dans ces conditions en reportant par exemple la valeur de Xj dans pi, on obtient : Pour tout j = ~n~ *((PjVnBi/)-1;'=.P i) (r' ln~ 8 2918204 (ci - m Blj - I 1 J * mj pj + j=1 n Si l'on précise que p- = (1/n) 1 (i=1à n) pi soit la moyenne des multiplicateurs liée à la contrainte des lignes, on est conduit aux deux relations suivantes : Pour tout j = i2\ *(pj- \n, Pour tout i pi = p- + ? * (ci -171 Bij - n * (17 1 p j û E ijBij )) i En effet : 1 * pi est égal à p -Dans ces conditions en reportant dans la relation : Cij = Bij + ( %2 ) * () on est conduit à la relation algébrique Cette formule d'ajustement permet de déduire la matrice des Cij de la matrice des Bij, par un calcul terme à terme ( ~n, ) Pour tout i pi = Cij = Bij + 1 J * (pj û EnBij) + 1 J * n m (ci û ~1Bij)- * (E;Ip j û E ijBij ) L'auteur a ainsi réussi à aboutir de manière tout à fait surprenante à un calcul de nature algébrique ne nécessitant pas le recours au calcul matriciel. La méthode algébrique autorise le traitement partiel de la matrice de référence qui dans beaucoup de cas peut suffire. La validation numérique de cette formule de traitement de signaux et 20 d'établissement des valeurs de définition de l'image recherchée sur le plan médical, apparaît ci après. La méthode classique de résolution d'un système linéaire implique l'inversion de la matrice des coefficients des relations entre les variables et la 9 multiplication par cette matrice inverse, du vecteur exprimant les seconds membres des relations. a) Il existe 12 relations entre les variables résultant de l'expression des dérivées partielles du LAGRANGIEN de la forme : Cijû/2-pi/2=Bij b) II existe 4 relations de contraintes relatives aux colonnes et 3 relations de contraintes relatives aux lignes, de la forme : Em : Cl/ = PJ n C'ZJ = et 1=~ La méthode classique nécessite alors l'inversion d'une matrice d'une taille égale à n * m + n + m soit dans notre cas 19x19, dont le temps de calcul est évidemment beaucoup plus élevé. On utilise, pour l'estimation des valeurs initiales Bij, une méthode de génération à partir des vecteurs contraintes eux-mêmes, comme cela est exposé plus après. Nous donnons ci-dessous un exemple réduit de génération des valeurs initiales Bij, ne comprenant, dans un but de clarté, que 3x4 termes. Soit donc la matrice de n lignes et m colonnes dans laquelle n=3, m=4 MATRICE INITIALE 1 2 3 4 lignes C 1 21,42 20,90 19,58 19,75 81,65 78 2 22,29 21,80 20,46 20,63 85,18 85 3 23,29 22,80 21,46 21,63 89,18 93 colonnes 67,00 65,50 61,50 62,00 256 P 70 67 59 60 256 10 Dans cette matrice, les valeurs estimées sont inscrites sur les trois lignes et les quatre colonnes, les contraintes de lignes sont inscrites dans la colonne C. Les contraintes de colonnes sont inscrites dans la dernière ligne P. Les valeurs initiales Bij ont été générées à partir des valeurs contraintes de ligne et de colonne en utilisant une moyenne arithmétique. Ainsi : 21,417=(78/4+72/3)/2 ; 20,9=(78/4+67/3)/2 ; 19,58=(78/4+60/3)/2 ; 19,75=(78/4+60/3)/2 ; 22,292=(85/4+70/3)/2 ; 21,8=(85/4+67/3)2 ; 20,46=(85/4+59/3)/2 ; 20,625=(85/4+60/3)/2 ; 23,292=(93/4+70/3)/2 ; 22,8=(93/4+67/3)/2 ; 21,46=(93/4+59/3)/2 ; 21,625=(93/4+60/3)/2 ; Ensuite, les valeurs initiales Bij ont été sommées par ligne et par colonne. Les sommes obtenues ont été portées dans le tableau. Ainsi, 21,417 + 20,9 + 19,58 + 19,75 = 81,647 ; 22,292 + 21,8 + 20,46 + 20,625 = 85,177 ; 23,292 + 22,8 + 21,46 + 21,625 = 89,177 ; 21,417 + 22,292 + 23,292 = 67 ; L'application de la formule générale d'ajustement se simplifie ici puisque 20 le total des contraintes de colonne (ou de ligne), 256, est égal à la somme des termes de colonne (ou de ligne), 256, et conduit à : EQUILIBRE APRES CALCULS 1 2 3 4 E2 C A 1 21,50 20,49 17,83 18,17 78 78 0 2 23,25 22,26 19,58 19,91 78 85 0 3 25,25 24,26 21,58 21,91 93 93 0 E 2 70 67 59 60 256 P 70 67 59 60 0 0 0 0 20,9+21,8+22,8=65,5; 19,58 + 20,46 + 21,46 = 61,5 ; 19,75 + 20,625 + 21,625 = 62 ; 11 Chaque valeur a été ajustée de la façon suivante : 21,5=21,417 +(78-81,647)/4 + (70-67)/3 ; 25,25=23,29 + (93-89,18)/4 + (70-67)/3 ; 22,3=21,80 + (85-85,18)/4 + (67-65,5)/3 ; 17,8=19,58 + ((78-81,65)/4 + (59-61,5)3 ; 21,6=21,46 + (93-89,18)/4 + (59-61,5)/3 ; 19,917=20,63 + (85-85,18)/4 + (60-62)/3 ; 23,25=22,292 + (85-85,177)/4 + (70-67)/3 ; 20,5=20,9 + (78-81,65)/4 + (67-65,5)/3 ; 24,3=22,8 + (93-89,18)/4 + (67-65,5)/3 ; 19,6= 20,46 + (85-85,18)/4 + (59-61,5)/3 ; 18,167=19,75 + (78-81,65)/4 + (60-62)/3 ; 21,917=21,63 + (93-89,18)/4 + (60-62)/3 ; On peut vérifier, en tenant compte de calculs réalisés par une simple calculatrice que la valeur des termes sommés verticalement ou horizontalement 10 satisfont les contraintes, mais aussi qu'il conduisent bien aux résultats recherchés. Dans cet exemple réduit, on a utilisé une moyenne arithmétique pour calculer les valeurs initiales Bij à partir des valeurs de contrainte de ligne et de colonne. II est également possible d'utiliser d'autres moyennes et par exemple une moyenne géométrique. Dans ce cas, chaque valeur Bij est égale à la racine carrée 15 du produit du terme homologue du vecteur contrainte de ligne par le terme homologue du vecteur contrainte de colonne. Cette moyenne augmente toutefois le temps de calcul, c'est pourquoi l'on préfère utiliser la moyenne arithmétique. Dans la pratique, on est amené à générer des matrices dont le nombre de lignes est égal au nombre de colonnes, à partir de vecteurs de tailles identiques 20 et dont la somme des termes est, pour des raisons physiques, égale à une même valeur pour chacun des vecteurs. Dans ces conditions, le troisième terme de la formule générale est nul. On peut alors générer directement la matrice résultat de taille n à partir de deux vecteurs générateurs tels que le premier corresponde au vecteur vertical représentatif des contraintes de lignes et le second, au vecteur 25 horizontal représentatif des contraintes de colonnes. Les vecteurs de génération directe sont tels que chaque terme du vecteur vertical est égal à la contrainte de ligne divisée par n, moins la somme des termes du vecteur divisée par 2n2. De manière symétrique, les termes du vecteur horizontal sont égaux à la contrainte de colonnes divisée par n, moins la somme 30 des termes du vecteur divisée par 2n2. Du point de vue du calcul, on divise par deux le temps de calcul. 12 Les étapes d'estimation (1) et (2) et d'ajustement (3) qui viennent d'être exposées en détail sont exécutées par exemple pour 18 matrices initiales, prises sous 18 paires d'angles décalées les unes des autres de 20 degrés et de moyenner les 18 matrices ajustées en une matrice de synthèse correspondant à l'image de synthèse recherchée. Dans cet exemple, le corps ne reçoit au total que 18 x (m + n) doses élémentaires d'irradiation pour obtenir une image de synthèse comportant (m x n) points, moyenne de 18 images prises sous 18 paires d'angles d'irradiation. Ainsi, le procédé décrit dans la demande FR 0606238 permet de diminuer le taux d'irradiation dans la proportion (m + n)/(m x n) des doses élémentaires d'irradiation nécessaires pour obtenir sous une paire d'angles, une image comportant (m x n) points. Pour un nombre de points égal à 144 000 000, la diminution du taux d'irradiation est de 99,7% en prenant 18 images. Une telle diminution autorise l'exploration d'une aire du corps à examiner de 300mm de côté sous une définition de 0,025mm. L'obtention d'une définition de 0,025mm permet de détecter des tumeurs cancéreuses avant que ces dernières n'atteignent la taille critique de vascularisation des cellules centrales de la tumeur, soit 1 mm. Sans la diminution du taux d'irradiation rapportée précédemment, l'obtention d'une image comportant 144 000 000 points sous une définition de 0,025mm conduirait à exposer le corps à examiner à une dose létale d'irradiation. Le procédé divulgué dans la demande FR 0606238 est mis en oeuvre notamment par deux dispositifs. Dans un premier cas, le dispositif est conçu pour explorer la tête d'une personne. Le détecteur s'étend sur 6cm dans le plan de coupe, c'est à dire perpendiculairement à l'axe de rotation du support ou du bâti mobile, et s'étend sur 25cm parallèlement à l'axe de rotation. En supposant que toute la largeur du détecteur est utile, la bande du détecteur qui est contenue dans le plan de coupe s'étend également sur 6cm. Si l'on souhaite disposer de 36 images pour réaliser une exploration sous 360 degrés, il faut répéter les étapes (1) à (3) du procédé pour 36 paires d'angles, décalées deux à deux de 10 degrés. La distance entre la source et le détecteur doit ainsi être égale à 68cm. Dans le deuxième cas, le dispositif est conçu pour explorer le thorax d'une personne. La largeur utile du détecteur s'étend sur 48cm, de sorte que la bande du détecteur qui est contenue dans le plan de coupe s'étend également sur 13 cm. L'acquisition des données se fait sur un détecteur plan ou sur un détecteur mobile en translation parallèlement à l'axe de rotation, par pas successifs sur des hauteurs de 10cm. On obtient alors des images élémentaires de 48cm xl0cm et l'on peut choisir le nombre de pauses, pour obtenir 18 ou 36 images, en effectuant à chaque étape une rotation de 20 ou 10 degrés. Si l'on souhaite disposer de 36 images pour réaliser une exploration sous 360 degrés, en répétant les étapes (1) à (3) du procédé pour 36 paires d'angles, décalées deux à deux de 10 degrés, la distance entre la source et le détecteur doit être égale à 550cm. Pour diminuer l'encombrement dû à une telle distance entre la source et le détecteur, il est prévu des modules de jonction mobiles en translation dans le plan de coupe pour être déplacés à travers une fenêtre et déboucher dans une cabine étanche aux rayons X. Dans ces deux dispositifs, la largeur d'exploration du corps est déterminée par la largeur utile du détecteur, laquelle détermine la distance entre la source et le détecteur pour un nombre de paires d'angles donné, c'est à dire pour une même qualité d'image sous une même définition. Plus la largeur d'exploration est importante, plus la distance entre la source et le détecteur est également importante, entraînant une augmentation de l'encombrement du dispositif et des arrangements conformes aux modules de jonction décrits précédemment. L'un des buts de l'invention est de modifier le procédé et les dispositifs qui viennent d'être décrits pour permettre leur application à toute largeur d'exploration du corps à examiner, en conservant un même nombre de paires d'images de base, c'est à dire une même qualité d'image sous une même définition et en s'affranchissant de la nécessité, dans l'état de la technique, d'augmenter l'encombrement du dispositif lorsque la largeur d'exploration passe par exemple de la tête au thorax d'une personne. A cet effet, l'invention a pour objet un procédé conforme à celui rappelé en introduction, caractérisé en ce que : - on fait tourner le support ou le bâti mobile en une ou plusieurs rotation(s) de façon telle à irradier ou à illuminer le corps à examiner sous 360 degrés, - on convertit en données les intensités détectées en un point du détecteur à chaque instant de la rotation du support ou du bâti mobile, et 14 - à l'aide de l'ordinateur dûment programmé, on effectue l'étape (1) de la façon suivante : (la) traiter les données correspondant à un premier et un deuxième arcs représentant chacun une fraction angulaire des 360 degrés, choisie en fonction d'une largeur d'exploration du corps à examiner dans le plan de coupe contenant le point considéré du détecteur, le premier et le deuxième arcs faisant entre eux un angle, de préférence égal à 90 degrés, (1b) construire, pour la paire d'arcs définie à l'étape (la), le vecteur générateur colonne avec les données correspondant au premier arc de la paire et le vecteur générateur ligne avec les données correspondant au deuxième arc de la paire, et l'on effectue les deux étapes suivantes : (4) répéter les étapes (1) à (3) pour un nombre N de paires d'arcs, choisi de façon à traiter l'ensemble des données acquises sous les 360 degrés, les premiers arcs, comme les deuxièmes arcs, des N paires présentant deux à deux un chevauchement angulaire fonction du nombre N de paires d'arcs et de la fraction angulaire des 360 degrés représentée par chaque arc, et (5) traiter les matrices ajustées, obtenues pour les paires d'arcs définies à l'étape (4), par un opérateur de rotation pour toutes les superposer sur une même paire d'arcs (0 -90 ) puis moyenner terme à terme les matrices ajustées et superposées pour aboutir à une matrice de synthèse exprimant une image des coefficients d'atténuation du corps examiné dans le plan de coupe. Pour une distance L entre la source et le détecteur, seulement limitée par la dimension physique de l'objet à examiner de sorte à ce que cet objet s'inscrive dans un cercle de rotation de diamètre égal ou supérieur à cette distance L, le procédé selon l'invention permet d'obtenir N images correspondant à des arcs de longueur I, alors que la somme de tous les arcs, soit N x I, est supérieure à la circonférence du cercle, soit H x L. Le chevauchement des arcs énoncé par l'étape (4) est défini à partir d'un décalage entre deux arcs successifs sur une longueur égale à la circonférence du cercle rapportée au nombre de paires d'arcs, soit (H x L) / N. 15 Ainsi, pour une exploration de 24cm de large de la tête d'une personne, une distance R entre le détecteur et l'axe de rotation égale à 23cm définissant un cercle de rotation de diamètre égal à 2R, la source pouvant être disposée sur le cercle de rotation ou à l'extérieur de ce cercle, permet d'acquérir 36 images de base chacune fondée sur une paire d'arcs élémentaires de 24cm, avec un décalage des arcs élémentaires égal à 4cm et un chevauchement deux à deux égal à 20cm, soit une proportion de chevauchement de 83%. Pour 36 paires d'arcs élémentaires deux à deux jointifs, entre premiers arcs comme entre deuxièmes arcs de deux paires successives, la distance R entre le détecteur et l'axe de rotation devrait être égale à 137,5cm, soit un cercle de rotation de 275cm de diamètre. De même, pour une exploration du thorax d'une personne sur une largeur de 48cm, une distance R entre le détecteur et l'axe de rotation égale à 52cm définissant un cercle de rotation de diamètre égal à 2R, la source pouvant être disposée sur le cercle de rotation ou à l'extérieur de ce cercle, permet d'acquérir 36 images de base chacune fondée sur une paire d'arcs élémentaires de 48cm, avec un décalage des arcs élémentaires égal à 9,1 cm et un chevauchement deux à deux égal à 38,9cm, soit une proportion de chevauchement de 81%. Pour 36 paires d'arcs élémentaires deux à deux jointifs, entre premiers arcs comme entre deuxièmes arcs de deux paires successives, la distance Rentre le détecteur et l'axe de rotation devrait être égale à 275cm, soit un cercle de rotation de 550cm de diamètre. Autrement dit, le chevauchement des arcs permet, pour une même largeur d'exploration et un même nombre d'images prises à l'issue de l'étape (4), soit encore un même nombre N de paires d'arcs, de diminuer la distance entre la source et le détecteur par rapport à ce qu'elle serait si les arcs étaient simplement jointifs, et par conséquent, d'utiliser une source de moindre puissance. D'où il résulte que le procédé selon l'invention permet comparativement de diminuer le taux d'irradiation du corps à examiner. Dans la mesure où la largeur d'exploration n'est plus déterminée par la largeur utile du détecteur, il est suffisant d'utiliser un détecteur ponctuel. Par comparaison avec l'utilisation d'un détecteur de plusieurs centimètres ou plusieurs dizaines de centimètres de large, l'utilisation d'un détecteur ponctuel permet d'annuler en pratique tout temps de pause entre deux prises d'images successives. On peut ainsi faire tourner le support ou le bâti à vitesse constante autour de l'axe 16 de rotation, ce qui là encore contribue à diminuer le taux d'irradiation du corps à examiner. L'irradiation ou l'illumination du corps à examiner sous 360 degrés est réalisée en une ou plusieurs rotations. Avec un unique couple, l'exposition en plusieurs rotations permet de marquer un arrêt entre deux rotations. On peut ainsi demander à une personne dont on doit explorer la cage thoracique, de retenir son souffle pendant chaque rotation et de reprendre sa respiration entre deux rotations. Avantageusement, on utilise deux ou trois couples source ù détecteur, disposés par exemple à 90 ou à 120 degrés respectivement. Dans ce cas, l'acquisition et le traitement des données conformément au procédé selon l'invention sont effectuées de la façon suivante : - on irradie ou illumine le corps à examiner sous 360 degrés à l'aide de plusieurs couples source - détecteur, de préférence régulièrement décalés angulairement d'un couple à l'autre autour de l'axe de rotation, - on convertit les intensités détectées en un point de chacun des détecteurs, à chaque instant de la rotation, du support ou du bâti mobile supportant le ou les couples source - détecteur, et - on traite l'ensemble des données acquises sous 360 degrés selon les étapes (1) à (5) du procédé pour obtenir une image du corps à examiner dans le plan de coupe perpendiculaire à l'axe de rotation et contenant le point considéré du détecteur de chacun des couples source - détecteur. Avec deux ou trois couples fonctionnant simultanément, on divise par deux ou par trois, le temps d'irradiation ou d'illumination sous 360 degrés. Là encore, on contribue à diminuer le taux d'irradiation. De surcroît, on diminue également le temps d'acquisition des données par le fonctionnement simultané des couples source -détecteur. L'invention s'étend à un dispositif spécialement conçu pour la mise en oeuvre d'un procédé conforme à l'exposé qui précède et qui comprend : - un support pour recevoir un corps à examiner, 17 - une source émettant un faisceau de rayons X ou lumineux suivant une direction de propagation pour irradier ou illuminer le corps à examiner, - un détecteur irradié ou illuminé par le faisceau pour détecter une intensité atténuée en fonction d'une traversée des rayons X ou lumineux à travers le corps à examiner, un axe de rotation par rapport auquel, le support est monté mobile pour tourner d'un angle de rotation par rapport à la source et au détecteur montés sur un bâti ou par rapport auquel la source et le détecteur sont montés sur un bâti mobile pour tourner autour d'un angle de rotation par rapport au support, - un convertisseur analogiqueûnumérique pour convertir les intensités détectées en données pour déterminer une atténuation par le corps à examiner des rayons X ou lumineux dans un plan de coupe perpendiculaire à l'axe de rotation, et - un ordinateur dûment programmé pour effectuer les étapes suivantes : (1) construire un vecteur générateur colonne et un vecteur générateur ligne dont les termes sont constitués par les données acquises sous deux angles de rotation, de préférence différents de 90 degrés, du support ou du bâti mobile autour de l'axe de rotation, (2) construire une matrice initiale (n,m) avec les termes des deux vecteurs générateurs, en affectant à chaque zone élémentaire un terme de ligne et de colonne (Bij) représentant un coefficient d'atténuation et défini par la demi-somme, du terme homologue (i) du vecteur générateur colonne divisé par le nombre (m) de termes du vecteur générateur ligne et du terme homologue (j) du vecteur générateur ligne divisé par le nombre (n) de termes du vecteur générateur colonne, (3) ajuster le coefficient d'atténuation en chaque zone élémentaire par une méthode des moindres carrés en tenant compte, des valeurs de bordure de ligne définies par la somme des termes (Bij) sur chacune des lignes de la matrice initiale et des valeurs de bordure de colonne définies par la somme des termes (Bij) sur chacune des colonnes de la matrice initiale, des termes des vecteurs générateurs comme contraintes de ligne ou de colonne et en utilisant la formule suivante : 18 / 1 * (ci û lm Bij) - (1) * m) '-' nm Cij = Bij + n, * (pi - 1,Bj)+ ( EUBij ) où, dans cette formule, Cij = la valeur recherchée Bij = la valeur estimée initialement (n) = le nombre de lignes de la matrice initiale (m) = le nombre de colonnes de la matrice initiale n i -1 Cij = pj pour toutes les valeurs de i, la contrainte de la colonne j m Cij = ci pour toutes les valeurs de j, la contrainte de la ligne i, j=1 pour aboutir à une matrice ajustée pour laquelle les valeurs de bordure de ligne et de colonne calculées à l'aide des valeurs ajustées (Cij) sont égales, pour chaque ligne et pour chaque colonne, respectivement aux termes des vecteurs générateurs ligne et colonne, caractérisé en ce que, - le support ou le bâti est mobile de façon telle à irradier ou illuminer le corps sous 360 degrés en une ou plusieurs rotation(s), - le convertisseur convertit en données les intensités détectées en un point du détecteur à chaque instant de la rotation du support ou du bâti mobile, et - l'ordinateur est dûment programmé pour, effectuer l'étape (1) de la façon suivante : (la) traiter les données correspondant à un premier et un deuxième arcs représentant chacun une fraction angulaire des 360 degrés, choisie en fonction d'une largeur d'exploration du corps à examiner dans le plan de coupe 19 contenant le point considéré du détecteur, le premier et le deuxième arcs faisant entre eux un angle., de préférence égal à 90 degrés, (1b) construire, pour la paire d'arcs définie à l'étape (la), le vecteur générateur colonne avec les données correspondant au premier arc de la paire et le vecteur générateur ligne avec les données correspondant au deuxième arc de la paire, et pour effectuer les deux étapes suivantes : (4) répéter les étapes (1) à (3) pour un nombre N de paires d'arcs, choisi de façon à traiter l'ensemble des données acquises sous les 360 degrés, les premiers arcs, comme les deuxièmes arcs, des N paires présentant deux à deux un chevauchement angulaire fonction du nombre N de paires d'arcs et de la fraction angulaire des 360 degrés représentée par chaque arc, et (5) traiter les matrices ajustées, obtenues pour les paires d'arcs définies à l'étape (4), par un opérateur de rotation pour toutes les superposer sur une même paire d'arcs (0 -90 ) puis moyenner terme à terme les matrices ajustées et superposées pour aboutir à une matrice de synthèse exprimant une image des coefficients d'atténuation du corps examiné (9) dans le plan de coupe. Comme indiqué précédemment, le chevauchement des arcs permet, pour une même largeur d'exploration et un même nombre d'images prises à l'issue de l'étape (4), soit encore un même nombre N de paires d'arcs, de diminuer la distance entre la source et le détecteur par rapport à ce qu'elle serait si les arcs étaient simplement jointifs. D'où il résulte que le dispositif selon l'invention présente comparativement un encombrement plus faible pour des largeurs d'exploration plus grandes que la largeur correspondant à des arcs jointifs, à même nombre de prises d'image. Avantageusement, le détecteur comprend une ligne de points disposés suivant un axe parallèle à l'axe de rotation. Plus particulièrement, le détecteur est une barrette disposée parallèlement à l'axe de rotation. Cet agencement permet d'explorer le corps à examiner suivant toute une longueur d'exploration, parallèlement à l'axe de rotation. En traitant les données acquises par chacun des points de la barrette de détection, on obtient une succession d'images du corps à examiner suivant une succession de plans de coupe contenant chacun un point 20 considéré de la barrette. Il est ensuite permis de construire des images en trois dimensions du corps à examiner. Un détecteur ponctuel ou ayant la forme d'une barrette disposée parallèlement à l'axe de rotation présente avantageusement une inertie en rotation comparativement plus faible qu'une plaque de détection. De surcroît, l'utilisation d'une source comparativement moins puissante permet également de diminuer son inertie en rotation. Par conséquent, le bâti supportant la source et le détecteur peut plus aisément être mis en rotation autour de l'axe de rotation et peut tourner à une vitesse de rotation comparativement plus élevée. Le dispositif selon l'invention permet ainsi de diminuer le temps d'acquisition imposé par la rotation du détecteur et de la source autour de l'axe de rotation sous 360 degrés. D'où il résulte là encore que le dispositif selon l'invention contribue à diminuer le taux d'irradiation. 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 un tableau rapportant une simulation de données acquises sous 360 degrés, degré par degré, dans un procédé selon l'invention. La figure 2 illustre de façon schématique les étapes (la), (lb) et (4) du procédé selon l'invention, limitées aux 3 première paires d'arcs présentant deux à deux un chevauchement. La figure 3 est un tableau rapportant les données du tableau de la figure 1, regroupées en sous-ensembles se chevauchant en correspondant aux premiers et aux deuxièmes arcs des 6 premières paires d'arcs, les 3 premières étant illustrées par la figure 2. La figure 4 montre de façon schématique un premier mode d'exécution d'un dispositif selon l'invention, dans lequel un couple source û détecteur est mobile autour d'un axe de rotation vertical. La figure 5 montre de façon schématique un deuxième mode d'exécution d'un dispositif selon l'invention, dans lequel deux couples source û détecteur sont mobiles autour d'un axe de rotation horizontal. 21 La figure 6 montre de façon schématique un troisième mode d'exécution de l'invention d'un dispositif selon l'invention dans lequel le support du corps à examiner est mobile autour d'un axe de rotation vertical. La figure 7 montre une variante du troisième mode d'exécution de l'invention, à quatre barrettes de détection alignées verticalement et irradiées par deux sources. La figure 8 montre une autre variante du troisième mode d'exécution d'un dispositif selon l'invention, prévu pour fonctionner à l'air libre. La figure 9 montre dans un plan de coupe horizontal, un quatrième mode d'exécution d'un dispositif selon l'invention, à trois couples source û détecteur. Pour l'exposé d'un mode d'exécution détaillé d'un procédé selon l'invention, on se contentera tout d'abord de faire un renvoi intégral à ce qui a été décrit précédemment, concernant tant les principes de la scannographie que le traitement des données par les étapes (1) et (2) d'estimation et par l'étape (3) d'ajustement. L'invention tient plus spécifiquement en ce que : - on fait tourner le support ou le bâti mobile en une ou plusieurs rotation(s) de façon telle à irradier ou à illuminer le corps à examiner sous 360 degrés, - on convertit en données les intensités détectées en un point du détecteur à chaque instant de la rotation du support ou du bâti mobile, et - à l'aide de l'ordinateur dûment programmé, on effectue l'étape (1) de la façon suivante : (la) traiter les données correspondant à un premier et un deuxième arcs représentant chacun une fraction angulaire des 360 degrés, choisie en fonction d'une largeur d'exploration du corps à examiner dans le plan de coupe contenant le point considéré du détecteur, le premier et le deuxième arcs faisant entre eux un angle, de préférence égal à 90 degrés, 22 (1 b) construire, pour la paire d'arcs définie à l'étape (1 a), le vecteur générateur colonne avec les données correspondant au premier arc de la paire et le vecteur générateur ligne avec les données correspondant au deuxième arc de la paire, et l'on effectue les deux étapes suivantes : (4) répéter les étapes (1) à (3) pour un nombre N de paires d'arcs, choisi de façon à traiter l'ensemble des données acquises sous les 360 degrés, les premiers arcs, comme les deuxièmes arcs, des N paires présentant deux à deux un chevauchement angulaire fonction du nombre N de paires d'arcs et de la fraction angulaire des 360 degrés représentée par chaque arc, et (5) traiter les matrices ajustées, obtenues pour les paires d'arcs définies à l'étape (4), par un opérateur de rotation pour toutes les superposer sur une même paire d'arcs (0 -90 ) puis moyenner terme à terme les matrices ajustées et superposées pour aboutir à une matrice de synthèse exprimant une image des coefficients d'atténuation du corps examiné dans le plan de coupe. L'acquisition des donnée s'effectue sous 360 . Une simulation est présentée ci-après et illustrée à l'aide des figures 1 à 3. L'utilisateur, par le procédé selon l'invention, peut augmenter le nombre d'images à traiter, par exemple obtenir 72 images au lieu de 36, avec les mêmes données. Dans ce nouvel algorithme, nous avons donc la possibilité d'effectuer un mode d'acquisition et deux modes de découpage des données (36 ou 72 dans l'exemple considéré). Si l'on veut travailler avec le même taux d'erreur, le rayon de rotation peut être divisé par deux et atteindre 76,4 cm. A titre de comparaison, si nous prenons par exemple une définition de 81 microns, on pourra choisir pour obtenir 20 images par faisceau et une largeur de 48cm, avec le procédé à paires d'arcs jointifs selon l'art antérieur exposé dans la demande FR 0606238, un rayon de 152,8 cm environ, mesurant la distance entre l'axe de rotation et le détecteur. Dans ces conditions, le temps d'acquisition sera de l'ordre de 118 secondes, pour une définition de 81 microns. Pour les parties moins importantes de l'organisme, par exemple la tête, on pourra, également selon le procédé exposé dans la demande 23 FR 0606238, se contenter d'images de 24cm et obtenir 40 images pour chaque faisceau.. Le chevauchement deux à deux des arcs, dans le procédé selon l'invention, procède d'un découpage des données acquises sous 360 degrés. II peut se faire par exemple par groupes de 10 données pour chaque arc, les groupes se chevauchant de sorte à augmenter le nombre d'images de référence et améliorer ainsi la qualité de l'information, ou plus simplement, diminuer le rayon R du cercle de rotation en autorisant ce chevauchement d'images. Les groupes ou sous-ensembles chevauchés seront ainsi constitués : 15 groupe Al : 1 à 10, groupe A2 : 6 à 15, groupe A3 : 11 à 20, groupe A4 : 16 à 25, groupe A5 : 21 à 30, groupe A6 : 26 à 35, groupe B1 : 91 à 100 groupe B2 : 96 à 105 groupe B3 : 101 à 110 groupe B4 : 106 à 115 groupe B5 : 111 à 120 groupe B6 : 116 à 125 Nous rappelons que l'on génère une image ajustée à partir de deux vecteurs orthogonaux pris dans les bitmaps homologues des images initiales. Il est donc possible de réaliser dans un plan 36 images ajustées dont on tirera une image moyenne qui constituera l'image tomodensitométrique. Nous aurons ainsi pour un 20 seul tour la possibilité de réaliser 72 images et l'on pourra soit améliorer la qualité de l'information pour le même rayon de rotation ; soit diminuer le rayon de moitié pour obtenir 36 images de dimension équivalente. La figure 2 illustre de façon schématique les étapes (la), (1b) et (4) du procédé à partir de 3 paires d'arcs présentant deux à deux un chevauchement. Plus 25 précisément, on a représenté 3 premières paires d'arcs élémentaires numérotées (Al,B1), (A2,B2) et (A3,B3). Dans chaque paire, les arcs (A) et (B) sont orthogonaux. Bien que la représentation ne soit pas à l'échelle, les arcs élémentaires (A) et (B) correspondent à un angle de 10 degrés, soit une fraction angulaire des 360 degrés égale à 1/36. Les premiers arcs (A), comme les 30 deuxièmes arcs (B), présentent deux à deux, c'est à dire d'une paire à la paire 24 suivante, un décalage correspondant à un angle de 5 degrés. Le nombre total N de paires est égal à 360/5 = 72. Le tableau de la figure 3 rapporte les données correspondant aux premiers et aux deuxièmes arcs de chacune des 3 premières paires d'arcs illustrées par le figure 2 et des trois suivantes (A4,B4), (A5,B5) et (A6,B6). Les données sont supposées être acquises pour toute position angulaire du détecteur, de degré en degré. Seulement les données acquises sous les 125 premiers degrés sont considérées. Ces données sont regroupées en sous-ensembles ayant chacun les données acquises sous 10 degrés consécutifs et chaque sous-ensemble de données correspond à un arc élémentaire représentant angulairement une fraction égale à 1/36 des 360 degrés. Les sousensembles se chevauchent en étant deux à deux décalés de 5 degrés consécutifs. Pour construire la première paires d'arcs élémentaires (A1,B1) on considère le sous-ensemble des données acquises entre let 10 degrés pour l'arc Al et le sous-ensemble des données acquises entre 91 et 100 degrés pour l'arc B1, orthogonal à l'arc Al. Pour construire la deuxième paire d'arcs élémentaires (A2,B2), on considère le sous-ensemble des données acquises entre 6 et 15 degrés pour l'arc A2 et le sous-ensemble des données acquises entre 96 et 105 degrés pour l'arc B2, orthogonal à l'arc A2. On procède de même pour les paires d'arcs (A3,B3) à (A6,B6) et jusqu'à la paire d'arcs (A72,B72). Pour chaque paire d'arcs (An,Bn), on construit le vecteur générateur colonne avec les données du sous-ensemble correspondant au premier arc (An) de la paire et le vecteur générateur ligne avec les données du sous-ensemble correspondant au deuxième arc (Bn) de la paire. Dans le tableau 2, les données acquises sont simulées aléatoirement entre 1 et 256. Les figures 2 et 3 dans l'illustration du procédé selon l'invention. Si l'on souhaite effectuer une image d'une région de 6 cm de large, de la tête, avec 72 images de base, on peut se contenter d'un cercle de rotation C dont le rayon est égal à 34,4 cm et d'une distance entre le foyer du tube à rayons X et le détecteur à peine supérieure à 68,8 cm, ce qui suppose un chevauchement de 5 degrés pour 2 arcs élémentaires de 10 degrés. Par comparaison, des arcs de 5 degrés seulement jointifs imposeraient une distance entre la source et le détecteur égale ou supérieure à 137,5 cm. De même si l'on souhaite effectuer une image élémentaire d'une région de 48 cm de large, par exemple du thorax, et 72 images de base 25 jointives, on serait conduit à un cercle de rotation de 11 m et à une distance entre la source et le détecteur égale ou supérieure à 11 m. Par contre si l'on accepte un chevauchement égale à 5 degrés pour des arcs de 10 degrés, ce diamètre peut être divisé par 2, ce qui conduit à une distance de 5,5 m entre le foyer du tube à rayons X et le détecteur. Les figures 4 et 5 illustrent de façon schématique un premier et un deuxième mode d'exécution d'un dispositif selon l'invention. Figure 4, un dispositif comprend : - un support 1 pour recevoir un corps à examiner 3, - une source 5 émettant un faisceau 7 de rayons X ou lumineux suivant une direction de propagation pour irradier ou illuminer le corps à examiner 3, la source 5 est fixée à un plaque de support 9 en plomb comportant une fente verticale 10, c'est à dire parallèle à l'axe de rotation, de faible hauteur pour limiter la quantité de rayons X à ce qui est strictement utile à l'irradiation d'un détecteur 11, - un détecteur 11 irradié ou illuminé par le faisceau 7 pour détecter une intensité atténuée en fonction d'une traversée des rayons X ou lumineux à travers le corps à examiner 3, le détecteur 11 étant fixé à une plaque de support 13 en plomb pour protéger l'électronique 12 embarquée avec le détecteur 11, la plaque de support 13 comportant une fente ouvrant sur le détecteur 11, - une potence verticale 15 permettant de supporter une console 17 mobile en rotation autour d'un axe de rotation S2, à laquelle le détecteur 11 et la source à rayons X 5 sont fixés, - la source 5 et le détecteur 11 étant montés sur la console 17 mobile pour tourner autour d'un angle de rotation par rapport au support 1, la console 17 étant mobile autour de l'axe de rotation SZ de façon telle à irradier ou illuminer le corps sous 360 degrés en une ou plusieurs rotation(s), un convertisseur analogiqueûnumérique 19, pour convertir les intensités détectées en données pour déterminer une atténuation par le corps à examiner des rayons X ou lumineux dans un plan de coupe perpendiculaire à l'axe de rotation Q, le convertisseur convertissant en données les intensités détectées en un point du détecteur 11 à chaque instant de la rotation de la console mobile 17, et 26 - un ordinateur 21 dûment programmé pour effectuer (la) à (5) du procédé selon l'invention. Le détecteur 11 a la forme d'une barrette disposée verticalement, c'est-à-dire parallèlement à l'axe de rotation SZ de la console mobile 17. On fait monter la personne sur le support 1 puis on fait tourner la console 17 de sorte à ce que la vitesse de rotation au niveau du détecteur 11 soit par exemple de 2,70cm par seconde, pour 27 microns de définition ou 8,10 cm par seconde, pour 81 microns de définition. Il est prévu dans ce mode d'exécution de disposer l'ensemble dans une cabine de petite dirnension, faite de parois en plomb 23 et comportant une porte en verre au plomb 25. Cette cabine de protection a une forme parallélépipédique, par exemple de hauteur 2 m et de largeur et profondeur 1 m. Compte-tenu que la source 5 utilisée est comparativement de plus faible puissance, les parois 23 comprennent avantageusement constituées des tissus au plomb. La figure 5 montre d'une manière classique un deuxième mode d'exécution d'un dispositif selon l'invention. Les mêmes références identiques les mêmes éléments décrits en rapport avec la figure 3. Une personne est étendue sur un lit 1 parallèlement à l'axe de rotation S2 horizontal et déplaçable horizontalement. Le tube à rayons X de la source 5 et le détecteur 11 sont montés sur une structure circulaire 17 que l'on fait tourner à une vitesse appropriée à la prise de vue à l'aide de moteurs 27 montés sur des supports 29. Dans ce cas là aussi, le détecteur 11 a la forme d'une barrette de détection et est placé parallèlement à l'axe de rotation û. Le dispositif comprend deux détecteurs 11 irradiés par deux sources à rayons X 5, disposés de façon orthogonale. Cet agencement permet de réaliser les 360 degrés d'exploration du corps à examiner en seulement un demi-tour du bâti mobile 17 au lieu du tour complet nécessaire lorsque l'on utilise un seul couple source û détecteur. Les figures 6 et 7 montrent un troisième mode d'exécution d'un dispositif selon l'invention, dans lequel l'axe de rotation S2 est vertical et la personne est reçue debout sur le support 1 mobile en rotation autour de l'axe de rotation S2.- Le patient 3 est placé sur le support rotatif mis en rotation par un moteur, non représenté, 27 commandé par une électronique destinée à contrôler la vitesse de rotation ainsi que la ou les séquences d'arrêts et démarrages de la rotation si l'on doit prévoir des arrêts pour permettre au patient de reprendre sa respiration entre deux rotations. Le patient 3 peut être tenu verticalement ou assis sur un siège transparent aux rayons X. On peut dans ce cas augmenter considérablement la hauteur de prise de vue en plaçant deux couples faisceaux - détecteurs 5, 11 à deux hauteurs différentes parallèlement à l'axe de rotation S2. Compte tenu de la dimension des détecteurs actuels, on obtiendra ainsi des prises de vue sur des hauteurs pouvant aller jusqu'à 92cm, sans laisser d'espace entre les différentes barrettes de détection 13, grâce à un dispositif approprié de connexions. les barrettes de détection, par exemple de 48 cm de hauteur, sont disposées exactement sur la même ligne verticale, parallèlement à l'axe de rotation, pour que le bas d'une barrette soit exactement au niveau du haut de la barrette suivante. Ceci peut aussi s'appliquer à un ensemble de 8 barrettes de détection verticales de 23 cm dont la hauteur totale est de l'ordre de 184 cm. Le dispositif peut être encore simplifié conformément au schéma de la figue 8. L'ensemble est à l'air libre ou protégé s'il est à l'extérieur par une simple paroi en toile plombée ou en matériau plus fin. On pourra obtenir des images de n'importe quelle partie de l'organisme en déplaçant de manière homologue le tube à rayons X 5 et le détecteur 11 grâce aux deux potences à de vérins verticaux 31. Le dispositif comprend : - une petite plate-forme rotative 1 sur laquelle le patient 3 peut se tenir debout, -une potence à vérins 31 portant un détecteur 11 inséré dans une plaque de plomb 13, laquelle protège l'électronique 12 liée au détecteur 11, et -une potence à vérins 31 portant un tube à rayons X 5 devant lequel est placée une simple fente de collimateur 10. Les deux potences à vérins 31 permettent de déplacer de manière homologue le tube à rayons X 5 et le détecteur 11. Un tel dispositif peut être transporté dans une camionnette. Le détecteur sera connecté à un système informatique local comportant un écran de grande dimension pour obtenir de 28 bonnes images. Ce système disposera si nécessaire d'une connexion sur Internet pour envoyer les images de base à un centre de calcul à distance, lequel effectuera le traitement des données à l'aide d'un ordinateur dûment programmé pour la mise en oeuvre du procédé selon l'invention, pour produire des images tomodensitométriques qui seront renvoyées au système local pour permettre leur visualisation, notamment pour des images à très haute définition. Ces images peuvent d'ailleurs être soumises, à la demande du médecin traitant, à des spécialistes connectés au réseau, dont l'avis sera particulièrement précieux pour des médecins installés dans des pays en voie de développement peu équipés. Dans le cas où il s'agirait d'un examen entraînant la production d'un grand nombre d'images tomodensitométriques, ces images pourraient être stockées sur des clefs USB de grande capacité et adressées par courrier rapide au médecin traitant, ou au médecin consultant. La mise en service du matériel d'acquisition et la connexion par Internet à un centre de calcul situé dans un pays développé bénéficierait ainsi à tous les pays du monde, en particulier aux pays en voie de développement qui pourront utiliser le dispositif simplifié dans les prochaines années. La figure 9 montre schématiquement et dans un plan de coupe horizontal, un quatrième mode d'exécution d'un dispositif selon l'invention, à trois couples source ù détecteur. Pour encore diminuer le temps d'acquisition, il est prévu de réaliser un dispositif comprenant, dans une cabine de protection de formehexagonale, trois couples 5, 11 tubes à rayons X ù détecteur de puissance limitée. Une rotation de 120 degrés suffit pour obtenir une image complète ,ce qui divise par trois le temps d'acquisition. Pour obtenir une image à une définition de 81 microns, les 3 paires de détecteur peuvent se déplacer verticalement et horizontalement pour se rapprocher de l'organe à examiner. Si l'on admet en effet, un rayon de 15 cm, pour l'observation de la tête, ce qui donne un diamètre de 30 cm, le cercle de rotation sera égal à 94,3 cm, une 30 rotation totale s'effectuerait en 11,6 secondes, mais grâce aux trois ensembles l'acquisition correspondant à un cercle complet peut s'effectuer 3,87 secondes et l'on peut obtenir ainsi en continuant à faire tourner le cercle, obtenir une image toutes les 3,87 secondes, en ajustant la définition pour que le temps de rafraîchissement de l'image soit inférieur ou égal au tiers du temps de la rotation. 25 29 Ceci peut être obtenu avec les performances actuelles de calcul par un PC de puissance suffisante. Avec le procédé et les dispositifs selon l'invention, nous avons ainsi examiné les possibilités d'utiliser des images chevauchées pour obtenir un nombre suffisant d'images, par exemple 36, en vue de générer des scans de qualité en ayant ainsi un temps d'acquisition réduit par exemple de moitié. Nous disposons actuellement d'un détecteur de 23cm, ce qui permet avec un tube à rayons X d'ouverture totale de 24 degrés, de placer l'objet à une distance de 52 cm (ou plus) pour que toute la hauteur du détecteur soit atteinte. Dans ces conditions il est possible de réaliser un appareil bénéficiant des avantages suivants : - le temps d'acquisition à 81 microns : 21 secondes, si la distance entre l'axe de rotation et le détecteur est de 27cm, ce qui laisse un espace suffisant même pour un patient corpulent, ou 71 secondes, dans les mêmes conditions mais pour une définition de 27 microns. - la quantité de rayons X reçue : particulièrement réduite. - la possibilité d'utiliser un tube à rayons X de faible puissance et donc de faible coût. Pour certaines parties de l'organisme, par exemple la tête, on peut réduire la distance entre l'axe de rotation et le détecteur par exemple à 15 cm ce qui réduit le temps d'acquisition à 11,1 secondes pour 81 microns ou à 33.3 secondes pour 27 microns, et la dose reçue de manière encore plus élevée en réduisant la puissance du faisceau par suite du rapprochement du détecteur et du tube à rayons X. Pour des positions rapprochées du tube et du détecteur, la distance au centre du détecteur est significativement inférieure à la distance aux extrémités et une correction doit être effectuée qui est fonction du carré des rapports des distances. La fonction de correction peut aisément être incluse dans le processus d'acquisition pour soumettre au calcul des matrices de valeurs corrigées, en 30 fonction de la taille du détecteur et de la distance relative du foyer du tube à rayons X au centre et du même foyer aux points sur le détecteur Les avantages découlant de l'utilisation d'images initiales chevauchées peuvent encore être mis en évidence sur l'exemple concret suivant : a) supposons un rayon de rotation de 25 cm b) la circonférence est de 158 cm (arrondi à 160 cm) c) pour obtenir 20 images dont la longueur d'arc soit de 50 cm il faudrait théoriquement une circonférence de 316 cm, un décalage de l'ordre de 8 cm est suffisant pour obtenir 20 images de 50 cm recouvrant la totalité du cercle de 160 cm avec u chevauchement pour la dernière prise de vue au delà des 160 cm Mais il est aussi possible d'obtenir 32 images avec les mêmes données de base en décidant un décalage de 5 cm au lieu de 8. On peut aller au delà car avec un décalage de 1 cm on pourrait produire 160 images distinctes et théoriquement en se limitant à une distance légèrement supérieure à 81 microns soit 100 microns on pourrait obtenir 16 000 images distinctes (ce qui est sans intérêt d'ailleurs à cause du temps de calcul prohibitif). Les procédures de base, au cours des prochaines années consisteront donc, à retenir pour l'essentiel la définition de 81 microns, ce qui donne un temps d'acquisition de 103 secondes pour l'ensemble de l'organisme et un temps de calcul d'une cinquantaine de minutes et de traiter par le procédé du zoom des zones plus petites de dimension de l'ordre de quelques dizaines de cm3 ; mais les techniques évoluant, le calculateur associé permettra d'offrir des temps se réduisant rapidement et des fonctionnalités nouvelles. Pour une définition de 243 microns le temps de calcul serait de 111 secondes. En effet le logiciel permet de calculer les images tomodensitométriques à toute définition multiple de la définition de l'image d'acquisition. L'invention permet de réaliser un scanner médical à haute définition dont le procédé d'acquisition des images de base est amélioré, ce qui a comme 31 conséquences de faciliter la réalisation d'un scanner simplifié pouvant être utilisé en particulier dans les pays en voie de développement. En effet, dans tous les pays en voie de développement ou dans les pays émergents les méthodes actuelles conduisent à des scanners médicaux dont le prix minimum est de l'ordre de 650 000{ et dont la meilleure définition est de l'ordre de 650 microns. L'auteur de l'invention a été conduit à mettre au point une technique d'acquisition et de traitement des données permettant la réalisation d'un scanner à haute définition, de l'ordre de 27 microns, mais dont le coût est largement inférieur aux scanners actuels. En résumé, il est possible d'améliorer considérablement les dispositifs déjà brevetés en recourant à une méthode de génération des images initiales par superposition partielle des images successives. Cette amélioration conduit a pouvoir rapprocher considérablement les tubes à rayons X des détecteurs ; il résulte de ce rapprochement plusieurs conséquences essentielles : 1/ la dose de rayons X reçue peut être considérablement réduite, puisqu'elle est fonction du carré de la distance entre le foyer du tube à rayons X et le détecteur. 2/ la puissance et le coût du tube à rayons X peuvent être réduits considérablement. 3/ l'acquisition peut se faire de manière continue ou pas à pas, par simple rotation du patient, ou par simple rotation de l'ensemble tube à rayon X détecteur. 4/ le temps d'acquisition étant réduit, à proportion de la distance entre le tube et le détecteur,, le nombre d'opérations réalisables dans un temps donné par un scanner peut être augmenté considérablement. 5/ la simplification du dispositif et la réduction des coûts qu'il entraîne autorise la réalisation d'un scanner simplifié, transportable facilement, utilisable 32 dans les pays en voie de développement à cause de son faible coût et de sa transportabilité. 6/ L'extension du réseau Internet touche maintenant un nombre considérable de pays, autorise la transmission des données, acquises localement, par Internet et le traitement par des centres de calcul puissants et évolutifs, dont les images tomodensitométriques produites peuvent être transmises soit par Internet, soit autrement (par exemple grâce à des clefs USB de grande capacité) au médecin traitant ou à des médecins consultants de haut niveau, susceptibles d'interpréter rapidement les résultats de l'imagerie médicale et de transmettre au médecin traitant des réponses qualifiées. 33 REVENDICATIONS 1. Procédé d'imagerie X ou infrarouge d'un corps dans lequel, un corps à examiner (9) étant reçu par un support (1), - on irradie ou illumine le corps à examiner (3) à l'aide d'une source (5) émettant un faisceau (7) de rayons X ou lumineux suivant une direction de propagation, - on détecte une intensité atténuée en fonction d'une traversée des rayons X ou lumineux à travers le corps à examiner (3), à l'aide d'un détecteur (11) irradié ou illuminé par le faisceau, - on fait tourner d'un angle de rotation le support (1) monté mobile autour d'un axe de rotation (û) par rapport à la source (5) et au détecteur (11) montés sur un bâti (23,31) ou tourner d'un angle de rotation la source (5) et le détecteur (11) montés sur un bâti (17) mobile autour d'un axe de rotation (S)) par rapport au support (1), - on convertit les intensités détectées en données pour déterminer une atténuation par le corps à examiner des rayons X ou lumineux dans un plan de coupe perpendiculaire à l'axe de rotation (û), à l'aide d'un convertisseur (19) analogiqueûnumérique, et - à l'aide d'un ordinateur (21) dûment programmé, on effectue les étapes suivantes : (1) construire un vecteur générateur colonne et un vecteur générateur ligne dont les termes sont constitués par les données acquises sous deux angles de rotation, de préférence différents de 90 degrés, du support ou du bâti mobile autour de l'axe de rotation, (2) construire une matrice initiale (n,m) avec les termes des deux vecteurs générateurs, en affectant à chaque zone élémentaire un terme de ligne et de colonne (Bij) représentant un coefficient d'atténuation et défini par la demi-somme, du terme homologue (i) du vecteur générateur colonne divisé par le nombre (m) de termes du vecteur générateur ligne et du terme homologue (j) du vecteur générateur ligne divisé par le nombre (n) de termes du vecteur générateur colonne, 34(3) ajuster le coefficient d'atténuation en chaque zone élémentaire par une méthode des moindres carrés en tenant compte, des valeurs de bordure de ligne définies par la somme des termes (Bij) sur chacune des lignes de la matrice initiale et des valeurs de bordure de colonne définies par la somme des termes (Bij) sur chacune des colonnes de la matrice initiale, des termes des vecteurs générateurs comme contraintes de ligne ou de colonne et en utilisant la formule suivante : Cij = Bij + n* (Pi ù 1 Bl>) + (m) * ( ci où, dans cette formule, Cij = la valeur recherchée Bij = la valeur estimée initialement (n) = le nombre de lignes de la matrice initiale (m) = le nombre de colonnes de la matrice initiale n i _ 1Cij = pj pour toutes les valeurs de i, la contrainte de la colonne j m Cij = ci pour toutes les valeurs de j, la contrainte de la ligne i, j=1 pour aboutir à une matrice ajustée pour laquelle les valeurs de bordure de ligne et de colonne calculées à l'aide des valeurs ajustées (Cij) sont égales, pour chaque ligne et pour chaque colonne, respectivement aux termes des vecteurs générateurs ligne et colonne, caractérisé en ce que, -on fait tourner le support (1) ou le bâti (17) mobile en une ou plusieurs rotation(s) de façon telle à irradier ou illuminer le corps (3) sous 360 degrés, - on convertit en données les intensités détectées en un point du détecteur (11) à chaque instant de la rotation du support (1) ou du bâti (17) mobile, et En" Bij) - J=l 1\ 35- à l'aide de l'ordinateur (21) dûment programmé, on effectue l'étape (1) de la façon suivante : (la) traiter les données correspondant à un premier (Al) et un deuxième (B1) arcs représentant chacun une fraction angulaire des 360 degrés, choisie en fonction d'une largeur d'exploration du corps à examiner dans le plan de coupe contenant le point considéré du détecteur (11), le premier (Al) et le deuxième arcs (B1) faisant entre eux un angle, de préférence égal à 90 degrés, (1b) construire, pour la paire d'arcs (Al,Bl) définie à l'étape (la), le vecteur générateur colonne avec les données correspondant au premier arc (Al) de la paire et le vecteur générateur ligne avec les données correspondant au deuxième arc (B1) de la paire, et l'on effectue les deux étapes suivantes : (4) répéter les étapes (1) à (3) pour un nombre N de paires d'arcs (An,Bn) choisi de façon à traiter l'ensemble des données acquises sous les 360 degrés, les premiers arcs (An), comme les deuxièmes arcs (Bn), des N paires présentant deux à deux un chevauchement angulaire fonction du nombre N de paires d'arcs et de la fraction angulaire des 360 degrés représentée par chaque arc, et (5) traiter les matrices ajustées, obtenues pour les paires d'arcs (An,Bn) définies à l'étape (4), par un opérateur de rotation pour toutes les superposer sur une même paire d'arcs (0 -90 ) puis moyenner terme à terme les matrices ajustées et superposées pour aboutir à une matrice de synthèse exprimant une image des coefficients d'atténuation du corps examiné (3) dans le plan de coupe. 2. Procédé selon la revendication 1, caractérisé en ce que, - on irradie ou illumine le corps à examiner (3) sous 360 degrés à l'aide de plusieurs couples source - détecteur (5,11), de préférence régulièrement décalés angulairement d'un couple à l'autre autour de l'axe de rotation (S)), 36-on convertit les intensités détectées en un point de chacun des détecteurs (11), à chaque instant de la rotation, du support (1) ou du bâti (17) mobile supportant le ou les couples source - détecteur (5,11), et - on traite l'ensemble des données acquises sous 360 degrés selon les étapes (1) à (5) du procédé pour obtenir une image du corps à examiner (3) dans le plan de coupe perpendiculaire à l'axe de rotation (S)) et contenant le point considéré du détecteur (11) de chacun des couples source -détecteur. 3. Dispositif d'imagerie X ou infrarouge d'un corps, spécialement conçu pour la mise en oeuvre d'un procédé selon la revendication 1 ou 2 en comprenant un support (1) pour recevoir un corps à examiner (3), - une source (5) émettant un faisceau (7) de rayons X ou lumineux suivant une direction de propagation pour irradier ou illuminer le corps à examiner (3), - un détecteur (11) irradié ou illuminé par le faisceau pour détecter une intensité atténuée en fonction d'une traversée des rayons X ou lumineux à travers le corps à examiner (3), - un axe de rotation (Q) par rapport auquel, le support (1) est monté mobile pour tourner d'un angle de rotation par rapport à la source (5) et au détecteur (11) montés sur un bâti (23,31) ou par rapport auquel la source (5) et le détecteur (11) sont montés sur un bâti (17) mobile pour tourner autour d'un angle de rotation par rapport au support (1), - un convertisseur (19) analogiqueûnumérique pour convertir les intensités détectées en données pour déterminer une atténuation par le corps à examiner des rayons X ou lumineux dans un plan de coupe perpendiculaire à l'axe de rotation (SI), et - un ordinateur (21) dûment programmé pour effectuer les étapes suivantes : (1) construire un vecteur générateur colonne et un vecteur générateur ligne dont les termes sont constitués par les données acquises sous deux angles de rotation, de préférence différents de 90 degrés, du support (1) ou du bâti (17) mobile autour de l'axe de rotation (Q), 37(2) construire une matrice initiale (n,m) avec les termes des deux vecteurs générateurs, en affectant à chaque zone élémentaire un terme de ligne et de colonne (Bij) représentant un coefficient d'atténuation et défini par la demi-somme, du terme homologue (i) du vecteur générateur colonne divisé par le nombre (m) de termes du vecteur générateur ligne et du terme homologue (j) du vecteur générateur ligne divisé par le nombre (n) de termes du vecteur générateur colonne, (3) ajuster le coefficient d'atténuation en chaque zone élémentaire par une méthode des moindres carrés en tenant compte, des valeurs de bordure de ligne définies par la somme des termes (Bij) sur chacune des lignes de la matrice initiale et des valeurs de bordure de colonne définies par la somme des termes (Bij) sur chacune des colonnes de la matrice initiale, des termes des vecteurs générateurs comme contraintes de ligne ou de colonne et en utilisant la formule suivante : i Cij = Bij + n / * (pj û Ii Bij) + m J * ( ci û Bij) - nm ( E",PiûEUBU) où, dans cette formule, Cij = la valeur recherchée Bij = la valeur estimée initialement (n) = le nombre de lignes de la matrice initiale (m) = le nombre de colonnes de la matrice initiale n_ 1Cij = pj pour toutes les valeurs de i, la contrainte de la colonne j m Ii j Cij = ci pour toutes les valeurs de j, la contrainte de la ligne i, =1 pour aboutir à une matrice ajustée pour laquelle les valeurs de bordure de ligne et de colonne calculées à l'aide des valeurs ajustées (Cij) sont égales, pour chaque ligne et pour chaque colonne, respectivement aux termes des vecteurs générateurs ligne et colonne, * 38caractérisé en ce que, - le support (1) ou le bâti (17) est mobile autour de l'axe de rotation (S2) de façon telle à irradier ou illuminer le corps (3) sous 360 degrés en une ou plusieurs rotation(s), - le convertisseur (19) convertit en données les intensités détectées en un point du détecteur (11) à chaque instant de la rotation du support (1) ou du bâti (17) mobile, et - l'ordinateur (21) est dûment programmé pour, effectuer l'étape (1) de la façon suivante : (la) traiter les données correspondant à un premier (Al) et un deuxième (B1) arcs représentant chacun une fraction angulaire des 360 degrés, choisie en fonction d'une largeur d'exploration du corps à examiner (3) dans le plan de coupe contenant le point considéré du détecteur (11), le premier (Al) et le deuxième (B1) arcs faisant entre eux un angle, de préférence égal à 90 degrés, (1b) construire, pour la paire d'arcs définie à l'étape (la), le vecteur générateur colonne avec les données correspondant au premier arc (Al) de la paire et le vecteur générateur ligne avec les données correspondant au deuxième arc (B1) de la paire, et pour effectuer les deux étapes suivantes : (4) répéter les étapes (1) à (3) pour un nombre N de paires d'arcs (An,Bn), choisi de façon à traiter l'ensemble des données acquises sous les 360 degrés, les premiers arcs (An), comme les deuxièmes arcs (Bn), des N paires présentant deux à deux un chevauchement angulaire fonction du nombre N de paires d'arcs et de la fraction angulaire des 360 degrés représentée par chaque arc, et (5) traiter les matrices ajustées, obtenues pour les paires d'arcs définies à l'étape (4), par un opérateur de rotation pour toutes les superposer sur une même paire d'arcs (0 -90 ) puis moyenner terme à terme les matrices ajustées et superposées pour aboutir à une matrice de synthèse exprimant une image des coefficients d'atténuation du corps examiné (3) dans le plan de coupe. 39. Dispositif selon la revendication 3, caractérisé en ce que le détecteur (11) a de préférence la forme d'une barrette de détection et comprend une ligne de points disposés suivant un axe parallèle à l'axe de rotation (Q). 5. Dispositif selon la revendication 3 ou 4, caractérisé en ce qu'il comprend une cabine (23,25) comportant des parois en tissu au plomb (23). 6. Ordinateur (21) dûment programmé pour la mise en oeuvre d'un procédé selon l'une des revendications 1 à 2. 7. Programme d'ordinateur pour la mise en oeuvre d'un procédé selon l'une des revendications 1 à 2 lorsqu'il est chargé dans un ordinateur (21). 40