Brevet FR2610752A1 — publié le 19880812

La présente invention a tout d'abord pour objet un procé- dé de représentation, sur un support plan déterminé et en perspective depuis un point d'observation déterminé, de l'image dtun terrain connu par les coordonnées des noeuds d'un réseau maillé déterminé, dans lequel on détermine par interpolation une surface passant par ces noeuds, on détermine un ensemble de points sur cette surface, et on allume les pixels de l'image à représenter qui sont asso- ciés à cet ensemble de points, pour visualiser l'image sur le support. L'invention s'applique en particulier à la présentation, au pilote d'un avion, d'une image synthétique du sol survolé, ou au pilote d'un véhicule terrestre, d'une image synthétique du site dans lequel il évolue, quelles que soient les conditions de visibilité. L'invention permet également la présentation d'images de sites réels ou fictifs à des fins de simulation. Dans le cas de l'aide aux tâches de pilotage et de navi- gation, les coordonnées des noeuds du réseau maillé sont stockées dans une mémoire remplie à partir d'une base de données disponible, par exemple, auprès de l'Institut Géographique National. Une telle base de donn6es contient les altitudes de noeuds obtenus en 6chantillonnant le ter- ra toutes les trois secondes d'arc, en latitude et en rain toutes les trois secondes d'arc, en latitude et en/ 26 10752 en c'est-à-dire environ 65 m en latitude et 92m/longitude, pour le territoire français. Le point d'observation, ainsi que le support plan sur lequel est représentée l'image,dépendent de la position de l'avion par rapport au terrain, ainsi que de son attitude de façon à ce que, par exemple, l'image représentée corresponde le mieux à ce que le pilote peut voir de son habitacle dans de bonnes conditions de visibilité. Un circuit de traitement, relié A la mémoire de stockage des coordonnées des noeuds du réseau et recevant des données en provenance des instru- ments de mesure de la position et de l'attitude de l'avion, détermine donc les pixels d'image à allumer et commande en conséquence un écran de visualisation dis- posé devant le pilote. On connait déjà un procédé du type défini ci-dessus, décrit dans la demande de brevet française no 2 524 177. Dans ce procédé, la surface déterminée par interpolation et passant par les noeuds comprend un ensemble de fa- cettes triangulaires qui épouse donc approximativement la forme du terrain, et l'ensemble de points de cette surface auquel seront associés les pixels d'image allumés comprend les noeuds d'un réseau maillé auxiliaire, dif- férent du réseau maillé précédent, et dépendant ici du point d'observation. Ce réseau maillé auxiliaire se pro- jette sur un plan horizontal selon un ensemble de cercles concentriques centrés sur la projection, sur le m8me plan horizontal, du point d'observation, et selon un ensemble de droites passant par cette projection. En d'autres ter- mes, la projection du réseau auxiliaire sur un plan hori- zontal constitue une sorte de toile d'araignée centrée sur la projection du point d'observation. La régularité de cette projection du réseau auxiliaire fait que la représentation des noeuds de ce réseau auxiliaire obtenue en allumant les pixels d'image qui leur sont associés 2 6 1 0 7 5 2 donne une image du relief de la surface à facettes, donc du terrain à représenter. Toutefois, si ce procédé permet une visualisation du relief du terrain, il ne donne que peu d'informations sur le site correspondant, dans la mesure o il ne permet de représenter ni le type de végétation recouvrant le terrain, ni les voies de communication comme les rivières, les routes et les voies ferrées, ni les agglomérations, par exemple. Pour éviter cet inconvénient, on connatt des procédés permettant la présentation synthétique de paysages rela- tivement réalistes, o figurent en relief les arbres, constructions et superstructures, représentés de plus, comme le terrain, avec des couleurs pour en améliorer la lisibilité. Ces procédés présentent cependant l'inconvé- nient d'être très lourds à mettre en oeuvre, car ils né- cessitent une modélisation, à l'aide de facettes colorées, de l'ensemble du paysage, ce qui, à l'heure actuelle, ne peut être fait que manuellement, et par un opérateur qualifié. De plus, lorsque cette coûteuse modélisation a été faite, sa mémorisation occupe une place importante en mémoire, ce qui limite son application à de faibles surfaces de terrain. La présente invention vise à pallier les inconvénients précédents en procurant un procédé du type défini ci- dessus, permettant la présentation d'informations relati- vement complètes tout en demeurant assez simple pour être mis en oeuvre de façon automatique, sans intervention manuelle. A cet effet, elle a pour objet un procédé du type défini ci-dessus, caractérisé par le fait que l'on détermine une carte, ou une photographie, planimétrique du terrain avec des informations de couleur, on numérise la carte pour ob- tenir des pixels de carte de couleurs déterminées, les points dudit ensemble de points étant déterminés pour qu'ils soient tous respectivement associés à tous les pixels de l'image, on associe à chacun de tous ces points, et donc A chacun des pixels d'image, un pixel de carte, et on allu- me chaque pixel d'image dans la couleur du pixel de carte associé. Dans le procédé de l'invention, on superpose à l'image en perspective du terrain, une carte planimétrique numérisée de ce terrain. Par carte planimétrique, on entend ici une carte de type classique, comme une carte routière ou une carte d'état major, mais aussi une photographie aérienne, ou une vue satellite, sous incidence verticale. De telles cartes, numérisées, sont disponibles par exemple auprès de l'Institut Géographique National. Naturellement, le mot "couleur" est à prendre ici au sens large, c'est-A-dire correspondant soit A une représentation multicolore, soit àA une représentation du type noir et blanc, comportant une échelle de niveaux de gris. Dans le cas o l'on utilise une carte de type classique, on obtient alors une représentation en perspective du ter- rain, "habillé" A l'aide des informations qui se trouvent sur la carte, a savoir d'une part les types de végétation, le cours des rivières, le tracé des routes, les aggloméra- tions et d'autre part les noms de ces rivières, routes et agglomérations, les limites départementales, et les fron- tières par exemple. Dans le cas o l'on utilise un plan sur lequel figurent des projets de voies routières ou ferrées, ou autres, on simule en perspective l'intégration future, dans le site, de ces réalisations projetées. Dans le cas o l'on utilise une photographie aérienne, obtenue A partir d'un avion ou d'un satellite, on visualise, de façon assez réaliste et en perspective, le paysage correspondant. En d'autres termes, on peut dire qu'on obtient une resti- tution en relief de l'image à deux dimensions que consti- tue la carte planimétrique, ce qui améliore en même temps sa lisibilité et celle du relief. Avantageusement, les pixels du support plan qui ne sont pas des pixels de l'image du terrain sont allumés d'une couleur figurant le ciel. Avantageusement encore, le point dudit ensemble de points, associé à un pixel d'image est le premier point d'inter- section de ladite surface avec le rayon qui part du point d'observation et qui passe par ce pixel d'image. Alors, les parties cachées du terrain, c'est-à-dire les parties non visibles depuis le point d'observation, se trouvent éliminées. Avantageusement encore, le pixel de carte associé à un point dudit ensemble de points est le pixel de carte de mêmes coordonnéesplanimétriques que ce point. La présente invention a également pour objet un système pour la mise en oeuvre du procédé précédent de représen- tation, sur un support plan déterminé et en perspective depuis un point d'observation déterminé, de l'image d'un terrain connu par les coordonnées des noeuds d'un réseau maillé déterminé, comprenant: - une première mémoire pour stocker les coordonnées desdits noeuds, 2 6 1 0 7 5 2 - des moyens pour engendrer un signal de données repré- sentatives du point d'observation et du support plan, - des moyens reliés à ladite première mémoire et auxdits moyens de génération, pour déterminer par interpolation une surface passant par lesdits noeuds, puis En ensemble de points sur cette surface, et les pixels de l'image à représenter qui sont associés à cet ensemble de pàints,et - des moyens de visualisation reliés auxdits moyens de détermination pour allumer lesdits pixels de l'image à représenter, système caractérisé par le fait que: - il comprend une deuxième mémoire pour stocker une carte, ou une photographie, planimétrique du terrain avec des in- formations de couleur, numérisée sous forme de pixels de carte de couleurs déterminées, - lesdits moyens de détermination sont reliés à ladite deuxième mémoire et agencés pour déterminer les points dudit ensemble pour qu'ils soient tous respectivement associés à tous les pixels de l'image, pour associer à chacun de tous ces points, donc à chacun des pixels d'ima- ge, un pixel de carte et pour allumer chaque pixel d'image dans la couleur du pixel de carte associé. La présente invention sera mieux comprise grâce à la description suivante de la mise en oeuvre préférée du procédé de l'invention, ainsi que de formes de réalisation du système de l'invention,faite en se référant aux dessins annexés, sur lesquels, - la figure 1 représente schématiquement un terrain survolé par un avion, ce terrain étant associé à une carte numérisée, - la figure 2 représente une vue détaillée de la carte numérisée de la figure 1, - la figure 3 représente une vue de la représentation synthétique, en perspective, du terrain de la figure 1, obtenue A l'aide du procédé de l'invention, - la figure 4 représente une vue de la représentation synthétique, en perspective, du relief seul du terrain de la figure 1, - la figure 5 est un diagramme illustrant les procédés d'interpolation pouvant être utilisés dans le procédé de l'invention, et, - la figure 6-est une représentation schématique, par blocs, d'un système mettant en oeuvre le procédé de l'invention. En référence à la figure i un procédé de représentation syn- thétique d'un terrain 2, pour le pilote d'un avion 8 survo- lant ce terrain 2, est décrit A titre d'exemple. Le point d'observation est le point P et un support plan 1, ici rectangulaire, figure par exemple la fenêtre de l'habi- tacle de l'avion, A travers laquelle peut être vu le ter- rain 2 survolé, si les conditions de visibilité sont bon- nes. Sur la figure 1, et pour des raisons évidentes de clar- té, le support plan 1 est considérablement agrandi. Le pro- cédé de représentation synthétique décrit se propose de fournir A l'observateur une image en perspective, synthé- tique, du terrain survolé, comparable & l'image vraie, mais indépendante des conditions de visibilité. Cette image syn- thétique, produite & partir d'informations concernant le terrain 2 stockées en mémoire, et d'informations en prove- nance des instruments de bord concernant la position et l'attitude de l'avion 8 par rapport au terrain 2, est pré- sentée A l'observateur de façon connue, par exemple A l'aide d'un écran de même taille que le support plan 1 ou de taille différente, ou encore par l'intermédiaire d'un dispositif connu de superposition avec l'image vraie. Les informations concernant le terrain 2 stockées en mémoire comprennent d'une part des informations alti- métriques et d'autre part des informations planimétriques, ces deux types d'informations étant fournies ici et pour le territoire français, par l'Institut Géographique - National. Les informations altimétriques comprennent les coordonnées de noeuds N d'un réseau maillé obtenu en échantillonnant le terrain toutes les trois secondes d'arc, c'est-à-dire, et comme cela a déjà été indiqué, environ 65 m en latitude et environ 92 m en longitude. Pour chacun de ces noeuds, on connaIt donc, dans un repère trirectangle OX, OY, OZ lié au sol, la coordonnée X de latitude, la coordonnée Y de longitude, et la coordonnée Z d'altitude. Les informations planimétriques comprennent une carte 5 du terrain, avec des informations de couleur, représentant ici et à titre d'exemple, le type de végétation, le cours des rivières, le tracé des routes, les agglomérations, ainsi que les noms de ces rivières, routes et aggloméra- tions, et d'autres indications de ce type. La carte 5 est numérisée, c'est-à-dire qu'elle a été divisée en une pluralité de surfaces élémentaires, ou éléments d'images, ou pixels, tels que les pixels 7 et 7' de la figure 1. La taille de ces pixels est, de façon connue, choisie assez petite pour que l'on puisse affecter à chacun d'entre eux une couleur uniforme, choisiedans un ensemble fini de couleurs quantifiées, sans perte notable d'information. De telles cartes sont disponibles auprès de l'Institut Géographique National. La figure 2, qui est une photogra- phie d'un écran cathodique, fournit un exemple de carte 2 6 1 0752 numérisée, représentant en particulier des routes et des agglomérations. Naturellement une telle carte peut être ai- sément mémorisée puisqu'elle est totalement définie par la connaissance des coordonnées planimétriques, c'est-A-dire de latitude X et de longitude Y, et de la couleur quantifiée de chacun des pixels de carte. Du fait que l'altitude du terrain 2 n'est connue que pour les noeuds N et que, comme on va le voir dans la suite, il sera utile de connattre cette altitude pour certains points ne coïncidant pas nécessairement avec ces noeuds N, on détermine une surface 3 connue partout, et épousant approximativement la forme du terrain 2. Pour cela, et de façon connue, on détermine une surface 3 par interpola- tion à partir des coordonnées des noeuds N. Sur la figure 1, on a représenté la surface 3 obtenue en remplaçant chaque portion de terrain comprise entre trois noeuds par une face plane triangulaire. Il existe naturellement d'autres types d'interpolation qui seront décrits dans la suite, après une première description d'ensemble du procédé. Après avoir défini le type d'interpolation adopté, on se trouve en mesure de déterminer, pour chaque pixel 4 du support plan 1 qui appartient à l'image du terrain 2, le point Q de la surface 3 associé à ce pixel 4, et que ce pixel 4 représente, c'est-à-dire le premier point d'inter- section Q de cette surface 3 avec le rayon qui part du point P et passe par le pixel 4. La démarche suivie pour déterminer le point Q sera décrite dans la suite. Les coordonnées planimétriques X et Y de ce point Q per- mettent de déterminer le pixel 7' de la carte 5 associé à ce point Q de la surface 3 c'est-à-dire, pratiquement, du terrain 2. 26 1 0752 Comme on le remarque sur la figure 1, une ligne horizontale, par exemple, de pixels d'image 4 correspond A une suite de points Q de la surface 3 qui dessine dans l'espace une cour- be gauche. Cette courbe gauche ne se projette évidemment pas sur la carte 5 selon une droite, mais selon une courbe quel- conque- Il en résulte que, en allumant chaque pixel d'ima- ge 4 dans la couleur du pixel de carte 7' qui lui est asso- cié, on obtient une représentation déformée de la carte 5, comme si cette carte était dessinée sur la surface 3 elle- même au lieu d'être dessinée sur un plan- l'ensemble est donc représenté en perspective. La figure 3, qui est une photographie d'un écran cathodique, montre le résultat obte- nu à partir de la carte de la figure 2 et des informations altimétriques du terrain correspondant. Sur la figure 3, qui montre la déformation des routes et des agglomérations de la figure 2, on voit que l'on a réalisé une superposi- tion de la carte 5 et de l'image en perspective du terrain 2. Sur la figure 3, afin d'améliorer encore la lisibilité de l'image obtenue, on a allumé chaque pixel d'image 4 dans la couleur du pixel de carte 7' associé, avec une intensité qui dépend de l'éclairement du terrain 2, au voisinage du point Q associé à ce pixel d'image 4, c'est-A-dire du point Q que représente ce pixel d'image 4. A cet effet, on choisit une direction d'éclairement du terrain, par exemple qui coïncide avec celle du soleil, et on calcule l'éclaire- ment au voisinage de chaque point Q- Cette technique, connue dans le cas o l'on affecte A la surface 3 représentée une couleur uniforme, et décrite plus en détails dans la suite, permet de restituer une bonne lisibilité au relief, comme le montre, dans ce cas, la figure 4. La figure 4, toujours une photographie d'un écran cathodique, montre un "effet de dunes" qui représente, en fait, le résultat obtenu dans le cas de la représentation du relief seul de la portion de terrain représentée sur la figure 3. Naturellement, cette modulation de l'intensité lumineuse de il chaque pixel, si elle améliore fortement la lisibilité de l'image obtenue, n'est pas toujours obligatoire. Par exem- pie, dans le cas o la carte utilisée est une photographie aérienne, celle-ci comporte généralement un ombrage natu- rel, qui suffit & donner une idée du relief. Ce relief ap- parait également du fait que certaines parties de la sur- face et du terrain étant cachées, certaines parties de la carte ne seront pas représentées. En effet, comme on asso- cie à chaque pixel d'image le premier point d'intersec- tion Q de la surface 3 avec le rayon qui part du point P et qui passe par le pixel 4, c'est-A-dire le point d'inter- section le plus proche du point P, les parties de la sur- face 3 non visible du point P ne sont pas représentées. Naturellement, lorsque le rayon qui part du point P et qui passe par un pixel du support plan 1 n'intercepte pas la surface 3, ou l'intercepte A une distance supérieure A un seuil définissant l'horizon, on considère que ce pixel n'est pas un pixel d'image 4 et on l'allume d'une couleur figurant le ciel, par exemple la couleur bleue. Ceci peut se produire si le rayon passe au-dessus de la ligne d'hori- zon, ou si les données mémorisées, tant altimétriques que planimétriques, ne concernent qu'une portion du terrain qui se trouve dans le champ de vision défini par le point P et le support plan 1. En référence A la figure 5, différents exemples de déter- mination, par interpolation, de la surface 3 sont décrits. Sur cette figure le point E représente un point de coor- données planimétriques X et Y, dont il faut déterminer l'altitude Z sachant qu'il se trouve dans la maille du réseau maillé définie par les quatre noeuds A, B, C, D de coordonnées (XA, YA' ZA), (XC, YA' ZB)' (XC' YC' ZC)' et (XA, YC' ZD), respectivement. Un premier type d'interpolation, assez grossière, consis- te à affecter au point E l'altitude Z du noeud le plus proche, ici le noeud D: Z = ZD Dans ce cas, le terrain 2 est assimilé à une surface 3 formée de rectangles horizontaux, de surface égale à celle d'une maille, et centrés sur chaque noeud. Un deuxième type d'interpolation est l'interpolation par facettes représentée sur la figure 1, qui revient à affec- ter au point E une altitude Z dépendant des altitudes des trois noeuds les plus proches, ici les noeuds A, D et C, selon la formule suivante, bilinéaire en X et Y Z =)Z +Z A +.. - Z (YA-YC) (Xc-XA)(Yc-YA) D (XC-XA) C Un troisième type d'interpolation est l'interpolation par surfaces non planes se projetant selon la maille ABCD, qui revient à affecter au point E une altitude Z dépendant des altitudes des quatre noeuds A, B, C et D de la maille, selon la formule suivante, également biliné- aire en X et Y: (X-Xc) (Y-YC) (X-X XA) (-YC) Z + (x -XC)(YA-Yc) (Xc-XA) (YA YC) (X-XA) (Y-YA) (X-Xc) (Y-CA) ZC + - ZD (Xc-XA) (YC-YÀ) (XA-Xc) (YC-YiA) Les trois types d'interpolation qui précèdent, ainsi que d'autres, sont bien connus de l'homme de métier, qui est apte à choisir celle qui lui parait la mieux adaptée, compte tenu du compromis à réaliser entre la précision de l'interpolation et la complexité des calculs qui en ré- sultent. Lorsque le type d'interpolation déterminant la surface 3 à laquelle on assimile le terrain 2 est définie, on calcule le point d'intersection Q associé à chaque pixel du support plan 1, c'est-à-dire à chaque pixel d'écran, en procédant de la façon suivante. On décompose le support plan 1 en une matrice de pixels qui correspondent aux pixels de l'écran grace auquel l'image sera visualisée. On procède colonne par colonne, en commançant à chaque fois par le pixel le plus bas de la colonne, qui est, en principe, celui qui est le plus proche de la surface 3 et pour lequel on est presque certain de trouver un point d'intersection Q avec la surface 3. Pour le pixel le plus bas d'une colonne, on utilise un point E, situé sur le rayon qui part du point P et qui passe par ce pixel, et dont on fait varier la position. Initialement, ce point E est sur le support plan 1. Ainsi, on est certain que le point E est au-dessus de la surface 3. Après avoir choisi un pas de progression D sur le rayon, on augmente la distance PE pas par pas, en regar- dant, à chaque fois, la position relative du point E par rapport à la surface 3. Dès que le point E se trouve en dessous de la surface 3, on procède à une recherche du point d'intersection par la méthode connue de dichotomie, c'est-à-dire que l'on réduit la distance PE de D/2, 26 1 0752 puis D/4, et ainsi de suite, c'est-à-dire par pas de valeur chaque fois égale à la moitié de celle du pas précédent, Jusqu'à ce que le point E passe au-dessus de la surface 3. Alors on augmente à nouveau la distance PE par pas de valeur chaque fois égale à la moitié de celle du pas précédent, c'est-à-dire de D/2n+l, D/2 n+2, et ainsi de suite, si la valeur D/2n était celle qui a fait passer le point E au dessus de la surface 3. On continue jusqu'à par exemple, une augmentation de distance égale à D/2P, qui fait passer le point E au-dessous de la sur- face 3, et ainsi de suite. Comme cela est connu, si l'on arrête lorsque la variation de distance entre deux points E successifs est inférieure à un seuil donné, par exemple D/2N et si l'on appelle E' le dernier point E obtenu, on est certain que la distance entre E et le point d'intersection Q est inférieure au seuil. En choisissant N assez grand, on peut donc affecter au point Q les coordonnées du point E' sans faire une erreur importante. En principe, comme il s'agit du pixel le plus bas de la colonne, on trouve un point d'intersection Q à distance finie. Si on ne trouve pas de point d'intersection alors que la distance PE dépasse un seuil correspondant sensi- blement à l'horizon, on considère que le pixel n'est pas un pixel 4 de l'image du terrain 2 et on l'allume de la couleur du ciel, ainsi que les pixels suivants de la colonn0. 3o Si, comme cela est en général le cas, on a trouvé le point d'intersection Q à distance finie, le pixel considéré appartient à l'image du terrain 2, et on l'allume dans la couleur du pixel 71 de la carte 5 qui lui est associé, c'est-à-dire du pixel 7' de même coordonnées planimétriques 2 6 1 0752 X et Y que le point d'intersection Q. On passe alors au pixel. suivant de la colonne, pour lequel on procède à la recherche du point d'intersection Q de la surface 3 et du rayon qui part du point P et qui passe par ce pixel en procédant sensiblement comme pour le pixel le plus bas, à ceci près que, d'une part, on part d'un point E' dont la distance initiale au support plan 1,au lieu d'être nulle, est égale à la distance PQ trouvée pour le pixel précédent, et que, d'autre part, le pas de progression sur:le rayon est choisi égal à D' avec D' inférieur à D. Comme, en principe, les points Q associés à deux pixels d'image 4 adjacents ne sont pas très éloi- gnés, on gagne ainsi en temps de calcul. Pour améliorer la lisibilité de l'image obtenue, on peut allumer chaque pixel d'image 4 avec une intensité variable, dépendant du point Q du terrain auquel il est associé, et qu'il représente. A titre d'exemple, on peut faire dépendre l'intensité lumineuse du pixel d'image 4 de la distance PQ, pour que cette intensité soit d'autant plus faible que la distance PQ est grande. Une loi inversement proportionnelle au carré de la distance convient bien. On peut également faire dépendre l'intensité lumineuse de l'altitude du point Q représenté, pour que cette intensité soit d'autant plus faible que l'altitude est élevée, ou encore en traitant cette altitude par paliers discontinus, pour faire apparattre des courbes de niveau. On peut également, et comme cela a déjà été signalé et représenté sur les figures 3 et 4, définir une direction d'éclairement à l'aide d'un vecteur normalisé, et faire 2 6 1 0 752 dépendre l'intensité lumineuse de l'éclairement du point Q représenté. A cet effet, on calcule la normale à la surface 3 en Q, éventuellement par une méthode approchée de type connu, et on effectue le produit scalaire de la normale et du vecteur d'éclairement. Ce produit scalaire dépend de l'angle entre les deux vecteurs. L'intensité lumineuse dépend du module du produit scalaire si celui-ci est négatif, et si ce produit est positif, l'intensité lumineuse est nulle car, bien que le point Q se trouve vu depuis le point d'observation P, il n'est pas éclairé par la source d'éclairement. On notera par ailleurs que, dans le procédé employé, le problème de l'élimination, sur l'image représentée, des parties cachées, c'est-à-dire des parties qui ne se trouvent pas vues depuis le point d'observation P, ne se pose pas dans la mesure o l'on choisit toujours le pre- mier point d'intersection Q sur le rayon en provenance du point P. Le procédé qui vient dt'être décrit est, par exemple, bis en oeuvre à l'aide du système de la figure 6. Ce système comprend principalement un circuit de traite- ment 13, relié à deux mémoires 11 et 14 et au système 12 de navigation inertielle de l'avion 8, et commandant un dispositif de visualisation 15. Le système 12 de navigation inertielle de l'avion est du type connu qui délivre un signal de données S représenta- tif de la position de l'avion et de son attitude donc du point d'observation P et du support plan 1. La mémoire 11 contient des données altimétriques concer- nant le terrain 2 survolé. Elle est pourvue d'une entrée 26 1 0 752 recevant un signal d'adressage SN et d'une sortie déli- vrant un signal ZN. Le signal d'adressage SN représente les coordonnées X et Y d'un noeud N du réseau maillé et le signal de sortie ZN l'altitude Z de ce noeud N. La mémoire 14 contient les données planimétriques du ter- rain 2 survolé, c'est-à-dire la carte numérisée 5. Elle comprend une entrée recevant un signal d'adressage SC et une sortie délivrant un signal PC. Le signal d'adressage SC représente les coordonnées X et Y d'un pixel de carte 7 et le signal de sortie PC représente la couleur de ce pixel de carte 7. Le dispositif de visualisation 15 est du type connu com- portant un écran, par exemple un écran cathodique, orga- nisé en une matrice de pixels, comportant quatre entrées recevant les signaux I, J, IT et PC' correspondant res- pectivement à la ligne, la colonne, l'intensité lumineuse et la couleur du pixel à allumer. Le circuit de traitement 13 est agencé pour déterminer par interpolation la surface 3, puis chaque point Q de cette surface 3 associé à chaque pixel d'image 4, c'est-à-dire chaque point Q d'intersection de la surface 3 avec le rayon qui part du point P et qui passe par ce pixel d'image 4, et enfin chaque pixel de carte 7' associé à chaque point Q, donc à chaque pixel d'image 4. Le circuit de traitement 13 comprend un circuit 136 qui, en réponse au signal S, calcule, dans le repère lié au * sol OX, OY, OZ les coordonnées du point P et du support 1. Il délivre un signal de données CPS. Un circuit de balayage 135 délivre les signaux I et J de façon à commander le balayage de l'écran du dispositif de 2 6 1 0752 visualisation 15 colonne par colonne, en commençant, dans chaque colonne par le pixel le plus bas. Un circuit 132 calcule, en réponse au signal CPS et aux signaux I et J, les coordonnées du point E, ou du point E', situé sur le rayon partant de P et passant par le pixel défini par les signaux I et J. Le circuit 132 délivre deux signaux SC et Z représentatifs le premier des coordonnées planimétriques X et Y et le second de l'al- titude Z du point E, ou du point E'. Un circuit 131 de calcul de la surface 3 par interpola- tion, relié à la mémoire 11 avec laquelle il échange les signaux SN et ZN, reçoit le signal SC et délivre un signal ZS qui représente l'altitude du point de la surface 3 ayant les coordonnées planimétriques X et Y du point E, ou du point E'. Un circuit de comparaison 133, recevant les signaux Z et ZS, compare leurs valeurs. Il est pourvu d'une sortie délivrant un signal binaire R. Le signal R est au niveau haut lorsque Z est supérieur à ZS, c'est-à-dire lorsque le point E, ou le point E' est au-dessus de la surface 3. Il commande alors le circuit de calcul 132 pour faire avancer le point E, ou le point Et, sur le rayon, c'est- à-dire pour augmenter la distance PE. Le signal R est au niveau bas lorsque Z est inférieur à ZS, c'est-à-dire lorsque le point E, ou le point E', est au dessous de la surface S. Il commande alors le circuit de calcul 132 pour faire reculer le point E, ou le point E', sur le rayon, c'est-à-dire pour réduire la distance PE. Le circuit 132 est agencé pour effectuer la dichotomie précédemment décrite et lorsque la distance entre deux points successifs est inférieure au seuil D/2 défini ci-dessus, il délivre un signal binaire EG de commande d'un circuit de porte 134, qui permet d'appliquer au dispositif de visualisation le signal PCt identique au signal PC représentant la couleur du pixel de carte 7' dont les coordonnées planimétriques X et Y sont définies par le signal SC, c'est-à-dire qui sont celles du point d'intersection Q. Ainsi, on peut dire que le circuit de porte 134 permet d'associer à chaque point d'intersection Q un pixel de carte 7'. De même, on peut dire que le circuit de calcul 132 et le circuit de comparaison 133 associent, à un pixel d'image 4, le point d'intersection Q de a surface 3 avec le rayon qui passe par ce pixel d'image 4. Le circuit 137 calcule ici l'intensité lumineuse avec laquelle il faut allumer le pixel défini par les signaux I et J pour tenir compte de l'éclairement de la surface 3, c'est-à-dire du terrain 2. A cet effet, il reçoit les signaux SC, Z et EGjde façon à disposer des coordonnées du point d'intersection Q et il est relié à la mémoire Il avec laquelle il échange les signaux SN et ZN pour calcu- ler la normale en Q à la surface 3. A partir de la direc- tion d'éclairement défini par un signal EC extérieur et que l'opérateur peut modifier, il calcule l'intensité lumineuse à affecter au pixel comme cela a été décrit précédemment. 3O Naturellement, le circuit de traitement 13 peut être réalisé également à l'aide d'un circuit de calcul unique qui accomplit séquentiellement les tâches des circuits 131 à 137. A titre d'exemple, les pages 20, 21 et 22 décrivent un organigramme illustrant le fonctionnement d'un tel circuit. 2 6 1 0 7 5 2 B3Eousle s.ur le colonnez de l'écran J=Jdeb,Jfi n Boucle sur les lignes 102- I = Ideb,lfin Calcul des coordonnées du point E de 1'écran 103 done Io reprte liS 4 l'écren Pas.ege dans le repère lié au sol,O/j- E('X,YZ) Non /,/ r 05'oJ iNon I=Ideb? Oui b N'x B /,.10. A 1 Calcul de la direction Calcul de la direction du rayon PE du rayon PE Calcul de l'intersection Q Calcul de l'intersection Q du rayon avec le terrain du rayon avec le terrain Non H orssecteur de visualiaion = iel PIX0lJ> est obtenu fo au point Q t I,lfin Non 1il Non 261 0752 i initial1setion " 1061 D!NT = j CIEL Fausse. Calcul des incréments de progression sur le rayon: DX = D*Dray*(X-XYI) 106 DY = D*DrauI*(Y-YVI) DZ = D*Dra.*(Z-ZYI) 1;".E' Y'=Y'+DE' "I Z'= Z+ DZ /, -% 0Non C' E' appartient-il au secteur visualis ? / E'est-il eu dessous du terrain? > Non \., %,,,-/5 Dichotomie sur les coordonnees de E et E' ci jusqu'à ce que: Dist(E',E) DT / |I$OT R Q (intersection) = E',1, 69 Calcul de DINT: DINT = Dist(P,Q) 26 1 0752 CIEL Peusse Calcul des incréments de progression sur le rayon: DX = D'*Dra.*-(X-X'fI) DY = D'*Dray*(YW-YVI)10 DZ = D'*Draj*g(Z-ZVI) Calcul des coordonnées X' Y',Z' du point E' situé sur le rauon tel aue: 10 Dist(P,E') = DIN E'appasrtient -il au secteur visualisl ? Position de E' par rapport au terrain (au dessCs ou au dcessous) DX = '* DiV, "., E" Y' Y'+ DY H Z" =Z' + DZ Non -s E"et E's0nt-ils de part et d'autre du terrain ? JCIEL vre Z I NTHi jusqu'è ce que: Dist(E', E") 2 6 1 0 75 2 Sur la page 20 les blocs 101 à 105 n'appellent pas de commentaires particuliers. Le bloc 106 est détaillé sur la page 21. Dans le bloc 1061, la variable DINT représente la distance entre point P et le point d'intersection Q. Dans les blocs 1061 et 1065 la variable CIEL est "fausse" si le pixel d'écran est un pixel 4 de l'image du terrain et "vraie" dans le cas contraire, o le pixel correspondant est allumé de la couleur du-ciel. Dans le bloc 1062, la variable Dray représente l'inverse de la norme du segment de rayon qui passe par le point d'observation P, de coordonnées XVI, YVI et ZVI et le point E de coordonnées X, X, Z. La variable D représente le pas de progression sur le rayon. Dans le bloc 1066, la distance DHOR est la distance au delà de laquelle on peut considérer qu'on a atteint la ligne d'horizon. Le bloc 107 est détaillé sur la page 22 à propos duquel des remarques similaires à celles qui ont été faites à propos du bloc 106 peuvent être faites. Dans le bloc 109 de la page 20, le tableau IPIX représente le tableau des pixels d'écran et la variable Nciel repré- sente la couleur choisie pour le ciel. On remarque que, lorsqu'un pixel d'écran n'est pas un pixel d'image, on affecte systématiquement la couleur du ciel à tous les pixels de la colonne qui sont disposés au- dessus, ce qui permet de gagner en temps de calcul sans risque d'une trop grande perte d'information. Dans le bloc 110, le tableau IMA représente les pixels de la carte 5 numérisée. Naturellement, la présente invention n'est pas limitée à la description qui vient d'en &tre faite. En particulier il est à la portée de l'homme de métier de modifier certaines étapes, de façon connue, pour réduire le temps de calcul. Par exemple, tant que le point E, ou le point E', se trouvent relativement loin du terrain 2, il n'est pas obligatoire de calculer, par interpolation bilinéaire, l'altitude du point de la surface 3 qui lui correspond. Une simple interpolation grossière peut suffire, le calcul de l'interpolation bilinéaire n'étant alors entre- pris que lorsque le point E est assez proche du terrain 2. Revendications 1. Procédé de représentation, sur un support plan déter- miné (1) et en perspective depuis un point d'observation déterminé (P), de l'image d'un terrain (2) connu par les coordonnées des noeuds (N) d'un réseau maillé déterminé, dans lequel on détermine par interpolation une surface (3) passant par ces noeuds (N), on détermine un ensemble de points (Q) sur cette surface (3), et on allume les pixels (4) de l'image à représenter, qui sont associés cet ensemble de points (Q), pour visualiser l'image sur le support (1), procédé caractérisé par le fait que ou un _hotographie, l'on détermine une carte,5) planimétrique du terrain (2) avec des informations de couleur, on numérise la carte pour obtenir des pixels (7, 7') de carte de couleurs déterminées, les points (Q) dudit ensemble de points étant déterminés pour qu'ils soient tous respectivement associés à tous les pixels (4) de l'image, on associe à chacun de tous ces points (Q), et donc à chacun des pixels d'image (4), un pixel de carte (7') et on allume chaque pixel d'image (4) dans la couleur du pixel de carte (7') associé. 2. Procédé selon la revendication 1, dans lequel les pixels du support plan (1) qui ne sont pas des pixels (4) de l'image du terrain (2) sont allumés d'une couleur figurant le ciel. 3. Procédé selon l'une des revendications i ou 2, dans lequel le point dudit ensemble de points, associé à un pixel d'image (4)> est le premier point d'intersection (Q) de ladite surface (3) avec le rayon qui part du point d'observation (p) et qui passe par ce pixel d'image (4). 4. Procédé selon l'une des revendications 1 à 3, dans lequel le pixel de carte associé à un point (Q) dudit ensemble de points est le pixel de carte (71) de mêmes coordonnées planimétriques que ce point (Q). 5. Procédé selon l'une des revendications 3 ou 4, dans lequel, pour déterminer ledit point d'intersection (Q): - on part d'un point (E; E') situé sur ledit rayon à une distance initiale, de valeur déterminée, dudit support plan (1), - on détermine la position relative dudit point (E; E') par rapport à ladite surface (S) et si ledit point (E; E') se trouve au-dessus de ladite surface (S), on augmente sa distance du support plan (1) par pas de valeur constante (D; D') jusqu'à ce que ledit point (E; E') se trouve au dessous de ladite surface (S), - on réduit alors ladite distance, par pas de valeur cha- que fois égale à la moitié de celle du pas précédent (D/2, D/4,..., D/2n) jusqu'à ce que ledit point (E; E') se trouve au-dessus de ladite surface (3), - on augmente alors ladite distance, par pas de valeur chaque fais égale à la moitié de celle du pas précédent (D/2n+l, D/2n+2,..., D/2p) jusqu'à ce que ledit point (E; E') se trouve au dessous de ladite surface (S), - on procède ainsi jusqu'à obtenir un point (E'; E") dont la distance au point précédent (E; E') est inférieure à un seuil déterminé (D/2N) et on affecte audit point d'in- tersection (Q) les coordonnées de ce point (E'; E"). 6. Procédé selon la revendication 5, dans lequel les pixels du support plan (1) sont organisés en une matrice ayant une pluralité de lignes et une pluralité de colonnes et on procède colonne par colonne, comme suit: on détermine si le pixel le plus bas de la colonne du support plan (1) est un pixel (4) de l'image du terrain (2) en recherchant le premier point d'intersection (Q) de ladite surface (3) avec le rayon qui part du point d'ob- servation (P) et passe par ce pixel, à partir d'une dis- tance initiale de valeur nulle entre le support plan (1) et ledit point (E) sur ce rayon, - si ce pixel n'est pas un pixel d'image (4) on l'allume d'une couleur figurant le ciel, ainsi que tous les pixels suivants de la colonne, - si ce pixel est un pixel d'image (4) on l'allume dans la couleur du pixel de carte (7') associé et on passe au pixel suivant de la colonne, pour lequel on recherche le premier point d'intersection (Q), en prenant cette fois comme distance initiale entre le support plan (1) et ledit point (E') sur le rayon, la distance entre le point d'observation (P) et le point d'intersection (Q) précé- demment déterminé. 7. Procédé selon l'une des revendications 1 à 6, dans lequel on allume chaque pixel d'image (4), dans la couleur du pixel de carte (7') associé, avec une intensité d'au- tant plus faible que la distance entre le point d'observa- tion (P) et le point (Q) dudit ensemble de points, asso- cié à ce pixel d'image (4), est grande. 8. Procédé selon l'une des revendications 1 à 6, dans lequel on allume chaque pixel d'image (4), dans la couleur du pixel de carte (7') associé, avec une intensité liée à l'altitude du point (Q) dudit ensemble de points, asso- cié à ce pixel d'image (4). 3o 9. Procédé selon l'une des revendications 1 à 6, dans lequel une direction d'éclairement du terrain (2) étant déterminée, on allume chaque pixel d'image (4), dans la couleur du pixel de carte (7') associé, avec une intensité qui dépend de l'éclairement de la surface (3) au voisinage du point (Q) dudit ensemble de points, associé à ce pixel dWimage (4). 10. Système pour la mise en oeuvre du procédé selon l'une des revendications 1 à 9 de représentation, sur un support plan déterminé (1) et en perspective depuis un point d'observation déterminé (P), de l'image d'un terrain (2) connu par les coordonnées des noeuds (N) d'un réseau maillé déterminé, comprenant: - une première mémoire (11) pour stocker les coordonnées desdits noeuds (N), - des moyens (12) pour engendrer un signal (S) de données représentatives du point d'observation (P) et du support -plan (1), - des moyens (13) reliés à ladite première mémoire (il) et aux dits moyens de génération (12), pour déterminer par interpolation une surface (3) passant par lesdits noeuds (N), puis un ensemble de points (Q) sur cette sur- face (3), et les pixels (4) de l'image à représenter qui sont associés à cet ensemble de points (Q), et - des moyens de visualisation (14) reliés auxdits moyens de détermination (13) pour allumer lesdits pixels (4) de l'image à représenter, système caractérisé par le fait que: - il comprend une deuxième mémoire (14) pour stocker une ou uei_ photographie, carteD5 planimetrique du terrain (2) avec des informa- tions de couleur, numérisée sous forme de pixels (7, 7') de carte de couleurs déterminées, - lesdits moyens de détermination (13) sont reliés à ladite deuxième mémoire (14) et agencés pour déterminer les points (Q) dudit ensemble pour qu'ils soient tous respectivement associés à tous les pixels (4) de l'image, pour associer à chacun de tous ces points (Q), donc à chacun des pixels d'image (4), un pixel de carte (7') et pour allumer chaque pixel d'image (4) dans la couleur du pixel de carte (7') associé. 11. Système selon la revendication 10, dans lequel les- dits moyens de détermination (13) comprennent: - des moyens (131) de calcul de ladite surface (3) reliés à ladite première mémoire (11), - des premiers moyens (132, 133), reliés auxdits moyens de génélation (12) et auxdits moyens de calcul (131), pour associer/un pixel d'image (4) le premier point d'inter- section (Q) de ladite surface (3) avec le rayon qui part du point d'observation (P),et passe par ce pixel d'ima- ge (4), et, - des deuxièmes moyens (134) reliés auxdits moyens d'as- sociation (132, 133) et à la deuxième mémoire (14), pour associer à chacun des points d'intersection (Q) le pixel de carte (7') sur lequel se trouve le point de même coor- données planimétriques que ce point d'intersection (Q).