Brevet FR2846121A1 — publié le 20040423

PROCEDE DE RECONSTRUCTION D'UNE IMAGE A PARTIR D'UN JEU DE PROJECTIONS PAR APPLICATION D'UNE TRANSFORMEE EN ONDELETTE5 DESCRIPTION Le sujet de cette invention est un procédé de reconstruction d'une image à partir de jeux de projections de cette image et en appliquant une transformation en ondelette. Les procédés de tomographie consistent à examiner un objet inanimé ou un être vivant par un réseau de détecteurs mobiles qui prend une série de vues en tournant autour de lui. Ces vues sont des projections de la propriété permettant d'exprimer 15 l'image (normalement une absorption d'un rayonnement traversant l'objet ou une scintillation d'un corps radioactif ingéré par l'objet), c'est-à-dire des sommes de la propriété le long de lignes traversant l'objet et définies par la collimation des détecteurs. Chaque 20 détecteur mesure une projection de l'image à chaque vue. Quand un nombre suffisant de vues et de projections a été pris, on procède à l'inversion des résultats pour obtenir la valeur de la propriété à chaque point de l'objet; cette inversion est 25 comparable à l'inversion d'un système d'équations de grandes dimensions et peut s'effectuer ouvertement par des méthodes algébriques, ou, plus fréquemment, par des méthodes analytiques par lesquelles on applique des opérations numériques successives aux projections sans 30 inverser directement le système. Un grand nombre de méthodes existent, parmi lesquelles on citera celle qui est décrite dans le brevet français 2 615 619 qui est le premier de cette équipe de recherche, et le brevet français 2 810 141 plus récent qui présente quelques ressemblances avec le procédé qui sera décrit ici. On 5 peut aussi citer un article (donnez la référence complète)... de Grass. Ces opérations aboutissent à ce qu'on appelle la rétroprojection des mesures, c'est-àdire au calcul de la valeur de la propriété prise à chacun des points de la ligne de la projection. Il peut être avantageux de travailler avec des résultats exprimés dans le domaine fréquentiel de Fourier, comme le prouve le deuxième document. Des transformées numériques d'autre nature ont aussi été employées. On se reportera aux figures 6, 7 et 8 pour une description concrète et schématique d'un procédé de prise de mesures tomographiques. Une source F de rayonnement et un système de détection 3 sont mobiles le long d'un cadre annulaire 2 à des positions 20 diamétralement opposées, et un rayonnement conique originaire de la source F atteint le système de détection 3 après avoir traversé l'objet 1 à étudier. La pièce essentielle du système de détection 3 est un réseau 4 bidimensionnel de détecteurs 5. Des 25 projections R de l'image tridimensionnelle de l'objet 1 sont mesurées par ceux des détecteurs 5 qui sont inclus dans un pourtour 15 de l'"ombre" de l'objet 1. Un nombre important de vues de ce genre est pris à autant d'angles 0 différents d'orientation du réseau 4 des 30 détecteurs 5. Il est fréquemment considéré, conventionnellement pour les calculs, un réseau fictif 4' de détecteurs 5' sur un plan de détection Pdet qui est parallèle au réseau 4 réel et passe par le centre O du cadre 2. Des coordonnées p et q sont définies pour identifier les détecteurs 13 d'après leurs lignes et 5 leurs colonnes. Des calculs de réarrangement, courants dans l'art, permettent de transposer les mesures d'un réseau 4 quelconque au réseau fictif 4' et d'appliquer les algorithmes de reconstruction à celui-ci. Le brevet français 2 615 619 sera rappelé 10 ici puisqu'il expliquait en détail quelles opérations numériques, notamment celles qui sont dites de filtrage et de rétroprojection, permettent d'obtenir l'image de l'objet à partir de ses projections; des méthodes différentes existent cependant. Le problème auquel on s'attaque ici est la réduction du temps ou du volume des calculs pour inverser le système des mesures. On sait qu'il s'agit d'une des limitations les plus graves dans les procédés de tomographie, et de 20 nombreux procédés nouveaux ont été conçus pour le même objectif que l'invention, dont celui du deuxième brevet cité. L'idée exposée ici est d'exploiter les propriétés particulières d'une transformée numérique, 25 dite décomposition en ondelette, des projections pour prévoir des emplacements de projections négligeables ou insignifiantes et de ne pas appliquer les calculs d'inversion à ces emplacements. Des perfectionnements offrent encore une plus grande réduction des calculs. Sous sa forme la plus générale, l'invention concerne un procédé de reconstruction d'une image comprenant un série de décompositions successives des jeux de projections en ondelette donnant des imagettes des jeux de projection, comprenant une approximation et des séries successives de détails homologues de chaque 5 jeu, puis de rétroprojections des imagettes des jeux de projection en imagettes de l'image à reconstruire et une combinaison des imagettes à reconstruire par une inversion de décomposition en ondelette donnant ladite image, caractérisé en ce qu'il comprend, dans chacune 10 des séries et successivement pour des imagettes de détails ayant un nombre croissant de points, une recherche de portions insignifiantes, jugées pour dépourvues de contenu, et une recherche de portions homothétiques des portions insignifiantes pour chacune 15 des imagettes de détails qui suivent dans la série; et en ce qu'il comprend une omission des rétroprojections pour toutes les portions insignifiantes et toutes, les portions homothétiques. L'art antérieur comprend un exemple de 20 reconstruction d'une image par une décomposition en ondelette (US - 5 953 388 - A), qui cependant est appliquée là pour reconstruire une portion seulement de l'image, en exploitant la propriété de "localité" de la décomposition, qui n'est guère sensible aux autres 25 portions de l'image, qu'on peut ainsi négliger dans les calculs. L'invention sera maintenant décrite plus concrètement et complètement en liaison aux figures suivantes: - la figure 1 illustre une décomposition d'image en ondelette, - les figures 2, 3 et 4 illustrent certains aspects de l'invention, - la figure 5 illustre un organigramme résumant une réalisation complète du procédé, - la figure 6 illustre un dispositif de principe du procédé, - et les figures 7 et 8 illustrent plus en détail le procédé et certaines rotations employées. Nous commencerons par exposer la 10 transformation d'un signal par une décomposition en ondelette. Plusieurs modèles d'ondelettes existent, qui ont en commun d'être comparables à un filtre passe-bas. Le signal est séparé en deux portions dont l'une, associée aux basses fréquences, peut être tenue pour 15 une approximation du signal, alors que l'autre, associée aux hautes fréquences, exprime plutôt ses détails. Une propriété des ondelettes est que les portions peuvent contenir chacun une moitié des points du signal si bien qu'il n'y a pas de perte 20 d'information par cette décomposition. La décomposition peut être faite dans le domaine direct d'expression du signal ou dans le domaine de Fourier. Dans le cas de projections d'un objet examiné sur un réseau généralement bidimensionnel de 25 détecteurs, les projections peuvent être groupées en jeux bidimensionnels d'après deux de leur coordonnées (en général p et q sur les axes d'un réseau de détecteurs). Toutefois, comme la figure 8 l'illustre, on considère le plus souvent des jeux de projections 30 réarrangées sur un réseau fictif de détecteurs. Dans l'exemple de la figure 8, on regroupe des projections Rx originaires d'un certain nombre de positions Fx successives de la source F du rayonnement, de manière que les projections aboutissant à une même colonne (à q constant) de détecteurs 5' du réseau fictif 4' soient 5 originaires d'une même position Fx, et aussi que les plans de projections soient tous parallèles jusqu'aux positions Fx: on a alors transformé le problème d'une géométrie conique du rayonnement à une géométrie en éventails parallèles dont la résolution et plus simple. 10 De plus, le réseau fictif 4' passe par le centre 0 de rotation et appartient donc au plan de détection Pdet, ce qui facilite aussi les calculs. L'invention pourrait encore être appliquée à des rencontructions de coupes de l'objet 1 à l'aide 15 d'un réseau monodimensionnel de détecteurs (tous à lamême coordonnée p4). On pourrait encore appliquer le principe exposé ci-dessus du réarrangement des projections d'une géométrie en éventail plan à géométrie parallèle. Le traitement des mesures réarrangées se fait en suivant les lignes de détecteurs 5 du réseau fictif 4', successivement pour tous les points. Un signal peut être décomposé successivement en ondelette pour donner plusieurs 25 niveaux de résultats. Les nouvelles décompositions ne concernent que la portion de rondelette qui donnait l'approximation du signal, la portion ou les portions qui donnaient les détails étant préservées. Prenons comme exemple d'objet d'une 30 décomposition en ondelette une image formée de cinq cercles dont un cercle extérieur et quatre cercles de diamètres différents, tous inscrits dans le premier. L'image de départ comprenait n x n points, et si on applique deux fois une décomposition en ondelette de cette image selon le principe précédent, le résultat est donné à la figure 1. La décomposition en ondelettes de l'image donne un ensemble d'imagettes, dont trois sont plus importantes que les autres, comprenant n/2 x n/2 points chacune, et correspondent aux détails horizontaux, aux 10 détails diagonaux et aux détails verticaux de l'image de départ à grande échelle; elles sont notées respectivement Dhl, Ddl et Dvl. Les détails horizontaux de l'image sont obtenus à partir des projections d'angles 0 proche de o ou A, les détails verticaux à 15 partir des projections d'angles 0 proche de n/2 OU 3C/2, et les détails diagonaux à partir des projections d'angles intermédiaires avec les conventions de la.figure 1. Le reste de l'image est composé de quatre imagettes comprenant n/4 x n/4 points chacune et dont 20 trois sont des imagettes de détails horizontaux, diagonaux et verticaux à plus petite échelle, notées Dh2, Dd2 et Dv2, alors que la dernière imagette est une approximation notée AA de l'image de départ. Si on appelle Q et v les fonctions de décomposition d'une 25 image ou d'une imagette en ondelette, la fonction 4 donnant l'approximation et la fonction v les détails, les fonctions à appliquer à l'image de départ pour obtenir la décomposition de la figure 1 sont données par le tableau I. AA:A X (X1), (X2) 2) 2 Dv2: (x1) 4 (x2) () (x2) 2 2 Dd2: (xl) 4 (x2) (X) 2 () 2 2 Dh2: (xl) 4 (x2) ' (2) X (2 2 2 Dvl: (x1) ' (X2) Ddl: ' (x1) ' (x2) Dhl: ' (x1) 4 (x2) L'invention consiste à effectuer la rétroprojection sur chacune des imagettes des jeux de projections décomposées en ondelette et à combiner les 10 imagettes rétroprojetées par une inversion de la disposition en ondelette pour obtenir l'image recherchée. La décomposition en ondelettes est propice à diverses simplifications qui accélèrent beaucoup les calculs. Ces simplifications sont faites entre la 15 décomposition et la combinaison. La première d'entre elles concerne des filiations qu'on peut établir entre des détails homologues à des échelles différentes. On considère pour cela des séries d'imagettes donnant des détails de 20 même nature. La figure 2 (qui illustre une décomposition d'une image ressemblant à celle de la figure 1, mais a trois niveaux de décomposition) illustre, pour les trois imagettes des détails horizontaux Dhl, Dh2 et Dh3, des portions homologues 25 J1, J2 et J3 qui occupent la même position et la même superficie relative sur chacune de ces imagettes et se déduisent ainsi l'une de l'autre par une homothétie géométrique dans leurs imagettes. On peut faire l'hypothèse que, pour la plupart des images rencontrées en pratique (à l'exception notamment des images texturées), si une portion telle que J3 a un contenu insignifiant, 5 c'est-à-dire qui ne révèle rien par rapport aux valeurs significatives de l'imagette, les portions homologues à plus grande échelle, ici J2 et Jl, seront elles aussi insignifiantes. Conformément à l'invention, on commence 10 donc, pour les imagettes des détails, à rétroprojeter les détails à l'échelle la plus petite, puis les détails aux échelles de plus en plus grandes. Un seuil numérique est appliqué aux coefficients de l'ondelette, c'est-à-dire aux valeurs que prend la transformée dans 15 l'imagette considérée. Une valeur inférieure à ce seuil donne une portion insignifiante, telle que J3. Cependant, les portions insignifiantes des imagettes ne sont pas reconstruites, c'est-à-dire que les calculs de rétroprojection ne sont pas effectués pour elles. On procède concrètement en construisant un masque numérique avant de rétroprojeter l'imagette. Sur les détails horizontaux Dh, le masque est construit pour la première fois pour l'imagette Dh3. Il prend une valeur égale à 0 pour les portions insignifiantes 25 telles que J3 et égale à 1 ailleurs. Les coefficients du masque se suivent dans un ordre déterminé, par exemple ligne après ligne. Les calculs de rétroprojection sont appliqués sur l'imagette considérée dans l'ordre des coefficients du masque. 30 Quand un coefficient est égal à 0, aucun calcul n'est effectué pour le point correspondant de l'imagette, auquel on attribue une valeur nulle dans l'imagette retroprojectée. En abordant l'image suivante des détails horizontaux (Dh2), le masque numérique permet de ne pas 5 considérer la portion J2, dont les points ne sont pas traités par l'unité de calcul qui - effectue la rétroprojection. Si par exemple le masque de l'imagette Dh3 possède un coefficient nul à la ligne i et la colonne j, on prévoit que les quatre points des lignes 10 2i et 2i+1, et des colonnes 2j et 2j+l de l'imagette Dh2 auront aussi des valeurs insignifiantes. Les calculs de rétroprojection ne seront pas effectués pour ces points. Un masque numérique est ainsi construit 15 pour chaque niveau de décomposition. Le masque décrivant l'imagette Dh2 comprendra par convention des coefficients égaux à O pour toute portion homologue à une portion aux coefficients nuls (comme J3) du masque de l'imagette correspondante à plus petite échelle; 20 pour déterminer les valeurs 1 ou O des autres points du masque de Dh2, on utilisera encore les comparaisons avec le seuil conventionnel. D'autres portions insignifiantes, ayant des valeurs nulles de masque, pourront apparaître. C'est ce qu'on a illustré sur les 25 imagettes Dvl, Dv2 et Dv3, des décompositions successives des détails verticaux. L'imagette Dv3 comprend une portion insignifiante K3 dont on retrouve les homologues K2 et Kl sur les décompositions à plus grande échelle. L'imagette Dv3 ne comprend, dans cet 30 exemple, aucune autre portion insignifiante, mais il a été possible d'en trouver trois autres, notées L2, M2 il et N2, sur l'imagette suivante Dv2. Au niveau de décomposition à plus grande échelle, celui de l'imagette Dvl, on ne s'occupera pas de rétroprojecter les points situés aux portions Ll, Ml et Ni homologues de L2, M2 et N2. On va maintenant examiner d'autres particularités du procédé de l'invention favorables, à une accélération des calculs. La première repose sur l'égalité entre la 10 transformée de Fourier d'une projection de l'image à un angle fixé et la transformée de Fourier de l'image sur une droite de même angle passant par l'origine. En se reportant à la figure 3, o un jeu de projections à été converti dans le domaine de Fourier 15 pour donner des projections de nature fréquentielle dans le système d'axes notés 4l et 42, le support du jeu de projection dans le domaine de Fourier comprend les valeurs comprises entre -V0 et v0 pour 41 comme pour 42. Une décomposition en ondelette, comme celle de la 20 figure 1, fait apparaître l'approximation encore notée AA dans les fréquences les plus basses, autour de l'origine 41 = 42 = 0, alors que les détails se trouvent de part et d'autre de cette approximation, à des fréquences de plus en plus grandes pour les détails 25 à grande échelle. On a représenté une droite passant par l'origine, en faisant un angle 8 avec l'axe horizontal 4l. Cette droite traverse l'approximation AA, ainsi que les imagettes des détails verticaux Dvl et Dv2. Mais on 30 est assuré que les projections faisant cet angle 0 seront tout à fait inutiles pour les calculs de rétroprojection des détails diagonaux Dd et horizontaux Dh puisque la droite d'angle 0 passe à l'écart de leurs imagettes. L'application de l'invention comprend alors 5 une sélection, pour chacune des catégories d'imagettes, des angles de projections qui seront utilisés dans les calculs. Le calcul sur le support des projections est élémentaire. La figure 4 reprend la division d'une image 10 décomposée en ondelette et transposée dans le domaine de Fourier. Une reconstruction parfaite peut être obtenue avec l'emploi d'un nombre déterminé de projections, dépendant de la fréquence (discrétisation) de reconstruction de l'image de l'objet. Ainsi, pour 15 reconstruire l'approximation AA, on peut démontrer qu'il suffira de choisir un nombre de projections, parmi celles qui ont été faites, correspondant à une fréquence maximale vl servant de rayon à un cercle circonscrivant la représentation fréquentielle de 20 l'approximation AA dans la décomposition. De même, les détails à petite échelle seront complètement rendus en utilisant un nombre de projections correspondant à la fréquence v2 dans le cercle ayant ce rayon, dans le système d'axes 41, 42, et circonscrivant les 25 représentations fréquentielles du groupe de détails Dv2, Dd2 et Dh2 à même échelle. Le support des projections dans le domaine de Fourier permet donc de connaître aisément les fréquences maximales nécessaires pour la rétroprojection parfaite des imagettes 30 respectives. La mise en oeuvre de l'invention pour exploiter ces deux particularités comprendra donc, à la rétroprojection de chacune des imagettes, une sélection des projections utiles pour cette rétroprojection, en 5 écartant les autres; et éventuellement une restriction du nombre de projections utiles effectivement exploitées par le calcul pour n'en garder que le nombre utile. La figure 5 est un organigramme de 10 l'ensemble du procédé décrit ici. REVENDICATIONS 1. Procédé de reconstruction d'une image à partir de jeux de projections de cette image, comprenant une 5 série de décompositions successives des jeux de projections en ondelette donnant des imagettes des jeux des projections, comprenant des images d'une approximation (AA) et des séries successives (Dd, Dh, Dv) de détails homologues de chaque jeu, puis des 10 rétroprojections des imagettes des jeux de projection en imagettes de l'image à reconstruire et une combinaison des imagettes à reconstruire par une inversion de décomposition donnant ladite image, caractérisé en ce qu'il comprend, dans chacune des 15 séries et successivement pour des imagettes de détails ayant un nombre croissant de points, une recherche de portions insignifiantes jugées pour dépourvues de contenu, et une recherche de portions homothétiques des portions insignifiantes dans chacune des imagettes de 20 détails qui suivent dans la série; et en ce qu'il comprend une omission des rétroprojections pour toutes les portions insignifiantes et toutes les portions homothétiques. 2. Procédé de reconstruction d'une image suivant la 25 revendication 1, caractérisé en ce qu'il comprend une sélection de régions d'angles (e) des jeux de projections qui sont utilisés dans les rétroprojections. 3. Procédé de reconstruction d'une image suivant la 30 revendication 2, caractérisé en ce qu'il comprend une sélection d'un nombre de projections qui sont utilisées dans la rétroprojections. 4. Procédé de reconstruction d'une image suivant la revendication 2, caractérisé en ce que la sélection des 5 régions d'angle est réalisée en fonction d'un support d'une transformée de Fourier des jeux de projections décomposés en ondelette. 5. Procédé de reconstruction d'une image suivant la revendication 3, caractérisé en ce que la sélection du 10 nombre de projections est réalisée en fonction de fréquences maximales d'un support d'une transformée de Fourier des jeux de projection décomposés en ondelette.