Brevet FR2906375A1 — publié le 20080328

L'invention concerne un système de reconstruction optique d'informations obtenues à partir d'un système d'exploration tomographique. Il concerne toutes les géométries d'acquisition 3D conduisant à des calculs qui peuvent se décomposer en une série d'étapes. L'invention est plus particulièrement utilisable pour l'exploitation des résultats d'exploration d'un scanner utilisé dans le domaine médical. ETAT DE LA TECHNIQUE Un scanner à rayon X est un instrument de diagnostic privilégié pour la détection de diverses pathologies et en particulier les tumeurs. La dernière génération de scanner, tel que celui schématisé en figure 1, est constitué d'une table T sur laquelle est allongé le patient P et d'un système d'acquisition circulaire SC. Les prises de vue sont faites par l'intermédiaire d'un ensemble constitué d'un émetteur de rayon X et d'un capteur associé. Chaque pixel du capteur reçoit une dose de radiation atténuée en fonction de la quantité et du type de tissu traversé. La table avance dans le système d'acquisition et les prises de vues sont faites en continu de sorte que du point de vue du patient, le mouvement est hélicoïdal. Le principe de la reconstruction tomographique est d'utiliser ces acquisitions, que l'on appelle données tomodensitométriques (TDM) brutes pour reconstruire le modèle tridimensionnel du patient. Ce modèle sera utile au (médecin) dans les différentes phases du traitement de la maladie diagnostic, planification, simulations 2906375 2 préopératoires, assistance au geste (thérapeutique) peropératoire. L'opération de base qui nous intéresse dans le cas de la reconstruction d'image tomographique est 5 l'opération de transformation de Radon inverse. Cette opération peut se faire de différentes manières. La méthode d'inversion par recombinaison filtrée est la plus courante (voir le document "Principles of Computerized Tomographic Imaging" de A. C. Kak, M. Slaney, ed. SIAM, 10 2001. ISBN : 0-89871-494-X). Elle consiste à appliquer un filtre rampe à chaque projection avant d'appliquer l'opérateur dit de rétroprojection. L'opérateur de transformée de Radon n'est plus suffisant pour des scanners hélicoïdaux. La complexité de 15 l'acquisition nécessite l'utilisation d'algorithmes plus puissant. L'ASSR (Advanced Single Slice Rebinning) en est un, consistant à faire précéder le calcul de la transformée de Radon inverse par une phase de recombinaison des données. Il fut introduit en 2000 par 20 Kachelriess (voir document "Advanced single-slice rebinning for tilted spiral cont-beam CT" de M. Kachelreiss and T. Fuchs, in Medical Physics, vol. 28, no. 6. Juin 2001). Il existe dans la littérature des traces 25 d'architectures optoélectroniques permettant des reconstructions des images explorées. La méthode Nishimura, publiée en 1978 (voir document "Implementation of the inverse Inverse Radon Transform by optical convolution" de M. Nichimura, D. Psaltis, F. Caimi and D. 30 Casasent, in Optics Communication, vol. 25, num. 3., pp : 2906375 3 301-304. June 1978.), permet le calcul optique de la transformée de Radon inverse par une méthode analytique. A l'époque de sa publication, la faible puissance des outils de calcul justifiait qu'une telle méthode soit 5 développée. De nos jours, le coût, la complexité et le peu de performances de ce système optique tel qu'il est présenté réduit son intérêt. La méthode Lu, décrite dans le brevet US5414623, propose un processeur optoélectronique complet capable de 10 reconstruire des images selon des méthodes algébriques itérative (SIRT et PIRT). Il est constitué de deux coeurs optiques : un projecteur et un rétroprojecteur. Le sinogramme de l'image est d'abord utilisé par le rétroprojecteur pour calculer une première image 15 reconstruite. Le projecteur permet ensuite de calculer le sinogramme de l'image reconstruite. Celui-ci est comparé avec le sinogramme initial. La différence entre Les deux est envoyée sur le rétroprojecteur pour donner la deuxième image reconstruite. Celle-ci est alors 20 additionné à la première. On opère ainsi de manière itérative jusqu'à obtenir un résultat satisfaisant (critère à définir). Ce type de méthodes d'inversion présente l'avantage d'être très robuste au bruit. En contrepartie, 25 la reconstruction nécessite plusieurs itérations et est donc lente. La partie limitante de cet algorithme de traitement, en terme de temps de calcul, est incontestablement le calcul des transformées de Radon inverses. L'algorithme 30 FBP présente une complexité de calcul en MNK où M représente le nombre de lignes de l'image reconstruite, N 2906375 4 le nombre de colonnes et K le nombre de projections dans le sinogramme. De par sa nature, l'algorithme n'est pas propice a un traitement parallèle, ni pipeline. Les techniques d'optimisations classiques ne permettent donc 5 pas de diminuer fortement le temps de calcul. L'invention apporte une solution permettant d'obtenir une diminution de temps de traitement. OBJET DE L'INVENTION L'invention concerne donc un système de 10 reconstruction optique d'informations obtenues à partir d'un système d'exploration hélicoïdal d'un objet. Ce système comporte: - une mémoire d'entrée destinée à stocker l'ensemble d'informations de trames d'explorations, 15 - une mémoire de trames destinée à recevoir des trames d'éléments d'information d'explorations nécessaires à la reconstruction d'un sinogramme parallèle en vue d'obtenir une image d'une coupe de l'objet, -une mémoire d'adresses contenant une table de 20 correspondance entre les adresses d'éléments d'information ou de groupes d'éléments d'information dans la mémoire de trames et les adresses de ces éléments ou de ces groupes d'éléments dans un sinogramme parallèle à obtenir, 25 - une unité de contrôle, - une mémoire de sortie destinée à recevoir successivement des éléments d'informations ou des groupes d'éléments d'informations appartenant tous à un même plan de reconstruction dans différentes trames d'explorations 30 successives et enregistrant ces éléments ou ces groupes 2906375 5 d'éléments d'information à des adresses correspondantes du sinogramme fournies par la mémoire d'adresses, - un système de filtrage et de rétroprojection optique, construisant une ou plusieurs trames images 5 planes à partir d'un sinogramme fourni par la mémoire de sortie, - un dispositif de stockage ou de visualisation recevant lesdites trames images planes du système de filtrage et de rétroprojection et les stockant ou 10 permettant leur visualisation. Selon un mode préféré de réalisation, le système de l'invention comporte une table de coefficients de correction contenant un coefficient de correction pour chaque élément ou pour chaque groupe d'éléments 15 d'informations. A chaque traitement d'un élément d'informations ou d'un groupe d'éléments d'informations, un coefficient correspondant est combiné avec l'adresse de l'élément d'information ou du groupe d'éléments d'informations fourni par la mémoire d'adresses pour 20 corriger la valeur de l'élément d'information ou d'un groupe d'éléments d'informations dans la mémoire de sortie. L'invention est applicable notamment à un système dans lequel lesdites explorations sont hélicoïdales et 25 donc, dans lequel les trames d'explorations fournies sont hélicoïdales. Avantageusement, le traitement de recombinaison permet d'obtenir un sinogramme dans la mémoire de sortie représenté par la formule suivante: 30 R,(n,m,P)=K(m,p).Rhd(i,j,k)=K(m,p).R,(f(n,m,p)) dans laquelle : 2906375 6 - Rä(n,m, p) est le sinogramme parallèle correspondant à la coupe de coordonnées z = zn. Les données dans le sinogramme sont représentées une matrice d'indices m et p, m représentant la variable d'angle et p la variable de 5 position radiale, - Rhe,(i, j,k) représente un élément image d'indice (j,k) sur la i-ème exploration hélicoïdale, - K(m,p) est une constante de correction variant en fonction de (m,p), 10 - f est une fonction de R3 vers R3 permettant de faire la correspondance entre les triplets (i,j,k) et (n,m,p) Les équations de correspondances sont : i=n+m+aresin(K,p) j = K2 p J1_(Kip)2 K3p cos (K4 (m + aresin (K, p))) (m + aresin (K,p)) ù K5 J1_(Kip)2 27z- 20 K1, K2, K3, K4 et K5 sont des constantes dépendantes de la géométrie du système de scanner. Selon une forme de réalisation de l'invention, le système de filtrage et de traitement de rétroprojection réalise un filtrage optique et un 25 traitement de rétroprojection optique. Selon une variante de réalisation, un dispositif de filtrage numérique reçoit les sinogrammes numériques de la mémoire de sortie et les fournit filtrés à un système de rétroprojection optique. 15 k=K2 2906375 7 Le système de filtrage et de traitement de rétroprojection optique peut alors comporter les différents composants optiques suivants alignés selon un même axe: 5 - un premier modulateur spatial de lumière éclairé par une source de lumière, ledit modulateur étant modulé spatialement par un signal de sinogramme de trame provenant de la mémoire de sortie, - une lentille sphérique qui reçoit un faisceau 10 lumineux modulé par le premier modulateur et qui réalise une transformée de Fourier du signal reçu, - un deuxième modulateur spatial de lumière qui est situé dans le plan de la transformée de Fourier de la lentille sphérique et qui module le faisceau qu'il reçoit 15 à l'aide d'un signal fixe dans le temps mais paramétrable, - une première optique asymétrique qui reçoit le faisceau modulé par le deuxième modulateur et qui fournit la transformée de Fourier inverse du signal transmis, 20 - un prisme de Dove qui tourne autour de l'axe du système et qui permet de faire tourner l'image, - une deuxième optique asymétrique qui constitue un objectif dont le grandissement n'est pas le même dans les deux directions, 25 - un dispositif de détection optique matriciel qui fournit une image plane. Selon la variante de réalisation indiquée précédemment, le système de filtrage et de traitement de rétroprojection optique comporte, d'une part, un 30 dispositif de filtrage numérique qui reçoit des signaux de sinogrammes de la mémoire de sortie et qui fournit des 2906375 8 signaux filtrés et, d'autre part, les différents composants optiques suivants alignés selon un même axe: - un troisième modulateur spatial de lumière éclairé par une source de lumière, ce modulateur étant 5 modulé spatialement par un signal filtré fourni. par le dispositif de filtrage numérique de sinogramme de trame provenant de la mémoire de sortie, - un système d'optique sphérique, cylindrique ou asphérique constituant un objectif dont le grandissement 10 est différent dans les deux directions, - un prisme de Dove tournant autour de l'axe du système et permettant de faire tourner l'image, - un dispositif de détection optique matriciel et fournissant une image plane. 15 De préférence, lesdits modulateurs sont des modulateurs linéaires. On peut prévoir que la première optique asymétrique est une lentille cylindrique dont la section droite est parallèle à la direction d'alignement des 20 modulateurs et en ce que la deuxième optique asymétrique est une lentille cylindrique dont la section droite est parallèle à la section droite du prisme, cette lentille cylindrique tournant en synchronisme avec le prisme. Alternativement, on peut également prévoir que la 25 première optique asymétrique est une lentille asphérique et en ce que la deuxième optique asymétrique est une lentille cylindrique dont la section droite est orthogonale à la section droite du prisme. Selon la forme de réalisation de l'invention, on 30 pourra prévoir que la source de lumière émet une lumière cohérente ou incohérente. 2906375 9 BREVE DESCRIPTION DES FIGURES Les différents objets et caractéristiques de l'invention apparaîtront plus clairement dans la description qui va suivre ainsi que dans les figures 5 annexées qui représentent: - la figure 1, une illustration d'un scanner permettant de pratiquer une exploration hélicoïdale d'un corps humain, - la figure 2a, un exemple de réalisation 10 simplifié du système de traitement optique d'informations selon l'invention, - la figure 2b, une variante de réalisation du système de la figure 2a, - la figure 3a, un exemple de réalisation plus 15 complet du système de traitement optique selon l'invention comportant des moyens de corrections, - la figure 3b, une variante de réalisation du système de la figure 3a, 20 - la figure 4, un exemple de réalisation d'un noyau optique des systèmes des figures 2a et 3a, la figure 5, un exemple de réalisation d'un noyau optique des systèmes des figures 2b et 25 3b - la figure 6, une variante de réalisation du noyau optique de la figure 4. DESCRIPTION DETAILLEE 30 En se reportant à la figure 2a, on va donc tout d'abord décrire un exemple de réalisation simplifié du 2906375 10 système de traitement d'informations selon l'invention. Ce système comporte un noyau de recombinaison RC et un noyau de traitement de rétroprojection BPC. Un scanner SC fournit de façon continue des 5 résultats d'explorations à une mémoire de stockage d'entrée Ml. A titre d'exemple, le scanner SC fournit des résultats d'explorations hélicoïdales, mais ces explorations pourraient être d'un autre type sans sortir du cadre l'invention. 10 Le noyau de recombinaison définit la coupe de l'objet que l'on cherche à reconstruire. Il s'agit donc de construire un sinogramme d'une coupe plane d'un objet à partir de plusieurs explorations hélicoïdales. On désigne par sinogramme un ensemble d'informations 15 numériques représentant une transformée de Radon d'explorations parallèles qui permettront d'obtenir une image d'une coupe plane transversale de l'objet. De préférence, on va construire le sinogramme d'une coupe dont le plan est orthogonal à l'axe de l'objet. 20 Mais cela n'est pas obligatoire, on pourrait envisager de construire une image d'une coupe oblique. Le noyaux de recombinaison peut donc être implémenté sous une forme systématique pour chaque case (ou adresse) du sinogramme en géométrie parallèle, une table 25 codée résidente dans la mémoire M3 donne l'adresse correspondante dans le jeu d'acquisition. En d'autres termes, les différentes trames hélicoïdales contenant les informations utiles à la reconstruction de la coupe recherchée sont transférés 30 dans la mémoire M2. Pour chaque élément du sinogramme parallèle, la mémoire M3 contient l'adresse de cet 2906375 11 élément dans les trames hélicoïdales. L'unité de contrôle FPGA inscrit, donc dans une mémoire M4, chaque élément mémoire à une adresse correspondant à sa position dans la coupe plane. En principe, dans un système de scanner 5 parfait, les éléments ayant même adresse dans différentes trames d'exploration hélicoïdales consécutives sont destinés à appartenir à la même image de coupe plane. Le sinogramme d'une coupe plane est ainsi construit dans la mémoire M4. 10 Les données numériques de cette image sont ensuite transmises à un noyau de filtrage numérique et de rétroprojection BPC. Selon l'exemple de réalisation de la figure 2a du système de l'invention, le noyau BPC comporte un système 15 de filtrage et de rétroprojection optique DC1 qui effectue par voie optique le filtrage et le traitement de rétroprojection. Il reçoit de la mémoire de sortie M4 les signaux numériques d'un sinogramme de l'image à reconstruire et fournit une image plane reconstituée à 20 une unité de stockage de coupes reconstruites M6 (ou une unité de visualisation d'image). Selon l'exemple de réalisation de la figure 2b, un système de filtrage numérique de type DSP traite les signaux numériques de la mémoire M4 et les fournit à un 25 système de rétroprojection optique OC2 qui effectue par voie optique le traitement de rétroprojection. Les figures 3a et 3b représentent des exemples de réalisation plus complets dans lesquels il est prévu un système de correction qui permet de corriger l'image de 30 la coupe plane construite lors de son enregistrement dans la mémoire M4. 2906375 12 Pour cela une mémoire M5 contient pour chaque élément d'image un coefficient de correction qui tient compte de la structure du scanner et des conditions d'explorations hélicoïdales. Les informations d'adresses 5 contenues dans la mémoire M3 permettent d'adresser la mémoire M5 en même temps que la mémoire M2. Les valeurs des trames hélicoïdales ainsi que les coeffic:_ents de correction sont transmis simultanément à l'opérateur opl. L'opérateur opl fournit une valeur de chaque élément 10 d'image corrigée par le coefficient de correction. De préférence, l'opérateur effectue une multiplication qui est le seul calcul effectué par le noyau de recombinaison RC. Le sinogramme parallèle reconstitué est donc stocké 15 dans la mémoire tampon M4 qui est de préférence une mémoire d'accès rapide. La construction d'un tel sinogramme peut s'exprimer par la formule suivante: Rii(n,m,p)=K(m,p).R.(i,j,k)=K(m,p).Rher(f(n,m,p) 20 dans laquelle : - R,,(n,m,p) est le sinogramme parallèle correspondant à la coupe de coordonnées Z = Z,,. Les données dans le sinogramme sont représentées une matrice d'indices m et 25 p, m représentant la variable d'angle et p la variable de position radiale. ù Rhe,(i, j,k) représente un élément image d'indice (j,k) sur la i-ème exploration helicoïdale. -K(m,p) est une constante de correction variant en 30 fonction de (m,p). 2906375 13 - f est une fonction de R3 vers R3 permettant de faire la correspondance entre les triplets (i,j,k) et (n,m,p) Les équations de correspondance sont : i=n+m+aresin(K,p) j=KZ p J1_(Kip)2 K3pcos(K4 (m+aresin(K,p))) ù (m+aresin(K,p)) Jlù(K,p) z 27r 10 KI, K2, K3, K4 et K5 sont des constantes dépendantes de la géométrie du système de scanner. Pour simplifier, on peut dire que : - la mémoire Ml contient l'ensemble des données Rh,, (i, j,k) . 15 - la mémoire M2 contient un extrait de Rh, / (i, j,k) pour i allant de i1 à i2. - la mémoire M3 contient la fonction f. La case de coordonnées (n,m,p) de la mémoire M3 contient le triplet (i,j,k) qui servira par la suite à adresser la mémoire M2. 20 En pratique, comme seul i dépend de n, la mémoire M3 ne contient que la fonction réduite : la case de coordonnées (m,p) contient le triplet (i',j,k) avec i'=i-n=m+aresin(K,p). - la mémoire M4 contient la matrice R11(n,m,p) pour un n donné. 25 - la mémoire M5 contient les valeurs de du coefficient de correction. La case de coordonnée (m,p) contient la valeur de K(m, p) . 5 k=K2 2906375 14 Dans la figure 3a, le filtrage et la rétroprojection sont réalisés par le bloc optique OC1 comme dans la figure 2a. La figure 3b représente une variante de réalisation 5 de la figure 3a dans laquelle le filtrage est réalisé par un filtre numérique DSP comme dans la figure 2b. Dans la description qui précède on a considéré que le noyau de recombinaison effectue une recombinaison élément par élément. Cependant, sans sortir du cadre de 10 l'invention, et selon un mode de réalisation préféré de l'invention, le noyau de recombinaison peut également effectuer un traitement de groupe d'éléments par groupe d'éléments. En se reportant à la figure 4, on va maintenant 15 décrire un exemple de réalisation du noyau de rétroprojection BPC. Ce noyau peut-être entièrement réalisé en optique. Il comporte un premier modulateur spatial de lumière Modl éclairé par une source de lumière SOl. Ce premier 20 modulateur est adressé par des signaux électriques provenant de la mémoire M4 des figures 2a à 3b. Le faisceau lumineux modulé par le modulateur Modl est transmis par une lentille sphérique LS1 à un deuxième modulateur spatial de lumière Mod2 qui est adressé par un 25 filtre rampe fixe dans le temps mais qui est paramétrable. Le rôle de la lentille sphérique LS1 est de calculer la transformée de Fourier des projections contenues dans la mémoire M4. Avantageusement, les modulateurs Modl et Mod2 sont 30 des modulateurs ligne . 2906375 15 L'organe de rétroprojection est constitué quant à lui d'une première lentille cylindrique LC1, d'un prisme de Dove SPI et d'une deuxième lentille cylindrique LC2. La première lentille cylindrique LC1 fournit la 5 transformée de Fourier inverse du spectre du sinogramme filtré. Le prisme de Dove SPI permet, lorsqu'il tourne selon la flèche indiquée sur la figure 4, de faire tourner l'image. Un système de lentille, tel que la lentille cylindrique LC2, constitue un objectif dont le 10 grandissement n'est pas le même dans les deux directions et tourne également en synchronisme avec le prisme. Le faisceau de sortie de la lentille cylindrique LC2 est transmis sur un dispositif de détection tel qu'une caméra CCD ou un détecteur d'un autre type, en vue soit 15 du stockage de l'image détectée dans une mémoire, soit de sa transmission à un dispositif de visualisation. La figure 5 représente une variante de réalisation de la figure 4 dans laquelle il est prévu un dispositif de filtrage numérique DSP qui reçoit de la mémoire M4 un 20 sinogramme et qui fournit un sinogramme filtré lequel permet d'adresser le modulateur Mod3. Cette variante de réalisation ne comporte pas alors de modulateur tel que le Mod2 de la figure 4. La figure 6 représente une variante de réalisation 25 du noyau optique de la figure 4. La première lentille cylindrique LC1 a été remplacée par une lentille asphérique LAS1. La deuxième lentille cylindrique LC2 ne tourne pas et est disposée selon un plan orthogonal à la section droite du prisme. 30 Le processeur optique ainsi décrit peut donc fonctionner de deux manières : 2906375 16 • Soit en lumière cohérente, avec la source SOl et le système de filtrage et de projection optique de la figure 4. • Soit en lumière incohérente, avec la source SO2 5 du système de la figure 5 où le système de filtrage est remplacé par un processeur numérique DSP. Le système optique permet alors simplement le calcul de la rétroprojection. En comparaison avec les systèmes existants, le 10 système de l'invention propose un traitement plus rapide. A titre d'exemple, les temps de calcul, pour une reconstruction complète d'une image ont permis d'obtenir les résultats suivants: • Recombinaison d'un sinogramme : 1.8 ms 15 • Filtrage des projections : 10 ms • Reconstruction de l'image: 45 ms Le fonctionnement du système a été validé en simulation. Les points suivants on aussi été étudiés en lumière incohérente: 20 • Discrétisation des niveaux à l'entrée - bruit de quantification : 8 bits suffisent à rendre l'écart quadratique moyen inférieur à 1 %. • Erreur de positionnement instantané du moteur : La précision quant à la position du moteur doit être 25 maîtrisé au dixième de degré près. • Vitesse de rotation du moteur : La précision quant à la vitesse de rotation du moteur doit être maîtrisé à 0,1% près. En tenant compte de tous ces paramètres, l'invention 30 permet d'obtenir une bonne qualité de reconstruction. 2906375 17 L'architecture selon l'invention permet ainsi un gain de temps d'un à deux ordres de grandeur par rapport à une structure purement électronique pour des images de grande taille. 5 REVENDICATIONS 1. Système de reconstruction optique d'informations obtenues à partir d'un système d'exploration d'un objet, caractérisé en ce qu'il comporte: - une mémoire d'entrée (Ml) destinée à stocker l'ensemble d'informations de trames d'explorations, - une mémoire de trames (M2) destinée à recevoir des trames d'éléments d'information d'explorations nécessaires à la reconstruction d'un sinogramme parallèle en vue d'obtenir une image d'une coupe de l'objet, - une mémoire d'adresses (M3) contenant une table de correspondance entre les adresses d'éléments d'information ou de groupes d'éléments d'information dans la mémoire de trames et les adresses de ces éléments ou de ces groupes d'éléments dans un sinogramme parallèle à obtenir, - une unité de contrôle (FPGA), - une mémoire de sortie (M4) destinée à 25 recevoir successivement des éléments d'informations ou des groupes d'éléments d'informations appartenant tous à un même plan de reconstruction dans différentes trames d'explorations successives et enregistrant ces 30 éléments ou ces groupes d'éléments d'information à des adresses correspondantes du 15 20 2906375 19 sinogramme fournies par la mémoire d'adresses (M3) , - un système de filtrage et de rétroprojection optique (CO1, CO2), construisant une ou plusieurs trames images planes à partir d'un sinogramme fourni par la mémoire de sortie, - un dispositif de stockage ou de visualisation (CA, M6) recevant lesdites trames images planes du système de filtrage et de rétroprojection (CO1, CO2) et les stockant ou permettant leur visualisation. 2. Système de reconstruction optique selon la revendication 1, caractérisé en ce qu'il comporte une table de coefficients de correction (M5) contenant un coefficient de correction pour chaque élément ou pour chaque groupe d'éléments d'informations ; à chaque traitement d'un élément d'informations ou d'un groupe d'éléments d'informations, un coefficient correspondant étant combiné avec l'adresse de 1' élément d'information ou du groupe d'éléments d'informations fournie par la mémoire d'adresses (M3) pour corriger la valeur de l'élément d'information ou d'un groupe d'éléments d'informations dans la mémoire de sortie (M4). 3. Système de reconstruction optique selon la revendication 2, caractérisé en ce que lesdites explorations sont hélicoïdales et en ce que lesdites trames d'explorations sont des trames d'explorations hélicoïdales. 5 10 2906375 20 4. Système de reconstruction optique selon la revendication 3, caractérisé en ce que le traitement de recombinaison permet d'obtenir un sinogramme dans la mémoire de sortie (M4) 5 représenté par la formule suivante: e,(r4m,p)=K(m,p).&(i, j,k)=K(14p).&i(f(n,m,P)) dans laquelle: -R,,(n,m,p) est le sinogramme parallèle 10 correspondant à la coupe de coordonnées Z = Z,, . Les données dans le sinogramme sont représentées une matrice d'indices m et p, m représentant la variable d'angle et p la variable de position radiale. 15 - Rhe, (i, j,k) représente un élément image d'indice (j,k) sur la i-ème exploration hélicoïdale. - K(m,p) est une constante de correction variant en fonction de (m,p). - f est une fonction de R3 vers R3 permettant de 20 faire la correspondance entre les triplets (i,j,k) et (n,m,p) Les équations de correspondances sont : i=n+m+aresin(K,p) j=Kz p J1_(Kip)2 K3pcos(K4 (m+aresin(K,p))) (m+aresin(K,p)) ù KS J1ù(Jcip)2 2,r 25 k=K2 2906375 21 K1, K2, K3, K4 et K5 sont des constantes dépendantes de la géométrie du système de scanner. 5. Système de reconstruction optique selon la revendication 4, caractérisé en ce que le système 5 de filtrage et de traitement de rétroprojection réalise un filtrage optique et un traitement de rétroprojection optique. 6. Système de reconstruction optique selon la revendication 4, caractérisé en ce qu'il comporte 10 un dispositif de filtrage numérique (DSP) recevant les sinogrammes numériques de la mémoire de sortie (M4) et les fournissant filtrés à un système de rétroprojection optique (002). 7. Système de reconstruction optique selon la 15 revendication 5, caractérisé en ce que le système de filtrage et de traitement de rétroprojection optique (OC1) comporte les différents composants optiques suivants alignés selon un même axe (XX'): - un premier modulateur spatial de lumière (Modl) 20 éclairé par une source de lumière (SO1), ledit modulateur étant modulé spatialement par un signal de sinogramme de trame provenant de la mémoire de sortie (M4), - une lentille sphérique (LS1) recevant un 25 faisceau lumineux modulé par le premier modulateur (Modl) et réalisant la transformée de Fourier du signal reçu, - un deuxième modulateur spatial de lumière (Mod2) situé dans le plan de la transformée de Fourier de 30 la lentille sphérique et modulant le faisceau 2906375 22 qu'il reçoit à l'aide d'un signal fixe dans le temps mais paramétrable, -une première optique asymétrique (LC1., LSA1) recevant le faisceau modulé par le deuxième 5 modulateur (Mod2) et fournissant la transformée de Fourier inverse du signal transmis, - un prisme de Dove tournant autour de l'axe (XX') du système et permettant de faire tourner l'image, - une deuxième optique asymétrique (LC2, LC3) 10 constituant un objectif dont le grandissement n'est pas le même dans les deux directions, - un dispositif de détection optique matriciel et fournissant une image plane. 8. Système de reconstruction optique selon la 15 revendication 6, caractérisé en ce que le système de filtrage et de traitement de rétroprojection optique comporte, d'une part, un dispositif de filtrage numérique (DSP) recevant des signaux de sinogrammes de la mémoire de sortie (M4) et 20 fournissant des signaux filtrés et, d'autre part, les différents composants optiques suivants alignés selon un même axe (XX'): -un troisième modulateur spatial de lumière (Mod3) éclairé par une source de lumière (SO2), 25 ledit modulateur étant modulé spatialement par un signal filtré fourni par un dispositif de filtrage numérique (DSP)de sinogramme de trame provenant de la mémoire de sortie (M4), - un système d'optique sphérique, cylindrique ou 30 asphérique (LS1, LC1, LC2) constituant un objectif 2906375 23 dont le grandissement est différent dans les deux directions, - un prisme de Dove tournant autour de l'axe (XX') du système et permettant de faire tourner :L'image, 5 - un dispositif de détection optique matriciel et fournissant une image plane. 9. Système de reconstruction optique selon l'une des revendications 7 ou 8, caractérisé en ce que lesdits modulateurs (Modl, Mod2, Mod3) sont des 10 modulateurs linéaires. 10. Système de reconstruction optique selon la revendication 9, caractérisé en ce que la première optique asymétrique est une lentille cylindrique (LC1) dont la section droite est parallèle à la 15 direction d'alignement des modulateurs et en ce que la deuxième optique asymétrique est une lentille cylindrique (LC2) dont la section droite est parallèle à la section droite du prisme (SPI), cette lentille cylindrique tournant en 20 synchronisme avec le prisme (SPI). 11. Système de reconstruction optique selon la revendication 9, caractérisé en ce que la première optique asymétrique est une lentille asphérique (LAS1) et en ce que la deuxième optique 25 asymétrique est une lentille cylindrique dont la section droite est orthogonale à la section droite du prisme (SPI). 12. Système de reconstruction optique selon la revendication 7, caractérisé en ce que la source 30 de lumière (SOI) est une source de lumière cohérente. 2906375 24 13. Système de reconstruction optique selon la revendication 8, caractérisé en ce que la source de lumière (SO2) est une source de lumière incohérente. 5