PROCEDE DE RECTIFICATION EN TEMPS REEL D'IMAGES DE SATELLITES METEOROLOGIQUES GEOSTATIONNAIRES

La présente invention est relative à un pro¬ cédé de rectification en temps réel d'images de satel- lites météorologiques géostationnaires.

Les satellites météorologiques géostation¬ naires actuellement en orbite (satellites qui tournent à environ 36 000 kms de la Terre au-dessus de l'Equateur et dont la rotation dure 24 heures de sorte qu'ils semblent immobiles dans le ciel) sont constitués par le satellite européen METEOSAT, les satellites américains GOES et le satellite japonais GMS, lequels, pour la prise d'images météorologiques, adoptent tous le même principe de ba¬ layage à savoir, construire des lignes de données par ro- tation du satellite autour de son propre axe longitudinal rendu parallèle à l'axe Nord-Sud (NS) de la Terre.

Aucun satellite géostationnaire ne peut être considéré comme étant en position fixe par rapport à la Terre : en fait l'orbite, l'orientation (attitude), la vitesse de rotation et le début du balayage s'écartent dans le temps par rapport aux valeurs idéales. Ceci donne des images qui sont déformées par rapport à une image de référence, que l'on suppose prise dans des conditions no¬ minales (à savoir, absolument stables) . Cette situation exige que la déformation des images soit déterminée avec une haute précision et que les données relatives aux images transmises par le satel¬ lite et reçues à partir de celui-ci soient corrigées avant utilisation. Les satellites évoqués plus haut peuvent ef¬ fectuer une telle correction seulement après qu'une image complète de la Terre ait été reçue. En particulier, pour un satellite tel que METEOSAT, ceci entraîne un retard qui peut aller jusqu'à 40 mn pour la partie de l'image correspondant à l'hémisphère Sud et jusqu'à 20 mn pour l'Europe, avant que l'image corrigée puisse être diffusée aux utilisateurs. Bien que de tels retards soient, en principe, acceptables pour les météorologistes utilisant les images pour les prévisions météorologiques classiques sur des périodes allant jusqu'à une semaine, ils ne sont pas du tout acceptables pour la description détaillée du temps atmosphérique conjointement avec les prévisions mé¬ téorologiques obtenues par extrapolation avec anticipa¬ tion pouvant aller jusqu'à 2 heures (dite "no casting" suivant la terminologie anglo-saxonne, cf. la définition citée par KA BROWNING dans son introduction au "Proceedings of the second international symposium on nowcasting, edt. B. BATTRICK and E. ROLFE, Norrkoeping, Sweden, 3-7 September, 1984) . L'état de la technique le plus proche est pré¬ cisément défini par le satellite européen METEOSAT qui est un satellite stabilisé par rotation autour de son axe et qui prend une image de la Terre chaque demi-heure. Un tel intervalle d'une demi-heure entre les images reçues est dénommé "créneau". Quatre images sont prises simulta¬ nément, une dans l'infrarouge, deux dans le spectre vi¬ sible et une dans la bande de vapeur d'eau. On fait ba¬ layer le télescope qui équipe le satellite de 18° dans la direction SN, ce qui donne un balayage complet de la Terre de 2500 lignes en 25 mn. Dans chacune de ces 2500 lignes, constituant l'image infrarouge thermique ou de vapeur d'eau, se trouvent 2500 éléments (ou points) d'images ou pels ("pixels" suivant la terminologie anglo- saxonne) . Toutefois, un seul canal visible a 2500 lignes de 5000 pixels, et, étant donné que deux canaux visibles agissent simultanément chacun regardant des lignes d'images alternées, c'est un total de 5000 lignes de don¬ nées correspondant à la bande visible qui sont dispo¬ nibles, c'est-à-dire 5000 x 5000 pixels pour l'image vi- sible. Les données d'images reçues par la station d'acquisition de données, de télécommande et de poursuite (DATTS) , sont envoyées au système de calcul terrestre ME¬ TEOSAT (MGCS) situé au Centre d'Opérations de 1'AGENCE SPATIALE EUROPEENNE à DARMSTADT (ESOC) .

Le satellite européen METEOSAT, en service de¬ puis le 9 décembre 1977 se trouve sur le méridien 0 (GREENWICH) , au-dessus de l'Equateur. Son système de prise d'images voit toujours la même portion, environ 1/3 de la surface de la Terre et l'image brute provient d'un radiomètre qui "tâte" la Terre en 25 mn. Un système de miroirs discrimine les rayonnements en trois bandes spec¬ trales et les conduit à trois détecteurs : l'un est sen¬ sible à la lumière visible (500 à 1000 nm) , l'autre au rayonnement infrarouge à 10500/12500 nm (rayonnement de chaleur), et le troisième, pour l'infrarouge aussi, mais de 5700 à 7100 nm (vapeur d'eau) . Le signal de chaque dé¬ tecteur est amplifié, filtré, changé d'analogique en di¬ gital et ensuite mémorisé (la Terre vue de 36 000 kms tient dans un angle de 18 ^* ).

La principale station terrestre du système METEOSAT est située au sud-est de Francfort, en Allemagne de l'Ouest. L'image brute reçue est transmise à partir de là par câble au Centre ESOC. Les images sont alors intro- duites dans un système de traitement où elles sont corri¬ gées géométriquement, calibrées et évaluées.

Les images dans le spectre visible consistent en 5000 lignes de 5000 points par ligne : c'est-à-dire, 25 millions de points d'image. Chaque point d'image peut prendre 64 valeurs de brillance possible. Les images in¬ frarouges sont formées de 2500 lignes de 2500 points, dans lesquelles les points sont codés en 256 valeurs. Les images dans la bande de vapeur d'eau consistent en 6,25 millions de points, comme pour les images dans la bande infrarouge, les points étant eux aussi codés en 256 va¬ leurs.

Les points mentionnés représentent une résolu¬ tion de 2,5 km en lumière visible et 5 km dans l'infrarouge, les valeurs optimales étant validées pour la surface directement sous le satellite (Golf de Gui¬ née) .

Les images brutes sont traitées et évaluées de plusieurs façons au Centre de DARMSTADT, afin de les rendre disponibles pour les recherches météorologiques. METEOSAT est alors utilisé comme émetteur-relais pour la distribution des images sélectionnées à l'intention des centres Météo, des Instituts de recherche et d'autres parties intéressées. Ces images traitées sont transmises au satellite à partir de la station terrestre de DARM¬ STADT et retransmises par le satellite aux fréquences ap¬ propriées.

D'autres renseignements très utiles concernant METEOSAT peuvent être trouvés dans la publication "Introduction to the METEOSAT operational System" ESA- BR32, September 1987, ainsi que dans l'article de WOLFF "An image geometry model for METEOSAT, Int. J. Remote sensing, 1985, vol. 6, n ° 10, pages 1599 à 1606.

Actuellement, les données d'images sont démul- tiplexées par un ordinateur dit intermédiaire (ce qui équivaut à l'expression anglo-saxonne "Front end Compu¬ ter") et stockées sous forme d'images brutes, dans une mémoire à disque, puis traitées successivement suivant trois étapes principales, qui sont illustrées à la figure 1 accompagnant la présente description, à savoir :

1) détermination des paramètres géométriques du système d'images ;

2) calcul de la matrice de déformation (à sa¬ voir, de l'écart de l'image réelle par rapport à l'image de référence) ; 3) rectification des images brutes évoquées plus haut suivant les indications de la matrice de défor¬ mation précitée.

Ces trois étapes correspondent à ce que la littérature scientifique désigne sous l'expression "traitement d'images géométriques", qui vise à recons ^¬ truire géométriquement l'image nominale.

L'image rectifiée est distribuée par le satel¬ lite METEOSAT à un grand nombre d'usagers comprenant les centres nationaux de prévision du temps ainsi que les amateurs.

Le traitement informatique de 1 ' image brute permet de déterminer les paramètres qui sont nécessaires pour connaître la déformation géométrique de l'image brute correspondant au créneau k-1. Cet ensemble de para¬ mètres, désigné sous la référence P _jç -i, est utilisé pour calculer la déformation géométrique sur un certain nombre de points de grille (26 x 26 dans l'application visée dans le cadre de la présente invention) . Les variables à considérer sont les suivantes : i) la vitesse de rotation du satellite autour de son axe longitudinal (dit axe de roulis) ; ii) l'orbite et l'orientation du satellite ; iii) les paramètres du satellite qui sont en relation avec le pointage (à savoir, le repérage dans l'espace) du détecteur ; iv) le centrage Est-Ouest (EO) du disque correspondant à la Terre dans l'image.

En ce qui concerne la vitesse de rotation propre du satellite, celle-ci est mesurée par le système de minutage dont est équipé le satellite.

L'orbite et l'orientation du satellite sont déterminés selon un traitement en différé (à savoir, en mode autonome) à partir des données de télémesure de ME- TEOSAT et prédites pour les jours suivants. Cette in- formation orbitale est utilisée sans modification. L'orientation basée sur la télémesure est utilisée comme valeur initiale, alors que les horizons infrarouges Sud et Nord sont utilisés pour améliorer la précision de l'orientation. Cette amélioration est nécessaire parce que METEOSAT n'est pas équipé, -pour des raisons écono¬ miques-, d'un système de mesure d'orientation sophistiqué comme c'est le cas d'un traceur de cartes stellaires, en sorte que 1'orientation mesurée ne présente pas la préci- sion exigée pour le traitement d'images géométriques.

Les horizons Est et Ouest dans l'image brute sont utilisés pour améliorer la position Sud-Nord de 1'image.

L'ensemble de paramètres ainsi affinés P^- _i est utilisé pour déterminer la déformation de l'image brute sur les points de grille évoqués plus haut. La dé¬ formation sur chacun de ces points est un vecteur repré¬ sentant la déformation de l'image en lignes et pixels. L'ensemble des valeurs correspondant à ces déformations est mis sous forme de matrice de déformation D^.-.^. Une méthode précise de détermination de cette matrice est dé¬ crite par WOLFF dans l'article "A simple approach to solve the METEOSAT image déformation problem based in ho¬ rizon extraction from image data and orbit information", Proc. Int. Symp. Spacecraft Flight Dynamics, DARMSTADT, FRG, 18-22 May, 1981, ESA SP-160,Aug. 1981, ainsi que dans l'article -également de WOLFF- déjà évoqué plus haut. Or, la matrice de déformation constitue l'entrée du programme de rectification qui effectue l'interpolation de la déformation sur la base des points de grille préci¬ tés, de : anière à déterminer la déformation pour chacun des pixels dans l'image rectifiée, et utilise cette dé¬ formation pour corriger le pixel correspondant à l'emplacement dans l'image brute. Le concept, d'ailleurs simple, d'un satellite stabilisé par rotation autour de son axe longitudinal, conjointement avec l'utilisation des données d'images elles-mêmes pour améliorer la précision des paramètres géométriques, fournit une erreur RMS (moyenne quadra¬ tique) de moins de 2 pixels. Cette erreur RMS a été dé¬ terminée en mettant certains points de contrôle en cor¬ rélation avec des données de référence (cf MUNDOCART Car- tographical Database, Petroconsultants -CES- Ltd. CAM- BRIDGE, UK) . La méthode de mesure a été décrite en détail par ADAMSON, J. , KERR, GW, JACOBS, GHP dans l'article "Rectification of quality assessment of Meteo- sat images", ESA Journal, vol. 12, pages 467 à 482, 1988) . Comme conséquence de cette méthode existante, le traitement d'images géométriques peut être mis en route seulement après que le disque complet de la Terre ait été reçu (à cet effet, sont utilisés tous les hori¬ zons infrarouges). L'image du créneau k-1 est traitée pendant la réception de l'image du créneau k, comme il¬ lustré à la figure 1. Le traitement d'images lui-même né¬ cessite environ 6 à 8 mn. Ensuite, les données sont dis¬ tribuées aux usagers de METEOSAT ce qui implique que les usagers dans l'hémisphère Sud doivent attendre approxima- tivement 40 mn avant de recevoir les images qui couvrent leur zone, alors que les usagers européens doivent at¬ tendre jusqu'à 20 mn avant de recevoir l'image correspon¬ dante.

La présente invention s'est donc donné pour but de pourvoir à un satellite météorologique géosta- tionnaire qui reponde mieux à la nécessité de la pratique que les satellites du même type et visant au même but an¬ térieurement connus, notamment en ce qu'il effectue la correction de la déformation des images en temps réel, et ce sans le besoin de détecteur d'attitude additionnel so¬ phistiqué à bord du satellite.

La présente invention a pour objet un procédé de rectification en temps réel de données d'images détec- tées à distance par un instrument de mesure de radiations tel qu'utilisé à bord d'un satellite geostationnaire sta¬ bilisé en rotation, la rectification étant obtenue par détermination d'un certain nombre de paramètres du satel¬ lite comprenant les paramètres de pointage de l'instrument de mesure de radiations, lequel procédé com¬ prend les étapes suivantes :

1) prédiction des paramètres de rectification par analyse statistique à partir d'un certain nombre d'images préalables ; 2) lecture d'un nombre limité de lignes d'images brutes pour faciliter la position précise de l'horizon Sud ou de l'horizon Nord, et ce en fonction de la direction de balayage de 1'image reçue ;

3) détermination de la vitesse de rotation du satellite à partir d'un nombre limité de lignes d'images brutes et prédiction des variations de la vitesse de ro¬ tation pour la partie restante de 1' image ;

4) redéfinition de l'attitude et des para¬ mètres de pointage de 1'instrument de mesure de radia- tions utilisant l'information d'horizon obtenue en temps réel ;

5) calcul informatique de la position de l'autre horizon et du centre de l'image basé sur l'horizon premièrement détecté et sur les données géomé- triques ;

6) création d'une matrice de déformation fai¬ sant suite immédiatement à la détermination du premier horizon ; et

7) rectification en temps réel des données de 1'image reçue immédiatement après. Dans un mode d'exécution préféré du procédé l'instrument de mesure de radiations est constitué par un radiomètre, du type de ceux dont sont équipés les satel¬ lites européens METEOSAT. L'invention sera bien comprise à l'aide de la description qui va suivre faite en référence à la figure 2 du dessin annexé qui est un schéma en bloc-diagramme illustrant de façon schématique le procédé selon l'invention et qui montre de façon claire la différence qui existe par rapport au procédé de l'Art antérieur il¬ lustré lui aussi schématiquement à la figure 1.

La figure 2 illustre graphiquement la circula¬ tion de données impliquées dans la rectification en temps réel de données d'images, objet de la présente invention. La correspondance entre les figures 1 et 2 est très claire et ne requiert pas de commentaires particuliers : en fait, les blocs 1*, 2*, 3* et 4* de la figure 2 correspondent aux blocs 1, 2, 3 et 4 de la figure 1, alors que les cylindres 5 à 7 et 5 à 7 correspondent à des mémoires à disques.

Etant donné que la plupart des perturbations influençant les paramètres du satellite sont de nature périodique, il est logique de s'attendre à que la matrice de déformation correspondant au créneau réel soit basée sur une analyse statistique de l'information établie pour les créneaux précédents. Toutefois, la Demanderesse a déjà fait l'expérience du fait qu'un tel système ne four¬ nit pas la précision requise et c'est pourquoi elle a dé¬ veloppé le procédé perfectionné de la présente Demande, décrit ci-après en référence à un balayage de la terre du Sud vers le Nord et qui toutefois, bien entendu, est aussi applicable aux instruments de mesure de radiations pour satellites balayant la terre du Nord vers le Sud.

La rectification en temps réel, selon l'invention, des images de satellites météorologiques géostationnaires utilise un mode unique de combinaison des paramètres qui peuvent être évalués, ce mode de com¬ binaison étant basé sur l'analyse statistique des para¬ mètres établis pour les créneaux précédents, et sur les paramètres qui peuvent être déterminés à partir des don¬ nées en entrée, cette combinaison seulement, de para¬ mètres évalués (par extrapolation) et mesurés, fournissant la précision exigée.

Le procédé d'évaluation en temps réel de l'ensemble des paramètres ^ est effectué pour l'image k directement au début du créneau k. A ce moment, les hori¬ zons Est, Ouest ou Nord, tels qu'ils sont utilisés dans le procédé de rectification selon l'invention, ne sont pas encore disponibles. Toutefois, compte tenu de l'observation faite que la diffusion des images rectifiées n'est pas utile avant de disposer de quelques images contenant au moins l'horizon Sud, qui peut ainsi être utilisé pour améliorer la précision des paramètres, la position des lignes de l'horizon Sud peut être mesurée, -comme d'ailleurs pour l'horizon Est-Ouest-, avant le début du procédé de recti¬ fication. En plus, on a l'avantage de pouvoir déterminer le décalage de la vitesse de rotation du satellite autour de son axe à partir d'une information de minutage dans les premières lignes d'images.

Cela étant, les seuls paramètres qui ne sont pas connus à ce moment dans le temps sont 1'orientation affinée et les paramètres du satellite de même que le po¬ sitionnement de l'horizon Nord à l'intérieur de l'image. Cependant, l'orientation peut être affinée uniquement en utilisant l'horizon Sud. De plus, en admet¬ tant l'hypothèse que les autres paramètres de pointage restent constants sur un créneau, on peut adopter pour eux les mêmes valeurs que celles correspondant à 1'image précédente. Au demeurant, l'horizon Nord peut être obtenu à partir de la position du satellite, du rayon polaire de la Terre et des paramètres de pointage, à savoir de repé¬ rage du détecteur du satellite.

Il est possible de tenir compte des variations de la vitesse de rotation du satellite autour de son axe à l'aide d'une équation différentielle du premier ordre qui fait partie d'un algorithme approprié connu sous le nom de KALMAN FILTER (cf GELB, A., "Applied optimal esti¬ mation", MIT press,1974). Un prototype du système de création de la ma¬ trice de déformation et d'évaluation des paramètres a été développé par la Demanderesse et testé pendant plusieurs mois, ce qui a permis d'en établir une fois pour toutes la précision et la qualité. Les tests ont été effectués en tenant compte, comme élément de comparaison, de METEOSAT.

Le procédé, objet de la présente invention, peut être appliqué aux prévisions météorologiques en temps réel comme requis par le contrôle des trafics aé- riens et routiers.