Brevet FR2655751A1 — publié le 19910614

L'invention concerne les procédés de reconstruction tomographique; elle concerne, plus particulièrement, la reconstruction tomographique à partir d'un nombre limité de projections, en deux et trois dimensions, avec une géométrie d'acquisition qui peut être conique, en éventail ou parallèle. La reconstruction tomographique est effectuée par interpolation des projections de mesure suivant les angles de prise de vue, par filtrage passe-haut des projections de mesure et des projections interpolées, par épandage sur le domaine de reconstruction, suivant la géométrie d'acquisition, des projections de mesure filtrée et de la partie négative des projections interpolées filtrées. Les projections seront fréquemment issues de mesures de la transmission de photons ou de l'émission d'isotopes radioactifs de l'objet à mesurer. Les projections de mesure représentent en transmission et en émission l'intégrale d'une fonction d'atténuation ou densité de radio-éléments suivant un ensemble de rayons. On peut cependant utiliser les techniques de reconstruction tomographique chaque fois que l'on dispose d'un ensemble de mesures représentant l'intégrale d'une fonction suivant des incidences particulières. Il faut rappeler au passage que les méthodes de reconstruction tomographiques se subdivisent en deux groupes. Celles qui utilisent un modèle discret de l'objet de reconstruction: l'objet est découpé en une suite de carrés ou de cubes auquel sont attribués les valeurs de la fonction d'intensité de l'objet, et une discrétisation du processus de mesure: l'intersection entre un rayon et un pixel ou voxel représente le poids affecté à ce pixel ou voxel. L'ensemble de ces poids représente la matrice de projection. L'objet est reconstruit par inversion de cette matrice, généralement de façon itérative. Ces méthodes algébriques sont particulièrement bien adaptées à des conditions difficiles de mesure, ce qui est le cas présent puisque l'on ne dispose pas d'un nombre suffisant de projections de mesures, compte tenu de la résolution attendue sur l'objet reconstruit.Elles nécessitent par contre un nombre important de calculs. Celles qui appartiennent au second groupe de procédés utilisent un modèle continu de l'objet à reconstruire: il existe une relation entre les projections de mesure et la fonction représentant l'objet, appelée théorème coupe-projection, qui indique que la transformée de Fourier des projections mesurées suivant un angle e est la transformée de l'image suivant cette même direction, ce qui permet de reconstruire l'image ou le volume à partir des projections de mesure.Cette méthode, appelée reconstruction analytique, peut revêtir plusieurs formes: l'objet est obtenu soit par inversion du domaine de Fourier préalablement remplit par les mesures, soit par filtrage-épandage des projections de mesure: les projections de mesure sont filtrées par un filtre passe haut (donc qui augmente le contraste des projections), puis les projections sont épandues, c'est à dire que l'intensité en chaque point de la reconstruction est la somme des valeurs des points associés sur chaque projection. Cette méthode de reconstruction est simple de conception et de mise en oeuvre pour un coût de calcul inférieur à la reconstruction algébrique. En trois dimensions la reconstruction analytique peut être soit une approximation de la méthode appliquée en deux dimensions, soit une inversion exacte du processus de génération des mesures. Le procédé de reconstruction tomographique suivant la présente invention peut s'appliquer aux deux groupes et elle est particulièrement efficace pour la méthode de reconstruction par filtrage-épandage. Elle trouve une application particulièrement importante dans le domaine de la reconstruction scanner, lorsque les conditions de mesure ne permettent pas l'acquisition d'un nombre important de projections : la méthode proposée permet alors la reconstruction d'un objet de meilleur qualité que le procédé standard, notamment pour des objets de fort contraste. Le procédé s'applique aussi bien à des reconstructions en deux et en trois dimensions. L'invention part de l'hypothèse qu'on peut considérer les erreurs de reconstruction de 1' objet comme le résultat du sous-échantillonnage angulaire des projections. Ces erreurs de reconstruction se caractérisent généralement par des motifs en étoile autour des parties de fortes intensités de l'objet, appelés aussi artéfacts d'épandage dans le cas particulier de la reconstruction par filtrage-épandage.Pour augmenter le nombre de projections, une interpolation est effectuée suivant la coordonnée angulaire décrivant les angles de vue successifs des mesures; la reconstruction à partir des projections de mesure et des projections interpolées réduira l'amplitude des erreurs décrites ci-dessus, mais introduira une nouvelle source d'erreur due à l'interpolation de mesures sous-échantillonnées, qui se traduira par un élargissement des structures de haute densité contenues dans l'objet; dans la présente invention on se propose de reconstruire un objet à partir des projections interpolées, d'annuler la partie négative de cet objet et de l'ajouter à la reconstruction à partir des projections de mesure, afin de réduire l'amplitude des erreurs de reconstruction, tout en ne contribuant pas à la construction de l'objet et donc sans élargir les structures de forte densité.Du fait de la linéarité de la reconstruction analytique, le procédé peut prendre plusieurs formes suivant la phase du processus de reconstruction où il est inséré: dans le cas particulier de la reconstruction par filtrage épandage, le procédé décrit plus haut peut s'appliquer en n'épandant que la partie négative des projections interpolées filtrées passe haut. De la même manière, d'éventuelles corrections de la valeur moyenne de l'objet reconstruit dont la partie positive a été annulée peuvent tout aussi bien s'appliquer aux projections interpolées filtrées passe-haut et modifiées dans la forme particulière du procédé présentée plus haut. On a déjà proposé (High quality ct image reconstruction from a small number of projections, Saito et al., ICASSP 1988, IEEE, PAGE 1272-1275) d'effectuer la reconstruction tomographique à partir d'un faible nombre de projections, par décomposition de la transformée de Radon (transformée des images en projections, en 2D) en valeurs singulières. Cependant cette méthode ne s'applique théoriquement qu'en géométrie parallèle et ne permet pas la reconstruction d'un objet en trois dimensions à partir d'acquisitions en géométrie conique. De plus elle requiert la décomposition des projections sur une base de fonctions particulières, opération qui engendre un coût de calcul important. L'interpolation des projections a été aussi proposée par plusieurs auteurs (Stark 1981), mais cependant ils ne proposent pas de modifier les projections interpolées, ce qui dans la cas de figure qui nous intéresse (reconstruction à partir d'un faible nombre de projections de structures de support réduit et de contraste important) introduira un élargissement des structures, comme nous l'avons précisé plus haut. La présente invention vise à fournir un procédé de reconstruction ne présentant pas les limitations ci-dessus, tout en ne nécessitant qu'un volume raisonnable de calculs arithmétiques. Dans ce but, l'invention propose notamment un procédé caractérisé en ce que: - on effectue sur les projections de mesure une interpolation suivant les angles de prise de vue, soit par transformée de Fourier de ces projections, introduction de zéros autour de la fréquence de Nyquist, puis transformée de Fourier inverse, soit par tout autre procédé d'interpolation d'un signal numérisé, qui peut être la multiplication des échantillons des projections pris suivant la coordonnée angulaire par la réponse impulsionnelle non-causale d'un filtre passe-bas, afin de calculer la valeur d'un ou de plusieurs échantillons que l'on intercalera entre les échantillons de mesure, - on applique le procédé de reconstruction tomographique sur les projections de mesure, - on ajoute une contribution des projections interpolées à l'objet reconstruit, qui correspond environ à la partie négative de la reconstruction tomographique réalisée à partir de ces projections interpolées. Pour cela on peut appliquer le procédé de reconstruction tomographique sur les projections interpolées, puis ajouter la partie négative de la reconstruction à partir des projections interpolées à la reconstruction à partir des projections de mesure pour reconstruire l'objet final, - on corrige éventuellement l'intensité de l'objet en retranchant le niveau moyen de la partie négative de la reconstruction à partir des projections interpolées. On peut appliquer le procédé suivant une forme équivalente à un autre endroit du processus de reconstruction, du fait de la nature particulière que peut prendre ce processus, et notamment sur les projections filtrées dans le cas de la reconstruction par filtrage-épandage. Dans ce cas le processus de reconstruction est le suivant:: - On interpole les projections de mesure, - on corrige l'amplitude des projections de mesure et interpolées, afin de tenir compte de la géométrie d'acquisition, on filtre les projections par un filtre passe-haut, - on élimine la partie positive des projections interpolées filtrées, - on ramène éventuellement la valeur moyenne des projections interpolées modifiées à zéro, - on épand les projections de mesure et les projections interpolées modifiées sur le volume de reconstruction, - on divise la valeur d'intensité obtenue par épandage par le nombre de projections de mesure. Ces opérations constituent la procédure de reconstruction tomographique analytique d'un objet par filtrage-épandage, qui s'applique en géométries parallèle et éventail. La procédure de reconstruction en trois dimensions peut être soit une approximation de la reconstruction en deux dimensions, soit comporter deux épandages des données issus du volume de Radon rempli par les projections de mesure (Grangeat,"procédé et dispositif d'imagerie tridimensionnelle à partir de mesures bidimentionnelles de l'atténuation de rayonnement", brevet no 8707134). Le procédé peut s'appliquer dans tous ces cas de figure. De multiples variantes du procédé peuvent êtres envisagées: la valeur positive de la reconstruction à partir des données interpolées peut être réduite fortement, sans être annulée.La correction de la valeur moyenne de l'intensité de l'objet reconstruit peut être locale, par exemple plan par plan pour un objet à trois dimensions, ou globale, c'est à dire que la valeur moyenne calculée sur toute la partie négative de l'objet interpolé corrige uniformément l'intensité de l'objet reconstruit. Le procédé appliqué sous sa forme modification des projections possède les mêmes variantes: la valeur moyenne de la partie négative des projections filtrées interpolées peut être définie sur chaque ligne de projection, sur la projection complète qui est par exemple une image en 3 dimensions, ou sur l'ensemble des projections interpolées.Cette correction n'est pas indispensable si l'objet contient peu de structures (objet fil aire par exemple) et dans le cas particulier de certains objets dont on ne recherche qu'une forme binaire (11 objet binaire est l'objet de reconstruction dont la densité est discrétisée en 1 et 0 à partir de la valeur d'un seuil). On observe sur la reconstruction d'un objet à partir des projections de mesure sous-échantillonnées angulairement et des projections interpolées, un élargissement des structures de haute densité et une réduction des erreurs de reconstruction: on peut schématiser ce phénomène de la façon suivante: la partie négative de l'objet reconstruit à partir des projections interpolées réduit l'amplitude des erreurs de reconstruction, tandis que la partie positive de cet objet interpolé participe à la construction de objet et , du fait de la non-consistance des projections interpolées, élargit les structures de haute densité : ce phénomène de filtrage passe-bas s'observe d'autant plus que la structure est éloignée du centre de l'objet et que par conséquent la fonction d'intensité des projections de mesure suivant la coordonnée angulaire possède de hautes fréquences.Le procédé présenté, en éliminant la partie positive de l'objet reconstruit à partir des projections interpolées, supprime l'élargissement des structures de fort contraste: ces structures seront sensiblement aussi bien définies que sur l'objet reconstruit à partir des mesures seules. Par contre le procédé réduit de façon importante les erreurs à l'extérieur de ces structures et engendrées par ces dernières. L'invention sera mieux comprise à la lecture de la description qui suit d'un mode particulier de mise en oeuvre de l'invention, donné à titre d'exemple non limitatif. Dans cet exemple nous utilisons la reconstruction par filtrage épandage étendue à la géométrie 3D conique et nous appliquons le procédé sous sa forme où les projections interpolées sont modifiées, bien que ce soit équivalent du point de vue coût de calcul à l'addition de la partie négative de l'objet reconstruit à partir des projections interpolées à la reconstruction réalisée à partir des mesures.La description se réfère aux dessins qui l'accompagnent, dans lesquels: - La Figure 1 est un graphe destiné à illustrer la reconstruction tomographique avec ajout de la partie négative de la reconstruction réalisée à partir des projections interpolées, à la reconstruction réalisée à partir des mesures , conformément à l'invention, - La figure 2 est un graphe destiné à illustrer la forme particulière que peut prendre le procédé lorsque la reconstruction s'effectue par filtrage épandage des mesures: on épand la partie négative des projections interpolées; - La figure 3 est un schéma synoptique montrant une constitution possible du dispositif de reconstruction. - La figure 4 est un dessin qui montre, à titre d'exemple, comment on modifie les projections interpolées filtrées, avant épandage, avec en (1) une projection interpolée, en (2) cette projection filtrée passe-haut et en (3) la suppression de la partie positive de cette projection filtrée; - Les figures 5A, 5B, 5C et 5D sont des graphiques destinés à montrer les modifications apportées à la reconstruction analytique par le procédé présenté, la Figure SA montrant un plan du fantôme 3D, la Figure 5B la reconstruction analytique à partir de 17 projections de mesure obtenues par simulation, la Figure 5C la reconstruction à partir de 17 projections de mesures simulées et de 16 projections interpolées et la Figure 5D la reconstruction à partir de 17 projections de mesure simulées et de 16 projections interpolées dont la partie positive avant épandage a été supprimée et le niveau de gris moyen ramené à zéro. La reconstruction d'un objet (image ou volume) est effectuée à partir d'un ensemble de projections qui représente l'intégrale de la fonction de densité de cet objet suivant un rayon. Dans l'exemple décrit, nous nous plaçons dans le cas particulier où on cherche à caractériser l'atténuation aux rayons X d'un objet en trois dimensions. La géométrie de mesure est conique, la source de rayons X décrit une trajectoire circulaire. Un dispositif appelé statif, fournit à un instant donné des projections rayon X mesurées suivant cette géométrie, au nombre de P=33 par exemple, avec un décalage angulaire entre chaque projection constant, et des mesures effectuées sur un demi tour( de o à a). De telles projections, répondant à ces critères sont disponibles par exemple auprès du Laboratoire de Traitement du Signal de l'Université de Rennesl. Chaque projection 2D est constituée de pixels de 255 niveaux de gris possibles, rangés suivant les lignes i et les colonnes j. Ces projections sont alors utilisées dans le processus de reconstruction schématisé par la Figure 1. On décrira ce processus en supposant qu'il est mis en oeuvre dans un dispositif de reconstruction ayant la constitution montrée en Figure 3, destiné à être relié à un système hôte, tel qu'un ordinateur assurant la gestion de l'ensemble des dispositifs de prise de projections (statif), de reconstruction de l'objet, de visualisation de traitement et d'archivage de la reconstruction. Le dispositif de reconstruction (Figure 3) comporte alors une mémoire vive principale (4) qui contient un microprogramme de calcul. La mémoire vive (4) constitue également la mémoire de travail. Elle est reliée par un bus de donnée (3) et un bus d'adresse (10) et un accès (2) au système hôte. Le dispositif comprend encore, reliés aux bus (2) et (10), une unité de calcul (5) et un circuit d'aiguillage (7). Ce dernier permet de relier le disque (6) au calculateur (5), soit au contrôleur (8) de commande d'un moniteur de visualisation (9). Le fonctionnement du dispositif est le suivant. Une première opération de chargement de projections de mesure, qui n'est pas modifiée par l'invention, consiste à ranger ces projections sur le disque (6) (mémoire de masse). Le système hôte (non représenté) commande la prise des mesures par le statif (1) et charge les projections de mesure dans le disque (6) par l'intermédiaire de l'accès (2) et des bus (3) et (10), puis il commande le début de la reconstruction à l'unité de calcul (5). Processus d'interpolation des projections de mesure: Cette unité de calcul pour réaliser l'opération d'interpolation des projections, recherche sur le disque l'ensemble des points, au nombre de P, repérés par les mêmes indice i et j suivant les lignes et les colonnes de chaque projection de mesure et vient les ranger en mémoire de travail dans un tampon. Cette opération sera répétée N fois si chaque projection de mesure contient N pixels. La séquence d'échantillons ainsi obtenue est complétée pour obtenir un nombre d'échantillons N égal à une puissance de 2 ; L'unité de calcul réalise une transformée de Fourier rapide sur cette séquence s(n), équivalente à la formulation suivante: où F(t) (complexe) est la transformée de Fourier discrète de s(n). Pour obtenir deux fois plus de projections par exemple, on double le nombre de points de ce signal en insérant autour de la fréquence de Nyquist des zéros: (t) = F(t) O 2 N N 2 2 F(t)=F(t-N) 3N 2 L'unité de calcul réalise alors une transformée de Fourier rapide sur cette nouvelle séquence: où s(n) est la séquence interpolée, stockée dans la mémoire de travail. L'unité de calcul range chaque échantillon impair dans une table à l'indice (i,j), chaque table étant référencée comme une projection interpolée d'angle milieu par rapport aux deux projections de mesure les plus proches, et rangé sur le disque (6). L'unité de calcul répète cette opération d'interpolation pour tous les points des projections, recherche sur le disque (6) les projections interpolées et complètes ces dernières par les échantillons interpolés s(2n#1). Processus de reconstruction du volume à partir des projections de mesure et des projections interpolées: L'algorithme de reconstruction en géométrie 3D conique dans l'exemple que nous donnons pour illustrer le procédé présenté dans l'invention, est du à L.A. Feldkamp et al, "Practical cone-beam algorithm", J. Opt. Soc. Am. A/VOL. 1, No. 6 Juin 1984. L'unité de calcul effectue un traitement de pondération, de filtrage et d'épandage de chaque projection pour reconstruire le volume. L'unité de calcul vient rechercher sur le disque la première projection (qui est une projection de mesure) et range celle-ci en mémoire de travail. l'Unité de calcul réalise une pondération de 1 'amplitude de chaque point de la projection en fonction de la géométrie d'acquisition. La première pondération en géométrie 3D conique est donnée par la rapport entre la distance de la source au centre du support de l'objet ( qui est aussi le centre de rotation de la source) et la distance de la source au point de la projection sur lequel on effectue la pondération.L'unité de calcul effectue une pondération de chaque point de la projection L'unité de calcul vient rechercher une séquence d'échantillons stockée sur le disque, qui représente la réponse impulsionnelle symétrique (cette réponse est pré-calculée), et la range en mémoire vive. Cette réponse impulsionnelle est donnée par la formule suivante: H(t,)= na2 rra -4 H(tk)= #&alpha;2(4k2-1) k#1, #2, #3 ou a désigne l'interval d'échantillonnage des projections. On limite le nombre de termes par exemple à 30.Le filtrage s'effectue suivant les lignes de projection: l'unité de calcul produit un échantillon de projection filtrée en multipliant terme à terme les échantillons d'une ligne de projection par la séquence du filtre dérivateur centrée sur l'échantillon à filtrer, et en additionnant le résultat de chaque multiplication. L'unité de calcul range au fur et à mesure les échantillons de la projection filtrée dans la mémoire de travail. L'unité de calcul effectue ensuite l'épandage de la projection filtrée sur le volume. Le volume est créé en mémoire de travail et la valeur de chaque point, appelé voxel, du maillage (i,j,k) de ce volume discret est positionné à zéro. Pour chaque point du volume (i,j,k) l'unité de calcul détermine en fonction de la géométrie d'acquisition le point sur le plan de projection, donné par les valeurs x et y. Les paramètres correspondant à cette géométrie (distance source/ centre de l'objet, distance source/plan de détection et rapport de la dimension du voxel et du pixel) sont codés dans le microprogramme de reconstruction. Le dernier paramètre de cette géométrie, qui est l'angle d'acquisition de la projection, est déterminée par l'unité de calcul en effectuant la division de n par le nombre de projections total et en multipliant cet incrément par l'ordre de la projection courante. La contribution des valeurs des pixels du plan de projection est déterminé par l'unité de calcul par interpolation linéaire entre les quatre voisins les plus proches du point de projection A(x,y), suivant la formulation: distance du pixel@ A(x,y) Valeur~voxel = ##valeur d'un pixel voisna A(x,y).disance entre deux pixels d'une même ligne La valeur obtenue est pondérée pour tenir compte de l'effet de zoom du faisceau conique, avant d'être attribué au voxel.L'unité de calcul pondère la valeur obtenue par le carré de la distance source/voxel divisée par la distance source/point de projection. Le volume provisoirement reconstruit est stocké dans la mémoire de travail. L'unité de calcul vient chercher une nouvelle projection sur le disque , qui cette fois ci sera une projection interpolée, suivant l'ordre croissant des angles de vue. La nouvelle projection prend la place de la première projection dans la mémoire de travail. L L'unité de calcul réitère le processus de pondération de la projection, ainsi que le processus de filtrage passe-haut de celle-ci, comme décrit plus haut. L'unité de calcul remplace chaque point de la projection filtrée, présente dans la mémoire de travail, soit par zéro ci cette valeur est positive, soit par la même valeur si celle-ci est négative. Pour chaque ligne j de la projection ainsi définie, l'unité de calcul détermine la moyenne des valeurs des pixels en divisant la somme des valeurs de chaque pixel par le nombre de pixels de la ligne. La valeur de chaque pixel de la ligne est alors corrigée en soustrayant la valeur de la moyenne déterminée précédemment. L'opération d'épandage de la projection interpolée filtrée et modifiée suit le même processus que décrit plus haut pour la première projection à ceci près que l'on détermine la contribution de la projection à la valeur du voxel courant (i,j,k), que l'unité de calcul vient cette fois ci rechercher en mémoire de travail (4) la valeur du voxel précédemment calculée, et que cette valeur, ajoutée à la contribution de la projection courante, et attribuée au voxel (i,j,k). Le volume de reconstruction est balayé complètement par l'unité de calcul et rangé en mémoire vive. Le traitement de la troisième projection, une projection de mesure, est analogue à celui de la première projection, excepté que la contribution de cette projection est ajoutée à la valeur de chaque voxel du volume reconstruit. Le traitement de la quatrième projection, une projection interpolée, suit le même processus que celui de la deuxième projection. P projections de mesure et P - 1 projections interpolées sont ainsi traitées. L'unité de calcul divise ensuite la valeur de chaque voxel du volume reconstruit par le nombre de projections de mesure P, et range le volume ainsi calculé sur le disque 6. L'unité de calcul transfère le volume reconstruit sur le microcalculateur de visualisation (8) par l'intermédiaire du contrôleur (7). Ce microcalculateur contrôle un moniteur de visualisation sur laquelle apparait une représentation tri-dimensionnelle du volume reconstruit. L'exemple de visualisation de la reconstruction analytique à partir des projections de mesure et des projections interpolées avec et sans modification, ainsi que de l'objet original illustré en Figure SA - 5C, montre que la reconstruction sans modification des projections interpolées élargit et réduit l'amplitude des structures de forte densité, ce qui n'est plus le cas avec le procédé présenté. La complexité du calcul reste faible: l'épandage des projections interpolées est aussi coûteux que celui des projections de mesure; la phase d'interpolation des projections de mesure par transformée de Fourier rapide est effectuée sur des données de dimensions inférieures d'une unité à la dimension de l'objet, et intervient très peu dans le coût de calcul total de la reconstruction. L'invention est susceptible de nombreuses variantes de réalisation et, en particulier, en ce qui concerne l'atténuation de la valeur positive des projections interpolées filtrées, qui ne nécessitent pas d'être exactement nulles. La correction de la valeur moyenne des projections filtrées modifiées n'est important que lorsque la densité de l'objet possède une valeur moyenne importante, et que l'on recherche une comparaison des niveaux de gris de l'objet de départ et de la reconstruction. On peut utiliser le procédé de reconstruction tomographique à partir d'un nombre limité de projections en dehors de son application aux techniques de rayons X. REVENDICATIONS 1. Procédé de reconstruction tomographique d'images ou de volumes à partir de projections de mesure obtenues par transmission ou émission d'ondes influencées par le passage dans un objet, caractérisé en ce que: - on effectue sur les projections de mesure une interpolation suivant les angles de prise de vue afin de calculer la valeur d'une ou de plusieurs nouvelles projections, - on applique la reconstruction tomographique sur les projections de mesure, - on ajoute une contribution des projections interpolées à l'objet reconstruit, qui correspond environ à la partie négative de la reconstruction tomographique réalisée à partir de ces projections interpolées. 2. Procédé selon la revendication 1, caractérisé en ce que l'élimination de la contribution positive des projections interpolées à l'objet reconstruit peut être obtenue en appliquant la reconstruction tomographique sur les projections interpolées et en ajoutant la partie négative de la reconstruction ainsi obtenue à la reconstruction réalisée à partir des projections de mesure. 3. Procédé selon la revendication 1, caractérisé en ce que l'élimination de la contribution positive des projections interpolées à l'objet reconstruit peut être obtenue en insérant cette suppression de la contribution positive des projections interpolées dans le processus de reconstruction tomographique. 4. Procédé selon les revendications 1 et 3, caractérisé en ce que l'élimination de la contribution positive des projections interpolées à l'objet reconstruit peut être obtenue, pour la méthode de reconstruction par filtrage-épandage, en éliminant la partie positive des projections interpolées filtrées, avant épandage. 5. Procédé selon la revendication 1, caractérisé en ce qu'on peut modifier l'intensité de l'objet reconstruit en retranchant le niveau moyen local de la partie négative de la reconstruction à partir des projections interpolées. Le calcul du niveau moyen local peut être effectué sur un plan de cette reconstruction, la correction s'effectuant sur le plan correspondant de l'objet reconstruit. 6. Procédé selon la revendication 1, caractérisé en ce qu'on peut modifier l'intensité de l'objet reconstruit en retranchant le niveau moyen global de la partie négative de la reconstruction à partir des projections interpolées. La calcul du niveau moyen global est effectué sur l'ensemble des points de cette reconstruction, la correction s'appliquant sur l'ensemble des points de 1 'objet reconstruit. 7. Procédé selon les revendications 1 et 3, caractérisé en ce que, pour la reconstruction par filtrage épandage, on peut modifier l'intensité de l'objet reconstruit en modifiant l'intensité des projections interpolées filtrées dont la partie positive est annulée. On retranche à ces projections la valeur moyenne calculée localement. Cette valeur moyenne peut être calculée sur une ligne d'une projection bidimensionnelle, la correction est alors appliquée sur cette même ligne. 8. Procédé selon les revendications 1 et 3, caractérisé en ce que, pour la reconstruction par filtrage-épandage, on peut modifier l'intensité de l'objet reconstruit en modifiant l'intensité des projections interpolées filtrées dont la partie positive est annulée. On retranche à ces projections la valeur moyenne calculée de façon globale. Cette valeur moyenne est calculée sur l'ensemble des points de chaque projection mono ou bidimensionnelle. La correction est alors appliquée sur l'ensemble des points de ces mêmes projections. 9. Procédé selon la revendication 1, caractérisé en ce que la méthode permettant la mise en oeuvre de la dite interpolation angulaire sur les projections de mesure, est décrite par les étapes suivantes: - on construit une séquence à partir des points repérés suivant les mêmes coordonnées sur chaque projection discrétisée, - on effectue une transformée de Fourier discrète de cette séquence, - on insère des zéros sur cette séquence transformée, autour de la fréquence de Nyquist, - on effectue une transformée de Fourier discrète inverse de cette nouvelle séquence, - on construit les projections interpolées à partir des échantillons qui ne correspondent pas aux mesures dans la séquence obtenue.