FR 2955190 A1 20110715 FR 1000140 A 20100114 La présente invention appartient au domaine de l'exploitation d'images obtenues par un capteur optique embarqué à bord d'un véhicule survolant la terre. Plus particulièrement, la présente invention concerne un procédé de détection de cibles dans une image ou une séquence d'images représentant sensiblement une même scène à la surface d'un corps céleste tel que la Terre. Par « cible », on entend tout événement ponctuel qui apparaît dans une image. L'invention trouve une application particulièrement avantageuse, bien que nullement limitative, dans la détection de cibles du type missile balistique, d'après des images acquises par un capteur matriciel optique embarqué dans un satellite en orbite Terrestre. L'invention pourrait également être appliquée à la détection de départs de feux de forêts par exemple. Les images acquises par de tels capteurs se présentent généralement sous la forme de matrices de pixels, agencés en lignes et en colonnes, chaque pixel représentant une portion d'une scène observée. La détection de cibles dans de telles images, en particulier des missiles balistiques, pose un problème lié à la dimension desdites cibles. En effet, du fait notamment de l'éloignement du satellite par rapport à la Terre, les dimensions réelles des cibles sont généralement de dimensions inférieures à celles des portions représentées par les pixels des images. Les dimensions des cibles, considérées par rapport aux dimensions d'une image, sont généralement pixelliques (c'est-à-dire à l'échelle du pixel), voire souspixelliques (c'est-à-dire inférieures au pixel). Il est possible de réduire les dimensions desdites portions au moyen d'un capteur optique haut de gamme. Toutefois, l'utilisation d'un capteur optique haut de gamme s'accompagne d'une augmentation de la complexité du système d'observation, et par conséquent d'une augmentation de son coût. De plus, la résolution de tels capteurs haut de gamme est elle-même limitée, et on comprend que d'autres moyens, en particulier des moyens de traitement des images acquises, doivent être mis en oeuvre pour améliorer la détection de cibles faiblement résolues, c'est-à-dire de dimensions pixelliques ou sous-pixelliques. Pour la détection de cibles dans des images, il est connu, en particulier de la publication scientifique « Detecting small moving objects using 2 temporal hypothesis testing », des auteurs A. Tzannes et D.H. Brooks (IEEE Transactions on Aerospace and Electronics Systems, vol. 38, n° 2, Avril 2002), de supprimer préalablement le fond des images. Le fond d'une image correspond à la partie essentiellement statique de la scène observée, c'est-à- dire dont les variations sont lentes par rapport au temps estimé d'apparition d'une cible. Suivant cet état de l'art, la détection s'effectue en comparant, pour chaque pixel, la radiométrie mesurée, à laquelle le fond a été préalablement retranché, à un seuil prédéfini qui est le même pour tous les pixels de l'image. Toutefois, un inconvénient d'une telle approche réside dans le fait que les performances de la détection, notamment exprimées sous la forme d'une probabilité de fausse alarme (c'est-à-dire la probabilité de détection d'une cible en l'absence de cible) dépendent du rapport signal sur bruit dans l'image. Le rapport signal sur bruit peut varier d'un pixel à l'autre, en fonction de la scène observée représentée par le fond de l'image. Par exemple, le niveau de signal mesuré pour un nuage sera généralement supérieur au niveau d'un signal mesuré en l'absence de nuage. De plus, le niveau de bruit dépend également de la scène observée. En particulier, le bruit photonique, introduit par le capteur optique, peut être modélisé par une loi Normale dont la variance serait égale au niveau du signal reçu de la portion de scène observée. On comprend donc que le niveau de bruit n'est pas nécessairement le même d'un pixel à l'autre. On comprend que la probabilité de fausse alarme ne sera pas la même suivant, par exemple, qu'on détecte une cible sur fond nuageux (c'est-à- dire avec un nuage entre la cible et la Terre) ou sur fond de Terre (c'est-à-dire sans nuage entre la cible et la Terre). Cet inconvénient est généralisable à tous les procédés connus de détection de cibles. Ceci est d'autant plus gênant en limite de nuage, puisque le rapport signal sur bruit, et donc les performances de détection, vont varier d'une image à l'autre. En particulier, des mouvements parasites de la ligne de visée du capteur au cours de l'acquisition d'images induisent une erreur, dite « de registration », et les scènes observées d'une image à l'autre ne sont pas exactement les mêmes. Un pixel en limite de nuage correspondra, d'une image à l'autre, à un fond nuageux ou à un fond de Terre, et on comprend que le passage d'un signal correspondant à un fond de Terre à un signal correspondant à un fond nuageux pourra entraîner une fausse alarme. La présente invention vise à proposer un procédé de détection de cibles dont les performances seront améliorées par rapport aux procédés connus. En particulier, la présente invention vise à proposer un procédé de détection dont les performances, notamment exprimées sous la forme de probabilité de fausse alarme, sont sensiblement les mêmes quel que soit le niveau de signal détecté dans une image. La présente invention vise également à proposer un procédé de détection dont certains modes de mise en oeuvre soient robustes en la présence de décalages résiduels entre images acquises successivement, décalages qui seraient dus notamment à des erreurs de registration induites par des mouvements parasites de la ligne de visée du capteur optique. La présente invention a pour objet un procédé de détection d'une cible dans au moins une image, dite « image de détection », acquise par un capteur optique embarqué dans un véhicule survolant un corps céleste, ladite image de détection étant constituée d'une pluralité de pixels, caractérisé en ce qu'il comporte les étapes : - d'estimation d'un fond de l'image de détection en chaque pixel à évaluer, - d'évaluation, en chaque pixel à évaluer, d'un test de détection par rapport à au moins un seuil dont la valeur en ce pixel est déterminée en fonction du fond estimé en ce pixel. De préférence, pour chaque pixel à évaluer, on estime des caractéristiques de l'erreur d'estimation du fond en ce pixel, puis la valeur de l'au moins un seuil en ce pixel est déterminée en fonction du fond estimé en ce pixel et des caractéristiques estimées de l'erreur d'estimation du fond en ce pixel. Suivant un premier exemple, les caractéristiques estimées de l'erreur d'estimation du fond en un pixel à évaluer sont au moins l'écart-type de l'erreur d'estimation en ce pixel. De préférence, pour chaque pixel à évaluer, la valeur de l'au moins un 4 seuil en ce pixel est exprimée dans un premier temps en fonction du fond estimé en ce pixel, puis la valeur de l'au moins un seuil en ce pixel est déterminée en remplaçant le fond estimé en ce pixel par le fond estimé en ce pixel augmentée d'une quantité représentative de l'erreur d'estimation du fond en ce pixel, cette quantité étant par exemple proportionnelle à l'écart-type estimé de l'erreur d'estimation du fond en ce pixel, et la constante de proportionnalité étant avantageusement la même en chaque pixel à évaluer. Suivant un autre exemple, les caractéristiques estimées de l'erreur d'estimation du fond en un pixel à évaluer sont au moins une fonction représentative d'une densité de probabilité de l'erreur d'estimation du fond en ce pixel. De préférence, on utilise l'estimation du fond pour modéliser, en chaque pixel à évaluer, un bruit instationnaire dont les caractéristiques dépendent d'un niveau du signal mesuré en ce pixel. Le bruit instationnaire est par exemple un bruit photonique, modélisé de préférence par une variable aléatoire gaussienne. II est utilisé en particulier pour modéliser le bruit de mesure servant à l'estimation du fond et pour normaliser la fonction de test de détection. De préférence, pour chaque pixel à évaluer, la valeur de l'au moins un seuil en ce pixel est déterminée en fonction du fond estimé en ce pixel, des caractéristiques estimées de l'erreur d'estimation du fond en ce pixel, et de caractéristiques supposées du signal émis par la cible. Les caractéristiques supposées du signal émis par la cible sont par exemple un niveau minimal EMIN de la cible à son apparition. De préférence, la fonction de test, utilisée par le test de détection, s'exprime en chaque pixel sous la forme d'une différence de distances entre la mesure en ce pixel et les valeurs que cette mesure devrait avoir selon au moins deux hypothèses respectivement de présence et d'absence de cible dans l'au moins une image de détection. De préférence, l'étape d'estimation du fond met en oeuvre un filtre de Kalman, de préférence de type « sans odeur », prenant en compte un modèle d'évolution du fond de l'image de détection. De préférence, le modèle d'évolution du fond de l'au moins une image de détection comporte une composante d'évolution spatiale et une composante d'évolution temporelle. Avantageusement, les deux composantes sont supposées indépendantes et la composante d'évolution spatiale rend compte d'une variation due aux erreurs de registration des images. Dans un mode de réalisation particulier, on considère un groupe de pixels autour du pixel à évaluer, et on effectue des calculs statistiques sur ce groupe de pixels afin de déterminer lesdites deux composantes. De préférence, on estime le fond en fonction d'une image ou d'une séquence d'images, dites « images d'estimation ». Avantageusement, les images d'estimation sont des images acquises préalablement à l'image de détection, séparées de ladite image de détection par un délai non nul, égal ou supérieur au temps supposé d'apparition d'une cible. De préférence, pour chaque pixel à évaluer, la valeur de l'au moins un seuil en ce pixel est déterminée de sorte que la probabilité de fausse alarme dudit test de détection en ce pixel est théoriquement égale ou inférieure à une valeur spécifiée, la valeur spécifiée étant avantageusement la même pour chaque pixel à évaluer. Par exemple, pour chaque pixel à évaluer, on exprime la probabilité de fausse alarme du test de détection en ce pixel en fonction de la valeur de l'au moins un seuil en ce pixel, du fond estimé en ce pixel et, le cas échéant, des caractéristiques de l'erreur d'estimation du fond en ce pixel et de caractéristiques supposées du signal émis par la cible, puis la valeur dudit au moins un seuil en ce pixel est déterminée de sorte que la probabilité de fausse alarme en ce pixel soit égale ou inférieure à la valeur spécifiée. Dans une réalisation particulière, pour chaque pixel à évaluer, on exprime la probabilité de fausse alarme du test de détection en ce pixel en fonction de la valeur d'au moins un seuil en ce pixel, du fond estimé en ce pixel, et de ladite fonction représentative d'une densité de probabilité de l'erreur d'estimation du fond, puis la valeur dudit au moins un seuil en ce pixel est déterminée en fonction dudit fond estimé en ce pixel et d'un tirage statistique d'un ensemble de réalisations de l'erreur d'estimation du fond en ce pixel suivant ladite densité de probabilité, de sorte que la probabilité de fausse alarme en ce pixel soit égale ou inférieure à une valeur spécifiée. De préférence, on détermine au moins deux seuils pour chaque pixel à évaluer, et le test de détection est un test séquentiel. De préférence, parmi les au moins deux seuils pour chaque pixel à évaluer, la valeur d'un premier seuil en ce pixel est déterminée de sorte que la probabilité de fausse alarme en ce pixel est théoriquement égale ou inférieure à une première valeur spécifiée, et la valeur d'un second seuil en ce pixel est déterminée de sorte que la probabilité de détection en ce pixel est théoriquement égale ou supérieure à une seconde valeur spécifiée. De préférence, le test de détection utilise une pluralité d'images de détection. Par exemple, on applique le procédé en considérant une image reconstituée dont la valeur en chaque pixel à évaluer est égale à la moyenne des valeurs des pixels correspondants dans la pluralité d'images de détection (après registration desdites images de détection). De préférence, le procédé comporte deux phases d'exécution : - une première phase, au cours de laquelle l'étape d'estimation du fond et l'étape d'évaluation du test de détection par rapport à des seuils déterminés en fonction du fond estimé sont exécutées tant qu'une cible n'a pas été détectée dans une image de détection, une seconde phase, exécutée lorsqu'une cible a été détectée dans une image de détection, au cours de laquelle l'étape d'estimation du fond n'est pas exécutée. La description suivante de modes de mise en oeuvre de l'invention est faite en se référant aux figures, dans lesquelles des références identiques désignent des éléments identiques ou analogues, qui représentent de manière non limitative : Figure 1 : un diagramme représentant schématiquement les étapes principales d'un procédé de détection de cibles selon l'invention, Figures 2, 3 et 4 : des diagrammes représentant schématiquement des modes particuliers de mise en oeuvre du procédé de détection, Figure 5 : un exemple d'évolution temporelle de la valeur d'un pixel d'une image, au cours de l'apparition d'une cible dans ledit pixel. La présente invention concerne un procédé de détection de cibles dans une image, dite « image de détection », acquise par un capteur embarqué 7 dans un véhicule survolant un corps céleste, de préférence la Terre. Dans le contexte de l'invention, une image est préférentiellement une image bidimensionnelle, constituée par exemple d'une matrice de pixels. Chaque pixel représente une portion d'une scène observée depuis le véhicule survolant la Terre. Une scène observée correspond à une zone à la surface de la Terre, dont certaines parties sont occultées, notamment par des nuages. Le capteur est de préférence un capteur optique. Par exemple, il s'agit d'un capteur qui mesure le rayonnement de la scène observée dans le domaine des longueurs d'ondes infrarouges. Le véhicule, embarquant le capteur optique, peut être de tout type adapté à l'acquisition d'images successives d'une même scène observée. Par exemple, le véhicule est un satellite en orbite autour de la Terre, de préférence une orbite géostationnaire. Suivant un autre exemple non limitatif, le véhicule est un aéronef, tel qu'un hélicoptère. Par « cible », on entend de manière générale tout événement modifiant le fond de l'image. Le fond de l'image correspond à la partie essentiellement statique de la scène observée, c'est-à-dire dont les variations sont lentes par rapport au temps d'apparition d'une cible (qui dépend du type de cible considéré). La présente invention trouve une application avantageuse dans la détection de cibles du type engin aérien, par exemple un aéronef (avion, hélicoptère) ou un missile en phase propulsée. La mise en oeuvre du procédé de détection de cibles selon l'invention est par exemple assurée par un dispositif de calcul (microcontrôleur, ordinateur muni d'un microprocesseur, etc.). Le dispositif de calcul comporte des moyens de mémorisation (disque dur magnétique, mémoire flash, disque optique, etc.) dans lesquels est mémorisé un produit programme d'ordinateur, sous la forme d'un ensemble d'instructions de code de programme à exécuter pour mettre en oeuvre les étapes du procédé selon l'invention, ainsi que des images acquises par le capteur optique. Suivant certains modes de réalisation, le dispositif de calcul comporte également des circuits intégrés spécialisés, de type ASIC, FPGA, etc., adaptés à mettre en oeuvre tout ou partie des étapes d'un ou des procédés selon l'invention. Le dispositif de calcul est de préférence embarqué avec le capteur optique dans le véhicule survolant la Terre dans le cas d'un traitement en temps réel, adapté par exemple à la détection de missiles balistiques. Dans le cas d'un traitement en temps différé, le dispositif de calcul est embarqué avec le capteur optique dans le véhicule, ou bien il est déporté, c'est-à-dire localisé dans une station terrestre ou dans un second véhicule survolant la Terre. Définitions et notations Les définitions et notations utilisées dans la suite de la description sont données ci-après. On comprend que le choix d'une convention particulière n'est pas limitatif de l'invention qui pourrait être décrite de manière équivalente en adoptant d'autres conventions sans que l'invention en soit modifiée. Une image acquise à un instant t est désignée par y(t). L'image y(t) est constituée d'une pluralité de pixels, dont les valeurs sont désignées par y(p,t), où p désigne l'indice d'un pixel de l'image. Par exemple, p est compris entre 1 et NLxNc, où NL désigne le nombre de lignes de pixels de l'image y(t) et Nc désigne le nombre de colonnes de ladite image. Un pixel y(p,t) correspond principalement à un signal x(p,t) du rayonnement mesuré de la portion de la scène observée, perturbé par un bruit b(p,t) : y(p,t) = x(p,t) + b(p,t) (1) En l'absence de cible dans le pixel d'indice p, le signal x(p,t) est égal au fond de la scène observée, désigné par p(p,t). Si une cible est présente, le signal x(p,t) peut être exprimé comme suit : x(p,t) = p(p,t) + c(p,t) (2) où c(p,t) est le signal mesuré correspondant au rayonnement de la cible. Le bruit b(p,t) a diverses origines. En particulier, il peut être modélisé comme étant la somme d'un bruit bE(p,t) d'origine électronique, et un bruit bp(p,t) d'origine photonique. Le bruit électronique bE(p,t) représente un bruit de mesure du capteur optique. Le bruit électronique est généralement modélisé par une variable aléatoire BE(p,t) gaussienne, qui suit une loi de probabilités Normale N(0,a2E), da valeur moyenne nulle et d'écart-type aE indépendant de la scène observée. Le bruit photonique bp(p,t) tient compte du caractère aléatoire de la mesure du nombre de photons. Le bruit photonique est généralement modélisé 9 par une variable aléatoire Bp(p,t) qui suit une loi de probabilités de Poisson. Etant donné le nombre important de photons reçus, et en application de la loi des grands nombres, cette loi de probabilité est parfois approchée par une loi Normale N(0,x(p,t)) de valeur moyenne nulle et de variance x(p,t). En modélisant le bruit photonique par une variable aléatoire gaussienne, et en considérant les réalisations du bruit électronique et du bruit photonique comme indépendantes, le bruit b(p,t) peut être modélisé par une variable aléatoire B(p,t) gaussienne : B(p,t) ù N(O,o2E + x(p,t)) (3) De manière plus générale, on peut considérer d'autres sources d'erreurs, par exemple liées aux erreurs résiduelles de calibration de biais et de gain, désignées respectivement par aB(p) et par aG(p), qui varient peu au cours du temps, mais qui peuvent être différentes d'un pixel à un autre. Dans ce cas, l'expression de la valeur y(p,t) du pixel d'indice p peut s'exprimer sous la forme En pratique, les erreurs de calibration aB(p) et aG(p) peuvent être estimées en tout pixel en mettant en oeuvre des procédés connus de l'homme de l'art. A des fins de clarté de la description, mais à titre nullement limitatif de l'invention, on considère dans la suite de la description que aB(p) est égal à zéro et que aG(p) est égal à un. Les équations données ci-après peuvent être généralisées à des valeurs quelconques de aB(p) et de aG(p), en partant des équations (4) et (5). Procédé de détection de cibles La figure 1 représente schématiquement le procédé de détection de cibles dans une image, dite « image de détection ». Le procédé de détection de cibles comporte principalement une étape 10 d'estimation du fond de l'image de détection, une étape 20 de détermination de seuil pour chaque pixel à évaluer, et une étape 30 d'évaluation d'un test de détection par rapport aux seuils déterminés, en chaque pixel à évaluer de l'image de détection. suivante : y(p,t) = ac(p).x(p,t) + b(p,t) (4) où le bruit b(p,t) peut être modélisé par une variable aléatoire gaussienne : (5) B(p,t) ù N(aB(p),a2E + ap.x(p,t)) Par « à évaluer », on entend que la détection de cibles peut être restreinte à un sous-ensemble de pixels d'une image de détection, ou qu'au contraire la détection peut être effectuée sur toute l'image de détection. Par exemple, dans le cas d'une détection de missile en phase propulsée, on peut limiter l'évaluation du test de détection à un sous-ensemble de pixels représentant une zone autour d'une base de tir de missiles. Il est à noter que l'ordre d'exécution des différentes étapes du procédé de détection n'est pas limitatif de l'invention, et qu'un ordre particulier ne constitue qu'une variante d'implémentation du procédé de détection. On note toutefois que, pour un pixel donné, le fond dudit pixel doit être estimé (étape 10) avant de déterminer la valeur du seuil associé (étape 20), seuil qui doit être déterminé avant d'évaluer le test de détection (étape 30) pour ledit pixel donné. Lorsque l'on considère tous les pixels à évaluer, il est possible par exemple d'estimer le fond de l'image de détection pour tous les pixels à évaluer, de déterminer tous les seuils associés, puis d'évaluer le test de détection pour tous les pixels à évaluer. Suivant un autre exemple, les étapes 10 d'estimation, 20 de détermination et 30 d'évaluation sont exécutées successivement pour chaque pixel à évaluer. Suivant encore un autre exemple, l'étape 10 d'estimation est exécutée pour tous les pixels à évaluer, puis les étapes 20 de détermination de seuil et 30 d'évaluation de test de détection sont exécutées successivement pour chaque pixel à évaluer. Il est à noter que le procédé de détection peut comporter d'autres étapes, comme par exemple une étape préalable, dite « de registration », au cours de laquelle les images acquises sont recalées les unes par rapport aux autres, pour compenser des décalages introduits par les mouvements parasites de la ligne de visée du capteur optique au cours de l'acquisition. L'étape de registration, non représentée sur les figures, peut être mise en oeuvre selon tout procédé de registration connu de l'homme de l'art. Tel que représenté sur la figure 2, les étapes du procédé de détection sont de préférence exécutées de façon récurrente, au moins tant qu'on n'a pas détecté de cible, afin de surveiller si une cible apparaît et, le cas échéant, signaler une alarme. Cela revient à acquérir 40 et à traiter régulièrement une 11 nouvelle image de détection, par exemple de façon sensiblement périodique. II est à noter que les étapes du procédé peuvent, lorsqu'elles sont récurrentes, avoir des périodicités différentes. En particulier, l'évaluation du test de détection (étape 30) se fait de préférence à chaque nouvelle image de détection, tandis que l'estimation du fond (étape 10) peut se faire à chaque nouvelle image de détection, ou toutes les deux, trois ou plus images de détection, suivant la vitesse de variation de la scène observée. a) Etape 10 d'estimation du fond de l'image de détection Tel que décrit précédemment, le fond de l'image de détection correspond à la partie de la scène observée dont les variations sont lentes par rapport au temps estimé d'apparition d'une cible. L'étape 10 vise à calculer une estimation pe(p,t) du fond p(p,t) pour chaque pixel à évaluer. Dans un mode préféré de mise en oeuvre de l'étape 10, le fond de l'image de détection est estimé d'après une image, ou une séquence d'images, dites « images d'estimation », acquises à différents instants dans le temps. Les images d'estimation représentent sensiblement la même scène observée que l'image de détection, et sont acquises par un capteur se trouvant dans le plan focal du capteur optique ayant servi à acquérir l'image de détection. De préférence, lesdites images d'estimation sont acquises par le même capteur optique que celui ayant servi à acquérir l'image de détection. Rien n'exclut que l'image de détection soit également une image d'estimation, c'est-à-dire que l'image de détection soit également utilisée au cours de l'étape 10 pour estimer son propre fond. Dans un mode préféré de mise en oeuvre du procédé de détection, 25 illustré sur la figure 3, les images d'estimation sont des images acquises antérieurement à l'acquisition de l'image de détection courante. De préférence, les images d'estimation sont séparées dans le temps de l'image de détection courante par un délai non nul, égal ou supérieur au temps d'apparition i5 d'une cible. 30 En pratique, l'apparition de la cible dans un ou plusieurs pixels d'une image est progressive, tel qu'illustré par la figure 5. Le temps d'apparition b d'une cible dépend du type de cible, et correspond au temps nécessaire pour que le signal c(p,t) soit égal ou supérieur à un niveau minimal, désigné par sM,N, prédéfini par exemple en fonction d'une cible de référence correspondant théoriquement à un pire cas pour la détection. Dans la figure 5, on a représenté trois intervalles de temps, respectivement Il, 12 et 13, ainsi que l'évolution des valeurs y(p,t) (désigné par 5 « y ») et x(p,t) (désigné par « x ») en fonction du temps. Les images acquises au cours de l'intervalle Il ne comportent pas de cible, les images acquises au cours de l'intervalle 13 comportent une cible, et l'intervalle 12 correspond à un intervalle de transition, au cours duquel la cible apparaît, de durée égale à 6. 10 On comprend qu'il est avantageux d'estimer le fond de l'image de détection à partir d'images dans lesquelles il n'y a pas de cible présente, c'est-à-dire principalement en fonction d'images acquises au cours de l'intervalle Il. Si on considère des images d'estimation acquises à des instants antérieurs d'au moins la durée b à celui d'acquisition de l'image de détection 15 courante, et si on n'a pas détecté de cible dans des images de détection acquises antérieurement au cours de l'intervalle de durée b, on assure théoriquement que l'estimation du fond est obtenue à partir d'images d'estimation dans lesquelles il n'y a ni cible présente, ni cible en cours d'apparition. 20 En pratique, le temps d'apparition b est par exemple préalablement déterminé comme étant le temps d'apparition d'une cible de référence au-dessus d'un nuage, une cible de référence étant une cible dont on connaît les caractéristiques (ergol utilisé, trajectoire en phase propulsée, etc.). Le temps d'apparition b est par exemple de l'ordre de quelques centaines de 25 millisecondes à quelques secondes. De préférence, tel qu'illustré par la figure 3, les images d'estimation sont, au moins en partie, des images de détection antérieures, c'est-à-dire des images pour lesquelles le test de détection a déjà été évalué. Par exemple, tel que représenté sur la figure 3, le procédé comporte 30 une étape 50 au cours de laquelle on retarde les images de détection avant de les utiliser pour estimer le fond à utiliser pour des images détection ultérieures. Les images de détection ainsi retardées sont par exemple stockées dans les moyens de mémorisation susmentionnés. La figure 4 représente un mode particulier de mise en oeuvre du procédé de détection, dans lequel la détection s'effectue principalement en deux phases. Au cours de la première phase, les étapes principales du procédé de détection, l'étape 10 d'estimation, l'étape 20 de détermination de seuils et l'étape 30 d'évaluation du test de détection, sont exécutées tant qu'une cible n'a pas été détectée. Lorsqu'une cible a été détectée dans une image de détection, la seconde phase consiste principalement à assurer le suivi d'une cible détectée. La différence principale entre la première et la seconde phase réside dans le fait que l'étape 10 d'estimation du fond n'est plus exécutée, afin d'éviter d'estimer le fond sur des images susceptibles de comporter une cible. Par exemple, on utilise la même estimation du fond que celle utilisée lorsque la cible a été détectée. De préférence, et tel que représenté sur la figure 5, l'étape 20 de détermination de seuils n'est plus exécutée au cours de la seconde phase. Par exemple, on utilise les mêmes les seuils que ceux utilisés lorsque la cible a été détectée. Si au cours de la seconde phase, la cible n'est plus détectée, le procédé de détection s'arrête ou rebascule sur la première phase. Suivant un mode préféré de mise en oeuvre du procédé de détection, non illustré par les figures, on estime également, au cours de l'étape 10 d'estimation, des caractéristiques statistiques de l'erreur d'estimation du fond. Suivant un premier exemple, on estime l'écart-type de l'erreur en chaque pixel, désigné par QEST(p,t). Suivant un autre exemple non limitatif, on estime une fonction représentative de la densité de probabilité de l'erreur d'estimation du fond. L'étape 10 d'estimation s'appuie de préférence sur un modèle d'évolution du fond de l'image, de préférence un modèle Markovien, par 30 exemple d'ordre 1, en chaque pixel d'indice p : x(p,t+At) = x(p,t) + w(p,t) (6) expression dans laquelle : - At est la période d'acquisition entre deux images consécutives d'une séquence d'images ; on considère, à titre nullement limitatif de l'invention, que la période d'acquisition est sensiblement constante, - w(p,t) est un terme d'évolution, modélisé par une variable aléatoire gaussienne de moyenne nulle et d'écart-type aEVOL(p,t)• Il est connu, en particulier de la publication scientifique précitée « Detecting small moving objects using temporal hypothesis testing », des auteurs A. Tzannes et D.H. Brooks (IEEE Transactions on Aerospace and Electronics Systems, vol. 38, n° 2, Avril 2002), qu'un tel modèle d'évolution temporelle est adapté notamment aux images de nuages dans le domaine infrarouge. De manière plus générale, d'autres modèles d'évolution peuvent être considérés, et le choix d'un modèle particulier ne constitue qu'une variante d'implémentation de l'invention. A partir du modèle d'évolution choisi, on utilise un estimateur qui fournit à la fois une estimation pe(p,t) du fond p(p,t) de l'image de détection, et une estimation des caractéristiques de l'erreur d'estimation, par exemple une estimation aeEsT(p,t) de l'écart-type de l'erreur d'estimation aEST(p,t) ou une fonction représentative de la densité de probabilité de l'erreur d'estimation. Un estimateur du type filtrage de Kalman est particulièrement adapté au cas d'un modèle Markovien d'évolution temporelle du fond. Un tel estimateur peut prendre en compte plusieurs images d'estimation et fournir une estimation du fond et, le cas échéant, une estimation des caractéristiques de l'erreur d'estimation. De préférence, l'estimation est faite au moyen d'un filtre de Kalman du type « sans odeur » (« Unscented Kalman Filter » ou UKF dans la littérature anglo-saxonne). Un tel filtre de Kalman sans odeur est connu par exemple du livre « The Unscented Kalman Filter », des auteurs E. Wan et R.V. der Merwe, publié aux éditions Wiley. De manière plus générale, l'étape 10 d'estimation peut mettre en oeuvre tout procédé d'estimation du fond d'images et, le cas échéant, d'estimation du bruit d'estimation, connu de l'homme de l'art. b) Etapes 20 de détermination de seuils et 30 d'évaluation du test de détection Selon l'invention, on définit une pluralité de seuils q(p,t), au moins un seuil par pixel, dont la valeur est déterminée en fonction de l'estimation pe(p,t) du fond p(p,t) de l'image de détection. Au cours de l'étape 30 d'évaluation, on évalue un test de détection en chaque pixel à évaluer, test de détection qui dépend du fond estimé pour le pixel évalué. Le test de détection est de préférence un test statistique d'hypothèses de type Neyman-Pearson, ou son extension sous forme de rapport de vraisemblance généralisée, dans lequel on considère au moins deux hypothèses Ho(p,t) et Ho(p,t) : - Ho(p,t) : il n'y a pas de cible à l'instant t dans le pixel d'indice p, - Ho(p,t) : une cible est présente à l'instant t dans le pixel d'indice p. Un test de détection du type test statistique d'hypothèses est généralement basé sur le choix d'une fonction de test, appelée fonction T dans la suite de l'exposé, construite à partir des mesures effectuées. L'évaluation du test de détection revient à calculer la valeur de la fonction de test T à partir des mesures considérées, et à comparer ladite valeur à un seuil. Suivant un exemple non limitatif, l'hypothèse Ho(p,t) est vérifiée si la valeur de la fonction de test T est inférieure au seuil, et l'hypothèse Ho(p,t) est vérifiée si la valeur de la fonction de test T est supérieure au seuil. Le calcul du seuil dépend avantageusement du test de détection. Plus particulièrement, il dépend de la fonction T de test choisie. D'un point de vue théorique, il est possible de définir un test de détection optimal, dit « test du rapport de vraisemblance généralisée », à partir d'un modèle du signal y(p,t) mesuré. Ce type de test est supposé connu de l'homme de l'art, et sa formulation théorique n'est pas rappelée. Toutefois, la mise en oeuvre du test du rapport de vraisemblance généralisée est difficile en pratique, du fait d'une complexité analytique élevée, qui est incompatible avec les capacités de calcul actuelles, en particulier les capacités des dispositifs de calcul actuellement embarqués dans des satellites, et/ou avec des contraintes d'évaluation en temps réel ou quasi-réel (en particulier dans le cas d'une détection de missiles balistiques). Avantageusement, on considère certaines approximations du test optimal de détection visant à réduire la complexité analytique de son évaluation. On décrit ci-après un exemple non limitatif de test de détection dont la complexité analytique est adaptée aux capacités des dispositifs de calcul actuellement embarqués dans des satellites. De préférence, le calcul des seuils q(p,t) tient compte de la présence du bruit photonique bp(p,t) dans le signal y(p,t), au moyen d'un modèle de bruit photonique. De manière plus générale, le calcul des seuils q(p,t) peut tenir compte de tout autre type de bruit instationnaire. Dans la suite de la description, on se place, à titre nullement limitatif, dans le cas où un modèle de bruit photonique est effectivement pris en compte. Ledit test de détection est obtenu en considérant tout d'abord une approximation du modèle théorique du bruit photonique (variable aléatoire qui suit une loi de Poisson), qu'on modélise par une variable aléatoire gaussienne N(0,x(p,t)), tel que décrit précédemment. Compte tenu des modèles choisis pour le bruit électronique et le bruit photonique, la valeur y(p,t) du pixel d'indice p peut être considérée comme la réalisation d'une variable aléatoire Y(p,t) gaussienne : Y(p,t) N(x(p,t), Q2E + x(p,t)) (7) où x(p,t) est égal à p(p,t) en l'absence de cible, à (p(p,t) + c(p,t)) en présence d'une cible. Le test de détection proposé est ensuite obtenu en approchant le test optimal formulé d'après ce modèle de y(p,t). Le test de détection proposé dépend principalement d'un niveau minimal SMIN du signal c(p,t) qui serait mesuré pour une cible, de l'écart-type o.E, 25 et du fond p(p,t) de l'image de détection. A titre nullement limitatif, on considère que le niveau minimal CMIN et l'écart-type QE sont connus a priori : le niveau minimal EMIN est par exemple prédéterminé en fonction d'une cible de référence correspondant théoriquement à un pire cas pour la détection, et l'écart-type aE est par 30 exemple prédéterminé par calibration du capteur optique. En pratique, le fond p(p,t) de l'image de détection est estimé au cours de l'étape 10 d'estimation, qui fournit une estimation pe(p,t) dudit fond, utilisée dans les différents calculs décrits ci-après. Dans un mode préféré de mise en oeuvre, la fonction T de test optimale au sens précédent peut être approximée comme suit : - si y(p,t) û pe(p,t) ? EMIN T(y(p,t),pe(p,t)) _ (Y(p,2)ûpe(p,t))2 QE +pe(p,t) (8) MIN La valeur de la fonction de test T en un pixel donné est en fait la différence des distances de la mesure en ce pixel aux valeurs que cette mesure devrait avoir suivant respectivement les au moins deux hypothèses du test, c'est-à-dire la valeur du fond en ce pixel s'il n'y a pas de cible, et la valeur du fond en ce pixel auquel s'ajoute le signal EMIN de la cible s'il y a une cible. On peut en effet écrire cette fonction de test T comme la différence : T(y(p, t), pe(p, t)) = Do (y(p, t), pe(p, t)) û D, (y(p, t), pe(p, t)) (10) où la distance Do pour l'hypothèse Ho(p,t) d'absence de cible est donnée par : Do (Y(p, t), pe(p, t)) _ (y(p,t) û pe(p, t))2 (11) QE +pe(p,t) et la distance DI pour l'hypothèse Ho(p,t) de présence de cible est donnée par : (y(p, t) û pe(p, t) - EMIN )2 si y(p, t) Dans cet exemple, la distance est avantageusement quadratique et pondérée par la valeur du signal de manière à tenir compte de la nature en partie photonique du bruit d'acquisition. On comprend que cet exemple n'est pas limitatif de l'invention, et qu'un autre type de distance peut être utilisé. Considérons un seuil q(p,t) tel que si T(y(p,t),pe(p,t)) est supérieur à ce seuil q(p,t), on conclut que l'hypothèse Ho(p,t) est la plus vraisemblable et dans le cas contraire, on conclut que c'est l'hypothèse Ho(p,t) qui est la plus vraisemblable. Dans le cadre de l'invention, on cherche principalement à maîtriser la - si y(p,t) û pe(p,t) A des fins de clarté de la description, les paramètres (p,t) sont temporairement omis. La probabilité de franchissement d'un seuil n par la fonction T, désignée par h(n,pe,x) peut alors s'exprimer : h(q,pe,x) = prob(T(Y,pe) q) (13) expression dans laquelle Y est distribué tel que décrit à l'expression (7). En remplaçant la fonction T par les expressions (8) et (9), la probabilité de franchissement d'un seuil n peut s'exprimer de la façon suivante : h(q,pe,x) = prob(Y'- gi(q,pe)) + prob(g2(q,pe) _ q(p, t) M(p, t))] (23) La fonction k peut s'exprimer sous la forme : k(ri(p,t),pe(p,t),QeEST (p,t)) = (pù ue)2 f h(q(p, t), pe(p, t), p) 1 QeEST (p, t) 2Tr exp _ 2QeÉST (p, t) Jdii (24) L'intégrale de l'expression (24) n'est pas calculable analytiquement, mais peut être approchée en utilisant toute technique d'intégration numérique connue de l'homme de l'art. Suivant un exemple non limitatif, cette intégrale est approchée en mettant en oeuvre une technique classique de Monte Carlo par pondération de vraisemblance. Pour des valeurs de pe(p,t) et QeEsT(p,t) données, on génère un nombre Ne d'échantillons pi (typiquement de l'ordre de 10000 échantillons) suivant la loi Normale N(pe(p,t), f2-Qe2EsT(p,t)) où f est un facteur d'amplification de creEsT(p,t) permettant, de manière connue, d'améliorer l'approximation de l'intégrale. Ladite intégrale peut ensuite être approximée comme suit : Ne k(q(p,t),pe(p,t),QeEST(p,t))Ewi •h(n(p,t),pe(p,t),) (25) i=1 expression dans laquelle : w ù ex J 1 1 (pe(p, t) -Pi )2 Ne p 2 _ 2• IQe/p,t (26) EST ) Comme précédemment, on peut déterminer des valeurs de la fonction k(q(p,t),pe(p,t),aeEST(p,t)) pour différentes valeurs du seuil q(p,t), et sélectionner le seuil ri(p,t) pour lequel on obtient une probabilité de fausse alarme sensiblement égale à PFAO. Lesdites valeurs de la fonction k(q(p,t),pe(p,t),QeEsT(p,t)) sont de préférence préalablement déterminées, pour différentes valeurs de pe(p,t) et de QeEsT(p,t), et mémorisées dans des moyens de mémorisation. Pour estimer des probabilités de fausse alarme autour de 1%, on considère de préférence une valeur de f comprise entre deux et dix. On comprend donc qu'en estimant des caractéristiques de l'erreur d'estimation et en déterminant les seuils en fonction des caractéristiques 22 estimées, il sera possible de mieux maîtriser la probabilité de fausse alarme. Par exemple, les seuils sont déterminés de sorte que la probabilité de fausse alarme est théoriquement la même en chaque pixel à évaluer, ou de sorte que la probabilité de fausse alarme est théoriquement inférieure à une valeur prédéfinie en chaque pixel à évaluer. Extension du procédé au cas de plusieurs images de détection De manière plus générale, on peut utiliser, au cours de l'étape 30 d'évaluation du test de détection, une pluralité de mesures, c'est-à-dire une pluralité d'images de détection, que l'on aura pris soin de recaler les unes par rapport aux autres, pour compenser des décalages introduits par les mouvements parasites de la ligne de visée du capteur optique au cours de l'acquisition. Dans la suite du paragraphe, on suppose que cette étape de registration est effectuée au préalable. Par exemple, on peut considérer, de manière non limitative, M (M > 1) images de détection acquises successivement avec une période At, formant un ensemble désigné par yM(t), comprenant les images y(t+mot), 0 m M-1. On désigne également par yM(p,t) l'ensemble des pixels y(p,t+mAt), 0 n) (33) La probabilité hM(ri,pe,xm) de franchissement d'un seuil peut également s'exprimer de la façon suivante : hM(q,pe,xm) = prob(ym gMi(n,pe)) + prob(gM2(ri,pe)) -` ym La probabilité d'avoir ym supérieur à une valeur p est donnée par : prob(ym'- p) = Q(p, xm, 62MOY) (37) En combinant les équations (30) et (33), on obtient : hM(n,pe,xm) = Q(gM,(q,pe), xm, Q2MOY) + Q(gM2(n,pe), xm, Q2MOY) (38) ù Q(pe + EMIN, xm, Q2MOY) Comme précédemment, on peut considérer que la probabilité de fausse alarme PFA(P,t) en un pixel p et un instant donné t est égale à PFA(p,t) = hM(n(p,t),pe(p,t),pm(p,t)), c'est-à-dire la probabilité que le test d'hypothèse déclenche l'alarme (la fonction Tm de test dépassant le seuil) alors qu'il n'y a pas de cible en ce pixel (xm(p,t) = pm(p,t) tel que définit à l'expression (31)). On peut par exemple déterminer la valeur du seuil q(p,t) pour laquelle la fonction hM(n(p,t),pe(p,t),pm(p,t)) vaut théoriquement PFAO. Tel que décrit précedemment, pm(p,t) est par exemple remplacé par l'estimation pe(p,t), ou par (pe(p,t) + k.oeEST(p,t)), ou par un calcul d'intégrale prenant en compte une densité de probabilité de l'erreur d'estimation du fond. Modèle d'évolution spatio-temporelle du fond Dans un mode préféré de mise en oeuvre, on raffine l'étape 10 d'estimation du fond en décomposant le terme d'évolution w(p,t) défini dans l'expression (6) en au moins deux composantes indépendantes : w(p,t) = WTEMP(P,t) + WSPATIAL(P,t) (39) WTEMP(P,t) est une composante d'évolution temporelle, par exemple distribuée suivant une loi Normale N(0,Q2TEMP(P,t)), et WSPATIAL(P,t) est une composante d'évolution spatiale, par exemple distribuée suivant une loi Normale N(0,Q-2SPATIAL(p,t)), de sorte que : G EVOL(P,t) - G TEMP(P,t)) + J SPATIAL(P,t) (40) La composante d'évolution temporelle modélise l'évolution temporelle de la scène observée. La composante d'évolution spatiale modélise l'erreur de registration entre images successives. Suivant un exemple non limitatif, les écart-types QTEMP(p,t) et QSPATIAL(p,t) sont estimés en considérant, pour un pixel donné, un groupe de neuf pixels centré sur ledit pixel donné. De manière plus générale, on peut considérer d'autres configurations de groupes de pixels, en particulier des groupes comportant un nombre plus important de pixels. La coordonnée pixel p est remplacée temporairement par des coordonnées (i,j) (ligne, colonne) conventionnelles. L'écart-type QTEMP(p,t) de l'évolution temporelle peut par exemple s'exprimer sous la forme : QTzEMP (I, i, t)_ - K p=1 q=1 2 1 1 1 P=1 q=1 Kù1E g EEpe(i+p,j+q,t) 9 EEpe(I+p,j+q,tùkAt) (41) k=1 p=-1q=-1 p=-lq=-1 qui représente l'évolution estimée du groupe de neuf pixels sur une fenêtre temporelle de K images (typiquement cinq à dix images). L'écart-type QSPATIAL(p,t) de l'évolution spatiale peut par exemple s'exprimer sous la forme : P=1 q=1 pe(i, j, t) ù 9 E E pe(i + p, j + q, t) p=-1q=-1 Q SPATIAL (I, J, t) - G REG (42) OÙ °REG désigne l'écart-type de l'erreur de registration, qui peut être estimé en mettant un procédé connu de l'homme de l'art, par exemple en fonction de 10 mesures de l'attitude du véhicule embarquant le capteur optique. Généralisation au test séquentiel De manière plus générale, dans des modes de mise en oeuvre du procédé de détection basés en particulier sur le traitement de plusieurs images de détection (cas « multi-mesures »), on peut considérer, pour chaque pixel à 15 évaluer, deux seuils n1(p,t) et q2(p,t) (avec g1(p,t) Emin). De préférence, la valeur du seuil r12(p,t) est déterminée de sorte que la probabilité de détection en ce seuil est égale ou supérieure à une valeur spécifiée. On comprend que, du fait de la prise en compte du fond dans le test de détection et dans les seuils associés à chaque pixel à évaluer, l'invention permet de mieux maîtriser les performances de la détection. En particulier, l'invention permet de maîtriser la probabilité de fausse alarme en chaque pixel à évaluer, notamment du fait d'une prise en compte, d'une part des caractéristiques de l'erreur d'estimation du fond, et d'autre part du bruit photonique présent dans les images acquises, dont les caractéristiques varient d'un pixel à l'autre. Par ailleurs, le procédé de détection améliore les performances de la détection, en particulier face à des décalages résiduels entre des images acquises à des instants différents, décalages qui sont dus notamment à des mouvements parasites de la ligne de visée du capteur optique au cours de l'acquisition. Les performances en présence de décalages résiduels sont notamment améliorées dans les modes de mise en oeuvre comportant l'une au moins des caractéristiques suivantes : - l'estimation du fond se fait sur une séquence d'images ; - l'estimation du fond se fait d'après un modèle d'évolution comportant une composante spatiale modélisant les erreurs de registration ; - le test de détection exploite une pluralité d'images de détection. Du fait que les performances de la détection sont améliorées, on comprend que l'on peut avantageusement obtenir un niveau de performances 25 donné avec un dispositif qui sera moins complexe et moins coûteux qu'un dispositif mettant en oeuvre un procédé de détection suivant l'état de l'art. REVENDICATIONS1 Procédé de détection d'une cible dans au moins une image, dite « image de détection », acquise par un capteur optique embarqué dans un véhicule survolant un corps céleste, ladite image de détection étant constituée d'une pluralité de pixels, caractérisé en ce qu'il comporte les étapes : - (10) d'estimation d'un fond de l'image de détection en chaque pixel à évaluer, - (30) d'évaluation, en chaque pixel à évaluer, d'un test de détection par rapport à au moins un seuil dont la valeur en ce pixel est déterminée (20) en fonction du fond estimé en ce pixel. 2 - Procédé selon la revendication 1, dans lequel, pour chaque pixel à évaluer, on estime des caractéristiques de l'erreur d'estimation du fond en ce pixel, et la valeur de l'au moins un seuil en ce pixel est déterminée (20) en fonction du fond estimé en ce pixel et desdites caractéristiques estimées de l'erreur d'estimation du fond en ce pixel. 3 - Procédé selon la revendication 2, dans lequel les caractéristiques estimées de l'erreur d'estimation du fond en un pixel à évaluer sont au moins un écart-type de l'erreur d'estimation en ce pixel. 4 - Procédé selon la revendication 3, dans lequel la valeur de l'au moins un seuil associé à un pixel à évaluer est déterminée (20) en remplaçant, dans la détermination (20) de l'au moins un seuil, le fond estimé par une fonction dudit fond estimé et dudit écart type estimé, de préférence en remplaçant le fond estimé par la somme dudit fond estimé et d'un terme proportionnel à l'écart-type estimé. 5 - Procédé selon la revendication 2, dans lequel les caractéristiques estimées de l'erreur d'estimation du fond en un pixel à évaluer sont au moins une fonction représentative d'une densité de probabilité de l'erreur d'estimation du fond en ce pixel. 6 - Procédé selon l'une des revendications précédentes, dans lequel on utilise le fond estimé pour modéliser, dans la détermination (20) de l'au moins un seuil en chaque pixel à évaluer, un bruit instationnaire dont les caractéristiques dépendent d'un niveau du signal mesuré en ce pixel, depréférence un bruit instationnaire photonique. 7 - Procédé selon l'une des revendications précédentes, dans lequel, pour chaque pixel à évaluer, la valeur de l'au moins un seuil en ce pixel est déterminée (20) en fonction du fond estimé en ce pixel, des caractéristiques estimées de l'erreur d'estimation du fond en ce pixel, et de caractéristiques supposées du signal émis par la cible. 8 - Procédé selon l'une des revendications précédentes, dans lequel la fonction de test utilisée par le test de détection s'exprime en chaque pixel sous la forme d'une différence des distances entre le signal mesuré en ce pixel et les valeurs que ce signal devrait avoir selon au moins deux hypothèses respectivement de présence et d'absence de cible dans l'au moins une image de détection. 9 - Procédé selon l'une des revendications précédentes, dans lequel l'étape (10) d'estimation du fond met en oeuvre un filtre de Kalman, de préférence de type « sans odeur », prenant en compte un modèle d'évolution du fond de l'image de détection. 10 - Procédé selon la revendication 9, dans lequel le modèle d'évolution du fond de l'au moins une image de détection comporte une composante d'évolution spatiale et une composante d'évolution temporelle. 11 - Procédé selon l'une des revendications précédentes, dans lequel on estime (10) le fond en fonction d'une image ou d'une séquence d'images, dites « images d'estimation ». 12 - Procédé selon la revendication 11, dans lequel les images d'estimation sont des images acquises préalablement à l'image de détection, séparées de ladite image de détection par un délai non nul, égal ou supérieur au temps d'apparition (Ô) d'une cible. 13 - Procédé selon l'une des revendications précédentes, dans lequel, pour chaque pixel à évaluer, la valeur de l'au moins un seuil en ce pixel est déterminée (20) de sorte que la probabilité de fausse alarme dudit test de détection en ce pixel est théoriquement égale ou inférieure à une valeur spécifiée, la valeur spécifiée étant de préférence la même pour chaque pixel à évaluer. 14 - Procédé selon l'une des revendications 1 à 12, dans lequel on détermine(20) au moins deux seuils pour chaque pixel à évaluer, et dans lequel le test de détection est un test séquentiel. 15 - Procédé selon la revendication 14, dans lequel, parmi les au moins deux seuils pour chaque pixel à évaluer, la valeur d'un premier seuil en ce pixel est déterminée (20) de sorte que la probabilité de fausse alarme en ce pixel est théoriquement égale ou inférieure à une première valeur spécifiée, et la valeur d'un second seuil en ce pixel est déterminée (20) de sorte que la probabilité de détection en ce pixel est théoriquement égale ou supérieure à une seconde valeur spécifiée. 16 - Procédé selon l'une des revendications précédentes, dans lequel le test de détection utilise une pluralité d'images de détection. 17 - Procédé selon l'une des revendications précédentes, comportant deux phases d'exécution : une première phase, au cours de laquelle l'étape (10) d'estimation du fond et l'étape (30) d'évaluation du test de détection par rapport à des seuils déterminés (20) en fonction du fond estimé sont exécutées tant qu'une cible n'a pas été détectée dans une image de détection, une seconde phase, exécutée lorsqu'une cible a été détectée dans une image de détection, au cours de laquelle l'étape (10) d'estimation du fond n'est pas exécutée.