Brevet FR2871911A1 — publié le 20051223

PROCEDE DESTINE A REDUIRE L'EXPOSITION A DES RADIATIONS, DES FAISCEAUX ULTRASONIQUES OU DES IMPULSIONS MAGNETIQUES DANS DES DISPOSITIFS D'IMAGERIE MEDICALE L'invention est relative à un procédé destiné à réduire l'exposition à des radiations, des faisceaux ultrasoniques ou des impulsions magnétiques dans un organisme, dans des dispositifs d'imagerie médicale, et comprend une méthode de traitement d'un signal ou d'un ensemble de signaux en vue de produire dans des conditions très améliorée des images numériques de l'intérieur d'un organisme, pour détecter et traiter le cas échéants les anomalies ou les singularités de l'organisme examiné. Différentes méthodes sont maintenant couramment utilisées dans ce but, en particulier, la radiographie, la scannographie, l'échographie, la résonance magnétique pour obtenir et traiter des informations concernant un ensemble de points ou de zones élémentaires, situés dans un organisme, en vue de produire des éléments d'information relatifs à chaque point, et susceptibles d'être exprimés par des images utilisables par des praticiens à la recherche d'anomalies ou de singularités. De manière classique, on obtient par un calcul approprié, une valeur, associée à chaque point ou zone du coefficient d'atténuation dans le cas du scanner médical, de réflexion dans le cas de l'échographie, ou de la densité de protons dans le cas de la résonance magnétique pour produire une image interne de l'organisme. L'exemple du scanner médical illustre le niveau de l'art en la matière. L'imagerie médicale est en effet un moyen fondamental pour détecter et traiter le cancer, et un certain nombre de maladies graves, ou d'observer les organes internes d'un patient. Mais il est apparut nécessaire à la communauté médicale, d'aller beaucoup plus loin en exprimant la demande de mise au point d'un procédé d'imagerie médicale, à haute définition à partir d'une amélioration significative d'obtention et de traitement des signaux et des images qu'ils expriment. II faut d'abord rappeler, à titre d'exemple, les principes de la scannographie et l'état actuel de l'art en la matière. La Scannographie (ou Tomodensitométrie), a été découverte par un ingénieur de la firme EMI, G. N. Hounsfield, en 1968. Le brevet (US 3 924 131, US 3 919 552) de 1972 s'intitule: A method and apparatus for examination of a body by radiation such as X or gammaradiation". 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 I x la valeur en un point X, On peut écrire la relation suivante: f A (x) dx = F 1=IoèF In=IoE Ln= - = foA(x) dx Dont on déduit par discrétisation: t Ln Io = J A (x) dx = A1X1 + A2X2 + +AnXn In Les valeurs successives, AI, A2i An correspondent aux valeurs de chaque segment définis par X,, X2i X,,. 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 (h2o), le coefficient d'atténuation de l'eau, on peut utiliser une échelle relative telle que: Bn = [An A(h2o)] * 1000 / A(h2o) 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ées, 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. 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 10 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 15 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 20 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, 25 - 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é. On arrive aujourd'hui à des volumes de chaque zone élémentaire de l'ordre du mm 3. C'est cependant loin de l'échelle microscopique puisque le nombre de 5 cellules vivantes est de l'ordre du milliard dans un mm 3. La détection précoce du cancer suppose un gain considérable de la définition. Mais le temps d'utilisation du système pour un patient déterminé ne saurait dépasser un seuil économique évident. Mais surtout l'accroissement du nombre de profils augmente la dose d'irradiation. On admet en effet qu'un nodule cancéreux accélère son développement lorsqu'il provoque une vascularisation associée, ce phénomène apparaît lorsqu'une taille critique est atteinte soit par exemple 50 microns. Dans les procédés classiques la dose d'irradiation et le temps de calcul sont multipliés par 8000, pour atteindre ce niveau de finesse. C'est la raison pour laquelle l'auteur s'est d'abord efforcé d'obtenir un gain de la définition et de trouver une solution qui tienne compte de la capacité croissante des puissances de calculs et de la progression beaucoup plus lente de l'appareillage à RAYONS X. A cet effet, l'invention a pour objet un procédé destiné à réduire l'exposition à des radiations, ou des faisceaux ultrasoniques ou des impulsions magnétiques dans un organisme, dans des dispositifs d'imagerie médicale, grâce à une modification du processus de traitement du signal en vue de l'obtention d'images représentatives de valeurs associée à chaque terme ou zone d'une surface plane ou non, qui consiste: à réaliser une première image par les procédés connus pour retenir une estimation aussi précise que possible de la valeur associée à chaque point ou zone dans une représentation de l'image, par une matrice rectangulaire, en évaluant par interpolation des valeurs entre deux points ou par répartition de la valeur d'un zone dans des micro zones, résultant de la division de la zone, pour aboutir à une matrice amplifiée à partir de la matrice initiale, donnant ainsi une image plus fine de l'image initiale à mesurer de manière précise, dans des balayages physiques selon des séries croisées de balayages explorant l'aire examinée, la variation entre les deux extrémités de chaque balayage de sorte à calculer les valeurs de bordures de la matrice représentative de l'image recherchée, une valeur de bordure étant la somme des termes d'une ligne ou d'une colonne, de la matrice représentative. à ajuster la valeur de chaque terme et donc de chaque valeur associé à un point de l'image recherchée en utilisant un ajustement par la méthode des moindres carrés tenant compte des contraintes relatives aux bordures, en vue d'obtenir une résolution plus fine de l'image. De préférence, on ajuste la valeur de chaque terme en utilisant la formule suivante: Cij=Bij+ n/*(pj EBZJ)+1 *(ci - Em Bij) J j-1 ( E,, P 1 E UBU) où, dans cette formule, Cij = la valeur recherchée Bij = la valeur estimée initialement (n) = le nombre de lignes de la matrice rectangle (m) = le nombre de colonnes de la matrice rectangle n Ei -1= pj pour toutes les valeurs de j, la contrainte de la colonne j Em Cij = ci pour toutes les valeurs de i, la contrainte de la ligne i j=1 et on applique, sur le même ensemble de données, découpé physiquement de manière différente à chaque fois, la même méthode, plusieurs fois pour obtenir une série d'évaluations de la valeur associée à chaque zone de sorte à obtenir la moyenne des mesures et estimer un écart type. On évalue chacun des termes Bij à partir d'un ensemble de balayages physiques espacés de sorte à produire une matrice réduite représentative de l'image, en répartissant chacun des termes de la matrice réduite dans un carré dont le nombre de lignes et de colonnes est égal à k, considéré comme coefficient d'amplification de la matrice initiale pour obtenir la matrice de référence, amplifiée k fois, et conduisant à une première image plus fine, destinées à être ajustée. L'invention s'étend à un dispositif de mise en oeuvre d'un procédé selon l'invention, caractérisé en ce qu'un faisceau de rayon X ou d'ultrasons est déplacé pas à pas de manière linéaire, grâce à un guidage par des actuateurs piézo électriques, de sorte à obtenir des coefficients relatifs à chaque profil correspondant à un pas donné pour réaliser un ensemble d'images, obtenues grâce à deux séries de profils, telles que dans chaque série les profils soient parallèles entre eux, et que les deux séries soient sécantes. De préférence, le faisceau pivote autour d'un axe se déplaçant lui- même linéairement grâce à des actuateurs piézo électriques. Des barrettes de détection multiples sont prévues pour réaliser simultanément plusieurs images parallèles susceptibles d'un traitement individuel ou en trois dimensions. Sont également prévues des plaques ou autres supports de faisceaux, de rayons X ou de lasers ou d'ultrasons, ou des champs magnétiques sur lesquels les actuateurs piézo électriques permettent d'orienter les faisceaux ou de les déplacer, en vue d'obtenir grâce à l'ensemble des détecteurs des informations concernant l'atténuation de chaque profil engendré par un pas du faisceau. Deux groupes de faisceaux peuvent fonctionner, l'un pour l'exploration, l'autre pour le traitement. Le dispositif peut comprendre également des plaques d'obturation perforées, lesquelles peuvent être occultées par deux réseaux croisés de fils ou lamelles déplaçables pour dégager une obturation ou plusieurs par déplacement des plaques ou lamelles. Ainsi la combinaison d'un appareillage à balayage fin, associé à des caches de très faible dimension, positionnables et commutables de manière ultra précise, et, d'un traitement de l'information permettant une évaluation acceptable des valeurs du coefficient d'atténuation de chaque micro zone élémentaire, permet de diviser par un facteur considérable le nombre de jets de faisceaux donc le nombre de profils à enregistrer, tout en portant la définition de un mm à 10 ou 20 microns, sans atteindre un niveau d'irradiation prohibitif. En effet pour une matrice image 1000x1000 la résolution par la 35 méthode standard nécessite un million de profils. Par la méthode exposée ci après il est possible dans un premier temps d'établir une image correspondant à une définition standard correspondant à une matrice 20x20 puis d'en déduire une deuxième image simplement extrapolée pour aboutir à une image 1000x1000, de balayer à nouveau le long de deux axes l'organisme étudié grâce à seulement 2000 profils pour obtenir une estimation convenable du coefficient d'atténuation de chacun des un million de points. Selon ce procédé il suffit d'effectuer 400 profils dans une première étape puis 2000, soit au total 2400 profils, ce qui entraîne 99.76%, de réduction du nombre de profils nécessaire. Le taux d'irradiation est dans ce cas divisé par 417, cette amélioration est encore plus significative pour des aires d'explorations importantes. Pour un examen d'une zone de 20x20 cm, avec une définition de 1 mm, telle que pratiquée actuellement, on obtient les coefficients d'atténuation relatifs à 200x200 point soit 40 000 points, ce qui nécessite l'établissement de 40 000 profils. Si l'on souhaite obtenir une définition de 10 microns, il faut disposer des valeurs des coefficients d'atténuation d'un nombre de points 10 000 fois plus nombreux, soit 400 millions de points, et infliger au patient un taux d'irradiation létal, mais on est ainsi conduit à des délais et des coûts virtuels prohibitifs quelle que soit d'ailleurs le procédé utilisé, échographique ou par résonance magnétique. Dans le procédé ici exposé, il est prévu de réaliser 40 000 profils pour obtenir l'image à faible définition, plus 40 000 profils supplémentaires pour le balayage croisé à haute définition soit au total 80 000 profils au lieu de 400 millions, soit une division par 5000 le nombre de profils et de l'irradiation associée. Mais il reste à résoudre le problème du traitement du signal et des valeurs qui lui sont associées. L'auteur a donc été conduit à mettre au point des méthodes de traitement du signal et de l'information, en vue de résoudre le problème particulier concernant l'évaluation des termes d'une matrice rectangulaire, représentatives des profils de découpage d'une aire déterminée explorée par un scanner lorsque l'on dispose d'estimations approchées relatives à chaque terme, mais d'informations précises concernant la somme de chaque ligne ou colonne de la matrice. Plusieurs méthodes ont été mises au point et explorées pour résoudre le problème rencontré dans différents contextes de l'imagerie. Les recherches ont conduit à une méthode originale, laquelle permet une simplification considérable du processus de calcul, en vue d'obtenir l'image recherchée, représentée par une matrice rectangle. L'ensemble des coefficients d'atténuation des point d'une aire quelconque située dans un plan d'examen, ou sur une surface, peut en effet être représentée par une telle matrice. Elle permet dans le cas de l'imagerie médicale de limiter le nombre de profils à réaliser. Cette méthode qui va ici être exposée permet en effet de produire des matrices de plus grande dimension par une répartition ou une extrapolation des valeurs de chaque terme ou zone d'une matrice initiale de faible définition pour obtenir des valeurs estimées dans des micro zones, résultant du découpage de chaque terme ou zone, de la matrice initiale. On obtient ainsi une matrice amplifiée. Cette matrice amplifiée est alors ajustée par un procédé de calcul permettant d'évaluer de manière fiable chaque terme si l'on peut établir les valeurs des éléments de bordure, un élément de bordure est la somme des termes d'une ligne ou d'une colonne. Dans ce cas, on peut passer par exemple d'une matrice 5x4 contenant 20 termes en une matrice 25x20, contenant 500 termes, puis calculer chacun des termes par le procédé exposé. Il y a bien alors, amplification et ajustement de la matrice initiale représentative des signaux obtenus pour la réalisation des différents profils. La description du procédé d'ajustement est exposée ci-dessous. Ce procédé joue un rôle particulier dans le calcul effectué après, le cas échéant, obtention d'une matrice de grande dimension à partir d'une matrice initiale plus petite. La description qui va suivre expose donc la méthode de calcul d'ajustement proprement dite selon l'invention 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 Bij la valeur estimée à la ligne i et à la colonne j, - si on appelle Cij 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 cj la somme des termes de la ligne i. L'estimation de Bij résulte soit du processus d'amplification matricielle, soit de toute autre méthode permettant une telle estimation, notamment à partir des techniques de l'ajustement linéaire ou polynomial. Dans le cas présent, on recherchera la solution des valeurs de en tenant compte des contraintes de lignes et de colonnes, c'est à dire le minimum de la fonction: E (Cij B02 Pour toutes valeurs de i et de j sous les contraintes: E Cij = pj pour toutes les valeurs de j E 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: L = - (Cij Bij)2 + 1 (Ln' (Cu - pj))+ Ln i_1 i (Ei 1 (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 2j et,sri, 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 dL/dCij 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 m Cij = ci, pour tout i l= Cij = Bij+ + i 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 égal à 25 et m, le nombre de colonnes, est égal à 30, la solution par l'algèbre linéaire consiste à traiter: - 750 variables Cij - 25 variables correspondant aux multiplicateurs de lignes, les p i - 30 variables correspondant aux multiplicateurs de colonnes, les ?L. j Un premier objectif est déjà atteint puisque seulement 55 profils doivent ê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 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: E 1By + 2 * XJ+,u i/2) =Pj Em Bai + 2 * pi+ 1 2 j / 2) = ci On peut déduire de ces relations: 11 *(2*(pj- 1Bij)- t 1,ui) n pi= m*(2(ci-)-E',=1Âj) Dans ces conditions en reportant par exemple la valeur de Xj dans pi, on obtient: Pour tout j = ? 1 J ((P1 - Ei Bij) n - i) n E=1P P Pour tout i pi (1\ n 2l +,n 1 m ni ( i) + m J (ci - L.,=1 Bij - n * ;=1 P j + * Si l'on précise que p- = (1/n) E (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 = ? \ n, *(pj- E:_1 Bi) - p - + I 2J * ( ci - >' Bij - m-' n E OU. )) En effet: m Ls* En pi est égal à p- Dans ces conditions en reportant dans la relation: Cij = Bij + () * (Xj+ pi) 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 Cij = Bij + *(pj InBij)+ *(ci Bij) - 1 * ( n m.1=1 nm EijBJ) 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 d'établissement des valeurs de définition de l'image recherchée sur le plan médical, apparaît ci après. Exemple d'application de la méthode, sur un modèle réduit Soit donc une matrice de n lignes et m colonnes dans laquelle n=3, m=4 TABLEAU 1 MATRICE INITIALE 1 2 3 4 lignes c 1 22 24 18 16 80 78 Pour tout i pi = p- 2871911 14 2 24 22 18 20 84 85 3 26 20 22 24 92 93 colonnes 72 66 58 60 256 P 70 67 59 60 256 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. L'application de la formule ci dessus se simplifie puisque le total des contraintes de colonnes (ou de lignes) est égal à la somme des termes et conduit à. TABLEAU D' EQUILIBRE APRES CALCULS 1 2 3 4 2 C A 1 20.83333 23.83333 17.8333 15.5 78 78 0 2 23.5833 22.58333 18.5833 20.25 85 85 0 3 25.5833 20.5333 22.583 24.25 93 93 0 2 70 67 59 60 256 P 70 67 59 60 A 0 0 0 0 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 satisfont les contraintes, mais aussi qu'il conduisent bien aux résultats recherchés. Equilibrage en utilisant la méthode de l'algèbre linéaire Cette méthode exprime directement les relations linéaires entre les variables Cij et Bij, et les variables 2^, j et pi. 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 multiplication par cette matrice inverse, du vecteur exprimant les seconds membres des relations. a) II existe 12 relations entre les variables résultant de l'expression des dérivées partielles du LAGRANGIEN de la forme: Cij Xj/2-pi/2=Bij b) Il existe 4 relations de contraintes relatives aux colonnes et 3 relations de contraintes relatives aux lignes, de la forme: Em1 c = pj En C'Z.J = et i=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é. L'analyse en deux étapes, correspondant à 6 séquences A) L'amplification matricielle Elle se déroule selon la séquence suivante: 1. Création d'une matrice initiale, traduisant une double série de profils croisés correspondant aux coefficients d'atténuation élémentaire d'une exploration tomodensito-métrique; nous reprenons par exemple le tableau N 1 Matrice initiale (page 11). 2. Mise en place d'un tableau de structure voir en annexe le tableau Na 2 intitulé AVANT AMPLIFICATION , figure 17. 3. Inscription dans les cases à proximité de chacune des valeurs divisée par 25 ce qui conduit au tableau N 3 après AMPLIFICATION , figure 18. B) L'ajustement se déroule ainsi: 4. Inscription des contraintes de lignes et de colonnes ce qui conduit au tableau N 4 AMPLI + CONTRAINTES , figure 19 dans lequel les sommes des lignes et des colonnes, et les contraintes correspondantes sont indiquées. Chaque contrainte résulte de la mesure physique de la valeur totale d'atténuation d'un profil obtenu par lecture de la valeur atteinte sur le détecteur d'extrémité du profil, selon un balayage particulier. 5. Calculs des valeurs par la formule indiquée auparavant, ce qui conduit au tableau N 5 AMPLI + CONTRAINTES + CALCULS , figure 20. 6. On peut alors traduire le tableau N 5, en ce qui concernent les valeurs des Cij, par une image colorée et graduée, constituant le tableau N 6-1, figure 21. On a ainsi abouti au résultat recherché par le praticien d'une image médicalement utilisable dont la définition est meilleure avec un nombre limité d'opérations, et donc une quantité réduite d'irradiation telle que nécessaire par les procédés existants. En réalité toute aire dans un plan, ou sur une surface, conduit à la création d'une matrice rectangle, cette matrice résulte simplement de l'établissement de deux séries de profils parallèles entre eux dans chaque série, que chaque série soit perpendiculaire ou non à l'autre. Elle peut aussi résulter des informations obtenues en coordonnées polaires, avec des séries profils passant par un pole donné, se recoupant le cas échéant avec des profils obtenus à partir d'un autre pole. Rappelons qu'il est toujours possible de ne traiter q'une partie de l'aire de référence sans aucune complication. Il est, en effet, toujours possible de limiter le calcul à une ou plusieurs zones de l'image (Etant donné la nature même du procédé de calcul). Les étapes du processus sont alors les suivantes: a) on met en oeuvre un examen avec un balayage à définition classique par exemple 1 mm b) On obtient une image appropriée c) On définit une aire d'exploration fine d) On met en oeuvre un balayage à haute définition par un dispositif approprié, qui suppose des modes de déplacement spécifiques e) On procède au traitement de l'information de l'aire choisie ce qui suppose: 1) L'estimation des valeurs d'atténuation des micro zones résultant de la division des zones élémentaires, de manière raisonnée, et ce à partir d'une image de référence obtenue par un examen avec un balayage à définition classique, 2) Le calcul des valeurs de termes de la matrice correspondant à l'aire déterminée, par la méthode d'ajustement exposée précédemment, sur la base du balayage à haute définition, 3) Le cas échéant on peut décider la réitération du processus, fondée sur un changement total ou partiel de l'échelle du découpage, ou de la zone explorée, ou le cas échéant sur la réitération du balayage sous des angles différents de la même zone pour obtenir plusieurs valeurs de l'atténuationde chaque micro zone, dont la moyenne procure une évaluation améliorée et peut être associée à un écart type. f) La création d'une base de données des valeurs d'atténuation de chaque micro zone dans les différentes phases de l'examen. 25 g) L'édition, partielle ou complète des résultats, et la visualisation qui en découle de manière plane ou en trois dimensions, en noir gris et blanc ou en couleurs. On peut, si besoin est, réaliser un traitement focalisé, utilisant les mêmes techniques de balayage fin et le cas échéant l'appareillage d'exploration, si la puissance nécessaire doit devenir plus élevée, ce qui est le cas. Il est alors possible d'utiliser un faisceau laser X, (ou tout faisceau de faible section) lequel sera utilisé à plusieurs reprises à un niveau approprié d'amplification de sorte à associer exploration et traitement en diminuant le risque d'irradiation des organes non concernés. En vue d'obtenir l'information nécessaire au traitement qui suivra on peut procéder de deux manières: 1 - Avec deux dispositifs associés a) On peut en effet réaliser une analyse tomodensitométrique avec un scanner existant, comportant une couronne, autorisant un nombre suffisant de balayages, pour permettre une définition de l'image de un mm. On conservera l'information dans la mémoire du système informatique, on évaluera ensuite la valeur des coefficients d'atténuation de chaque micro zone résultant d'une division raisonnée dans chaque zone, et en limitant de manière précise l'exploration supplémentaire. b) Puis on effectuera un balayage à haute définition, dans un deuxième appareil spécialisé de sorte à obtenir, dans une aire déterminée la valeur des coefficients d'atténuation, le long de rayons du faisceau beaucoup plus rapprochés. La configuration de ce deuxième appareil peut être identique à celle du premier, mais les balayages peuvent être limités à des zones angulaires plus petites, l'angle est alors fonction du nombre de pas à réaliser, pour obtenir un ensemble de mesures croisées suffisantes. Le deuxième appareil peut être plus simplement rectangulaire pour matérialiser physiquement une matrice. Dans ce cas deux faisceaux se déplacent perpendiculairement pour établir directement une matrice rectangle qui sera soumise au traitement mathématique de redressement dont la nature a été exposée. c) Le problème est aussi un problème de repérage pour que les zones puissent être traitées par superposition exacte. Deux solutions peuvent alors être mises en oeuvre - Une solution par analyse des images obtenues par les deux dispositifs, l'un et l'autre rotatif en utilisant un logiciel permettant une comparaison précise des coefficients d'atténuations dans chaque image de sorte à permettre une superposition des images obtenues, à la même échelle. - Une autre solution en positionnant par exemple sur la peau des repères physiques, tels que grain de plomb ou d'autres matériaux ayant des propriétés faciles à détecter. Ces qrains sont l'équivalent de bornes topographiques. 2 - Avec un dispositif intégré Ce dispositif permet un fonctionnement à double balayage par 15 exemple: a) un balayage à déplacement de faisceau, avec un pas de l'ordre de 1 mm, pour obtenir une image telle que celle produite par un scanner de type actuel b) séquentiellement (ou simultanément) un balayage à déplacement très limité du faisceau (avec un pas de un dixième à un centième de mm) pour sur la même image produire les deux catégories d'information sur la base de la même topographie Un dispositif décrit ci après autorisera l'utilisation de l'un et ou l'autre 25 des procédés de balayages au prix d'une certaine complication dans la mise en oeuvre des deux balayages, surtout lorsqu'il sont effectués simultanément. On pourra, dans ce cas aussi, compléter le calage des images par repérage physique, en plaçant sur la peau, ou même dans l'organisme des matériaux à très forte absorption. 10 En effet on peut injecter des molécules susceptibles de se concentrer, dans les cancers ou dans les zones de vascularisation associées. Dans le cas du dispositif intégré il suffirait de réaliser un système physique rectangulaire stable, supportant les plaques d'obturation et de détection, associé à deux faisceaux mobiles, se déplaçant perpendiculairement, en s'inclinant, l'un et l'autre indépendamment pour couvrir chacun une zone large de par exemple 90 d'inclinaison. Ou en se déplaçant tout en gardant une position angulaire donnée, la combinaison des déplacements pouvant être ajustée selon des protocoles appropriés au cas à traiter. L'objectif final est de disposer de deux ensembles de données: a) Une matrice de référence à balayage large, permettant d'identifier chaque zone b) Une matrice, résultant des informations produites par balayage fin, susceptible d'être traitée, par les méthodes mathématiques exposées au début de la description du procédé, pour évaluer la valeur des coefficients d'atténuation dans chaque micro zone. La planche N 1 intitulée Schéma d'ensemble montre un exemple de réalisation. Dans cette planche, figure 1, l'accélérateur de particules ou faisceau 1 peut se déplacer. Le plateau de positionnement 2 est déplaçable en hauteur et translatable horizontalement 9. La plaque d'absorption 3 est composée de matériaux lourds pour éviter l'effet des rayonnements. La plaque support des barrettes de détection 4 supporte un ensemble de barrettes permettant la lecture des coefficients d'atténuation en l'absence ou en présence du corps ou de l'objet à examiner. L'obturateur 5 est composé d'un ensemble d'orifices à position et ouvertures ajustables. Ce schéma représente une coupe de l'ensemble du système physique. La table 2 sur laquelle repose le patient est placée à l'intérieur d'un cadre. Le faisceau 1 vertical, traverse une plaque d'obturateurs 5 puis le patient, puis la table, pour atteindre la plaque 4 supportant les détecteurs, laquelle est reliée au dispositif informatique. Nous avons choisi de représenter le cas nouveau, d'un cadre rectangulaire, associé à au moins deux faisceaux inclinables et déplaçables. La Planche N 2 intitulée Faisceaux et obturateurs, vue perspective est une projection de l'ensemble montrant les balayages verticaux et horizontaux vus en perspective, dans une variante, dans laquelle le cadre est posé sur la table et le patient installé sur un matelas; pour lui permettre de pénétrer dans le cadre. Dans cette planche, figure 2, 1 désigne un faisceau vertical, 10 un faisceau horizontal, 2 une table de positionnement, 8 une potence supportant le faisceau horizontal, 4 une plaque de détection horizontale, 6 une plaque de détection verticale et 7 une grille d'obturation horizontale. Dans cette variante le faisceau ne traverse pas la table, cette disposition a pour objectif d'éliminer dans le traitement de l'information, le problème de ladite table. II est bien entendu possible d'utiliser des faisceaux à puissance normale pour limiter les risques, dus au rayonnement. Dans ce cas le traitement supposera que l'on dispose sur le même site des deux types de faisceaux, ou d'un laser X fonctionnant à plusieurs niveaux d'amplification. Pour que l'organisme reçoive le faisceau, de manière précise, on place dans le trajet une plaque de plomb dans laquelle une fenêtre est ouverte; dont la 20 dimension et la découpe est réalisée en fonction de chaque protocole d'analyse ou de traitement envisagé. Dans le cas de notre système on installe au surplus une plaque contenant des orifices très petits de un à deux centièmes de millimètres; ces orifices peuvent être obturés préférablement par déplacement de fils très fin. On constitue ainsi une grille de deux ensembles superposés de fils, horizontaux 12 et verticaux 13, figure 10, se croisant à angle droit. Ces deux ensembles obturent simultanément les orifices. Pour une définition de un mm, les fils sont remplacés par des lamelles de largeur 0.8 à 1 mm, voir planche 5 Système d'obturation , figures 10 à 13. Le déplacement de deux fils, commandés en abscisse et en ordonnée dégage un seul orifice; les déplacements simultanés de deux sous ensembles dégagent une aire rectangulaire (voir la planche 5, figure 10). La figure 9 montre des fils 13 translatés horizontalement. La figure 10 montre une ouverture 14 après une translation verticale des fils 13. Ces déplacements sont obtenus de manière précise par l'utilisation d'actuateurs pièzo électriques. Le dégagement d'une aire précise est effectué pendant un temps très court, de sorte à découvrir les micro zones à examiner ou traiter en une ou quelques milli secondes. Un boîtier 15 commande les déplacement des fils, figure 11. Voir aussi la planche N 3 intitulée Faisceau et obturation . Dans cette planche, figure 3, les obturateurs 5 du faisceau 1 sont composés d'une série de lucarnes fonctionnement simultanément ou indépendamment, groupées par blocs et déplaçables. Différentes formes d'obturateurs sont représentées par les figures 4 à 8. La figure 8 montre plus particulièrement une fente à section variable. Pour éviter les difficultés inhérentes à de tels procédés, il sera préférable d'utiliser un laser à double fonction de positionnement et de visée, 15 analogue à ceux utilisés pour des mesures de distance et de positionnement topographique. Le schéma et la description sont donnés sur la planche N 4 intitulée Laser de repérage . Un tel laser est utilisable lorsque l'on traite une fraction de l'organisme transparent comme l'oeil, ou si l'on veut examiner une surface comme la peau ou un organe interne pendant une intervention chirurgicale. L'utilisation de plaques d'obturation et de détection et de véritables grilles à fonctionnement très rapides est destinée à permettre un cadrage de l'aire à examiner sans déplacement physique de l'objet ou du patient. Dans la planche 4, figure 9, le laser de repérage 16 occupe la position exacte du faisceau X. Le cadre portant le système d'obturation 5 est dans la trajectoire du rayon. La table support 2 supporte aussi le cadre d'obturation et celui de détection. Le matelas de repos 17 s'ajuste sur le niveau du cadre détection. L'armoire électrique 18 relaie toutes les informations au système informatique. La Planche N 6 intitulée Analyse du balayage montre la méthode de 30 balayage utilisée. Figure 14, le faisceau 1 se déplace le long d'une plaque de balayage 19, il effectue une translation pas à pas. Dans le cas de la méthode classique, à chaque pas il s'incline progressivement de sorte à ce que les faisceaux successifs convergent vers un point central O. Chaque balayage permet au faisceau de traverser la plaque d'obturation 5, puis la zone à examiner 20, enfin d'atteindre le détecteur 4 dans une zone alignée dans l'axe défini par la position de l'articulation de faisceau avec la plaque et par l'angle de tir. Le nombre de profils à obtenir doit être égal au nombre de points d'examens élémentaires, pour déterminer avec certitude le coefficient d'atténuation de chaque zone. Une bonne approximation peut être obtenue, par une séquence de balayages fins tels que l'on obtienne un nombre suffisants d'estimations de chaque micro zone pour calculer une moyenne et un écart type convenables, pour chaque terme de la matrice amplifiée. Dans le cas du balayage fin, l'inclinaison reste constante, pour l'un 1 et l'autre 10 des faisceaux, ou pour le même faisceau, si, figure 14, les plaques de support 19,21 du faisceau, de l'obturateur 5,22 et du détecteur 4,6 passent de la position horizontale 19,5,4 à la position verticale 21,22,6 en effectuant un déplacement pour procéder à deux séries d'analyse permettant aux deux balayages de se recouper dans toute la zone à examiner. Le système fonctionne donc en quatre étapes pour le balayage fin 1. Définition de l'aire à examiner par exemple 4x6 mm, 2. Positionnement de cette aire en définissant les abscisses et les ordonnées des quatre sommets, 3. Allumage du faisceau approprié, accélérateur de particules pour le traitement, ou de rayons X pour l'analyse, ou Laser X 25 fonctionnant à plusieurs niveaux, un niveau faible pour l'exploration, et plus élevé pour le traitement, 4. Désobturation rapide par déplacement selon deux axes de coordonnées des fils d'obturation. L'information provenant des détecteurs est transmise en temps réel au système informatique, dans le cas de l'utilisation du même dispositif pour les deux phases de l'opération, qui supposent: - un balayage avec un pas de un mm pour l'ensemble de l'examen - puis un balayage avec par exemple un pas de vingt microns Le premier balayage sera réalisé en déplaçant et inclinant le faisceau selon un protocole standard, les obturateurs s'ouvrant pour permettre l'examen de zones prédéfinies dans différentes conditions standard. Pour une aire de 10x10 cm, on réalisera 10'000 profils pour calculer de manière précise les valeurs de 10'000 points par inversion d'une matrice 100x100; ou dans des conditions beaucoup plus rapide si l'on se contente des valeurs moyennes associées à des écarts types, pour une série de mesures traitées par la méthode d'ajustement exposée précédemment. Cette procédure limite en particulier l'irradiation. Pour le deuxième balayage, les deux faisceaux pourront préférablement garder leurs inclinaisons et se déplacer l'un et l'autre dans un champ limité pour permettre l'exploration d'une aire par exemple rectangulaire de 4x6 mm, on est conduit à générer (pour une définition de 20 microns) une matrice rectangle de 200x300, par plan de coupe soit 60 000 points, dont le traitement numérique restera dans des limites convenables, compte tenu de l'algorithme déjà exposé, tout en limitant l'irradiation. La planche 7 intitulée Balayage spécial présente un intérêt particulier, elle montre qu'il est possible d'utiliser une seule plaque de balayage 19, par exemple horizontale, et un seul plan de détection 4 ici horizontal. Figure 15, on a représenté un trapèze de sommets A, B, C et D. Les références 23 et 24 désignent respectivement des lignes horizontales et verticales. Dans ces conditions, figure 16, l'objet 20 est balayé, dans le plan tomographique, par exemple sous deux angles a et f3 tels que les balayages se recoupent en un nombre suffisant de points pour obtenir une évaluation pour chaque point éloigné du plus proche de 1 mm, et obtenir ainsi une évaluation précise. Ce travail sera complété par au moins deux balayages par translation avec des pas de 20 microns environ ou par rotation, avec une fraction de degré par pas, pour obtenir les données nécessaires au traitement par la méthode déjà exposée. Les angles des deux groupes de balayage peuvent être ou non égaux à 90 . La complication résultant de l'utilisation d'axes de coordonnées non orthogonaux est négligeable sur le plan du calcul. Nous rappelons que le dispositif peut comprendre un faisceau LASER A RAYONS X dont le rayon va traverser l'objet ou le patient à examiner. Après la traversée un détecteur permet de lire la puissance résiduelle du faisceau en vue de mesurer le coefficient d'atténuation pendant la traversée de l'objet ou du patient examiné. Ce faisceau progresse point par point devant une batterie d'obturateurs, pouvant être constituée par des petits caches déplaçables en face d'orifices de très faible section. II est possible d'utilise un faisceau suffisamment large pour couvrir une batterie d'obturateurs, et affecter une barrette homologue de détecteurs. Si l'on admet que certains déplacements naturels (par exemple la respiration) puissent déplacer des éléments d'organe à une vitesse de 1 cm par seconde; et si l'on veut obtenir une précision de l'image de 10 à 20 microns, le temps de fonctionnement de l'obturateur doit être de l'ordre de la milli seconde, et le dispositif doit créer une série de flashes successifs de très courte durée. Si donc le faisceau doit être établi et interrompu très rapidement (dans le cas notamment de l'utilisation d'un laser lumineux ou à rayons X), ou si l'obturation doit être modifiée rapidement le fonctionnement de la commande (électrique ou mécanique) pendant une période très courte se heurte aux limites des systèmes électromagnétiques, sur le plan technique ou des systèmes électroniques, sur le plan des coûts. Dans un brevet européen n 99401358 déposé par Norbert BEYRARD, concernant un contacteur à commande piézo électrique l'auteur décrit un dispositif de commande d'un contacteur disjoncteur dont les déplacements sont effectués par un actuateur piézoélectrique. Les essais les plus récents montrent: - Que la coupure du courant peut s'effectuer en un temps de l'ordre de la milli seconde. - Que la coupure de courant peut grâce à un circuit électronique approprié, s'effectuer lorsque le courant passe, par exemple en courant alternatif, au point ou l'intensité est nulle. On peut sans extrapoler prévoir un dispositif de commande comportant: - un alternateur (ou un moduleur) délivrant un courant de fréquence égal par exemple à 200 ou 400 Hertz. - un circuit électronique permettant que la mise en circuit s'effectue au moment d'un premier passage du courant à intensité nulle, et que la coupure s'effectue au passage suivant. Dans ces conditions, le dispositif sera allumé, ou désobturé, pendant une demi période égale à 1/800ème de seconde, soit 1.25 milli seconde, ou de 1/400ème de seconde ou 2.5 milli secondes. Dans le dispositif mis au point par Norbert BEYRARD et qui comprend un circuit électronique de commande, il est possible d'enclencher, puis de déclencher le circuit électrique pour que la période de fonctionnement soit ajustée de manière précise, tout en maintenant le fait que l'ouverture ou la fermeture du circuit s'effectue au moment où l'intensité du courant alternatif (sinusoïdal) passe par le point d'intensité nulle. Il est ainsi possible de provoquer l'ouverture et la fermeture d'un circuit à un moment précis dans le temps. Si donc le déplacement de l'objet ou de la zone à examiner s'effectue à 1 cm par seconde, il pourra être examiné au cours d'une période de 1 à 2 milli seconde, correspondant à un déplacement de l'ordre de 10 à 20 microns. Pour tenir compte des aléas le faisceau serait d'un diamètre de 100 microns (ce qui est normalement le cas pour certains laser X). II s'agit là d'un point essentiel qui autorise la recherche d'une définition de l'ordre de 10 à 20 microns au lieu de 1 mm actuellement, soit un gain de 50 à 30 100 de la définition par rapport à la technique actuelle. 10 Le gain économique qui en résulte peut s'exprimer simplement. Un tel dispositif de coupure, tel que décrit dans le brevet précité, ne devrait coûter qu'une dizaine d'euros, pour une fabrication en série ayant éventuellement d'autres objets; la coupure électronique par exemple par un thyristor suppose un système de coupure dont le prix est environ cent fois plus élevé. Si l'on dispose de systèmes ayant une telle précision, la quantité d'informations produite par le détecteur est multipliée dans des proportions considérables et dépasserait rapidement la limite logiciels et des calculateurs généralement associés aux scanners actuels. Les calculateurs associés au nouveau dispositif devraient traiter un ensemble d'équations linéaires et de variables 125 000 fois plus important pour un gain de définition de 50. Prenons un objet dont la section est de 20 cm sur 20 cm pour une zone de 1 mm. On est conduit à calculer 200 x 200 coefficient d'atténuation par section, soit 40.000 points. En dépit des astuces de programmation le temps et le coût des calculs sont encore aujourd'hui considérables. Si l'on améliore la définition dans un rapport de 50 il faut multiplier dans une proportion considérable (de l'ordre de 125 000) la vitesse, et donc la puissance de calcul actuelle. Si l'on considère que la Loi de Moore continuera à s'appliquer dans le futur, que la puissance de calcul double donc tous les 18 mois, il faudrait attendre 25 ans pour que les systèmes deviennent opérationnels. C'est la raison pour laquelle il est nécessaire d'utiliser un autre procédé tel que décrit précédemment, en vue d'accélérer le processus. Le processus de calcul est alors le suivant: Au cours d'une première étape on s'efforce d'estimer les coefficients d'atténuation de zones réduites. Admettons que l'on explore une zone limitée prise à l'intérieur d'une section explorée de manière classique, et que cette zone ait par exemple 4 mm sur 3mm. Dans une telle zone nous avons 12 zones, pour lesquels différents coefficients d'atténuation ont été calculés, de manière classique ou non. La représentation de valeur de chaque zone est exprimée dans le tableau 1 Matrice initiale page 2. Chaque zone est alors découpée en vingt cinq (5x5), c'est la raison pour laquelle dans le tableau N 2 AVANT AMPLIFICATION la valeur du coefficient d'atténuation apparaît au centre de chaque zone. Le coefficient relatif à la zone initiale N 1 amplifiée est égal à1125, de 22 soit 0.88, on voit, sur le tableau 3 Matrice après amplification la matrice amplifiée, qui montre les coefficients d'atténuation de chaque micro zone après répartition. Les valeurs des contraintes sont indiquées au tableau N 4 Après ampli + contraintes , et complètent les informations produites au tableau N 3. Nous avons testé la méthode d'ajustement décrite auparavant sur la 15 base du tableau 4, et vérifié son fonctionnement en établissant le tableau 5 ampli + contraintes + calculs , qui montre que la somme des coefficients d'atténuation estimée est bien égales à la valeur de la contrainte. L'examen de ce tableau 5 permet de suivre le processus de calcul: 1. Les lignes numérotées indiquent la valeur estimée de chaque 20 terme, elles sont les mêmes que celles homologues du tableau 4 (Les Bij) 2. La colonne 21 indique la somme des valeurs estimées du tableau 3. La colonne 22 les contraintes de lignes 25 4. La colonne 23 l'écart entre les valeurs de 22 et 21 De la même manière en bas du tableau apparaissent les différences entre contraintes de colonnes et somme des valeurs estimées pour chaque colonne. Enfin, sur toutes les lignes en gras situées au dessus des valeurs des lignes numérotées, apparaissent les valeurs calculées grâce à une macro instruction qui lie les valeurs (Cij) résultant du calcul, aux valeurs indiquées sur les lignes numérotées (les Bij). Par exemple: Si l'on estime les coefficients indiqués en ligne 3, tirés de la matrice d'origine et ceux apparaissant dans la ligne au dessus, on vérifie sur la colonne 21 que la somme est bien égale à 13,987, soit une valeur très proche de celle de la contrainte égale à 14 qui apparaît dans la colonne 22. Ce calcul effectué pour l'ensemble des coefficients du tableau montre que la méthode ainsi découverte permet le calcul dans des conditions d'extrême rapidité ; puisqu'une simple formule algébrique se substitue aux formules matricielles résultant normalement de ce processus de calcul. Les calculateurs actuels peuvent traiter le problème correspondant et même aller bien au delà en appliquant des processus de calcul similaires à ceux exposés précédemment. 1. Lorsque l'on procède au calcul d'équilibrage avec un tableau obtenu par répartition des valeurs de chaque zone en par exemple 25 micro zones, on peut se trouver loin des valeurs de contraintes de ligne ou de colonnes de deux manières: 2. Si l'on ne dispose pas des valeurs de bordures, celles ci peuvent être évaluées par répartition des valeurs connues sur la ligne ou la colonne, grâce à l'étape de définition standard (1mm). 3. On peut procéder alors à un calcul pas à pas avec des pas réduits à 20 microns et des temps de fonctionnement de l'ordre de la milli seconde. Cette méthode mettra en évidence les écarts entre estimation et observation. Le cas échéant, on effectuera plusieurs balayages fins sous des angles différents, on disposera ainsi de plusieurs évaluations des Cij, conduisant à une amélioration de la valeur recherchée. 30 Si les écarts entre estimation et observation sont inférieurs à un seuil fixé à l'avance, on admettra par exemple un écart relatif de 20 %, et on procédera sans risque au calcul par les méthodes exposées ci dessus. Si les écarts par contre sont plus élevés et s'ils restent inférieurs à un deuxième seuil, par exemple +1- 40% d'écart relatif, on diminuera le taux d'amplification en passant par exemple d'une division de chaque zone en 25 micro zones à une division de chaque zone en 9 micro zones. Si en dépit de cette nouvelle répartition des écarts persistent à l'intérieur de l'aire considérée, on procédera à un traitement selon les modalités classiques pour calculer les coefficients d'atténuation avec certitude, ou de manière plus économique par la méthode utilisant une séquence de balayages. Il est donc possible aujourd'hui d'envisager de traiter avec des calculateurs puissants mais disponibles, des mini zones de 1 cm x 1 cm, soit 250 000 micro zones de 20 microns de coté, au lieu de 1 mm (la résolution actuelle), Cette technique de résolution conduit à une amélioration de la définition d'un facteur 50, et ce en tenant compte des limites physiques de la dimension de l'obturateur et temporelles du fonctionnement du faisceau ou l'interruption de l'obturateur. L'exemple des tableaux 3, 4 et 5 proches de ce qui pourrait être envisagé conduirait à mesurer 35 valeurs de bordures du faisceau pour 300 valeurs estimées des coefficients d'atténuation. 9. Le gain sur le plan économique et de santé est encore plus spectaculaire avec des tableaux plus importants. On est cependant encore loin de la situation idéale permettant 30 d'atteindre le niveau d'une micro zone de 1 micron, c'est à dire la taille d'une cellule. Chacune des micro zones de 20 microns contient en effet 8 000 cellules. 4. 5. 6. 7. 8. Dans l'immédiat, ces calculs pourraient atteindre la limite des calculateurs associés au scanner, pour des examens de zones importantes, par exemple 20 x 20 cm. Il est à nouveau rappelé, qu'il est possible d'effectuer dans la même zone, plusieurs balayages, par exemple une dizaine, qui donneront pour chaque point, une dizaine de valeurs voisines dont on pourra ainsi calculer la valeur moyenne et l'écart type, ce qui donne dans la plupart des cas, une estimation tout à fait convenable et évite le recours au calcul matriciel qui entraîne des temps de calcul prohibitifs. Pour être utilisée, l'information doit être éditée de manière appropriée. On peut procéder de manière habituelle en éditant des cartes faisant apparaître pour chaque coupe des zones plus ou moins grises en fonction du coefficient d'atténuation, et ce de deux manières: - de manière standard (ou par la méthode des balayages 15 multiples) pour l'ensemble de la coupe avec une résolution de 1 mm - de manière complémentaire pour chaque mini zone traitée plus en détail avec différents tons de gris, ou différentes couleurs. On prévoira d'agrandir éventuellement les mini zones pour faciliter la 20 lecture. On prévoira avantageusement de tracer des lignes d'iso atténuation (courbes rejoignant des micro zones de même niveau d'atténuation) pour voir apparaître les anomalies particulières. Grâce à un traitement coloré on peut envisager superposer les lignes 25 d'iso atténuation sur les fonds à gris variables. La planche 6-1 en annexe montre un exemple saisissant puisque outre lacoloration, des lignes de niveaux d'atténuation sont établies, ce qui affine encore le résultat de l'investigation. Il faut aussi tenir compte du balayage axial, c'est à dire du déplacement 30 du plan de coupe, qui ne nécessite plus aujourd'hui de véritable déplacement physique, puisqu'on peut disposer des barrettes comprenant un nombre suffisant de détecteurs sur une longueur de plusieurs centimètres. Le problème est aussi un problème de traitement de l'information puisque le déplacement axial devrait être du même ordre de grandeur que le balayage dans un plan de coupe, ce qui conduit à des calculs qui même largement simplifiés par les méthodes exposées, ne doivent pas devenir prohibitifs. Lorsque pour des raisons d'exploration fine on sera conduit à rapprocher les coupes on pourra réaliser une projection à trois dimensions des informations, particulièrement pour les micro zones, tout aussi bien que pour des zones correspondantes à la définition actuelle. L'auteur est donc conduit à prévoir un système opérant ainsi: II faut ainsi: a) Faire fonctionner un scanner avec une résolution de 1 mm pour obtenir des mesures relatives à des zones de même dimension qu'actuellement; le calcul peut être effectué de deux manières, soit par l'utilisation de l'algèbre linéaire, ce qui produira une valeur précise du coefficient d'atténuation de chaque zone, soit en utilisant la méthode exposée initialement qui suppose une (ou plusieurs) dizaine de balayages pour obtenir une valeur moyenne et un écart type pour chaque point. II est d'ailleurs possible de procéder à un nombre plus élevé de balayage pour atteindre une précision suffisante pour tout point b) Définir des zones de dimension plus limitée pour un examen plus fin (par exemple de 1 cm de coté, comportant 100 mini zones par zone, soit une amplification de coefficient de 10, ou se contenter de seulement 25 mini zones par zone) c) Répartir dans chaque mini zone la valeur obtenue pour chaque zone par exemple en 25 micro zones (voir tableau 3, 4 et 5). Estimer ainsi les valeurs des coefficients d'atténuation de chaque micro zones et les valeurs des bordures d) Mesurer pas à pas les valeurs de bordures 20 30 e) Procéder à un traitement mathématique et numérique fondé sur l'utilisation de la méthode décrite au début de ce texte et montrée dans le tableau 5 f) Procéder éventuellement en deux ou plusieurs étapes, pour les opérations précédentes pour affiner progressivement la définition Procéder à un examen comparatif des valeurs des coefficients d'atténuation estimés par la méthode de répartition proportionnelle et ceux calculés par la méthode exposée. Si une différence importante apparaît ou une micro zone précise ou pour un ensemble de micro zones contiguës, on réduira le nombre de micro zones par zone et l'on procédera à un nouveau calcul de répartition et d'estimation. Si la difficulté persiste on procédera pour l'aire considérée à un calcul des coefficients d'atténuation par les méthodes classiques h) Editer le document de référence (définition 1 mm) et les documents obtenus par haute définition (définition 20 à 100 microns) et ce dans des présentations nouvelles telles que décrites ou envisagées ci dessus. Le procédé dont le fonctionnement est ainsi décrit autorise l'obtention d'une définition pouvant être 100 fois plus fine que la résolution actuelle, en utilisant des techniques et des outils de calcul à notre portée. Trois images ont été éditées tableaux 6-2, 6-3 et 6-4, respectivement figures 22, 23 et 24, montrant pour une aire réduite des images correspondant aux définitions de 20, 12.5, 10 microns. Dans le dernier cas l'image traduit un processus de calculs relatif à la résolution d'un problème d'algèbre linéaire de 40400 équations et variables, en ne nécessitant que 400 balayages secondaires, plus quatre balayages primaires soit une réduction de 99% de l'irradiation, laquelle irradiation diminue relativement lorsque la dimension de l'aire primaire croît. Comme ceci a déjà été exprimé l'examen d'une section de 20x20cm avec une définition de 10 microns, détermine 400 millions de points soit autant de profils et nécessite l'inversion d'une matrice 20millex20mille par le procédé classique; et seulement 80 mille profils soit cinq mille fois moins en deux étapes et g) 10 un calcul algébrique à notre portée, notamment pour des zones réduites choisies de manière raisonnée ce que le procédé autorise sans aucune complication. Les progrès techniques en matière de mécanique de déplacement de positionnement de commutation électrique et enfin de calcul permettront sans changer les fondements du système décrit d'atteindre une résolution de un micron soit la détection de cellules anormales de manière individuelle. D'ores et déjà on peut utiliser les techniques exposées par l'auteur dans un brevet français N 00401532, relatif à un moteur piézo électrique à micro crémaillère. Ce procédé consiste à agir par des barrettes pièzo électriques à travers une pointe très fine sur un élément mobile dont la surface est striée de manière très fine de sorte à constituer une micro crémaillère. On obtient ainsi des déplacements avec par exemple des pas de 100 microns. Mais on peut aussi utiliser des actuateurs piézo électriques de la société CEDRAT TECHNOLOGIES en France ou NANOMOTION en Israël pour obtenir une série de déplacements avec des pas inférieurs, tout en disposant d'une possibilité d'effort élevé. A titre d'exemple il est envisageable d'utiliser des actuateurs piézo électriques de CEDRAT TECHNOLOGIES connus sous l'appellation APAs. Dans le cas présent, il est aussi possible de transmettre un mouvement rotatif grâce à des dispositifs tels que décrits par CEDRAT TECHNOLOGIES, sous l'appellation RPMs. Les méthodes de traitement des signaux pour obtenir des images numérisées, susceptibles de bénéficier des méthodes de traitement exposées au début peuvent à l'évidence être utilisées dans d'autres domaines que l'imagerie médicale en particulier en ce qui concerne: - la cartographie des sols, des sous sols, et des pressions atmosphériques - la photographie et la cinématographie numérique - le traitement et la transmission d'images notamment dans le domaine de la télévision En ce qui concerne l'imagerie médicale obtenue par échographie ultrasonore ou résonance magnétique nucléaire, le procédé d'amélioration de l'image est de même nature, car comme il a été expliqué plus haut la matrice représentative d'une image ou d'une aire déterminée peut être constituée de données associées à des séries de profils centrés sur un ou deux pôles, ou exprimés en coordonnées orthogonales ou obliques. Le procédé d'amplification et ou d'ajustement matriciel conduit de la même manière à une amélioration sensible de la définition. Le traitement par ultrasons peut aussi être associé aux explorations. REVENDICATIONS 1. Procédé destiné à réduire l'exposition à des radiations (1), ou des faisceaux ultrasoniques ou des impulsions magnétiques dans un organisme, dans des dispositifs d'imagerie médicale, grâce à une modification du processus de traitement du signal en vue de l'obtention d'images représentatives de valeurs associée à chaque terme ou zone d'une surface plane ou non, qui consiste: - à réaliser une première image par les procédés connus pour retenir une estimation aussi précise que possible de la valeur associée à chaque point ou zone dans une représentation de l'image, par une matrice rectangulaire, en évaluant par interpolation des valeurs entre deux points ou par répartition de la valeur d'un zone dans des micro zones, résultant de la division de la zone, pour aboutir à une matrice amplifiée à partir de la matrice initiale, donnant ainsi une image plus fine de l'image initiale - à mesurer de manière précise, dans des balayages physiques (19,21) selon des séries croisées de balayages explorant l'aire examinée (20), la variation entre les deux extrémités de chaque balayage de sorte à calculer les valeurs de bordures de la matrice représentative de l'image recherchée, une valeur de bordure étant la somme des termes d'une ligne ou d'une colonne, de la matrice représentative. - à ajuster la valeur de chaque terme et donc de chaque valeur associé à un point de l'image recherchée en utilisant un ajustement par la méthode des moindres carrés tenant compte des contraintes relatives aux bordures, en vue d'obtenir une résolution plus fine de l'image. 2. Procédé exposé selon la revendication N 1, caractérisé en ce que: - on ajuste la valeur de chaque terme en utilisant la formule suivante: r 1 [nm 1 Cij=Bij+ n*(pj El Bij)+m*(ci Em Bij)-;_' J E;,Pj EO) ( * où, dans cette formule, Cij = la valeur recherchée Bij = la valeur estimée initialement (n) = le nombre de lignes de la matrice rectangle (m) = le nombre de colonnes de la matrice rectangle En_ 1Ci/ = pj pour toutes les valeurs de j, la contrainte de la colonne j Em Cij = ci pour toutes les valeurs de i, la contrainte de la ligne i j=1 et l'on applique, sur le même ensemble de données, découpé physiquement de manière différente à chaque fois, la même méthode, plusieurs fois pour obtenir une série d'évaluations de la valeur associée à chaque zone de sorte à obtenir la moyenne des mesures et estimer un écart type. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que l'on évalue chacun des termes BU à partir, d'un ensemble de balayages physiques espacés de sorte à produire une matrice réduite représentative de l'image, en répartissant chacun des termes de la matrice réduite dans un carré dont le nombre de lignes et de colonnes est égal à k, considéré comme coefficient d'amplification de la matrice initiale pour obtenir la matrice de référence, amplifiée k fois, et conduisant à une première image plus fine, destinées à être ajustée. 4. Dispositif de mise en oeuvre d'un procédé selon les revendications 1, 2 ou 3, caractérisé en ce qu'un faisceau de rayon X (1,10) ou d'ultrasons est déplacé (19,21) pas à pas de manière linéaire, grâce à un guidage par des actuateurs piézo électriques, de sorte à obtenir des coefficients relatifs à chaque profil correspondant à un pas donné pour réaliser un ensemble d'images, obtenues grâce à deux séries de profils (23,24) , telles que dans chaque série les profils soient parallèles entre eux, et que les deux séries soient sécantes. 5. Dispositif selon la revendication 4, caractérisé en ce que le faisceau (1) pivote autour d'un axe se déplaçant (19,21) lui même linéairement grâce à des actuateurs piézo électriques. 6. Dispositif selon la revendication 4 ou 5, caractérisé en ce qu'il comprend des barrettes de détection multiples (4) pour réaliser simultanément plusieurs images parallèles susceptibles d'un traitement individuel ou en trois dimensions. 7. Dispositif selon la revendication 6, caractérisé en ce qu'il comprend des plaques (19,21) ou autres supports de faisceaux, de rayons X ou de lasers ou d'ultrasons, ou des champs magnétiques sur lesquels les actuateurs piézo électriques permettent d'orienter (a,I) les faisceaux ou de les déplacer, en vue d'obtenir grâce à l'ensemble des détecteurs (4,6), des informations concernant l'atténuation de chaque profil engendré par un pas du faisceau. 8. Dispositif selon les revendications 4, 5, 6 ou 7, caractérisé en ce qu'il comprend deux groupes de faisceaux fonctionnant, l'un pour l'exploration, l'autre pour le traitement. 9. Dispositif selon l'une des revendications 4 à 8, caractérisé en ce qu'il comprend des plaques d'obturation (5,7,22), perforées, lesquelles peuvent être occultées par deux réseaux croisés de fils (12,13) ou lamelles déplaçables pour dégager une obturation (14) ou plusieurs par déplacement des plaques ou lamelles.