Brevet FR3092565A1 — publié le 20200814

Description Titre de l'invention : Systèmes informatiques pour l'acquisition d'images satellitaires avec prise en compte de l'incertitude météo- rologique [pool L'invention concerne le domaine de la télédétection par satellite. Plus précisément, elle concerne des systèmes informatiques pour l'acquisition d'images satellitaires. Elle concerne aussi des dispositifs de calcul de plan mission pour satellite ainsi que des satellites d'observation de la Terre. [0002] L'espace circumterrestre contient nombre de satellites d'observation de la Terre qui permettent notamment de surveiller des paramètres du système Terre (océans, atmosphères, ressources naturelles, impacts des activités humaines, etc.), obtenir des cartographies ou fournir des informations qui permettent d'assurer la défense des États. Pour cela, ces systèmes se basent sur un ensemble de données d'imagerie satellitaire, c'est-à-dire un ensemble de mesures réalisées à l'aide d'instruments embarqués au niveau des satellites. [0003] Du fait de la banalisation de l'observation de la Terre, les demandes d'imagerie sa- tellitaire se multiplient de jour en jour. Or, certaines des ressources des satellites sont limitées, comme le nombre de prises de vues réalisables durant le survol d'une zone du globe ou encore les capacités de transmission des données vers le sol. On doit donc planifier l'utilisation de ces ressources tout en répondant aux demandes d'imagerie satellitaire. [0004] Dans le domaine de l'observation optique de la Terre, l'imagerie satellitaire dépend des conditions météorologiques. En effet, des nuages présents dans la zone d'acquisition peuvent masquer l'objet ou le territoire que l'on cherche à observer. Or, il est bien connu que les prévisions météorologiques présentent beaucoup d'incertitudes, quelle que soit la modélisation des phénomènes météorologiques. [0005] Afin de réduire ces incertitudes, il est connu de réaliser une prévision dite d'ensemble. Ainsi, plutôt que de simuler un seul jeu de données d'entrée avec le modèle, ce qui constitue une prévision déterministe, on calcule plusieurs jeux de données d'entrée déduit des incertitudes que l'on a estimées. On dispose alors d'un ensemble de prévisions que l'on espère représentatif de l'incertitude de la prévision déterministe. [0006] Toutefois, la prévision d'ensemble nécessite beaucoup de temps de calcul. [0007] Un tel caractère ne répond pas aux exigences toujours plus grandes de rapidité et de précision requises par les clients dans des systèmes de télédétection par satellite. [0008] La présente invention vise donc à pallier les inconvénients précités. 2 [0009] Pour cela, un premier aspect de l'invention concerne un système informatique pour l'acquisition d'images satellitaires. [0010] Enfin, un deuxième aspect de l'invention concerne un dispositif dc calcul dc plan mission pour satellite. [0011] Ainsi, l'invention concerne un système informatique pour l'acquisition d'images sa- tellitaires associées à au moins une zone prédéterminée de la Terre préalablement divisée selon une grille à mailles. Le système comprend : - au moins un satellite d'observation de la Terre comprenant, un système d'imagerie optique pour acquérir au moins une image satellitaire à au moins une date d'acquisition, un premier émetteur-récepteur satellite pour recevoir un signal de commande, et un premier processeur pour commander le système d'imagerie optique à partir du signal de commande, - au moins un dispositif de calcul de plan mission pour satellite comprenant, un deuxième émetteur-récepteur satellite pour émettre le signal de commande, une mémoire pour stocker, pour chaque maille, un historique de valeurs de variables météorologiques qui comprend des observations météorologiques et/ou des prévisions météorologiques qui ont été observées et/ou qui ont été prévues à une date qui précède la date d'acquisition, un deuxième processeur pour o déterminer, pour chaque maille, à partir de l'historique de valeurs de variables météorologiques, une pluralité de groupes de valeurs qui comprennent chacun, au moins : une première valeur de variable météorologique associée à une première date d'historique, une deuxième valeur de variable météorologique associée à une deuxième date d'historique postérieure ou antérieure à la première date d'historique, et une différence de temps d'historique entre la première date d'historique et la deuxième date d'historique, o calculer, pour tout ou partie de la zone prédéterminée de la Terre, à partir de la pluralité de groupes de valeurs, une fonction de probabilité conditionnelle qui décrit la probabilité d'observer une deuxième valeur de variable météorologique comprise dans un intervalle prédéterminé à une deuxième date d'observation sous la condition : d'observer une première valeur de variable météorologique à une première date d'observation qui est antérieure à la deuxième date d'observation, et/ou de prévoir une troisième valeur de variable météorologique à une date de prévision située autour de la deuxième d'observation, à savoir, à une date qui précède, suit et/ou égale à la deuxième date d'observation, et sous la condition supplémentaire d'observer une différence de temps prédé- 3 terminée entre : la première date d'observation et la deuxième date d'observation, et/ou la première date d'observation et la date de prévision, o sélectionner une pluralité de mailles de la grille, prévues pour être survolées par le satellite à la date d'acquisition, à partir des caractéristiques orbitales du satellite, de sorte à obtenir des mailles candidates, o obtenir, pour chaque maille candidate, au moins une observation météorologique courante associée à une date de planification qui précède la date d'acquisition et/ou au moins une prévision météorologique courante associée à une date située autour de la date d'acquisition, o calculer un plan mission du satellite à partir des mailles candidates et de la fonction de probabilité conditionnelle, et o insérer le plan mission dans le signal de commande. [0012] L'invention couvre également un dispositif de calcul de plan mission pour satellite. Le dispositif comprend : - un premier émetteur-récepteur satellite pour émettre un signal de commande qui commande un système d'imagerie optique d'un satellite d'observation de la Terre pour réaliser l'acquisition d'au moins une image satellitaire à au moins une date d'acquisition, - une mémoire pour stocker un historique de valeurs de variables météorologiques qui comprend des observations météorologiques associées à au moins une zone prédéterminée de la Terre préalablement divisée selon une grille à mailles, les observations météorologiques et/ou des prévisions météorologiques qui ont été observées et/ou qui ont été prévues à une date qui précède la date d'acquisition, - un deuxième processeur pour o déterminer, pour chaque maille, à partir de l'historique de valeurs de variables météorologiques, une pluralité de groupes de valeurs qui comprennent chacun, au moins : une première valeur de variable météorologique associée à une première date d'historique, une deuxième valeur de variable météorologique associée à une deuxième date d'historique postérieure ou antérieure à la première date d'historique, et une différence de temps d'historique entre la première date d'historique et la deuxième date d'historique, o calculer, pour tout ou partie de la zone prédéterminée de la Terre, à partir de la pluralité de groupes de valeurs, une fonction de probabilité conditionnelle qui décrit la probabilité d'observer une deuxième valeur de variable météorologique comprise dans un intervalle prédéterminé à une deuxième date d'observation sous la condition : d'observer une première valeur de variable météorologique à une première date d'observation qui est antérieure à la deuxième date d'observation, et/ou de prévoir une troisième valeur de variable météorologique à une date de prévision située autour de la deuxième d'observation, à savoir, à une date qui précède, suit et/ou égale à la deuxième date d'observation, et sous la condition supplémentaire d'observer une différence de temps prédéterminée entre : la première date d'observation et la deuxième date d'observation, et/ou la première date d'observation et la date de prévision, o sélectionner une pluralité de mailles de la grille, prévues pour être survolées par le satellite à la date d'acquisition, à partir des caractéristiques orbitales du satellite, de sorte à obtenir des mailles candidates, o obtenir, pour chaque maille candidate, au moins une observation météorologique courante associée à une date de planification qui précède la date d'acquisition et/ou au moins une prévision météorologique courante associée à une date située autour de la date d'acquisition, o calculer un plan mission du satellite à partir des mailles candidates et de la fonction de probabilité conditionnelle, et o insérer le plan mission dans le signal de commande. [0013] Dans une première mise en oeuvre, la pluralité de groupes de valeurs comprend, en outre, des valeurs de variables météorologiques qui sont associées à des mailles voisines de la maille courante. [0014] Dans une deuxième mise en oeuvre, pendant le calcul du plan mission du satellite, le deuxième processeur est en outre configuré pour : - sélectionner tout ou partie des mailles candidates à partir du résultat de la fonction de probabilité conditionnelle associée et des observations météorologiques courantes et/ou des prévisions météorologiques courantes, de sorte à obtenir des mailles sélectionnées, et - calculer le plan mission du satellite à partir des mailles sélectionnées. [0015] Dans un exemple de la deuxième mise en oeuvre, les mailles sélectionnées sont celles dont la fonction de probabilité conditionnelle sort, pour la date d'acquisition, une probabilité supérieure ou inférieure à une première probabilité prédéterminée. [0016] Dans une troisième mise en oeuvre, le deuxième processeur calcule la fonction de probabilité conditionnelle à partir de l'estimation d'une densité de probabilité de tout ou partie des groupes de valeurs. [0017] Dans une quatrième mise en oeuvre, le deuxième processeur calcule la fonction de probabilité conditionnelle grâce à un modèle d'apprentissage automatique préentraîné à partir de tout ou partie des groupes de valeurs. [0018] D'autres caractéristiques et avantages de l'invention seront mieux compris à la lecture de la description qui va suivre et en référence aux dessins annexés, donnés à titre illustratif et nullement limitatif. [0019] [fig.l]La figure 1 représente un système informatique selon l'invention. [0020] [fig.2]La figure 2 représente une autre vue du système informatique selon l' invention. [0021] [fig.3]La figure 3 représente une zone prédéterminée de la Terre. [0022] [fig.4]La figure 4 représente une grille à mailles qui couvre la zone prédéterminée de la Terre. [0023] Pour des raisons de clarté, les éléments représentés ne sont pas nécessairement re- présentés à la même échelle, les uns par rapport aux autres, sauf mention contraire. [0024] Dans le cadre de la description, on entend par système d'observation de la Terre, un système comprenant des constellations satellitaires telles que Landsat, Spot ou encore Pléiades. En outre, on précise que l'invention concerne les systèmes satellitaires optiques, pour lesquels l'acquisition d'images de la Terre dépend d'un éclairement solaire. Ces systèmes utilisent des capteurs passifs qui permettent d'enregistrer le rayonnement solaire réfléchi par le sol, dans diverses longueurs d'onde. [0025] Le principe général de l'invention est basé sur le constat que l'incertitude liée aux conditions météorologiques peut mener à l'élaboration de plans mission, pour satellite, peu performants qui couvrent des zones d'acquisitions présentant une faible probabilité d'obtenir des prises de vues acceptables à la date souhaitée d'acquisition. Ainsi, l'invention propose de caractériser cette incertitude et de l'utiliser pour améliorer l'élaboration des plans mission pour satellite. Dans l'invention, on construit un modèle des incertitudes à partir d'observations et/ou de prévisions météorologiques passées. Ensuite, on exploite le modèle, au moment de la planification, en présentant en entrée du modèle les observations et/ou les prévisions courantes. De cette manière, on peut identifier les zones d'acquisitions pour lesquelles il est plus probable d'obtenir des prises de vues acceptables à la date souhaitée d'acquisition. [0026] La figure 1 illustre selon l'invention, un système informatique 100 pour l'acquisition d'images satellitaires associées à au moins une zone prédéterminée de la Terre. [0027] Dans l'exemple de la figure 1, on considère que le système informatique 100 reçoit une requête en imagerie satellitaire similaire aux requêtes reçues et traitées dans les centres mission pour satellite. Toutefois, on envisage également que le système informatique 100 reçoive plusieurs requêtes concurrentes. [0028] Dans l'invention, la requête définit au moins une zone prédéterminée de la Terre dans laquelle l'acquisition en imagerie satellitaire doit être réalisée. Dans un exemple, on définit la zone prédéterminée de la Terre par un ou plusieurs polygones dont les sommets correspondent à des positions géographiques du globe terrestre. Dans une mise en oeuvre particulière, la requête est prévue pour être associée à au moins un 6 système d'observation de la Terre prédéterminé qui est configure pour survoler la zone prédéterminée de la Terre. [0029] La figure 3 illustre un exemple de zone prédéterminée de la Terre 10. [0030] Dans l'invention, la zone prédéterminée de la Terre 10 est divisée selon une grille composée de mailles. On entend par grille, par exemple, un quadrillage horizontal et vertical. On entend par maille, l'un quelconque des quadrilatères élémentaires qui composent la grille. Par exemple, chaque maille de la grille présente une forme carrée ou voisine d'un carré, par exemple, une forme rectangulaire. La figure 4 illustre un exemple de grille à mailles 20 superposée à la zone prédéterminée de la Terre 10. Dans un premier exemple, la grille selon l'invention se présente sous la forme d'un fichier image dans lequel chaque maille associée à la zone prédéterminée est géoréférencée. Dans un deuxième exemple, la grille selon l'invention se présente sous la forme d'un fichier texte dans lequel chaque ligne décrit au moins une maille de la grille avec des informations de géoréférencement associées à la zone prédéterminée. [0031] De retour à la figure 1, le système informatique 100 comprend au moins un satellite d'observation de la Terre 110 et au moins un dispositif de calcul de plan mission pour satellite 120. En outre, le système informatique 100 est agencé selon une architecture qui repose sur un réseau de communication satellitaire 130 de type connu. [0032] Dans l'exemple de la figure 2, le satellite d'observation de la Terre 110 comprend un système d'imagerie optique 111, un premier émetteur-récepteur satellite 112 et un premier processeur 113. [0033] Le système d'imagerie optique 1 1 1 est prévu pour acquérir au moins une image sa- tellitaire. Dans un exemple, le système d'imagerie optique Ill comprend au moins un instrument optique qui fonctionne selon un mode du type choisi parmi (en anglais) : push-broom », «push-frame », « step-and-store ». [0034] Dans l'invention, le système d'imagerie optique 111 est prévu pour acquérir l'image satellitaire à au moins une date d'acquisition. On entend par date, l'indication de l'heure, du jour, du mois et de l'année. [0035] Dans un exemple, on peut considérer que la date d'acquisition est prévue pour le « dd/MM/YYYY HH:mm» où dd correspond au jour, MM correspond au mois, YYYY correspond à l'année concernée, HH correspond à l'heure et mm correspond aux minutes. Dans un autre exemple, on peut considérer que l'acquisition satellitah-e est réalisée sur plusieurs jours aux dates suivantes «dd/MM/YYYY Hamm», « (dd-F2)/MM/YYYY (HH-FI):mm», « (dd+5)/MM/YYYY HH:(nun-F15)» et « (dd-F6/MM/YYYY (HH+2):(mm-3)» Dans l'exemple précédent, la date « (dd-F2)/MM/YYYY (HH-FI):mm» succède à la date «dd/MM/YYYY Hlimm» de deux jours plus une heure. La date «(dd-F5)/MM/YYYY HH:(mm-F15)» succède à la date «MI/MM/YYYY HH:rmn» de cinq jours moins quinze minutes. Et la date 7 (dd+6/MM/YYYY (HH-F2):(mm-3)» succède à la date «dd/MM/YYYY HH:mm » de six jours plus deux heures et moins trois minutes. Bien sûr, on peut envisager d'autres intervalles de temps d'acquisition plus court ou plus long selon les besoins du client. [0036] Le premier émetteur-récepteur satellite 112 est prévu pour recevoir un signal de commande. Dans un exemple, le signal de commande est un signal radio satellitaire qui utilise un protocole de communication satellitaire de type connu. [0037] Le premier processeur 113 est prévu pour commander le système d'imagerie optique à partir du signal de commande. [0038] Dans l'exemple de la figure 2, le dispositif de calcul de plan mission pour satellite 120 comprend un deuxième émetteur-récepteur satellite 121, une mémoire 122 et un deuxième processeur 123. [0039] Le deuxième émetteur-récepteur satellite 121 est prévu pour émettre le signal de commande. [0040] La mémoire 122 est prévue pour stocker, pour chaque maille de la grille 20, un historique de valeurs de variables météorologiques qui comprend des observations météorologiques et/ou des prévisions météorologiques passées que l'on a observées et/ou prévues à unc date qui précède la date d'acquisition. [0041] Dans l'invention, les valeurs stockées dans la mémoire 122 sont associées à des variables météorologiques connues telles que la nébulosité, la température, le niveau de précipitation, les caractéristiques du vent, l'humidité relative ou le niveau d'ensoleillement. Dans un exemple, l'historique de valeurs de variables météorologiques peut être contenu dans un fichier de données ERA-Interim produit par le CEPMMT (Centre Européen pour les Prévisions Métérologiques à Moyen Terme, ou ECMWF en anglais pour European Centre for Medium-Range Weather Forecasts). Un tel fichier, selon l'invention, peut comprendre au moins une valeur de variable météorologique pour chacune des mailles de la grille 20, et ce, sur une ou plusieurs dates d'historique au sein d'une ou plusieurs années antérieures à la date d'acquisition. [0042] Le deuxième processeur 123 est prévu pour traiter l'historique de valeurs de variables météorologiques pour former le signal de commande. [0043] Tout d'abord, le deuxième processeur 123 est prévu pour calculer pour chaque maille de la grille 20, à partir de l'historique de valeurs de variables météorologiques, une pluralité de groupes de valeurs qui comprennent chacun, au moins : - une première valeur de variable météorologique associée à une première date d'historique, - une deuxième valeur de variable météorologique associée à une deuxième date d'historique postérieure ou antérieure à la première date d'historique, et - une différence de temps d'historique entre la première date d'historique et la deuxième date d'historique. 8 [0044] Ainsi, pour chaque maille de la grille 20, le deuxième processeur 123 calcule au moins un triplet de valeurs. En pratique, le deuxième processeur 123 identifie au moins une première valeur de variable météorologique associée à une première date d'historique. Ensuite, le deuxième processeur 123 identifie au moins une deuxième valeur de variable météorologique associée à une deuxième date d'historique qui est postérieure ou antérieure à la première date d'historique. Puis, le deuxième processeur 123 calcule la différence de temps entre la première date d'historique et la deuxième date d'observation. Par la suite, le deuxième processeur 123 associe les première et deuxième valeurs météorologiques avec la différence de temps au sein d'un groupe de valeurs. Enfin, le deuxième processeur 123 sauvegarde le groupe de valeurs dans la mémoire 122. [0045] Prenons un exemple pour illustrer cet aspect de l'invention. Dans cet exemple, on suppose que l'une des mailles de la grille 20 est associée à un historique de valeurs de néhulosité exprimé en pourcentage de ciel clair. De manière classique, on entend par néhulosité la part de ciel couvert par tout ou partie des nuages au niveau de ladite maille de la grille 20. Comme indiqué plus haut, on peut obtenir les valeurs de nébulosités par observation et/ou par prévision. Dans cet exemple, la maille en question de la grille 20 comprend les valeurs de nébulosité suivantes : 100 % (date 1), 97 % (date 2), 99 % (date 3), 80 % (date 4), 77 % (date 5), 79 % (date 6), pour lesquelles les dates date 1, date 2, date 3, date 4, date 5 et date 6 peuvent être exprimées sous la même forme que celle de la date d'acquisition, comme décrit plus haut. De retour à l'exemple, le deuxième processeur 123 peut déterminer les groupes de valeurs suivants en considérant date_1 comme la première date d'historique : {100 %, 97 %, date_2-date_1}, {100 %, 99 %, date_3-date_1 }, 1100 %, 80 %, date_4-date_11, {100 %, 77 %, date_5-date_1 } et 1100 %, 79 %, date_6-date_11. Lorsque l'on considère date_3 comme la première date d'historique, le deuxième processeur 123 peut déterminer les groupes de valeurs suivants : {99 %, 80 %, date_4-date_3}, {99 %, 77 %, date_5-date_3} et {99 %, 79 %, date_6-date_3 J. [0046] Toutefois, le deuxième processeur 123 est également configuré pour calculer des groupes de valeurs qui comprennent plus de trois valeurs. Dans ce cas, le groupe de valeurs associé à une maille courante de la grille 20 peut comprendre, en outre, des valeurs de variables météorologiques qui sont associées à des mailles voisines de la maille courante de la grille 20. On entend par maille voisine, une maille de la grille 20 qui est située sur la même ligne, la même colonne et/ou autour de la maille courante. Comme indiqué plus haut, on peut obtenir les valeurs de variables météorologiques associées aux mailles voisines par observation et/ou par prévision. [0047] De retour à la figure 2, le deuxième processeur 123 est prévu pour calculer une fonction de probabilité conditionnelle spécifique à tout ou partie de la zone prédé- 9 terminée de la Terre 10, à partir de la pluralité de groupes de valeurs. La fonction de probabilité conditionnelle selon l'invention permet de relier des observations météorologiques et/ou des provisions météorologiques, observées et/ou réalisées à des dates différentes, et ce, en fonction d'une différence de temps entre les dates d'observation et/ou les dates de prévision. [0048] Dans l'invention, la fonction de probabilité conditionnelle décrit la probabilité d'observer une deuxième valeur de variable météorologique comprise dans un intervalle prédéterminé à une deuxième date d'observation sous la condition : - d'observer une première valeur de variable météorologique à une première date d'observation qui est antérieure à la deuxième date d'observation, et/ou - de prévoir une troisième valeur de variable météorologique à une date de prévision située autour de la deuxième d'observation, à savoir, à une date qui précède, suit et/ou égale à la deuxième date d'observation. En outre, le deuxième processeur 123 calcule la fonction de probabilité conditionnelle sous la condition supplémentaire d'observer une différence de temps prédéterminée entre : - la première date d'observation et la deuxième date d'observation, et/ou - la première date d'observation et la date de prévision. [0049] Dans un exemple, la fonction de probabilité conditionnelle selon l'invention permet de répondre aux questions du type : - quelle est la probabilité d'observer dans une maille de la grille 20, une valeur de nébulosité supérieure à 82 % dans trois heures, sachant que j'y observe maintenant la valeur 77 %? Dans ce cas, l'intervalle prédéterminé peut être compris entre 82% et 100%. - quelle est la probabilité d'observer dans une maille de la grille 20, une valeur de nébulosité supérieure à 70 % et inférieure à 80% dans cinq heures, sachant qu'une prévision réalisée il y a deux heures y prévoyait la valeur 10 %? Dans ce cas, l'intervalle prédéterminé peut être compris entre 70% et 80%. - quelle est la probabilité d'observer dans une maille de la grille 20, une valeur de nébulosité supérieure à 60 % dans quatre heures, sachant qu'une prévision réalisée pour dans quatre heures y prévoit la valeur 35 % ? Dans ce cas, l'intervalle prédéterminé peut être compris entre 60% et 100%. - quelle est la probabilité d'observer dans une maille de la grille 20, une valeur de nébulosité inférieure à 20 % dans deux heures, sachant que j'y observe maintenant la valeur 30 % et qu'une prévision réalisée pour dans deux heures y prévoit la valeur 10 %? Dans ce cas, l'intervalle prédéterminé peut être compris entre 0% et 20%. [0050] Dans une mise en oeuvre particulière, dans laquelle les groupes de valeurs com- prennent plus de trois valeurs, la fonction de probabilité conditionnelle décrit la pro- habilité d'observer une deuxième valeur de variable météorologique comprise dans un intervalle prédéterminé à une deuxième date d'observation sous la condition supplémentaire : - d'observer au moins une quatrième valeur de variable météorologique dans une maille voisine de la maille courante, à une cinquième date d'observation qui est antérieure à la deuxième date d'observation, et/ou - de prévoir au moins une cinquième valeur de variable météorologique dans la maille voisine de la maille courante, à une date de prévision située autour de la deuxième d'observation, à savoir, à une date qui précède, suit et/ou égale à la deuxième date d'observation. [0051] Dans un exemple de cette mise en oeuvre particulière, la fonction de probabilité conditionnelle selon l'invention permet de répondre aux questions du type : - quelle est la probabilité d'observer dans une maille de la grille 20, une valeur de nébulosité supérieure à 82 % dans trois heures, sachant que j'y observe maintenant la valeur 77 % et que j'observe maintenant la valeur 52 % dans une maille voisine ? Dans ce cas, l'intervalle prédéterminé peut être compris entre 82% et 100%. - quelle est la probabilité d'observer dans une maille de la grille 20, une valeur de nébulosité supérieure à 70 % et inférieure à 85% dans cinq heures, sachant qu'une prévision réalisée il y a deux heures y prévoyait la valeur 10 % et qu'une prévision réalisée il y a deux heures prévoyait la valeur 10 % dans une maille voisine ? Dans ce cas, l'intervalle prédéterminé peut être compris entre 70% et 85%. - quelle est la probabilité d'observer dans une maille de la grille 20, une valeur de nébulosité supérieure à 60 % dans quatre heures, sachant qu'une prévision réalisée pour dans quatre heures y prévoit la valeur 35 % et qu'une prévision réalisée pour dans quatre heures y prévoit la valeur 35 % dans une maille voisine ? Dans ce cas, l'intervalle prédéterminé peut être compris entre 60% et 100%. - quelle est la probabilité d'observer dans une maille de la grille 20, une valeur de nébulosité inférieure à 20 % dans deux heures, sachant que j'y observe maintenant la valeur 30 % et qu'une prévision réalisée pour dans deux heures y prévoit la valeur 10 % et que j'observe maintenant la valeur 27 % dans une maille voisine ? Dans ce cas, l'intervalle prédéterminé peut être compris entre 0% et 20%. - quelle est la probabilité d'observer dans une maille de la grille 20, une valeur de nébulosité inférieure à 20 % et supérieure à 10 % dans deux heures, sachant que j'y observe maintenant la valeur 30 % et qu'une prévision réalisée pour dans deux hem-es y prévoit la valeur 10 % et qu'une prévision réalisée pour dans deux heures prévoit la valeur 15 % dans une maille voisine ? Dans ce cas, l'intervalle prédéterminé peut être compris entre 10% et 20%. [0052] Dans un exemple de l'invention, le deuxième processeur 123 calcule la fonction de 11 probabilité conditionnelle à partir de l'estimation d'une densité de probabilité de tout ou partie des groupes de valeurs. On entend par densité de probabilité, une fonction qui pondère mathématiquement l'éventail des valeurs possibles que peut prendre un paramètre. Dans cet exemple, le deuxième processeur 123 estime la densité de probabilité en utilisant des techniques connues, telles que l'estimation par les histogrammes, estimation par les histogrammes mobiles ou l'estimation à noyau (par ex., la méthode de Parzen-Rosenblatt). [0053] Dans un autre exemple de l'invention, le deuxième processeur 123 calcule la fonction de probabilité conditionnelle grâce à un modèle d'apprentissage automatique préentraîné à partir de tout ou partie des groupes de valeurs. Dans cet exemple, on entraîne un modèle d'apprentissage automatique à partir d'un algorithme de classification supervisée. Ainsi entraîné, le modèle d'apprentissage automatique précntraîné pourra prédire, pour chaque maille de la grille 20, une classe représentative d'une valeur de probabilité conditionnelle associée à la deuxième valeur de variable météorologique à la date deuxième date d'observation en fonction, au moins : - de l'observation météorologique effectuée à la première date d'observation et/ou de la prévision effectuée autour de la deuxième date d'observation, et - de la différence de temps entre la première date d'observation et la deuxième date d'observation et/ou la différence de temps entre la première date d'observation et la date de prévision. [0054] Dans ce cas, la phase d'entraînement nécessite de former une pluralité de vecteurs qui comprennent chacun un groupe de valeurs tel que décrit plus haut. Ensuite, on calcule, pour chaque groupe de valeur, la valeur de probabilité conditionnelle, tel que décrit plus haut. Dans un exemple, la probabilité conditionnelle est calculée à partir d'une densité de probabilité de tout ou partie des groupes de valeurs, comme décrit plus haut. Puis, on définit une classe indicative de la valeur de probabilité conditionnelle calculée. Enfin, on entraîne le modèle d'apprentissage automatique en présentant respectivement les vecteurs et les classes associées, à l'entrée et à la sortie du modèle d'apprentissage automatique, de sorte à obtenir un modèle d'apprentissage automatique préentraîné. Dans une mise en oeuvre particulière, on constitue le modèle d'apprentissage automatique à partir d'un réseau de neurones artificiels de type perceptron multicouche ou un réseau de neurones de type convolutif. [0055] De retour à la figure 2, le deuxième processeur 123 sélectionne une pluralité de mailles de la grille 20, que le satellite 110 prévoit de survoler à la date d'acquisition, à partir des caractéristiques orbitales du satellite à la date d'acquisition, de sorte à obtenir des mailles candidates. [0056] On entend par caractéristiques orbitales d'un satellite, les caractéristiques conven- tionnelles telles que l'altitude, le nombre de révolutions par jour ou l'inclinaison de 12 l'orbite. On peut également définir d'autres paramètres du système conjointement aux caractéristiques orbitales d'un satellite, comme la liste des stations sols permettant la redescente des images réalisées par le satellite, la résolution du système d'imagerie, les bandes spectrales supportées par le système d'imagerie du satellite, la fauchée de l'instrument optique ou encore les limites d'angle de dépointage autorisées pour la prise de vues. [0057] Puis, le deuxième processeur 123 obtient, pour chaque maille candidate, au moins une observation météorologique courante associée à une date de planification qui précède la date d'acquisition et/ou au moins une prévision météorologique courante associée à une date située autour de la date d'acquisition, à savoir, qui précède, suit ou est égal à la date d'acquisition. Dans un exemple, la date de planification correspond à une date de commencement de planification d'une mission d'acquisition. En pratique, cette date correspond à la date à laquelle le centre mission pour satellite reçoit ou prend en charge la requête d'imagerie satellitaire formulée par les clients. Dans l'invention, on obtient l'observation météorologique courante et/ou la prévision météorologique courante à partir des méthodes connues. [0058] Dans une mise en oeuvre particulière, comme indiqué plus haut, le deuxième processeur 123 est également prévu pour obtenir, pour chaque maille voisine de la maille candidate courante, des observations et/ou les prévisions courantes. [0059] Par la suite, le deuxième processeur 123 est prévu pour calculer un plan mission du satellite 110 à partir des mailles candidates et de la fonction de probabilité conditionnelle. [0060] On entend par plan mission (en anglais, mission plan), le plan de programmation du satellite considéré que l'on télécharge/envoie au satellite via le centre de contrôle et qui permet d'effectuer les prises de vue de la zone prédéterminée de la Terre 10 suivant un ordre et des caractéristiques de prise de vue fournies dans le plan. Il est connu que le plan mission prend en compte des contraintes cinématiques liées à la plateforme et intègre également les plans de vidages sur les différentes stations de réception. On peut obtenir un tel plan mission par un générateur de plan mission de type connu, par exemple, comme décrit dans la demande FR3039728. [0061] Enfin, de retour à la figure 2, le deuxième processeur 123 est prévu pour insérer le plan mission dans le signal de commande. [0062] Dans une mise en oeuvre de l'invention, pendant le calcul du plan mission du satellite, le deuxième processeur 123 est en outre prévu pour sélectionner tout ou partie des mailles candidates à partir du résultat de la fonction de probabilité conditionnelle associée et des observations météorologiques courantes et/ou des prévisions météorologiques courantes, de sorte à obtenir des mailles sélectionnées. Ensuite, le deuxième processeur 123 calcule le plan mission du satellite à partir des mailles sélectionnées. 13 [0063] Dans cette mise en oeuvre, on peut écarter de l'acquisition satellitai rc les mailles candidates pour lesquelles il est peu probable d'obtenir des prises de vues acceptables. [0064] Par exemple, les mailles sélectionnées sont celles dont la fonction de probabilité conditionnelle sort, pour la date d'acquisition, une probabilité supérieure ou inférieure à une première probabilité prédéterminée. [0065] Toutefois, d'autres critères peuvent en outre être utilisés pour réaliser la sélection. Par exemple, on peut utiliser une information de priorité de la requête en imagerie satellitaire ou une information de surface des mailles candidates. [0066] Dans un mode particulier de réalisation de l'invention, des instructions de programmes d'ordinateur mettent en oeuvre les différentes opérations réalisées par les premier et deuxième processeurs 113, 123. Par conséquent, l'invention vise aussi un programme avec un code de programme d'ordinateur fixé sur un support de stockage non transitoire, et qui est susceptible d'exécuter les différentes opérations réalisées par les premier et deuxième processeurs 113, 123 lorsque le programme d'ordinateur est chargé dans l'ordinateur ou exécuté dans l'ordinateur. [0067] On a décrit et illustré l'invention dans la présente description détaillée et dans les figures. Toutefois, la présente invention ne se limite pas aux formes de réalisation présentées. Ainsi, d'autres variantes et modes de réalisation peuvent être déduits et mis en oeuvre par la personne du métier à la lecture de la présente description et des figures annexées. 14 [Revendication 1]