DOMAINE TECHNIQUE

La présente invention appartient au domaine de l'imagerie terrestre par des engins spatiaux, tels que des satellites d'observation, et concerne plus particulièrement un procédé de planification d'acquisition d'images de zones terrestres par un engin spatial ou une constellation de tels engins spatiaux. L'invention trouve une application particulièrement avantageuse, bien que nullement limitative, dans le cas de satellites d'observation fonctionnant selon un mode d'acquisition d'images par bandes, mode désigné par la suite dans le texte par le terme « peigne » connu de l'homme de l'art (mode encore dit « pushbroom » en référence à la littérature anglo-saxonne).

ÉTAT DE LA TECHNIQUE

Les missions d'observation terrestre effectuées par un engin spatial consistent à acquérir des images de zones terrestres, c'est-à-dire situées à la surface de la Terre, en réponse à des requêtes de clients. En particulier, un tel engin spatial suit une orbite défilante autour de la Terre de sorte à être en mesure d'acquérir lesdites zones terrestres lorsqu'il les survole pendant une durée prédéterminée. En outre, le développement de l'agilité d'un tel engin spatial permet d'accroître le nombre de zones terrestres imagées, ces dernières pouvant désormais être situées de part et d'autre de la dite orbite ou bien acquises sous différents angles au cours de plusieurs orbites. Ainsi, à tout instant le long de ladite orbite défilante correspond une ou plusieurs opportunités d'acquisition d'images de zones terrestres différentes.

La réception desdites requêtes par l'engin spatial s'effectue de manière régulière, généralement sur une base journalière. Actuellement, les requêtes pour de telles acquisitions de zones terrestres sont en nombre croissant, car non plus limitées uniquement aux seuls secteurs industriels historiquement liés au secteur de l'imagerie spatiale. Par exemple, et à titre non limitatif, le secteur de l'agriculture fait désormais un usage important d'observations terrestres dans le but d'optimiser l'exploitation de surfaces agraires. Dès lors, le nombre de requêtes à traiter par un tel engin spatial ne cesse d'augmenter, ce qui est aussi le cas de leur complexité dans la mesure où ces requêtes comportent des contraintes bien précises associées à la zone terrestre devant être acquise, comme par exemple des conditions de localisation, d'exposition lumineuse, ou encore de nécessité d'acquisitions multiples dans le but d'obtenir des images stéréo, tri-stéréo ou multibandes. A cela s'ajoute aussi le fait que les requêtes sont différentiables selon leur niveau de priorité.

La gestion des contraintes associées auxdites requêtes doit aussi être envisagée en parallèle de la gestion de contraintes d'exploitation associées audit engin spatial, qu'elles soient de type cumulatives (taille mémoire, puissance électrique consommée, temps maximum de fonctionnement d'instruments embarqués), ou bien encore de type locales (durée minimale entre deux acquisitions successives).

On comprend donc que l'objectif consistant à satisfaire un ensemble de telles requêtes est un problème très fortement contraint nécessitant la planification dans le temps, selon l'orbite dudit engin spatial, des prises de vues à effectuer. Plus largement, une telle problématique s'inscrit dans la classe des problèmes d'optimisation dit du « voyageur de commerce », bien connus pour être très difficiles à résoudre en un temps raisonnable.

Le problème est particulièrement critique dans le cas où la demande de l'utilisateur doit être satisfaite dans un délai court. L'augmentation du nombre de satellites donne la possibilité d'avoir des accès satellites sur une zone dans des délais beaucoup plus courts. Cependant cette capacité ne peut être pleinement utilisée que si la mise à jour du plan d'acquisition peut également se faire dans des délais courts.

S'il est néanmoins connu qu'un tel problème de planification admet des solutions d'approximation de bonne qualité grâce à des techniques heuristiques d'optimisation (par exemple de type gloutonnes), leurs temps de calcul n'apportent pas satisfaction eu égard au volume de requêtes à traiter en conditions opérationnelles.

Plus récemment, des techniques de résolution basées sur des hypothèses de simplifications importantes, notamment concernant une gestion découplée des contraintes locales et cumulatives dudit engin, ont permis d'apporter en temps raisonnable des solutions à la gestion desdites contraintes locales , ce qui est une première étape avant la résolution globale dudit problème de planification. Si certaines de ces techniques permettent effectivement de prendre en compte les contraintes cumulatives, aucune d'entre elles ne permet toutefois de prendre en compte l'acquisition d'images stéréo, tri-stéréo ou multibandes, notamment dans le cas d'un satellite d'observation fonctionnant en mode peigne. Par conséquent, ces techniques sont limitatives et non optimales quant au nombre de requêtes qu'elles peuvent satisfaire.

EXPOSÉ DE L'INVENTION

La présente invention a pour objectif de remédier à tout ou partie des inconvénients de l'art antérieur, notamment ceux exposés ci-avant, en proposant une solution qui permette d'avoir un procédé de planification optimal quant à la satisfaction de contraintes locales et cumulatives d'un engin spatial et qui permette de prendre en compte des zones terrestres nécessitant d'être acquises une pluralité de fois ou bien encore une seule fois.

A cet effet, l'invention concerne un procédé de planification d'acquisition d'images de zones terrestres Zi ,...,Z _N , par un engin spatial en mission autour de la Terre suivant une orbite défilante prédéterminée, comme par exemple un satellite d'observation, chaque zone terrestre Z, étant associée à une requête R, de sorte qu'à ladite zone Z, correspond un intervalle temporel d'accessibilité visuelle Ti le long de ladite orbite défilante, 1 < i < N. En outre, ledit procédé comporte les étapes successives suivantes :

- une étape de détermination pour chaque intervalle T, d'opportunités d'acquisition discrètes de la zone Z, de sorte qu'à chaque dite opportunité d'acquisition discrète sont associées une date de début incluse dans ledit intervalle T,, une durée d'exécution, une contrainte locale cinématique ainsi qu'un poids,

- une étape de regroupement des opportunités d'acquisition discrètes des zones Z, dans un ensemble D,

- une étape de classement des opportunités d'acquisition discrètes de l'ensemble D, les opportunités d'acquisition étant classées par date de début croissante et désignées par D ¹ ,

- une étape d'évaluation de compatibilité cinématique entre lesdites opportunités d'acquisition discrètes de l'ensemble D, prises deux à deux selon leurs classements respectifs,

- une étape de détermination d'une séquence optimale d'opportunités d'acquisition discrètes ayant un poids maximal, et compatible cinématiquement.

Dans des modes particuliers de mise en œuvre, le procédé de planification d'acquisition d'images de zones terrestres comporte l'une ou plusieurs des caractéristiques suivantes, prises isolément ou selon toutes les combinaisons techniquement possibles.

Dans un mode particulier de mise en œuvre, le procédé comporte une étape de prise en compte de requêtes d'acquisitions multiples, consécutive à l'étape d'évaluation de compatibilité cinématique et antérieure à l'étape de détermination d'une séquence optimale, dans laquelle les poids respectifs des opportunités d'acquisition discrètes associées à des zones Z, devant être acquises une pluralité de fois sont mis à jour et dans laquelle une variable mustHit force l'enchaînement desdites opportunités d'acquisition discrètes afin de satisfaire lesdites requêtes d'acquisitions multiples.

Dans un mode particulier de mise en œuvre, l'étape de détermination d'une séquence optimale comporte la mise à jour dynamique d'une variable pathMesh stockant les zones Z, associées aux opportunités d'acquisition discrètes de ladite séquence, de sorte à empêcher la duplication d'opportunités d'acquisition discrètes de zones devant être acquises une seule fois.

Dans un mode particulier de mise en œuvre, l'étape de détermination d'une séquence optimale comporte pour chaque opportunité d'acquisition discrète D' de l'ensemble D une sélection partielle d'opportunités d'acquisition discrètes précédant D' dans D, et compatibles cinématiquement avec D', parmi toutes les opportunités d'acquisition discrètes précédant D' dans D, et compatibles cinématiquement avec D'.

Dans un mode particulier de mise en œuvre, le procédé comporte une étape de vérification de respect d'au moins une contrainte cumulative imposée audit engin spatial, ladite étape de vérification étant consécutive à l'étape de détermination d'une séquence optimale.

Dans un mode particulier de mise en œuvre, ladite au moins une contrainte cumulative comporte une contrainte de taille mémoire.

Dans un mode particulier de mise en œuvre, ladite au moins une contrainte cumulative comporte une contrainte de puissance électrique consommée.

Dans un mode particulier de mise en œuvre, ladite au moins une contrainte cumulative comporte une contrainte de temps maximum de fonctionnement d'instruments embarqués dans ledit engin.

Dans un mode particulier de mise en œuvre, l'étape de vérification de respect de contraintes cumulatives comporte une mise à jour dynamique de ladite séquence optimale par suppression d'opportunités d'acquisition discrètes, tant que lesdites contraintes cumulatives ne sont pas satisfaites.

Dans un mode particulier de mise en œuvre, les dates de début des opportunités d'acquisition discrètes de chaque intervalle T, sont échantillonnées temporellement selon un pas de temps constant.

Dans un mode particulier de mise en œuvre, les pas de temps respectifs des intervalles Ti , ...,T _N sont égaux entre eux.

Dans un mode particulier de mise en œuvre, les pas de temps respectifs d'au moins deux intervalles T, et T _j sont différents.

Dans un mode particulier de mise en œuvre, le pas de temps de chaque intervalle T, est ajusté de sorte que ledit intervalle T, est associé à 20 opportunités d'acquisition discrètes.

Dans un mode particulier de mise en œuvre, la contrainte locale cinématique d'une opportunité d'acquisition discrète comporte une consigne d'attitude dudit engin spatial au début et à la fin de ladite opportunité d'acquisition discrète.

Dans un mode particulier de mise en œuvre, le poids de chaque opportunité d'acquisition discrète associée à une zone Z, dépend de la priorité de la requête R, associée à ladite zone Z,.

Dans un mode particulier de mise en œuvre, chaque zone Z, est associée à un critère qualité w, dépendant des caractéristiques géométriques de ladite zone Z, ainsi que des conditions météorologiques prévues au cours de l'intervalle temporel d'accessibilité visuelle Ti.

Dans un mode particulier de mise en œuvre, le poids de chaque opportunité d'acquisition discrète associée à une zone Z, dépend de la priorité de la requête R, de sorte à être un vecteur dont la seule composante non nulle est égale au critère qualité w,, la position de ladite composante non nulle dans ledit vecteur dépendant de la priorité de ladite requête R,.

Dans un mode particulier de mise en œuvre, ladite composante non nulle dépend du critère qualité w, et d'un critère d'intérêt d'acquisition associé à ladite opportunité d'acquisition discrète.

Dans un mode particulier de mise en œuvre, l'étape de détermination d'une séquence optimale comporte une comparaison des poids respectifs des opportunités d'acquisition discrètes compatibles cinématiquement entre elles.

Dans un mode particulier de mise en œuvre, le procédé comporte une étape d'agrégation de zones terrestres, antérieure à l'étape de détermination d'opportunités d'acquisition discrètes.

Dans un mode particulier de mise en œuvre, ledit engin spatial est un satellite d'observation fonctionnant en mode peigne.

PRÉSENTATION DES FIGURES

Les caractéristiques et avantages de l'invention seront mieux appréciés grâce à la description qui suit, description qui expose les caractéristiques de l'invention au travers de modes de mise en œuvre préférés, qui n'en sont nullement limitatifs.

La description s'appuie sur les figures annexées qui représentent :

- Figure 1 : une représentation d'un organigramme d'un exemple de mise en œuvre d'un procédé de planification d'acquisition d'images de zones terrestres Z-i, ...,Z _N par un engin spatial.

- Figure 2 : une représentation d'un mode préférentiel de mise en œuvre du procédé de la figure 1 au cours duquel le procédé comporte une étape de vérification de respect d'au moins une contrainte cumulative imposée audit engin spatial.

- Figure 3 : une représentation d'un mode particulier de mise en œuvre du procédé de la figure 1 au cours duquel le procédé comporte une étape d'agrégation de zones terrestres. - Figure 4 : une représentation d'un mode préférentiel de mise en œuvre du procédé de la figure 1 au cours duquel sont prises en compte les zones terrestres devant être acquises une unique fois ou bien une pluralité de fois.

DESCRIPTION DÉTAILLÉE D'UN MODE DE RÉALISATION DE

L'INVENTION

La présente invention appartient au domaine de la planification d'acquisition d'images de zones terrestres par un engin spatial en mission autour de la Terre suivant une orbite défilante prédéterminée.

La figure 1 représente un organigramme d'un exemple de mise en œuvre d'un procédé de planification d'acquisition d'images de zones terrestres ZI , .. . ,ZN par un engin spatial.

Par « acquisition d'images », on entend ici la mesure d'un rayonnement électromagnétique reçu desdites zones terrestres Zi, ...,Z _N . A cet effet, l'engin spatial comporte une optique destinée à réaliser une pluralité d'acquisitions d'images, ladite optique comportant en outre des capteurs adaptés à mesurer ledit rayonnement électromagnétique. Lesdites images sont finalement converties au format numérique et stockées au sein d'une mémoire de capacité limitée embarquée par ledit engin spatial.

La suite de la description vise plus spécifiquement, mais de manière non limitative, un procédé de planification d'acquisition d'images de zones terrestres par un satellite d'observation agile en orbite défilante (LEO ou MEO, acronymes des expressions anglais respectives « Low Earth Orbit » et « Médium Earth Orbit »). Par exemple, ledit satellite d'observation agile suit une orbite héliosynchrone quasi-polaire à altitude constante. Il comporte un instrument adapté à acquérir des images en mode peigne par bandes successives selon différentes directions d'acquisition.

Rien n'exclut cependant, suivant d'autres exemples non détaillés, de considérer d'autres types d'engins spatiaux (navette spatiale, sonde, ...) adaptés à acquérir des images de zones terrestres, ainsi que d'autres types d'instruments ou d'autres modes d'acquisition d'image, par exemple avec des détecteurs matriciels en mode dit « step and stare ».

Les zones terrestres Z-i , ...,Z _N sont des surfaces au sol, d'aires respectives finies, éventuellement contigûes et intersectent la fauchée dudit satellite d'observation. L'agilité du satellite autorise les zones terrestres acquises à être positionnées d'une part, le long de la projection de ladite orbite défilante sur la surface terrestre, et d'autre part de par et d'autre de cette projection sur la surface terrestre.

Dans la suite de la description, chaque zone terrestre Z, est associée à une requête R, émise par un utilisateur. Chaque requête R, pouvant être elle- même associée à différents paramètres :

- une priorité p,, p, étant par exemple un entier naturel, qui caractérise la nécessité de satisfaire cette requête R, avant éventuellement d'autres requêtes. Chaque priorité p, dépend de plusieurs facteurs comme par exemple l'importance et le coût de la requête R, de laquelle elle dépend,

- des contraintes d'acquisition dépendant de plusieurs facteurs, comme par exemple l'attitude du satellite (qui détermine l'orientation de la direction de visée), ou bien encore les conditions météorologiques nécessaires à une acquisition de bonne qualité et qui doivent être satisfaites aux moments où le satellite est potentiellement en mesure d'acquérir la zone Z, associée à ladite requête R,,

- un type d'acquisition (mono, stéréo, tri-stéréo). En effet, certaines zones Zi nécessitent d'être acquises une pluralité de fois, éventuellement selon un angle d'acquisition différent comme par exemple dans le cas de la stéréo ou de la tri-stéréo.

Il est à noter que d'une manière générale, c'est-à-dire au moment où une requête R, est signifiée par un utilisateur, la zone Z, associée à ladite requête R, est de forme quelconque. C'est pourquoi, et ceci dans un but d'optimisation de fonctionnement, une campagne d'acquisition d'images par le satellite est en règle générale précédée d'une étape de mise en forme consistant à découper ladite zone Z, en zones élémentaires adaptées à être acquises par ledit satellite. Par exemple, dans le cas d'une prise de vue en mode peigne, lesdites zones Z, sont découpées en bande de largeur inférieure ou égale à celle de la fauchée du satellite. Pour la suite de la description, on considère que les zones Z, associées aux requêtes R, sont de formes respectives adaptées à être acquises par l'optique dudit satellite. Autrement dit, ladite étape de mise en forme est soit inutile eu égard aux formes initiales respectives desdites zones Zi, soit considérée comme ayant été préalablement effectuée.

Par ailleurs, à chaque zone terrestre Z, correspond un intervalle temporel d'accessibilité visuelle T, le long de ladite orbite défilante. Ledit intervalle T, comporte une date de début ainsi qu'une date de fin, et correspond à une durée temporelle pendant laquelle le satellite peut acquérir une image de la zone Z, tout en respectant les contraintes d'acquisition ainsi que le type d'acquisition de la requête R, à laquelle est associée ladite zone Z,.

Pour la suite de la description, on adopte aussi la convention que les composantes d'un vecteur sont lues de la gauche vers la droite en étant indexées à partir de 1 . Ainsi un vecteur V à N composantes est noté V = (V ¹ , V ² , V ^N ). De la même manière, un élément positionné à la i-ème ligne et j- ème colonne d'une matrice M est noté M'' ^j .

Le procédé de planification d'acquisitions d'images comporte plusieurs étapes successives. Dans son principe général, ledit procédé consiste à déterminer tout d'abord des séquences d'acquisition comportant des points le long de l'orbite défilante du satellite, de sorte qu'en réalisant des acquisitions d'images en les points respectifs desdites séquences, le plus grand nombre de requêtes d'utilisateurs soient satisfaites. Ensuite, ledit procédé a pour objectif de déterminer parmi lesdites séquences d'acquisition celle satisfaisant un critère d'optimalité défini ci-après ainsi que des contraintes inhérentes au fonctionnement dudit satellite.

A cet effet, le procédé comporte dans un premier temps une étape 100 de détermination, pour chaque intervalle T,, d'opportunités d'acquisition discrètes a , ... de la zone Z, de sorte qu'à chaque dite opportunité d'acquisition discrète sont associées une date de début incluse dans ledit intervalle T,, une durée d'exécution, une contrainte locale cinématique ainsi qu'un poids.

Une opportunité d'acquisition discrète a{ est un point le long de ladite orbite défilante qui désigne la possibilité pour le satellite d'acquérir une image d'une zone Z,, à partir d'une date de début ainsi que pendant une durée d'exécution comprises dans l'intervalle T, de ladite zone Z,, de sorte que soient respectées les contraintes d'acquisition ainsi que le type d'acquisition de la requête R, associée à cette zone Z,.

Dans la suite de la description, une séquence d'acquisition désigne une suite d'opportunités d'acquisition discrètes.

Un intervalle T, comporte une ou plusieurs opportunités d'acquisition discrètes qui discrétisent ce dernier selon leurs dates de début respectives. En effet, l'agilité du satellite autorise le fait qu'une zone Z, puisse être acquise potentiellement en des temps distincts de son intervalle TY

Dans un mode particulier de mise en œuvre de l'étape 100, les dates de début des opportunités d'acquisition discrètes de chaque intervalle T, sont échantillonnées temporellement selon un pas de temps constant. Il est à noter que plus ledit pas de temps constant est choisi petit, plus il y a d'opportunités d'acquisition au sein de chaque intervalle Tj, augmentant ainsi le nombre de séquences d'acquisition possibles. En contrepartie, un pas de temps petit augmente la complexité de calcul du procédé.

Dans un mode particulier de mise en œuvre de l'étape 100, les pas de temps respectifs d'au moins deux intervalles T, et T _j sont différents. Une telle configuration est avantageuse lorsque par exemple les priorités respectives des requêtes R, et R _j associées auxdits intervalles T, et T _j sont différentes. En effet, un moyen de satisfaire une requête de priorité élevée consiste à associer ladite requête à un grand nombre d'opportunités d'acquisition discrètes de sorte à augmenter le nombre de points le long de l'orbite en lesquels le satellite effectue des acquisitions de la zone associée à ladite requête, Autrement dit, plus la priorité d'une requête est élevée, plus il est avantageux de discrétiser finement l'intervalle d'accessibilité associé à cette dernière au moyen d'un pas de temps petit ceci afin d'avoir un nombre important d'opportunités d'acquisition discrètes au sein dudit intervalle d'accessibilité.

Dans un mode particulier de mise en œuvre de l'étape 100, les pas de temps respectifs des intervalles Ti ,...,T _N sont égaux entre eux. De cette manière, un seul pas de temps est choisi, ce qui est avantageux en terme de temps de calcul lorsque les priorités respectives des requêtes Ri ,...,R _N associées aux intervalles Ti, ...,T _N sont toutes identiques.

Dans un mode particulier de mise en œuvre de l'étape 100, le pas de temps de chaque intervalle T, est ajusté de sorte que ledit intervalle T, est associé à un nombre prédéfini d'acquisition discrètes. Par exemple, un nombre d'acquisitions discrètes égal à 20 est avantageux dans la mesure où la pratique montre qu'il s'agit d'un bon compromis entre l'optimisation du procédé du point de vue du temps de calcul et la possibilité de construire suffisamment de séquences d'acquisition en vue de satisfaire au mieux un grand nombre de requêtes.

Par ailleurs, la durée d'exécution correspond au temps nécessaire à l'optique du satellite pour réaliser l'acquisition sous l'hypothèse que l'attitude dudit satellite le permet dès ladite date de début de l'acquisition. Ainsi, la durée d'exécution d'une opportunité d'acquisition discrète ne comprend pas le temps nécessaire à la modification de l'attitude du satellite de sorte que l'optique de ce dernier soit configurée pour réaliser ladite acquisition ou bien une éventuelle autre acquisition ultérieure. A titre d'exemple nullement limitatif, et correspondant aux configurations les plus usitées, le temps d'exécution d'une opportunité d'acquisition est compris entre 5 secondes et 15 secondes.

La contrainte locale cinématique d'une opportunité d'acquisition discrète correspond à la configuration dans laquelle doit se trouver le satellite, et donc l'orientation de l'axe de visée de l'instrument, au début et à la fin d'une opportunité d'acquisition discrète afin de satisfaire la requête à laquelle est associée ladite opportunité d'acquisition discrète.

Par exemple, la contrainte locale cinématique d'une opportunité d'acquisition discrète comporte une consigne d'attitude dudit engin spatial au début et à la fin de ladite opportunité d'acquisition discrète. Par consigne d'attitude, on entend ici l'orientation du satellite lors d'une acquisition afin de satisfaire la requête associée à ladite acquisition.

Le poids c[a{] d'une opportunité d'acquisition discrète a{ correspond à une mesure générale de l'intérêt d'acquérir une zone terrestre en ladite opportunité d'acquisition discrète au regard de l'ensemble des requêtes reçues par le satellite.

Par exemple, le poids de chaque opportunité d'acquisition discrète associée à une zone Z, dépend de la priorité p, de la requête R, associée à ladite zone Z _h de sorte que plus ladite priorité est élevée, plus les poids respectifs desdites opportunités d'acquisition discrètes sont importants.

Le poids d'une opportunité d'acquisition discrète relève aussi avantageusement de caractéristiques de la zone terrestre à laquelle ladite opportunité est associée. Ainsi, dans un mode particulier de mise en œuvre de l'étape 100, chaque zone Z, est associée à un critère qualité w[Z,] dépendant des caractéristiques géométriques de ladite zone Z, ainsi que des conditions météorologiques prévues au cours de l'intervalle temporel d'accessibilité visuelle Ti. Ledit critère qualité w[Z,] est une quantité scalaire caractérisant, par exemple, le fait que la zone Zi est de taille restreinte et donc plus facile à acquérir qu'une zone de plus grande taille, mais aussi que les conditions météorologiques rencontrées pendant l'intervalle temporel d'accessibilité visuelle T, associé à Z, sont sensiblement identiques à celles exigées par les contraintes d'acquisition de la requête reliée audit intervalle T,.

Dans un mode particulier de mise en œuvre de l'étape 100, le poids de chaque opportunité d'acquisition discrète associée à une zone Z, dépend de la priorité p, de la requête R, de sorte à être un vecteur dont la seule composante non nulle est égale audit critère qualité w[Zj], la position de ladite composante non nulle dans ledit vecteur étant la priorité p, de ladite requête R, _. Dans ce mode de mise en œuvre, la taille dudit vecteur est égale au nombre de priorités distinctes utilisées pour classifier les requêtes reçues par le satellite lors d'une mission. En outre, on adopte la convention qu'une priorité p, est de plus grande importance qu'une autre priorité P _j si et seulement si p, < P _j . Par exemple, une opportunité d'acquisition discrète a{ d'une zone Z, ayant un critère qualité w[Zi], et étant associée à une requête R, de priorité p, égale à 2, parmi un ensemble de requêtes dont les priorités vont de 1 à 4, a pour poids

c[a{] = (c ¹ ^ ], c ² [a{], c ³ [a{], c ⁴ [a ]) = (0, w[Z _f ], 0, 0).

Une telle représentation des poids respectifs des opportunités d'acquisition discrètes est avantageuse puisqu'elle permet de comparer les poids entre eux au moyen d'une relation d'ordre naturelle notée > opérant de la façon suivante : on note c [a] > c [b] au sens que le poids d'une opportunité d'acquisition discrète a est supérieur au poids d'une opportunité d'acquisition discrète b si et seulement si

{cHa] > c ¹ ^]} V {(cHa] = c ¹ ^]) Λ (c ² [a] > c ² [b])} V

les opérateurs Λ et v représentant respectivement les opérateurs booléens « et » et « ou », et l'opérateur > étant l'opérateur classique de comparaison entre scalaires. On note que la relation d'ordre > opère en comparant des vecteurs composante par composante, et est ainsi un ordre lexicographique.

Une telle représentation est aussi avantageuse car elle permet de définir de manière naturelle le poids d'une séquence d'acquisition comportant une pluralité d'opportunités d'acquisition discrètes comme la somme des poids desdites opportunités d'acquisition discrètes, la sommation vectorielle se faisant composante par composante.

Dans un mode particulier de mise en œuvre de l'étape 1 00, ladite composante non nulle du poids d'une opportunité d'acquisition discrète a{ dépend du critère qualité w[Z,] ainsi que d'un critère d'intérêt d'acquisition w[a ] associé à ladite opportunité d'acquisition discrète a{ . Ledit critère d'intérêt d'acquisition w[a ] est un scalaire qui mesure l'intérêt d'acquérir plus particulièrement l'opportunité a{ plutôt qu'une autre opportunité associée elle aussi à la même zone Z,. Par exemple, et à titre nullement limitatif, le critère d'intérêt d'acquisition dépend de l'angle d'acquisition de chaque opportunité d'une même zone Z\. En outre, dans ce mode de mise en œuvre, le poids d'une opportunité d'acquisition discrète est mis à jour de sorte qu'en reprenant l'exemple ci-avant :

c[a{] = (0, w[Zi] + w[a{], 0, 0).

Une telle manière de procéder est avantageuse dans la mesure où le poids d'une opportunité d'acquisition discrète tient à la fois compte des caractéristiques de la zone et de la requête à laquelle est associée ladite opportunité d'acquisition discrète, mais aussi des différences entre les opportunités d'acquisition discrètes d'une même zone. Ainsi, le poids c[a{] mesure de façon exhaustive l'intérêt et la variabilité existant entre chaque opportunité d'acquisition discrète.

Le procédé comporte ensuite une étape 200 de regroupement des opportunités d'acquisition discrètes des différentes zones Z, dans un ensemble D.

Chaque élément de l'ensemble D est une opportunité d'acquisition discrète à laquelle est associée une date de début, une durée d'exécution, une contrainte locale cinématique ainsi qu'un poids. L'ensemble D, tel que décrit dans cette étape 200, regroupe ainsi toutes les opportunités susceptibles de faire partie d'une ou plusieurs séquences d'acquisition en vue de satisfaire un maximum de requêtes.

Le procédé comporte ensuite une étape 300 de classement des opportunités d'acquisition discrètes de l'ensemble D par date de début croissante.

Ladite étape 300 consiste donc à ordonner les éléments de l'ensemble

D. Il est possible de classer les éléments de D selon d'autres caractéristiques des opportunités d'acquisition discrètes. Par exemple, ces dernières peuvent être classées par poids croissants (au sens de la relation d'ordre >), par durées d'exécution croissantes ou bien encore selon les requêtes auxquelles elles sont associées. Toutefois, dans le présent mode de mise en œuvre, le choix s'est avantageusement porté sur un classement par date de début croissante pour des raisons exposées ci-après.

Il est à noter que tout sous ensemble de l'ensemble D ainsi trié constitue une séquence d'acquisition potentielle, c'est-à-dire une suite d'opportunités d'acquisition discrètes respectant la chronologie du mouvement du satellite le long de son orbite défilante. Il convient néanmoins de s'assurer que le satellite puisse effectivement réaliser les unes après les autres les opportunités d'acquisitions discrètes de ladite séquence d'acquisition potentielle.

A cet effet, le procédé comporte ensuite une étape 400 d'évaluation de compatibilité cinématique entre lesdites opportunités d'acquisition discrètes de l'ensemble D, prises deux à deux selon leurs classements respectifs. Par compatibilité cinématique, on entend ici que deux opportunités d'acquisition discrètes D' et D ^j de l'ensemble D, telles que ladite opportunité D' ait une date de début antérieure à la date de début de ladite opportunité D ^j (donc telles que i < j), sont chaînables par le satellite le long de son orbite. Cela est théoriquement possible lorsque le temps de manœuvre nécessaire au satellite afin de se rendre de l'opportunité D' à l'opportunité D ^j est inférieur à la durée séparant la date de fin de l'opportunité D' de la date de début de l'opportunité D ^j .

Par extension, une séquence d'acquisition A qui est un sous ensemble de l'ensemble D est dite compatible cinématiquement lorsque ses éléments sont deux à deux compatibles cinématiquement, c'est-à-dire lorsque A ^j est compatible cinématiquement avec A ^l+1 .

Classiquement, les éléments de l'ensemble D peuvent être vus comme les sommets d'un graphe G, ces sommets étant reliés entre eux au moyen d'arrêtés lorsqu'ils sont compatibles cinématiquement entre eux. La structure de l'ensemble D, obtenue à l'issue de l'étape 300, conditionne avantageusement la structure du graphe G dans la mesure où :

- le classement des éléments de D par date de début croissante fournit une orientation naturelle aux arêtes du graphe de G dans le sens des temps croissants : le graphe G est donc orienté,

- une opportunité D' n'est pas compatible cinématiquement avec une opportunité d'acquisition discrète D ^j dont la date de début est antérieure à celle de ladite opportunité D' (donc telles que i > j). Cela revient à dire que les éléments de toute séquence d'acquisition compatible cinématiquement sont deux à deux distincts. Cela garantit l'impossibilité du satellite à passer deux fois au même endroit le long de son orbite : le graphe G est donc aussi acyclique.

Dès lors, on comprend que chercher à planifier l'acquisition d'images de zones terrestres revient à chercher des séquences d'acquisition compatibles cinématiquement au sein dudit graphe G. Plus particulièrement, dans l'objectif de satisfaire le plus grand nombre de requêtes possible, cela revient à chercher une séquence d'acquisition compatible cinématiquement et ayant le poids le plus élevé au sein du graphe G. Il s'agit là d'un problème d'optimisation, encore connu sous l'expression de « problème de recherche du plus long chemin » (« longest path problem » en anglais).

Le fait que le graphe G soit un graphe orienté acyclique est avantageux dans la mesure où il est connu que dans ce cas ledit problème d'optimisation est de complexité polynomiale en nombre de sommets du graphe G, alors qu'il est de complexité exponentielle pour un graphe quelconque.

Dans le présent mode de mise en œuvre, le graphe G est représenté au moyen d'une matrice carrée M de booléens, dite matrice d'adjacence, dont le terme situé à la ligne i et à la colonne j est :

- M'' ^j = 1 si D' et D ^j sont compatibles cinématiquement,

- M ^l,j = 0 sinon.

Les conditions de compatibilité entre deux opportunités d'acquisition discrètes impliquent que M est une matrice triangulaire supérieure stricte.

Il est à noter que la mise en œuvre des étapes 100 à 400 est de préférence effectuée au sol, et en amont de la mission d'observation du satellite, c'est-à-dire éventuellement longtemps avant la résolution du problème d'optimisation. Une telle mise en œuvre est avantageuse dans la mesure où les tâches effectuées pendant les étapes 100 à 400, et plus particulièrement celle de détermination des opportunités d'acquisition discrètes de l'étape 100 ainsi que celle de détermination du graphe G au cours de l'étape 400, sont facilement exécutables par des moyens informatiques selon un schéma de programmation parallèle, permettant ainsi un gain de temps. En outre, une première mise en œuvre des étapes 100 à 400 permet avantageusement de s'affranchir d'une nouvelle détermination complète du graphe G lorsqu'on souhaite par exemple inclure à la volée dans le plan de mission du satellite une requête supplémentaire, le graphe G étant en effet rapidement et facilement mis à jour. De la même manière, la mise à jour de priorités et/ou de critères qualité de zones terrestres et/ou de poids d'opportunités d'acquisition discrètes peut être effectuée au dernier moment avant la résolution du problème d'optimisation.

Dans une étape 500, consécutive à l'étape 400, on détermine une séquence optimale d'opportunités d'acquisition discrètes ayant un poids maximal, et compatible cinématiquement.

Ladite séquence optimale est cherchée au sein du graphe G en parcourant la partie située au-dessus de la diagonale de la matrice M selon une méthode de calcul de plus long chemin, et fournit un plan d'acquisition théorique suivi par le satellite au cours de sa mission d'observation.

Pour ce faire, on procède en deux temps. Dans un premier temps, on détermine pour chaque opportunité d'acquisition discrète D' de l'ensemble D un poids objectif obj(D'), qui est le poids maximal des séquences compatibles cinématiquement et menant à ladite opportunité d'acquisition discrète D'. Pour cela, on introduit un ensemble P(D') de toutes les opportunités d'acquisition discrètes précédant D' dans l'ensemble D, et compatibles cinématiquement avec D'. Alors le poids objectif obj(D') est égal au poids de D' auquel s'ajoute le maximum des poids objectifs des opportunités d'acquisition discrètes de P(D'), c'est-à-dire

obj(O ^[ ) = c[O ^[ ] + max _{sep(Di )} (oô;(s)).

En faisant désormais l'hypothèse que l'ensemble D est de cardinal Q(D), et en introduisant une variable previous{D) qui stocke l'opportunité d'acquisition discrète de l'ensemble P(D') maximisant le poids de la séquence optimale menant à D', la détermination du poids objectif de tous les éléments de l'ensemble D se fait dynamiquement au moyen d'un premier algorithme s'écrivant :

Pour i allant de 1 à Q(D) ; obj(Û) = c[Û] Fin de Pour ;

Pour i allant de 1 à Q(D) ; previous(D') = - 1 Fin de Pour

Pour i allant de 1 à Q(D) :

Pour j allant de (i+ 1) à Q(D) :

Si IV(' ^J = 1 :

Si obj(Û) + c[Û] > obj(O) :

previous(O) = D ¹ ;

obj(Û) = obj(Û) + c[0] ;

Fin de Si ; Fin de Si ;

Fin de Pour ;

Fin de Pour ; II est à noter que la recherche de la séquence optimale par la méthode de calcul du plus long chemin, telle que présentée dans le premier algorithme, est avantageusement adaptée au fait que les opportunités d'acquisition discrètes de D sont classées par dates de début croissantes au cours de l'étape 300. En effet, un tel classement correspond à un ordre topologique sur le graphe G qui facilite la mise en place des étapes séquentielles dudit premier algorithme.

Dans un deuxième temps, ladite séquence optimale est déterminée en sélectionnant parmi les opportunités d'acquisition discrètes de D, celle de poids objectif maximal, puis en remontant itérativement les prédécesseurs de ladite opportunité de poids objectif maximal au moyen des informations stockées dans la variable previous.

Comme noté ci-dessus, le classement des opportunités d'acquisition discrètes de D fournit avantageusement un ordre topologique sur le graphe G, qui dès lors est orienté acyclique, et par conséquent toute séquence cinématiquement compatible obtenue en sortie de l'étape 500 est optimale au sens de la maximalité du poids.

Dans un mode particulier de mise en œuvre du procédé, l'étape 500 comporte pour chaque opportunité d'acquisition discrète D' de l'ensemble D une sélection partielle d'opportunités d'acquisition discrètes précédant D' dans l'ensemble D, et compatibles cinématiquement avec D', parmi toutes les opportunités d'acquisition discrètes précédant D' dans l'ensemble D, et compatibles cinématiquement avec D'.

Une telle manière de procéder revient à restreindre l'ensemble P(D') pour chaque opportunité d'acquisition discrète D' de l'ensemble D. Ainsi, ledit ensemble P(D') peut être avantageusement remplacé par un ensemble P _C (D') défini comme un sous-ensemble de ce dernier, et comportant les opportunités d'acquisition discrètes précédant D' dans l'ensemble D, compatibles cinématiquement avec D' ainsi que proches de D'. Par proches de D', on entend ici que les éléments de P _C (D') occupent respectivement une position dans D de sorte que la différence entre cette position et celle de D' reste bornée par une constante prédéfinie.

Dès lors, le poids objectif obj(D') s'obtient par la formule

obj(O ^[ ) = c[O ^[ ] + max _sePc(Di) (oô; (s)).

A cet effet, on introduit une variable maxNextVertice(D') qui contient pour chaque opportunité d'acquisition D' le nombre maximum de sommets du graphe G visités au cours de la recherche de la séquence optimale. Le scalaire maxNextVertice(D') constitue ainsi un majorant du cardinal de l'ensemble P _C (D')- Dès lors, le premier algorithme devient :

Pour i allant de 1 à Q(D) ; obj(Û) = c[Û] Fin de Pour ;

Pour i allant de 1 à Q(D) ; previous(D') = - 1 Fin de Pour

Pour i allant de 1 à Q(D) :

Pour j allant de (i+ 1) à maxNextVertice(Û) :

Si M ^J = 1 :

Si obj(Û) + c[Û] > obj(O) :

previous(O) = D ¹ ;

obj(Û) = obj(Û) + c[0] ;

Fin de Si ;

Fin de Si ;

Fin de Pour ;

Fin de Pour ; Limiter l'ensemble P(D') à l'ensemble P _C (D') est avantageux car l'exploration des sommets de G par la méthode de calcul de plus long chemin est réduite en temps de calcul. Par exemple, et à titre nullement limitatif, lorsque l'ensemble P _C (D') contient uniquement des opportunités d'acquisition discrètes séparées de soixante secondes au plus de D ¹ le long de l'orbite défilante, le temps de calcul est réduit de 60% à 70% et la séquence optimale reste inchangée par rapport au cas où P _C (D') n'est pas utilisé.

La figure 2 représente un mode préférentiel de mise en œuvre du procédé de la figure 1 au cours duquel le procédé comporte, à la suite de l'étape 500, une étape 600 de vérification de respect d'au moins une contrainte cumulative imposée audit satellite.

Une contrainte cumulative imposée au satellite est une contrainte qui prend en compte l'ensemble des choix faits dans la planification de la mission d'observation avant que le satellite n'atteigne sur son orbite une opportunité d'acquisition discrète. Autrement dit, et en opposition aux contraintes locales cinématiques associées à chacune des opportunités d'acquisition discrètes, une contrainte cumulative a un caractère global relativement à l'historique de mission du satellite.

Par exemple, et à titre nullement limitatif, ladite au moins une contrainte cumulative comporte une contrainte de taille mémoire. En effet, en mission d'observation, le satellite stocke sous format numérique les acquisitions qu'il réalise dans une mémoire de type connu en soi et de taille limitée.

Alternativement ou en complément de l'étape 600, ladite au moins une contrainte cumulative comporte une contrainte de puissance électrique consommée. En effet, en mission d'observation, le satellite consomme pour se déplacer de l'énergie électrique stockée dans des batteries de capacités respectives limitées. Il en est de même pour le fonctionnement de l'optique du satellite. Il est à noter que lesdites batteries sont rechargeables au moyen généralement de panneaux solaires.

Alternativement ou en complément de l'étape 600, ladite au moins une contrainte cumulative comporte une contrainte de temps maximum de fonctionnement d'instruments embarqués dans ledit engin.

Rien n'exclut, suivant d'autres exemples non détaillés ici, d'avoir d'autres contraintes cumulatives appliquées au satellite.

La vérification de ladite au moins une contrainte cumulative consiste à valider que la séquence optimale d'acquisition obtenue en sortie de l'étape 500 est compatible avec le fonctionnement général du satellite. A cet effet, cette vérification est mise en œuvre au cours de l'étape 600 au moyen, par exemple, d'un processus itératif consistant à supprimer une à une de ladite séquence optimale les opportunités d'acquisition discrètes par ordre de priorité croissant tant que ladite au moins une contrainte cumulative n'est pas satisfaite. De cette manière, le procédé produit en sortie de l'étape 600 une séquence qui est une sous séquence de la séquence optimale obtenue en sortie de l'étape 500, et qui est compatible avec les contraintes cumulatives associées au satellite. Cette séquence obtenue en sortie de l'étape 600 est le plan d'acquisition effectivement suivi par le satellite au cours de sa mission.

La figure 3 représente un mode particulier de mise en œuvre du procédé de la figure 1 au cours duquel le procédé peut comporter une étape 50 d'agrégation de zones terrestres, antérieure à l'étape 100 de détermination d'opportunités d'acquisition discrètes.

Par agrégation de zones terrestres, on entend ici le fait de regrouper certaines zones terrestres Z-i , ...,Z _N , au sein d'une ou plusieurs zones dites englobantes, lesdites zones regroupées devant être proches géographiquement les unes des autres, de sorte à maximiser le nombre de zones acquises au cours d'une même opportunité d'acquisition discrète.

Par exemple, une zone engloblante peut contenir plusieurs zones terrestres alignées dans la direction d'acquisition du satellite le long de son orbite.

Dans un autre exemple, éventuellement couplé au précédent, une zone englobante peut contenir plusieurs zones terrestres alignées dans une direction transverse à l'orbite du satellite. Dans ce cas, la taille de ladite zone englobante est bornée par la taille de la fauchée du satellite.

Dans la mise en œuvre de l'étape 50, chaque zone englobante est par exemple associée à une priorité égale à la priorité de plus haute importance parmi les priorités de toutes les zones regroupées au sein de ladite zone englobante. De plus, le critère qualité d'une zone englobante est par exemple égal à la somme des critères qualité respectifs des zones regroupées au sein de ladite zone englobante.

A l'issue de l'étape 50, chaque zone englobante est vue comme une zone terrestre à part entière de sorte que les étapes 100 à 600 ne sont pas modifiées.

Il est à noter que le regroupement de zones terrestres au sein de zones englobantes se révèle pertinent lorsque plusieurs zones terrestres sont contigues, de sorte à pouvoir être toutes acquises en une seule fois. Une telle manière de procéder est avantageuse car elle permet au satellite d'économiser le temps de manœuvre nécessaire autrement pour se déplacer entre les opportunités d'acquisition discrètes associées respectivement à chacune desdites zones terrestres regroupées. Cette économie de temps peut avantageusement être utilisée pour acquérir davantage de zones terrestres qui ne peuvent quant à elles pas être regroupées au sein de zones englobantes.

La figure 4 représente un mode préférentiel de mise en œuvre du procédé de la figure 1 au cours duquel l'étape 500 de détermination d'une séquence optimale comporte avantageusement la mise à jour dynamique d'une variable pathMesh stockant les zones Z, associées aux opportunités d'acquisition discrètes de ladite séquence optimale, de sorte à empêcher la duplication d'opportunités d'acquisition discrètes de zones devant être acquises une seule fois.

A cet effet, en notant mesh(D') la fonction qui renvoie la zone terrestre à laquelle est associée l'opportunité d'acquisition discrète D', le premier algorithme se réécrit de la façon suivante :

Pour i allant de 1 à Q(D) ; obj(Û) = c[Û] Fin de Pour ;

Pour i allant de 1 à Q(D) ; previous(D') = - 1 Fin de Pour ;

Pour i allant de 1 à Q(D) :

Pour j = maxNextVertice(Û) :

Si M ^J = 1 :

Si obj(Û) + c[0] > obj(Û) :

Si mesh(Û) g pathMesh(d) ;

previous(Û) = Û ;

obj(O) = obj(Û) + c[Û] ;

pathMesh(Û) = {pathMesh(Û) mesh(Û)} ;

Fin de Si ;

Fin de Si ;

Fin de Si ;

Fin de Pour ;

Fin de Pour ; En stockant dans la variable pathMesh l'historique des zones terrestres visitées lors de la détermination de la séquence optimale, le procédé est avantageusement adapté à empêcher qu'une zone terrestre Z, ayant un critère qualité élevé, associée à une requête de haute priorité ainsi que devant être acquise une unique fois, soit finalement acquise une pluralité de fois. En effet, dans l'hypothèse où une telle zone Z, est associée à plusieurs opportunités d'acquisition discrètes compatibles cinématiquement entre elles, la méthode de calcul de plus long chemin cherchant à maximiser le poids de la séquence optimale, il est tout à fait possible que plusieurs desdites opportunités d'acquisition discrètes associées à Z, soient inclues dans ladite séquence optimale. Il est donc avantageux d'introduire la variable pathMesh tel qu'explicité ci-dessus afin de discriminer ce type de configurations.

Les caractéristiques détaillées précédemment permettent d'obtenir une séquence optimale garantissant une prise de vue unique pour chaque zone terrestre associée à une requête monovue. Reste alors à considérer également le cas de zones terrestres associées à des requêtes stéréo/tri- stéréo.

C'est pourquoi, dans ledit mode préférentiel de mise en œuvre illustré par la figure 4, le procédé comporte aussi une étape 450 de prise en compte de requêtes d'acquisitions multiples, consécutive à l'étape 400 d'évaluation de compatibilité cinématique et antérieure à l'étape 500 de détermination d'une séquence optimale, dans laquelle les poids respectifs des opportunités d'acquisition discrètes associées à des zones Z, devant être acquises une pluralité de fois sont mis à jour et dans laquelle une variable mustHit force l'enchaînement desdites opportunités d'acquisition discrètes afin de satisfaire lesdites requêtes d'acquisitions multiples.

Il est à noter que les opportunités d'acquisition discrètes associées à une zone Z, devant être acquise une pluralité de fois appartiennent au même intervalle temporel d'accessibilité visuelle T, et sont en outre distinctes les unes des autres. Lesdites opportunités ont donc des dates de début différentes et permettent ainsi par exemple d'acquérir ladite zone Z, sous un angle différent.

Le principe général de l'étape 450 consiste à augmenter artificiellement du double ou du triple le poids des opportunités d'acquisition discrètes associées à des zones Z, devant être acquises respectivement en mode stéréo ou tri-stéréo. On comprend par ailleurs que l'agilité du satellite permet d'envisager de satisfaire des requêtes monovues entre des opportunités d'acquisition associées à des zones Z, devant être acquises plusieurs fois.

Pour cela, on détermine dans un premier temps, pour chaque requête correspondant à une demande d'acquisitions multiples, l'ensemble des séquences d'acquisition compatibles cinématiquement, comprenant exactement le nombre d'acquisitions requis de sorte à satisfaire ladite requête (ce nombre est de 2 et 3 respectivement dans les cas stéréo et tri-stéréo), ainsi qu'éventuellement des opportunités d'acquisition discrètes associées à d'autres requêtes monovues tel qu'évoqué ci-dessus. Une telle séquence est dite séquence multi-acquisition et est construite de sorte que ses premier et dernier éléments sont des opportunités d'acquisition discrètes appartenant à une même zone terrestre et associées à la requête dont dépend ladite séquence multi-acquisition. Autrement dit, lesdits premier et dernier éléments ne sont pas associés à des requêtes monovues.

Dans un deuxième temps, pour chaque séquence multi-acquisition telle que déterminée ci-dessus :

- on détermine un critère de chaînage. Ledit critère de chaînage est un scalaire qui correspond à une mesure de l'intérêt de satisfaire la requête à laquelle est associée ladite séquence multi-acquisition selon effectivement ladite séquence multi-acquisition plutôt qu'une autre. Le critère de chaînage est donc une échelle de comparaison des séquences multi-acquisition associées à une même requête,

- on détermine un poids de chaînage. Dans un mode particulier de mise en œuvre, ledit poids de chaînage dépend de la priorité de la requête, du critère qualité w[Z,] de la zone Z, associée à ladite requête ainsi que dudit critère de chaînage. Par exemple, la prise en compte dudit critère de chaînage dans le poids de chaînage se fait en additionnant ledit critère de chaînage au critère qualité w[Z,] de manière identique à ce qui est décrit au cours de l'étape 100 concernant l'utilisation du critère d'intérêt d'acquisition w[a ]. Dès lors, le poids de chaînage de la séquence multi-acquisition est un vecteur dont une seule composante est non nulle,

- le graphe G est mis à jour selon un graphe G _mu iti- D'une part, le graphe G _mu iti comporte des sommets correspondant respectivement aux opportunités d'acquisition discrètes associées à des requêtes monovues et ayant respectivement un poids tel que déterminé dans l'étape 100. D'autre part, le graphe G _mu iti comporte des sommets correspondant respectivement à chaque opportunité d'acquisition discrète de ladite séquence multi-acquisition et ayant tous un poids nul, sauf le sommet associé à la première composante de ladite séquence et dont le poids est égal au poids de chaînage,

- on crée une variable mustHit qui à chaque opportunité d'acquisition discrète de ladite séquence multi-acquisition associe l'opportunité d'acquisition discrète qui lui succède dans le temps au sein de ladite séquence multi-acquisition. Par exemple, à titre nullement limitatif, si ladite séquence multi-acquisition est associée à une requête tri-stéréo d'une zone Z, et s'écrit (a- , af, af), alors on a mustHit(al) = af et mustHit(af) = af . Dans un autre exemple, et à titre nullement limitatif, si ladite séquence multi-acquisition est associée à une requête tri-stéréo d'une zone Z, et s'écrit af, a , af), où aj est une opportunité associée à une requête monovue d'une zone Z _j différente de alors on a mustHit(al) = af , mustHit(af ) = aj et mustHit(aj) = af .

Il est à noter que dans le cas où une requête correspondant à une demande d'acquisitions multiples est associée à plusieurs séquences multi- acquisition possibles, le graphe G _mu iti comporte des sommets correspondant à des opportunités d'acquisition discrètes identiques en termes de date de début, mais néanmoins différentes les unes des autres en termes de poids. Par exemple, et à titre nullement limitatif, supposons qu'une requête stéréo soit associée à un intervalle d'accessibilité discrétisé selon quatre opportunités d'acquisition discrètes (a} , a , af , af de sorte que seules les séquences multi- acquisition (a- , a ) et (a- , a ) sont cinématiquement compatibles et adaptées à satisfaire ladite requête. Alors le sommet a} va apparaître deux fois dans le graphe G _mU | _t i, à savoir une première fois avec un poids égal au poids de chaînage de la séquence (a}, af et de sorte que = af , ainsi qu'une deuxième fois avec un poids égal au poids de chaînage de la séquence (a- , a ) et de sorte que = af. Les sommets af, af et af n'apparaissent quant à eux qu'une seule fois dans le graphe G _mu iti et ont un poids nul.

Le graphe G _mu iti est représenté au moyen d'une matrice M _mu iti selon un mode de mise en œuvre identique à celui utilisé pour représenter le graphe G au moyen de la matrice M dans l'étape 400. On comprend donc, du fait que certaines opportunités d'acquisition discrètes puissent apparaître plusieurs fois dans le graphe G _mu iti tel que décrit ci-avant, que ladite matrice M _mu iti a une taille supérieure à la matrice M. En outre le graphe G _mu iti est celui utilisé lors de l'étape 500 de recherche d'une séquence optimale lorsque l'ensemble des requêtes transmises au satellite comporte au moins une requête d'acquisitions multiples.

Par ailleurs, la variable mustHit est aussi adaptée à forcer l'exécution d'une séquence multi-acquisition une fois que le premier élément de cette dernière est effectivement intégré à la séquence optimale au cours de l'étape 500. Une telle façon de procéder est avantageuse car elle permet au procédé de prendre en compte dans leur entièreté les requêtes associées à des demandes d'acquisitions multiples sans pour autant privilégier ces dernières au détriment de requêtes monovues comportant des opportunités d'acquisition discrètes de poids respectifs supérieurs des opportunités que comprennent lesdites requêtes d'acquisitions multiples.

La variable mustHit est aussi avantageusement utilisée au cours du calcul du plus long chemin de l'étape 500 afin de vérifier que lorsqu'une opportunité d'acquisition discrète est ajoutée à la séquence optimale en construction, cette opportunité d'acquisition discrète est bien stockée dans la variable mustHit appliquée aux éléments de la séquence optimale précédant ladite opportunité ajoutée.

Dans le cas d'une constellation de satellites, il est avantageux de pouvoir répondre rapidement à un client si une requête peut-être réalisée et sous quel délai. Dans ce cas le procédé peut être :

a) réception de la requête du client,

b) identification du prochain satellite qui peut faire cette acquisition, sans prendre en compte le plan d'acquisition en cours de ce satellite

c) génération d'un nouveau plan pour ledit satellite avec les contraintes suivantes :

d . toutes les acquisitions programmées dans le plan courant doivent le rester (sauf si l'on accepte d'annuler des acquisitions prévues), par contre leur conditions d'acquisition peuvent changer.

c2. insertion de la nouvelle requête

d) en fonction du résultat de l'optimisation :

d1 . si l'acquisition est incluse dans le nouveau plan alors le plan modifié est transmis au satellite

d2. sinon on réitère sur le satellite suivant

Il n'est pas nécessaire de prendre les nouvelles demandes client une par une, cela peut se faire de façon périodique, par exemple toutes les 5 minutes. Moyennant un choix adapté des priorités, il est possible de forcer la prise en compte de la nouvelle requête.

De part sa rapidité le calcul du plan peut être fait fréquemment. En outre, la mise à jour de la solution est moins chère qu'un calcul complet en partant des résultats du calcul précédent, et notamment la matrice M représentant le graphe obtenu à la fin de l'étape 400 ou le graphe G _mU | _t associé à la matrice M _mu iti à la fin de l'étape 450

Ainsi il est possible de minimiser le temps de calcul pour les étapes antérieures à l'étape 500, qui sont les plus gourmandes en temps de calcul.