FR 2890451 A1 20070309 FR 0509000 A 20050902 DOMAINE TECHNIQUE GENERAL La présente invention concerne un procédé de détection d'un objet situé dans un champ optique d'un système optronique. Elle concerne également un dispositif de détection mettant en oeuvre le procédé. ETAT DE L'ART Comme le montre la figure 1, la détection de missiles, dans le cadre de I'autoprotection des avions 10 de transport ou de chasse, se situe dans le cadre général de la détection d'un objet 15 quasi ponctuel par un système optronique 11 monté sur les avions 10. En général, le système optronique 11 est un imageur, composé d'une 15 optique 112 et d'un détecteur matriciel 111. Pour être efficace, la détection des missiles 15 doit être omnidirectionnelle. Il en résulte tout d'abord que pour diminuer le nombre de systèmes optroniques 11, un champ 16 de l'optique 112 est le plus grand possible. La résolution du système en subit les conséquences: du point de vue de l'imagerie, la scène est en générale sous résolue. Il en résulte ensuite que la détection doit s'opérer sur tous types de fonds, à savoir des fonds uniformes de type ciel clair , des fonds de ciel structurés ou nuageux, ou des fonds de sols urbains, ruraux ou naturels. Les missiles 15 devant être détectés sont des missiles sol-air ou air-air. Ce sont des objets quasi ponctuels par rapport à l'avion 10 ou à l'environnement observé dans le champ 16, et fortement mobiles dans le référentiel terrestre. Leur vitesse est de l'ordre de Mach 2,5. Le profil temporel de l'émission infrarouge de ces missiles est 30 constitué de plusieurs phases, qui correspondent à différentes phases de vol. Dans le cas des missiles SATCP (Sol Air Très Courte Portée), la figure 2 montre que le profil temporel de l'émission comporte: - une phase 1 d'éjection du missile 15 d'un tube 14 (pour les missiles 15 SATCP tirés par un fantassin 13), - une courte phase 2 de vol balistique (quelques centaines de ms), - une mise à feu 3 d'un réacteur du missile (ou over shoot selon la terminologie anglo-saxonne généralement utilisée par l'homme du métier), - une phase 4 de propulsion ou de croisière, elle-même décomposée en plusieurs phases suivant les modèles de missile. À la suite de son éjection 1, le missile 15 possède une vitesse de l'ordre de 50 m/s. Puis le missile 15 subit une forte accélération (de 30 à 50 10 g) pendant 1 à 2 secondes (phase 3), pour le mener à une vitesse de l'ordre de 500 à 700 m/s (phase 4). Cette vitesse est maintenue plus ou moins constante pendant 5 à 15 secondes. Avec des dispositifs SATCP, quelques secondes seulement séparent l'éjection 1 du missile 15 de son impact avec l'avion 10, si ce dernier ne 15 déclenche pas de dispositifs de contre-mesure à temps. Par conséquent, il est nécessaire de détecter le missile 15 le plus tôt possible pour que l'avion 10 puisse déclencher ses dispositifs de contremesure. Actuellement, les systèmes optroniques ne détectent pas le missile au 20 plus tôt et ne sont pas optimaux. PRESENTATION DE L'INVENTION L'invention propose de pallier au moins un de ces inconvénients. A cet effet, on propose selon l'invention un procédé de détection d'un objet situé dans un champ optique d'un système optronique, caractérisé en 25 ce qu'il comporte les étapes selon lesquelles: - des moyens de traitement déterminent un spectre spatio-temporel à partir d'un signal issu du système optronique; - un filtre spatio-temporel crée des zones de recherche dans le spectre spatio-temporel pour éliminer les basses fréquences spatiales et les basses 30 fréquences temporelles et; -des moyens de détection détectent dans les zones de recherche des zones d'alarmes dont une magnitude de spectre est supérieure à un seuil, ces zones d'alarmes correspondant à une détection possible de l'objet. 2890451 3 L'invention est avantageusement complétée par les caractéristiques suivantes, prises seules ou en une quelconque de leur combinaison techniquement possible: - le filtre spatio-temporel crée des zones de recherche dans le spectre 5 spatio-temporel en effectuant les étapes suivantes: - un filtrage d'un signal issu du système optronique par un filtre morphologique non linéaire; - un filtrage du signal issu du filtre morphologique non linéaire par un filtre spatial passe-bas; - un filtrage du signal issu du filtre spatial passe-bas par un filtre morphologique non linéaire; - l'étape de filtrage du signal issu du système optronique par le filtre morphologique non linéaire et l'étape de filtrage du signal issu du filtre morphologique non linéaire par le filtre spatial passe-bas sont itérées une pluralité de fois; - le nombre d'itérations est fonction de la sélectivité en fréquence du filtre; - le filtre effectue un filtrage de détection spatiale ponctuelle avant le filtrage par le filtre morphologique non linéaire; - le filtre effectue un filtrage antidistorsion de repliement avant le filtrage par 20 le filtre morphologique non linéaire; - la détermination du spectre spatio- temporel se fait pour chaque signal issu d'un pixel du système optronique; -l'objet est un missile sol-air à très courte portée éjecté à partir d'un dispositif d'éjection; - le système optronique est monté sur un véhicule; - un volume de données composant le signal pour la détermination du spectre spatio-temporel est fonction de la vitesse du véhicule dans un repère terrestre; et - le volume de données est proportionnel à la racine cubique de la vitesse 30 du véhicule dans un repère terrestre. L'invention concerne également un dispositif de détection mettant en oeuvre le procédé. L'invention présente de nombreux avantages. 2890451 4 La technique de détection selon l'invention détecte notamment la phase 1 d'éjection, qui est un événement bref (de l'ordre de 50 ms) et limité dans l'espace. Il s'agit donc d'un événement de haute fréquence du point de vue spatio-temporel. Cette stratégie de détection de l'éjection (ou early- warning selon la terminologie anglo-saxonne généralement utilisée par l'homme du métier) permet de gagner un temps précieux, et de déclencher plus rapidement les dispositifs de contre-mesure de l'avion. La sécurité de l'avion est donc grandement augmentée. La détection de l'objet s'opère sur tous types de fonds, à savoir des 10 fonds uniformes de type ciel clair , des fonds de ciel structurés ou nuageux, ou des fonds de sols urbains, ruraux ou naturels. La technique de détection selon l'invention a un taux de fausses alarmes faible. 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 schématiquement un véhicule du type avion comportant un système optronique dans le champ duquel est 20 situé un dispositif d'éjection d'objet; - la figure 2, déjà commentée, représente schématiquement un profil temporel de l'émission infrarouge d'un missile; - les figures 3 et 4 représentent schématiquement un spectre spatio-temporel d'une image acquise par un système optronique et un spectre spatio-temporel d'une séquence d'images en une dimension acquises par un système optronique; - les figures 5 et 6 représentent schématiquement un spectre spatio-temporel d'une éjection d'un missile acquise par un système optronique; - la figure 7 représente schématiquement des zones de recherche dans un spectre spatio-temporel regroupant les spectres des figures 4 et 6; - la figure 8 est un équivalent de la figure 7 avec plusieurs vitesses de défilement du système optronique; 2890451 5 - les figures 9 et 10 représentent schématiquement un spectre spatio-temporel d'une éjection d'un missile dans différents cas d'échantillonnage du signal; - la figure 11 montre l'évolution d'un spectre spatio-temporel d'une éjection d'un missile en fonction de la vitesse de défilement du système optronique; - la figure 12 représente schématiquement les étapes principales d'un procédé selon l'invention; -la figure 13A représente schématiquement une partie d'un domaine spatio-temporel détecté avec une éjection d'un missile et un point fixe sans 10 appliquer un procédé selon l'invention; - la figure 13B représente schématiquement une partie d'un domaine spatio-temporel détecté avec une éjection d'un missile et deux points mobiles sans appliquer un procédé selon l'invention;et - la figure 14 représente schématiquement le domaine de la figure 13B en 15 appliquant un procédé selon l'invention. Dans l'ensemble des figures, les éléments similaires portent des références numériques identiques. DESCRIPTION DETAILLEE La présente description s'applique avantageusement à un exemple selon la figure 1, à savoir la détection d'un missile 15 pour I'autoprotection d'un avion 10 de transport ou de chasse. Le missile 15 est bien entendu situé dans un champ optique 16 d'un système optronique 11 monté sur l'avion 10. La présente invention s'applique très avantageusement à un missile 15 sol-air à très courte portée tiré à partir d'un dispositif d'émission 14 transporté par un fantassin 13. On comprend cependant que l'invention s'applique également à toute technique de détection d'un objet situé dans un champ optique d'un système optronique, monté sur un véhicule ou un immeuble. Un dispositif de détection selon l'invention comporte principalement le 30 système optronique 11. Le système optronique 11 est un imageur, composé d'une optique 112 et d'un détecteur matriciel 111. 2890451 6 Le champ 16 de l'optique 112 est le plus grand possible, préférentiellement supérieur à 21t stéradian (sr). Le détecteur matriciel 111 comporte une matrice présentant classiquement un ensemble de pixels de détection. Le dispositif de détection comporte de plus des moyens de traitement 12 aptes à déterminer un spectre spatio-temporel à partir d'un signal issu du système optronique 11. Une telle étape de détermination est décrite plus en détail dans la suite de la présente description. Il comporte de plus un filtre spatio-temporel 17 apte à filtrer le spectre 10 spatio-temporel. Une telle étape de filtrage est décrite plus en détail dans la suite de la présente description. Il comporte en outre des moyens de détection 18 aptes à détecter des zones d'alarmes dont une magnitude de spectre est supérieure à un seuil, ces zones alarmes correspondant à une éjection possible de l'objet. Une telle étape de détection est décrite plus en détail dans la suite de la présente description. Lors d'une première étape principale, le système optronique 11 acquiert une séquence d'images. Les images sont prises dans le champ optique 16 du système optronique 11. Lors d'une deuxième étape principale, les moyens de traitement 12 déterminent un spectre spatio-temporel à partir d'un signal issu du système optronique 11 et représentant la séquence d'images. On explique ici pour des raisons de clarté la détermination du spectre spatio-temporel pour des images à une dimension, à savoir un signal issu d'un pixel du système optronique 11. On comprend que la détermination du spectre spatio-temporel se fait en pratique pour chaque signal issu d'un pixel du système optronique 11. Le calcul suivant établit le spectre spatio-temporel d'une séquence d'images S, obtenue par translation d'une même image 1D de position r au 30 cours du temps t, à la vitesse v. S(F,t)=1(ï--i.t) où I est une image définie par le niveau électrique issu du détecteur (dans le plan des réels 932), et où r est le vecteur position dans le détecteur supposé fixe, et où v est le vecteur vitesse de défilement de l'image dans le référentiel du détecteur, résultant du mouvement du porteur du détecteur. Cette équation suppose qu'il existe une échelle suffisamment petite pour que la vitesse v puisse être considérée comme uniforme, et suffisamment grande pour que le déplacement entre deux images soit petit devant cette échelle (pour limiter les effets de bord). La transformation de Fourier de S s'écrit; S(k'w)-ff1(k -vt)'e iki ra> r.d2k-.dt Le changement de variable (':,t) F---> (6,0= (i - t., .t, donne: k s(k,w)= fff 1(l). " e-r(w+ti a}' . d i É dt On obtient finalement: Sk,w)=8(w+k I(k) Les figures 3 et 4 (la figure 4 étant une projection de la figure 3 dans le plan (k co)) montrent que le spectre 31 obtenu pour la séquence S est donc non nul en magnitude (la magnitude est le troisième axe du repère et est perpendiculaire au plan (k,w)) seulement le long d'une droite 33 de pente -v dans le plan (k w). L'enveloppe en magnitude est définie par le spectre 32 de l'image à une dimension pris en la valeur (w/v) (à savoir l'image du spectre par une affinité). Le spectre 32 est nul sur les points de la droite situés au delà du point (Qmax optique, Kmax optique) où Kmax optique est la fréquence de coupure de l'optique du système optronique et S2max optique est égale à v multipliée par Kmax optique. Par ailleurs, l'éjection 1 du missile 15 est un signal de courte durée qui est, bien que ce ne soit pas tout à fait exact, supposé être entraîné par le mouvement du fond de scène. Cette durée brève conduit à un élargissement du spectre temporel de ce signal. L'émission radiométrique du missile 15 peut être modélisée par 30 l'équation: 1=10(.x, v)x f(t) où lo est une fonction dans le plan des réels 93' décrivant la tache optique, c'est à dire la fonction obtenue à partir du détecteur de l'image lorsque le système observe un point à l'infini; et où f(t) est le profil temporel de l'émission (fonction du temps t), c'est à dire une fonction positive qui vaut 0 avant le tir du missile et 1 lorsque celui-ci est le plus visible. x et y désignent les coordonnées dans le détecteur supposé fixe. Il est possible d'écrire: So(i,t)f(t) Io (r -'.t) Le calcul de la transformée de Fourier donne: So(k,w)-JJJ/(t) IO(r VÉt)Ée-rkï e-iur.d'rÉdt Le changement de variable (F't) H (û, t) _ - t, t) donne: So(k'w)-fif.f(t) Io(t) e tk. ïi e d'ûdt Finalement, on obtient: So(k,=.i(w+ i;)'Io(k) Comme le montrent les figures 5 et 6 (de même que précédemment, la figure 6 est une projection de la figure 5 dans le plan (k,w)), contrairement au signal issu du fond, le signal provenant de l'éjection du missile 15 présente un spectre 41 qui est étendu dans l'espace des fréquences a). La figure 7 représente une mise en correspondance du spectre de la figure 4 et du spectre de la figure 6. La mise en correspondance de la figure 7 permet de montrer que le filtre 17 crée, lors d'une troisième étape principale, des zones de recherche 71 dans le spectre spatio-temporel déterminé, pour éliminer les basses fréquences spatiales et les basses fréquences temporelles. Le filtre peut ainsi éliminer notamment la zone correspondant à un spectre 31 du signal issu du fond du paysage en défilement. Il n'est pas possible de détecter une éjection dans la zone du spectre correspondant à la droite 33 ou le spectre 31, car il y aurait un nombre de fausses alarmes trop important. 2890451 9 Les zones de recherche 71 correspondent aux fréquences qui ont les valeurs absolues les plus élevées. Les zones 72 correspondant aux zones du spectre 41 d'où a été soustrait le spectre 31 fréquences de valeur absolue intermédiaire ne sont pas choisies pour la détection. Dans la pratique en effet, il y a dans une séquence d'images, non pas une seule vitesse v, mais un ensemble de vitesses. Dans le champ 16 du système optique 11, les vitesses peuvent prendre un grand nombre de valeurs. Ce phénomène est particulièrement vrai dans le cas des images grands champs 16. La figure 8 montre un spectre (toujours dans le cas à une dimension) pour plusieurs vitesses. On constate que seules les zones de recherche 71 correspondant à des fréquences supérieures aux droites 33 (ces droites 33 étant des cônes dans le cas général) d'équations: w k vmaX permettent de détecter le signal d'éjection dans tous les cas de figure. Les zones 81 correspondent à des zones mixtes où le paysage et l'éjection sont superposés. C'est donc le module de la vitesse maximale qui est la grandeur dimensionnante du filtre, notamment dans le cas où l'on souhaite n'utiliser 20 qu'un seul et même filtre pour toutes les vitesses de l'image. Toutefois l'estimation de la vitesse maximale peut être variable en fonction de la position dans l'image, et être plus ou moins précise si, par exemple, il est possible de connaître le mouvement du porteur, et de tenir compte de cette connaissance pour estimer les mouvements dans l'image. Les développements qui précèdent s'appliquent à un domaine continu, spatialement et temporellement. Les développements qui suivent visent à décrire les effets de l'échantillonnage dans ces deux dimensions. La figure 9 montre un spectre, toujours dans le cas 1 D spatiale, issu d'un échantillonnage vérifiant, dans l'espace et le temps, les conditions de Shannon. La référence 93 montre la pulsation d'échantillonnage temporel et 2890451 10 la référence 94 montre la pulsation d'échantillonnage spatial. On constate qu'on peut déterminer comme précédemment des zones de recherche 71. Dans, la pratique, l'échantillonnage spatial reste invariable. Par contre, la qualité de l'échantillonnage temporel est directement liée à la vitesse de défilement maximale observée dans l'image. La figure 10 montre un spectre, toujours dans le cas 1 D spatiale, issu d'un échantillonnage vérifiant les conditions de Shannon dans l'espace, mais ne les vérifiant pas dans le temps. On constate donc classiquement un repliement du spectre en fréquences. Les zones de recherche 71 sont réduites et décalées vers les basses fréquences. La conséquence du repliement spectral est donc de diminuer la surface des zones de recherche. Il faut par conséquent un filtre 17 plus sélectif en fréquence. La réponse impulsionnelle du filtre est alors moins bien localisée. Ceci nous amène donc à considérer la résolution en fréquence du filtre 17. La résolution en fréquence est directement liée à la dimension de la portion spatio-temporelle considérée pour la détection. Comme le montre la figure 11, lorsque la vitesse v de défilement augmente, la surface du spectre 41 correspondant au signal à détecter diminue. Il apparaît sur la figure 11 que la surface accessible dans le plan (k,w) pour la détection (il s'agira bien entendu d'un volume dans le cas général à deux dimensions spatiales) diminue lorsque la vitesse V augmente. Pour une détection efficace, il faut donc une résolution spatiotemporelle d'autant plus fine que la vitesse est élevée. Dans le cas général on peut évaluer le volume Vol de l'ellipsoïde du spectre 41 dans le plan (k,w) par: Vol = 7r É K11,ax É wmax cos(arctan V Vol = 7r Kmax. wmax ce qui se simplifie en + V' Or cette résolution est donnée par l'inverse d'une dimension 0 typique du volume de données utilisées pour appliquer la détection. La dimension 0 est proportionnelle à Vol''3 Pour maintenir des performances comparables lorsque V varie, il 5 faudrait donc que la dimension 0 varie en sens contraire des dimensions caractéristiques des dimensions du volume Vol calculé plus haut à savoir: La dimension typique 0 du volume de donnée nécessaire pour conserver une qualité de détection varie donc comme vil'. En 10 conséquence, les dimensions de la réponse impulsionnelle associée du filtre suivent la même évolution. Les développements qui précèdent s'attachent à la création des zones de recherche 71. Lors d'une quatrième étape principale, les moyens de détection 18 15 détectent dans les zones de recherche 71 des zones d'alarmes dont une magnitude de spectre est supérieure à un seuil. Ces zones alarmes correspondant à une éjection possible de l'objet. Le dispositif de détection peut donc alors donner une instruction pour le déclenchement des dispositifs de contre-mesure. Les développements qui suivent s'attachent à un mode de mise en ceuvre d'un procédé selon l'invention. La figure 12 représente schématiquement les étapes principales d'un tel procédé selon l'invention. Lors d'une étape 120, le système optronique 11 acquiert une séquence d'images. Les images sont prises dans le champ optique 16 du système optronique 11. Le spectre spatio-temporel 120 est déterminé lors de l'étape 120. Lors d'une étape 121, optionnelle, le filtre 17 effectue un filtrage de détection spatiale ponctuelle. L'étape 122 est avantageusement mise en 2890451 12 oeuvre pour éliminer les objets manifestement trop volumineux pour être des missiles 15. Lors d'une étape 122, optionnelle également, le filtre 17 effectue un filtrage antidistorsion de repliement. L'étape 122 est avantageusement mise 5 en oeuvre dans des conditions de sous-échantillonnage. Un compteur N est ensuite mis à zéro. Lors d'une étape 123, le filtre 17 crée des zones de recherche 71 dans le spectre spatio-temporel 41 en effectuant tout d'abord un filtrage par un filtre morphologique non linéaire. Ce filtre correspond à un filtrage de détection temporelle ponctuelle. Ce filtrage diminue le niveau du signal issu du fond de scène, mais pas celui provenant de l'éjection. Lors d'une étape 124, le filtre 17 effectue ensuite un filtrage du signal issu du filtre morphologique non linéaire par un filtre spatial passe-bas. Ce filtre correspond à un grossissement de pixel. Le fait de grossir la taille du pixel permet de diminuer la proportion de signal hautes fréquences dans le signal sortant du pixel. Ce filtre est avantageusement mis en oeuvre dans les cas où les conditions de faible défilement de la scène ne sont pas réunies. A l'issue de l'étape 124, le compteur N est incrémenté d'un incrément. 20 Lors d'une étape 125, le filtre détermine si la valeur du compteur est inférieure à une valeur maximale Nmax. Si la réponse est oui, alors le filtre 17 retourne à l'étape 123 déjà décrite. Si la réponse est non, le filtre passe alors à une étape 126. On comprend que l'étape 123 et l'étape 124 peuvent être itérées une pluralité de fois. Le nombre d'itérations et le nombre Nmax sont liés à la sélectivité en fréquence du filtre. Plus la vitesse de défilement est importante, plus le filtre doit être sélectif, notamment en fréquence, et plus le nombre d'itérations doit être important. En contrepartie, plus le nombre d'itérations est important, plus la dimension de la réponse du filtre est importante également. Le filtre est d'autant moins localisé qu'il y a d'itérations. Lors de l'étape 126, le filtre 17 effectue enfin un filtrage du signal issu du filtre spatial passe-bas par un filtre morphologique non linéaire. En effet, lors de l'étape 124, il est possible qu'une partie basse fréquence temporelle soit apparue. Il est alors utile de réappliquer le traitement temporel pour ne sélectionner que les composantes vraisemblables du signal. En effet, les basses fréquences temporelles ne sauraient provenir d'une éjection, dont la durée est brève, quelle que soit l'échelle spatiale à laquelle elle est considérée. Lors d'une étape 127, les moyens de détection 18 détectent des signaux d'alarmes dont une magnitude est supérieure à un seuil. Ces zones alarmes correspondant à une éjection possible de l'objet. On fixe la valeur du seuil en fonction du taux de fausses alarmes que l'on veut obtenir. La figure 13A représente schématiquement une partie d'un domaine spatiotemporel détecté avec une éjection 140 d'un missile et un point fixe 141 sans appliquer un procédé selon l'invention. La figure 13B représente schématiquement une partie d'un domaine 15 spatio-temporel détecté avec une éjection 140 d'un missile et deux points 142 et 143 mobiles sans appliquer un procédé selon l'invention. Dans les cas des figures 13A et 13B, le dispositif de détection ne détecte pas l'éjection 140 de l'objet, à causes de la détection également des autres points 141, 142 et 143. Il y a indifférenciation entre les phénomènes. La figure 14 représente schématiquement le domaine de la figure 13B en appliquant un procédé selon l'invention. Pour l'application du procédé, le filtre ponctuel spatial de l'étape 121 est le laplacien de la tache optique correspondant à un pixel, cette tache étant une gaussienne. L'élément structurant du filtre morphologique non linéaire de l'étape 123 est d'une durée de cinq échantillons. Le filtre passe-bas de l'étape 124 est une fonction constante égale à 1 sur sept pixels. Le paramètre Nmax est fixé à 2, il y a donc trois filtres morphologiques non linéaires temporels. Dans ce cas, l'éjection de l'objet est bien détectée par le dispositif, 30 puisque les autres phénomènes ne sont plus eux détectés. On a donc un taux de fausses alarmes faible. REVENDICATIONS 1. Procédé de détection d'un objet (15) situé dans un champ optique (16) d'un système optronique (11), caractérisé en ce qu'il comporte les étapes 5 selon lesquelles: - des moyens de traitement (12) déterminent un spectre spatio-temporel (41) à partir d'un signal issu du système optronique (11) ; - un filtre spatio-temporel (17) crée des zones de recherche (71) dans le spectre spatio-temporel (41) pour éliminer les basses fréquences spatiales 10 et les basses fréquences temporelles et; - des moyens de détection (18) détectent dans les zones de recherche (71) des zones d'alarmes dont une magnitude de spectre est supérieure à un seuil, ces zones d'alarmes correspondant à une détection possible de l'objet (15). 2. Procédé selon la revendication précédente, dans lequel le filtre (17) spatio-temporel crée des zones de recherche (71) dans le spectre spatiotemporel (41) en effectuant les étapes suivantes: - un filtrage (123) d'un signal issu du système optronique (11) par un filtre 20 morphologique non linéaire; - un filtrage (124) du signal issu du filtre morphologique non linéaire par un filtre spatial passe-bas; - un filtrage (126) du signal issu du filtre spatial passe-bas par un filtre morphologique non linéaire. 3. Procédé selon la revendication précédente, dans lequel l'étape de filtrage (123) du signal issu du système optronique (11) par le filtre morphologique non linéaire et l'étape de filtrage (124) du signal issu du filtre morphologique non linéaire par le filtre spatial passe-bas sont itérées une pluralité de fois. 4. Procédé selon la revendication précédente, dans lequel le nombre d'itérations est fonction de la sélectivité en fréquence du filtre (17). 1 5 5. Procédé selon l'une des revendications 2 à 4, dans lequel le filtre effectue un filtrage (121) de détection spatiale ponctuelle avant le filtrage par le filtre morphologique non linéaire. 6. Procédé selon l'une des revendications 2 à 5, dans lequel le filtre effectue un filtrage (122) antidistorsion de repliement avant le filtrage par le filtre morphologique non linéaire. 7. Procédé selon l'une des revendications 2 à 6, dans lequel la 10 détermination du spectre spatio-temporel se fait pour chaque signal issu d'un pixel du système optronique (11). 8. Procédé selon l'une des revendications précédentes, dans lequel l'objet est un missile (15) sol-air à très courte portée éjecté à partir d'un dispositif 15 d'éjection (14). 9. Procédé selon l'une des revendications précédentes, dans lequel le système optronique (11) est monté sur un véhicule (10). 10. Procédé selon la revendication précédente, dans lequel un volume de données (0) composant le signal pour la détermination du spectre spatiotemporel est fonction de la vitesse (V) du véhicule dans un repère terrestre. 11. Procédé selon la revendication précédente, dans lequel le volume de 25 données (0) est proportionnel à la racine cubique de la vitesse (V113) du véhicule dans un repère terrestre. 12. Dispositif de détection d'un objet (15), comportant un système optronique (11), l'objet (15) étant situé dans un champ optique (16) du 30 système optronique, caractérisé en ce qu'il comporte: - des moyens de traitement (12) aptes à déterminer un spectre spatio-temporel (41) à partir d'un signal issu du système optronique (11) ; - un filtre spatiotemporel (17) apte à créer des zones de recherche (71) dans le spectre spatio-temporel (41) pour éliminer les basses fréquences spatiales et les basses fréquences temporelles et; - des moyens de détection (18) aptes à détecter dans les zones de recherche (71) des zones d'alarmes dont une magnitude de spectre est supérieure à un seuil, ces zones alarmes correspondant à une détection possible de l'objet (15).