Brevet FR2863079A1 — publié le 20050603

PROCEDE ET DISPOSITIF DE RECONSTITUTION D'IMAGE SUR DES BARRETTES DE DETECTEURS. DOMAINE TECHNIQUE GENERAL. La présente invention concerne un procédé de reconstitution d'image acquise grâce à des barrettes de détecteurs montées sur un projectile. ETAT DE L'ART. Des projectiles peuvent être autoguidés sur une cible en mouvement lorsque cette dernière dégage de la chaleur par exemple la chaleur de propulseurs ou de moteurs - grâce à un dispositif de détection de points chauds embarqué sur le projectile. Une fois le point chaud repéré, le projectile ajuste sa trajectoire sur le point chaud pour atteindre sa cible. La figure 1 montre schématiquement une coupe longitudinale du plan focal d'un dispositif de détection de points chauds de l'état de la technique. Comme le montre la figure 1, le dispositif de détection de points chauds comporte classiquement quatre barrettes 1 de détection de forme rectangulaire et placées les unes par rapport aux autres pour former une croix grecque régulière centrée sur l'axe longitudinal 2 du dispositif. Les barrettes 1 sont placées dans le plan focal du dispositif de détection. Leur longueur est égale au champ de l'optique, et leur largeur approximativement au diamètre de la tache image (réponse percussionnelle) de l'optique. Les barrettes 1 sont fixes par rapport au projectile symbolisé par le repère ortho normal 3 lui même centré sur l'axe 2. Comme les barrettes 1 sont fixes par rapport au repère 3, pour pouvoir observer la scène, on fait tourner l'axe optique 6 du dispositif autour de l'axe 2 selon la flèche référencée par 5, l'axe optique 6 étant à une distance R de l'axe 2. La trajectoire de l'axe 6 dans le plan focal est représentée par le cercle 4. Lors de cette rotation 5 de l'axe optique 6 autour de l'axe 2, les barrettes 1 détectent un point chaud. La détection et le positionnement de ce point chaud par rapport au projectile permettent au projectile de se diriger vers le point chaud et par conséquent la cible. Les dispositifs de détection et les procédés de traitement de signal permettant le positionnement de points chauds précédents présentent cependant des inconvénients. En effet, la technique antérieure ne permet de traiter qu'un seul point 5 chaud. Ils ne permettent pas de reconstituer une image de la scène autour du projectile. Par conséquent, la fiabilité des projectiles comportant de tels dispositifs et fonctionnant avec de tels procédés de détection et de 10 traitement de signal n'est pas optimale. PRESENTATION DE L'INVENTION L'invention propose de pallier ces inconvénients. Un des buts de l'invention est de proposer un procédé permettant de reconstituer une image grâce aux signaux captés sur des barrettes de 15 détecteurs montés sur un projectile. A cet effet, l'invention propose un procédé de reconstitution d'une image lm acquise grâce à des barrettes de détecteurs montés sur un projectile, caractérisé en ce qu'on reconstitue le vecteur représentatif de l'image lm à partir du vecteur représentatif des signaux B issus des barrettes 1 en résolvant le système linéaire tel que: LÉIm=B où L est la matrice reprenant les signaux issus de chaque barreau pour tous les pixels de l'image. L'invention est avantageusement complétée par les caractéristiques 25 suivantes, prises seules ou en une quelconque de leur combinaison techniquement possible: - le procédé comporte les étapes selon lesquelles: - on numérote les barrettes et on construit le vecteur B en juxtaposant les signaux reçus de chaque barrette (1) en fonction du numéro de 0 barrette-; - - on divise en pixels le plan dans lequel sont situées les barrettes, les pixels étant situés à l'intersection de lignes et de colonnes et on construit le vecteur lm en juxtaposant les lignes du plan ainsi divisé ; et - on construit chaque coefficient L; de la matrice L en reportant en colonne pour chaque pixel courant j l'intensité du signal numérisé issu de chaque barrette i pour un point chaud centré sur le pixel j, chaque ligne de L correspondant à un barreau et chaque colonne à un pixel; - la résolution du système linéaire utilise une méthode du type minimisation d'une fonctionnelle quadratique; - la méthode de résolution est une méthode des moindres carrés et en ce que la fonctionnelle est du type 2 É Im É L É L É lm Im L É B où Im,L désignent les transposées de lm et de L; - la méthode des moindres carrés comporte une étape: d'initialisation selon laquelle on annule l'image sur tout le champ des barrettes; - de définition d'un seuil initial; et - de bouclage consistant à effectuer: une descente par gradient conjugué, - un seuillage de l'image et une mise à zéro des signaux de l'image dont l'intensité est en dessous du seuil; et - un abaissement du seuil; la boucle s'arrêtant au bout d'un certain nombre d'itérations ou lorsque 20 l'image reconstituée est suffisamment résolue pour être exploitée par des moyens de traitement et de guidage du projectile; - la résolution du système linéaire utilise une méthode du type minimisation ou maximisation d'une fonctionnelle sous contrainte; - la méthode de résolution est une maximisation de l'entropie de l'image 25 avec la contrainte L Im = B; - la méthode de résolution utilise une technique de reconstruction algébrique multiplicative MART ; L'invention concerne également un dispositif de reconstitution d'image comportant des barrettes de détection et des moyens de traitement aptes à 30 mettre en oeuvre un procédé selon l'invention. PRESENTATION DES FIGURES D'autres caractéristiques, buts et avantages de l'invention ressortiront de la description qui suit, qui est purement illustrative et non limitative, et qui doit être lue en regard des dessins annexés sur lesquels: - la figure 1, déjà commentée, représente une coupe longitudinale d'un plan focal d'un dispositif de détection selon l'art antérieur; - la figure 2 représente schématiquement une coupe longitudinale d'un plan focal d'un dispositif de détection selon l'invention; - la figure 3 représente schématiquement le trajet d'une barrette devant une 10 scène; - la figure 4 représente schématiquement l'échantillonnage du plan des fréquences le long de directions perpendiculaires aux barrettes; - la figure 5A montre l'influence du nombre des barrettes et de la taille du champ de l'image sur l'échantillonnage des fréquences; - la figure 5B représente schématiquement un mode de réalisation d'un dispositif comportant sept barrettes; - les figures 6A et 6B représentent schématiquement une méthode de construction du vecteur lm, - les figures 7 et 8 représentent schématiquement une méthode de 20 construction de la matrice L; - les figures 9A et 9D représentent schématiquement images reconstituées par un procédé possible, au fur et à mesure des itérations du procédé ; - la figure 10 représente une mire; et - la figure 11 représente une image de la mire de la figure 10 reconstituée 25 par un autre procédé possible après une dizaine d'itérations. DESCRIPTION DETAILLEE. POSITION DU PROBLEME. Dans toute la suite de la description, on ne considère pas un quelconque balayage de la scène par le système d'auto guidage et de reconstitution d'image. Le mouvement relatif du projectile (notamment la rotation sur lui-même) lors du vol n'est donc pas pris en compte. La figure 2, dont les éléments similaires à la figure.-1 _portent des références numériques identiques, représente schématiquement une coupe longitudinale d'un plan focal d'un dispositif mettant en oeuvre un procédé selon l'invention. Ainsi, un dispositif de détection mettant en oeuvre un procédé selon l'invention comporte au moins trois barrettes 1 de détection de forme rectangulaire et placées les unes par rapport aux autres selon au moins trois directions non colinéaires entre elles. Dans toute la description suivante, sauf mention explicite du contraire, les détecteurs sont orientés radialement, vers l'axe 2 du dispositif, et vers les sommets d'un polygone régulier. Sur la figure 2, les barrettes sont ainsi disposées entre elles en forme de Y régulier centré sur l'axe longitudinal 2 du dispositif. Les barrettes 1 sont placées dans le plan focal du dispositif de détection. Leur longueur est égale au champ de l'optique, et leur largeur approximativement au diamètre de la tache image (réponse percussionnelle) de l'optique. Les barrettes 1 sont fixes par rapport au projectile symbolisé par le repère ortho normal 3 lui même centré sur l'axe 2. L'axe optique 6 du dispositif tourne autour de l'axe 2 selon la flèche référencée par 5 grâce à une toupie, l'axe optique 6 étant à une distance R de l'axe 2. La trajectoire de l'axe 6 dans le plan focal est représentée par le cercle 4. Lors de cette rotation 5 de l'axe optique 6 autour de l'axe 2, les barrettes 1 acquièrent un signal. Dans la suite de la présente description, on suppose que la position 25 angulaire de l'axe optique 6 est connue au fur et à mesure de l'acquisition du signal sur les barrettes 1. Les angles 0 de rotation sont donnés par des moyens de mesure reliée à la toupie. Lors du déplacement de la toupie, l'axe optique 6 effectue un mouvement de translation circulaire devant les barrettes 1. On peut bien entendu également dire, en changeant de référentiel, que les barrettes 1 effectuent un mouvement de translation circulaire par rapport à l'axe optique 6. La figure 3 montre une barrette 1 passant devant l'image 30 de la scène si l'on prend comme référentiel les axes (Ox, Oy) centré sur l'image 30 selon une direction 31 de déplacement en translation. L'image 30 de la scène est acquise pendant une étape 32, c'est-à-dire tout au long de la translation selon la direction 31. On obtient donc un signal 33. On note L l'opérateur qui permet à partir d'une image 30 d'obtenir les signaux issus des barrettes 1. On note lm la fonction de x et de y dans le plan focal qui représente 10 l'image de la scène que l'on cherche à obtenir. La fonction lm représente donc l'image 30 sur la figure 3. On note B les signaux obtenus sur les barrettes 1. La fonction B représente donc les signaux 33 de la figure 3. On rappelle que l'on cherche à déterminer lm à partir de B, c'est-à-dire 15 reconstituer une image de la scène à partir d'un signal reçu par des barrettes. L'équation du problème que l'on cherche à résoudre peut être ainsi écrite: (1) L.Im=B. On voit donc que la résolution de l'équation (1) peut s'effectuer au moyen d'un calcul numérique. Le fait d'effectuer une intégration sur toute la longueur de la barrette 1 s'apparente à une technique d'acquisition de mesures utilisée en tomographie. La tomographie est utilisée dans le domaine médical notamment. Cette caractéristique d'acquisition des mesures permet d'utiliser une méthode de reconstruction d'image tomographique. Tout d'abord, l'équation (1) met en valeur l'aspect linéaire du problème et suppose la linéarité des barrettes. De plus, la résolution numérique de ce problème nécessite une discrétisation du plan focal. On a en effet besoin d'un procédé discret de reconstruction tomographique, qui s'appuie par définition et contrairement aux méthodes analytiques, sur une modélisation discrète de l'image à reconstruire et des mesures. Enfin, on verra dans la suite de la présente description que le pouvoir séparateur intrinsèque du traitement du signal lors de la mise en oeuvre 5 d'un procédé est limité par le nombre et la dimension des barrettes. ECHANTILLONNAGE DANS L'ESPACE DES FREQUENCES. D'un point de vue théorique, les signaux issus des barrettes sont directement liés à la transformée de Fourier de l'image analysée. De part la linéarité de l'équation (1), un procédé selon l'invention peut utiliser une étape qui s'apparente à une transformée de Radon. On rappelle en effet qu'on effectue dans un procédé selon l'invention une intégration 32 de l'intensité du signal acquis sur les barrettes, et ce sur toute la longueur de chaque barrette 1 lors de sa translation 31. On obtient alors ensuite une reconstitution d'une image correspondant à un tel signal 33 acquis. Si on applique à ce signal 33 une transformée 34 de Fourier à une dimension notée TF, dans la suite de la présente description, alors on obtient sensiblement un échantillonnage de l'espace des fréquences 300 sur une droite 310. L'espace des fréquences fx et fy représente la transformée de Fourier à deux dimensions de l'image 30 - notée TF2 dans la suite de la présente description. Plus précisément, dans l'hypothèse de barrettes 1 de longueur infinie, on peut écrire, sans tenir compte des caractéristiques des barrettes 1 représentées usuellement par les gains et les décalages ou offsets selon la terminologie de l'homme du métier: (2) TF, (Sk) = TF2 (I)(n (Bk)) x TF, (Det) avec les notations suivantes: I (image) désigne l'éclairement (W/m2) du plan focal; TF, désigne la transformée de Fourier à une dimension; TF2 désigne la transformée de Fourier à deux dimensions; Sk désigne le signal issu du détecteur k en fonction du. déplacement de l'image; n(Bk) est la direction perpendiculaire au grand côté du détecteur k, à savoir la direction 31 dans le cas de la barrette 1 de la figure 3; et Det désigne la fonction de transfert de la forme du détecteur, typiquement une fonction rectangle dont la largeur est la largeur du détecteur. Le paramètre de calcul est le déplacement de l'image suivant n(Bk) . La figure 4 montre que l'image 300 dans l'espace de Fourier d'une image acquise par des barrettes dans l'espace réel est d'autant mieux échantillonnée que le nombre de barrettes dans l'espace réel est important. On rappelle que la direction d'échantillonnage dans l'espace des fréquences est perpendiculaire à la direction d'extension de la barrette dans l'espace réel. On comprend grâce à la figure 4 qu'il est préférable d'avoir un nombre important de barrettes pour avoir une meilleure reconstitution. Dans l'exemple de la figure 4, on a cinq directions d'échantillonnage 41, 42, 43, 44 et 45. Cela signifie que cinq barrettes sont placées dans le plan focal. La figure 4 montre également que la largeur D des barrettes intervient comme une convolution de l'image dans l'espace des fréquences par un disque 50 dont le diamètre est YD. L'influence de la largeur des barrettes 1 sur le nombre de barrettes à disposer dans le plan focal est détaillée dans la suite de la présente description. Le dispositif de mise en oeuvre du procédé selon l'invention comporte 25 en général préférentiellement sept barrettes. PRINCIPE DE RESOLUTION. On rappelle que pour pouvoir résoudre l'équation (1), il faut effectuer une discrétisation du plan focal. La figure 6A montre ainsi schématiquement que l'on divise l'image 60 30 acquise par les barrettes 1 en pixels référencés par 600. On définit ainsi des lignes, que l'on référence par 61, 62, 63 et 64 sur la figure 6A_. On effectue de plus une numérotation de chaque pixel de l'image 60. Chaque pixel 600 de l'image 60 porte ainsi un numéro. La figure 6B montre que l'on représente alors l'image 600 comme un vecteur 65 en mettant les unes après les autres les lignes 61 à 64 par 5 exemple. De même, on construit un vecteur B représentatif de B de l'équation (1) en juxtaposant les signaux reçus de chaque barrette 1 après avoir numéroté les barrettes. La juxtaposition des signaux dans le vecteur B est effectuée en fonction des numéros des barrettes. La figure 7 montre que l'opérateur L de l'équation (1) est représenté par sa matrice. La figure 7 montre que chaque ligne de la matrice L correspond à une barrette courante i du plan focal. Chaque colonne de la matrice L correspond à un pixel j. On reporte alors, en colonne et pour chaque pixel courant jdu plan focal, l'intensité du signal numérisé issu de chaque barrette i pour un point chaud centré sur le pixel j lors d'un tour de toupie. L'évolution de cette intensité sur une barrette i lors d'un tour de toupie est schématisée à la figure 8. On obtient alors la matrice de L. On rappelle que l'équation (1) s'écrit: (1) LÉIm=B On comprend alors qu'avec la connaissance de la matrice de L et du vecteur de B, le problème de reconstitution d'image est ramené à la 25 résolution d'un système linéaire représentatif de l'équation (1). D'une manière générale, le système de l'équation (1) n'est pas inversible et admet une infinité de solutions, du fait que la matrice L est rectangulaire. Plusieurs méthodes de résolution sont possibles pour faciliter la 30 résolution de cette équation. Les performances d'un dispositif mettant en oeuvre un procédé selon l'invention seront alors conditionnées, outre les caractéristiques physiques du dispositif, par les étapes de résolution utilisées pour le système linéaire, ainsi que par la pertinence du paramétrage du problème pour la résolution. Dans toutes les simulations numériques qui vont suivre, les dimensions des barrettes sont (en pm) 700 x 100 et le plan représenté est un carré de côté 750 pm. Ces dimensions sont données à titre indicatif non limitatif et dépendent des dimensions du champ du système optique. PREMIER EXEMPLE DE PROCEDE. Un premier mode de réalisation possible d'un procédé selon l'invention consiste à minimiser une fonctionnelle quadratique au lieu de résoudre 10 directement l'équation (1). De préférence, on minimise l'erreur quadratique entre les mesures et l'image attendue. Cette méthode est appelée méthode des moindres carrés par l'homme du métier. Plusieurs fonctionnelles différentes peuvent être utilisées dans une 15 méthode des moindres carrés. De préférence, on cherche à minimiser la fonctionnelle définie par: (3) !ÉIm LÉLÉIm ImÉLÉB où Imt,L désignent les transposées de lm et de L. Ainsi, selon le procédé possible, on définit lors d'une étape 20 d'initialisation que l'image de départ est nulle dans tout le champ des détecteurs. On définit un seuil initial. Le seuil dépend des applications et des barrettes d'acquisition. Dans notre exemple, il peut être fixé à'/2. On reprend ensuite en boucle la séquence suivante: - descente par gradient conjugué connu de l'homme du métier, seuillage de l'image et mise à zéro des signaux de l'image dont l'intensité est en dessous du seuil abaissement du seuil. La boucle s'arrête au bout d'un certain nombre d'itérations ou lorsque l'image reconstituée est suffisamment résolue pour être exploitée par les moyens de traitement et de guidage du projectile. Bien entendu, d'autres étapes de descente de gradients peuvent être utilisées dans un procédé selon l'invention. Il peut notamment s'agir d'une méthode de descente de gradients simple ou à pas optimal par exemple. Le fait d'effectuer un seuillage de l'image permet d'éviter la formation de signaux parasites sur l'image et d'assurer une meilleure convergence du calcul. Le fait d'abaisser le seuil à chaque itération du calcul permet de maîtriser la limite des itérations. On peut ainsi définir le seuil s par (4) s= 1/2n où n est le numéro de l'itération. Les figures 9A à 9D montrent les résultats des quatre premières itérations d'un procédé utilisant la méthode des moindres carrés. Dans cet exemple, l'image est divisée en 1600 pixels. Le plan focal comporte deux points chauds 91 et 92 de même intensité et visibles notamment sur la figure 9D. Le calcul prend en compte le traitement du signal, notamment l'élimination d'une composante basse fréquence. La figure 9A est l'image reconstituée après une itération avec un seuil égal à 1/2. La figure 9B est l'image reconstituée après deux itérations avec un seuil égal à'/. La figure 9C est l'image reconstituée après trois itérations avec un seuil égal à 1/6. La figure 9D est l'image reconstituée après quatre itérations avec un 25 seuil égal à 1/8. On voit ainsi que le procédé de calcul permet une convergence rapide vers une image reconstituée de la scène. DEUXIEME EXEMPLE DE PROCEDE. Un deuxième mode de réalisation possible d'un procédé selon l'invention consiste à utiliser une méthode de résolution du type maximisation ou de minimisation d'une fonctionnelle sous contrainte. Cette contrainte permet de choisir une solution acceptable parmi une infinité de solutions possibles à l'équation (1). Il peut en effet y avoir une infinité de solutions lors d'une sousdétermination dans les mesures. Une telle sous-détermination peut se produire lorsque les mesures sont d'un nombre limité, à savoir lors d'une acquisition rapide de la scène par les barrettes, ou lorsque les mesures sous certaines incidences ne sont pas accessibles ou exploitables. Un procédé de maximisation ou de minimisation d'une fonctionnelle sous contrainte peut bien entendu être utilisé dans tous les cas d'acquisition. Dans un procédé de résolution sous contrainte, on utilise de préférence une méthode de résolution du type Technique de reconstruction algébrique ART (Algebraic Reconstruction Technique) ou Technique de reconstruction algébrique multiplicative MART (Multiplicative Algebraic Reconstruction Technique) selon la terminologie généralement utilisée par l'homme du métier. De préférence, on utilise une méthode de type MART. Cette méthode de résolution peut s'interpréter comme une minimisation ou une maximisation sous contrainte avec un critère entropique. Un procédé selon l'invention fait en effet converger la solution vers un maximum d'entropie de l'image avec la contrainte donnée par l'équation (1) . Dans une méthode de type MART, on corrige à chaque itération un ensemble de mesures pour obtenir rapidement une solution. Cependant, cet ensemble de mesures est uniquement situé au niveau d'une barrette donnée. On ne corrige donc pas à chaque itération toutes les mesures sur l'image définie sur la figure 6A, ce qui assure la rapidité de la convergence. L'ensemble de mesures à corriger est ensuite déplacé à chaque itération en fonction du déplacement de la barrette considérée par rapport à 30 la scène. On renouvelle bien entendu les étapes du procédé précédemment décrit pour chacune des barrettes. Les principes de la méthode MART sont décrits ci-après. Pour chaque pixel k compris dans la barrette considérée, la correction de l'image lors de l'itération n+1 du procédé s'effectue sensiblement par une multiplication de l'image jn obtenue au pixel k après correction lors d'une itération n par le rapport du signal Bn reçu par une barrette donnée en ce pixel k lors de l'itération n d'une part au signal Rn reconstitué sur l'ensemble de la barrette donnée lors de l'itération n d'autre part. On a donc une expression de la forme: t/(k) _ (k) (5)1 n+ 1 j, nBn Rn où jn+ est l'image au pixel k à l'itération n+1; 10 j(k) est l'image au pixel k à l'itération n n L,,, est un opérateur valant 0 dans les calculs faisant intervenir des coefficients de la matrice L ne concernant pas un pixel compris dans la barrette considérée, et 1 dans les calculs faisant intervenir des coefficients de la matrice L concernant un pixel compris dans la barrette considérée; Ânk est un coefficient permettant de régler la convergence des itérations en répartissant la correction sur l'ensemble des pixels de la barrette considérée; et R. est défini par une relation du type: 6) R = I Lnj I,, l dans laquelle l'indice j de la sommation indique que l'on somme les Ln] sur l'ensemble des pixels j de la barrette considérée. Rn représente donc le signal que donne l'image reconstituée sur la barrette considérée. Il est entendu que lorsque n désigne un indice de matrice, il décrit les 25 indices de chaque matrice cycliquement, modulo la taille de la matrice. L'indice n n'est donc pas limité par la taille des matrices dans les itérations. D'après (6), l'équation (5) peut également s'écrire: (7) n+ +1 = j(k) . Bn Lnj'InJ) J /Znk permet d'effectuer un moyennage de la correction sur l'ensemble des pixels de la barrette. On assure ainsi une conservation de l'énergie lors 5 des différentes itérations et on permet la convergence du procédé itératif. On peut définir /lnk par l'expression: Les itérations sont arrêtées lorsqu'un nombre d'itérations est atteint ou 10 lorsque l'image reconstituée est suffisamment résolue pour être exploitée par les moyens de traitement et de guidage du projectile. La figure 10 montre ainsi l'image d'une mire comportant des motifs remarquables 100. La figure 11 montre l'image reconstituée de la mire de la figure 10 15 après une dizaine d'itérations du procédé. On y distingue clairement les images 110 reconstituées des motifs 100. D'autres méthodes de résolutions du système de l'équation (1) sont bien entendu possibles. NOMBRE DE BARRETTES DANS LE PLAN FOCAL. On rappelle que du fait de la taille limitée du champ de l'image, et si l'on appelle D le diamètre du champ circulaire dans l'espace réel, les images subissent une convolution par un disque 50 dont le diamètre est en 1/D dans l'espace des fréquences. Prenons l'exemple d'une sinusoïde dans l'espace réel, cette sinusoïde étant vue par le disque de diamètre D. Alors l'image de cette sinusoïde tronquée est un disque de diamètre 11D dans l'espace des fréquences. Or, la figure 4 montre clairement que les fréquences élevées ont moins de chances d'être perçues par l'échantillonnage donné par les droites 41 à 45 par exemple, puisqu'elles s'éloignent du centre du repère (Ofx, Ofy). Ainsi, plus les fréquences sont élevées, moins le disque 50 de convolution à de chances d'être perçu par l'échantillonnage. En référence à la figure 5A, le disque 50 est détecté, dans le pire des 10 cas, lorsque il est tangent à deux droites en même temps. En première approximation, on a la relation: (8) tan(-) 2D 1 f 1 ax où O est l'angle entre deux barrettes successives; et f max est la fréquence max que l'on détecte. On a, en approximant la tangente à son angle: (9) 0= 1 On souhaite avoir: (10)fmaxD pour pouvoir observer plus d'une période de la sinusoïde dans l'espace réel. 1 doit donc être inférieur à 1. f max.D Or, comme on a par définition: (11) 0 = 2z, où n est le nombre de barrettes. n D'après (10) et (11), 2 doit être inférieur à 1. n Par conséquent, n est préférentiellement supérieur ou égal à 7. Le dispositif de mise en oeuvre d'un procédé selon l'invention comporte ainsi préférentiellement sept barrettes 1 également réparties, comme le montre la figure 5B. f max D 16 REVENDICATIONS. 1. Procédé de reconstitution d'une image (lm) acquise grâce à des barrettes (1) de détecteurs montés sur un projectile, caractérisé en ce qu'on reconstitue le vecteur représentatif de l'image (lm) à partir du vecteur représentatif des signaux (B) issus des barrettes (1) en résolvant le système linéaire tel que: LÉIm=B où L est la matrice reprenant les signaux issus de chaque barreau pour tous 10 les pixels de l'image. 2. Procédé selon la revendication 1, caractérisé en ce qu'il comporte les étapes selon lesquelles: - on numérote les barrettes et on construit le vecteur B en juxtaposant 15 les signaux reçus de chaque barrette (1) en fonction du numéro de barrette; - on divise en pixels le plan dans lequel sont situées les barrettes, les pixels étant situés à l'intersection de lignes et de colonnes et on construit le vecteur lm en juxtaposant les lignes du plan ainsi divisé ; et - on construit chaque coefficient L,1 de la matrice L en reportant en colonne pour chaque pixel courant j l'intensité du signal numérisé issu de chaque barrette i pour un point chaud centré sur le pixel j, chaque ligne de L correspondant à un barreau et chaque colonne à un pixel. 3. Procédé selon l'une des revendications 1 ou 2, caractérisé en ce que la résolution du système linéaire utilise une méthode du type minimisation d'une fonctionnelle quadratique. 4. Procédé selon la revendication 3, caractérisé en ce que la méthode 30 de résolution est une méthode des moindres carrés et en ce que la fonctionnelle est du type I ÉImÉL ÉLÉIm ImÉL.B où Im,L désignent les transposées de lm et de L. 5. Procédé selon la revendication 4, caractérisé en ce que la méthode des moindres carrés comporte une étape: d'initialisation selon laquelle on annule l'image sur tout le champ des barrettes; - de définition d'un seuil initial; et - de bouclage consistant à effectuer: - une descente par gradient conjugué, un seuillage de l'image et une mise à zéro des signaux de l'image dont l'intensité est en dessous du seuil; et - un abaissement du seuil; la boucle s'arrêtant au bout d'un certain nombre d'itérations ou lorsque l'image reconstituée est suffisamment résolue pour être exploitée par des 15 moyens de traitement et de guidage du projectile. 6. Procédé selon l'une des revendications 1 ou 2, caractérisé en ce que la résolution du système linéaire utilise une méthode du type minimisation ou maximisation d'une fonctionnelle sous contrainte. 7. Procédé selon la revendication 6, caractérisé en ce que la méthode de résolution est une maximisation de l'entropie de l'image avec la contrainte L É Im = B. 8. Procédé selon la revendication 7, caractérisé en ce que la méthode de résolution utilise une technique de reconstruction algébrique multiplicative MART . 9. Dispositif de reconstitution d'image, caractérisé en ce qu'il comporte 30 des barrettes de détection et des moyens de traitement aptes à mettre en oeuvre le procédé selon l'une des revendications là 8. 10. Dispositif selon la revendication 9, caractérisé en ce qu'il comporte sept barrettes disposées radialement par rapport à l'axe du dispositif et non colinéaires entre elles.