FR 2924513 A1 20090605 FR 0708458 A 20071204 PROCEDE DE TRAITEMENT D'UNE IMAGE NUMERIQUE DEFOCALISEE ET DISPOSITIF CORRESPONDANT. L'invention porte sur le traitement des images défocalisées, en particulier le traitement des images radar dans le cadre de l'imagerie produite selon la technique dite SAR pour � Synthetic Array Radar � selon la dénomination anglo-saxonne couramment employée par l'homme du métier. Classiquement, une antenne embarquée à bord d'un porteur émet un signal vers une zone choisie, et récupère les échos émis par des objets localisés dans cette zone. C'est l'analyse de ces échos qui permet alors de reconstituer une image de la zone visée. La technique de l'imagerie SAR permet de synthétiser une antenne de n'importe quelle longueur de façon à obtenir des images radar de haute résolution. Cette technique consiste à scanner le terrain selon une direction choisie en émettant une série d'impulsions tout en dirigeant le faisceau du radar vers la zone dont on souhaite avoir l'image. L'émission de signaux dure un temps d'éclairement Te. Cependant, en raison de l'effet Doppler, les fréquences des échos renvoyés par la zone éclairée sont décalées en fonction du taux de réflexion des objets contenus dans la zone visée. L'effet Doppler est défini par le décalage existant entre la fréquence de l'onde émise et celle de l'onde reçue lorsque l'émetteur et le récepteur sont en mouvement l'un par rapport à l'autre. Ce décalage fréquentiel ainsi que l'imprécision des capteurs inertiels, c'est-à-dire les capteurs embarqués à bord du porteur, entraîne une perte de netteté (ou défocalisation) au sein de l'image. Cette perte de netteté se traduit par l'apparition d'une erreur de phase au sein de la phase de l'image radar. Des méthodes telle que la méthode dite � Map Drift � (également connue sous le nom de � Multilook registration �) décrite dans le manuel � Air and Spaceborne Radar Systems, An Introduction, de Philippe Lacomme, Jean-Philippe Hardange, Jean-Claude Marchais et Eric Normant, Ed. William Andrew publishing, Scitech Publishing. INC., p. 297-301) ou encore la méthode de gradient de phase dite PGA (ou encore sa variante dite PPP pour � Prominent Point Processing �), présentées aux pages 303-311 du manuel précité permettent de mesurer cette erreur de phase afin de traiter le flou de l'image. Cependant ces méthodes sont incompatibles avec des images obtenues avec des temps d'exposition Te très longs (par exemple, au moins supérieur à 10s). Pour estimer l'erreur de phase de l'image défocalisée, la méthode appelée � Map Drift � propose d'établir un modèle polynomiale a priori de cette erreur de phase, de décomposer l'image défocalisée en K imagettes, K correspondant à l'ordre du modèle polynomiale, puis de calculer l'erreur de phase de l'image défocalisée en traitant chaque imagette formée. Cependant, il n'est pas possible de modéliser l'erreur de phase d'une image numérique de haute résolution à l'aide d'un polynôme. La méthode dite � Map Drift � est donc limitée uniquement à des images de faible résolution, obtenue à l'aide d'un temps d'éclairement inférieur à 10s. L'invention vise notamment à apporter une solution à ce problème. Un but de l'invention est de proposer un procédé (et un dispositif correspondant) apte à traiter la défocalisation d'images numériques quel que soit le temps d'éclairement utilisé pour former l'image, et en particulier pour un temps d'éclairement considéré comme long, voire très long. A cet effet, l'invention propose un procédé de traitement d'une image numérique défocalisée, obtenue par l'éclairement d'une zone pendant un temps donné, comprenant une détermination de l'erreur de phase de ladite image numérique défocalisée, ladite détermination comprenant : - une formation à partir de ladite image numérique défocalisée, de N imagettes, N étant un entier, chaque imagette ayant une résolution N fois inférieure à celle de l'image numérique une fois traitée, et - un calcul de ladite erreur de phase de l'image numérique défocalisée. Selon une caractéristique générale de cet aspect, le nombre N d'imagettes formées est fonction de la résolution de ladite image numérique défocalisée. En d'autres termes, on détermine l'erreur de phase de l'image numérique défocalisée à partir d'imagette de résolution inférieure à celle de l'image à traiter. Ceci a notamment pour avantage de permettre la détermination de l'erreur de phase de toute image numérique défocalisée, quelque soit sa résolution et la loi régissant la variation de son erreur de phase. Selon un mode de mise en oeuvre, la détermination de ladite erreur de phase comprend pour chaque imagette, une localisation puis une extraction d'au moins un point remarquable dit ponctuel fort, une estimation pour au moins un desdits ponctuels forts extraits, de son erreur de phase partielle, ledit calcul de l'erreur de phase de l'image numérique défocalisée étant effectué à partir de l'ensemble des erreurs de phase partielles. Chaque erreur de phase partielle peut comprendre un premier terme traduisant la formation de l'imagette considérée et un deuxième terme traduisant la défocalisation de l'image numérique. De préférence, pour chaque imagette : - on localise et extrait au moins deux ponctuels forts, - on estime son erreur de phase partielle pour chaque ponctuel fort localisé et extrait, - on calcule une erreur de phase partielle moyenne en effectuant une moyenne des erreurs de phase partielles estimées pour l'ensemble des ponctuels forts, et - on calcule l'erreur de phase de l'image numérique défocalisée à partir des erreurs de phase partielles moyennes. Le procédé peut comprendre en outre une correction de l'image numérique défocalisée en multipliant le signal temporel représentatif de ladite image numérique défocalisée par le conjugué du signal représentatif de l'erreur de phase déterminée de ladite image numérique défocalisée. Le premier terme de chaque erreur de phase partielle peut comprendre un coefficient de phase traduisant le décalage entre une imagette considérée et une imagette de référence, et une constante de phase. Ledit coefficient de phase peut être déterminé par une intercorrélation de l'imagette considérée et de l'imagette de référence. Ladite constante de phase peut être déterminée en considérant la continuité temporelle de l'erreur de phase partielle entre deux imagettes successives. Selon un mode de mise en oeuvre, on peut itérer l'étape de détermination de l'erreur de phase de l'image numérique défocalisée au plus N fois, en décrémentant le nombre d'imagettes formées à chaque répétition, chaque nouvelle itération étant effectuée sur l'image numérique défocalisée corrigée avec l'erreur de phase déterminée lors de la précédente itération. En variante, chaque imagette formée peut partiellement recouvrir au moins une imagette adjacente. Selon un autre aspect, il est proposé un dispositif de traitement d'une image numérique défocalisée obtenue à l'aide de moyens d'éclairement d'une zone pendant un temps donné, comprenant un moyen de détermination de l'erreur de phase de l'image numérique défocalisée Le moyen de détermination comporte : - un moyen de formation apte à former à partir de ladite image numérique défocalisée, N imagettes, N étant un entier, chaque imagette ayant une résolution N fois inférieure à celle de l'image numérique défocalisée, - un moyen de calcul couplé audit moyen de formation et apte à déterminer l'erreur de phase de l'image numérique défocalisée, à partir des N imagettes. Selon une caractéristique générale de cet aspect, le nombre N d'imagettes formées par le moyen de formation est fonction de la résolution de ladite image numérique défocalisée. Selon un mode de réalisation, lequel le moyen de calcul peut comprendre : -un sous-moyen de détection de ponctuel apte à localiser et à extraire au moins un point remarquable dit ponctuel fort au sein des imagettes, - un sous-moyen d'estimation apte à estimer pour chaque imagette, l'erreur de phase partielle de chaque ponctuel fort extrait, et - un sous-moyen de calcul apte à calculer l'erreur de phase de l'image numérique défocalisée à partir de l'ensemble des erreurs de phase partielles estimées. Le sous-moyen d'estimation peut comprendre un premier module apte à estimer un premier terme traduisant la formation de l'imagette considérée et un deuxième module apte à estimer un deuxième terme traduisant la défocalisation de l'image numérique. Le deuxième module peut être apte à estimer le deuxième terme en effectuant une intercorrélation entre une imagette et une imagette de référence. Le sous-moyen de détection peut être apte à localiser et à extraire au moins deux ponctuels forts. Le sous-moyen de détermination peut alors comprendre en outre un module de moyennage apte à calculer pour chaque imagette, une erreur de phase partielle moyenne à partir des erreurs de phase partielles estimées pour l'ensemble des ponctuels forts. Ledit sous-moyen de calcul principal est alors apte à calculer l'erreur de phase de l'image numérique à partir des erreurs de phase partielles moyennes. Selon un mode de réalisation, le dispositif peut comprendre en outre un moyen de correction de l'image numérique défocalisée, monté en sortie du moyen de détermination et comprenant un multiplieur apte à multiplier un signal temporel représentatif de l'image numérique défocalisée par le conjugué d'un signal représentatif de l'erreur de phase de ladite image numérique défocalisée. Selon un mode de réalisation, le dispositif peut comprendre en outre une boucle de retour couplant la sortie du moyen de correction et une entrée du moyen de formation des imagettes via un module d'ajustement de la valeur du nombre d'imagettes. Selon un autre aspect, il est proposé une utilisation d'un dispositif tel que décrit ci-avant, au sein d'un radar comprenant une antenne formant ledit moyen d'éclairement, l'image numérique défocalisée étant alors une image de type radar. D'autres avantages et caractéristiques de l'invention apparaîtront à l'examen de la description détaillée de modes de réalisation et de mise en oeuvre de l'invention nullement limitatifs, et des dessins annexés sur lesquels La figure 1 illustre le lobe principal d'un faisceau émis par une antenne radar ; La figure 2 illustre un mode de mise en oeuvre d'un procédé de traitement d'une image numérique selon l'invention ; La figure 3 illustre plus en détail un mode de mise en oeuvre d'une étape de détermination de l'erreur de phase de l'image SAR et de sa correction ; La figure 4 illustre plus en détail un mode de mise en oeuvre d'une étape de formation d'imagettes ; La figure 5 illustre différents signaux mis en oeuvre entre la réception d'un écho et la formation des imagettes ; La figure 6 illustre un exemple d'erreur de phase d'une image défocalisée pour un temps d'éclairement Te ; la figure 7 illustre un mode de réalisation d'un dispositif de traitement d'image selon l'invention ; la figure 8 illustre un mode de réalisation d'un moyen de formation d'imagettes selon l'invention ; la figure 9 illustre plus en détail un mode de réalisation d'un moyen de détermination d'une erreur de phase selon l'invention ; - la figure 10 illustre plus en détail un mode de réalisation du module de calcul secondaire et du sous-moyen de calcul, selon l'invention ; et la figure 11 illustre plus en détail un mode de réalisation d'un moyen de correction selon l'invention. On se réfère à la figure 1 qui représente de façon schématique le lobe principal LBP d'un faisceau émis par une antenne radar d'un porteur AV, par exemple un avion. La trajectoire TJ de l'avion AV est symbolisée par une flèche. Dans la réalité, le faisceau comporte également des lobes secondaires. Ces derniers ne sont pas représentés à des fins de simplification du schéma. On se place dans cet exemple dans le cadre d'un mode télescope, c'est-à-dire que l'antenne est asservie à la zone éclairée, ici une zone ZN, quelque soit la position du porteur. Le faisceau est émis pendant un temps d'éclairement Te supérieur à 30s. Le faisceau éclaire la zone ZN à l'intérieur de laquelle ont été représentés, à titre d'exemple, deux ponctuels pct0 et pctl. Les ponctuels sont des points particuliers de la zone ZN d'où sont renvoyés des échos vers l'antenne radar. Typiquement, les ponctuels correspondent à des arrêtes ou coins de bâtiments, ou encore à des zones particulièrement réfléchissantes tels que des trièdres métalliques disposés au sol, sur la zone éclairée, dans le but de renvoyer des échos aisément identifiables. On localise chaque ponctuel pct0 et pctl à l'aide de leur gisement selon l'axe de l'azimut, respectivement GO et G1. On rappelle que l'azimut est l'angle horizontal entre la direction d'un objet (ici un ponctuel) et une direction de référence (ici la trajectoire TJ de l'avion AV). Les ponctuels pct0, pctl représentés ici sont des ponctuels dits forts, c'est-à-dire que les échos renvoyés par ces éléments sont remarquables. Sur l'image radar représentant la zone ZN, l'existence de ponctuels forts se traduit visuellement par l'apparition de points particulièrement brillants. Le spectre fréquentiel des ponctuels forts possède des raies de forte amplitude en comparaison aux raies des autres ponctuels. Il est donc particulièrement intéressant de privilégier le traitement de ces ponctuels forts. Comme indiqué en préambule, une image radar directement formée à partir des informations issues de l'antenne, est défocalisée. Autrement dit, les contours des objets représentés sur l'image sont flous. On rappelle que les principales raisons de cette défocalisation sont l'effet Doppler et l'imprécision des capteurs inertiels embarqués à bord du porteur. Le signal temporel d'un écho renvoyé par un ponctuel et en particulier par un ponctuel fort est caractérisé entre autres, par sa phase qui traduit la défocalisation de l'image radar. Le signal temporel s(t) d'un écho peut être modélisé par : s(t)=A.rect[eJ.exp[i(2xx fo xt+O(t))], où s(t) est le signal temporel représentatif de l'écho, A est l'amplitude du signal, rect[] est la fonction rectangle, fo est la fréquence de Doppler du ponctuel considéré, et - 0(t) est la loi de phase de défocalisation du ponctuel considéré, encore appelée erreur de phase. On rappelle que la fréquence de Doppler correspond à l'écart de fréquence entre l'écho et le signal émis par l'antenne. A présent, un mode de mise en oeuvre d'un procédé de traitement d'image défocalisée, va être décrit. Ce mode de mise en oeuvre comprend la détermination à partir de l'étude de la loi de phase d'un ponctuel fort d'une image traitée, dans le but ici de corriger cette défocalisation. L'organigramme de la figure 2 présente différentes étapes de ce mode de mise en oeuvre. Tout d'abord, lors d'une première étape 1, on reçoit les échos renvoyés par les ponctuels appartenant à la zone éclairée par le radar. On considère ici que l'image radar est de haute résolution et par conséquent, que le temps d'éclairement Te est relativement long (par exemple supérieur à 10s), cas pour lequel ce mode de mise en oeuvre est particulièrement adapté. Chaque écho est alors numérisé 2 et on forme à partir des échos numérisés une image radar 3 dite SAR en référence à la technique SAR décrite en préambule. Cette image est de haute résolution : la résolution peut par exemple être comprise entre 15 cm et 50 cm. Cette image SAR est défocalisée pour les raisons évoquées ci-dessus. On effectue une transformée de Fourier inverse de l'image SAR venant d'être formée 4 de manière à élaborer le spectre de l'image. On détermine alors l'erreur de phase de l'image SAR 5, à partir de son spectre. Une fois cette erreur de phase connue, on corrige la défocalisation 6 de l'image SAR, de façon que les contours des différents objets représentés soient nets. Pour déterminer l'erreur de phase de l'image SAR considérée, le mode de mise en oeuvre illustré sur la figure 2 propose de s'attacher aux ponctuels forts de l'image (ayant une brillance particulièrement remarquable au sein de l'image radar considérée). En d'autres termes, on détermine l'erreur de phase de l'image SAR défocalisée en étudiant l'erreur de phase d'un ponctuel fort, voire de plusieurs ponctuels forts. Cependant, au lieu de calculer l'erreur de phase d'un ponctuel directement à partir de l'image SAR défocalisée, on forme à partir de cette image SAR, des imagettes ayant une résolution inférieure à celle de l'image SAR initiale. Puis on détermine des erreurs de phase partielles du ou des ponctuels forts au sein de ces imagettes, avant de calculer l'erreur de phase de l'image SAR à partir des erreurs de phase partielles des ponctuels. Plus précisément, comme indiqué sur la figure 3 pour déterminer l'erreur de phase globale, on commence par former N imagettes à partir de l'image SAR traitée (N étant un entier). Ces imagettes diffèrent de l'image SAR traitée, de par leur résolution, N fois inférieure à la résolution qu'aura l'image numérique une fois traitée. Le choix de la valeur de N est effectué de manière que la défocalisation de chaque imagette ne soit pas trop importante. En effet, une trop forte défocalisation implique un étalement de l'énergie, et par conséquent une plus grande difficulté à détecter les points brillants au sein des imagettes. Typiquement, une défocalisation acceptable correspond à une erreur de phase de l'imagette inférieure à 200°. Mais avant tout, la valeur de l'entier N est fonction de la résolution de l'image SAR. On rappelle que la résolution angulaire est l'angle minimum entre deux objets situés à une même distance du radar pour qu'ils ne soient pas confondus. La résolution dépend de la distance entre l'antenne radar et l'objet éclairé. Ces imagettes correspondent en fait à des images obtenues avec une portion du temps d'éclairement Te, soit Te/N lorsque N imagettes sont formées. Chacune de ces portions est définie par son centre tk, tel que : Nù11 Te tk=\kù 2 xN-, où: - N est le nombre d'imagettes formées, - Te est le temps d'éclairement utilisé pour former l'image SAR, et - k est un entier variant entre 0 et N-1. En variante, chaque imagette peut recouvrir au moins une imagette adjacente. En d'autres termes, la durée de la portion du temps d'éclairement générant une imagette de rang k (k étant un entier) recouvre la durée de la portion du temps d'éclairement ayant généré l'imagette de rang k-1. Le chevauchage des imagettes entre elles, permet notamment d'atténuer les artéfacts liés au processus de formation des imagettes. Par exemple les imagettes peuvent être formées selon le mode de mise en oeuvre illustré sur la figure 4. A l'issue de la réalisation de la transformée de Fourier inverse du spectre de l'image SAR, 4, on applique des fonctions rectangle décalées de Te/N de façon à obtenir des signaux correspondant à des images résultant d'un éclairement égal Te/N, 510. Enfin on effectue une transformée de Fourier de chaque signal issu de l'application d'une fonction rectangle, 511. La figure 5 illustre les différents signaux mis en jeux entre la réception d'un écho, et la formation des imagettes. Un exemple d'écho ECH est représenté sur le repère A. Cet écho est numérisé : on obtient alors un signal discret s(nT), repère B. On forme alors l'image SAR et on détermine le spectre S(w) du signal numérisé s(nT), w correspondant à une pulsation repère C. On effectue une transformée de Fourier inverse sur le spectre S(w) de façon à trouver le signal s(t), repère D. Suivent ensuite les signaux liés à la formation des imagettes proprement dit. Sur le signal s(t), on applique les fonctions rectangle avec un décalage temporel choisi. Dans cet exemple, on considère que N=2. Par conséquent, on applique deux fonctions rectangle sur le signal s(t), respectivement comprises entre [0, T/2] et [Te/2, Te]. Deux signaux so(t) et si(t) sont issus de l'application des deux fonctions rectangle, repère E. On effectue une transformée de Fourier sur chacun des deux signaux so(t) et si(t) de façon à obtenir les deux imagettes IMGO(f) et IMGI(f) proprement dites, repère F. On se réfère à nouveau à la figure 3. Une fois les imagettes formées, on s'attache à localiser puis à extraire au moins un ponctuel fort au sein de chacune des imagettes 52. La localisation d'un ponctuel fort peut se faire de manière classique à l'aide de méthode de détection de contraste, de façon à repérer les raies de plus forte amplitude au sein de l'imagette considérée. Ensuite, on extrait une fenêtre autour du ponctuel fort localisé et on calcule sa transformée de Fourier inverse. Dans la suite de la description, sauf indication contraire, on considère que l'on n'extrait qu'un seul ponctuel fort. Mais il est possible en variante d'extraire plusieurs ponctuels forts auxquels on applique le même processus. On calcule alors l'erreur de phase de l'image défocalisée à l'aide de N erreurs de phase partielles moyennes. Une erreur de phase partielle moyenne correspond de la moyenne des erreurs de phase partielles obtenues pour l'ensemble des ponctuels forts, sur une imagette. En considérant les erreurs de phase partielles moyennes, on s'affranchit des déformations parasites du signal, liées aux erreurs de mesure ou autres. Soit O(t) l'erreur de phase d'un ponctuel fort de la zone éclairée, et par conséquent l'erreur de phase de l'image SAR défocalisée, compte tenue de l'hypothèse d'extraction d'un seul ponctuel fort, mentionnée ci-dessus. On note 0k(t)la composante de cette erreur de phase sur une imagette k (05ksN-1), appelée erreur partielle du ponctuel fort. II vient alors : 0k (t) _ 0(t) x rect t ù tk Te' ' où Te' =Te/N. Si l'origine des temps de chaque imagette est recentrée sur son propre milieu, l'erreur de phase du ponctuel fort dans ce nouveau repère Wk (t) peut alors s'écrire : ~k (t) = 0k (t + tk) = o(t + tk) x rect[ Te' ] Par conséquent, l'erreur de phase du ponctuel fort peut s'écrire : 0(t) = 0k (t) = yrk (t ù tk) . k k En d'autres termes, la valeur de l'erreur de phase d'un ponctuel, ici un ponctuel fort, est égale à la somme de ses composantes (ou erreurs de phase partielles) sur chaque imagette. Un exemple d'une erreur de phase q5(t)de type parabolique et de ses composantes Ok(t), est représenté sur la figure 6. L'image SAR est ici découpée en quatre imagettes IMGO, IMGI, IMG2 et IMG3, chacune étant respectivement centrée sur les instants to, t2 et t3. Chaque imagette correspond à une image obtenue pour un temps d'éclairement Te/4. Les erreurs de phase partielles du ponctuel fort sont référencées 00(0 , çbl (t) , 02 (t) , 03 (t) respectivement pour les imagettes IMGO, IMG1, IMG2 et IMG3. L'erreur de phase 0(t) correspond à la succession des erreurs de phases partielles du ponctuel fort, mises bout à bout. Comme on peut le voir sur la figure 6, la dernière valeur d'une erreur de phase partielle correspond à la première valeur de l'erreur de phase partielle lui succédant, en raison de la continuité temporelle de l'erreur de phase de l'image. On se réfère à nouveau à la figure 3. Si l'on effectue un développement de Taylor de la composante de l'erreur de phase sur la k1ème imagette vk (t) , k appartenant à [O, N-1], il vient : Yak (t) ù Y1k (0) + Y/k' (0) X t + CSI//k (t) , ou Tl T2 Svk(t) ne contient que des ordres >_ 2, et t E [ùTe'/2,Te'/2]. Le premier terme Ti traduit la formation de l'imagette considérée, de rang k ici. Le second terme T2 traduit quant à lui, la défocalisation de l'image SAR proprement dite. Ainsi, si l'on considère l'ensemble des imagettes des rangs 0 à N-1, les erreurs de phase partielles du ponctuel fort forment le système d'équation suivant : Wo (t) = vo (0) + teo' (0) x t + 6'vo (t) ~l (t) = Yil (0) + lvl' (0) X t + sçvl (t) (S). ~N 1 (t) = l/'N-1 (0) + YIN-l'(0) X t + UIVN-1 (t) Par conséquent, pour déterminer l'erreur de phase de l'image SAR défocalisée, il faut résoudre le système de N équations (S) comportant 3N inconnues, chaque équation correspondant à l'expression de l'erreur de phase partielle du ponctuel fort pour une imagette donnée. On s'attache tout d'abord à estimer pour chaque imagette, le deuxième terme T2 traduisant la délocalisation de l'image SAR, 53. Pour chaque imagette, le terme 8i'k (t) s'écrit : 8yrk(t)-8(t)-2xgxA Comme expliqué plus en détail ci-dessous, les termes v'k (0) traduisent le décalage entre deux imagettes : ce décalage peut être déterminé, à une constante près, par la mesure de la position du pic de corrélation résultant de la corrélation entre ces deux imagettes. On choisi tout d'abord une imagette de référence, par exemple l'imagette centrale à laquelle on attribue l'indice kfef. Puis on effectue une corrélation de chaque imagette avec l'imagette de référence, de façon à mesurer la position du pic de corrélation traduisant le décalage entre l'imagette en cours et l'imagette de référence. Sur la figure 6, ce décalage correspond à la pente de la droite tangente à la portion de parabole 0k (t) en son centre tk. Pour calculer la position du pic de corrélation, on commence tout d'abord à calculer un signal c(t), tel que : C(t)=Sk(t)XSk,.ef*(t), ou encore, si l'on considère les imagettes IMGk et IMGkref, c(t) = FFT-' (IMGk) x FFT' (IMGkref) , Où FFT est l'opérateur de la transformée de Fourier. La corrélation est donnée ensuite par : C(f) = FFT [c(t)] . Dans tous les cas, le terme de ((phase du signal c(t) vaut : Vk(t)ùWkref(t), soit en développant : ~k (t) ù Y'kref (t) = Wk (0) ù ykref (0)]+ bek' (0) ù yfkref�°)]X t + ['5k (t) ù 8Vkref (t)] La position du pic de corrélation est donnée par le coefficient des termes de premier ordre. L'expression [8yfk(t)ù8yfkref(t)] ne contient que des termes d'ordre supérieur ou égal à deux. Elle influe donc sur la largeur du pic de corrélation et non sur sa position. L'expression [yrk (0) û yrl.ef(0) j est une constante : elle ne contient donc pas de terme du premier ordre. Par conséquent, la position Ak, kref du pic de corrélation entre l'imagette IMGk et IMGkref est liée aux termes yrk'(0), yrkref'(0) par : vk'(0)-Vkref'(0)=2XrcXAk,kref OÙ Ak,kref est exprimé en fréquence. En effectuant une corrélation entre chaque imagette IMGk, k étant différent de kref, et l'imagette de référence IMGkref, on réalise en tout N-1 corrélations. A partir des différents termes calculés précédemment, on peut calculer l'erreur de phase de l'image SAR défocalisée. Ce calcul est effectué en tenant compte de deux hypothèses concernant l'erreur de phase pour t=0, ainsi que d'un critère de continuité entre les erreurs de phases partielles. Tout d'abord, on considère que la phase absolue de l'erreur de phase n'a pas d'influence sur la qualité de la correction de la défocalisation de l'image. On fixe alors : 0(0) = 0. Par ailleurs, on souhaite qu'en compensant la défocalisation de l'image à l'aide de l'erreur de phase calculée, il n'y ait pas de décalage de points dans l'image, ce qui se produit lorsque l'erreur de phase calculée pour corriger la défocalisation de l'image comporte une partie linéaire. Par conséquent, on pose : aO(t) at En outre, comme on peut le voir sur la figure 6, il existe une continuité de l'erreur de phase q5(t) sur l'ensemble du temps d'éclairement Te. II vient alors les N-1 équations suivantes : =0. ,=0 r Te" _ ( Te' " ~Vo 12 ) 2 \ \ 1 (-- 2 /1 Vl 'J = ~Vz2e'J C CTe'~ = Cù Te'1 iVN-2 ù VINù1 2 2 J On peut alors construire l'erreur de phase de l'image SAR défocalisée en mettant bout à bout les erreurs de phases partielles tout en tenant compte du critère de continuité ci-dessus. On calcule ainsi l'erreur de phase de l'image SAR défocalisée à une constante et une phase linéaire près. En effet, on sait calculer l'expression suivante appelée dans la suite du texte Ok(t) : Ok(t)= ,/k(t+tk)ùlek(0)-v1'kref (0)xt, pour t[-Te'/2,Te'/2]. La mise � bout-à-bout � les erreurs de phases partielles pour le calcul de l'erreur de phase donne alors : O(t)=Ok(t+tk)xrect tùtk k Te Soit, en utilisant le critère de continuité précité : O(t) = E Bk (t) ù E iV'kref (0) x (t ù tk) x rectl t Te tk + cste = 9(t) ù yr'kref (0) x t + cste avec t e [û Te'/ 2, Te'/ 2]. Finalement, après un éventuel filtrage pour lisser l'erreur de phase calculée, on choisit la valeur de la constante et la phase linéaire yr'kref (0) , de façon à obtenir la valeur de la constante indéfinie de l'erreur de phase. On peut alors corriger l'image SAR défocalisée 56, en multipliant la transformée de Fourier inverse de celle-ci avec le conjugué de l'erreur de phase calculée. k k II est possible de réitérer le processus (51-57) sur l'image venant d'être corrigée après avoir décrémenté la valeur de N tant que N>1 avant la décrémentation. De cette façon, on affine au maximum la correction apportée. Dans le cas où N=1, l'image n'est quasiment plus défocalisée. II est alors possible de simplifier l'étape de détermination de l'erreur de phase de l'image 5 en appliquant directement une méthode dite PPP ou PGA sur l'image radar elle-même. On va à présent décrire un mode de réalisation d'un dispositif de traitement d'image selon l'invention, pouvant être réalisé de façon logicielle. Ce dispositif peut être incorporé par exemple au sein d'un système radar (exemple pris dans la suite du texte) ou encore dans un sonar. On se réfère à la figure 7. La référence DIS désigne un dispositif de traitement d'une image numérique. Ce dispositif est couplé à une antenne ANT, ici une antenne radar. Cette antenne ANT émet un signal d'éclairement SEM pendant une durée Te vers une zone au sol, de façon à pouvoir former une image radar de la zone éclairée. Certains points de la zone (arrêtes de bâtiment par exemple) renvoient des échos plus intenses vers l'antenne ANT. A des fins de simplification, on considérera que la zone éclairée ne renvoie qu'un seul écho ECH. Cet écho ECH est reçu par l'antenne ANT et numérisé par un module de numérisation NUM. Ce dernier délivre un signal discret s(nT) à un module MSAR de formation d'image radar par exemple de type SAR, ayant pour spectre fréquentiel S(w), w étant une pulsation. Un module FFT-' applique une transformation de Fourier inverse à l'image SAR formée S(w) et délivre en sortie un signal temporel s(t) représentatif de l'image SAR. Le module FFT"' est couplé à un moyen de détermination de l'erreur de phase de l'image SAR formée MDET. Celui-ci comporte un moyen de formation d'imagettes MFI et un moyen de calcul de l'erreur de phase MCAL. Plus précisément, le signal s(t) est délivré au moyen de formation d'imagettes MFI, apte à former N imagettes à partir du signal s(t). Ces N imagettes ont une résolution N fois inférieure à celle de l'image radar formée par le module MSAR, comme expliqué ci-avant. On considère dans cet exemple que le moyen de formation d'imagettes MFI délivre deux imagettes IMGO(f) et IMGI(f). Sauf indication contraire, on considérera dans cet exemple que seules deux imagettes sont formées. Ces deux imagettes IMGO(f) et IMG1(f) sont délivrées à un moyen de calcul MCAL, apte à déterminer l'erreur de phase 0(t) de l'image SAR, cette erreur causant une défocalisation de l'image SAR. L'erreur de phase 0(t) déterminée est délivrée à un moyen de correction MCOR couplé en sortie du moyen de détermination MDET et capable de corriger le signal s(t) représentatif de l'image SAR. A l'aide d'une boucle de retour BR reliant la sortie du moyen de correction MCOR et une entrée du moyen de formation des imagettes MFI, le signal corrigé sfoc(t) est renvoyé en entrée du moyen de formation des imagettes MFI via un module d'ajustement MAJU du nombre d'imagettes N. En effet de façon à affiner la correction de la défocalisation, le signal corrigé sfoc(t) est retraité de façon à obtenir une image radar le plus net possible. On diminue alors le nombre d'imagettes formées à chaque itération, jusqu'à ce que N soit égal à 1. La figure 8 illustre un mode de réalisation d'un moyen de formation des imagettes MFI. Dans ce mode de réalisation, le moyen de formation MFI reçoit le signal s(t) représentatif de l'image radar de type SAR. Un module RECT applique N (ici N=2) fonctions rectangle décalées de Te/N sur le signal s(t) de façon à former N (ici N=2) signaux temporels so(t) et si(t). Un module FFT permet d'appliquer une transformée de Fourier à ces deux signaux pour former les imagettes proprement dit, respectivement IMGO(f) et IMG1(f). La figure 9 illustre un mode de réalisation du moyen de calcul MCAL de l'erreur de phase q(t). Ce moyen comprend un sous-moyen de détection de ponctuels forts DTP, un sous-moyen d'estimation SMEST de l'erreur de phase partielle de chaque imagette, et un sous-moyen de calcul SMCAL, apte à calculer l'erreur de phase de l'image O(t) à partir des erreurs de phases partielles estimées. Les imagettes IMGO(f) et IMGI(f) sont reçues par le sous-moyen de détection de ponctuels forts DTP. Celui-ci comprend un module de localisation LCL des ponctuels au sein des imagettes. Par exemple, un algorithme de détection de contraste peut être implémenté pour la localisation des ponctuels forts. Une fois localisés, les ponctuels forts sont extraits à l'aide d'une fenêtre appliquée par un module de fenêtrage MFN couplé en sortie du moyen de localisation LCL. Dans cet exemple, sauf indication contraire, on considère qu'un seul ponctuel fort est détecté par imagette, respectivement IMGO(f0) et IMGI (fO). Ces ponctuels forts IMGO(f0) et IMGI(f0) sont délivrés à un premier module MPGA apte à déterminer les deuxièmes termes 8vo(t) et 8y'1(t) de l'erreur de phase partielle de chaque imagette à l'aide par exemple d'un algorithme dit PPP ou PGA bien connus de l'homme du métier. En parallèle, les ponctuels forts IMGO(f0) et IMGI(fl) sont délivrés à un module d'intercorrélation CORL (deuxième module) apte à mesurer la position des pics de corrélation Ao,ref A1,ref , issues de la corrélation entre chaque imagette et une imagette de référence indicée kref. Ces positions des pics de corrélation Ao,ref A1.ref sont délivrés à un multiplieur par 27c afin de déterminer le premier terme de chaque erreur de phase partielle à une constante et une partie linéaire près, Vo'(0) ûVik,,ef'(0) et vl�0) ûtgkref '(0). Chaque élément issu du module MPGA et du multiplieur 2it, est délivré en entrée d'un module de calcul secondaire CALS. Celui-ci estime l'erreur de phase partielle de chaque imagette pour le ponctuel fort extrait IMGO(f0) et IMGI(f0). Les erreurs de phase partielles estimées sont transmises en entrée du sous-moyen de calcul SMCAL qui calcule alors l'erreur de phase de l'image numérique défocalisée, 0(0. Un module de moyennage pourrait être inséré entre le module de calcul secondaire CALS et le sous-moyen de calcul SMCAL, dans le cas où on localise et extrait plusieurs ponctuels forts. Dans cette hypothèse, le sous-moyen de calcul SMCAL calcule l'erreur de phase de l'image numérique défocalisée à partir de N erreurs de phase partielles moyennes. Chaque erreur de phase partielle moyenne correspond à une moyenne des erreurs de phase partielles estimées pour l'ensemble des ponctuels forts sur une imagette. La figure 10 représente plus en détail un mode de réalisation du module de calcul secondaire CALS et du sous-moyen de calcul SMCAL. Le module de calcul secondaire CALS comprend deux intégrateurs INT de façon à élaborer la partie linéaire de chaque erreur de phase partielle à l'aide des termes vo'(0) û vkref'(0) et v,'(0) û yrkref'(0) . Les termes linéaires obtenus sont respectivement sommées au deuxième terme correspondant 8vo (t) et 8yr, (t) à l'aide de deux sommateurs SOM. Ces derniers délivrent alors les erreurs partielles de chaque imagettes à une constante et un terme linéaire près, Oo (t) et 01(0. Un moyen de construction MCST incorporé dans le sous-moyen de calcul SMCAL, élabore l'erreur de phase de l'image à une constante et une partie linéaire près O(t) en mettant bout à bout les erreurs de phase partielles Oo (t) et O, (t) . Pour ce faire, le moyen de construction MCST tient compte d'un critère de continuité temporelle CNT entre les erreurs de phases partielles. Un module de réglage REG couplé en sortie du moyen de construction MCST, fixe la constante et la partie linéaire de l'erreur de phase 0(t) de façon à délivrer l'erreur de phase finale de l'image q(t) au moyen de correction MCOR. La figure 11 représente un mode de réalisation d'un moyen de correction MCOR. Dans cet exemple, il comprend un moyen CONJ apte à déterminer le conjugué 0` (t) de l'erreur de phase 0(0. Un multiplieur MUL multiplie le conjugué obtenu 0*(t) par le signal s(t) représentatif de l'image défocalisée, de façon à la corriger. Le dispositif décrit ci- avant peut être incorporé au sein d'un système radar ou encore d'un système sonar. Dans toute la description, il a été question d'images monodimensionnelles. L'homme du métier saura très facilement adapter l'invention à des images bi-dimensionnelles. REVENDICATIONS 1-Procédé de traitement d'une image numérique défocalisée, obtenue par l'éclairement d'une zone (ZN) pendant un temps donné (Te), comprenant une détermination (5) de l'erreur de phase de ladite image numérique défocalisée, ladite détermination comprenant : - une formation (51) à partir de ladite image numérique défocalisée, de N imagettes, N étant un entier, chaque imagette ayant une résolution N fois inférieure à celle de l'image numérique une fois traitée, et - un calcul (55) de ladite erreur de phase de l'image numérique défocalisée, caractérisé par le fait que le nombre N d'imagettes formées est fonction de la résolution de ladite image numérique défocalisée. 2-Procédé selon la revendication précédente, dans lequel la détermination (5) de ladite erreur de phase comprend pour chaque imagette, une localisation puis une extraction (52) d'au moins un point remarquable dit ponctuel fort, une estimation pour au moins un desdits ponctuels forts extraits, de son erreur de phase partielle (53, 54), ledit calcul de l'erreur de phase de l'image numérique défocalisée étant effectué à partir de l'ensemble des erreurs de phase partielles. 3- Procédé selon la revendication précédente, dans lequel chaque erreur de phase partielle comprend un premier terme traduisant la formation de l'imagette considérée et un deuxième terme traduisant la défocalisation de l'image numérique. 4-Procédé selon l'une des revendications précédentes, dans lequel pour chaque imagette : - on localise et extrait au moins deux ponctuels forts, - on estime son erreur de phase partielle pour chaque ponctuel fort localisé et extrait,- on calcule une erreur de phase partielle moyenne en effectuant une moyenne des erreurs de phase partielles estimées pour l'ensemble des ponctuels forts, et - on calcule l'erreur de phase de l'image numérique défocalisée à partir des erreurs de phase partielles moyennes. 5-Procédé selon l'une des revendications précédentes, comprenant en outre une correction (56) de l'image numérique défocalisée en multipliant le signal temporel représentatif de ladite image numérique défocalisée par le conjugué du signal représentatif de l'erreur de phase déterminée de ladite image numérique défocalisée. 6-Procédé selon l'une des revendications 3 à 5, dans lequel le premier terme de chaque erreur de phase partielle comprend un coefficient de phase traduisant le décalage entre une imagette considérée et une imagette de référence, et une constante de phase, ledit coefficient de phase étant déterminé par une intercorrélation de l'imagette considérée et de l'imagette de référence, et ladite constante de phase étant déterminée en considérant la continuité temporelle de l'erreur de phase partielle entre deux imagettes successives. 7-Procédé selon l'une de revendications 5 à 6, dans lequel on itère l'étape de détermination de l'erreur de phase de l'image numérique défocalisée au plus N fois, en décrémentant (57) le nombre d'imagettes formées à chaque répétition, chaque nouvelle itération étant effectuée sur l'image numérique défocalisée corrigée avec l'erreur de phase déterminée lors de la précédente itération. 8- Procédé selon l'une des revendications précédentes, chaque imagette formée recouvre partiellement au moins une imagette adjacente. 9- Dispositif (DIS) de traitement d'une image numérique défocalisée obtenue à l'aide de moyens d'éclairement d'une zone pendant un temps donné, comprenant un moyen de détermination de l'erreur de phase de l'image numérique défocalisée qui comporte :- un moyen de formation (MFI) apte à former à partir de ladite image numérique défocalisée, N imagettes, N étant un entier, chaque imagette ayant une résolution N fois inférieure à celle de l'image numérique défocalisée, - un moyen de calcul (MCAL) couplé audit moyen de formation et apte à déterminer l'erreur de phase de l'image numérique défocalisée, à partir des N imagettes, caractérisé par le fait que le nombre N d'imagettes formées par le moyen de formation est fonction de la résolution de ladite image numérique défocalisée. 10- Dispositif de traitement selon la revendication précédente, dans lequel le moyen de calcul comprend : -un sous-moyen de détection (DTP) de ponctuel apte à localiser et à extraire au moins un point remarquable dit ponctuel fort au sein des imagettes, - un sous-moyen d'estimation (SMEST) apte à estimer pour chaque imagette, l'erreur de phase partielle de chaque ponctuel fort extrait, et - un sous-moyen de calcul (SMCAL) apte à calculer l'erreur de phase de l'image numérique défocalisée à partir de l'ensemble des erreurs de phase partielles estimées. 11-Dispositif selon l'une des revendications 9 à 10, dans lequel le sous-moyen d'estimation (SMEST) comprend un premier module (MPGA) apte à estimer un premier terme traduisant la formation de l'imagette considérée et un deuxième module (CORL) apte à estimer un deuxième terme traduisant la défocalisation de l'image numérique. 12- Dispositif selon la revendication précédente, dans lequel le deuxième module (CORL) est apte à estimer le deuxième terme en effectuant une intercorrélation entre une imagette et une imagette de référence. 13-Dispositif selon l'une des revendications 10 à 12, dans lequel le sous-moyen de détection (DTP) est apte à localiser et à extraire au moinsdeux ponctuels forts, dans lequel le sous-moyen de détermination comprend en outre un module de moyennage apte à calculer pour chaque imagette, une erreur de phase partielle moyenne à partir des erreurs de phase partielles estimées pour l'ensemble des ponctuels forts, ledit sous-moyen de calcul (SMCAL) étant alors apte à calculer l'erreur de phase de l'image numérique à partir des erreurs de phase partielles moyennes. 14-Dispositif selon l'une des revendications 10 à 13, comprenant en outre un moyen de correction de l'image numérique défocalisée, monté en sortie du moyen de détermination (MCOR) et comprenant un multiplieur (MUL) apte à multiplier un signal temporel représentatif de l'image numérique défocalisée par le conjugué d'un signal représentatif de l'erreur de phase de ladite image numérique défocalisée. 15-Dispositif selon la revendication précédente, comprenant en outre une boucle de retour (BR) couplant la sortie du moyen de correction (MCOR) et une entrée du moyen de formation (MFI) des imagettes via un module d'ajustement (MAJU) de la valeur du nombre d'imagettes. 16-Utilisation d'un dispositif selon l'une des revendications 9 à 15, au sein d'un radar comprenant une antenne (ANT) formant ledit moyen d'éclairement, l'image numérique défocalisée étant alors une image de type radar.