Brevet FR2693011A1 — publié le 19931231

PROCEDE ET DISPOSITIF DE SURVEILLANCE D'UNE SCENE TRIDIMENSIONNELLE, METTANT EN OEUVRE DES CAPTEURS D'IMAGERIE La présente invention concerne la surveillance d'une scène éclairée, tridimensionnelle, ayant une structure géométrique connue ou déterminable, permettant notamment de détecter des anomalies de taille supérieure à une valeur déterminée (représentant par exemple des intrusions) ou, dans un mode évolué de réalisation, de localiser un mobile se déplaçant à l'intérieur de la scène. On connaît déjà de nombreux procédé de surveillance d'une scène éclairée à trois dimensions, l'éclairement étant dû à une source naturelle, telle que le soleil, ou à une source artificielle de rayonnement électromagnétique ou d'ultrason Beaucoup de ces procédés font une comparaison entre l'image de la scène à un instant donné et l'image formée précédemment Ces procédés ont notamment l'inconvénient de mal discriminer les modifications locales dues à une intrusion de celles dues à une modification locale d'éclairement De façon générale, ces procédés et dispositifs ont pour la plupart une sélectivité insuffisante pour permettre de différencier, avec une probabilité élevée, des anomalies caractérisées par leur taille, et notamment leur saillie par rapport à la scène d'origine, donc de différencier par exemple un rongeur et un être humain. La présente invention vise à fournir un procédé et un dispositif de surveillance ou de contrôle d'une scène, répondant mieux que ceux antérieurement connus aux exigences de la pratique, notamment en ce qu'ils assurent une sélectivité élevée en ne mettant en oeuvre que des moyens matériel simples. Dans ce but l'invention propose notamment un procédé de surveillance d'une scène éclairée à trois dimensions, suivant lequel: on constitue un modèle géométrique tridimensionnel (dit 3 D) de référence de ladite scène associant, à chaque point de la scène, des points homologues dans le champ d'au moins deux capteurs placés à distance l'un de l'autre et sensibles à un paramètre radiométrique (tel que la brillance) de la scène; on synthétise périodiquement la scène réelle sur le modèle géométrique tridimensionnel, à partir de chacune des représentations qui en sont fournies par les capteurs; et, à partir de là, on compare les scènes synthétisées entre elles pour faire apparaître les écarts dus à une modification du modèle 3 D qui ont une taille supérieure à une valeur déterminée. On voit que le modèle de référence est purement géométrique et n'intervient que pour permettre d'établir la correspondance entre points homologues dans les champs des différents capteurs La comparaison peut s'effectuer entre des images prises à des instants identiques, donc dans des conditions d'éclairement qui sont les mêmes Etant donné que la comparaison est effectuée entre des scènes synthétisées, une modification des conditions extérieures modifiant les valeurs radiométriques (brillance ou couleur) d'un point sera sans effet sur la comparaison, à condition que les capteurs aient des caractéristiques de réponse similaires. Il faut incidemment noter que le terme "scène" ne doit pas être interprété comme limité au cas d'une zone continue, correspondant à la zone de recouvrement des champs des capteurs La scène peut être constituée de portions discrètes dans les champs, seules jugées importantes et à surveiller. La comparaison des scènes synthétisées peut être effectuée par détermination des coefficients d'intercorrélation entre les deux scènes synthétisées, sur des voisinages déterminés des points de la scène, et par comparaison de ces coefficients à un seuil qui sera évidemment toujours inférieur à l'unité. Le modèle géométrique 3 D de la scène peut être constitué à partir des images de la scène, dans un état de référence, fournies par les capteurs, mémorisées soit pixel par pixel soit, plus généralement, à l'aide d'éléments géométriques tels que des facettes. Le modèle géométrique peut également être établi directement sous forme mathématique à partir d'une connaissance du terrain fournie par des observations antérieures ou à l'aide de cartes et les modèles de représentation géométrique en deux dimensions par les capteurs peuvent en être déduites à partir de la connaissance de la position des capteurs par rapport à la scène. Dans la plupart des cas, on sera amené à donner aux capteurs un écartement mutuel inférieur à la distance entre les capteurs et les points de la scène qui en sont les plus proches: le rapport entre l'écartement des capteurs et la distance sera un paramètre prépondérant pour sélectionner la saillie minimale que doit avoir un objet par rapport à la scène d'origine pour être détecté. Les opérations de détection définies ci-dessus peuvent être répétées, par exemple à la cadence de prise de vue des capteurs, qui seront généralement des caméras Ainsi on peut suivre les évolutions d'un objet dans la scène. L'invention propose également un dispositif de surveillance permettant de mettre en oeuvre le procédé ci- dessus défini Ce dispositif comprend plusieurs capteurs sensibles à un paramètre radiométrique, dans des longueurs d'onde comparables, comprenant des moyens de mémorisation d'un modèle purement géométrique tridimensionnel de référence d'une scène à surveiller, associant, à chaque point de la scène, des points homologues dans le champ des capteurs; des moyens permettant, à partir du modèle géométrique et des signaux représentatifs de la valeur radiométrique de chaque pixel d'image fourni par les capteurs, de synthétiser la scène réelle par projection sur le modèle géométrique tridimensionnel; et des moyens de comparaison de voisinages de dimension déterminée des scènes synthétisées, permettant de déterminer les écarts entre eux. L'invention est susceptible de nombreux modes de mise en oeuvre différents, permettant de l'adapter aussi bien au cas o l'on peut disposer ou obtenir un modèle géométrique tridimensionnel sous forme d'un modèle numérique de terrain qu'au cas o il est nécessaire de déterminer une transformation établissant, au moins de façon approchée, une correspondance entre chaque point du champ d'un des capteurs avec le point homologue du champ de l'autre capteur. L'invention sera mieux comprise à la lecture de la description qui suit de modes particuliers de réalisation, donnés à titre d'exemple non limitatifs La description se réfère aux dessins qui l'accompagnent, dans lesquels: la figure 1 montre schématiquement la disposition de deux capteurs par rapport à une scène et l'effet d'une anomalie en relief sur la scène, dans le cas d'un dispositif de surveillance que l'on peut qualifier de "fond de terrain"; la figure 2 est un schéma destiné à faire apparaître l'influence de différents paramètres sur la détection; la figure 3 est un organigramme d'un mode de réalisation du procédé; la figure 4 est un schéma de principe montrant une constitution possible d'un dispositif; les figures 5 A et 5 B montrent des dispositions particulières, avantageuses du point de vue du traitement, des images fournies par les caméras constituant les capteurs d'un dispositif; la figure 6 montre une architecture possible de la chaîne de traitement d'un dispositif. Avant de décrire de façon détaillée les étapes d'un procédé suivant l'invention de détection d'intrusions, on définira l'approche générale mise en oeuvre conformément à l'invention, dans le cas particulier d'une scène constituée par un terrain dont il est possible de représenter la structure par un modèle géométrique sous forme d'un modèle numérique de terrain (en abrégé MNT), donnant les coordonnées de chaque point du terrain dans un système d'axes fixe x, y et z. La figure 1 montre un dispositif n'ayant que deux capteurs 10 et 12, constitués par exemple par des caméras sensibles dans le domaine visible, qui ont des champs angulaires qui se recouvrent dans la scène à observer Dans l'image fournie par un capteur d'ordre i (i étant égal à 1 ou 2 dans le cas considéré) un point du modèle numérique de terrain ayant des coordonnées (x,y,z) dans le système d'axe référentiel choisi correspondra à un point de coordonnées (ai, pi) et de radiométrie Si dans un référentiel à deux dimensions qui est associé cette fois au capteur respectif. Il est possible de déterminer, par le calcul ou des mesures préalables, les transformations permettant de passer des coordonnées et de la radiométrie dans un référentiel lié à la scène aux coordonnées et à la radiométrie dans un référentiel lié à chaque capteur. Connaissant le modèle géométrique 3 D de la scène à observer et la répartition des valeurs radiométriques données par un des capteurs, on peut construire une scène de synthèse associée à ce capteur, en réaffectant à chaque point connu dans la scène, c'est-à-dire à chaque point de coordonnées (x,y,z), la valeur radiométrique fournie par le capteur. On constitue ainsi, pour chaque capteur, une scène de synthèse en associant la géométrie réelle, fournie par le MNT, et la valeur radiométrique, soit directement fournie par le capteur, soit soumise à un traitement préalable, par exemple de filtrage. Si la géométrie 3 D de la scène réelle n'a pas subi de modification dans une zone déterminée, les scènes synthétisées resteront quasi- identiques l'une à l'autre pour cette zone, même s'il y a eu évolution de la radiométrie réelle Si même, du fait par exemple d'une réponse différente des deux capteurs et des variations d'éclairement, les scènes synthétisées ne restent pas strictement identiques, les scènes synthétisées resteront dans tous les cas étroitement corrélées, en l'absence de modification de la géométrie du terrain. Si en revanche la scène observée a subi des modifications géométriques par rapport au modèle connu, par suite d'une intrusion, la valeur radiométrique correspondant à un point de la perturbation, vu d'un des capteurs, sera projetée lors de la synthèse sur une position géométrique du point d'o elle n'est pas issue, étant donné que le dispositif ne connaît que le modèle géométrique n'ayant pas cette perturbation Il en sera de même pour la valeur radiométrique fournie, pour le même point de la perturbation, par l'autre capteur En conséquence les scènes de synthèse élaborées par deux capteurs, chacun à partir des informations radiométriques qu'il reçoit et du modèle géométrique de référence, divergeront Les écarts entre les deux scènes de synthèse seront fonction de la hauteur et du volume de la perturbation. Cette influence de la hauteur de la perturbation apparaît immédiatement si l'on considère le cas particulier d'une perturbation n'ayant qu'une très faible épaisseur: dans ce cas il y a uniquement modification des valeurs radiométriques dans une zone particulière de la scène réelle. Cette scène restera géométriquement identique au modèle de référence La perturbation, comme une perturbation qui serait due à une modification locale d'éclairement (passage d'une ombre portée par exemple) ne sera pas détectée comme une intrusion. Le dispositif peut en conséquence être robuste aux variations de conditions d'éclairage et être insensible au passage d'objets de petite taille sur le sol. On verra plus loin que la taille, et notamment la hauteur, minimale détectée dépend notamment de la résolution des capteurs et de la distance entre eux. On voit également sur la figure 1 que l'anomalie 14 montrée en tirets se traduit par une modification de l'image fournie par le capteur 10 sur la largeur 1, ce qui peut se traduire comme l'adjonction d'une ombre sur l'image La même anomalie se traduit pour l'image fournie par le capteur 12, par une modification sur la largeur 1 ': les scènes de synthèse ne coïncideront donc plus sur l'ensemble des "ombres" Cela reste vrai quelle que soit la nature de l'éclairement, qu'il soit optique, électromagnétique à fréquence élevée ou même ultrasonore. La sensibilité de détection du dispositif dépend de la résolution au niveau o sont susceptibles de se produire des intrusions Si on désigne par exemple par D l'intervalle entre deux capteurs, supposés identiques, et par Ei = Pl Pi+l la résolution au niveau du terrain, pour le point i et le capteur 12, une anomalie sera détectée à la seule condition qu'elle déborde des domaines de résolution angulaire des deux capteurs, associés à un même emplacement du modèle de terrain En conséquence, l'anomalie indiquée en 14, ne sera pas détectée, alors que les anomalies 142 et 143 peuvent l'être En conséquence, on peut déterminer les conditions à remplir pour détecter une anomalie dès qu'elle dépasse une hauteur déterminée au-dessus du sol. Si en particulier le terrain est loin des capteurs et n'est pas trop perturbé, on peut considérer que la hauteur minimale d par rapport au terrain que doit présenter un obstacle pour être détecté est reliée aux autres paramètres par D/(L-d) = (L pix)/d o L est la distance moyenne des capteurs au terrain, D est la distance entre les capteurs, et pix désigne la résolution angulaire correspondant à la résolution au sol E 1 Dans le cas d'un dispositif que l'on peut qualifier de "fond de terre" destiné à détecter des intrus au sol, tout en évitant les fausses alarmes dues à des animaux de petite taille, L sera toujours très supérieur à d On peut alors écrire: L 2 = D d/pix. A titre d'exemple, on peut adopter, pour détecter des objets de taille supérieure à 10 centimètres, d= 20 cm. Le résultat peut alors être atteint avec pix = 10-3 rad L= 20 m D= 2 m. Quelle que soit la solution adoptée, on voit qu'on utilise la connaissance de la structure géométrique de la scène pour donner à la détection un caractère de comparaison géométrique, ce qui augmente les performances La connaissance que l'on a a priori élimine l'influence des éléments connus et favorise la détection des modifications géométriques dues aux intrusions. Le dispositif peut être complété par des capteurs d'autres natures, ou mettre en oeuvre plusieurs systèmes associés travaillant dans des longueurs d'onde différentes. De plus il peut mettre en oeuvre, en plus d'une détection faisant appel à une comparaison entre valeurs radiométriques fournies par deux capteurs ou davantage (par exemple en infra-rouge), une détection temporelle utilisant une comparaison entre images successives. Enfin, une fois un objet détecté, il est possible de le faire disparaître en incorporant ses caractéristiques dans le modèle géométrique de la zone observée, soit par une nouvelle prise de vue et une analyse des nouvelles vues, soit par voie purement mathématique. Le procédé, tel qu'il vient d'être décrit jusqu'ici, implique la connaissance du modèle géométrique 3 D de la scène observée pour élaborer les scènes de synthèse. Mais il suffit, pour mettre en oeuvre le procédé, de disposer d'une correspondance géométrique entre les pixels des images des différents capteurs Cette correspondance géométrique est établie pour la scène de référence afin de déceler ensuite les perturbations Toute perturbation géométriques de la scène de référence provoque une différence radiométrique entre les points associés par la correspondance dans la scène de référence Comme on le verra plus loin, cette correspondance peut être établie à partir de l'identification dans l'image de points d'amer dans les champs des différents capteurs Enfin, un modèle peut être établi en assimilant le terrain à une mosaïque de facettes ayant chacune une forme géométrique simple. Dans ces conditions, à partir de l'approche exposée ci-dessus, divers procédé de détection sont susceptibles d'être mis en oeuvre On mentionnera maintenant deux solutions, qui diffèrent par la méthode de mise en correspondance géométrique du modèle et de l'image synthétisée. La première méthode suppose qu'un modèle numérique du terrain a été réalisé et mémorisé Ce modèle numérique peut être constitué à partir de photographies aériennes, de cartes, ou même de mesures sur le terrain A partir de la connaissance de ce modèle mathématique, on peut déterminer par le calcul, pour une position donnée de chacun des capteurs, au nombre de deux au moins pour chaque scène, les transformations qui font correspondre, à chaque point de coordonnées (x,y,z) du modèle numérique, les coordonnées ai, pi de l'image de ce point dans le repère associé au capteur d'ordre i. En d'autres termes, on peut écrire, en désignant par T les transformées déterminables ai = Tad(x,y,z) Pl = Tp (x,y,z) ( 1) Le signal radiométrique S(ai,pi) fourni en chaque point de l'image du capteur i est proportionnel au flux réémis par l'élément de surface associé aux points (x,y,z) du modèle, le coefficient de proportionnalité étant sensiblement constant sur un voisinage de quelques pixels. Le schéma de détection, dans le cas de deux capteurs 10 et 12 constitués par deux caméras, est alors celui montré en figure 3 A un point (x,y,z) diffusant un flux S, on peut associer les points image (al, pl) de la caméra 1 et (a 2,p 2) de la caméra 2 On peut également en déduire les signaux Sl et 52. On peut donc ainsi créer deux objets de synthèse correspondant chacun à une des caméras, constitués par des points ayant chacun: des coordonnées (x,y,z) qui ne varient pas dans le temps, puisqu'elles correspondent à la géométrie réelle de la scène, des radiométries qui fluctuent suivant les conditions d'éclairement. L'opération de détection s'effectue par comparaison des signaux Sl et 52 associés à un même point (x,y,z) et cela pour chacun des points (x,y,z). Cette comparaison peut s'effectuer par diverses opérations mathématiques Une solution donnant de bons résultats consiste à calculer le coefficient C(p) d'intercorrélation normalisé local entre les deux modèles de synthèse, sous la forme: z Sl(Tal(P'),TPI(P')52 (Ta 2 (p'),T 2 (P')) p' V(p) C(p)= lE Sl(Tal(P'),Tl(p,))2 X E 52 (Ta 2 (P'),T 2 (P'))2 l 1/2 p' EV/ (p) o V(p) désigne un voisinage des points P respectifs, en général sur 3 x 3 ou 5 x 5 pixels dans le cas d'un réseau carré; p' désigne un point courant de V(p), de coordonnées le (x' y, z') - Le coefficient d'intercorrélation normalisé est alors égal à 1 si, pour un voisinage V(p) donné, on a 51 (T*(p'), Tpl(p')) = k 52 (T 2 (p'), T 92 (p')) Cette égalité est respectée si la scène observée n'est pas perturbée. Cette méthode permet de s'affranchir des conséquences de différences de gains entre les capteurs et des effets dus à des conditions angulaires distinctes de prise de vue, ayant une influence sur la rétrodiffusion, vers les capteurs. Dans la mesure o la scène observée est parfaitement diffusante, il suffit d'une comparaison pixel à pixel des deux images. Dans le cas d'une intrusion, le modèle de synthèse est reconstruit avec une radiométrie qui n'est plus celle associée à la géométrie de la zone observée On associe en effet, à un point de coordonnées (x,y,z) masqué par l'intrus, une radiométrie qui n'est pas issue du point de la scène de référence Et cela est vrai pour tous les points masqués par la perturbation. La détection est déterminée en émettant un signal d'alarme lorsque le coefficient de corrélation est inférieur à un seuil s, qui est évidemment choisi toujours inférieur à 1. Une telle technique de détection a l'avantage d'être peu sensible aux variations de dynamique et notamment d'éclairement, ce qui constitue un facteur de robustesse important et permet d'omettre tout étalonnage radiométrique. Une seconde approche doit être utilisée lorsque l'on ne connaît pas initialement un modèle numérique de terrain: elle permet de revenir aux conditions initiales du cas o un MNT est disponible. Une solution consiste à utiliser des points remarquables dans la scène, facilement repérables dans chaque image et dont on connaît les coordonnées (x,y,z) dans un même repère géographique, constituant des points d'amer On désignera par úXk,Yk,Zk} les coordonnées du point d'amer d'ordre k et par: {akl,pkl} les coordonnées des images du point d'ordre k dans le repère lié à la caméra 1, {ak 2,pk 2} les coordonnées des images dans le repère lié à la caméra 2. On recherche alors la transformation permettant de passer des coordonnées dans la scène réelle aux coordonnées dans l'image: -ai = Tai (x,y,z) p = Tgi (x,y,z). Suivant la complexité de la scène, on sera amené à adopter des équations de formes plus ou moins complexes pour les transformations T Le nombre de points d'amer nécessaire pour donner une représentation suffisamment exacte augmentera avec la complexité de la transformation Lorsque la scène n'est pas trop tourmentée, il est possible d'adopter des transformations Ta et T 'de forme bilinéaire: Ta(x,y,z) = (a OOO)+(al OO)x+(a Ol O)y+(a OO 1)z+(all O)xy+ (al Ol)xz+(a 011)yz TP(x,y,z) = (POO 00)+(pl OO)x+(P 010)y+(p O Ol)z+(P 10)xy+ (P 101)xz+( p Oll)yz Le calcul des coefficients a OOO, p, p Ol s'effectue en écrivant, pour l'ensemble des points d'amer, les équations: ak = Ta (xk, Yk, Zk) pki = Ti (xk, Yk, Zk) On obtient ainsi des systèmes d'équations linéaires qui peuvent être résolus sans difficulté: dans le cas ci- dessus sept coefficients sont à déterminer, ce qui exige au moins sept points d'amer Un nombre plus important peut être utilisé, de façon à retenir des coefficients représentant un optimum. Une fois les transformations Ta et T, connues, on se retrouve dans le cas précédent, puisque l'on connaît l'équivalent d'une MNT. Une autre solution consiste à utiliser une grille de transformation, projetée sur la scène réelle et dont on repère la représentation dans chaque image Cette grille peut notamment être constituée par un corroyage, obtenu par exemple par balayage à l'aide d'un laser et repérage, à chaque impulsion laser, du point de l'image représentant l'impact sur le terrain. On décrira maintenant, en faisant référence à la figure 4, la constitution de principe d'un dispositif de surveillance, destiné à la détection automatique d'intrusion dans des zones connues La figure 4 montre deux capteurs, constitués par des caméras 10 et 12 dont les champs présentent une zone de recouvrement 16 dont une partie ou la totalité constitue la scène surveillée L'éclairage peut être naturel; de nuit ou pour améliorer la sensibilité de détection à courte distance, l'éclairage naturel peut être remplacé ou complété par une source artificielle, visible ou infra-rougre Dans le cas de détection de personne, il suffit d'une résolution au sol de 10 x 10 cm et d'une cadence de deux détections par secondes, qui correspond à un déplacement de 1,5 m pour une personne à 10 km/h Les capteurs peuvent être des caméras infra-rouge à CCD ayant une matrice de 600 x 700 pixels Les deux caméras 10 et 12 sont synchronisées, par exemple en commandant la caméra 12 à partir du signal de synchronisation de la caméra 10. Les signaux vidéo provenant des caméras sont appliqués à une unité de traitement 18 qui peut être regardée comme comprenant: deux cartes d'acquisition 20 et 22 permettant de récupérer la partie utile des images des caméras, une carte processeur de signal 24, une carte de mémoire vive 26, et, éventuellement, une carte de communication 28 permettant de connecter l'unité de traitement à un réseau téléphonique. Bien que l'on puisse disposer les caméras 10 et 12 dans n'importe quelle orientation relative, il est cependant avantageux, pour réduire la complexité du traitement et la capacité de mémoire requise, de respecter des contraintes qui dans la pratique ne sont pas gênantes et dont certaines seront maintenant exposées. On rappelera tout d'abord que les pixels correspondants de deux images fournies d'une même scène par deux caméras sont sur deux droites, dites épipolaires, passant chacune par le centre optique O 1 et 02 de l'image respective: dans chaque image, le pixel correspondant au même point P de la scène est dans l'image -constitué par une droite épipolaire du plan Q P 02 Les épipolaires ont la propriété que tout pixel de l'épipolaire dans une image a un pixel homologue dans l'épipolaire correspondante de l'autre image si les repères de référence des caméras restent inchangés. Une configuration intéressante est de placer les caméras de façon que leurs repères (directions ligne et trame) soient parallèles et que les points Q 1 et 02 soient situés sur une même droite parallèle à la direction de ligne. Deux épipolaires telles que 30, et 302 sont alors des droites parallèles à la direction de ligne (figure 5 A) Les traitements peuvent alors s'effectuer ligne à ligne suivant les épipolaires: il suffit de mémoriser dans la carte 26 les lignes correspondant aux épipolaires des zones observées et un nombre de lignes en-dessus et au-dessous correspondant à l'étendue du voisinage choisi pour le calcul de corrélation. Lorsqu'il n'est pas possible ou peu commode d'adopter la disposition de la figure 5 A, celle de la figure B est encore avantageuse Les épipolaires 32, et 322 sont alors des droites inclinées, mais faisant le même angle par rapport aux lignes des deux images. Il faut alors effectuer un traitement supplémentaire, en comparaison de la figure 5 A Ou bien on sélectionne les points homologues des épipolaires et des parallèles aux épipolaires appartenant au voisinage choisi, ou bien on augmente la taille de la mémoire 26. La figure 6 montre une constitution possible de la partie traitement d'un dispositif dans le cas d'un modèle de terrain à l'aide d'une grille mémorisée On supposera que les centres optiques des caméras 10 et 12 sont disposés au voisinage d'une même ligne droite parallèle au sens de balayage dans la direction ligne Dans ce cas la mémoire 26 (figure 4) peut comporter deux espaces mémoire 26, et 262 correspondant à un paquet de quelques lignes de pixels seulement, centré sur la ligne d'ordre i dont les pixels sont voisins Le processeur peut être regardé comme ayant deux voies de traitement 24, et 242 effectuant les projections et recevant, sur des entrées respectives 34, et 342, les modèles géométriques sous forme numérique. Le processeur comporte également une unité 36 de comparaison des projections, généralement par détermination d'un facteur d'intercorrélation sur chaque voisinage et de comparaison à un seuil. La sortie de l'unité 36 est appliquée à un bloc 38 de transfert des images vers un appareil d'enregistrement et éventuellement un appareil de visualisation 40 Le processeur 36 comporte également une sortie 40 commandant le passage d'un paquet de lignes centré sur la ligne i à un paquet centré sur la ligne i+ 1 une fois la comparaison terminée En d'autres termes, le processeur pilote de façon asynchrone la mémorisation des paquets de lignes. L'architecture montrée en figure 6 permet de limiter la quantité de mémoire vive à écriture et lecture rapide. L' invention est susceptible de nombreuses variantes de réalisation et de nombreuses applications autres que celles qui ont été spécifiquement évoquées. Par exemple le dispositif peut être utilisé pour localiser un robot àl'intérieur de la scène, le robot constituant une anomalie. Dans ce cas, il peut être avantageux de placer les capteurs sur un organe de déplacement (chariot de pont roulant par exemple) muni de moyens de mesure permettant de disposer en permanence des coordonnées des capteurs par rapport au terrain Il est ainsi possible, à partir d'une mémoire centrale stockant un modèle numérique de terrain correspondant à la totalité du domaine d'évolution du robot, de mettre en permanence à jour la mémoire du processeur de traitement pour qu'elle corresponde à l'emplacement des capteurs Le dispositif suivant l'invention a l'avantage d'être exempt de la plupart des défauts des systèmes de visualisation montés sur les robots eux-mêmes, notamment la limitation du champ de vision en distance et l'impossibilité d'observer le domaine de déplacement du robot derrière des obstacles. REVENDICATIONS 1 Procédé de surveillance d'une scène éclairée tridimensionnelle, ayant une structure géométrique connue ou déterminable, permettant de détecter des anomalies, caractérisé en ce que: on constitue un modèle géométrique tridimensionnel de référence de ladite scène associant, à chaque scène, des points homologues dans le champ d'au moins deux capteurs placés à distance l'un de l'autre et sensibles à un paramètre radiométrique (tel que la brillance) de la scène; on synthétise périodiquement la scène réelle sur le modèle géométrique tridimensionnel, à partir de chacune des représentations qui en sont fournies par les capteurs; et, à partir de là, on compare les scènes synthétisées entre elles pour faire apparaître les écarts dus à une modification du modèle 3 D qui ont une taille supérieure à une valeur déterminée. 2 Procédé selon la revendication 1, caractérisé en ce qu'on compare les scènes synthétisées par détermination des coefficients d'intercorrélation entre les deux scènes synthétisées, sur des voisinages déterminés des points de la scène, et en ce qu'on compare les coefficients obtenus à un seuil. 3 Procédé selon la revendication 1 ou 2, caractérisé en ce qu'on constitue le modèle géométrique de la scène à partir des images de la scène, dans un état de référence, fournies par les capteurs, mémorisées pixel par pixel. 4 Procédé selon la revendication 1 ou 2, caractérisé en ce qu'on établit le modèle géométrique sous forme mathématique à partir d'une connaissance du terrain fournie par des observations antérieures ou à l'aide de cartes. Procédé selon la revendication 1 ou 2, caractérisé en ce qu'on constitue le modèle géométrique de la scène à partir de points d'amer de coordonnées (x,y,z) connues et d'une équation de forme prédéterminée. 6 Procédé selon la revendication 1 ou 2, caractérisé en ce que l'on constitue le modèle géométrique à partir d'une grille de corroyage, réalisée sur la scène par balayage à l'aide d'un laser et enregistrée sous forme du repérage par chaque capteur, à chaque impulsion laser, du point de l'image représentant l'impact sur le terrain. 7 Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'on compare de plus des images successives fournies par les capteurs pour réaliser une détection temporelle. 8 Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que les capteurs ont un écartement mutuel inférieur à la distance entre les capteurs et les points les plus proches de la scène. 9 Dispositif de surveillance d'une scène éclairée tridimensionnelle ayant une structure géométrique connue ou déterminable, permettant de détecter des anomalies en relief sur la scène, comprenant au moins deux capteurs ( 10,12) sensibles à un paramètre radiométrique, dans des longueurs d'onde comparables, caractérisé en ce qu'il comprend également des moyens de mémorisation d'un modèle purement géométrique tridimensionnel de référence de la scène à surveiller, associant, à chaque point de la scène, des points homologues dans le champ des capteurs; des moyens permettant, à partir du modèle géométrique et des signaux représentatifs de la valeur radiométrique de chaque pixel d'image fourni par les capteurs, de synthétiser la scène réelle par projection sur le modèle géométrique tridimensionnel; et des moyens de comparaison de voisinage de dimension déterminée des scènes synthétisées, permettant de déterminer les écarts entre eux. Dispositif selon la revendication 9, caractérisé en ce que lesdits capteurs comportent de plus des capteurs travaillant dans des plages de longueur d'onde différentes des deux dits capteurs. 11 Dispositif selon la revendication 10 de surveillance destiné à localiser un robot mobile dans la scène, caractérisé en ce que lesdits capteurs sont portés par des moyens de déplacement par rapport à la scène et sont associés à une mémoire à mise à jour du modèle géométrique mémorisé en fonction de la position des capteurs par rapport à la scène.