Brevet FR2670039A1 — publié le 19920605

PROCEDE ET DISPOSITIF DE RECONSTRUCTION D'IMAGES TRIDIMENSIONNELLES D'UN OBJET EN UTILISANT DEUX TRAJECTOIRES CIRCULAIRES D'ACQUISITION DESCRIPTION L'invention concerne un procédé et un dispo- sitif de reconstruction d'images tridimensionnelles d'un objet à L'aide d'un réseau bidimensionnel de capteurs parcourant deux trajectoires circulaires d'acquisition des mesures ou, de manière équivalente, deux réseaux bidimensionnels de capteurs parcourant chacun une trajectoire circulaire Elle peut encore être élargie à un nombre quelconque de trajectoires. L'image de l'objet est définie par des valeurs prises par une fonction sur chacun de ses points La fonction est d'autre part une propriété d'un rayonnement (rayons X ou gamma entre autres) de forme conique ayant un point focal et qui passe à travers l'objet; chaque rayon est reçu par un des capteurs du réseau bidimensionnel et représente donc la somme de la fonction sur tous les points de l'objet appartenant à ce rayon Un traitement approprié des sommes sur tous les rayons, et ceci pour un nombre suffisant de mesures suivant des incidences différentes autour de l'objet, permet de reconstituer l'image de l'objet. On se contente en pratique de considérer des nombres finis de rayons et de points de l'objet selon des discrétisations ou des maillages. Cette invention constitue un perfectionnement d'une invention précédente, décrite dans la demande de brevet européen'EP-A-O 292 402, mais on peut envisa- ger d'employer l'invention actuelle dans d'autres circonstantes ou en utilisant d'autres méthodes mathé- matiques de traitement des données. Les méthodes envisageables utilisent en particulier ce qu'on appelle la transformée de Radon de la fonction à mesurer, qui est définie comme la somme de la fonction sur chacun des plans, appelés plans de Radon, qui passent à travers l'objet considéré ou, mieux encore, la dérivée première de cette trans- formée La contribution des points de ces plans qui n'appartiennent pas à l'objet est considérée comme inexistante, ce qui est valable dans le cas d'un rayon- nement qui traverse un milieu gazeux ambiant sans être atténué Ici encore, une discrétisation est effectuée pour n'effectuer les calculs que sur un nombre fini de plans. La dérivée première de la transformée de Radon est définie comme la dérivée de la transformée de Radon en fonction de la variable p définissant la distance du plan considéré à une origine Pour chaque plan, elle correspond à la somme sur ce plan de la dérivée première de la transformée de Radon dans la direction perpendiculaire au plan l'invention décrite dans le brevet EP-A-0 290 402 montre que pour un plan passant par une position du point focal du cône et rencontrant le réseau bidimensionnel de capteurs, il est possible de calculer à partir des mesures la valeur exacte de la dérivée première de la transformée de Radon sur ce plan Ce calcul, suivant les formules exposées dans le brevet EP-A-0 292 402, met en oeuvre de façon préférentielle une pondération de correction d'éLoignement du point focal, deux filtrages correspon- dant aux calculs des dérivés premières respectivement suivant les lignes et les colonnes du réseau de cap- teurs, deux sommations suivant la droite d'intersection entre le plan à traiter et le réseau bidimensionnel de capteurs, puis une combinaison linéaire et une normalisation des résultats Cette sommation met en oeuvre Les interpo Lations nécessaires car ces droites d'intersection passent entre Les rangées de capteurs ou Les coupent. Dans Le reste de La demande, on entend par somme La somme pondérée des valeurs mesurées (dans Le cas de La transformée de Radon) aussi bien que La combinaison Linéaire des sommes des va Leurs pondé- rées et filtrées (dans Le cas de La dérivée première de La transformée de Radon) obtenues selon Les Lignes et Les colonnes du réseau. La somme de la fonction sur les points des p Lans de Radon est faci Le à obtenir à condition que ces plans aient une intersection avec le réseau bi- dimensionnel de capteurs et qu'ils passent par le point foca L unique que visent les capteurs IL suffit d'effectuer La somme des valeurs mesurées par chaque capteur situé à l'intersection, avec Les interpolations nécessaires car Les intersections des plans de Radon passent entre des rangées de capteurs ou Les coupent. Une fois que Les valeurs de la fonction sur les plans de Radon ont été calculées, il existe des formules d'inversion, exposées dans La demande de brevet men- tionnée ci-dessus, qui permettent de parvenir aux valeurs de la fonction sur Les points du maillage de l'objet. Il faut toutefois revenir sur les conditions qui permettent d'obtenir un nombre de plans de Radon suffisant pour permettre une description satisfaisante de L'objet Chaque plan de Radon peut être défini par ce qu'on appelle son point caractéristique, c'est- à-dire le point de projection sur ce plan d'une origine choisie arbitrairement Ce point caractéristique, noté C sur la figure 1, peut être défini par ses coordonnées sphériques P, 4 et O de rayon, longitude et colatitude respectivement à partir de l'origine 0 Le plan de Radon P passant par Le point caractéris- tique C peut être défini par Le rayon p et le vecteur unitaire de direction OC. Les valeurs de la fonction sur les plans de Radon ne peuvent effectivement être calculées que pour les plans de Radon qui coupent la trajectoire parcourue par le point foca L du rayonnement Ce point focal est concrètement, dans le cas d'une fonction d'atténuation, une source ponctuelle de rayons X, gamma, etc C'est la conception qui est décrite dans la demande de brevet européenne mentionnée Les mêmes conditions géométriques existent en tomographie d'émis- sion, quand la fonction à mesurer est l'activité émise par le corps Le point foca L n'a alors pas d'existence physique et correspond simplement au point de conver- gence vers lequel tous les collimateurs utilisés devant le réseau bidimensionnel sont dirigés. Il faut, pour obtenir un nombre suffisant de plans de Radon, effectuer plusieurs mesures avec des positions différentes du point focal. Considérons une trajectoire circulaire T parcourue par le point focal (ou la source) S, le réseau plan de capteurs Pdét parcourant la même trajec- toire que le point focal S ou éventuellement une tra- jectoire concentrique de rayon différent Supposons encore que l'origine O qui sert à définir les points caractéristiques C soit confondue avec le centre de la trajectoire T Le volume enveloppant les points caractéristiques correspond à un tore To produit par la rotation d'une surface sphérique de diamètre OS autour de l'axe de rotation de la trajectoire T En effet, les plans de Radon qui passent par le point focal S ont leurs points caractéristiques répartis sur la surface sphérique de diamètre OS car l'angle SCO est droit Le tore To est appelé le volume carac- téristique des mesures, qui dépend donc de la forme de la trajectoire T et de sa position par rapport à l'origine O. Il est par ailleurs suffisant, pour obtenir une description comp Lète de l'objet par les procédés utilisant Les sommes de la fonction sur les plans de Radon, de pouvoir disposer des points caractéris- tiques appartenant à un volume caractéristique de l'objet Ce volume caractéristique de l'objet est toujours inclus dans une sphère V centrée sur l'origine O et qui englobe l'objet M On peut donc assurer que si le volume caractéristique de l'objet est inclus dans le volume caractéristique des mesures, la recons- truction de l'image de l'objet M sera possible Dans le cas d'une trajectoire circulaire T, cette condition n'est malheureusement pas remplie car une zone d'ombre subsiste pour les points caractéristiques de la sphère V qui n'appartiennent pas au tore To et dont les plans de Radon ne coupent pas la trajectoire T. La demande de brevet précédente contient la constatation que le choix de certaines trajectoires plus compliquées permet d'agrandir le volume caracté- ristique des mesures pour remplir complètement le volume caractéristique de l'objet Un cas intéressant est celui de la superposition de deux trajectoires circulaires concentriques et appartenant à des plans différents, par exemple perpendiculaires, car la super- position des deux tores concentriques permet d'éliminer la zone d'ombre autour de l'origine Un dispositif adapté à de telles mesures comprend un rail ou une glissière matérialisant la trajectoire T, sur lequel un émetteur de rayons X et le réseau bidimensionnel de capteurs circulent de manière à être toujours diamé- tralement opposés, et qui peut basculer autour d'un axe fixe Aucune indication n'avait toutefois été donnée sur la façon de reconstruire effectivement l'image de l'objet M à partir des deux séries de mesures et de leur synthèse La présente invention a précisément pour sujet un procédé adapté à la recons- truction d'images au moyen de mesures prises suivant deux trajectoires circulaires, sans qu'aucune condition ne soit exigée sur Les relations géométriques entre ces trajectoires. Pour récapituler, L'invention concerne un procédé de reconstruction d'images tridimensionnelles d'un objet défini par des valeurs prises par une fonc- tion sur des points de l'objet, la fonction étant une propriété d'un rayonnement conique ayant un point focal et passant à travers l'objet, dans lequel la fonction est calculée par l'intermédiaire de sommes de la fonction sur des plans passant par au moins un point de l'objet et définis dans un repère d'objet, comportant au moins deux séries de mesures, chaque somme de la fonction étant calculée, sur lesdits plans, à partir d'au moins une des séries de mesures, chacune des séries de mesures étant menée avec un réseau bi- dimensionne L de capteurs de rayonnement, orientés vers le point focal, qui est déplacé autour de l'objet sur une trajectoire circulaire respective, caractérisé en ce que, pour chaque série de mesures, une série des sommes de la fonction est calculée sur certains au moins desdits plans qui sont sécants ou tangents à la trajectoire associée à ladite série de mesures et sont alors définis dans un repère de trajectoire associé à Ladite trajectoire, et en ce que les séries de sommes de la fonction sont regroupées en effectuant des ca Lculs pour pouvoir définir les plans dans le repère d'objet. Ainsi, chaque série de mesures de la fonction peut être calculée d'abord avec les paramètres (inci- dence de la mesure, distance entre le point focal et le réseau, orientation du réseau par rapport à l'axe de la trajectoire) spécifiques à la série de mesures correspondante avant de se placer dans des conditions homogènes au moment du regroupement des mesures. Avantageusement, Les repères de trajectoire et d'objet ont un axe de même direction, ce qui est toujours possible dans Le cas o il n'y a que deux trajectoires, ou, encore, les repères de trajectoire ont un axe confondu avec l'axe de La trajectoire circu- laire correspondante; de préférence, l'origine des repères de trajectoire est confondue avec le centre du cercle. Il est encore avantageux que les groupes de coordonnées dans les repères de trajectoire et le repère d'objet soient des coordonnées sphériques et que le calcul de changement de repère comporte une conversion temporaire des coordonnées sphériques en coordonnées cartésiennes. Il est encore possible que le repère d'objet coîncide avec un des repères de trajectoire. Un cas particulier important est celui de trajectoires concentriques. Le repère d'objet aura de préférence un axe suivant une direction du maillage décrivant L'image de l'objet. Le dispositif selon l'invention comprend au moins un réseau bidimensionnel de capteurs, les réseaux étant déplacés autour d'un objet à examiner suivant au moins deux trajectoires circulaires et reliés à des moyens d'acquisition de mesures et de traitement des mesures appropriés à la mise en oeuvre du procédé de l'invention. On va maintenant décrire l'invention à l'aide des figures suivantes annexées à titre illustratif et non limitatif: la figure 1, déjà décrite, résume les conditions géométriques à La source de l'invention; la figure 2 représente un dispositif per- mettant la mise en oeuvre pratique de l'invention; la figure 3 illustre les moyens géomé- triques mis en oeuvre; la figure 4 est une explication d'une étape d'interpolation; et la figure 5 représente un autre dispositif possible. Dans le cas o l'objet à examiner est un corps humain irradié par une source externe, le patient peut être allongé sur une table 1 Un rail circulaire 2 entoure la table 1 et peut être mis en rotation autour d'un axe vertical à l'aide de deux charnières 3 qui le relient à des points fixes de l'installation. Le rail circulaire 2 est en fait une crémaillère qui porte en des positions diamétralement opposées une source 4 émettrice d'un faisceau conique 6 de rayons X et un écran 5 porteur d'un réseau bidimensionnel de capteurs La source 4 et l'écran 5 sont déplacés le long du rail 2 de conserve sous l'action de moteurs non représentés commandés automatiquement et qui actionnent des pignons engrenant avec la crémaillère. Les trajectoires circulaires différentes sont obtenues par la rotation du rail 2 Il est encore possible d'avoir un rail fixe et une table 1 mobile entre des positions déterminées à l'avance Il existe bien d'autres dispositifs qui peuvent être utilisés, et dont quelques-uns sont décrits dans la demande de brevet antérieure. La figure 3 représente le mode de mise en oeuvre du procédé dans le cas le plus général o deux trajectoires circulaires T 1 et T 2 ont des axes non concourants et non parallèles et dont les directions sont quelconques AB désigne le segment le plus court pour joindre ces axes Les repères de trajectoire R 1 et R 2 ont pour origine les centres 01 et 02 des trajectoires T 1 et T 2, les axes Z 1 et Z 2 sont perpendi- culaires aux trajectoires respectives et les axes Y 1 et Y 2 sont parallèles On définit en outre un repère d'objet R 3 dans lequel les reconstructions de l'image de l'objet sont effectuées L'axe Y 3 de ce repère est parallèle aux axes Y 1 et Y 2 pour faciliter les calculs Aucune autre condition n'est requise pour que les calculs qui vont suivre soient applicables. En pratique, les origines 01, 02 et 03 seront générale- ment confondues et le repère R 3 sera identique à un des autres repères R 1 ou R 2, mais comme ces cas simpli- ficateurs ne sont pas considérés ici par souci de généralité, il faut donc en outre définir deux repères intermédiaires R'1 et R'2 d'origine 03 et dont les axes sont parallèles à ceux des repères de trajectoire R 1 et R 2 respectivement On désigne par e 1 et f 2 les angles entre l'axe Z 3 d'une part, les axes Z 1 ou Z'1 et Z 2 ou Z'2 d'autre part. Soit maintenant un plan de Radon P sur lequel le calcul de somme a été effectué et pour lequel un changement de repère est nécessaire pour effectuer le regroupement ou la fusion des séries de mesures. Les projections orthogonales des trois origines 01, 02 et 03 sur le plan P sont désignées par M 1, M 2 et M 3 Le point M 3 est le point caractéristique du plan de Radon dans le repère d'objet R 3 et donc essentiel pour les calculs d'inversion de la transformée de Radon qui mènent à la reconstruction de l'image de l'objet Il n'a cependant aucune signification dans les repères de trajectoire R 1 et R 2 o le même plan de Radon P est défini par les points caractéristiques M 1 ou M 2. Les points M 1 à M 3 sont définis généralement par leurs coordonnées sphériques car ces coordonnées sont les plus pratiques pour définir des maillages réguliers sur des volumes mesurés à l'aide de trajec- toires circulaires ainsi que pour définir les points caractéristiques Le passage en coordonnées carté- siennes n'est en principe effectué qu'à l'étape finale, pour définir des images de l'objet par des plans de coupe Une conversion temporaire en coordonnées carté- siennes est cependant effectuée dans le procédé décrit ici Si i est l'indice d'un des repères R 1, R 2, ou R 3, on dispose des formules suivantes: Xi = pi sin( Oi) cos(;i) Yi = p sin( Oi) sin(;i) Zi = Pl cos(Oi) pour assurer la correspondance avec des coordonnées cartésiennes qui seront utilisées dans les calculs de changement de repère. Les formules suivantes: X 1 = X'1 ( 1 + (cxl X'I+cyl Y'I+czl Z'I)) X '12 +y'12 +z '12 Y 1 = Y'1 ( 1 + (cxl X'1 +cyl Y'1 +czl Z'1) X'12 +y'12 +z 12 Z 1 = Z'1 ( 1 + (cxl X'1 +cyl Y'1 +czl Z')) X'1 2 +y '12 +Z '12 XI 1 = X 3 cos(C 1) + Z 3 sin (^ 1)Y Y 1 = Y 3 Z'1 = -X 3 sin(p 1) + Z 3 cos(fl) o cxl, cyl et czl sont les coordonnées cartésiennes de l'origine 03 dans le repère R 1 de trajectoire, permettent de retrouver les coordonnées X 1, Y 1 et 1 1 Z 1 du point M 1dans le repère de trajectoire R 1 en fonction des coordonnées X 3, Y 3 et Z 3 du point M 3 dans le repère d'objet R 3 Des transformations inverses seraient également possibles Le calcul pour le point M 2 est identique à un changement d'indice près. On calcule ensuite les coordonnées sphériques 01, Pl et 41 du point M 1 dans le repère R 1 par les formules suivantes, o "atan" signifie arc tangente: si Z 1 > O: 01 = atan(IX 12 + Y 12 / Z 1) Pl = VX 12 + Y 12 + Z 12 si Z 1 O: si Y 1 > O: gl = acos(X 1 / /X 12 + Y 12) si Y 1 0: (Pl = O si X 1 + 1 = T si X 1 = O: 1 = O. si Z 1 4 O si Y 1 > O: I 1 = Ir+ acos(X 1/ X 12 + Y 12) si Y 1 O: 01 = r si X 1 ( 1 = T si X 1 = O: 1 = O. Cependant, le point M 1 ainsi obtenu n'est normalement pas un point d'échantillonnage du volume de Radon correspondant à la première trajectoire T 1, mais un point intermédiaire entre les points du maillage retenu On effectue alors une interpolation du maillage dans un volume pour affecter une valeur au point M 1 (ou Mi dans le cas général) à partir de huit points environnants Mil à Mi 8, dont les valeurs de la somme de la fonction sur le plan de Radon associé sont Vi,j Les coefficients d'interpolation sont dpi, d Oi et d si Pi,j a le même signe que P 3; sinon, Vci,j = -Vi,j. On rappelle que les coordonnées sphériques sont défi- nies avec -oo Pour chacun des points du volume caractéris- tique de l'objet défini dans le repère d'objet, le regroupement consiste, par exemp Le sur un point de L'intersection des vo Lumes caractéristiques de chacune des trajectoires, à effectuer une moyenne des sommes de fonction associées à chaque trajectoire ou, pour un point appartenant exclusivement à l'un des volumes caractéristiques des trajectoires, à considérer La somme de La fonction calculée dans Le repère de cette trajectoire, ces volumes caractéristiques étant définis par rapport à L'origine du repère de l'objet. Les ca Lcu Ls d'inversion, conformément à la demande de brevet antérieur par exemple, sont effec- tués quand des sommes de valeurs ont été réunies pour tous les points correspondant au volume caractéristique de l'objet Les mêmes étapes de filtrage, de pondéra- tion et de rétroprojection sont effectuées. Le procédé peut aussi être employé avec des modélisations différentes ou des repères diffé- rents, et on peut utiliser plus de deux trajectoires circulaires. Par ail Leurs, on a considéré avant tout L'utilisation de ce procédé pour l'inversion de la transformée de Radon proprement dite ou de sa dérivée première On peut encore procéder sensiblement de La même façon avec une transformée de Fourier de tous Les rayons de l'espace de Radon qui ne coupent pas La zone d'ombre, et ceci pour chacun des volumes carac- téristiques de mesure La synthèse des mesures ainsi transformées dans le troisième repère peut être menée sans changement pour peu que les repères R 1, R 2 et R 3 aient leurs centres confondus. Au Lieu de la transformée de Radon ou de ses dérivées premières, on peut encore procéder de la même façon avec la transformée de Hilbert de ces transformées. Le filtrage de la dérivée première pourrait également être effectué avant de réaliser L'inter- polation de la figure 4, ce qui permettrait de ne pas déterminer les orientations Par contre, il fau- drait effectuer le filtrage des dérivées une fois pour chaque volume R 1 ou R 2. Si les mesures prises autour des trajectoires T 1 et T 2 sont en nombre insuffisant, plutôt que de recourir à une interpolation de réarrangement, exposée dans la demande de brevet précédente, o les plans de Radon représentatifs des points M 1 ou M 2 ne sont pas mesurés directement mais à partir d'interpolations réa Lisées sur des mesures correspondant à deux posi- tions voisines du point focal S, il est préférable d'incorporer l'opération de réarrangement dans la phase de regroupement décrite dans la présente demande, afin de n'avoir à réaliser qu'une seule interpolation. Comme il a déjà été signa Lé, l'invention peut aussi être utilisée en tomographie d'émission. On peut alors plus facilement accepter des trajectoires circulaires incomplètes Comme la figure 5 le montre, un autre dispositif comporte un premier rail 21 entou- rant la table 1 et appartenant à un plan transversal par rapport au patient Un second rail 22 est situé à une extrémité de la table 1 et appartient à un plan longitudinal du patient; contrairement au précédent, il est interrompu pour livrer passage à la table 1 et au patient Une telle disposition est intéressante par exemple pour l'examen du cerveau. Un écran 51 parcourt une trajectoire T 1 en cercle complet autour du rail 21 et un écran 52 parcourt une trajectoire T 2 en arc de cercle autour du rail 22 Comme toutefois il n'y a pas de source en tomographie d'émission, les points de la trajectoire T 2 opposés à la partie manquante du rail 22 peuvent servir de support à des mesures. L'invention peut servir pour les applications usuelles d'imagerie médicale ou de contrôle non des- tructif de pièces. REVENDICATIONS 1 Procédé de reconstruction d'images tri- dimensionnelles d'un objet (M) défini par des valeurs prises par une fonction sur des points de l'objet, la fonction étant une propriété d'un rayonnement conique ayant un point focal (S) et passant à travers l'objet (M), dans lequel la fonction est calculée par l'intermédiaire de sommes de la fonction sur des plans (P) passant par au moins un point de l'objet et définis dans un repère d'objet (R 3), comportant au moins deux séries de mesures, chaque somme de la fonction étant calculée, sur lesdits plans, à partir d'au moins une des séries de mesures, chacune des séries de mesures étant menée avec un réseau bidimen- sionnel de capteurs ( 5) de rayonnement orientés vers le point focal, qui est déplacé autour de l'objet sur une trajectoire circulaire (T 1, T 2) respective, caractérisé en ce que, pour chaque série de mesures, une série des sommes de la fonction est calculée sur certains au moins desdits plans, qui sont sécants ou tangents à la trajectoire associée à ladite série de mesures et sont alors définis dans un repère de trajectoire associé à ladite trajectoire, et en ce que les séries de sommes de la fonction sont regroupées en effectuant des calculs pour pouvoir définir les plans dans le repère d'objet. 2 Procédé de reconstruction d'images tri- dimensionnelles suivant la revendication 1, caractérisé en ce que les repères (R 1, R 2, R 3) de trajectoire et d'objet ont un axe de même direction (Y 1, Y 2, Y 3). 3 Procédé de reconstruction d'images tri- dimensionnelles suivant l'une quelconque des revendica- tions 1 ou 2, caractérisé en ce que les repères de trajectoire ont un axe (Z 1, Z 2) perpendiculaire à la trajectoire respective (T 1, T 2). 4 Procédé de reconstruction d'images tri- dimensionnelles suivant l'une quelconque des revendica- tions I à 3, caractérisé en ce que les groupes de coordonnées dans les repères de trajectoire et le repère d'objet sont des coordonnées sphériques, et en ce que les calculs de changement de repère comportent une conversion temporaire des coordonnées sphériques (p,,) en coordonnées cartésiennes (X, Y, Z). 5 Procédé de reconstruction d'images tri- dimensionnelles suivant l'une quelconque des revendica- tions 1 à 4, caractérisé en ce que le repère d'objet coïncide avec un des repères de trajectoire. 6 Procédé de reconstruction d'images tri- dimensionnelles suivant l'une quelconque des revendica- tions 1 à 5, caractérisé en ce que les trajectoires sont concentriques. 7 Procédé de reconstruction d'images tri- dimensionnelles suivant l'une quelconque des revendica- tions 1 à 6, caractérisé en ce que les plans sont définis dans les repères de trajectoire et d'objet par des points caractéristiques appartenant à ces plans, les points caractéristiques associés aux repères de trajectoire et d'objet étant différents, et en ce que les opérations de regroupement consistent à faire une conversion entre les coordonnées dans le repère de trajectoire (X 1, Y 1, Z 1) du point caractéris- tique (M 1) associé audit repère de trajectoire (R 1) et les coordonnées (X 3, Y 3, Z 3) dans le repère d'objet (R 3) du point caractéristique (M 3) associé au repère d'objet. 8 Procédé de reconstruction d'images tri- dimensionnelles suivant l'une quelconque des revendica- tions 1 à 7, caractérisé en ce que le rayonnement est un rayonnement d'atténuation émis par une source ( 4) ponctuelle au point focal, tournant autour de l'objet suivant une trajectoire circulaire de source. 9 Procédé de reconstruction d'images tri- dimensionnelles suivant l'une quelconque des revendica- tions 1 à 7, caractérisé en ce que le rayonnement est un rayonnement d'émission provenant de l'objet, le point focal étant un point vers lequel convergent les directions de mesure des capteurs. Procédé de reconstruction d'images tri- dimensionnelles suivant l'une quelconque des revendica- tions 1 à 9, caractérisé en ce que le repère d'objet a un axe correspondant à une direction du maillage décrivant l'image tridimensionnelle de l'objet. 11 Dispositif comprenant au moins un réseau bidimensionnel de capteurs, les réseaux étant déplacés autour d'un objet à examiner suivant au moins deux trajectoires circulaires et reliés à des moyens d'ac- quisition de mesures et de traitement des mesures appropriés à la mise en oeuvre d'un procédé conforme à l'une quelconque des revendications précédentes.