FR 2665601 A1 19920207 FR 9009743 A 19900731 i Procédé et dispositif de caractérisation et de localisation en temps réel de singularités d'une image numérisée, notamment pour la reconnaissance de formes dans un traitement d'analyse de scène La présente invention concerne un procédé et un disposi- tif pour caractériser et localiser en temps réel certaines singularités d'une image numérisée. L'invention sera principalement décrite dans le cadre d'une application à l'analyse d'une scène dans laquelle on cherche à reconnaître des formes caractéristiques correspondant à des con- tours à forte courbure, par exemple les angles d'un contour polygo- nal ou des régions de petite taille. Une fois ces singularités caractérisées (c'est-à-dire discri- minés du reste de l'image, notamment des points de contour voisins) et localisées (c'est-à-dire référencées par leurs coordonnées dans l'image), elles vont constituer des " points caractéristiques " corres- pondant à des " primitives " de l'image traitée qui pourront se prêter à un certain nombre de traitements tels que par exemple des seg- mentations d'image (consistant à faire croître une région à partir des " points graines " que forment les points caractéristiques préa- lablement reconnus) ou pour des applications faisant appel à des techniques de mise en correspondance de points caractéristiques. il en est ainsi, par exemple, lorsque l'on cherche à locali- ser dans un environnement naturel des objets de fabrication humaine (route, ponts, bâtiments, voies ferrées, canaux, etc) sur des images infrarouge délivrées par une caméra embarquée à bord d'un aéronef, de tels objets apparaissant en général à l'observation sous une forme polygonale Comme application typique, on peut citer le recalage de navigation de l'aéronef. Cette application à la reconnaissance de formes n'est cependant pas limitative de la présente invention, qui peut égale- ment être utilisée pour d'autres applications nécessitant l'extraction de points caractéristiques, par exemple des applications de stéréo- vision, d'analyse de mouvement, etc. De façon générale, si l'on considère l'image comme étant la restriction à N 2 d'une fonction mathématique F(x,y) à deux varia- bles et à valeurs réelles, ci-après appelée " fonction de niveau de gris ", le principe du procédé consiste à élaborer, à partir de cette fonction F (fonction image brute), une autre fonction permettant de caractériser et de localiser de façon fiable des points caractéristiques de l'image Cette autre fonction sera désignée sous le terme d'" attribut ". Il est en effet nécessaire de passer par l'intermédiaire d'un tel attribut, car si l'on se contentait de considérer les points tels que la valeur de la fonction de niveau de gris F (x,y) soit extrémale dans un voisinage de taille plus ou moins importante, il faudrait que les points caractéristiques recherchés correspondent à de simples pics de niveau de gris de l'image, ce qui constitue une hypothèse particulièrement restrictive En outre, un tel test direct serait extrê- mement sensible aux pointes de bruit de l'image. Ayant donc défini un attribut et calculé la valeur de celui-ci pour tous les pixels de l'image (étape de caractérisation), on retrouve les points caractéristiques (étape de localisation) par appli- cation d'un critère simple à cet attribut, par exemple un test de mniximalité dans l'ensemble de l'image, dans une sous-image (test local) ou bien dans un voisinage d'un pixel central donné. Un des premiers types d'attribut proposés a été un attri- but statistique, c'est-à-dire un attribut constitué d'un opérateur réa- lisant en chaque point de l'image une estimation de la variance locale dans des voisinages orientés centrés sur le point traité On pourra se référer à cet égard aux travaux de H P Moravec exposés dans Towards Automatic Visual Obstacle Avoidance, in Proceedings of the International Joint Conference on Artificial Intelligence (IJCAI), 1977 et de M J Hannah, dans Bootstrap Stereo, in Pro- ceedings of the Image Understanding Workshop, DARPA Confer- ence, 1980 Ces techniques sont caractérisées par un volume de cal- culs très pénalisant pour une mise en oeuvre matérielle. Une autre approche consiste à utiliser un attribut non plus statistique mais différentiel, c'est-à-dire dans lequel on réalise tout d'abord une approximation du gradient de la fonction du niveau de gris de l'image (c'est-à-dire que l'on calcule le vecteur formé en chaque point par les dérivées partielles de la fonction F), puis à ana- lyser la fonction vectorielle ainsi obtenue. On va ainsi définir pour chaque pixel une valeur com- plexe (c'est-à-dire, sous une autre formulation, un vecteur) conte- nant pour chaque point deux informations, à savoir: une mesure de la transition locale des niveaux de gris au voisinage de ce point, représentée par la norme G(x,y) du vecteur gradient Q(x,y), et une estimation de la direction dans laquelle s'effec- tuer cette transition, représentée par l'argument (D(x,y) du vecteur gradient Q(x,y); en cas de présence effective d'un contour au voisinage de ce point, cette direction de gradient sera donc perpendiculaire à la direction du contour. On peut alors, à partir de la mesure de la transition locale du niveau de gris, extraire les seuls points de contour, c'est-à- dire ne conserver que les points correspondants à des maxima locaux de cette fonction (c'est-à-dire des maxima dans un voisinage V donné, par exemple un voisinage 3 x 3), en ne gardant ainsi que les " lignes de crête " de la fonction de niveau de gris de l'image. Cette condition peut être formulée de la manière sui- vante: un point M(xy) donné est un point de contour si et seule- ment si la relation suivante est vérifiée: G(x,y) > G(x',y') V (x',y') e V(x,y) N D(x,y), D(x,y) désignant la droite d'orientation c( désignant un voisinage donné de (xy). L'information de gradient ne peut cependant pas être uti- lisée directement pour la recherche des points caractéristiques En effet, un point dont l'amplitude de gradient est localement maximale est nécessairement reconnu comme étant un point de contour, et il est donc délicat de le distinguer de ses voisins situés sur le même contour. Certaines techniques ont été proposées pour rechercher, après avoir ainsi détecté les points de contour, ceux pour lesquels la variation de la direction du gradient est localement maximale (angle, point de rebroussement, etc) De telles techniques sont notamment décrites par P R Beaudet dans Rotational Invariant Image Operators, in Proceedings of the International Joint Confer- ence on Pattern Recognition (IJCPR), 1978, par L Kitchen et A. Rosenfeld dans Gray Level Corner Detection, in Pattern Recognition Letters, Volume 1, 1982, par D A Zuniga et R Haralick dans Cor- ner Detection Using the Facet Model, in Proceedings of the IEEE Conference on Vision and Pattern Recognition (CVPR), 1983 et par R.L Dreschler et H H Nagel dans On the Selection of Critical Points and Local Curvature Extrema of Regions Boundaries for Interframe Matching, in Image Sequence Processing and Dynamic Scene Analysis, T S Huang editor, NATO ASI Series, Volume F 2, Springer Verlag, 1983, Ces techniques sont cependant relativement lourdes à mettre en oeuvre, car elles sont toutes basées sur l'utilisation et la combinaison des dérivées partielles d'ordre 1 et 2 de la fonction de niveau de gris de l'image; il est donc nécessaire de connaître et donc de calculer ces dérivées pour tous les pixels de l'image. Une autre technique d'analyse a été proposée par J Y. Dufour et H Waldburger dans Recalage d'images par association de primitives, construction et analyses d'histogrammes multidimen- sionnels, in Actes du congrès GRETSI, Antibes 1989. Dans cette technique, on recherche des points maximi- sant un critère de radialité locale du champ vectoriel de gradient (ce qui revient plus ou moins à rechercher les centres de courbure du contour aux endroits o cette courbure est la plus forte) On notera que ces points recherchés, qui sont à l'intersection des supports des vecteurs de leurs voisinages, sont des points proches d'un contour mais non situés sur celui-ci. Cette dernière technique, comme les précédentes, se caractérise cependant par une implantation coûteuse, en termes de nombre de circuits à utiliser. L'un des buts de l'invention est donc de proposer un trai- tement particulier, avec une architecture simplifiée adaptée, per- mettant de réaliser en temps réel ce traitement de caractérisation et de localisation des singularités d'une image numérisée. On verra en particulier que les traitements effectués ont tous un caractère local, c'est-à-dire que l'analyse d'un pixel donné est effectuée exclusivement en fonction des pixels situés autour de lui dans un voisinage limité Ceci permet d'aboutir à une relative simplicité du traitement à effectuer, à la différence des procédés de l'art antérieur qui nécessitaient généralement de considérer la glo- balité de l'image, obligeant ainsi à prévoir une architecture relative- ment complexe et coûteuse en circuits (mémoires de trame de grande capacité, volume de calculs important, etc). On notera également, lorsque l'on exposera en détail l'invention, le caractère adaptatif des traitements vis-à-vis du contenu de l'image; cette caractéristique permet d'améliorer nota- blement la qualité de ces traitements. À cet effet, la présente invention propose un procédé pour caractériser et localiser des points caractéristiques d'une image numérisée, notamment pour la reconnaissance de formes dans un traitement d'analyse de scène, ces points caractéristiques étant des points de contour à forte courbure tels que des coins ou des régions de petite taille, ladite image étant formée d'une trame bidimension- nelle de pixels présentant chacun un niveau de gris déterminé, caractérisé en ce qu'il comprend les étapes consistant à: (a) approximer, pour chaque pixel à analyser, les déri- vées partielles d'ordre 2 de la fonction de niveau de gris de l'image (F(x,y)), (b) déterminer à partir de ces dérivées un attribut (X 2) représentatif de la caractéristique recherchée et affecter cet attribut au pixel ainsi analysé, et (c) rechercher parmi les pixels de l'image ainsi analysés les points maximisant ledit attribut. Très avantageusement, ladite approximation est réalisée en convoluant la fonction de niveau de gris de l'image (F(x,y)) avec des fonctions correspondant aux dérivées partielles d'ordre 2 d'une fonction de régularisation, cette fonction de régularisation permet- tant d'atténuer le bruit présent dans l'image. L'invention a également pour objet un dispositif constitué des moyens permettant de mettre en oeuvre ces fonctions, et dans lequel la caractérisation et la localisation des points caractéristiques sont alors avantageusement opérées en temps réel à la cadence vidéo. Avantageusement, ledit attribut est la seconde valeur propre de la matrice desdites dérivées partielles: Y(x,y) = I Fxx + Fxy _ l (F Fyy)2 + 4 Fxy 2 l 1/2, avec Fxx(x,y) = a 2 F/ax 2, Fyy(x,y) = a 2 F/ay 2 et Fxy(x,y) = a 2 F/axay. De façon également avantageuse, ladite fonction de régularisation est une fonction gaussienne, notamment une fonction gaussienne centrée: G(x,y) = ( 2 X I EI 1/2)-l exp l -1/2 (x,y) E-1 (x,y)t l, o E est la matrice de covariance de la gaussienne. On va maintenant décrire un mode de mise en oeuvre de l'invention, en référence aux dessins annexés. La figure 1 est un schéma par blocs illustrant une archi- tecture préférentielle permettant de mettre en oeuvre en temps réel le traitement de l'invention. La figure 2 illustre schématiquement une première forme de réalisation du circuit convolueur de la figure 1. La figure 3 illustre une seconde forme de réalisation, plus élaborée, du circuit convolueur de la figure 1. La figure 4 illustre schématiquement une forme de réali- sation du premier circuit de fonction de la figure 1. La figure 5 illustre schématiquement une forme de réali- sation du second circuit de fonction de la figure 1. Présentation générale du traitement Le procédé consiste à traiter une image brute formée d'une trame bidimensionnelle de pixels présentant chacun un niveau de gris déterminé, par exemple l'image délivrée par une caméra vidéo telle qu'une caméra infrarouge air-sol embarquée à bord d'un aéronef. A partir de cette image brute, on va tout d'abord détermi- ner pour chaque pixel les dérivées d'ordre deux de la fonction de niveau de gris F(x,y) de l'image analysée (ou une approximation de ces dérivées) On suppose bien entendu que la fonction de niveau de gris F(x,y) est deux fois différentiable. Ces dérivées secondes peuvent être présentées sous forme matricielle (hessienne de la fonction): ( Fxx(x,y) Fxy(x,y)) H(x,y) = t Fxy(x,y) Fyy(x,y) les notations utilisées étant les suivantes: Fxx(x,y) = a 2 F/ax 2, Fyy(x,y) = a 2 F/Dy 2, et Fxy(x,y) = a 2 Ff/axay. L'idée de départ de la présente invention consiste à utili- ser un attribut basé sur l'approximation en chaque point de la fonc- tion constituée par la seconde valeur propre de la hessienne H(x,y) de la fonction F(x,y), qui est: 2 (x,y) = I Fxx + Fxy l (Fxx Fyy)2 + 4 Fxy 2 l 1/2 Le choix de la seconde valeur propre correspond à un pla- cement dans une direction orthogonale à celle o Fxx est maximale, c'est-à-dire dans un repère relatif o la variation en x est la plus marquee. L'approximation directe des dérivées partielles de la fonction image se révèle cependant le plus souvent insuffisante pour obtenir des résultats satisfaisants dans les applications habituelle- ment envisagées. En effet, outre la grande sensibilité au bruit de tels opé- rateurs, le voisinage utilisé pour l'approximation est très restreint, et l'information extraite possède donc un caractère local trop mar- qué. Pour pallier cet inconvénient, la présente invention pro- pose d'opérer non pas sur l'image originelle (fonction F(x,y)), mais sur une image modifiée (fonction F'(x,y)) obtenue en convoluant l'image originelle avec une fonction " de régularisation " G(x,y) appropriée. On aura ainsi: F'(x,y) = (F * G) (x,y) La fonction de régularisation doit bien entendu être posi- tive, intégrable et au moins deux fois continûment différentiable. Dans le cas général, on peut adopter par exemple pour cette fonction G une fonction gaussienne centrée: G(x,y) = ( 2 z I Xy 1/2)-l exp l - 1/2 (x,y) -1 (xy)t l X étant la matrice de covariance de cette fonction (qui peut être diagonale ou non). On notera que ce type de fonction n'est pas limitatif de l'invention, et que l'on peut utiliser également dans le même but d'autres fonctions de régularisation, qui peuvent être des fonctions gaussiennes ou non, et des fonctions gaussiennes isotropes ou non. Si la fonction G est convenablement choisie (ce qui est le cas de la fonction gaussienne ci-dessus), l'opération de convolution apporte deux avantages. En premier lieu, elle " régularise " l'image, c'est-à-dire qu'elle en atténue le bruit et qu'elle en arrondit les contours angu- leux. En second lieu, elle facilite les opérations de dérivation (très lourdes du point de vue du volume des calculs, donc de la com- plexité des circuits à utiliser); en effet, on peut remarquer que: an/Diiayn-i (F*G) = F * an/a Diayn-i G. La fonction l'2 que l'on utilisera comme attribut sera alors définie par: 12 (xy) = I F * LGI -l (F * Tl G)2 + (F * T 2 G)2 l 1/2, avec les notations suivantes: LG = Gxx + G, (c'est-à-dire le laplacien de la fonction G) TIG = Gx-Gyy T 2 G = 2 Gxy Avec un tel attribut, les points pour lesquels la valeur absolue de 2 '2 est élevée sont des points situés au voisinage d'un contour de l'image F' présentant une courbure. En outre, la valeur absolue de ?J 2 est d'autant plus importante que le contraste local est important et que la courbure est accentuée, ce qui correspond typiquement aux propriétés des points caractéristiques habituellement recherchés. De plus, les points tels que 212 est positive sont situés dans la concavité de la courbure locale et la présence d'un point anguleux du contour est caractérisée, si sa courbure est suffisam- ment accentuée, par la présence d'un maximum local de V'2 en une position proche de celle du centre de courbure, ou du centre de la région si le contour délimite une région de taille suffisamment faible. On voit ainsi que le point caractéristique déterminé grâce à cet attribut est proche du contour mais extérieur à celui-ci, ce qui permet avantageusement de le distinguer des points de contour pro- prement dits La même image peut donc conserver après traitement à la fois les points de contour et les points caractéristiques, ces deux types de points n'étant pas confondus. Cet attribut est donc parfaitement adapté à la caractéri- sation et à la localisation ponctuelle de coins et de régions de petite taille dans une image. Architecture de mise en oeuvre en temps réel du traitement On va maintenant décrire, en référence aux schémas par blocs des figures, une architecture de circuit propre à réaliser un temps réel les traitements du processus exposé ci-dessus. En effet, les formulations des traitements exposées plus haut ont été choisies en vue de permettre ce traitement en temps réel en tenant compte des cadences très élevées, qui peuvent attein- dre jusqu'à 20 M Hz pour la cadence de pixel. Cette architecture, illustrée schématiquement dans son ensemble figure 1, comporte essentiellement trois blocs. Le premier bloc, référencé 10, exécute à partir de la fonc- tion de niveau de gris F(x,y) reçue en entrée à la cadence de l'image vidéo les trois convolutions: L(x,y) = (F * G) (x,y), T 1 (x,y) = (F * TG) (x,y), et T 2 (xy) = (F * T 2 G) (x,y) Le second bloc, référencé 20, est un bloc permettant de calculer, à partir de deux valeurs x et y en entrée, la moyenne qua- dratique (x 2 +y 2)1/2 de ces deux valeurs; ces deux valeurs x et y seront ici les résultats de convolution T 1 et T 2 délivrés par le convo- lueur 10. Le troisième bloc de fonction, référencé 30, est un bloc permettant de calculer, à partir de deux valeurs x et y, le terme I x I -a y, a étant une constante; ces deux termes x et y seront ici, respectivement, le résultat de convolution L délivré par le circuit 10 et la moyenne quadratique (T 12 +T 22)1/2, délivrée par le premier cir- cuit de fonction 20. Le résultat délivré par ce circuit 30 est donc: I L Ia l(T 12 +T 22)1/2 l c'est-à-dire l'attribut V 2 explicité plus haut. On va maintenant décrire en détail chacun des blocs. Le bloc 10 exécutant la convolution a été représenté sur les figures 2 et 3 sous deux formes différentes. Ce circuit, dans l'une ou l'autre de ses formes, est réalisé à partir de deux circuits universels VLSI développés dans le cadre du programme européen EUR KA (projet " MIP " n EU 34: Modu- lar Image Processing), à savoir les circuits Video Memory (fonction mémoire) et Linear Filter (fitre linéaire). Le circuit MIP Video Memory est un circuit mémoire destiné à l'agencement des données pour le circuit MIP Linear Fil- ter I 1 permet de mémoriser quatre lignes vidéo de 1024 pixels cha- cune (taille maximale), chaque pixel pouvant être codé sur 8 bits de niveau de gris Il peut ainsi délivrer, à la cadence vidéo, une colonne de cinq pixels (les quatre pixels mémorisés plus le pixel courant) à un filtre linéaire placé en aval. La longueur des lignes est programmable par commande externe, avec une taille maximale de 1024 pixels. Quant au circuit MIP Linear Filter (filtre linéaire), il s'agit d'un circuit dédié permettant d'effectuer la convolution d'une image E avec un masque K, conformément à la relation: C(n,m) = X E(n+l,m+j) K(ij) Ce circuit présente les caractéristiques fonctionnelles sui- vantes: voisinage de traitement: 5 x 10, deux modes de fonctionnement possibles: mode " réel " (convolution unique) et mode " complexe " (deux convolutions simultanées avec deux masques différents), format vidéo programmable (retour ligne, retour trame), cadence vidéo maximale: 20 M Hz, entrée 5 pixels de 8 bits, sortie 16 bits (en mode réel) ou 24 bits (en mode complexe), possibilité de post-traitements intégrés: * ajustement des sorties par une transformation du type: S(n,m) = a C(n, m) 2 b + c, avec a, b et c programmables, * seuillage: les valeurs inférieures à un seuil donné peuvent être forcées à zéro, les valeurs supérieures étant conservées en l'état, ou forcées à 1 selon le mode de seuillage, * calcul d'histogramme, et * recherche du minimum et du maximum sur les valeurs résultats. Dans la première architecture envisagée pour la mise en oeuvre de l'invention, illustrée figure 2, le circuit convolueur 10 com- porte une mémoire vidéo 11 alimentant en parallèle deux filtres linéaires 12 et 13. Le premier filtre linéaire 12 opère en mode réel et permet le calcul du laplacien L Le filtre linéaire 13 opère en mode complexe (deux convolutions simultanées avec des masques différents) et délivre en sortie parallèlement Tl et T 2 Les coefficients du masque sont appliqués aux filtres linéaires par une commande externe (non symbolisée), en même temps que les autres commandes de paramé- trage de ce circuit. Cette architecture est relativement simple du point de vue du nombre de circuits utilisés, mais on peut noter que son emploi est très limité en raison de la taille réduite ( 5 x 5) des noyaux de convolution utilisables. Dans le cas de la figure 3, on a associé en cascade deux mémoires vidéo 11 et 11 ' de manière à disposer d'un voisinage 9 x 9 (la ligne courante, plus les huit lignes précédentes, chargées dans les mémoires 11 et 11 '), procurant ainsi un noyau de convolution plus large. De la même façon, on met en cascade deux groupes de deux filtres linéaires 12, 12 ' (calcul du laplacien) et 13, 13 ' (calcul de Tl et T 2). Les deux blocs de fonction 20 et 30, illustrés isolément figures 4 et 5, ont une architecture semblable. Les opérations non linéaires qu'ils impliquent peuvent être effectuées intégralement par deux circuits MIP du type Func- tion Module (module de fonction) 21, 31 associés à des RA Ms respec- tives 22, 32 En effet, le module MIP Function Module permet d'approximer à la cadence vidéo toute fonction continue à deux variables. À cet effet, la RAM qui lui est associée contient les valeurs de la fonction sur un échantillonnage de points (X,,Yj), avec O valeur de la fonction par une interpolation bilinéaire ou linéaire. Ses caractéristiques sont les suivantes: stockage des valeurs de la fonction sur une grille 128 x 128, cadence vidéo maximale: 20 M Hz, entrées: 2 x 12 bits, sortie: 12 bits. Des signaux de commande externes (non symbolisés) per- mettent de charger dans les RA Ms respectives 22, 32 les valeurs de la fonction non linéaire que l'on veut réaliser (fonction (x 2 +y 2)1/2 dans le cas du circuit 20, fonction I x I -a y dans le cadre du circuit ). REVENDICATIONS 1 Un procédé pour caractériser et localiser des points caractéristiques d'une image numérisée, notamment pour la recon- naissance de formes dans un traitement d'analyse de scène, ces points caractéristiques étant des points de contour à forte courbure tels que des coins ou des régions de petite taille, ladite image étant formée d'une trame bidimensionnelle de pixels présentant chacun un niveau de gris déterminé, procédé caractérisé en ce qu'il comprend les étapes con- sistant à: (a) approximer, pour chaque pixel à analyser, les déri- vées partielles d'ordre 2 de la fonction de niveau de gris de l'image (F(x,y)), (b) déterminer à partir de ces dérivées un attribut (X 2) représentatif de la caractéristique recherchée et affecter cet attribut au pixel ainsi analysé, et (c) rechercher parmi les pixels de l'image ainsi analysés les points maximisant ledit attribut. 2 Le procédé de la revendication 1, dans lequel ladite ap- proximation est réalisée en convoluant la fonction de niveau de gris de l'image (F(x,y)) avec des fonctions correspondant aux dérivées partielles d'ordre 2 d'une fonction de régularisation, cette fonction de régularisation permettant d'atténuer le bruit présent dans l'image. 3 Le procédé de l'une des revendications 1 ou 2, dans lequel ledit attribut (X 2) est la seconde valeur propre de la matrice des dérivées secondes: X 2 (x,y) = I Fxx + Fxy l (Fx Fyy)2 + 4 FXY 2 l 1/2, avec: Fxx(x,y) = a 2 F/ax 2, Fyy(x,y) = a 2 F/ay, et Fxy(xy) = a 2 F/axay. 4 Le procédé de la revendication 2, dans lequel ladite fonction de régularisation est une fonction gaussienne. Le procédé de la revendication 2, dans lequel ladite fonction de régularisation est une fonction gaussienne centrée: G(x,y) = ( 2 7 z I EI 1/2)-l exp l -1/2 (x,y) -1 (x,y)t l, ú étant la matrice de covariance de la fonction gaus- sienne. 6 Un dispositif de caractérisation et de localisation de points caractéristiques d'une image numérisée, notamment pour la reconnaissance de formes dans un traitement d'analyse de scène, ces points caractéristiques étant des points de contour à forte courbure tels que des coins ou régions de petite taille, ladite image étant for- mée d'une trame bidimensionnelle de pixels présentant chacun un niveau de gris déterminé, dispositif caractérisé en ce qu'il comprend: des moyens dérivateurs ( 10) approximant, pour cha- que pixel à analyser, les dérivées partielles d'ordre 2 de la fonction de niveau de gris de l'image (F(x,y)) des moyens caractérisateurs ( 23, 30), déterminant à partir de ces dérivées un attribut (X 2) représentatif de la caractéristique recherchée et affectant cet attri- but au pixel ainsi analysé, et des moyens localisateurs, recherchant parmi les pixels de l'image ainsi analysés les points maximi- sant ledit attribut. 7 Le dispositif de la revendication 6, dans lequel les moyens dérivateurs ( 10) comprennent des moyens convolueurs effec- tuant une convolution de la fonction de niveau de gris de l'image (F(x, y)) avec des fonctions correspondant aux dérivées partielles d'ordre 2 d'une fonction de régularisation, cette fonction de régulari- sation permettant d'atténuer le bruit présent dans l'image. 8 Le dispositif de l'une des revendications 6 ou 7, dans lequel la caractérisation et la localisation des points caractéristiques sont opérées en temps réel à la cadence vidéo.