FR 3151452 A1 20250124 FR 2307360 A 20230717 Il est décrit un procédé de communication par liaison optique entre au moins un satellite (10) et une pluralité de stations émettrices et/ou réceptrices au sol (20), comprenant : l’obtention de données de couvertures nuageuses des stations au sol, et de données relatives à une orbite du satellite, la détermination, pour chaque station au sol (20), de prévisions sur la possibilité d’établir une liaison optique ininterrompue entre le satellite et la station au sol, pour au moins deux horizons de temps déterminés, à court terme et à long terme, et la sélection d’au moins une des stations au sol (20) pour l’établissement d’une communication par liaison optique montante ou descendante avec le satellite, en fonction de la durée prédite la plus longue pour laquelle il est possible d’établir une liaison optique ininterrompue entre la station et le satellite. Figure de l’abrégé : Figure 1 Procédé de gestion d’un système satellitaire de télécommunication par lien optique aidée par des modèles météorologiques pour une continuité de service. La présente divulgation concerne un procédé d’établissement de communications par liaison optique entre au moins un satellite et un ensemble de stations au sol. Il est connu de transférer des données entre une ou plusieurs stations au sol et un satellite, en vue de retransmettre ces données à des récepteurs au sol, par une liaison radiofréquence. Il peut toutefois être avantageux de remplacer une liaison radiofréquence par une liaison optique, qui présente notamment l’avantage d’un débit bien plus élevé, et également d’une précision plus importante dans la direction de la liaison, ce qui garantit une plus grande sécurité. En effet, il est difficile d’intercepter une liaison optique sans perturber cette liaison, et donc rendre cette intrusion détectable. En revanche, l’utilisation d’un lien optique entre un satellite et une station au sol présente l’inconvénient de risques de coupures accrus du fait du passage d’un nuage sur la ligne de visée entre la station et le satellite. Ce risque de coupure est significativement plus élevé que sur des liaisons radiofréquences. En effet, dans les liaisons radiofréquences, la pluie peut dégrader le signal mais ne le coupe pas entièrement, à la différence du passage d’un nuage pour un lien optique. Le passage d’un nuage étant en outre bien plus fréquent, ce risque de coupure est significativement plus élevé que sur des liaisons radiofréquences. Par conséquent, pour pouvoir assurer un service continu, il est nécessaire d’utiliser plusieurs stations au sol en parallèle et de déterminer à chaque instant la station qui sera choisie pour permettre d’établir une liaison optique avec le satellite. On connaît du document EP 2873172 un système de communication optique comprenant une constellation de satellites, chaque satellite comprenant des télescopes optiques inter-satellites et des télescopes optiques de liaison montante/descendante pour des communications avec des sites terrestres, et dans lequel le système de communication optique est configuré pour que lorsque la constellation de satellites passe devant un site terrestre donné, un ou plusieurs des télescopes de liaison montante/descendante d’un satellite suivent au moins deux télescopes optiques terrestres du site terrestre donné et commutent entre des télescopes optiques terrestres de manière à envoyer des données à celui ayant la visibilité directe la plus claire du satellite donné. Ce procédé permet donc de tenir compte de la nébulosité pour l’établissement d’une communication entre un site terrestre et un satellite. On connaît également du document AU2019363341 un système de gestion d’une station optique terrestre comprenant un dispositif de surveillance d’une couverture nuageuse basse, un dispositif de détermination d’une couverture nuageuse haute à partir d’un satellite météorologique, et un dispositif qui superpose les données sur la couverture nuageuse basse et la couverture nuageuse haute et qui contrôle la station optique terrestre à partir de ces informations. Résumé La présente divulgation vient améliorer la situation. En particulier, la présente invention vise à améliorer la qualité de service en choisissant de manière plus efficace quelle station au sol sera chargée d’établir la communication optique avec le satellite. En particulier, un but de la présente divulgation est de permettre le maintien en continu d’une communication par liaison optique entre une pluralité de stations au sol et un satellite. Il est proposé un procédé de communication par liaison optique entre au moins un satellite et une pluralité de stations émettrices et/ou réceptrices au sol, le procédé étant mis en œuvre par un calculateur et comprenant : l’obtention de données actualisées de couvertures nuageuses propres à chaque station au sol, l’obtention de données relatives à une orbite du satellite, la détermination, pour chaque station au sol, compte tenu des données relatives à l’orbite du satellite et par l’application d’au moins un modèle de risque, de prévisions sur la possibilité d’établir une liaison optique ininterrompue entre le satellite et la station au sol, en fonction des données actualisées de couverture nuageuse propre à chaque station au sol, pour au moins deux horizons de temps déterminés, à court terme et à long terme, ledit modèle de risque étant établi sur la base de données d’historiques de couvertures nuageuses des stations au sol, et la sélection d’au moins une des stations au sol pour l’établissement d’une communication par liaison optique montante ou descendante avec le satellite, cette station au sol étant sélectionnée sur la base de la durée prédite la plus longue pour laquelle il est possible d’établir une liaison optique ininterrompue entre la station et le satellite, à partir d’un instant déterminé d’établissement de la liaison optique, en fonction des prévisions sur chacun des deux horizons de temps. Dans des modes de réalisation, ledit modèle de risque fournit des prévisions sur la possibilité d’établir la liaison optique ininterrompue, associées à un indice de confiance, la sélection de ladite station étant réalisée également en fonction d’une durée d’établissement du futur lien et de l’indice de confiance, un poids plus ou moins faible étant respectivement attribué à la durée d’établissement du futur lien selon que l’indice de confiance est respectivement plus ou moins important. Dans des modes de réalisation, l’orbite du satellite est défilante et impose des créneaux temporels de visibilité propres à chaque station au sol, où les communications sont envisageables, ledit modèle de risque étant appliqué pour un ensemble de stations au sol dont les créneaux temporels de visibilité se recouvrent. Dans des modes de réalisation, la sélection comprend plusieurs stations au sol pour l’établissement de communications par liaison optique descendantes avec le satellite, ces stations au sol étant sélectionnées sur la base de durées prédites cumulées pour lesquelles il est possible d’établir des liaisons optiques ininterrompues entre les stations au sol et le satellite, à partir d’instants déterminés d’établissement de liaisons optiques propres à chaque station au sol et en fonction des prévisions, propres à chaque station au sol, sur deux horizons de temps. Dans des modes de réalisation, l’étape de détermination de prévisions utilise deux modèles de risque distincts, du type réseau de neurones, dont l’un est entrainé sur la base de données d’historiques de couvertures nuageuses locales des stations au sol pour le calcul de prévisions à court terme et l’autre est entraîné sur la base de données d’historiques de couvertures nuageuses globales des stations au sol pour le calcul de prévisions à long terme. Dans des modes de réalisation, le procédé comprend en outre une étape de mise à jour des données d’historiques de couvertures nuageuses par une intégration des données actualisées de couverture nuageuses et une étape de mise à jour dudit modèle de risque. Dans des modes de réalisation, la sélection d’au moins une station au sol comprend la détermination d’une séquence de communication comprenant une suite d’au moins deux stations au sol utilisées successivement pour l’établissement de la communication optique montante ou descendante avec le satellite, la séquence étant déterminée de façon à satisfaire au moins un premier critère portant sur l’absence de coupure de la communication optique au cours de la séquence. Dans des modes de réalisation, la séquence est déterminée de façon à satisfaire en outre un deuxième critère portant sur le nombre de changements de stations au cours de la séquence. Dans des modes de réalisation, les données actualisées de couverture nuageuse obtenues comprennent des données acquises par des capteurs de nébulosité au sol, des données d’observation météorologique obtenues par au moins un satellite, et des données obtenues par la mise en œuvre de modèles numériques de prédiction. Dans des modes de réalisation, la pluralité de stations au sol comprend une station dite courante en communication par liaison optique avec le satellite lors de la mise en œuvre du procédé, et l’étape de détermination comprend : une première détermination d’une possibilité d’établir une liaison optique entre le satellite et la station courante pendant au moins un horizon de temps déterminé, et ; si une liaison optique entre le satellite et la station courante est prédite comme impossible à l’horizon de temps: la mise en œuvre d’une deuxième détermination d’une possibilité d’établir une liaison optique entre le satellite et les autres stations au sol pendant au moins deux horizons de temps déterminés, cette deuxième détermination étant suivie de la sélection de l’une parmi les autres stations au sol pour l’établissement de la communication par liaison optique. Dans des modes de réalisation, l’instant d’établissement de la liaison optique avec la station sélectionnée est déterminé comme un instant antérieur à l’horizon de temps pour lequel une liaison optique est prédite comme impossible entre le satellite et la station courante. Dans des modes de réalisation, au moins l’un des horizons de temps considérés lors de la deuxième détermination se termine ultérieurement à l’horizon de temps considéré pour la première détermination, et le procédé comprenant la mise en œuvre d’un premier modèle entrainé pour la première détermination relative à la station courante, et d’un deuxième modèle entrainé pour la deuxième détermination relative aux autres stations. Dans des modes de réalisation, le procédé comprend : l’obtention de données d’observation de nébulosité acquises par un capteur pour la station courante, et la mise en œuvre de la première détermination sur la base desdites données d’observation de nébulosité, et l’obtention de données d’observation météorologiques obtenues par au moins un satellite et couvrant une région incluant les autres stations au sol, et la mise en œuvre de la deuxième détermination sur la base desdites données d’observation météorologiques. Dans des modes de réalisation, les données relatives à la couverture nuageuse au voisinage d’une station au sol comprennent une série temporelle de données relatives à la couverture nuageuse au voisinage d’une station au sol acquises à une fréquence déterminée. Selon un autre objet, il est décrit un produit programme d’ordinateur, comprenant des instructions de code pour la mise en œuvre du procédé selon la description qui précède, lorsqu’il est mis en œuvre par un calculateur. Selon un autre objet, il est décrit un système de communication par liaison optique entre un satellite et une pluralité de stations émettrices et/ou réceptrices au sol, le système comprenant au moins un terminal de contrôle comprenant un calculateur et une mémoire, le terminal de contrôle étant adapté pour communiquer avec les stations au sol, caractérisé en ce qu’il est configuré pour mettre en œuvre le procédé selon la description qui précède. Le procédé décrit ci-avant permet avantageusement de sélectionner une station émettrice ou réceptrice au sol pour l’établissement d’une communication par liaison optique, respectivement montante ou descendante, avec un satellite, la sélection de la station au sol étant réalisée sur la base d'estimations, à court terme et à long terme, de probabilités de coupure de la liaison optique causée par une couverture nuageuse entre l’ensemble des stations au sol considérées et le satellite. De façon avantageuse, il est possible, non seulement de maintenir une liaison de communication continue entre le satellite et une des stations au sol, en sélectionnant efficacement une station pour laquelle aucune couverture nuageuse n’est prévue, tout en minimisant le nombre de transitions d’une station à l’autre. Un autre avantage de la présente invention est que le calcul d’une disponibilité d’une station pour un satellite peut être réalisé pour différents horizons de temps, à long terme et à court terme et permet ainsi d’ajuster avec précision les transitions d’une station à l’autre ainsi que le nombre de transitions nécessaires. Le calcul de disponibilité est ici une estimation de la probabilité de disponibilité à une échéance donnée. Dans des modes de réalisation, la prédiction peut avantageusement être réalisée sur la base de données locales d’observation du ciel par un capteur au sol et sur la base de données à plus grande échelle, acquises par exemple par satellite, ce qui permet de tenir compte de deux échelles d’observations différentes et ainsi d’obtenir deux prédictions à long terme et à court terme avec leur précisions respectives. Alternativement, des données de prédiction de nébulosité, obtenues par des modèles numériques, peuvent également être utilisées pour déterminer la possibilité d’établir la communication optique entre le satellite et les stations au sol. D’autres caractéristiques, détails et avantages apparaîtront à la lecture de la description détaillée ci-après, et à l’analyse des dessins annexés, sur lesquels : Fig. 1 représente schématiquement un système de communication optique selon un exemple de réalisation. Fig. 2 représente schématiquement les principales étapes d’un procédé de communication selon un exemple de réalisation. Fig. 3 représente schématiquement un exemple de sélection d’une séquence de stations au sol pour l’établissement d’une communication optique. Fig. 4 représente schématiquement un exemple de mise en œuvre de la sélection d’une séquence de stations au sol utilisant deux modèles entrainés. Fig. 5 représente un exemple de donnée de couverture nuageuse acquise par un capteur de nébulosité. Fig. 6 représente schématiquement un exemple de modèle entrainé pour prédire la possibilité d’établir une liaison optique entre le satellite et une station au sol en fonction de données de couvertures nuageuses. Fig. 7 représente un exemple de structure de réseau de neurones pouvant être utilisé pour prédire la possibilité d’établir une liaison optique ininterrompue en fonction de données de couverture nuageuse. Fig. 8 illustre des résultats obtenus par la mise en œuvre du procédé selon un exemple de réalisation. Il est maintenant fait référence à la , qui représente un système 1 de communication par liaison optique entre au moins un satellite 10 et une pluralité de stations au sol 20. Les stations au sol 20 peuvent être des stations émettrices, la communication étant alors une communication dite « montante » dans laquelle des données sont envoyées par les stations au sol vers un satellite, ou des stations réceptrices, la communication étant alors dite « descendante », dans laquelle des données sont envoyées par un satellite aux stations au sol. Chaque station au sol peut être à la fois émettrice et réceptrice, les rôles d’émetteur et de récepteur étant variables en fonction des besoins de communication. Dans des exemples de réalisation, la communication par liaison optique peut comprendre une communication par liaison optique montante entre une ou plusieurs premières stations émettrices au sol et le satellite, puis une communication par liaison optique descendante entre le satellite et une ou plusieurs deuxièmes stations réceptrices au sol. En variante, et comme représenté sur la , la communication peut comprendre une communication par liaison optique montante entre des stations émettrices au sol et le satellite, puis une communication descendante par un autre moyen, typiquement par communication radiofréquence, entre le satellite et des stations au sol ou des terminaux utilisateurs 40. Selon encore une autre variante, on peut avoir une communication montante radiofréquence, et une communication descendante par liaison optique entre le satellite et les stations au sol. Les stations au sol 20 sont avantageusement en communication entre elles par une liaison terrestre, par exemple une communication filaire de type fibre optique. Le procédé selon l’invention peut être mis en œuvre pour une communication par liaison optique avec un ou plusieurs satellites 10 géostationnaires ou via des constellations défilantes. Les satellites utilisés pour établir le lien optique de connexion peuvent avoir une orbite du type GEO (géostationnaire), MEO (medium earth orbit) ou LEO (low earth orbit). Dans le cas d’un satellite non-géostationnaire les fenêtres de visibilité de chaque station au sol par chaque satellite sont fonction des positions des satellites en orbites et précisément connues et déterminées à l’avance. Un satellite défilant peut notamment contrôler son pointage en continu vers une station au sol, pendant une période déterminée de visibilité entre le satellite et cette station au sol. Le pointage d’un satellite vers une station au sol est par exemple réalisé par un contrôle de l’attitude du satellite. On peut aussi envisager que la tête de communication optique à bord du satellite s’oriente par rapport au corps du satellite, pour pointer vers une station au sol. Le pointage d’un satellite est modifié successivement pour pointer successivement vers différentes stations au sol. La tête de communication optique comprend par exemple un télescope optique. En tout état de cause, les informations relatives à l’orbite du ou des satellites, et donc à la position relative du ou des satellites par rapport aux stations au sol, sont connues. Le système de communication comprend par exemple en outre un terminal de contrôle 30, comprenant au moins un calculateur 31 et une mémoire 32, le terminal de contrôle 30 étant adapté pour communiquer avec les stations au sol et pour mettre en œuvre le procédé décrit ci-après. Le calculateur 31 peut comprendre un ou plusieurs processeurs, microprocesseurs, microcontrôleurs, processeur graphique, etc. La mémoire 32 comprend par exemple une mémoire non volatile permettant de stocker des instructions de code qui sont exécutées par le calculateur pour la mise en œuvre du procédé. De plus, la mémoire stocke par exemple au moins un modèle préalablement entrainé pour prédire la possibilité d’établir une liaison optique ininterrompue entre le satellite et au moins une station au sol, comme décrit plus en détails ci-après. Dans des exemples de réalisation, le système 1 de communication par liaison optique comprend en outre au moins un capteur de nébulosité 50 positionné au voisinage d’au moins une des stations au sol 20. Un capteur de nébulosité 50 est un capteur au sol, orienté vers le ciel, et adapté pour acquérir des données relatives à la couverture nuageuse instantanée, c’est-à-dire au moment de l’observation, dans la portion du ciel se trouvant dans le champ de vision du capteur. Par exemple, un tel capteur peut comprendre un capteur infrarouge associé à un calculateur configuré pour déterminer, à partir des informations du capteur infrarouge, une information d’épaisseur optique représentée par la couverture nuageuse. En référence à la , on a représenté un exemple de donnée acquise par un capteur de nébulosité. Le champ de vision du capteur correspond à la zone circulaire au centre de la figure. A chaque pixel contenu dans cette zone circulaire, correspond une direction, définie par des valeurs d’azimut et d’élévation, par rapport au capteur. De plus, chaque pixel est associé à une épaisseur optique qui dépend de l’épaisseur de la couche nuageuse présente sur la ligne de visée correspondant à cet azimut et cette élévation. Sur la , les épaisseurs optiques sont représentées en nuances de gris, mais l’information relative à l’épaisseur optique en chaque pixel peut prendre la forme d’une valeur scalaire comprise entre deux bornes minimale et maximale prédéterminées. Dans le cas où le système 1 comprend un ou plusieurs capteurs de nébulosité 50, le terminal de contrôle 30 est avantageusement en communication avec ce ou ces capteurs, par exemple par une connexion filaire, cette communication pouvant avoir lieu en temps réel ou à une fréquence déterminée. En outre, le terminal de contrôle 30 peut également recevoir des données de couverture nuageuse acquises par un satellite d’observation météorologique 2, ou des données de couverture nuageuse correspondant à une prévision météorologique générée par un modèle numérique (non représenté). Chaque station au sol peut également être en communication avec un capteur de proximité, tel qu’un capteur de nébulosité au sol, et un terminal de traitement local dédié. Ce terminal de traitement local permet alors d’estimer les probabilités de faisabilité de la liaison optique avec le satellite à différents horizons de temps. Ces probabilités sont transmises à un autre terminal de traitement global qui aura pour charge de décider et sélectionner la station au sol qui devra établir une liaison optique avec le satellite pendant une période de temps future. Cela permet avantageusement de minimiser les volumes de données transmis entre les stations au sol, leur capteur de proximité et le centre de gestion des stations. Dans le cas d’un modèle de risque pour évaluer une probabilité locale à court terme d’établissement de la liaison optique, comme ceux basés sur les réseau de neurones, le modèle de risque de chaque station au sol sera par exemple défini en fonction des données d’historiques de couverture nuageuse locale propres à chaque station au sol. On pourrait aussi envisager de partager les données d’historique de couverture nuageuse locale, générées par plusieurs capteurs de proximité pour définir un modèle de risque partagé entre plusieurs stations, pour établir les prévisions à court terme. En référence aux figures 2, 3 et 4, on va maintenant décrire un procédé de communication par liaison optique entre le satellite et les stations au sol. Ce procédé permet de sélectionner au moins une station au sol 20 parmi une pluralité de stations pour établir une communication par liaison optique ininterrompue avec le satellite 10, sur la base de données de couverture nuageuse des stations au sol. Dans des exemples de réalisation, le procédé permet de sélectionner une séquence S de plusieurs stations au sol 20 successivement utilisées pour l’établissement de la communication optique, de sorte à assurer une communication continue malgré une couverture nuageuse de certaines stations, cette couverture étant variable dans le temps. Le procédé comprend l’obtention 100 de données relatives à la couverture nuageuse d’un ensemble de stations au sol 20. Ces données peuvent par exemple comprendre, pour une station au sol considérée, une indication de la couverture nuageuse pour une pluralité de directions à partir de la station, chaque direction étant déterminée par des valeurs d’azimut et d’élévation. De plus, ces données peuvent comprendre, pour une station au sol considérée, et le cas échéant pour au moins une direction déterminée à partir de cette station, une séquence temporelle de données décrivant l’évolution temporelle de la couverture nuageuse, la séquence temporelle couvrant une durée déterminée et comprenant des données acquises à une fréquence déterminée. Cette fréquence peut par exemple être comprise entre une acquisition toutes les trente secondes et une acquisition toutes les dix minutes. De plus, la durée déterminée est par exemple comprise entre quelques minutes et plusieurs heures. Dans des exemples de réalisation, des données relatives à la couverture nuageuse des stations au sol sont acquises en continu à une fréquence déterminée, et transmises au terminal de contrôle 30 en temps réel ou à intervalle de temps fixe. Les données relatives à la couverture nuageuse comprennent avantageusement des données, du type image, acquises par des capteurs de nébulosité situés au voisinage des stations au sol considérées. Par « au voisinage », on entend que la distance entre le capteur de nébulosité et la station au sol est suffisamment faible pour que les informations relatives à la couverture nuageuse du capteur soient également applicables à la station. Typiquement, cette distance est inférieure à 100 mètres, avantageusement inférieure à 50 mètres, très avantageusement comprise entre 0 et 10 mètres. Comme indiqué précédemment, un tel capteur 50 fournit, pour une pluralité de directions par rapport au capteur, une valeur d’épaisseur optique. L’épaisseur optique fournie par le capteur permet de déterminer directement l’atténuation d’un signal optique causée par les nuages par la formule suivante : Où OD est l’épaisseur optique et où « attenuation » est l’atténuation évaluée en décibels. On considère typiquement qu’une communication par liaison optique est possible lorsque l’atténuation due aux nuages est inférieure à un seuil de 3 dB. Lorsque l’atténuation est supérieure, on peut considérer la liaison optique comme impossible. Toutefois, en augmentant la puissance d’émission, on pourrait compenser une perte plus grande. Par ailleurs, il y a d’autres perturbations à prendre en compte, comme par exemple les turbulences, non liées à l’épaisseur optique qui pourrait justifier le fait de prendre davantage de marge. En complément, ou alternativement, les données relatives à la couverture nuageuse comprennent des données, du type image, d’observation de couverture nuageuse acquises par au moins un satellite 2 observant une région de la Terre intégrant une ou plusieurs stations au sol de l’ensemble considéré. Les données obtenues par un satellite d’observation météorologique 2 peuvent également comprendre des données du type images infrarouges qui sont traitées pour calculer une épaisseur optique représentée par la couverture nuageuse, qui est comme précédemment convertie en valeur d’atténuation causée par les nuages. Un satellite d'observation a par exemple une résolution de l’ordre de 0.05 degrés voire inférieure, ce qui équivaut à des pixels de taille 5km voire de l'ordre du kilomètre, à l’équateur, les pixels représentant des zones plus étendues lorsqu'on s’éloigne de l’équateur. Dans des exemples de réalisation, les données obtenues comprennent à la fois des données acquises par des capteurs 50 de nébulosité voisins des stations au sol et des données obtenues par un satellite d’observation météorologique 2. En effet, les capteurs de nébulosité sont en mesure de fournir des données locales, particulièrement pertinentes à court terme, tandis que des données de satellite sont des données régionales permettant de déterminer des tendances à plus long-terme, qui sont ainsi complémentaires. En outre les capteurs de nébulosité fournissant des données locales, permettent une évaluation plus précise de la couverture sur la ligne de visée entre la station au sol et le satellite. En complément, ou alternativement, les données relatives à la couverture nuageuse peuvent également comprendre des données de prévision de couverture nuageuse d’une zone géographique donnée, ces données étant obtenues par application de modèles numériques de simulation, tels que les modèles utilisés pour les prévisions météorologiques. Par ailleurs, le procédé comprend également l’obtention 200, par le terminal de contrôle 30, de données relatives à une orbite du satellite considéré, ces données permettant de déterminer, compte-tenu des localisations géographiques des stations au sol, un sous-ensemble de stations au sol avec lesquelles le satellite est susceptible d’établir une communication par liaison optique en l’absence de couverture nuageuse. Le satellite mémorise par exemple les données accumulées sur plusieurs orbites différentes correspondant à plusieurs périodes orbitales. Par exemple dans le cas d’un satellite LEO, le satellite se décale régulièrement, par rapport à la terre, d’une orbite à l’autre. On peut par exemple vider les données stockées à bord du satellite une ou deux fois par jour, ce qui revient par exemple à établir une liaison optique environ toutes les 15 orbites pour un satellite à environ 600km d’altitude, en considérant par exemple une dizaine de stations au sol candidates pendant cette journée. Ce sous-ensemble de stations peut être variable dans le temps dans le cas d’un satellite défilant, c’est-à-dire non géostationnaire, mais dans ce cas, l’orbite du satellite étant connue, la composition du sous-ensemble de stations avec lesquelles une communication par liaison optique est possible, en l’absence de nuage, est connue pour la position actuelle du satellite comme pour les positions futures. L’enjeu est alors, parmi l’ensemble des stations au sol candidates, de sélectionner suffisamment de stations au sol avec le choix du créneau de communication associé lors du passage du satellite, afin de garantir une durée de communication totale sur l’horizon par exemple pour garantir le vidage depuis le satellite d’un volume de donnée préétabli. Le planificateur s’assure alors que les stations sont disponibles, et réserve suffisamment de créneaux de communication où la prédiction de faisabilité est bonne pour garantir le vidage. Une surréservation peut être envisagée pour pallier les cas où la prédiction serait erronée et où la communication prévue ne se ferait pas faute d’un lien bloqué par les nuages. Ainsi des stations au sol peuvent être sélectionnées sur la base de durées prédites cumulées pour lesquelles il est possible d’établir des liaisons optiques ininterrompues entre des stations au sol et un satellite, à partir d’instants déterminés d’établissement de liaisons optiques fonction de créneaux de visibilité et en fonction de prévisions de couvertures nuageuses, en prenant en compte deux horizons de temps. Par exemple, les données relatives à l’orbite du satellite peuvent être enregistrées dans la mémoire 32 et le calculateur 31 peut y accéder à l’étape 200 pour déterminer, parmi l’ensemble de stations pour lesquelles des données relatives à la couverture nuageuse sont obtenues à l’étape 100, un sous-ensemble de stations avec lesquelles le satellite peut établir une communication par liaison optique pendant un intervalle de temps déterminé. L’intervalle de temps est par exemple un intervalle de temps d’une durée déterminée à partir du moment auquel est mis en œuvre le procédé. Le calculateur 31 détermine ensuite au cours d’une étape 300, à partir des données de couverture nuageuse obtenues, des prévisions relatives à la possibilité d’établir une liaison optique ininterrompue entre le satellite et au moins une station au sol, pour au moins deux horizons de temps déterminés, respectivement à court-terme et long-terme. Les termes « court-terme » et « long-terme » sont interprétés de façon relative, c’est-à-dire que le premier horizon de temps est plus court que le second. A titre d’exemple non limitatif, le premier horizon de temps peut être compris entre 0 et 10 minutes, par exemple entre 30 secondes et 5 minutes. Le deuxième horizon de temps, à long terme, peut être compris par exemple entre 30 minutes et 3 heures, voire entre 30 et 90 minutes. Dans des exemples de réalisation, cette étape est mise en œuvre pour chaque station au sol du sous-ensemble de stations déterminé à partir des informations obtenues sur la position du satellite. Dans la suite, le terme « horizon de temps » est utilisé pour décrire un intervalle de temps compris entre un instant actuel, correspondant typiquement à l’instant de mise en œuvre du procédé, et un instant futur, l’intervalle de temps étant de durée déterminée. On comprend donc que l’horizon de temps « long-terme » inclut, de préférence strictement, l’horizon de temps « court-terme ». Le calculateur 31 est donc configuré pour prédire le risque qu’il existe au moins une interruption de la liaison optique entre le satellite 10 et la station au sol 20 considérée pendant cet intervalle de temps. Dans des exemples de réalisation, le calculateur 31 est configuré pour déterminer une probabilité qu’il existe au moins une interruption pendant l’horizon de temps considéré, et appliquer un seuillage à la probabilité pour obtenir une indication binaire sur la possibilité, ou non, d’établir la liaison optique en continu pendant l’horizon de temps. Cette étape 300 de détermination des prévisions relatives à la possibilité d’établir une liaison optique ininterrompue est mise en œuvre par application, sur les données de couverture nuageuse obtenues à l’étape 100, d’au moins un modèle M de risque, par exemple par un réseau de neurones préalablement entrainé sur des données d’historiques de couvertures nuageuses pour un ensemble de stations au sol. Dans des exemples de réalisation, les données d’historiques de couvertures nuageuses utilisées pour l’entrainement du modèle comprennent des données d’historiques de couvertures nuageuses de l’ensemble des stations au sol pour lesquelles des données de couverture nuageuse actuelle sont obtenues à l’étape 100. En variante, ces données d’historique peuvent concerner des stations au sol différentes de celles pour lesquelles les données relatives à la couverture nuageuse sont obtenues à l’étape 100. Les données d’historiques de couvertures nuageuses utilisées pour l’entrainement du modèle M de risque comprennent avantageusement, pour une station au sol considérée, des séquences temporelles de données comprenant une indication de couverture nuageuse ou d’atténuation déterminée, à une fréquence d’échantillonnage déterminée, par exemple comprise entre quelques secondes et quelques minutes, par exemple de l’ordre de 30 secondes à 2 minutes. Dans des exemples de réalisation, ces séquences temporelles de données comprennent des indications de couverture nuageuse ou d’atténuation, déterminée pour une pluralité de directions à partir de la station au sol considérée. Par ailleurs, comme indiqué précédemment à propos de l’étape 100 d’obtention de données relatives à la couverture nuageuse, ces données d’historique de couvertures nuageuses peuvent être des données obtenues par des capteurs de nébulosité au sol voisins des stations considérées, et/ou des données d’observation météorologiques obtenues par satellite. Dans des exemples de réalisation, l’étape 300 de détermination des prévisions relatives à la possibilité d’établir une liaison optique ininterrompue est mise en œuvre par application de deux modèles distincts, par exemple par des réseaux de neurones, dans lequel un premier modèle est configuré pour le calcul de prévisions à court-terme, c’est-à-dire au premier horizon de temps, et le deuxième modèle est configuré pour le calcul de prévisions à long-terme, c’est-à-dire pour au moins un deuxième horizon de temps ultérieur au premier. Ainsi, les horizons de chacun des deux modèles permettent avantageusement deux usages spécifiques. Le premier modèle à court terme permet de déterminer le risque de coupure par exemple sur un lien en cours d’utilisation et fournit une plus grande précision, alors que le deuxième modèle, à long terme, a pour but par exemple de faire une première sélection d’une nouvelle station, lorsqu’un risque de coupure important apparaît sur la station courante. Le réseau qui utilise des informations à court terme et des informations à long terme, permet avantageusement de minimiser le nombre de changements de stations au sol courantes. Ces modèles de risque peuvent être établis à partir des mêmes données d’historiques de couverture nuageuse. Alternativement, le modèle de risque configuré pour le calcul de prévisions à court-terme peut être établi à partir de données d’historiques locales, c’est-à-dire obtenues par des capteurs de nébulosité au sol, tandis que le modèle de risque configuré pour le calcul de prévisions à long-terme peut être établi à partir de données d’historiques régionales, c’est-à-dire des données d’observation de satellite. En référence à la , l’étape 300 de prévision de la possibilité d’établir une liaison optique ininterrompue entre le satellite et des stations au sol permet d’obtenir, pour chaque station au sol, représentée par les termes OGS1…OGSn, une indication relative à la possibilité d’établir une liaison avec le satellite de façon ininterrompue pour les au moins deux horizons de temps considéré. Sur la , on a représenté schématiquement un nombre n de stations au sol, et un ensemble d’horizons de temps H1,…Hj. Sur les figures 3 et 4, on a représenté en blanc les moments où une communication par liaison optique est possible entre le satellite 10 et la station 20 considérée, et en noir les moments où cette communication est impossible. L’étape 300 de détermination des prévisions relatives à la possibilité d’établir une liaison optique ininterrompue permet donc de déterminer par exemple : qu’une liaison optique ininterrompue est possible pour un horizon de temps court-terme H1 pour la première station, que ce n’est pas le cas pour la deuxième station, qu’une liaison optique ininterrompue est également possible pour un horizon de temps long-terme H2 pour la dernière station. Le procédé comprend ensuite une étape 400 de sélection d’au moins une station au sol 20 pour l’établissement d’une communication par liaison optique avec le satellite 10, à partir des informations de prévisions relatives à la possibilité d’établir une liaison optique ininterrompue obtenues à l’étape 300 pour chacun des horizons de temps considérés. La station est sélectionnée sur la base de la durée la plus longue prévue, en considérant les prévisions à long terme, pour laquelle il est possible d’établir une liaison optique ininterrompue entre la station et le satellite, à partir d’un instant déterminé d’établissement de la liaison optique. Cet instant d’établissement de la liaison optique peut être un instant actuel, notamment dans le cas où la station considérée est une station dite « courante », pour laquelle une communication par liaison optique est déjà établie avec le satellite et pour laquelle le risque de coupure est proche, selon les prévisions à court terme. Alternativement, cet instant d’établissement peut être un instant ultérieur, qui se situe entre l’instant actuel et un horizon de temps pour lequel un risque de coupure est établi pour la station courante ou, en d’autres termes, auquel la liaison optique est prédite comme impossible, selon les prévisions à court terme. Ainsi, dans des exemples de réalisation, l’étape 400 de sélection d’une station au sol 20 pour l’établissement d’une communication par liaison optique avec le satellite 10 comprend la détermination d’une séquence S de communication comprenant une suite d’au moins deux stations au sol 20 utilisées successivement pour l’établissement d’une communication par liaison optique avec le satellite 10, la séquence S étant déterminée de façon à satisfaire un critère portant sur l’absence de coupure de la communication au cours de la séquence. Ainsi, la séquence peut par exemple comprendre la détermination d’une première station au sol avec laquelle la communication est établie, pour laquelle aucune coupure n’est prévue pendant un premier horizon de temps, puis une deuxième station pour laquelle aucune coupure n’est prévue entre la fin du premier horizon de temps et la fin d’un deuxième horizon de temps, ultérieur au premier. De plus, et afin de limiter les transitions d’une station à l’autre, la séquence S peut également être déterminée de façon à respecter un deuxième critère portant sur le nombre de changement de stations au cours de la séquence, ce deuxième critère étant par exemple de minimiser le nombre de changement de stations, ou de respecter un nombre maximum de changements de stations, au cours de la séquence. Dans ce cas, le premier critère relatif à l’absence de coupure au cours de la séquence est prioritaire sur le deuxième critère. Par exemple, toujours en référence à la , une séquence S, représentée par des flèches, est déterminée dans laquelle : la première station au sol OGS1 est utilisée pour la communication jusqu’à un instant correspondant à un risque de coupure de la communication, et à partir de cet instant, la communication est basculée vers la deuxième station OGS2 pour laquelle aucun risque de coupure n’est identifié à partir de cet instant et jusqu’à un horizon de temps ultérieur. On remarque que la troisième station OGSn qui est représentée n’est pas sélectionnée car, bien qu’aucun risque de coupure ne soit prévu pour cette station au moment où le risque de coupure est établi pour la première station, l’horizon de temps pour lequel la liaison est ininterrompue est plus court que pour la deuxième station, ce qui implique plus de changements de stations. En référence à la , dans un exemple de réalisation, l’étape 300 de détermination de prévisions sur la possibilité d’établir une liaison optique ininterrompue entre le satellite 10 et chacune d’une pluralité de stations au sol 20 comprend une première détermination 310 de la possibilité d’établir une liaison optique ininterrompue entre le satellite et une seule station au sol, qui est la station courante avec laquelle une liaison optique est déjà établie, pour au moins un horizon de temps déterminé. Sur la , on a ainsi représenté par les flèches « p » cette première détermination, qui est mise en œuvre de façon itérative à une fréquence déterminée (représentée par l’intervalle delta sur la ), tant que la liaison optique avec la station courante est prédite comme possible sur l’horizon de temps considéré. De préférence, la fréquence à laquelle cette première détermination est mise en œuvre correspond à un intervalle de temps entre deux itérations qui est inférieur à l’horizon de temps H p considéré pour la prédiction. Ainsi sur la , on observe que la quatrième itération de cette détermination pour la station courante est la dernière, car cette itération détermine un risque de coupure de la liaison optique dans l’horizon de temps considéré. Dans ce cas, quand une liaison optique continue est prédite comme impossible pour l’horizon de temps considéré, l’étape de détermination comprend une deuxième détermination 320 d’une possibilité d’établir une communication par liaison optique entre le satellite et les autres stations au sol pour au moins deux horizons de temps considéré, comprenant un horizon de temps court-terme et un horizon de temps long-terme. On revient ensuite à l’étape 400 de sélection, parmi lesdites autres stations, d’au moins une station 20 pour l’établissement de la liaison optique, la bascule entre la station courante et la station sélectionnée ayant lieu avant la fin de l’horizon de temps H p considéré pour la première détermination et de préférence dès que la deuxième station a été sélectionnée. Également dans ce cas, l’étape 400 de sélection d’une station au sol 20 pour l’établissement d’une communication par liaison optique avec le satellite 10 peut comprendre la sélection d’une séquence S de stations satisfaisant le ou les critères mentionnés ci-avant. De retour à l’exemple de la , on a indiqué par des flèches « dm » les mises en œuvre de la deuxième détermination pour les autres stations que la station courante initialement. On observe que les stations n°2 et n représentées ont toutes les deux un horizon de temps à court-terme correspondant à une liaison optique possible au moment de cette deuxième détermination, mais que la station n°2 présente également un horizon long-terme correspondant à une liaison optique possible, contrairement à la station n. Par conséquent, c’est la deuxième station qui est sélectionnée. On ne se limite pas à une station en liaison optique avec un satellite. Le procédé peut être appliqué à N stations en service parmi M stations potentielles (avec N De plus, une fois que cette station est sélectionnée, le procédé revient à la première détermination 310 pour la station courante, qui est itérée à nouveau jusqu’à ce qu’un risque de coupure soit établi pour la station considérée, comme cela est représenté par la répétition des flèches « p » pour la deuxième station OGS2. Bien que les étapes 100 d’obtention de données relatives à la couverture nuageuse et 200 d’obtention de données relatives à l’orbite du satellite ne soient pas représentées sur les figures 3 et 4, on comprend que l’ensemble des stations au sol considérées pour l’établissement potentiel d’une liaison optique est mis à jour au cours du temps, en fonction des données relatives à la position du satellite obtenues à l’étape 200, et que les données de couvertures nuageuses pour ces stations sont également mises à jour par des répétitions de l’étape 100 d’obtention de données relatives à la couverture nuageuse. Dans un exemple de réalisation, la première détermination 310 relative à la station courante est mise en œuvre par application d’un premier modèle, ce premier modèle ayant été préalablement entrainé sur des données d’historiques de couverture nuageuse d’un ensemble de station, pour prédire la possibilité d’établir une liaison optique ininterrompue entre le satellite et une même station donnée pendant un horizon de temps déterminé. Cet horizon de temps pouvant être court-terme, et étant plus cours que l’horizon de temps le plus long considéré lors de la deuxième détermination 320 relatives aux autres stations, les données d’historiques de couverture nuageuse utilisées pour l’apprentissage du modèle, puis les données de couverture nuageuse utilisées pour l’inférence du modèle, peuvent être des données locales obtenues par un capteur de nébulosité voisin de la station considérée. Cet exemple n’exclut cependant pas que les données d’apprentissage et les données utilisées pour l’inférence comprennent également des données d’observation satellite et/ou des données obtenues par simulation numérique. Dans ce cas, comme la première détermination 310 est mise en œuvre pour une station au sol courante pour laquelle une communication par liaison optique est déjà établie avec le satellite, la base de données d’apprentissage peut comprendre exclusivement des séquences de données couvrant une durée supérieure ou égale à l’horizon de temps considéré, et comprenant des données acquises à une fréquence déterminée durant cette durée, dont les premières données de la séquence correspondent à une absence de couverture nuageuse, c’est-à-dire à une liaison optique possible entre le satellite et la station. Par ailleurs, la deuxième détermination 320 relative aux autres stations et la sélection d’une deuxième station pour poursuivre la communication peuvent être mises en œuvre par application d’un deuxième modèle, entrainé pour réaliser ces tâches. Puisque la deuxième détermination est mise en œuvre en considérant au moins deux horizons de temps, dont au moins l’un se termine de préférence ultérieurement à l’horizon de temps considéré pour la première détermination, les données d’historiques de couverture nuageuse utilisées pour l’apprentissage du modèle, puis les données de couverture nuageuse utilisées pour l’inférence du modèle, peuvent être des données régionales obtenues par un satellite d’observation météorologique. Ces données peuvent être complétées par des données acquises par des capteurs de nébulosité et/ou des données obtenues par simulation numérique. Dans le cas des capteurs locaux, tels que les capteurs de nébulosité, le volume des données générées est proportionnel au nombre de capteurs pour les stations concernées, contrairement au premier modèle qui utilise les données satellite. On exploite ainsi les données générées au niveau de la station au sol pour anticiper son basculement et les données générées au niveau de toutes les stations libres candidates pour servir de relai. Alternativement, les données d’historiques de couverture nuageuse et les données de couverture nuageuse utilisées respectivement pour l’apprentissage et l’inférence du modèle peuvent être des données de capteurs de nébulosité, ce qui permet d’utiliser les mêmes données pour les deux modèles. Alternativement, la première détermination 310 relative à la station courante, la deuxième détermination 320 relatives aux autres stations et la sélection 400 d’une station au sol 20 pour l’établissement d’une communication par liaison optique peuvent être mises en œuvre par un même modèle qui peut avoir fait l’objet d’un premier entrainement par exemple pour la deuxième détermination 320 et la sélection 400 d’une station au sol 20 pour l’établissement d’une communication par liaison optique, puis d’un apprentissage par renforcement pour mettre en œuvre la première détermination 310 avec précision. L’apprentissage par renforcement permet à un agent d’apprendre à opérer le système. Ainsi un tel agent, pourrait décider, à chaque pas de temps, de rester ou de changer de station, un tel agent étant également à même de choisir la prochaine station. Un tel apprentissage inclurait notamment le risque de coupure qui correspondrait à un cas de pénalisation très fort pour l’agent en question. Dans ce cas un même agent traduirait un modèle de risque sur les deux horizons de temps, à court terme et à long terme. Le procédé peut également comprendre la mise à jour des données d’historiques de couverture nuageuse au fur et à mesure que de nouvelles données sont obtenues, ainsi que la mise à jour du ou des modèles présentés ci-avant, sur la base de ces données et de l’écart entre les prévisions du modèle et la réalité des données de couverture nuageuse observée après avoir mis en œuvre le procédé. En référence aux figures 6 et 7, on va maintenant décrire un exemple de mise en œuvre de modèles utilisés pour la première détermination 310 relative à la station courante d’une part, et pour la deuxième détermination et la sélection d’une station d’autre part. Les données utilisées en entrée de chaque modèle sont des séquences de données, par exemple du type image, acquises à une fréquence 1 d’acquisition par exemple toutes les 30 secondes, par des capteurs de nébulosité voisins de stations au sol respectives. Les séquences de données couvrent par exemple une période de 5 minutes. Les pixels des images comportent par exemple des couleurs représentatives de la nébulosité observée. Chaque pixel représente par exemple un couple (azimut, élévation). A chaque pixel correspond par exemple une épaisseur optique (un scalaire entre 0 et 5) qui dépend de l’épaisseur du nuage présent sur chaque pixel. Les données sont par exemple prétraitées pour supprimer des zones inutiles des images, c’est-à-dire au moins les bords de l’image hors de la zone circulaire visible dans l’exemple de la . Un prétraitement peut aussi être réalisé par exemple pour extraire de la séquence de données une succession d’échantillons correspondant, pour au moins une direction donnée par rapport au capteur, à la série temporelle des données acquises par le capteur pour cette direction. Ce prétraitement peut par exemple être mis en œuvre par utilisation du format de données TFRecord de la suite logicielle TensorFlow, permettant de stocker de manière optimisée les échantillons d’apprentissage. Les données satellite d’observation de la terre sont par exemple du type image dans le spectre visible ou invisible, comme les rayons infra-rouges. Les images satellites sont traitées par des procédés connus de génération d’images satellite. Les images satellites fournies sont par exemple centrées autour du point d’intérêt qu’est la station. Les images satellite utilisées se présentent par exemple sous la forme d’une matrice. Les données, par exemple du type image matricielle, utilisées pour l’entraînement des réseaux de neurones, pour la mise à jour des réseaux de neurones et pour les prédictions réalisées avec les réseaux de neurones, sont par exemple du même type. Les données actualisées sont par exemple utilisées pour calculer une inférence d’un réseau et réaliser une prédiction. Les données d’archive permettent de calculer des inférences dans le passé, afin de vérifier des prédictions. La descente de gradient et les algorithmes associés sont par exemple utilisés pour la phase d'apprentissage. Pour ce faire, TensorFlow implémente par exemple la différenciation automatique, qui utilise le calcul pour déterminer les gradients, en générale, pour calculer le gradient de l'erreur. L’erreur calculée en sortie est par exemple rétro-propagée dans le réseau de neurones, lors de la phase d’apprentissage ou lors d’une mise à jour du réseau de neurones. Les deux modèles utilisés peuvent être par exemple mis en œuvre par des réseaux de neurones, par exemple de type convolutifs, recevant en entrée les séquences d’images prétraitées mentionnées ci-avant, référencé « I », se présentant par exemple sous la forme d’une matrice et produisant en sortie une matrice de scalaires compris entre 0 et 1, où chaque élément de la matrice correspond à une même position de pixel des images d’entrée, donc par exemple à une même direction par rapport au capteur de nébulosité, et déterminant si le lien est faisable ou non sur des intervalles de temps de [t, t+k] minutes avec minutes. Chaque intervalle de temps correspond donc à un horizon de temps considéré par le modèle. Un exemple de structure de réseau de neurones pouvant être utilisé est représenté sur la . Ce réseau comprend une succession de couches convolutives (dont le nombre de canaux d’entrée « ch » et la dimension du noyau de convolution « kernel » sont indiqués sous chaque couche) suivies d’une couche d’activation impliquant par exemple la fonction ReLu. Le nombre de canaux en sortie du réseau indiqué par « frame_out » correspond au nombre d’intervalles que nous souhaitons prédire en sortie du réseau. Le réseau est donc configuré de sorte qu’il reste « frame_out » couches (ou « channels » en anglais) en sortie du réseau et ainsi avoir le format souhaité en sortie. Par exemple si on a les intervalles [0, 5min], [0,10min] et [0, 15min] en sortie, le frame_out sera de 3. Le réseau peut également comprendre une couche de rognage (« crop ») permettant de réduire la sortie aux pixels correspondant à une direction supérieure à un seuil d’élévation donné. La sortie du réseau correspond à la sortie d’une fonction sigmoïde. Selon cet exemple, les entrées du réseau comprennent donc 11 matrices de nombres réels, représentant l’épaisseur optique sur chaque pixel d’une image acquise par un capteur de nébulosité, de taille 460x640. Les 11 matrices correspondent aux images acquises toutes les 30 secondes pendant cinq minutes entre les cinq dernières minutes et l’instant courant où la prévision est faire. Les sorties du réseau comprennent un nombre de matrices égale au nombre d’horizons de temps considérés, par exemple 8 selon l’exemple ci-dessus, chaque matrice comprenant des éléments scalaires entre 0 et 1 estimant la probabilité que la ligne de visée selon la direction correspondant à l’élément de la matrice soit continûment dépourvue de nuage sur l’horizon de temps considéré. Comme indiqué ci-avant, la base de données d’apprentissage pour l’entrainement du premier modèle pourrait par exemple être réduite à des données pour lesquelles le lien optique est initialement possible. De plus, pour le premier modèle qui prédit un risque qu’une liaison optique soit impossible sur au moins un horizon de temps considéré, un seuil est par exemple appliqué sur la sortie scalaire pour lever une alerte, qui déclenche l’appel au deuxième modèle. En revanche, pour le deuxième modèle, une exécution de ce modèle est faite pour chaque autre station que la station courante, et la station choisie peut être celle minimisant cette sortie scalaire, qui ne fait donc pas l’objet d’un seuillage. En référence à la , on a représenté les résultats obtenus par l’exemple présenté ci-avant, les courbes représentant la probabilité P qu’une communication par liaison optique soit interrompue au cours d’une simulation d’une durée de 3ans, avec la mise en œuvre du procédé décrit ci-avant (courbe A) pour choisir une séquence de stations au sol, comparée à un modèle persistant simple (courbe B) et à une séquence obtenue avec une connaissance parfaite des données de nébulosité future (courbe C). Les courbes sont représentées en fonction du nombre Nb de station au sol. Pour cette simulation on a tiré des dates aléatoirement et vérifié si le dm provoquait une coupure ou pas, les dates tirées étant générées sur 3 ans de données réelles. On constate que les performances obtenues sont meilleures que celles du modèle persistant et s’améliorent avec le nombre de stations, pour être assez proches des performances théoriques optimales correspondant à la courbe C. Procédé de communication par liaison optique entre au moins un satellite (10) et une pluralité de stations émettrices et/ou réceptrices au sol (20), le procédé étant mis en œuvre par un calculateur et comprenant : l’obtention (100) de données actualisées de couvertures nuageuses propres à chaque station au sol, l’obtention (200) de données relatives à une orbite du satellite, la détermination (300), pour chaque station au sol (20), compte tenu des données relatives à l’orbite du satellite et par l’application d’au moins un modèle de risque, de prévisions sur la possibilité d’établir une liaison optique ininterrompue entre le satellite et la station au sol, en fonction des données actualisées de couverture nuageuse propre à chaque station au sol, pour au moins deux horizons de temps déterminés, à court terme et à long terme, ledit modèle de risque étant établi sur la base de données d’historiques de couvertures nuageuses des stations au sol, et la sélection (400) d’au moins une des stations au sol pour l’établissement d’une communication par liaison optique montante ou descendante avec le satellite, cette station au sol étant sélectionnée sur la base de la durée prédite la plus longue pour laquelle il est possible d’établir une liaison optique ininterrompue entre la station et le satellite, à partir d’un instant déterminé d’établissement de la liaison optique, en fonction des prévisions sur chacun des deux horizons de temps. Procédé selon la revendication 1, dans lequel ledit modèle de risque fournit des prévisions sur la possibilité d’établir la liaison optique ininterrompue, associées à un indice de confiance, la sélection de ladite station étant réalisée également en fonction d’une durée d’établissement du futur lien et de l’indice de confiance, un poids plus ou moins faible étant respectivement attribué à la durée d’établissement du futur lien selon que l’indice de confiance est respectivement plus ou moins important. Procédé selon la revendication 1 ou 2, dans lequel l’orbite du satellite est défilante et impose des créneaux temporels de visibilité propres à chaque station au sol, où les communications sont envisageables, ledit modèle de risque étant appliqué pour un ensemble de stations au sol dont les créneaux temporels de visibilité se recouvrent. Procédé selon la revendication 3, dans lequel la sélection comprend plusieurs stations au sol (20) pour l’établissement de communications par liaison optique descendantes avec le satellite (10), ces stations au sol étant sélectionnées sur la base de durées prédites cumulées pour lesquelles il est possible d’établir des liaisons optiques ininterrompues entre les stations au sol et le satellite, à partir d’instants déterminés d’établissement de liaisons optiques propres à chaque station au sol et en fonction des prévisions, propres à chaque station au sol, sur deux horizons de temps. Procédé selon l’une des revendications précédentes, dans lequel l’étape de détermination (300) de prévisions utilise deux modèles de risque distincts, du type réseau de neurones, dont l’un est entrainé sur la base de données d’historiques de couvertures nuageuses locales des stations au sol pour le calcul de prévisions à court terme et l’autre est entraîné sur la base de données d’historiques de couvertures nuageuses globales des stations au sol pour le calcul de prévisions à long terme. Procédé selon l’une des revendications précédentes, comprenant en outre une étape de mise à jour des données d’historiques de couvertures nuageuses par une intégration des données actualisées de couverture nuageuses et une étape de mise à jour dudit modèle de risque. Procédé selon l’une des revendications précédentes, dans lequel la sélection (400) d’au moins une station au sol comprend la détermination d’une séquence (S) de communication comprenant une suite d’au moins deux stations au sol (20) utilisées successivement pour l’établissement de la communication optique montante ou descendante avec le satellite (10), la séquence étant déterminée de façon à satisfaire au moins un premier critère portant sur l’absence de coupure de la communication optique au cours de la séquence. Procédé selon la revendication précédente, dans lequel la séquence (S) est déterminée de façon à satisfaire en outre un deuxième critère portant sur le nombre de changements de stations au cours de la séquence. Procédé selon l’une des revendications précédentes, dans lequel les données actualisées de couverture nuageuse obtenues comprennent des données acquises par des capteurs de nébulosité (50) au sol, des données d’observation météorologique obtenues par au moins un satellite, et des données obtenues par la mise en œuvre de modèles numériques de prédiction. Procédé selon l’une des revendications précédentes, dans lequel la pluralité de stations au sol (20) comprend une station dite courante en communication par liaison optique avec le satellite lors de la mise en œuvre du procédé, et l’étape de détermination comprend : une première détermination (310) d’une possibilité d’établir une liaison optique entre le satellite et la station courante pendant au moins un horizon de temps déterminé, et ; si une liaison optique entre le satellite et la station courante est prédite comme impossible à l’horizon de temps: la mise en œuvre d’une deuxième détermination (320) d’une possibilité d’établir une liaison optique entre le satellite et les autres stations au sol pendant au moins deux horizons de temps déterminés, cette deuxième détermination étant suivie de la sélection de l’une parmi les autres stations au sol pour l’établissement de la communication par liaison optique. Procédé selon la revendication précédente, dans lequel l’instant d’établissement de la liaison optique avec la station sélectionnée est déterminé comme un instant antérieur à l’horizon de temps pour lequel une liaison optique est prédite comme impossible entre le satellite et la station courante. Procédé selon la revendication 10 ou 11, dans lequel au moins l’un des horizons de temps considérés lors de la deuxième détermination se termine ultérieurement à l’horizon de temps considéré pour la première détermination, et le procédé comprenant la mise en œuvre d’un premier modèle entrainé pour la première détermination relative à la station courante, et d’un deuxième modèle entrainé pour la deuxième détermination relative aux autres stations. Procédé selon l’une des revendications 10 à 12, comprenant l’obtention de données d’observation de nébulosité acquises par un capteur pour la station courante, et la mise en œuvre de la première détermination (310) sur la base desdites données d’observation de nébulosité, et l’obtention de données d’observation météorologiques obtenues par au moins un satellite et couvrant une région incluant les autres stations au sol, et la mise en œuvre (320) de la deuxième détermination sur la base desdites données d’observation météorologiques. Procédé selon l’une des revendications précédentes, dans lequel les données relatives à la couverture nuageuse au voisinage d’une station au sol comprennent une série temporelle de données relatives à la couverture nuageuse au voisinage d’une station au sol acquises à une fréquence déterminée. Produit programme d’ordinateur, comprenant des instructions de code pour la mise en œuvre du procédé selon l’une des revendications qui précèdent, lorsqu’il est mis en œuvre par un calculateur (31). Système de communication par liaison optique entre un satellite (10) et une pluralité de stations émettrices et/ou réceptrices au sol (20), le système comprenant au moins un terminal de contrôle (30) comprenant un calculateur (31) et une mémoire (32), le terminal de contrôle étant adapté pour communiquer avec les stations au sol, caractérisé en ce qu’il est configuré pour mettre en œuvre le procédé selon l’une des revendications 1 à 14.