Brevet FR2648304A1 — publié le 19901214

PROCEDE DE DETERMINATION D'UN ESPACE A PARTIR D'UN ESPACE DISCRET CONNU POUR LA RECONSTRUCTION D'IMAGES BI OU TRIDIMENSIONNELLES, DISPOSITIF DE MISE EN OEUVRE ET APPLICATION DU PROCEDE DESCRIPTION La présente invention a pour objet un procédé de détermination d'un espace à partir d'un espace discret connu pour la reconstruction d'images bi ou tridimensionnelles. Elle a également pour objet un dispositif de mise en oeuvre et une application du procédé. Dans des applications telles que l'imagerie médicale permettant de reconstituer des images du corps humain ou d'imagerie industrielle pour contrôle non destructif de pièces l'on fait appel à des techniques de reconstruction d'images. Ces reconstructions d'images se font à partir de projections obtenues par un dispositif d'acquisition émettant et captant un rayonnement d'irradiation. Ces projections sont obtenues à la suite d'une série de mesures de l'atténuation du rayonnement à travers l'objet a contrôler, l'incidence du rayonnement avant été modifiée entre chaque série de mesures. Les dispositifs d'acquisition sont constitués d'une source de rayonnement X ou r ou à neutrons. Ils comportent également un ou plusieurs détecteurs. Lorsqu'il n'y a qu'un seul détecteur, la source se déplace linéairement et émet à chaque déplacement pour obtenir une projection à partir d'une série de mesures; la géométrie du dispositif d'acquisition est dite parallèle. Lorsqu'il y a une barrette de détecteurs, la source est fixe, la série de mesure est obtenue à partir d'une seule émission, la géométrie d'acquisition est dite en éventail. En présence d'une matrice de détecteurs, la géométrie est dite conique. Généralement, l'on cherche à avoir un grand nombre de projections afin d'obtenir la meilleure reconstitution possible d'images de L'objet ainsi traversé par Le rayonnement. Pour ceLa, l'on fait tourner la source par petits incréments d'angles autour de l'objet afin d'obtenir des séries de mesures sous incidences différentes. D'une manière générale, la reconstruction d'images (images bidimensionnelles) ou de volumes (images tridimensionnelles) à partir de projections fait appel à deux grandes classes de méthodes: - la classe des méthodes analytiques, - la classe des méthodes algébriques. L'analyse de ces méthodes analytiques et algébriques montre que, outre les opérations de filtrage, réarrangement, masquage, etc..., les opérations de base sont La reprojection et La rétroprojection. - La re2roiection (ou projection) est une opération qui permet d'obtenir, à partir d'une image ou d'un volume, les projections relatives à ces espaces dans une géométrie donnée. Dans le cas de l'imagerie X, chaque point de La projection correspond à une estimation de l'atténuation énergétique d'un faisceau de rayons X sur un trajet linéaire dans l'espace considéré supposé connu. Cette opération de reprojection est utilisée en particulier dans les techniques itératives (méthode algébrique) et dans certaines méthodes de reconstruction de volumes par synthèse et inversion de la transformée de RADON. L'on pourra se référer à la thèse de P. Grangeat soutenue à l'Ecole Nationale Supérieure des Télécommunications le 30 juin 1987 ayant pour titre "Analyse d'un système d'imagerie 3D par reconstruction à partir de radiographies X en géométrie conique". - La Erétroprojecion (ou épandage) est le traitement dual de la reprojection. Elle consiste à redistribuer l'information issue des projections aux pixels images (ou aux voxels volumes) appartenant à une droite d'acquisition. Cette opération est utilisée dans l'ensemble des méthodes de reconstruction (images ou volumes) à partir de projections. Elle est parfois associée, dans le cadre des méthodes analytiques en particulier, à des opérations de filtrage et de réarrangement de données. Ces deux opérations (reprojection et rétroprojection) sont étroitement liées avec la géométrie du système d'acquisition ayant engendré les données initiales. Cette géométrie comme cela a été dit précédemment peut itre: - parallèle, - en éventail (pas linéaire constant ou pas angulaire constant), - conique (en 3-D). Le volume de données a traiter, et la complexité des procédés de reconstruction, en font un problème nécessitant la mise en oeuvre de moyens de traitement performants. La solution naturelle qui consiste à utiliser des processeurs vectoriels non dédiés (standards) pour mettre en oeuvre ces procédés conduit à un rapport prix/performance élevé. Cette solution assure une certaine souplesse (car ils sont programmables souvent en language évolué), mais ceci au détriment de la performance et du prix. Par aittlleurs, il existe à l'heure actuelle des développements de machines spécifiques. Des processeurs de rétroprojection bidimensionneLle spécifiques ont été réaLisés dans le cadre de systèmes tomographiques. Par exemple, la société STAR Technologie diffuse un rétroprojecteur (géométrie éventail) associé à un processeur vectoriel. L'ensemble des produits existant sur le marché sont, de par leur conception, liés à une géométrie d'acquisition, et n'ont pas en conséquence été conçus pour permettre une mise en oeuvre efficace des opérations de reprojection et rétroprojection, dans les différentes géométries. L'on rappelle que les procédés permettant de réaliser des reprojections et des rétroprojections se caractérisent par un ensemble de phases itératives imbriquées: phase sur les projections -phase sur les plans images uniquement dans le cas 3D -phase sur les lignes-images --phase sur les éléments de la ligne- la phase sur les projection peut être imbriquée n'importe o avant la phase sur les éléments. La phase de plus bas niveau porte sur La correspondance des éléments de ligne et est qualifiée de phase klémentaire. Cette phase met en correspondance deux ensembles monodimensionnels (que L'on peut assimiler à des lignes) dont l'un est connu. Les opérations internes à cette phase permettent de calculer une contribution à chaque élément de l'ensemble monodimensionnel inconnu. Cette phase est répétée de 3 4. l'ordre de N (respectivement de l'ordre de N fois) lors de la reconstruction d'une image de taille N (respectivement d'un volume de taille N) à partir de N projections. La présente invention consiste en un procédé et en un dispositif de détermination d'un espace à partir d'un espace connu discret pour la reconstruction d'image ou de volume permettant de remédier aux inconvénients apportés par les méthodes et les dispositifs usuels de reconstructions d'images. En effet le procédé conforme à l'invention permet d'être utilisé quelle que soit la géométrie du dispositif d'acquisition, parallèle 2D, 30, éventail 2D, conique 3D. Le dispositif de mise en oeuvre du procédé est conçu pour être utilisé quel que soit le dispositif d'acquisition contrairement à l'art antérieur. La présente invention a donc pour objet un procédé de détermination d'un espace à partir d'un espace discret connu par projection ou rétroprojection, pour la reconstruction d'une image bi-dimensionnelle ou tridimensionnelle d'un corps à analyser, à partir d'un ensemble de projections obtenu, au moyen d'un dispositif d'acquisition de géométrie définie comportant une source de rayonnement et un ensemble de détecteurs, l'un des espaces étant défini par un ensemble de points élémentaires repérés respectivement suivant trois axes par des indices (l, $, n) correspondant aux coordonnées suivant ces axes; l'autre espace étant défini par un ensemble de points élémentaires que l'on repère également respectivement par des indices (i,j,k) suivant trois axes de cet espace, ces indices (i, j, k) correspondant aux coordonnées suivant ces axes, la détermination de l'espace inconnu est obtenue par le calcul de 'a valeur numérique de chaque point de l'espace inconnu à partir. de la détermination des contributiorsapportées par les valeurs numériques de points correspondants de l'espace connu, le procédé étant caractérisé en ce qu'il comporte une phase élémentaire consistant à déterminer de façon itérative pour chaque point d'un ensemble monodimensionnel de points de l'espace inconnu portés par un axe E donné de cet espace, la contribution apportée par le point correspondant et éventuellement par les points voisins d'un ensemble monodimensionnel de points de l'espace connu portés par un axe F de cet espace, cette phase comportant de façon itérative pour chaque point les étapes suivantes: - déterminer l'indice associé à La coordonnée sur l'axe F du point correspondant de l'espace connu, en prenant la partie entière ADR de l'expression généralisée a (e) avec V1 + eV2 a (e) =-à------- 0 V3 + eV4 e étant l'indice d'itération correspondant à la coordonnée sur l'axe E du point considéré de l'espace inconnu et V1, V2, V3, V4 étant des paramètres constants dans la phase élémentaire, l'ensemble de ces paramètres étant déterminé préalablement à cette phase en fonction de la géométrie d'acquisition choisie; du type d'opération (reprojection ou rétroprojection) effectuée et des ensembles monodimensionnels considérés, - déterminer la contribution C à apporter à la valeur numérique du poinrt considéré de l'espace inconnu par interpolation sur les valeurs numériques du point correspondant d'indice ADR et éventuellement des points voisins, - pondérer cette contribution par un facteur f déterminé par la relation suivante: V5 f = --- - -- - (V6 + eV7) V5, V6, V7 étant des paramètres constants dans la phase élémentaire, l'ensemble de ces paramètres étant déterminé préalablement à cette phase en fonction de la géométrie d'acquisition choisie, du type d'opération et des ensembles monodimensionnels considérés, - additionner cette contribution pondérée à la valeur numérique du point considéré de l'espace inconnu, déterminee lors de phases élémentaires précédentes ayant mis en correspondance le point de l'espace inconnu à des points d'autres espaces monodimensionnels connus, l'ensemble des contributions obtenues à la fin du procédé de détermination de l'espace inconnu constituant la valeur finale de ce point. Selon une première variante du procédé de détermination d'un espace à partir d'un espace discret connu, la phase élémentaire représente une opération de rétroprojection en géométrie parallèle, l'espace connu étant un ensemble de données de projection indicé par l,m,n; m et n étant fixes, l'espace inconnu étant l'espace image indicé par i, j, k, (j et k) ou (i et k) étant fixés; l'indice recherché a (e) est soit l'indice l (i) lorsque i et k sont fixes, l'axe E étant O indicé par i et l'axe F par l, soit l'indice l (j) lorsque j et k sont fixes, l'axe E étant indicé par j et l'axe F par l; et l'on choisit alors: V3 = 1 V4 = 0 et pour i fixé à i k=O, m=O, n=n, (j) étant 0O' recherché Vl = L + i AI O O V2 = aJ À = j e= j ou pour j fixe j, k=0, m=O, n=n, L(i) étant 0 0nno 0 recherché V = L + j a.J O O V2 = AI e? = j, ei les termes LO, AI, J étant prédéterminés, l étant un nombre entier réel de partie entière l. Selon une deuxième variante de réalisation du procédé de détermination d'un espace à partir d'un espace discret connu, la phase élémentaire représente une opération de reprojection en géométrie parallèle, l'espace connu étant l'espace image indicé par L, m, n; m et n ou l et m étant fixés, l'espace inconnu étant l'ensemble de données de projection indicée i, j, k ou j et k sont fixes; l'indice a (e) recherché est l'indice lo(i) ou n o(i), l'axe E est alors indicé par i ou l'axe F est indicé par l ou n et l'on choisit alors: 1 5 m fixé à 0, n à no, j à jO k à O, on cherche l (i) o V3 = 1 V4 = 0 V1 = h (n, j) 2 0 jo) V2 = AL (j I 0 e = i ou m fixé à 0, h à l, j à jo, k à 0; on cherche n (i) V3 = 1 V4 = 0 Vl = h (l, j) 2 0 O) V2 = a L (J) e = i 2 Les termes h1, h, eL et L étant 1 2 I 2 prédéterminés et n et l sont des nombres réels de 0 O partie entière respectivement n et l. Selon une troisième variante de réalisation du procédé de détermination d'un espace à partir d'un espace discret connu, la phase élémentaire représente une opération de rétroprojection en géométrie en éventail à pas constant, l'espace connu est l'ensemble de données de projection indicees par l,m,n; m et n étant fixés, l'espace inconnu étant l'espace image indice par i, j, k, (i et k) ou (j et k) étant fixes, l'indice recherché a (e) est soit l'indice lo( i) soit l (j) l'axe E étant indicé par i ou par j et l'axe F par l et l'on choisit alors: k=O, m=O, n=n0, j= Jo' lot(i) étant recherché, Vl = h (j, n) 3 O' 0 V2 = h (n 4 0 V3 = h (j n Jo' no V4 = h (n) e = i 6 O ou k=O, m=O, n=n, i=i l (j) étant recherché, Vl = h'3 n V2 = h' (n) 4 0 V3 = h' (i n 0 V4 = h' (n) 6 0 À0 = j. les termes h, h, h5 h, h', h'4, h' h' étant prédéterminés et l est un nombre réel de partie entière L. Selon une quatrième variante de réalisation du procédé de détermination d'un espace à partir d'un espace discret connu, la phase élémentaire représente une opération de reprojection en géométrie en éventait, l'espace connu étant l'espace image indicé par L, m, n; m et n ou L et m étant fixes, l'espace inconnu étant l'ensemble de données de projection indicees par i, j, k ou j et k sont fixes, l'indice recherché a (e) est l'indice l (i) ou n (i) l'axe E étant indice par i et O O l'axe F par L ou n et l'on choisit alors: m fixe à 0, n à n, j à j, k à 0, l () étant recherche, Vl = h (n, j)j 7 0 0 V2 = h (nO, j) 9 0 V4 = h (j) 0 e = i ou:m fixé à 0, l à lo, j à k à O, 0 étant recherché, Vl = h' (lo j) V2 = h' (L O j) 8 0 V3 =h' (j) 9 0 V3 = h' (j) 0 e = i Les termes h h8, h, h h', h'8 ' 7 h8, h9, 1 O' h' h 9' h' étant prédéterminés et n, l étant des nombres O' réels de partie entière respectivement n et I. Selon une cinquième variante de réalisation du procédé de détermination d'un espace à partir d'un espace discret connu, la phase élémentaire représente une opération de rétroprojection en géométrie conique; l'espace inconnu étant un espace image indicé par i, j, k (j et k étant fixés), l'espace connu étant l'espace de projection indicé par l, m, n (n étant fixé et m fixé puiscue fonction de j et k), l'indice recherché a (e) est l'indice l Mi) l'axe E est indicé par i et O 0 l'axe F par l, et l'on choisit alors: V3 = 1 V4 = 0 C-Qo Vl = Zo tbQ C.Z V2 = les termes Z0, C, Q, d Z et 4Q étant prédéterminés et l étant un nombre réel dont la partie entière est l. O Selon une sixième variante de réalisation du procédé de détermination d'un espace à partir n espace discret connu, la phase élémentaire repres. :e i1 une reprojection en géométrie conique, l'espace inconnu étant l'espace de projection indicé par i, j, k, (j et k étant fixes), et l'espace connu étant l'espace image indicé par l, m, n (m étant fixé, n étant fixé puisque fonction de j et k); l'indice recherché a (e) est l'indice l Ci) l'axe E est indicé par i et l'axe F par L et L'on choisit aLors: V3 = 1 V4 = O Vi = Qc=zo 4z CaQ V2 = --- ,Z w les termes Z0, Q0, C, AZ et.Q étant prédéterminés et l étant un nombre réel dont la parie entière est l. Selon une septième variante de réalisation du procédé de détermination d'un *espace & partir d'un espace discret connu, pour une géométrie paralLèle et conique dans Le cas ou La phase élémentaire représente 2 une rétroprojection et pour une géométrie parallèle lorsque cette phase représente une reprojection, le facteur de pondération f est déterminé en prenant: V7 = O V6 = 1 et V5 est prédéterminé. Selon une huitième variante de réalisation du procédé de détermination d'un espace à partir d'un espace discret connu dans l'étape de pondération pour une géométrie éventail dans le cas d'une rétroprojection, le facteur de pondération f est déterminé en prenant: V5 prédéterminé à partir des paramètres du dispositif d'acquisition V6 = V3 V7 = V4 Seon une neuvième variante de réalisation du Selon une neuvième variante de réalisation du procédé de détermination d'un espace à partir d'un espace discret connu, dans l'étape de pondération pour une géométrie en éventail et conique dans le cas de la reprojection, le facteur de pondération f est déterminé en prenant: V5 = 1 V6 = 1 V? = 0 L'invention consiste en outre en un dispositif de mise en oeuvre du procédé caractérisé en ce Qu'il comporte outre le dispositif d'acquisition: - des moyens de génération d'adresse recevant les paramètres V1-V4, délivrant l'adresse ADR qu est la partie entière de a (e); O - des moyens de mémorisation des valeurs numériques des points de l'espace connu, adressées par les moyens de génération, dé.livrant l'ensemble de " +1) valeurs numériques des points de l'espace connu de V(ADDR-b) à V(ADR+d) (f et îentiers positifs; + =) correspondant a l'indice ADR, i étant l'ordre d'interpolation; - des moyens d'interpolation recevant les valeurs numériques V(ADR-I) à V(ADR+ Y) continu- - mémoire et adressables par l'indice ADR et déterminant a la fin de l'interpolation une contribution non pondérée à la valeur numérique de chaque point correspondant de l'espace inconnu; - des moyens de pondération recevant chaque contribution non pondérée et les paramètres V5, V6, V7 pour délivrer chaque contribution pondérée-; - des moyens de sommation pour additionner à la valeur numérique du point en cours de l'espace inconnu, la contribution pondérée relative à ce point en cours; des moyens de mémorisation contenant la valeur en cours. L'invention a également pour objet une application du procédé à la reconstruction par synthèse et inversion de La transformée de Radon. La reconstruction par synthèse et inversion de- la transformée comporte selon le procédé trois étapes: - une première étape consistant en une reprojection bidimensionnelle en géométrie parallèle permettant d'obtenir une transformée de Radon intermédiaire; - une deuxième étape consistant en une rétroprojection bidimensionnelle en géométrie parallèle permettant d'obtenir des projections réarrangées; - une troisième étape consistant en une rétroprojection bidimensionnelle à géométrie parallèle permettant d'obtenir le volume final recherché. D'autres caractéristiques et avantages de l'invention ressortiront mieux à la lecture de la description qui va suivre, donnée à titre illustratif et non limitatif. La description se réfère aux figures annexées dans lesquelles: - la figure I représente les étapes principales du procédé sous forme d'une succession de blocs, chaque bloc correspondant à une étape particulière, conformément à l'invention, - la figure 2 représente une figure géométrique permettant de rappeler le principe d'acquisition 2D en géométrie parallèle, - la figure 3 représente une figure géométrique permettant de rappeler le principe d'acquisition 2D en géométrie éventail; - la figure 4 représente une figure géométrique permettant de rappeler le principe d'acquisition 3D en géométrie conique, - la figure 5 représente le cas d'une rétroprojection avec un dispositif d'acquisition à géométrie parallèle, - la figure 6 représente une reprojection avec un dispositif à géométrie parallèle, - la figure 7 représente une rétroprojection avec un dispositif à géométrie en éventail, - la figure 8 représente une reprojection avec un dispositif à géométrie en éventail, - la figure 9 représente le cas de rétroprojection et reprojection 3D en géométrie parallèle, - la figure 10 représente le schéma d'adressage dans la phase élémentaire pour une rétroprojection 3D en géométrie conique, - la figure 11 représente le schéma du dispositif de mise en oeuvre du procéde conforme à l'invention, - la figure 12 montre l'espace volume à reconstruire de coordonnées x, y, z et l'espace projectin connu de coordonnées p, q, e, - la figure 13 montre l'espace de Radon de coordonnées p, e, y et le repère O, i,], k Lié à l'espace du volume à reconstruire, - la figure 14 montre l'espace des projections réarrangées de coordonnées r, z,; - la figure 15 montre une reprojection parallèle d'un plan radiographique, - la figure 16 montre une rétroprojection paralèle d'une partie d'un volume, - la figure 17 montre une rétroprojection parallèle d'une partie du volume de rétroprojection réarangées. La figure 1 permet de montrer sous forme d'une succession de blocs, les différentes étapes du procédé selon l'invention. Tout d'abord l'on rappelle quelques définitions. L'on entend par espace discret connu un espace bidimensionnel (2D) ou tridimensionnel (3D) numérisé, c'est-à-dire un espace défini par un ensemble de points élémentaires correspondant à des pixels en 2D ou à des voxels en 3D. Les coordonnées des points sont numérisées dans l'espace image ou volume et dans l'espace projection. Ces points sont référencés par leurs indices i, j, k dans l'espace inconnu, l, m, n dans l'espace connu. L'on rappelle également que dans l'opération de reprojection, l'espace image ou volume est connu et l'espace projection est à calculer. Dans l'opération de rétroprojection, l'espace projection est connu et l'espace image ou volume est à calculer. Selon le procédé objet de l'invention, la phase élémentaire telle qu'elle a été précédemment définie, comporte un ensemble d'étapes qui vont permettre de déterminer la contribution à la valeur d'un point de l'espace inconnu à partir de la valeur du point correspondant (et de son voisinage) de l'espace connu, de façon unifiée pour les opérations de reprojection et de rétroprojection, pour un grand nombre de géométriesd'acquisition. Cette phase comporte pour cela une étape de détermination d'un indice ou adresse dans l'espace connu, suivi d'une détermination de la contribution par interpolation sur les données avoisinant cette adresse et enfin une étape de pondération de cette contribution. L'ensemble de ces étapes se renouvelle par itération d'un indice apparaissant dans la détermination de l'adresse jusqu'à ce que cet indice ait atteint sa borne supérieure N, c'est-à-dire Le nombre de points élémentaires caractérisant le nombre de points suivant un axe de L'espace connu (taille de l'espace). Détermination de L'adresse (étape I). D'une manière générale cette étape consiste, en fonction de L'adresse d'un point de l'espace à calculer, à déterminer l'adresse du point correspondant dans l'espace connu. Le calcuL d'adresse est étroitement lié à La géométrie d'acquisition et fait intervenir Les formuLes affines de changement de repère. Cependant conformément à l'invention, La détermination de L'adresse dans la phase élémentaire est réduite à la détermination de la partie entière ADR de la valeur a (e) par la relation suivante: V1+eV2 a (e) = -à- 0 V3+eV4 dans laquelle les paramètres V1, V2, V3, V4 sont prédéterminés en fonction de la géométrie d'acquisition et du type d'opération effectuée (rétroprojection ou reprojection) et des ensembles monodimensionnels considérés. e est l'indice d'itération correspondant à L'adresse,-d'un point de L'espace inconnu. Détermination de la contribution C (étape II). Dans cette étape, L'on utilise Les points appartenant à l'espace connu, avoisinant l'adresse ADR précédemment calculée pour effectuer une interpolation. L'interpolation se fait à partir de procédés classiques. Il peut s'agir d'interpolation d'ordre zéro ou du plus proche voisin ou d'ordre supérieur. Pour une interpolation d'ordre supérieur, pour une adresse fournie par l'étape de calcut d'adresse, plusieurs données de l'espace connu au voisinage de cette adresse doivent être extraites. Les coefficients d'interpolation utilisés dans cette étape sont déterminés généralement à partir de la partie fractionnaire de a (e). Pondération de la contribution (étape III) La pondération consiste à appliquer un facteur de pondération f dont les paramètres V5, V6, V7 sont prédéterminés et fonction de la géométrie d'acquisition du type d'opérration et des ensembles monodimensionnels considérés. Addition de la contribution ondérée _à l'adres se en cours (étape IV) Cette étape consiste à ajouter la contribution pondérée au point considéré de l'espace à calculer soit à calculer V(e)+f.C, V(e) étant la valeur à l'adresse e du point de l'espace inconnu. Itération de l'indicee Cette dernière étape consiste à incrémenter l'indice e et à recommencer l'ensemble des étapes précédentes jusqu'à atteinte de la valeur N (N correspondant au nombre de points élémentaires sur une ligne de cet espace). La suite de la description va permettre de détailler les différentes étapes selon différentes applications correspondant chacune à un choix particulier de type de géométrie d'acquisition. Tout d'abord l'on rappelle au moyen des figures 2, 3 et 4 respectivement le principe d'une géométrie parallèle 2D, éventail 2D, et conique 3D. S est la source réelle ou virtuelle de rayonnement. Dans les figures 2 et 3, un point M de l'espace objet est caractérisé par ses coordonnées M(x, y") ou M(pM, q), (x, y) dans un repère M( M, YM M M M (0,l,j), et (pM qM) dans un repère (O,u,v) qui est le m M repère lié à la projection d'angle 0. Les points de la projection sur l'axe E sont notes A(p A), p éA tant la coordonnée liée au repère (O',u). Dans Le cas éventail figure 3, le rayon R1 caractérise la distance sourcecentre 0 du repère (O,u,v) et R2 la distance centre du repère au détecteur 0. Les coordonnées sont échantillonnées dans l'espace objet par: x(i) = X +i àX, i et j étant des nombres entiers y(j) = Y +j Y o à X et AY étant les pas d'échantillonnage. Les coordonnées sur l'axe E de projection est donnée par la relation: p(l) = P +l P, l étant un nombre entier o tP étant le pas d'échantillonnage. Un point M est repéré dans le repère 0,i,j par x et y et dans le repère 0, u,v par p et q avec 20. M M M M p = x cos*+y sin$ M M M q = -x sin +y cos$ Mi M M Dans le cas d'une géométrie d'acquisition conique (3D) telle que représentée sur la figure 4, L'ensemble DET formé par les points A relatifs aux points M d'une colonne de volume, forme une colonne dans la radiographie numérique parallèle à l'axe de rotation k et indépendamment de La valeur de l'angle 0. La phase élémentaire permet également dans ce cas de géométrie conique de mettre en relation des ensembles d'éléments (de l'espace connu et de l'espace à reconstruire) appartenant à des droites parallèles à l'axe de rotation t. Le procédé objet de l'invention s'applique tout à fait, à la reprojection et à la rétroprojection en géométrie conique et en géométrie parallèle puisque la phase élémentaire comporte des itérations sur Les éléments de l'espace à reconstruire ne nécessitant qu'une opération pour déterminer l'adresse (dans l'ensemble de l'espace connu) de l'élément correspondant. La figure 5 représente le cas de la rétroprojection pour une géométrieparallèle 2D. Dans ce cas les paramètres V1, V2, V3, V4 servant à la détermination de l'adresse sont précalculés. L'espace image inconnu est échantillonné comme: x(i) = X + i&X o y(i) = Y + jaY z(k) = 0 d'o k=O L'espace projection connu est échantillonné comme: p(l) = P + laP O q(m) = 0 donc m=O e(n) Les paramètres V1, V2, V3, V4 sont précalculés de la manière suivante: m=0, n fixé à n, k=O l est fonction de i j est donc fixé à j O V3 = 1 V4 = 0 e=i e = i Vl = L + j J O iO V2= =I ou m=O, n fixé à nM, k=O l est fonction de j i est fixé à i V3 = 10 V3 =1 V4 = 0 e = j Vl = L + i AI 0 0 V2 = AJ et avec tI = A'.cos e(n) aJ =.in (nO) dp 0 X0 yr P0 L = --.cos e(n) + --. sin(n) Op 0 j AP 0 AP X0, Yo, &P, 0P, sin C(no), cos *(n) A x, &Y étant o O connus. Dans le cas de la figure 5, l'espace de projection est indicée par L, m=0, et n=nO L'espace inconnu est l'espace image indicé i, j, k; k=0, i ou j sont fixés; l'adresse recherchée est l (i si j est fixe à j0' ou l 0(j) si i est fixé à i. La figure 6 représente le cas de la reprojection pour une géométrie parattllle 2D. Dans ce cas les coordonnées sont échantillonnées de la façon suivante: x(l) = X + l.x y(n) = Y + n aY zC(m) = 0 d'une part pour l'espace image connu. et p(i) = P + i P q(k) = 0 d'ou k=O *(i) d'une part pour l'espace projection inconnu. Les paramètres V1, V2, V3 et V4 utiliséspour la détermination sont précalculés de la manière suivante: m=0, k=0, j=j n fixé à n e = i et l ti) est recherchée V3 = 1 V4 = O P sin,(j) X 0 00 vi =---------- - (Y +n Y). ---------- 4xcose(jO) &X.cos(jO) ax &P V2 = ax cos O(j) O ou m=0, k=O, j=j I fixé à l e = i et n (i) est recherchée V3 = i V4 =0 P (X0+L X)cos(j o)) YO V1 = (à--- -) - _ Ysinj Aynsin (j) AY P et V2 = t) OYsinnvO X, Y, &P, P, sin (j0), cos0(j), ax, AY etant connus. Dans le cas de la figure 6, l'espace connu est l'esDace image indicé par l, m, n, m=O, l ou n étant fixés; l'espace inconnu est-L'espace projection indicé par i, j, k, k=0 et j=j0, L'adresse recherchée est L (i) si n est fixé à n, ou n (i) si l est fixé à O O O L La figure 7 représente le cas de la rétroprojection en géométrie éventail. Dans ce cas ces coordonnees des espaces sont échantillonnés comme image, inconnu espace projection connu xx(i)=X +i 4X p(L) = P +l ap O O y 1 O0 i P V2 = -- cX cos (n) + ----- X sin (n a P 0 R1+R2 0 V3 = R1 - X sin 0(n) + y(j) cos 0(n) V4 = - X sin e(n) ou m=O, k=O, n fixé à n, e=j et i fixé à i. On cherche L (j) O V1 = -- Y sin *(n) + x(i)cos*(n)- P O 0 O O O P0 - -- - Rl+Y cos(n)-x(i)sin(n RI+R2 0 0 0 0 P V2 = y--sin(n) ---- ycos (n)} P 0 O R1+R2 0 V3 = R1 - x(i) sin *(n) + Y cos (n) a 0 0 0 V4 = cos *(n) L'ensembLe des coefficients servant à précalculer ces paramètres étant connu. Dans le cas de la figure 7, l'espace de projection connu est indicé par L, m, n. m=O, n=no, l'espace inconnu est l'espace image indicé par i, j, k, k=0, i ou j fixes, l'adresse recherchée est l (i) si j est fixé à j ou l (j) si i est fixé à i O 0 O La figure 8 représente le cas de la reprojection en géométrie en éventail. Dans ce cas les coordonnées des espaces sont échantillonnées de la façon suivante: espace image inconnu x(e) = X + l 4X y(n) = Y + n aY z(n) = 0 donnant n=O espace projection connu p(i) = P + i AP q(k) = 0 donnant k=O + (j) Dans ce cas les paramètres V1, V2, V3 et V4 servant à la détermination de L'adresse sont précalculés de la manière suivante: m=0, k=O, j=J, n=n, e=i O O l (i) est recherché. [R1+y(n) cos +(j - X sin j 0 0 0 Vl = p on R1 + R2 y(n) sin 0(j - X cos 0(j). R +y(n) cos *(jo) X0 sin *(yo)] V2 =. P R1 + R2 p àx PObx - sin 0 (jo) + 6X cos 0iO) RI+R2 &P 4X sin *Cj) O V4 - R1 + R2 ou bien m=0, k=0, j=j e=i, l fixe à l0 n (i) est recherché. vl [R-x(l)sin(i0)-Y0cos(j0) xCL0)cos(j+y0sin AP V2 = ----- (R1 - x(l) sin4(j)+Y cos *(j)) R1+R2 0 0 0 0 P., O cos *(j)) V3 =.Y (sin (j) - ----- 0 0 R1+R2 *&.P V4 = -âY - cos (j) R1+R2 0 l'ensemble des coefficients servant à prétalculer ces paramétres étant connu. Dans le cas de La figure 8, l'espace connu est l'espace image indicé par l, m, n, m=O, l ou n étant fixé; l'espace inconnu. est l'espace projection indicé par i, j, k, k=O j=j; l'adresse recherchée est o I (i) si n est fixé à n ou n (i) à L est fixé à L La figure 9 représente le cas de rétroprojections et reprojections en géométrie parallèle 3D. Un plan d'altitude k se projette toujours o dans la projection d'angle 0 (quel que soit *) à une altitude constante v. La rétroprojections ou la reprojections dans cette géométrie est une succession de rétroprojectionS et reprojections 2D en géométrie parallèle. La détermination de l'adresse dans ce cas est telle que déjà décrite à propos des cas de rétroprojection et reprojection en géométrie paralèle 2D. La figure 10 représente le cas d'une rétroprojection en géométrie conique 3D. Dans ce cas la partie de L'espace à reconstruire figé est pour k=k et jOj qui donne x(k) et y(j) figé. La partie de L'espace connu figé est pour n=n qui fixe 0 à e(n) et m=m' fixé, puisque m est 0 o fonction de k et j. L'adresse recherchée est l (i) ; Les paramètres Vl, V2, V3, V4 servant à.la détermination de l'adresse sont précalculés de la façon suivante: V3 = 1 V4 = 0 Z C-Q V = ----- &Q CaZ V2 = - aQ L'ensemble des coefficients intervenant dans ces calculs est connu. En effet, les coordonnées sont échantillonnées comme suit: x(k) X +k-AX o y(j) = Y +j by zCi) = Z +k Z (n) p(m) = P + m aP q(L) = + L LQ R1 + R2 C = -----______ R2 + t (koJo,n) M 0 0 0 l dépend de k0, j et n n O' o l (k j,nO)=X(k)sinC(n)-y(j)Cos(n) M O 0 O O Dans le cas d'une reprojection e- geométrie conique 3D, l'adresse recherchée est l (i). La partie de l'espace à reconstruire qui est figée est pour k=k et j=j fixant ainsi * à *(k) et 0 0 p à p(j). O La première dimension que l'on fixe dans l'espace connu définit un plan m=m qui fixe l'abscisse x à x(m). Les paramètres Vl, V2, V3, V4 servant à la détermination de l'adresse sont précalculés de la façon suivante: Q C-Z O O V1 = ---- a z V2 = - 4z R +l (h,ff,n') R2 M o'mO avec C = --------- R1 + R2 P(j) [R2+x(m o)sin*(k 0?-y0oseo(k.) n' = ----------------_________________ i Y ER1+R2)sin$(k)+p(j)cos4(k) (R1+R2) LO(m o)cos(k) + Y osin4(k o) Y (Rl+R2) sin *(k) + P(jo)cOsC(k 0 et avec l (k n')x(m M O' m n)X(mo.sin*(h)-y(n')cos(hO) V3 = 1 V4 = 0 R1: distance source - centre du repère objet. R2: distance centre du repère objet-détecteur. L'ensemble des autres coefficients intervenant dans ce précalcul est connu. L'espace a été échantillonné de la façon suivante: x(m)=X +m nX y(n)=Y +n lY zCl)=Z +l AZ 0 (k) p(j)=P +j & P q(i)=Q+i aQ Conformément au procédé de l'invention, l'étape de pondération consiste à appliquer un facteur de pondération f tel que: V5 f = ---2-- (V6+eV7) les paramètres V5, V6, V7 étant fonction de la géométrie d'acquisition et du type de projection effectuée. Dans le cas de rétroprojection en géométrie parallèle 2D, conique, parallèle et conique 3D et de reprojection parallèle 2D et 3D, V5 est choisi par l'utilisateur V6 = I V7 = 0 Dans le cas de rétroprojection en éventail, V5 est choisi égal à (R1+R2), V6=V3, V7=V4. Dans le cas de reprojection éventail 2D et conique: V5 = 1 V6 = 1 V7 = 0 Sur la figure 11, l'on a représenté le schéma d'un dispositif permettant de mettre en oeuvre l'invention. Ce dispositif a une architecture pipe- line. Ce dispositif comporte des moyens de génération d'adresse 10, des moyens de détermination des contributions par interpolation 20. Le dispositif comporte également des moyens de pondération 30, des moyens de sommation 40. Le dispositif comporte en outre des moyens de stockage des données de l'espace connu 50; et des moyens de stockage des données de l'espace inconnu 60 et 70. Les moyens de génération d'adresse reçoivent les paramètres V1, V2, V3, V4 et sont aptes à effectuer le calcul de l'adresse ADR à déterminer. Le calcul est réalisé par des opérateurs incrémentaux constituant le générateur qui effectuent à chaque cycle une opération du type incrémental. Les moyens de mémorisation 50 sont adressés par les moyens de génération d'adresse 10 et délivrent les (o +1) valeurs consécutives V(ADR- È) à V(ADR+ Y) avec f+ï =, P et ' entiers positifs. Les moyens d'interpolation reçoivent les %+L) valeurs V(ADR-e) à V(ADR+ i) et déterminent à chaque cycle la contribution non pondérée. Le coefficient COEF utilisé pour effectuer l'interpolation est par exemple donné par le générateur d'adresse, ce coefficient étant généralement la partie fractionnaire de l'adresse. Le nombre cd est égal à 1 pour une interpolation linéaire. Les moyens de pondération 30 réalisent le produit de la contribution non pondérée avec le facteur f. Pour ce calcul, les moyens de pondération reçoivent les paramètres V5, V6, V7 et la donnée issue de l'interpolateur 20. En particulier, pour l'opération de rétroprojection éventail 2D, les paramètres V4 et V7 peuvent être issus du générateur d'adresse 10. L'invention concerne également une application du procédé à la reconstruction par synthèse et inversion de la transformée de Radon qui met en jeu qatre espaces. Les moyens de sommation 40 permettent d'ajouter la contribution ainsi pondérée au voxel (ou pixel) traité. La zone de donnée de sortie peut être, comme représenté sur la figure 11, constituée de deux bancs mémoire 60 et 70 travaillant en ping-pong. La figure 12 montre l'espace volume à reconstruire de coordonnées x, y, z et l'espace projection connu de coordonnées p, q,. La figure 13 montre l'espace de - Radon de coordonnées p, e, f et le repère e, i, Y, k lié à l'espace du volume à reconstruire. La figure 14 montre l'espace des projections réarrangées de coordonnées r, z, c. Le but de cette méthodde est de passer de L'espace projection connu à l'espace à reconstruire et donc inconnu. La reconstruction s'effectue en 3 étapes: La première étape consiste en une reprojection bidimensionnelle en géométrie parallèle de l'ensemble des radiographies prises une à une comme des Images formant ainsi un espace bidimensionnel image connu; L'espace inconnu associé à cette image est la projection de coordonnées d, *=*(m) comme le montre La figure 15. Une série de reprojections bidimensionnelles sur La totalité des radiographies permet d'obtenir un volume de données intermédiaire de coordonnées d, &0. Une opération de rééchantiLLonnage sur ce volume d,/, * permet d'obtenir aisément l'espace de Radon de coordonnées P, O, y. Après filtrage des données de l'espace de Radon, la deuxième (figure 16) étape est mise en oeuvre et on considère l'espace connu comme étant l'espace de Radon de coordonnées P, e, 9 et l'espace inconnu l'epace des projections réarrangées de coordonnées r, z, I. Une rétroprojection bidimensionnelle en géométrie parallèle est effectuée en prenant comme espace de projection connu chaque projection O, e à Y= f(m)= constante et comme espace bidimensionnel à reconstruire l'image de coordonnées r, z à À='?(m)= constante. En répétant ainsi cette opération pour l'ensemble des valeurs de d(m), l'espace ddes projections réarrangées de coordonnées r, z, óPest obtenu. Après normalisation des données de projections réarrangées r, z,, la (figure 17) troisième étape est mise en oeuvre et on considère l'espace connu comme étant l'espace des projections réarrangées de coordonnées r, z, f et l'espace inconnu l'espace volume à reconstruire de coordonnées x, y, z. Une rétroprojection bidimensionnelle en géométrie parallèle est effectuée en prenant comme espace de projection connu chaque projection r, à zz(K)= constante et comme espace bidimensionnelle à reconstruire l'image de coordonnées x, y à z=z(k)= constante. En répétant ainsi cette opération pour l'ensemble des valeurs de z(k), l'espace volume à reconstruire de coordonnées x, y, z est obtenu. Le tableau A suivant permet de synthétiser l'application décrite. TABLEAU A Espace connu Espace recherché projection ler ETAPE Radon intermédiaire radiographieques --------------- d, t, * (m) p. q, *(m) reprojection 2D géométrie parallèle Radon Intermédiaire -------------- Radon P, e, d, ô, Rééchantillonnage 2ème ETAPE Radon P,e,(mO) ------------- Projections réarangées Rétroprojection r. z, (m) 2D géométrie parallèle 3ème ETAPE Projections réar- -------------- Volume final recherché z,, zCk) Rétroprojection x, y, z(k) 0 0 2D géométrie parallèle REVENDICATIONS 1. Procédé de détermination d'un espace à partir d'un espace discret connu par projection ou rétroprojection, pour la reconstruction d'une image bidimensionnelle ou tridimensionnelle d'un corps à analyser, à partir d'un ensemble de projections obtenu, au moyen d'un dispositif d'acquisition de géométrie définie comportant une source de rayonnement et un ensemble de détecteurs, l'un des espaces étant défini par un ensemble de points élémentaires repérés respectivement suivant trois axes par des indices (l, m, n) correspondant aux coordonnées suivant ces axes; l'autre espace étant défini par un ensemble de points élémentaires que l'on repère également respectivement par des indices (i,j,k) suivant trois axes de cet espace, ces indices (i, j, k) correspondant aux coordonnées suivant ces axes, la détermination de l'espace inconnu est obtenue par le calcul de la valeur numérique de chaque point de l'espace inconnu à partir de la détermination des contributiors apportées par les valeurs numériques de points correspondants de l'espace connu, le procédé étant caractérisé en ce qu'il comporte une phase élémentaire consistant à déterminer de façon itérative pour chaque point d'un ensemble monodimensionnel de points de l'espace inconnu portés par un axe E donné de cet espace, la contribution apportée par le point correspondant et éventuellement par les points voisins d'un ensemble monodimensionnel de points de l'espace connu portés par un axe F de cet espace, cette phase comportant de façon itérative pour chaque point les étapes suivantes: - déterminer l'indice associé à la coordonnée sur l'axe F du point correspondant de l'espace connu, en prenant la partie entière ADR de l'expression généralisée a (e), avec Vl + eV2 a (e) = --------- 0 V3 + eV4 e étant l'indice d'itération correspondant à la coordonnée sur l'axe E du point considéré de l'espace inconnu et V1, V2, V3, V4 étant des paramètres constants dans La phase élémentaire, l'ensemble de ces paramètres étant déterminé préalablement à cette phase en fonction de la géométrie d'acquisition choisie; du type d'opération (reprojection ou rétroprojection) effectuée et des ensembles monodimensionnels considérés, - déterminer la contribution C à apporter à la valeur numérique du point considéré de l'espace inconnu par interpolation sur les valeurs numériques du point correspondant d'indice ADR et éventuellement des points voisins, - pondérer cette contribution par un facteur f déterminé par la relation suivante: V5 f = ---------- (V6 + eV7) VS, V6, V7 étant des paramètres constants dans la phase élémentaire, l'ensemble de ces paramètres étant déterminé préalablement à cette phase en fonction de la géométrie d'acquisition choisie, du type d'opération et des ensembles monodimensionnels considérés, - additionner cette contribution pondérée à la valeur numérique du point considéré de l'espace inconnu, détermirée lors de phases élémentaires précédentes ayant mis en correspondance le point de l'espace inconnu à des points d'autres espaces monodimensionnels connus, l'ensemble des contributions obtenues à la fin du procédé de détermination de L'espace inconnu constituant la valeur finale de ce point. 2. Procédé de détermination d'un espace à partir d'un espace discret connu, selon la revendication 1 dans lequel la phase élémentaire représente une opération de rétroprojection en géométrie parallèle, l'espace connu étant un ensemble de données de projection indicé par l,m,n; m et n étant fixes, l'espace inconnu étant l'espace image indicé par i, j, k, (j et k) ou (i et k) étant fixes; l'indice recherché a (e) est soit l'indice l (4) lorsque i et k sont fixés, l'axe E étant indicé pr et l'axe F par l, soit l'indice l (j) lorsque j et k u sont fixés, l'axe E étant indicé par j et l'axe F par L; et l'on choisit alors: V3 = 1 V4 = O et pour i fixé à i, k=0, m=0, n=n, L (j) étant O lo(J recherch() V1 = L + i LI O O V2 = J e=j ou pour j fixe à j, k=O, m=O, n=n, l (i) étant 0 0 recherché v1 = L t j AJ o o V2 = I e=i e = i les termes L, p I, AJ étant prédéterminés, l étant un nombre entier réel de partie entière l. 3. Procédé de détermination d'un espace à partir d'un espace discret connu, selon la revendication 1 dans lequel la phase élémentaire représente une opération de reprojection en géométrie parallèle, l'espace connu étant l'espace image indicé par l, m, n; m et n ou l et m étant fixés, l'espace inconnu étant l'ensemble de données de projection indicée i, j, k ou j et k sont fixés; l'indice a (e) O recherché est l'indice l (i) ou n (i), l'axe E est 0 0 alors indicé par i ou l'axe F est indicé par I ou n et l'on choisit alors: m fixé à O, n à no, j à jo, k à O, on cherche lo(i) V3 = 1 V4 = 0 V1 = h (n, j) V2 = bL (j) 1 0 e = i ou m fixé à 0, h à l, j à jo k à O; on cherche n (i) O' 0 V3 = 1 V4 = O Vl ='h (l io V2 = L (J e = i 2 0 Les termes h1, h, bL et àL étant I2 1 2 prédéterminés et n et l0 sont des nombres réels de partie entière respectivement n et l. 4. Procédé de détermination d'un espace à partir d'un espace discret connu, selon la revendication 1 dans lequel la phase élémentaire représente une opération de rétroprojection en géométrie en éventail à pas constant, l'espace connu est l'ensemble de données de projection indicées par l,m,n; m et n étant fixés, l'espace inconnu étant l'espace image indicé par i, j, k, (i et k) ou (j et k) étant fixés, l'indice recherché a (e) est soit-l'indice o l (i) soit l (j) l'axe E étant indicé par i ou par j et 0 o l'axe F par l et l'on choisit alors: k=O, m=O, n=n =Jo, l (i) étant recherché, o Vl = h (jo n) 335V2 = h (n0 V2 = h (n) 4 0 V3 = h (j n ) V4 = h (n) 6 0 e = i ou:k=O, m=O, n -n, i0 (j) étant recherché, V1 = hl (i n) 3 0' 0 V2 = h' (n) no) 4 O V3 = h' (i, n V4 = h' (n) 6 0 e = j. les termes h3 h4 h, h6, h', h' étant 3 45 6 3' 4 5' 6 prédéterminés et l est un nombre réel de partie o entière l. 5. Procédé de détermination d'un espace à partir d'un espace discret connu, selon la revendication 1 dans lequel la phase élémentaire représente une opération de reprojection en géométrie en éventail, l'espace connu étant l'espace image indicé par L, m, n; m et n ou l et m étant fixés, l'espace inconnu étant l'ensemble de données de projection indicées par i, j, k ou j et k sont fixes, l'indice recherché a (e) est l'indice L (i) ou n (i) l'axe E 0 0 - O étant indice par i et l'axe F par l ou n et L'on choisit alors: m fixé à 0, n à n, j à j k à 0, l (i) étant 0 'O recherché, Vl = h (n, j) 7 0 O V2 = h (n, j 8 0 0o V3 = h (j) 9 0 V4 = h (j o e = i ou.,m fixé à 0, l à L, j à j k à O, n (i) étant O'i. J o' ' 0 recherché, Vl = h' (l, j) V2 = h' (l, j 8 0 0 V3 = h' (j) 9 0 V3 -- h' 1 j 0 O e =i Les termes h7, h8, h h10, h'7, h' h' 9 7'8' h9, h' étant prédéterminés et n, L étant des nombres 0O réels de partie entière respectivement n et l. 6. Procédé de détermination d'un espace à partir d'un espace discret connu, selon la revendication 1 dans lequel la phase élémentaire représente une opération de rétroprojection en géométrie conique; l'espace inconnu étant un espace image indicé par i, j, k (j et k étant fixés), l'espace connu étant l'espace de projection indicé par l, m, n (n étant fixé et m fixé puisnue fonction de j et k), l'indice recherché a (e) est l'indice L i) l'axe E est 0 O indicé par i et l'axe F par l, et l'on choisit alors: v3 = 1 V4 = 0 V = Z C-Qo V1 = Zoà C Q V2 = Q les termes Z, C, Q, AZ et eQ étant prédéterminés et 0 0- O étant un nombre réel dont la partie entière est l. 7. Procédé de détermination d'un espace à partir d'un espace discret connu, selon. la revendication 1 dans lequel la phase élémentaire représente une reprojection en géométrie conique, l'espace inconnu étant l'espace de projection indicé par i, j, k, (j et k étant fixés), et l'espace connu étant l'espace image indicé par l, m, n (m étant fixé, n étant fixé puisque fonction de j et k); l'indice recherché a (e) est l'indice l (i) l'axe E est indicé par i et l'axe F par l et l'on choisit alors: V3 = 1 V4 = 0 V = QC-Z o z V2 CAO V2 = --- les termes Z0, OO, C, bZ et SQ étant prédéterminés et l étant un nombre réel dont la parie entière est l. 8. Procédé de détermination d'un espace à partir d'un espace discret connu, selon les revendications 1 et 2, ou 1 et 3 ou 1 et 6, caractérisé en ce qu'une géométrie parallèle et conique dans le cas ou la phase élémentaire représente une rétroprojection et pour une géométrie parallèle lorsque cette phase représente une reprojection, le facteur de pondération f est déterminé en prenant: V7 = 0 V6 = 1 et V5 est prédéterminé. 9. Procédé de détermination d'un espace à. partir d'un espace discret connu selon les revendications 1 et 4, caractérisé en ce que dans l'étape de pondération pour une géométrie éventail dans le cas d'une rétroprojection, le facteur de pondération f est déterminé en prenant: VS prédéterminé à partir des paramètres du dispositif d'acquisition V6 = V3 V7 = V4 10. Procédé de détermination d'un espace à partir d'un espace discret connu, selon les revendications 1 et 5 ou 1 et 7, caractérisé en ce que dans l'étape de pondération pour une géométrie en éventail et conique dans le cas de la reprojection, le facteur de pondération f est déterminé en prenant: V5 = 1 V6 = 1 V7 = 0 11. Dispositif de mise en oeuvre du procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comporte outre le dispositif d'acquisition: - des moyens de génération d'adresse recevant les paramètres V1-V4, délivrant l'adresse ADR qui est la partie entière de a (e); - des moyens de mémorisation des valeurs numériques des points de L'espace connu, adressées par les moyens de génération, délivrant l'ensemble de (K+1) valeurs numériques des points de l'espace connu de V(ADDR-p) à V(ADR+ ) ( et' entiers positifs; p+' = () correspondant à l'indice ADR, o étant l'ordre d'interpolation; - des moyens d'interpolation recevant les valeurs numériques V(ADR-e) à V(ADR+) continues en mémoire et adressablespar l'indice ADR et déterminant à la fin de l'interpolation une contribution non pondérée à la valeur numérique de chaque point correspondant de l'espace inconnu; - des moyens de pondération recevant chaque contribution non pondérée et les paramètres V5, V6, V7 pour délivrer chaque contribution pondérée; - des moyens de sommation pour additionner à la valeur numérique du point en cours de l'espace inconnu, la contribution pondérée relative à ce point en cours; - des moyens de mémorisation contenant la valeur en cours. 12. Application du procédé selon l'une quelconque des revendications 1 à 10, à la reconstruction par synthèse et inversion de la transformée de Radon. 13. Application du procédé selon la revendication 12, caractérisé en ce que la reconstruction comporte trois étapes: - une première étape consistant en une reprojection bidimensionnelle en géométrie parallèle permettant d'obtenir une transformée de Radon intermédiaire; - une deuxième étape consistant en une rétroprojection bidimensionnelle en géométrie parallèle permettant d'obtenir des projections réarrangées; une troisième étape consistant en une rétroprojection bidimensionnelle à géométrie parallèle permettant d'obtenir le volume final recherché.