Brevet FR2699781A1 — publié le 19940624

Procédé de détection d'apparition d'objets ponctuels dans une image La présente invention concerne un procédé de détection d'apparition d'objets ponctuels dans une image numérisée sous la forme d'une pluralité de pixels ayant des luminances respectives Les objets à détecter dans l'image sont ponctuels et ont donc une taille sensiblement de l'ordre d'un pixel, typiquement quelques pixels. Il est connu dans le domaine de l'imagerie infrarouge de représenter une image par un modèle à base d'un polynôme à deux variables, et dont des valeurs prises par ce polynôme sont une approximation des luminances respectives des pixels de l'image. L'invention vise à modéliser une image de pixels à des fins de détection d'objets ponctuels, ou cibles, dans l'image. A cette fin, un procédé de détection d'apparition d'objets ponctuels dans une image numérisée sous la forme d'une pluralité de pixels ayant des luminances respectives, est caractérisé en ce qu'il comprend pour chacun de blocs élémentaires de pixels centré sur l'un respectif desdits pixels de l'image, les étapes suivantes: une approximation polynomiale des luminances de pixels d'un anneau dudit bloc, ledit anneau étant composé des pixels dudit bloc situés à la périphérie d'un trou composé des autres pixels dudit bloc et localisé au centre dudit bloc, une prédiction des luminances desdits pixels du trou en fonction de ladite approximation polynomiale afin de produire des luminances prédites, une soustraction desdites luminances prédites aux luminances des pixels dudit trou, en des luminances de différence, et une détection d'objet ponctuel dans ledit bloc lorsque le rapport de l'une desdites luminances de différence relatives audit trou et d'une moyenne de luminances de différence sur l'ensemble de ladite image, est supérieur à un seuil prédéterminé. L'approximation polynomiale selon l'invention consiste en la détermination de filtres de convolution à partir de bases polynomiales orthonormées relatives respectivement à des dimensions de l'image, et en des convolutions desdits pixels de l'anneau avec lesdits filtres de convolution en des coefficients respectifs d'un polynôme. D'autres caractéristiques et avantages de la présente invention apparaîtront plus clairement à la lecture de la description suivante d'une réalisation préférée de l'invention, en référence aux dessins annexés correspondants dans lesquels: la figure 1 montre un bloc-diagramme schématique d'un système de pistage infrarouge; la figure 2 est un bloc-diagramme d'un module de détection selon l'invention inclus dans le système de pistage infrarouge montré à la figure 1; la figure 3 est un bloc-diagramme détaillé d'une réalisation de circuits d'estimation et de prédiction inclus dans le module de détection; la figure 4 est une représentation tridimensionnelle de luminances associées à un bloc de pixels de dimension ( 2 M + 1)( 2 N + 1) = 3 x 3; la figure 5 montre un bloc de pixels de dimension 3 x 3 avec un trou à un pixel; et I la figure 6 est un signal d'image de différence de luminances. En référence à la figure 1, un système de pistage infrarouge comprend typiquement une barrette de détecteurs infrarouges 1 disposée verticalement et associée à un dispositif de focalisation la, un module de détection d'apparition d'objets ponctuels 2 et un module de pistage 3. La barrette de détecteurs 1 et le dispositif de focalisation la sont montés sur un support tournant lb pour être entraînés en rotation à une vitesse angulaire constante autour d'un axe vertical YY' Des sorties de la barrette de détecteurs 1 sont appliquées à une entrée d'un dispositif d'acquisition d'images infrarouges de l'espace environnant 2 a pour en fournir cycliquement une image de pixels à un rythme suffisamment rapide. L'espace environnant est ainsi analysé périodiquement pour un angle de gisement de 360 degrés Le dispositif d'acquisition 2 a inclut, en outre, des moyens de balayage pour sélectionner tour à tour chacune desdites sorties des détecteurs respectifs constituant la barrette 1 Ces moyens de balayage sont synchronisés avec la commande de rotation de la barrette 1 incluse dans le support la. En sortie, le dispositif d'acquisition 2 a fournit un signal vidéo, ici numérique, représentant une image infrarouge glissante de l'espace environnant, par rafraîchissement cyclique. Cette image glissante consiste en une pluralité de pixels Chaque pixel est associé à une surface élémentaire d'image et une intensité de rayonnement infrarouge détectée à un instant donné par un détecteur de la barrette 1. Pour chaque tour complet de l'ensemble 1-la, le signal vidéo numérisé sortant du dispositif d'acquisition 2 a supporte un certain nombre de luminances Pij de pixels, par exemple codées en mots s à 8 bits, chacun de ces mots de luminance Pi,j étant représentatif du quotient de la luminosité émise par une portion élémentaire respective (pixel) de l'image de l'espace environnant par l'aire de cette portion. Ce signal vidéo numérisé sortant est appliqué à une entrée du module de détection 2 Ce dernier détecte dans l'image de l'espace environnant des cibles ou objets fixes ou mobiles quasi-ponctuels, et produit des alarmes en réponse aux cibles détectées Ces cibles quasi-ponctuelles sont des portions très réduites de l'image ayant des dimensions typiquement de 1 à ( 4 x 4) pixels Les alarmes produites par le module 2 sont appliquées à une entrée du module de pistage 3 Ce dernier suit l'évolution temporelle des cibles dans l'image afin "d'accrocher" ces cibles, par exemple dans un contexte de poursuite d'une cible se déplaçant sur le sol. L'invention est orientée vers le module de détection 2 Un tel module de détection 2 peut être intégré dans un système différent de celui représenté à la figure 2, tel qu'un simple système de visualisation de cible dans une image. En référence à la figure 2, un module de détection d'apparition d'objets ponctuels 2 selon l'invention comprend un circuit d'estimation sur anneau 20, un circuit de prédiction sur trou 21, un comparateur 22, un circuit de calcul de moyenne 23 et un circuit de détection 24. Préalablement à la description du module de détection d'apparition est présenté un modèle mathématique connu d'approximation polynomiale d'une image de pixels Ce modèle polynomial, notamment utilisé dans le domaine de l'imagerie, consiste en la représentation d'une image par un polynôme à deux variables qui est l'approximation des luminances de s pixels de l'image situés au voisinage d'un pixel central de coordonnées discrètes (uo, vo), et qui s'écrit sous la forme: L(u O + u, v O + v) = ko,O + k 1,o u + kol V + k 1,1 u v + k 2,o O U 2 + ko,2 v 2 + k 2,1 u 2 v + kl,2 u v 2 + dans laquelle équation kr s, o r et S sont des entiers, désigne un coefficient du polynôme, L désigne une fonction d'intensité, ici la fonction de luminance, et (u,v) dénote un couple de coordonnées, abscisse et ordonnée, de pixel dans un bloc d'image. Le polynôme L constitue un développement en série de Taylor au voisinage du pixel central de coordonnées (u 0, vo). En référence à la figure 5, l'approximation polynomiale est réalisée sur chacun des blocs élémentaires B Ei,j, carrés ou rectangulaires, centrés sur un pixel respectif Pi,j de l'image infrarouge. Les blocs élémentaires ont une dimension ( 2 M+l)x( 2 N+ l), o M et N sont des entiers Plus précisément, cette approximation polynomiale est utilisée selon l'invention pour estimer des luminances de pixels formant un anneau dans chaque bloc élémentaire B Ei,j Le polynôme résultant de cette approximation polynomiale des luminances des pixels de l'anneau sert à prédire les luminances ponctuelles dans un trou de ce bloc élémentaire B Ei j le trou comportant un ou plusieurs pixels et étant localisé au centre du bloc élémentaire B Ei,j. L'anneau du bloc B Eij est constitué par tous les pixels du bloc B Ei,j situés à la périphérie du trou, dans lequel est effectuée une détection d'objet ponctuel, ou cible En pratique, plusieurs dimensions de bloc élémentaire B Ei,j et de trou peuvent être choisies en fonction notamment de la taille probable de la cible à détecter dans l'image A titre d'exemple, des blocs élémentaires B Ei,j de pixels ont des dimension de ( 5 x 5) à ( 7 x 7) et des trous associés ont des dimensions de 1 et ( 4 x 4) pixels. Le modèle mathématique de représentation polynomiale de l'anneau d'un bloc élémentaire de pixels B Ei,j est le suivant Une image de pixels peut être représentée par un plan sur lequel des points discrets (pixels) sont définis, et auxquels sont associées des luminances respectives L'image est naturellement supportée par un plan à deux dimensions Pour les deux variables dimensionnelles u et v de l'image sont définies respectivement deux bases polynômiales orthonormées de Tschebycheff Le choix des dimensions de ces deux bases résulte d'un compromis entre une bonne localisation des pixels et une faible sensibilité aux bruits, notamment de texture. Des vecteurs des deux bases polynomiales orthornomées respectives de dimension E et F peuvent être définis de la façon suivante: Pa (u) = 1 Pb (V) = 1 Pl u P ((u) =v Pa(u) =i Pa(u) = M 2/ MO u Pb(v) =y 2/ YO V Pa(u) = x 3 ( 4 / 2)u Pb(v) = v 3 (Y 4 / Y 2)v etc etc 7 2699781 et ainsi de suite, par récurrence en considérant que dans chaque base, les polynômes doivent consituer une base orthogonale La norme suivante, sous forme de produit scalaire discret, est utilisée pour définir chacune des deux bases: P Pa P ≥ Pe(U) Pa (u), uv f f 1 Pb Pb ≥ Lpb(v) P ( ( 2) v V o e et f sont respectivement compris entre O et (E- 1), et O et (F-i), en considérant que les coefficients des termes de plus haut degré dans les polynômes respectifs de plus haut degré des deux bases sont égaux à 1, c'est-à- dire que ces polynômes sont normalisés. La définition de ces deux bases orthonormées de vecteurs polynomiaux s'effectue relativement à des intervalles U et V de tailles ( 2 M + 1) et ( 2 N+l) correspondant aux dimensions d'un bloc élémentaire de pixels B Ei,j. Les coefficients kk et Yk' o k est un entier, cités précédemment et définissant les coefficients des polynômes sont alors donnés par: U=+M V=+N k Yuk k= 1 YVku ( 3) u =-M v =-N Il est possible de définir, à partir de ce qui précède, des masques dérivateurs, permettant de créer le polynôme de luminances f(uo + u, vo + v) qui est l'approximation des luminances respectives des pixels dans l'anneau d'un bloc B Ei,j. Ces masques Kr s, de dimension égale à la dimension d'un bloc élémentaire B Ei, j, telle que ( 7 x 7), sont donnés par: Are F, Sa (U) pf(v) Krslu,vl = peau ( 4) e+f bases polynomiales respectivement. Le polynôme de luminance est alors défini par: L(uo+u, vo+v) = ko,0 + kl,o O u + ko,l v + kl,lu v + k 2, u 2 + ko,2 v 2 + ( 5 a) avec krs= Krs lu,vl I(uo + u,v O + v) ( 5 b) (u,v) e UXV I(u,v) désignant la luminance d'un pixel de coordonnées (u,v). Les deux dernières relations ( 5 a) et ( 5 b) sont à la base de la présente invention L'approximation d'une image par un polynôme à deux variables est connue en soi L'invention propose d'utiliser un tel modèle d'approximation-polynomiale d'image dans le cadre d'une détection de cible. Selon la présente invention, l'approximation polynomiale de chaque bloc élémentaire de pixel est limitée aux seuls pixels de l'anneau du bloc, et les luminances ponctuelles dans le trou dudit bloc élémentaire B Eij sont prédites à partir de cette approximation polynomiale Ces deux étapes sont matérialisées par les circuits 20 et 21 dans la figure 2. En revenant aux équations précédentes ( 5 a) et ( 5 b), il apparaît que le polynôme de luminance f est entièrement connu à partir des coefficients krs, r et S étant des entiers Dans la relation ( 5 a), les variables u et v désignent les abscisse et ordonnée d'un pixel, et f représente la luminance de ce pixel. Comme souligné par la relation ( 5 b), les coefficients krs sont entièrement définis à partir de convolutions spatiales des masques Krs et des luminances des pixels du bloc considéré B Ei-j Ces masques constituent des filtres linéaires de convolution spatiale En pratique, seuls quelques coefficients krs du polynôme de luminance L sont nécessaires pour estimer convenablement une luminance de pixel Comme montré par la relation ( 5 b), chaque coefficient du polynôme de luminance nécessite pour son calcul la définition d'un masque, ou filtre, de convolution Krs L'approximation polynomiale étant réalisée selon l'invention, à partir des seuls pixels de l'anneau, dans la relation 5 (b), le calcul des coefficients du polynôme de luminance ne s'effectue qu'en considérant les pixels de coordonnées (u,v) appartenant à l'anneau. La figure 3 illustre une réalisation pratique des circuits 20 et 21 de la figure 2, sous la forme de six filtres de convolution 200 à 205 associés aux masques Kr,s, et un sommateur 206, tels que O polynomiale. En revenant à la figure 2, l'image prédite Ip et l'image d'origine infrarouge Io, constituée de(s) luminance(s) réelle(s) de(s) pixel(s) du trou, sont fournies à deux entrées respectives du comparateur 22 Ce dernier produit en sortie une image de différence, notée ID, représentative de la différence de luminance entre image d'origine Io et image prédite Ip. Ainsi est obtenu le contraste entre un(des) luminances réelle(s) du(des) pixel(s) du trou et un(des) luminance(s) de ce trou prédite(s) à partir du modèle d'approximation polynomiale. La figure 6 montre un signal d'image de différence entre luminances prédites et réèlles de pixels de trou pour les différents blocs élémentaires B Eij d'une image Il apparaît que ce signal de différence de luminance fluctue autour d'une différence moyenne de luminance notée D. En se reportant à la figure 2, l'image de différence ID est appliquée à une entrée du circuit de calcul de moyenne 23 qui produit la différence moyenne de luminance D Cette différence moyenne de luminance D est représentative du contraste moyen entre les luminances réelle(s) et prédites pour tous les trous dans les blocs élémentaires B Eij de l'image Cette différence moyenne de luminances D est appliquée à une première entrée du circuit de détection 24 dont une seconde entrée reçoit l'image de différence ID = (Io Ip) Pour chaque trou de chaque bloc élémentaire de l'image I(, le(s) pixel(s) de l'image de différence ID correspondante est(sont) comparé(s) avec la différence moyenne D des luminances prédites et réelles Une alarme est produite par le circuit 24 lorsque l'un (des) pixel(s) d'un trou de l'image de différence ID a(ont) une(des) luminance(s) très supérieure(s) à la différence moyenne de luminance D. En pratique, il est choisi un coefficient multiplicateur C qui est multiplié à la différence moyenne de luminance D pour fournir un produit (c D) qui est comparé à la luminance de chaque trou d'image de différence dans le circuit 24; en d'autres termes cette comparaison est équivalente à celle du rapport ID/D et du coefficient C Une augmentation de ce coefficient C abaisse le taux de fausses alarmes Le choix du coefficient C résulte d'un compromis entre il un faible taux de fausses alarmes et un faible taux de cibles non détectées. A titre d'exemple, la figure 4 illustre une détection de cible pour un pixel central d'un trou. Sur cette figure, deux luminances d'un pixel Pij d'un trou respectivement indiquées en trait continu et en trait discontinu sont associées à des niveaux réel et prédit par approximation polynomiale L'écart entre ces deux luminances déclenche une alarme dans le circuit de détection 24 s'il est supérieur au produit de la différence moyenne des luminances réelles et prédites D sur l'ensemble de l'image par le coefficient C. lb 13 2699781 REVENDICATIONS 1 Procédé de détection d'apparition d'objets ponctuels dans une image numérisée sous la forme d'une pluralité de pixels ayant des luminances respectives, caractérisé en ce qu'il comprend pour chacun de blocs élémentaires de pixels (B Eij) centré sur l'un respectif desdits pixels de l'image, les étapes suivantes: une approximation polynomiale ( 20) des luminances de pixels d'un anneau dudit bloc, ledit anneau étant composé des pixels dudit bloc (B Eij) situés à la périphérie d'un trou composé des autres pixels dudit bloc et localisé au centre dudit bloc (B Ei, j), une prédiction ( 21) des luminances desdits pixels du trou en fonction de ladite approximation polynomiale afin de produire des luminances prédites, une soustraction ( 22) desdites luminances prédites (Ip) aux luminances des pixels (Io) dudit trou, en des luminances de différence (ID), et une détection ( 23, 24) d'objet ponctuel dans ledit bloc (B Eij) lorsque le rapport de l'une desdites luminances de différence (ID) relatives audit trou et d'une moyenne de luminances de différence (D) sur l'ensemble de ladite image, est supérieur à un seuil prédéterminé (C). 2 Procédé conforme à la revendication 1, caractérisé en ce que ladite approximation polynomiale consiste en la détermination de filtres de convolution (Kr,s) à partir de bases polynomiales orthonormées (papb) relatives respectivement à des dimensions (u,v) de l'image, et en des convolutions desdits pixels de l'anneau avec lesdits filtres de 14 2699781 convolution (Krs) en des coefficients respectifs (kr,s) d'un polynôme (L). 3 Procédé conforme à la revendication 2, caractérisé en ce lesdites luminances prédites sont des valeurs prises par ledit polynôme (L) pour des valeurs de coordonnées desdits pixels du trou.