Brevet FR2798804A1 — publié le 20010323

I Procédé, appareil pour reconstruire des données d'image acquises par un système d'imagerie à rayons X par tomosynthèse La présente invention concerne la tomosynthèse, et plus par- ticulièrement un procédé et un appareil qui utilisent un système de to- mosynthèse pour saisir des images en projection 2D d'un objet et qui reconstruisent ensuite une représentation 3D de l'objet en utilisant un algorithme de reconstruction de tomographie informatisée volumétri- que à faisceau conique. La tomosynthèse circulaire permet de construire une image tridimensionnelle (3D) d'un objet à partir d'un ensemble fini d'images radiographiques de projection bidimensionnelle (2D). La figure 1 illustre la géométrie d'un système typique de tomosynthèse circulaire. Le système comprend une source 1 de rayons X et un détecteur 2 2D de rayons X, qui est normalement un écran. Pendant l'acquisition des données, la source 1 de rayons X comme le détecteur 2 tournent sui- vant des trajectoires circulaires et le détecteur 2 saisit un ensemble de radiographies de projection, à des emplacements de source discrets. Les trajectoires circulaires sont dans des plans qui sont parallèles en- tre eux, et les déplacements de la source 1 et du détecteur 2 se font en sens contraire suivant ces trajectoires. Une fois obtenues les radiographies de projection, on les nu- mérise, puis elles subissent un mouvement de translation les unes par rapport aux autres pour les superposer de telle manière que les images de structures dans le plan de tomosynthèse se recouvrent exactement. Les images des structures situées à l'extérieur du plan de tomosynthèse ne se recouvrent pas exactement, ce qui se traduit par un rendu flou de ces structures, dépendant de la profondeur. En faisant varier la quanti- té de translation relative des radiographies de projection, on peut faire varier l'emplacement du plan de tomosynthèse. A chaque fois que l'on modifie le plan de tomosynthèse, les données d'image correspondant aux structures chevauchantes sont superposées et l'on obtient une ima- ge 2D de la structure dans le plan de tomosynthèse. Une fois obtenu un jeu complet d'images 2D de l'objet, on génère une image 3D de l'objet à partir du jeu d'images 2D. La tomosynthèse circulaire peut être très utile pour les appli- cations cliniques qui nécessitent une résolution spatiale élevée dans le sens de la coupe coronale à travers le patient, en particulier parce qu'il est habituellement difficile d'obtenir une haute résolution dans le sens de la coupe coronale avec la tomographie informatisée (TI) à rayons X. De plus, l'appareil de tomosynthèse est généralement simple et peu coûteux comparé à la TI à rayons X. Cependant, un défaut des systè- mes de tomosynthèse circulaire est que la qualité d'image des images 3D reconstruites est habituellement limitée en raison du flou des struc- tures qui se trouvent à l'extérieur du plan de tomosynthèse. De nom- breux algorithmes de reconstruction ont été développés pour améliorer la qualité d'image et pour réduire ce flou hors-plan. Par exemple, on a utilisé dans ce but des algorithmes basés sur des techniques d'inver- sion de matrices. Une technique bien connue pour générer des reconstructions 3D utilisant des images de projection en 2D est celle appelée tomogra- phie informatisée volumétrique à faisceau conique (TIV à faisceau co- nique). La figure 2 illustre la géométrie d'un système TIV à faisceau conique. La source 4 de rayons X projette un c6ne de rayons X sur l'objet 5. Les rayons X tombent sur un détecteur numérique 6 qui nu- mérise les données. Le détecteur 6 est dans un plan qui n'est pas paral- lèle au plan dans lequel se trouve la source 4. Au lieu de cela, le dé- tecteur 6 est situé dans un plan qui est orthogonal au plan dans lequel se trouve la source. Comme dans le cas de la tomosynthèse circulaire, un déplace- ment relatif est généré entre la géométrie détecteur 6/source 4 et l'ob- jet 5. Ceci peut être réalisé en faisant tourner la source 4 et le détec- teur 6 ou en faisant tourner l'objet 5. Les images de l'objet sont re- cueillies en certaines positions de la source 4 et du détecteur 6 ou en certaines positions de l'objet 5, en fonction de la manière dont est gé- néré le déplacement relatif. Au fur et à mesure que sont générées des cartes 2D de l'intensité d'image mesurée par le détecteur 6, ces cartes 2D sont filtrées en utilisant des techniques connues de filtrage d'ima- ges puis sont reconstruites en utilisant un ou plusieurs algorithmes connus pour obtenir une représentation 3D de l'objet. Les images 3D reconstruites en utilisant des systèmes TIV à faisceau conique ont une meilleure qualité d'image que les images re- construites obtenues en utilisant des systèmes typiques de tomosynthè- se circulaire. Il serait souhaitable d'améliorer la qualité des images obtenues en utilisant des systèmes de tomosynthèse circulaire, plutôt que de remplacer les systèmes de tomosynthèse circulaire par des sys- tèmes TIV à faisceau conique. Remplacer les systèmes de tomosynthè- se circulaire par des systèmes TIV à faisceau conique peut ne même pas être possible dans beaucoup de cas, en raison des coûts et d'autres facteurs. Par conséquent, il existe un besoin pour un procédé et un ap- pareil qui améliorent la qualité des images reconstruites à partir de données d'images acquises en utilisant des systèmes de tomosynthèse circulaire. La présente invention propose un procédé et un appareil qui utilisent un système de tomosynthèse circulaire pour recueillir des données radiographiques de projections de rayons X en 2D d'un objet analysé. Les données recueillies sont ensuite transformées pour avoir une forme appropriée à la reconstruction par tomographie informatisée volumétrique à faisceau conique (TIV à faisceau conique). Une fois les données transformées, on utilise généralement un algorithme de reconstruction de TIV à faisceau conique pour recon- struire une image 3D de l'objet. La présente invention propose un appareil pour générer une image tridimensionnelle (3D) à partir d'une pluralité d'images de pro- jection à deux dimensions (2D). Cet appareil comprend un système de tomosynthèse circulaire, comprenant une source de rayons X et un dé- tecteur, la source de rayons X projetant des rayons X vers un objet à analyser, le détecteur détectant des rayons X émis par la source de rayons X, le détecteur générant des signaux électriques en réponse aux rayons X qui le frappent, les signaux électriques correspondant aux données d'image acquises par le système de tomosynthèse circulaire; et un dispositif logique configuré pour recevoir les signaux électriques générés par le détecteur, ledit dispositif logique étant configuré pour transformer les données d'image de tomosynthèse circulaire en don- nées d'image de tomographie informatisée volumétrique à faisceau co- nique (TIV à faisceau conique). De préférence, le dispositif logique est un ordinateur, ledit ordinateur étant configuré pour exécuter un algorithme de transforma- tion, l'algorithme de transformation transformant les données d'image de tomosynthèse circulaire en données d'image TIV à faisceau coni- que, ledit ordinateur étant en outre configuré pour exécuter un algo- rithme de reconstruction TIV à faisceau conique qui opère sur les don- nées d'image transformées pour générer l'image 3D. Le détecteur du système de tomosynthèse circulaire est dispo- sé dans un plan qui est parallèle à un plan de trajectoire de la source de rayons X du système de tomosynthèse circulaire. L'algorithme de transformation transforme les données d'image recueillies par le détec- teur du système de tomosynthèse circulaire en un plan de détecteur vertical virtuel qui est orthogonal au plan dans lequel se trouve le dé- tecteur de tomosynthèse. Le plan de détecteur vertical virtuel corres- pond à un plan dans lequel un détecteur de système TIV à faisceau co- nique serait disposé par rapport à une source de rayons X de système TIV à faisceau conique. De préférence, lorsque les données d'image sont transformées en plan de détecteur vertical virtuel, une kème rangée du détecteur du système de tomosynthèse circulaire est transformée en une rangée du plan de détecteur vertical virtuel à une hauteur Zvk le long d'un axe des z. Une fois que les données d'image ont été transformées en données d'image dans le plan de détecteur vertical virtuel, les données d'image transformées peuvent être traitées par l'ordinateur pour corri- ger toute non-uniformité de l'espacement des détecteurs élémentaires dans les directions suivant les axes y et z des détecteurs élémentaires du plan de détecteur vertical virtuel. De préférence, l'ordinateur traite les données d'image trans- formées pour corriger toute non-uniformité de l'espacement des détec- teurs élémentaires en exécutant un algorithme de correction. L'algorithme de correction peut déterminer par interpolation une valeur d'intensité en utilisant les valeurs d'intensité du détecteur élémentaire intersecté et de la pluralité de détecteurs élémentaires en- tourant le détecteur élémentaire intersecté pour obtenir une valeur d'intensité moyenne qui est attribuée au détecteur élémentaire corres- pondant dans le plan de détecteur vertical virtuel. De préférence, l'algorithme de correction est constitué de l'algorithme de reconstruction TIV à faisceau conique. Ce dernier peut être un algorithme de Feldkamp. De préférence, l'algorithme de reconstruction TIV à faisceau conique est un algorithme de Feldkamp modifié, l'algorithme de Feld- kamp modifié comprenant des segments de code de filtrage et de rétro- projection qui ont été modifiés pour réaliser les fonctions de l'algo- rithme de correction. La présente invention propose également un procédé pour gé- nérer une image 3D à partir d'une pluralité d'images de projection 2D. Ce procédé comprend les étapes consistant à: utiliser un système de tomosynthèse circulaire, tel que précé- demment décrit, pour acquérir les images de projection 2D; transformer les données d'image de tomosynthèse circulaire en données d'image de tomographie informatisée volumétrique à fais- ceau conique (TIV à faisceau conique); et exécuter une reconstruction de TIV à faisceau conique sur les données d'image transformées pour générer l'image 3D. De préférence, après l'étape de transformation des données d'image et avant l'étape d'exécution d'une reconstruction TIV à fais- ceau conique, on exécute une étape de correction sur les données d'image transformées, l'étape de correction traitant les données d'ima- ge transformées pour corriger toute non-uniformité de l'espacement des détecteurs élémentaires dans les directions suivant les axes y et z des détecteurs élémentaires du plan de détecteur vertical virtuel. De préférence, l'étape de correction corrige toute non-unifor- mité de l'espacement des détecteurs élémentaires en exécutant les éta- pes décrites plus loin. De préférence, l'étape de correction est exécutée au cours de l'étape d'exécution d'une reconstruction TIV à faisceau conique. L'algorithme de reconstruction TIV à faisceau conique peut être un algorithme de Feldkamp. La présente invention propose enfin un programme d'ordina- teur pour générer une image 3D à partir d'une pluralité d'images de projection 2D, les images de projection 2D étant acquises par un systè- me de tomosynthèse circulaire tel que précédemment décrit, le pro- gramme d'ordinateur étant enregistré sur un support lisible par ordina- teur. De préférence, le programme comprend un premier segment de code transformant les données d'image de tomosynthèse circulaire en données d'image de tomographie infor- matisée volumétrique à faisceau conique (TIV à faisceau conique); et un deuxième segment de code exécutant une reconstruction TIV à faisceau conique sur les données d'image transformées pour gé- nérer l'image 3D. Le premier segment de code transforme une kème rangée du détecteur du système de tomosynthèse circulaire en une rangée du plan de détecteur vertical virtuel à une hauteur Zv,k le long d'un axe des z. Le programme peut comprendre en outre un troisième seg- ment de code, le troisième segment de code exécutant une opération de correction sur les données d'image transformées, l'opération de correc- tion traitant les données d'image transformées pour corriger toute non- uniformité de l'espacement des détecteurs élémentaires dans les direc- tions suivant les axes y et z des détecteurs élémentaires du plan de dé- tecteur vertical virtuel. De préférence, le troisième segment de code corrige toute non-uniformité de l'espacement des détecteurs élémentaires en exécu- tant les étapes décrites plus loin. Le troisième segment de code peut déterminer par interpola- tion la valeur d'intensité à attribuer en utilisant les valeurs d'intensité du détecteur élémentaire intersecté et de la pluralité de détecteurs élé- mentaires entourant le détecteur élémentaire intersecté pour obtenir une valeur d'intensité moyenne à attribuer au détecteur élémentaire correspondant dans le plan de détecteur vertical virtuel. La présente invention sera mieux comprise à la lecture de la description détaillée suivante faite en référence aux dessins d'accom- pagnement, dans lesquels: la figure 1 est une illustration de la géométrie d'un système de tomosynthèse typique de l'art antérieur; la figure 2 est une illustration de la géométrie d'un système typique de tomographie informatisée volumétrique à faisceau conique (TIV à faisceau conique) de l'art antérieur; la figure 3 est un schéma fonctionnel de l'appareil de l'inven- tion en conformité avec un mode de réalisation; la figure 4 est un organigramme illustrant le procédé de la présente invention en conformité avec un mode de réalisation de l'in- vention pour reconstruire une représentation 3D d'un objet analysé par l'appareil représenté en figure 3; la figure 5 illustre le processus de transformation réalisé par le procédé représenté par l'organigramme de la figure 4 pour transfor- mer les données recueillies par un détecteur horizontal d'un système de tomosynthèse circulaire en des données de plan de détecteur verti- cal virtuel. La figure 3 est un schéma fonctionnel de l'appareil de l'in- vention en conformité avec un mode de réalisation de l'invention. L'appareil est un système de tomosynthèse circulaire tel que celui évo- qué plus haut en référence à la figure 1. Toutefois, contrairement aux systèmes typiques de tomosynthèse circulaire, l'algorithme de recon- struction 14 exécuté par un dispositif logique 10 est un algorithme de reconstruction de TIV à faisceau conique. L'algorithme de transforma- tion 13 exécuté par le dispositif logique 10 transforme les données d'images acquises par le système de tomosynthèse circulaire pour leur donner une forme qui est adaptée à la reconstruction par application de l'algorithme 14 de TIV à faisceau conique. L'algorithme de transfor- mation 13 et l'algorithme de reconstruction 14 seront décrits plus loin en détail, en référence aux figures 4 et 5. Le dispositif logique 10 comprend typiquement un ordinateur qui est en communication avec un dispositif d'affichage 1 1 et avec un composant 12 de stockage de données du système. Le composant 12 de stockage de données du système stocke les données utilisées par l'ordi- nateur 10 pour réaliser la reconstruction. Le composant 12 de stockage de données du système stocke également les données d'image qui ont été acquises par le système de tomosynthèse. En outre, le composant 12 de stockage de données du système stocke le code correspondant aux algorithmes 13 et 14 et les données d'image qui ont été reconstrui- tes par application de ces algorithmes. L'ordinateur 10 est capable de faire afficher les images reconstruites sur le dispositif d'affichage 11. L'appareil de l'invention comprend un détecteur 15, qui est de préférence un détecteur numérique ou intensificateur d'image. L'appa- reil de l'invention comprend une source 19 de rayons X pour émettre des rayons X à travers l'objet à analyser 18, qui peut être un patient humain, par exemple. Le détecteur 15 reçoit des rayons X qui passent à travers l'objet 18 et génère des signaux de tension proportionnels à l'intensité des rayons X. Le détecteur 15 et la source 19 sont respecti- vement couplés électriquement à des unités de commande 16 et 17. Les unités de commande 16 et 17 sont couplées électriquement à l'ordina- teur 10. L'ordinateur 10 envoie des instructions aux unités de comman- de 16 et 17. En réponse à ces instructions provenant de l'ordinateur 10, les unités de commande 16 et 17 déplacent le détecteur 15 et la source 19 suivant des trajectoires circulaires, de la manière décrite plus haut en référence à la figure 1. La manière dont les unités de commande 16 et 17 génèrent le déplacement du détecteur 15 et de la source 19 est bien connue des personnes du métier. Comme établi ci-dessus, l'appareil de l'invention utilise de préférence un système typique de tomosynthèse circulaire pour acquérir des données d'image à rayons X. Par conséquent, pour plus de brièveté, il ne sera pas donné ici de description détaillée des modes de fonctionnement du système de tomosynthèse circulaire et du processus d'acquisition de données. Les instructions reçues par les unités de commande 16 et 17 instruisent les unités de commande 16 et 17 sur la manière dont l'ac- quisition de données doit être réalisée. En des positions particulières de la source 19 sur la trajectoire circulaire, des échantillons d'image sont acquis par le détecteur 15. Le détecteur 15 comprend, par exem- ple, un détecteur numérique ou intensificateur d'image. Typiquement, le détecteur 15 génère des signaux de tension analogiques en réponse aux rayons X frappant ce dernier, et si le détecteur 15 est un intensifi- cateur d'image, un convertisseur analogique-numérique (non représen- té) convertit ces signaux de tension analogiques en signaux de tension numériques. Les signaux de tension numériques sont ensuite envoyés à l'ordinateur 10 pour être traités, comme décrit en détail ci-dessous. Il faut remarquer qu'au lieu de déplacer le détecteur 15 et la source 19, on peut déplacer l'objet à analyser 18 afin de générer le mouvement de déplacement relatif entre le détecteur 15, la source 19 et l'objet 18. Par ailleurs, les personnes de l'art comprendront qu'une grande variété de mécanismes appropriés est disponible pour générer le mouvement du détecteur 15 et de la source 19. De même, les person- nes de l'art comprendront qu'une grande variété de mécanismes appro- priés est disponible pour générer le mouvement de l'objet à analyser 18. Bien que la figure 3 illustre des unités de commande distinc- tes 16 et 17 pour commander le déplacement respectivement du détec- teur 15 et de la source 19, il n'est pas nécessaire d'utiliser deux unités de commande. Une unité de commande unique peut être utilisée pour commander le mouvement soit du détecteur 15 soit de la source 19, et le détecteur 15 et la source 19 peuvent être couplés mécaniquement par un mécanisme de liaison (non représenté). Ce mécanisme de liai- son maintient la synchronisation du mouvement du détecteur 15 et de la source 19. Il est à noter également que l'invention de système de tomo- synthèse circulaire n'est limitée à aucun type particulier ni à aucune configuration particulière. Par ailleurs, la présente invention n'est li- mitée à aucun ordinateur particulier pour exécuter les tâches de traite- ment de l'invention. Le terme "ordinateur", tel qu'il est utilisé ici, est destiné à désigner toute machine capable d'exécuter les calculs néces- saires pour exécuter les tâches de l'invention. Le terme "ordinateur" est destiné à désigner toute machine capable d'accepter une entrée structurée et de traiter cette entrée selon des règles prescrites pour produire une sortie. L'ordinateur 10 comprend un microprocesseur (non représen- té) qui exécute toutes les tâches de calcul. Le composant 12 de stocka- ge de données du système peut être un composé intégré ou séparé, c'est-à-dire qu'il peut se trouver à l'intérieur du circuit intégré compre- nant le microprocesseur ou bien à l'extérieur de ce circuit intégré. Il est à noter que des dispositifs de calcul autres qu'un microprocesseur peuvent exécuter les calculs de l'invention, comme, par exemple, un circuit intégré personnalisé (ASIC) ou un processeur de traitement nu- mérique du signal (DSP). En outre, bien que le procédé de la présente invention est de préférence exécuté par une combinaison matériel/logiciel (c'est-à-dire, le microprocesseur exécutant les algorithmes de transformation et de reconstruction), il est à noter que le procédé de la présente invention peut être entièrement exécuté par du matériel ("hardware"), comme le comprendront les personnes du métier. Le procédé de l'invention conforme au mode de réalisation préféré est illustré par l'organigramme de la figure 4. Le bloc 21 re- présente le processus d'acquisition de données d'image. Une fois les données d'image acquises, elles sont stockées dans le composant 12 de stockage de données du système, comme représenté au niveau du bloc 23. L'ordinateur 10 extrait ensuite les données d'image du composant 12 de stockage de données, image par image, comme indiqué dans le bloc 25. En même temps qu'il extrait les données du composant 12 de stockage de données, l'ordinateur 10 transforme les données pour leur donner une forme adaptée à la reconstruction TIV à faisceau conique, comme indiqué dans le bloc 27. Pendant le processus de transformation représenté par le bloc 27, les données sont transformées, du plan de détecteur horizontal de la géométrie de tomosynthèse circulaire, elles deviennent des données d'un plan de détecteur vertical virtuel de la géométrie TIV à faisceau conique. Au cours de ce processus de transformation, on corrige égale- ment les données pour la non-uniformité de l'espacement des détec- teurs élémentaires dans le plan vertical de détecteur. L'ordinateur 10 réalise alors une reconstruction TIV à faisceau conique, comme indi- qué par le bloc 29. L'étape de transformation des données et de correc- tion des données transformées pour la non-uniformité de l'espacement des détecteurs élémentaires dans le plan vertical virtuel de détecteur sera décrite plus loin en détail en se référant à la figure 5. La reconstruction TIV à faisceau conique est bien connue dans l'industrie. Un algorithme bien connu pour réaliser une recon- struction à faisceau conique est l'algorithme de reconstruction à fais- ceau conique de Feldkamp. D'autres types d'algorithmes comportant des modifications de l'algorithme de Feldkamp avec segments de code de filtrage et de rétroprojection, sont également connus dans le métier et adaptés pour réaliser une reconstruction TIV à faisceau conique re- présentée par le bloc 29 de la figure 4. L'algorithme de Feldkamp est décrit dans un article intitulé "Practical Cone-Beam Algorithm" (algo- rithme de faisceau conique pratique), de L. Feldkamp, L. Davis et J. Kress, J. Opt. Soc. Am., A/Vol. 1, n 6, Juin 1984. * Cet algorithme est également décrit dans un manuel intitulé "Principles of Computerized Tomographic Imaging" (principes d'ima- gerie tomographique informatisée), d'A. Kak et M. Slaney. Les person- nes du métier comprendront que l'invention n'est limitée à aucun algo- rithme particulier pour réaliser la reconstruction TIV à faisceau coni- que. Etant donné que l'algorithme de Feldkamp et autres algorithmes de reconstruction TIV sont bien connus dans l'industrie, la manière dont ce type d'algorithme de reconstruction TIV est appliqué aux don- nées recueillies selon l'invention ne sera pas décrite ici pour plus de brièveté. Les personnes du métier comprendront comment on applique ce type d'algorithme. Les systèmes typiques de TIV à faisceau conique utilisent une source de rayons X qui projette un cône de rayons X à travers l'objet à analyser. Afin de capter ces rayons X, le détecteur, dans un système typique de TIV à faisceau conique, est placé dans un plan qui est orthogonal au plan dans lequel un détecteur est placé dans un sys- tème de tomosynthèse circulaire. Les algorithmes de reconstruction TIV à faisceau conique actuellement disponibles sont conçus pour ma- nipuler des données qui ont été acquises par un détecteur positionné de cette manière. Le procédé de l'invention permet à des données acqui- ses par un détecteur de tomosynthèse circulaire d'être transformées en données qui semblent avoir été recueillies par un détecteur TIV à fais- ceau conique. Ce procédé de transformation, qui comprend le proces- sus de correction des données transformées pour la non-uniformité de l'espacement des détecteurs élémentaires dans le plan de détecteur ver- tical virtuel, va maintenant être décrit en référence à la figure 5. A titre d'exemple, et sans caractère limitatif, la source 19 de rayons X est en une position particulière (x,y,z) = (x0,0,0) sur la tra- jectoire circulaire 36. Le détecteur 15 est placé dans un plan qui est parallèle au plan dans lequel on fait tourner la source 19. Le détecteur se trouve à une hauteur z = zt. Le détecteur 15 recueille les données d'intensité d'émission de l'objet balayé (non représenté). L'invention transforme les données recueillies par le détecteur horizontal 15 en données d'un plan de détecteur vertical virtuel, qui est le plan désigné par la référence 35 en figure 5. Chaque rangée 38 du détecteur hori- zontal 15 est transformée en rangée 39 du plan 35 de détecteur vertical virtuel. Les rangées 38 du détecteur horizontal 15 et les rangées 39 du plan 35 de détecteur vertical virtuel sont parallèles à l'axe des y. Pour un détecteur horizontal 15 ayant un espacement At des détecteurs élémentaires, la kême rangée du détecteur horizontal 15 est transformée enune rangée du plan 35 de détecteur vertical virtuel à la hauteur Zv,k, définie par l'équation Zv,k = ztD/(D - kAt), Equation 1 o D est la distance séparant l'emplacement de la source 19 de rayons X du plan 35 de détecteur vertical virtuel suivant l'axe des x. L'espa- cement vertical (direction z), Azk, entre la kème et la (k-l)ème rangées sur le détecteur vertical, est défini par l'équation: Az,k = Zv,k - Zv,k1 Equation 2 = (zv,k/(D-(k-1)At)) At Equation 3 L'espacement des détecteurs élémentaires, Ay,k, des données de rangée (dans la direction y) dans le plan 35 de détecteur vertical virtuel est défini par l'équation Ay,k = (Zv,k/Zt) Equation 4 Il est à noter que l'espacement des détecteurs élémentaires du détecteur horizontal 15 n'a pas besoin d'être égal dans les directions x et y, comme le comprendront les personnes du métier. Une fois les données transformées dans le plan 35 de détec- teur vertical virtuel en partant du détecteur horizontal 15 selon ces équations, on corrige les données pour la non-uniformité de l'espace- ment des détecteurs élémentaires dans le plan 35 de détecteur vertical virtuel. On peut alors utiliser des algorithmes de reconstruction TIV à faisceau conique, tels que, par exemple, l'algorithme de Feldkamp, pour reconstruire une image 3D de l'objet. Lorsque les données du détecteur horizontal 15 sont transfor- mées en données du plan 35 de détecteur vertical virtuel, l'espacement des détecteurs élémentaires du plan 35 de détecteur vertical virtuel n'est pas parfaitement uniforme. Les algorithmes de reconstruction TIV actuel- lement utilisés sont conçus pour manipuler des données recueillies par des détecteurs ayant un espacement uniforme des détecteurs élémen- taires. Par conséquent, les données d'image doivent être traitées pour compenser la non-uniformité de l'espacement des détecteurs élémentai- res avant d'utiliser un algorithme de reconstruction TIV "standard", tel que l'algorithme de Feldkamp, pour reconstruire l'image 3D. L'ordinateur 10 exécute la partie de l'étape 27 destinée à cor- riger la non-uniformité de l'espacement des détecteurs élémentaires en traçant une ligne à partir de chaque détecteur élémentaire du plan 35 de détecteur vertical virtuel correspondant aux données transformées, jusqu'à la source 19 en passant à travers le détecteur horizontal 15. L'ordinateur 10 détermine alors quel détecteur élémentaire du détec- teur horizontal 15 est coupé par la ligne. Les valeurs d'intensité du dé- tecteur élémentaire intersecté et d'une pluralité de détecteurs élémen- taires entourant le détecteur élémentaire intersecté sont ensuite utili- sées par l'ordinateur 10 pour déterminer par interpolation une valeur d'intensité à attribuer au détecteur élémentaire correspondant du plan de détecteur vertical virtuel. Le processus d'interpolation peut être aussi simple que d'effec- tuer la moyenne des valeurs d'intensité du groupe de détecteurs élémentaires pour obtenir une valeur d'intensité moyenne, qui est ensuite attribuée au détecteur élémentaire correspondant dans le plan 35 de détecteur vertical virtuel. Les personnes du métier comprendront que plusieurs méthodes d'interpolation différentes peuvent être appropriées. Il est à noter que les tâches représentées par le bloc 27 de la figure 4 peuvent être exécutées à l'aide d'un code qui est distinct de l'algorithme de reconstruction TIV à faisceau conique. Les personnes du métier comprendront de quelle manière on peut concevoir le code nécessaire pour exécuter les tâches représentées par le bloc 27. Il est à noter également qu'un algorithme existant de reconstruction TIV à faisceau conique peut être modifié pour exécuter les tâches représen- tées par le bloc 27, en plus de l'exécution de la tâche normale de re- construction TIV à faisceau conique. Ceci est accompli en modifiant les parties de filtrage et de rétroprojection de l'algorithme de recon- struction TIV à faisceau conique pour prendre en compte la non-uni- formité de l'espacement des détecteurs élémentaires dans le plan 35 de détecteur vertical virtuel. Les personnes du métier comprendront de quelle manière on peut modifier un algorithme existant de reconstruc- tion TIV à faisceau conique de cette façon pour exécuter toutes ces fonctions. Les personnes du métier comprendront que l'invention n'est pas limitée aux modes de réalisation évoqués ci-dessus en référence aux figures 3 à 5. Les personnes du métier comprendront que des mo- difications peuvent être apportées aux modes de réalisation évoqués cidessus qui sont dans la portée de l'invention. Il est bien entendu, également, que bien que la présente invention ait été décrite en réfé- rence à l'analyse d'un patient, l'invention peut aussi être appliquée à des utilisations industrielles, comme, par exemple, l'analyse de cartes à circuits imprimés. REVENDICATIONS 1. Appareil pour générer une image tridimensionnelle (3D) à partir d'une pluralité d'images de projection à deux dimensions (2D), caractérisé en ce qu'il comprend: un système de tomosynthèse circulaire, comprenant une sour- ce (19) de rayons X et un détecteur (15), la source de rayons X projetant des rayons X vers un objet à analyser (18), le détecteur détectant des rayons X émis par la source de rayons X, le détecteur générant des signaux électriques en réponse aux rayons X qui le frappent, les si- gnaux électriques correspondant aux données d'image acquises par le système de tomosynthèse circulaire; et un dispositif logique (10) configuré pour recevoir les signaux électriques générés par le détecteur, ledit dispositif logique (10) étant configuré pour transformer les données d'image de tomosynthèse cir- culaire en données d'image de tomographie informatisée volumétrique à faisceau conique (TIV à faisceau conique). 2. Appareil selon la revendication 1, caractérisé en ce que le- dit dispositif logique (10) est un ordinateur, ledit ordinateur (10) étant configuré pour exécuter un algorithme de transformation (13), l'algo- rithme de transformation transformant les données d'image de tomo- synthèse circulaire en données d'image TIV à faisceau conique, ledit ordinateur étant en outre configuré pour exécuter un algorithme (14) de reconstruction TIV à faisceau conique qui opère sur les données d'image transformées pour générer l'image 3D. 3. Appareil selon la revendication 2, caractérisé en ce que le détecteur (15) du système de tomosynthèse circulaire est disposé dans un plan qui est parallèle à un plan de trajectoire de la source de rayons X du système de tomosynthèse circulaire, dans lequel l'algorithme de transformation (13) transforme les données d'image recueillies par le détecteur du système de tomosynthèse circulaire en un plan (35) de dé- tecteur vertical virtuel qui est orthogonal au plan dans lequel se trouve le détecteur (15) de tomosynthèse, le plan de détecteur vertical virtuel correspondant à un plan dans lequel un détecteur de système TIV à faisceau conique serait disposé par rapport à une source de rayons X de système TIV à faisceau conique. 4. Appareil selon la revendication 3, caractérisé en ce que, lorsque les données d'image sont transformées en plan (35) de détec- teur vertical virtuel, une kème rangée du détecteur (15) du système de tomosynthèse circulaire est transformée en une rangée du plan (35) de détecteur vertical virtuel à une hauteur Zvk le long d'un axe des z, Zvk étant définie par l'équation: Zv,k = ztD/(D - kAt), dans laquelle zt correspond à la hauteur le long de l'axe des z du plan dans lequel le détecteur (15) du système de tomosynthèse circulaire est placé, dans laquelle D est la distance séparant l'emplacement de la source (19) de rayons X du système de tomosynthèse circulaire du plan (35) de détecteur vertical virtuel suivant l'axe des x, l'axe des x étant orthogonal à l'axe des z, et dans laquelle At correspond à un espace- ment des détecteurs élémentaires du détecteur (15) du système de to- mosynthèse circulaire, et o l'espacement vertical suivant l'axe des z, Azk, entre la knme et la (k-l)ème rangées des données d'image du plan (35) de détecteur vertical virtuel est défini par Az,k = Zv,k - Zv,k-1 = (Zvk/(D-(k- 1)At)) At et o l'espacement des détecteurs élémentaires, Ayk des données d'image dans le plan (35) de détecteur vertical virtuel suivant l'axe des y est défini par l'équation: Ay,k = (Zv,k/Zt) a 5. Appareil selon la revendication 4, caractérisé en ce que, une fois que les données d'image ont été transformées en données d'image dans le plan (35) de détecteur vertical virtuel, les données d'image transformées sont traitées par l'ordinateur (10) pour corriger toute non-uniformité de l'espacement des détecteurs élémentaires dans les directions suivant les axes y et z des détecteurs élémentaires du plan de détecteur vertical virtuel. 6. Appareil selon la revendication 5, caractérisé en ce que l'ordinateur (10) traite les données d'image transformées pour corriger toute nonuniformité de l'espacement des détecteurs élémentaires en exécutant un algorithme de correction qui réalise les étapes consistant àh: tracer une ligne à partir de chaque détecteur élémentaire du plan (35) de détecteur vertical virtuel correspondant aux données d'image transformées, jusqu'à la source (19) de rayons X du système de tomosynthèse circulaire en passant à travers le détecteur (15) du système de tomosynthèse circulaire; déterminer quel détecteur élémentaire du détecteur (15) du système de tomosynthèse circulaire est coupé par la ligne; et utiliser les valeurs d'intensité du détecteur élémentaire inter- secté et d'une pluralité de détecteurs élémentaires entourant le détec- teur élémentaire intersecté dans le détecteur (15) du système de tomo- synthèse circulaire pour déterminer par interpolation une valeur d'in- tensité à attribuer à chaque détecteur élémentaire correspondant du plan (35) de détecteur vertical virtuel. 7. Appareil selon la revendication 6, caractérisé en ce que l'algorithme de correction détermine par interpolation une valeur d'in- tensité en utilisant les valeurs d'intensité du détecteur élémentaire in- tersecté et de la pluralité de détecteurs élémentaires entourant le dé- tecteur élémentaire intersecté pour obtenir une valeur d'intensité moyenne qui est attribuée au détecteur élémentaire correspondant dans le plan (35) de détecteur vertical virtuel. 8. Appareil selon la revendication 7, caractérisé en ce que l'algorithme de correction est constitué de l'algorithme de reconstruc- tion TIV à faisceau conique. 9. Appareil selon la revendication 7, caractérisé en ce que l'algorithme de reconstruction TIV à faisceau conique est un algorith- me de Feldkamp. 10. Appareil selon la revendication 8, caractérisé en ce que l'algorithme de correction est constitué de l'algorithme de reconstruc- tion TIV à faisceau conique, l'algorithme de reconstruction TIV à fais- ceau conique étant un algorithme de Feldkamp modifié, l'algorithme de Feldkamp modifié comprenant des segments de code de filtrage et de rétroprojection qui ont été modifiés pour réaliser les fonctions de l'algorithme de correction. 1 1. Procédé pour générer une image 3D à partir d'une plurali- té d'images de projection 2D, caractérisé en ce qu'il comprend les éta- pes consistant à: utiliser un système de tomosynthèse circulaire pour acquérir les images de projection 2D, le système comprenant une source (19) de rayons X et un détecteur (15), la source de rayons X projetant des rayons X à travers un objet à analyser (18), le détecteur détectant les rayons X qui passent à travers l'objet, le détecteur générant des si- gnaux électriques en réponse aux rayons X qui le frappent, les signaux électriques correspondant aux données d'image acquises par le système de tomosynthèse circulaire; et transformer les données d'image de tomosynthèse circulaire en données d'image de tomographie informatisée volumétrique à fais- ceau conique (TIV à faisceau conique); et exécuter une reconstruction de TIV à faisceau conique sur les données d'image transformées pour générer l'image 3D. 12. Procédé selon la revendication 11, caractérisé en ce que le détecteur (15) du système de tomosynthèse circulaire est disposé dans un plan qui est parallèle à un plan de trajectoire de la source (19) de rayons X, dans lequel l'étape de transformation transforme les don- nées d'image recueillies dans le plan dans lequel se trouve le détecteur du système de tomosynthèse circulaire en un plan de détecteur vertical virtuel qui est orthogonal au plan de trajectoire de la source de rayons X, le plan de détecteur vertical virtuel correspondant à un plan dans lequel un détecteur de système TIV à faisceau conique serait disposé par rapport à une source de rayons X de système TIV à faisceau coni- que. 13. Procédé selon la revendication 12, caractérisé en ce que, au cours de l'étape de transformation, une kème rangée du détecteur du système de tomosynthèse circulaire est transformée en une rangée du plan (35) de détecteur vertical virtuel à une hauteur Zvk le long d'un axe des z, Zvk étant définie par l'équation Zvk = ztD/(D - kt), dans laquelle zt correspond à la hauteur le long de l'axe des z du plan dans lequel le détecteur du système de tomosynthèse circulaire est pla- cé, dans laquelle D est la distance séparant l'emplacement de la source (19) de rayons X du système de tomosynthèse circulaire du plan (35) de détecteur vertical virtuel suivant l'axe des x, l'axe des x étant or- thogonal à l'axe des z, et dans laquelle At correspond à un espacement des détecteurs élémentaires du détecteur du système de tomosynthèse circulaire, et o l'espacement vertical suivant l'axe des z, Azk, entre la kème et la (k-l)ème rangées des données d'image du plan de détecteur vertical virtuel est défini par Az,k = Zv,k - Zv,k-1 = (zvk/(D-(k-1)At)) At et o l'espacement des détecteurs élémentaires, AYk, des données d'image dans le plan (35) de détecteur vertical virtuel suivant l'axe des y est défini par l'équation: AYk = (Zv,k/Zt) ' At- 14. Procédé selon la revendication 13, caractérisé en ce qu'après l'étape de transformation des données d'image et avant l'étape d'exécution d'une reconstruction TIV à faisceau conique, on exécute une étape de correction sur les données d'image transformées, l'étape de correction traitant les données d'image transformées pour corriger toute non-uniformité de l'espacement des détecteurs élémentaires dans les directions suivant les axes y et z des détecteurs élémentaires du plan (35) de détecteur vertical virtuel. 15. Procédé selon la revendication 14, caractérisé en ce que l'étape de correction corrige toute non-uniformité de l'espacement des détecteurs élémentaires en exécutant les étapes consistant à: tracer une ligne à partir de chaque détecteur élémentaire du plan (35) de détecteur vertical virtuel correspondant aux données d'image transformées, jusqu'à la source (19) de rayons X du système de tomosynthèse circulaire en passant à travers le détecteur (15) du système de tomosynthèse circulaire; déterminer quel détecteur élémentaire du détecteur (15) du système de tomosynthèse circulaire est coupé par la ligne; et utiliser les valeurs d'intensité du détecteur élémentaire inter- secté et d'une pluralité de détecteurs élémentaires entourant le détec- teur élémentaire intersecté dans le détecteur (15) du système de tomo- synthèse circulaire pour déterminer par interpolation une valeur d'in- tensité à attribuer à chaque détecteur élémentaire correspondant du plan (35) de détecteur vertical virtuel. 16. Procédé selon la revendication 15, caractérisé en ce qu'au cours de l'étape d'utilisation des valeurs d'intensité, la valeur d'intensité à attribuer est déterminée par interpolation en utilisant les valeurs d'intensité du détecteur élémentaire intersecté et de la pluralité de détecteurs élémentaires entourant le détecteur élémentaire intersec- té pour obtenir une valeur d'intensité moyenne à attribuer au détecteur élémentaire correspondant dans le plan (35) de détecteur vertical vir- tuel. 17. Procédé selon la revendication 16, caractérisé en ce que le procédé est exécuté par un ordinateur (10), l'ordinateur étant confi- guré pour exécuter un algorithme de transformation (13), l'algorithme de transformation transformant les données d'image de tomosynthèse circulaire en données d'image de tomographie informatisée volumétri- que à faisceau conique, ledit ordinateur (10) étant configuré pour exé- cuter un algorithme de reconstruction de tomographie informatisée vo- lumétrique (TIV) à faisceau conique qui opère sur les données d'image transformées pour générer l'image 3D. 18. Procédé selon la revendication 17, caractérisé en ce que l'étape de correction est exécutée au cours de l'étape d'exécution d'une reconstruction TIV à faisceau conique. 19. Procédé selon la revendication 17, caractérisé en ce que l'étape de reconstruction TIV à faisceau conique est exécutée par un algorithme de reconstruction TIV à faisceau conique. 20. Procédé selon la revendication 19, caractérisé en ce que l'algorithme de reconstruction TIV à faisceau conique est un algorith- me de Feldkamp. 21. Procédé selon la revendication 19, caractérisé en ce que l'algorithme de reconstruction TIV à faisceau conique est un algorith- me de Feldkamp modifié, l'algorithme de Feldkamp modifié compre- nant des segments de code de filtrage et de rétroprojection qui ont été modifiés pour réaliser les fonctions de l'étape de correction.