Brevet FR2562371A1 — publié le 19851004

PROCEDE DE RECONSTRUCTION D'UNE IMAGE A HAUTE RESOLUTION, PAR TOMODENSITOMETRIE L'invention concerne un procédé de reconstruction d'une image à haute résolution, par tomodensitométrie, ayant pour but d'amé- liorer la qualité de l'image reconstruite sans alourdir la structure technologique des moyens de mesure de l'absorption de rayonnement pénétrant et sans modifier sensiblement le processus de prise des données. Les tomodensitomètres dits de "troisième génération" se composent essentiellement d'une source de rayonnement pénétrant (des rayons X) émettant un faisceau plat en éventail, d'une rangée courbe de détecteurs placés en regard de cette source pour mesurer la fraction non absorbée du rayonnement dans le plan de la coupe à imager, de moyens pour faire tourner la source et les détecteurs dans le plan de cette coupe et d'un ordinateur pour traiter les données fournies par les détecteurs et reconstruire une image à partir de ces données. La résolution spatiale de l'image reconstruite dépend notamment du nombre des détecteurs embrassant l'éventail du faisceau de rayons X. Pour une ouverture de faisceau donnée, plus le nombre de détecteurs est élevé et plus la résolution spatiale est bonne. Cependant, l'ensemble multi-détecteurs est l'un des organes les plus coûteux du tomodensitomètre et l'un de ceux qui est le plus difficile à fabriquer avec une bonne fiabilité. Par conséquent l'augmentation du nombre de détecteurs sur un secteur donné correspondant à l'éventail, se heurte à des limites technologiques. L'un des buts de l'invention est d'améliorer sensiblement la qualité de l'image sans modifier la structure technologique de l'ensemble source-détecteurs et sans modifier sensiblement la pro- cédure de prise des données proprement dite. Un autre but de l'invention est aussi de réduire le coût de construction d'un tomodensitomètre ayant une résolution spatiale comparable à celle des tomodensitomètres existants, cette réduc- tion de coût étant obtenue en utilisant un ensemble multidétecteurs comportant deux fois moins de détecteurs que les machines exis- tantes. Dans cet esprit, l'invention concerne essentiellement un pro- cédé de reconstruction d'images par tomodensitomètrie, du type consistant à faire tourner un ensemble source-détecteurs dans le plan d'une coupe à imager, ladite source émettant un faisceau de rayonnement pénétrant en forme d'éventail et à relever et mémo- riser une succession de vues pour différents angles/5 de rotation de l'ensemble source-détecteurs par rapport à un système dcaxes de référence fixe dudit plan, les détecteurs étant agencés pour con- verger vers le foyer de ladite source et chaque vue étant échantil- lonnée pour des angles e représentatifs des détecteurs et repérés a partir de la droite joignant le centre de rotation de l'ensemble sourcedétecteurs et ledit foyer, le décalage angulaire entre deux détecteurs adjacents étant AY constant et le décalage angulaire entre vues étant constant, caractérisé en ce qu'il consiste, le décalage angulaire entre ladite droite et le détecteur le plus proche de celle-ci étant réglé à j4A, - à relever et mémoriser une série de vues réparties sur un tour complet dudit ensemble source-détecteurs; - à élaborer pour chaque vue des valeurs intermédiaires représentant un échantillonnage entrelacé avec celui desdits détecteurs, chaque valeur intermédiaire, d'angle Y donné, étant obtenu par une combi- naison d'au moins deux valeurs d'absorption linéaires déduites de valeurs de mesure prélevées dans des vues situées de part et d'autre d'une vue fictive dont Pl'image par la transformation T(O4, y) = ( 67+ (2c - 1)-ff, -/) contient le point représentatif de ladite valeur intermédiaire dans l'espace [riYJavec 6; =/3+ y et k entier relatif, lesdites valeurs de mesure desdites vues étant prélevées sur le détecteur d'angle opposé à celui de ladite valeur intermédiaire recherchée;et - à appliquer un algorithme de reconstruction, connu en soi, com- prenant une convolution et une rétroprojection aux vues ainsi complétées par lesdites valeurs intermédiaires. L'invention sera mieux comprise à la lumière de la description qui va suivre d'un procédé conforme à son principe et des éléments constitutifs essentiels d'un tomodensitomètre fonctionnant selon ce procédé, faite en référence aux dessins annexés dans lesquels: - la figure 1 est une représentation schématique des principaux éléments constitutifs d'un tomodensitomètre mettant en oeuvre l'invention; - la figure 2 est un représentation dans l'espace [x 6Ides vues prises par un tomodensitomètre classique fonctionnant en mode pulsé pour la reconstruction d'une image; la figure 3 est une représentation dans le même espace[tQ3]des vues prises par un tomodensitomètre du même type dans lequel on a effectué l'une des modifications nécessaires à la rmise en oeuvre de l'invention, à savoir un décalage statique entre la source et les détecteurs; - la figure 4 illustre, dans l'espace r"]june transformation appli- cable aux vues de la figure 3 pour une mise en oeuvre du procédé de l'invention; - la figure 5 illustre dans l'espace L9Q] le résultat de cette transformation et montre le processus d'interpolation permettant de compléter les vues; - la figure 6 illustre la mise en oeuvre d'un algorithme conforme à l'invention, pour le calcul d'une valeur intermédiaire de la première vue V1 - la figure 7 illustre une phase de calcul du même algorithme;et - la figure 8 est une représentation dans l'espace rld'une vue, dans le cas d'un tomodensitormètre fonctionnant en mode continu. La suite de la description en référence aux figures I à 7, se place dans le cas particulier d'un tomodensitomètre fonctionnant en mode pulsé. Ce genre de tomodensitomètre est bien connu. Dans la plupart des cas, le mode pulsé se concrétise par le fait que le faisceau en -5ventail n'est 5mis que pendant un intervalle de temps très bref au moment o l'ensemble source-détecteurs se trouve à un angle/d prédéterminé. Dans ce cas, il suffit que la lecture de tous les détecteurs soit achevée lorsque l'ensemble source-détecteurs parvient à l'angle/, suivant. Cependant, une autre façon de réaliser ledit mode pulsé peut consister à générer en permanence le faisceau en éventail mais à "lire" tous les détecteurs simultanément à chaque angle/. En se reportant à la figure 1, on a représenté les éléments constitutifs essentiels d'un tomodensitomètre. Cet appareil com- porte un ensemble source-détecteurs mobile en rotation dans un plan de coupe à imager. Le centre de rotation est O et le plan de la coupe sera considéré conmme étant celui de la figure 1. On associe à ce plan de coupe un système d'axes de référence fixe XOY, orthonormés. La source est symbolisée sur le dessin par son foyer F. De façon bien connue, il s'agit d'une source de rayonnement péné- trant et plus particulièrement des rayons X, émettant un faisceau 11 en forme d'éventail. Dans la suite du texte, on parlera simplement de la source F. Le bottier 13 renfermant tous les détecteurs di a une forme courbe pour que ces détecteurs soient agencés suivant un arc de cerdce CC' géométriquement centré sur le foyer F. On sait que le boîtier peut aussi être rectiligne mais renfermer des détecteurs agencés pour converger vers le foyer. Dans l'exemple décrit,le faisceau en éventail 11 embrasse Parc de cercle CC' occupé par les faces d'entrée des détecteurs di. Ainsi, Pimage reconstituée sera à l'intérieur du cercle 16 de centre O, inscrit dans les limites de l'angle CFC'. On peut utiliser un ensemble comportant 1024 détec- teurs d.. Le bottier 13 et la source F sont solidaires d'un même support tournant (centre de rotation O) et fixes Pl'un par rapport à l'autre de sorte qu'on peut repérer la rotation de l'ensemble source- détecteurs par l'angle/3 entre l'axe OY et la droite OF. Les valeurs d'atténuation du faisceau de rayons X pour des angles/3 prédé- terminés -sont déduites des valeurs de mesure prélevées aux sorties des détecteurs et mémorisées dans une mémoire M d'un ordinateur. Pour chaque angle// prédéterminé, un ordre de lecture des détec- teurs diest émis, par exemple au moyen d'un codeur incrémental 15 agencé pour lire des informations d'angle solidaires de l'ensemble tournant source-détecteurs et pour délivrer (liaison de pilotage 17) des ordres de commande de lecture de l'état de tous les détecteurs. L'ensemble des valeurs d'atténuation linéaire déduites des lectures des détecteurs pour un angle fi donné est appelée "vue". On peut donc considérer que chaque vue est échantillonnée de par la structure même de l'ensemble de détecteurs, pour des angles > représentatifs des détecteurs. Chaque angle ' indiquant la position d'un détecteur di est repéré à partir de la droite OF. De ce fait, on appelera Ly le décalage angulaire constant entre les centres des faces d'entrée de deux détecteurs adjacents. Si on considère pour chaque angle tg d'une vue, l'axe OX perpendiculaire en P au rayon d'angle Y frappant un détecteur quelconque, on démontre facilement que l'angle 1 entre l'axe OX du repère fixe et l'axe OX est tel que: ceci dans tous les cas de figures, les angles 0,P et Y étant orientés. On a encore: ' = - D'après cette dernière relation, on voit que les différentes vues (pour des angles P différents) ont une représentation extrême- ment simple dans l'espace Éa, 3J. Comme le montre la figure 2, ces vues sont portées par des * peir isti droites parallèles à4a première bissectrice et décalées de âf le long de l'axe des 9. On peut porter sur ces droites des points représentant les coordonnées, dans l'espace, des valeurs échan- tillonnées effectivement acquises à la lecture des différentes vues. Ainsi, si on suppose que le nombre de détecteurs est pair, soit 2n et que la droite OF passe exactement entre les deux détecteurs du centre d et d ceux-ci apparaissent décalés de | 2- de part et d'autre de l'axe. Le but de l'invention est d'obtenir pour chaque vue des valeurs échantillonnées intermédiaires dont les coordonnées dans l'espace [r l]sont figurées par des étoiles X- sur la figure 2. Autrement dit, le but de l'invention est d'obtenir des vues com- portant 4n valeurs échantillonnées avec un ensemble de 2n détec- teurs seulement et d'améliorer ainsi la résolution spatiale de l'image reconstruite. Selon une particularité de l'invention, on provoque avant tout un décalage angulaire permanent et prédéterminé entre l'ensemble des détecteurs (c'est-à-dire le boîtier 13) et la droite OF. Selon l'exemple choisi, le décalage angulaire entre la droite OF et le détecteur le plus proche (do) est choisi égal à + 4. Cela signifie que l'ensemble des détecteurs est décalé "à gauche" en considérant la figure 1, de façon que le centre de la face crd'entrée du détecteur do soit décalé angulairement de +.- par rapport à OF et que le détecteur d1 soit décalé de --2 par rapport à cette même droite. La figure 3 montre la configuration des vues dans l'espace ] après ce décalage. Pour la mise en oeuvre de l'invention, on relève une série Vi de N vues réparties sur un tour complet de l'ensemble source- détecteurs, avec N = 27 Les valeurs intermédiaires recherchées si sont également repérées sur les différentes vues de la figure 3, chaque emplacement de valeur intermédiaire dans l'espace [,Xe étant équidistant des emplacements de deux détecteurs réels voisins. Soit la transformation: T (O., ï) = (9-tr, -Y) On peut démontrer que la fonction d'absorption d'une vue est invariante par cette transformation T. Or cette transformation T peut se décompose dans l'espace [6 en une translation - Vr suivie d'une symétrie par rapport à. Cette transformation est illustrée en soi à la figure 4. Si on applique ladite transformation à I. a la vue Vi, le résultat est la vue V'.. Si on applique ainsi la transformation T à la série des vues Vi, on constate que les valeurs échantillonnées des "vues" V'. trans- I formées ont la même ordonnée à que des valeurs intermédiaires si recherchées. On peut donc compléter les vues Vi par des inter- polations à partir de valeurs d'absorption linéaires correspondant à des valeurs échantillonnées des "vues" V'i, à > constant, dans 'espace[,, comme indiqué par les flèches sur la figure 5. Selon l'exemple, chaque valeur intermédiaire d'angle {donné est obtenue par une combinaison de deux valeurs d'absorption linéaire déduites des valeurs échantillonnées de même ordonnée >, prélevées dans des "vues" V'. situées de part et d'autre de la valeur intermédiaire I recherchée. La combinaison peut simplement consister en une interpolation linéaire entre les valeurs d'absorption linéaire corres- pondantes de deux "vues" adjacentes Vyi, comme indiqué figure 5. Autrement dit, chaque valeur intermédiaire est obtenue en faisant la moyenne entre ces deux valeurs d'absorption linéaire-, précitées, en les pondérant respectivement par l'inverse de l'écart (en prc- jection sur l'axe des 9, figure 5) entre leur position respective et la position de la valeur intermédiaire recherchée. Pour certaines "vues" V'i nécessaires aux calculs d'interpolation suggérés par la figure 5, le transformation T fait appel à une vue Vi située en dehors de l'intervalle L0, 2 W] des vues réellement acquises. Cependant, ces dernières peuvent être théoriquement prolongées dans l'espace [0, Y1 de -omà +i en raison de la périodicité par rapport à &9, de la fonction d'atténuation linéaire au-delà d'un tour complet (en raison de la condition N), ce qui peut se traduire ainsi: V =Vj>ii =Vji=j ModuloN Autrement dit, la transformation T peut être complétée de la façon suivante: avec k entier relatif. Il existe toujours une valeur de k telle que la vue V. sur laquelle est appliquée la translation + (2 k - 1) soit bien l'une des vues acquises Vi. La Iigure 5 illustre la justification du processus de l'invention. Ce mede de calcul des valeurs intermédiaires doit raaintenant être explicité car les Ivues' V'i ne représentent qu'un réarrangement des vues Vi et le calcul de chaque valeur intermé- diaire peut se faire en fait à partir de chaque "emplacement" de valeur intermédiaire recherchée, dans l'espace [,2. C'est ce qu'il- lustre la figure 6 o on n'a représenté que la première vue V1 ainsi que deux "vues" V':., résultat de la transformation T appliquée aux vues VN et VIN 1 (N étant pair), lesquelles sont aussi repré- 2 2 +, sentées. Les vues N et VN sont situées de part et d'autre d'une vue fictive VFo dont la transformée VF'o par la transformation T, passe par les coordonnées de la valeur intermédiaire so de la vue V1 dans l'espace[G,0. Selon l'invention, il suffit donc pour chaque valeur intermédiaire de chaque vue de rechercher la "vue fictive" correspondante dont la transformée par la transformation T passe par les coordonnées de ladite valeur intermédiaire, de sélectionner 1 les deux vues réelles adjacentes Vi, Vi + I situées de part et d'autre de ladite vue fictive et de prélever dans ces deux vues les valeurs d'absorption d'angle -' ( létant l'ordonnée de so dans l'espace [z,J) qui correspondent à des valeurs de mesure-d'un même détec- teur réel et d'effectuer une interpolation linéaire entre ces valeurs d'absorption linéaire. Il est à noter que, pour simplifier, on a supposé ci-dessus que le nombre N de vues acquises Vi était pair; le processus de l'invention s'applique néanmoins tel quel si N est impair. On va maintenant décrire Pun des algorithmes de calcul possibles, constituant une mise en oeuvre actuellement préférée du procédé décrit ci-dessus. Dans l'exemple décrit, on se place dans le cas d'un décalage des détecteurs "à gauche" tel qu'illustré sur les dessins et dans le cas d'un mode de fonctionnement pulsé o, comme on sait, les valeurs échantillonnées des vues Vi sont portées par des droites parallèles à la première bissectrice de l'espace [j4,. On va d'abord décrire en référence à la figure 7, l'algorithme permettant de compléter la vue V1. L'équation de la vue V1 est: n h r a v On cherche à calculer la valeur intermédiaire en si. Soit x. l'écart entre l'intersection de la vue à compléter avec l'axe (en l'occurence l'origine, dans le cas de la vue V1) et l'interconnection de la transformée de la vue fictive (passant par s.) avec le même axe. On a: xj = 2' sj (A) Le point d'intersection de la vue fictive correspondante avec l'axe des est donc: yj = "+xj (B) Le numéro de la vue réelle KG situé "à gauche" de cette vue y. fictive est donc déterminé par la valeur entière du rapport,. De façon plus précise: Zai G(i) Et Yi) + 1, modulo N (C) et le numéro de la vue KD située "à droite" de la vue fictive est donc: KD(j) = KG(j) + 1, modulo N (D) Il reste maintenant à déterminer "l'ordonnée". du détecteur à prendre en compte dans chaque vue KG et KD, soit, d'après la transformation T: Ydj = ysj. Or si on se reporte à la figure 2: s6. =A et la transformée KG' de KG. Par conséquent, on peut déduire de r. J la pondération à apporter aux valeurs d'atténuation linéaire corres- pondant aux valeurs de mesure prélevées sur le détecteur d des vues KG et KD. Soient hKD (dcj) et hKG (di), les valeurs d'atténuation linéaire correspondant au détecteur i dans les vues réelles KD(j) et KG(j), et soit: tor R(j) = Do; R0)- /3 on vérifie que la valeur d'atténuation linéaire h (sj) à attribuer à l'emplacement s ide la vue V1 (autrement dit la valeur intermédiaire s. elle-même) s'exprime: 1S h(sj) = R(j) [hKD (d) - hKG (d1)J + hKG (cI) (F) ! La programmation des expressions A à F ci-dessus est à la portée de l'homme du métier tout comme la conception d'un processeur de calcul spécialisé couplé à la mémoire tampon M. Pour compléter les vues suivantes Vm (c'est-à-dire V2, V3, 0... VN) on a, au lieu de la relation (A): xj (m - 1)= 2 V soit: x; = 2ysj + (m - I)A (A') L'expression A' remplace A dans l'algorithme de réar- rangement et d'interpolation mais les autres relations B à F restent inchangées. Il est aussi intéressant de remarquer que d'une vue à l'autre R(j) reste inchangé et que KG(j) est incrémenté de 1. Par conséquent, il est possible dans le logiciel de calcul, de précalculer r et KG uniquement pour la première vue et de procéder ensuite par incrémentations de KG. Par ailleurs, on vérifie que le même algorithme est valable pour un décalage des détecteurs "à droite", de sorte que le procédé décrit s'applique pour un décalage de l'ensemble de détecteurs de:. L'algorithme qui précède a été développé avec l'hypothèse d'un fonctionnement en mode pulsé. Il est facile de l'adapter à un fonctionnement en mode continuo Selon un mode de fonctionnement continu usuel, on suppose que la source S émet en permanence et que la lecture des valeurs de mesure aux sorties des détecteurs, pour chaque vue, se fait par scrutation progressive d'une extrémité à l'autre de l'ensemble des détecteurs, par exemple depuis le dé- tecteur d.n jusqu'au détecteur dn - l Deux paramètres supplémentaires doivent donc être pris en considération: - la vitesse de rotation, v, de l'ensemble source-détecteurs - l'intervalle de temps e séparant l'acquisition des valeurs échantil- lonnées correspondant à deux détecteurs voisins, en considérant une vitesse de scrutation constante. La figure 8 illustre la façon dont la représentation d'une vue dans l'espace [ é,J] évolue si on prend en compte les paramètres v et; La vue portée par la droite en trait interrompu correspond au mode pulsé tandis que la même vue portée par la droite en trait plein correspond au mode continu. -*Les valeurs échantillonnées se déplacent parallèlement à l'axe 9 puisque P =/+ y-; /,3 variant pendant l'acquisition. Les deux droites concourrent en d n puisque la scrutation des détecteurs commence par ce détecteur et la valeur échantillonnée suivante, d n + se trouve déplacée en d' n d'une distance v.V dans l'espace,f,], parallèlement à l'axe des. Soit a l'angle entre tga +v. = ô. soit = ga Les vues sont donc portées par des droites parallèles à la droite d'équation = (i) Par conséquent, pour la vue Vt la relation (A) devient: xj = 2 (1 +K) /sj et pour une vue Vm quelconque, la relation (A') devient: xj 2 (Il & nouvelles valeurs de xj définies ci-dessus. On a vérifié que le procédé décrit s'applique aussi sans changement quelque soit le type de scrutation en mode continu. Dans ce cas, les valeurs échantillonnées des vues ne sont plus obligatoirement portées par des droites mais par des courbes quel- conques auxquelles on peut appliquer la même transformation T pour en déduire les valeurs intermédiaires recherchées à partir d'inter- polations analogues à celles décrites plus haut. En se reportant à nouveau à la figure 1, on a décrit sous forme de schémabloc fonctionnel les parties essentielles de l'ordinateur associé aux détecteurs d.. Les valeurs échantillonnées lues aux sorties des détecteurs d. sont transformées en informations numé- riques grâce à un convertisseur analogique-numérique 18 puis trai- tées dans un préprocesseur 19 effectuant les opérations habituelles de calibration et de transformation logarithmique. Les informations ainsi transformées sont les valeurs d'atténuation linéaires corres- pondantes et ces informations numériques sont mémorisées dans la mémoire M. Pour des raisons qui apparaîtront clairement plus loin, la mémoire M se présente sous forme d'une matrice d'unités de mémorisation comprenant 4 n colonnes et N lignes (nombre de vues Vi) puisque P'acquisition des vues Vi se fait sur un tour complet et avec.: Ap=! - Le remplissage de la mémoire M par les informations repré- sentatives des vues Vi se fait à raison d'une unité de mémorisation sur deux; ces unités sont représentées en grisé sur la figure 1. Le remplissage de la mémoire M, tel que représenté, correspond à un décalage "à gauche" des détecteurs. La mémoire M est associée à des moyens de lecture 20, qui organisent le transfert des infor- mations vers un processeur 21, connu en soi, et dont le r8le est d'appliquer l'algorithme de reconstruction mentionné plus haut. Cependant la lecture d'une ligne de la mémoire M n'est commandée que lorsque toutes les unités de mémorisation de cette ligne qui n'ont pas été "remplies" pendant la phase de mémorisation des vues Vi, auront reçu des informations numériques représentatives des valeurs intermédiaires, calculées par mise en oeuvre de l'algorithme décrit ci-dessus. Pour cela, un processeur spécialisé ou un sous- programme de l'ordinateur effectue les opérations de calcul - d'adresse dans la mémoire M et d'interpolation linéaire avant de commander l'inscription des valeurs intermédiaires résultantes dans les unités de mémorisation laissées provisoirement "vides" lors de l'inscription des vues réelles dans cette même mémoire M. Ce processeur ou ce sous-programme est symbolisé à la figure 1 par le bloc fonctionnel 22 agencé pour calculer les adresses des valeurs d'atténuation linéaires utiles, lire ces valeurs, effectuer le calcul d'interpolation combinant ces valeurs et réinscrire le résultat à l'adresse convenable d'une valeur intermédiaire d'une vue à com- pléter. Il pilote aussi les moyens de lecture 20 (liaison fonctionnelle 23) gérant le transfert des informations depuis la mémoire M vers le processeur de reconstruction 21, le transfert d'une "ligne" de la mémoire n'étant autorisé que lorsque toutes les unités de mémori- sation correspondant aux valeurs intermédiaires ont été mémorisées. REVENDICATIONS 1. Procédé de reconstruction d'image par. tomodensitomètrie, du type consistant à faire tourner un ensemble source-détecteurs dans le plan d'une coupe à imager, ladite source émettant un faisceau de rayonnement pénétrant en forme d'éventail et à relever et mémoriser une succession de vues pour différents angles'5 de rotation de l'ensemble source-détecteurs par rapport à un système d'axes de référence fixe (XOY) dudit plan, lesdits détecteurs étant agencés pour converger vers le foyer (F) de ladite source et chaque vue étant échantillonnée pour des angles Y représentatifs des détecteurs et repérés à partir de la droite (OF) joignant le centre de rotation (0) de l'ensemble source-détecteurs et ledit foyer (F), le décalage angulaire entre deux détecteurs adjacents étant A' constant et le décalage angulaire entre vues étant constant, carac- térisé en ce qu'il consiste, le décalage angulaire entre ladite droite (OF) et le détecteur le plus proche de celle-ci étant réglé à - à relever et mémoriser une série de vues réparties sur un tour complet dudit ensemble source-détecteurs, - à élaborer pour chaque vue un nombre de valeurs intermédiaires représentant un échantillonnage entrelacé avec celui desdits détec- teurs, chaque valeur intermédiaire d'angle Y donné étant obtenue par une combinaison d'au moins deux valeurs d'absorption linéaire déduites de valeurs de mesure prélevées dans des vues situées de part et d'autre d'une vue fictive dont l'image par la transformation T (G$S) = ( entier relatif, lesdites valeurs de mesure desdites vues étant préle- vées sur le détecteur d'angle opposé à celui de ladite valeur intermédiaire recherchée, et - à appliquer un algorithme de reconstruction, connu en soi, compre- nant une convolution et une rétroprojection aux vues ainsi com- plétées par lesdites valeurs intermédiaires. 2. Procédé de reconstruction selon la revendication 1, carac- térisé en ce que ladite combinaison consiste, pour chaque valeur intermédiaire, à faire la moyenne entre les deux valeurs d'absorption linéaire précitées en les pondérant respectivement par l'inverse de l'écart, en projection sur l'axe dans l'espace[e,'{, entre leurs positions respectives résultant de ladite transformation T et la posi- tion de la valeur intermédiaire considérée.