FR 3039728 A1 20170203 FR 1557215 A 20150728 3039728 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 15 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 20 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 25 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 30 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 2 3039728 terrestre devant être acquise, comme par exemple des conditions de localisation, d'exposition lumineuse, ou encore de nécessité d'acquisitions 5 10 20 25 30 3 3039728 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 5 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). 10 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 15 ê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 4 3039728 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 5 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 10 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 Z1,...,ZN, par un engin spatial en 15 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 Ri de sorte qu'à ladite zone Z correspond un intervalle temporel d'accessibilité visuelle Ti le long de ladite orbite défilante, 1 5. i N. En outre, ledit procédé comporte les étapes successives suivantes : 20 une étape de détermination pour chaque intervalle Ti d'opportunités d'acquisition discrètes de la zone Zi de sorte qu'à chaque dite opportunité d'acquisition discrète sont associées une date de début incluse dans ledit intervalle Ti, une durée d'exécution, une contrainte locale cinématique ainsi qu'un poids, 25 une étape de regroupement des opportunités d'acquisition discrètes des zones Zi 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', 30 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 5 3039728 d'acquisition discrètes ayant un poids maximal, et compatible cinématique ment. Dans des modes particuliers de mise en oeuvre, le procédé de planification d'acquisition d'images de zones terrestres comporte l'une ou 5 plusieurs des caractéristiques suivantes, prises isolément ou selon toutes les combinaisons techniquement possibles. Dans un mode particulier de mise en oeuvre, 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 10 détermination d'une séquence optimale, dans laquelle les poids respectifs des opportunités d'acquisition discrètes associées à des zones Zi 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. 15 Dans un mode particulier de mise en oeuvre, l'étape de détermination d'une séquence optimale comporte la mise à jour dynamique d'une variable pathMesh stockant les zones Zi 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. 20 Dans un mode particulier de mise en oeuvre, l'étape de détermination d'une séquence optimale comporte pour chaque opportunité d'acquisition discrète Di de l'ensemble D une sélection partielle d'opportunités d'acquisition discrètes précédant Di dans D, et compatibles cinématiquement avec Di, parmi toutes les opportunités d'acquisition discrètes précédant Di dans D, et 25 compatibles cinématiquement avec D. Dans un mode particulier de mise en oeuvre, 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. 30 Dans un mode particulier de mise en oeuvre, ladite au moins une contrainte cumulative comporte une contrainte de taille mémoire. Dans un mode particulier de mise en oeuvre, ladite au moins une contrainte cumulative comporte une contrainte de puissance électrique 6 3039728 consommée. Dans un mode particulier de mise en oeuvre, ladite au moins une contrainte cumulative comporte une contrainte de temps maximum de fonctionnement d'instruments embarqués dans ledit engin. 5 Dans un mode particulier de mise en oeuvre, 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 oeuvre, les dates de début des 10 opportunités d'acquisition discrètes de chaque intervalle Ti sont échantillonnées temporellement selon un pas de temps constant. Dans un mode particulier de mise en oeuvre, les pas de temps respectifs des intervalles T1,...,TN sont égaux entre eux. Dans un mode particulier de mise en oeuvre, les pas de temps 15 respectifs d'au moins deux intervalles Ti et Ti sont différents. Dans un mode particulier de mise en oeuvre, le pas de temps de chaque intervalle Ti est ajusté de sorte que ledit intervalle Ti est associé à 20 opportunités d'acquisition discrètes. Dans un mode particulier de mise en oeuvre, la contrainte locale 20 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 oeuvre, le poids de chaque opportunité d'acquisition discrète associée à une zone Zi dépend de la priorité 25 de la requête Ri associée à ladite zone Z. Dans un mode particulier de mise en oeuvre, chaque zone Zi est associée à un critère qualité wi dépendant des caractéristiques géométriques de ladite zone Zi ainsi que des conditions météorologiques prévues au cours de l'intervalle temporel d'accessibilité visuelle Ti. 30 Dans un mode particulier de mise en oeuvre, le poids de chaque opportunité d'acquisition discrète associée à une zone Zi dépend de la priorité de la requête 1=1; de sorte à être un vecteur dont la seule composante non nulle est égale au critère qualité wi, la position de ladite composante non nulle dans 7 3039728 ledit vecteur dépendant de la priorité de ladite requête R1. Dans un mode particulier de mise en oeuvre, ladite composante non nulle dépend du critère qualité %Ni et d'un critère d'intérêt d'acquisition associé à ladite opportunité d'acquisition discrète. 5 Dans un mode particulier de mise en oeuvre, 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 oeuvre, le procédé comporte une étape d'agrégation de zones terrestres, antérieure à l'étape de détermination 10 d'opportunités d'acquisition discrètes. Dans un mode particulier de mise en oeuvre, 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 15 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 oeuvre 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 20 mise en oeuvre d'un procédé de planification d'acquisition d'images de zones terrestres Z1,...,4 par un engin spatial. Figure 2: une représentation d'un mode préférentiel de mise en oeuvre 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 25 contrainte cumulative imposée audit engin spatial. Figure 3: une représentation d'un mode particulier de mise en oeuvre 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 30 oeuvre 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 8 3039728 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. 5 La figure 1 représente un organigramme d'un exemple de mise en oeuvre d'un procédé de planification d'acquisition d'images de zones terrestres Z1,...,4 par un engin spatial. Par « acquisition d'images », on entend ici la mesure d'un rayonnement électromagnétique reçu desdites zones terrestres Z1,...,4. A cet 10 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. 15 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 « Medium Earth Orbit »). Par exemple, ledit satellite d'observation agile suit 20 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,...) 25 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 Zi,...,ZN sont des surfaces au sol, d'aires respectives finies, éventuellement contigües et intersectent la fauchée dudit 30 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. 9 3039728 Dans la suite de la description, chaque zone terrestre Zi est associée à une requête Ri émise par un utilisateur. Chaque requête Ri pouvant être elle-même associée à différents paramètres : une priorité p,, pl étant par exemple un entier naturel, qui 5 caractérise la nécessité de satisfaire cette requête Ri avant éventuellement d'autres requêtes. Chaque priorité pi dépend de plusieurs facteurs comme par exemple l'importance et le coût de la requête Ri de laquelle elle dépend, des contraintes d'acquisition dépendant de plusieurs facteurs, 10 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 Zi associée à ladite 15 requête Ri, 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, 20 est à noter que d'une manière générale, c'est-à-dire au moment où une requête Ri est signifiée par un utilisateur, la zone Zi associée à ladite requête Ri 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 25 consistant à découper ladite zone Zi 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 Zi 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 Zi 30 associées aux requêtes Ri 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. Il 10 3039728 Par ailleurs, à chaque zone terrestre Zi correspond un intervalle temporel d'accessibilité visuelle Ti le long de ladite orbite défilante. Ledit intervalle Ti 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 5 la zone Zi tout en respectant les contraintes d'acquisition ainsi que le type d'acquisition de la requête 1=1; à 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 = (V1, 10 V2, ..., VN). De la même manière, un élément positionné à la i-ème ligne et j- ème colonne d'une matrice M est noté Mil. 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 15 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 20 fonctionnement dudit satellite. A cet effet, le procédé comporte dans un premier temps une étape 100 de détermination, pour chaque intervalle Tb d'opportunités d'acquisition discrètes ail, al, ... de la zone Zi de sorte qu'à chaque dite opportunité d'acquisition discrète sont associées une date de début incluse dans ledit 25 intervalle Tb une durée d'exécution, une contrainte locale cinématique ainsi qu'un poids. Une opportunité d'acquisition discrète al 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 Zi, à partir d'une date de début ainsi que pendant une durée 30 d'exécution comprises dans l'intervalle Ti de ladite zone Zi, de sorte que soient respectées les contraintes d'acquisition ainsi que le type d'acquisition de la requête Ri associée à cette zone Z. 11 3039728 Dans la suite de la description, une séquence d'acquisition désigne une suite d'opportunités d'acquisition discrètes. Un intervalle Ti comporte une ou plusieurs opportunités d'acquisition discrètes qui discrétisent ce dernier selon leurs dates de début respectives. En 5 effet, l'agilité du satellite autorise le fait qu'une zone Zi puisse être acquise potentiellement en des temps distincts de son intervalle Ti. Dans un mode particulier de mise en oeuvre de l'étape 100, les dates de début des opportunités d'acquisition discrètes de chaque intervalle Ti sont échantillonnées temporellement selon un pas de temps constant. Il est à noter 10 que plus ledit pas de temps constant est choisi petit, plus il y a d'opportunités d'acquisition au sein de chaque intervalle Th 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 oeuvre de l'étape 100, les pas de 15 temps respectifs d'au moins deux intervalles Ti et Ti sont différents. Une telle configuration est avantageuse lorsque par exemple les priorités respectives des requêtes Ri et Ri associées auxdits intervalles Ti et Ti 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 20 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 25 d'acquisition discrètes au sein dudit intervalle d'accessibilité. Dans un mode particulier de mise en oeuvre de l'étape 100, les pas de temps respectifs des intervalles T1,...,TN 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 R1,...,RN 30 associées aux intervalles T1,...,TN sont toutes identiques. Dans un mode particulier de mise en oeuvre de l'étape 100, le pas de temps de chaque intervalle Ti est ajusté de sorte que ledit intervalle Ti est associé à un nombre prédéfini d'acquisition discrètes. Par exemple, un nombre 12 303 9 72 8 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 5 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 10 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. 15 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. 20 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. 25 Le poids c[al] d'une opportunité d'acquisition discrète al 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 30 associée à une zone Zi dépend de la priorité pi de la requête 1=1; associée à ladite zone Zi, de sorte que plus ladite priorité est élevée, plus les poids respectifs desdites opportunités d'acquisition discrètes sont importants. 13 3039728 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 oeuvre de l'étape 100, chaque zone Zi est associée à un critère qualité w[Zi] dépendant 5 des caractéristiques géométriques de ladite zone Zi 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 10 météorologiques rencontrées pendant l'intervalle temporel d'accessibilité visuelle Ti associé à Z sont sensiblement identiques à celles exigées par les contraintes d'acquisition de la requête reliée audit intervalle Ti. Dans un mode particulier de mise en oeuvre de l'étape 100, le poids de chaque opportunité d'acquisition discrète associée à une zone Zi dépend de la 15 priorité pi de la requête 1=1; de sorte à être un vecteur dont la seule composante non nulle est égale audit critère qualité w[Zi], la position de ladite composante non nulle dans ledit vecteur étant la priorité pi de ladite requête R. Dans ce mode de mise en oeuvre, 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 20 mission. En outre, on adopte la convention qu'une priorité pi est de plus grande importance qu'une autre priorité pi si et seulement si pi opérant de la façon suivante : on note c [a] > c[b] au sens que le poids d'une opportunité 30 d'acquisition discrète a est supérieur au poids d'une opportunité d'acquisition discrète b si et seulement si {cd[a] > cd[b]} {(c1 [a] = cl [b]) A (c2 [a] > c2 [b])} , 14 3039728 les opérateurs A et y 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. 5 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. 10 Dans un mode particulier de mise en oeuvre de l'étape 100, ladite composante non nulle du poids d'une opportunité d'acquisition discrète al dépend du critère qualité w[Zi] ainsi que d'un critère d'intérêt d'acquisition w[al associé à ladite opportunité d'acquisition discrète al. Ledit critère d'intérêt d'acquisition w[ah est un scalaire qui mesure l'intérêt d'acquérir plus 15 particulièrement l'opportunité al 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 oeuvre, le poids d'une opportunité d'acquisition discrète est mis à jour de sorte qu'en reprenant 20 l'exemple ci-avant : c[al] = (0, w[Zi] + w[ctij, 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 25 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[al] 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 30 opportunités d'acquisition discrètes des différentes zones Zi dans un ensemble D. 15 303 9 72 8 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 5 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. 10 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 15 elles sont associées. Toutefois, dans le présent mode de mise en oeuvre, 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 20 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. 25 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 Di et Di de l'ensemble D, telles que ladite opportunité Di 30 ait une date de début antérieure à la date de début de ladite opportunité Dl (donc telles que i j). Cela revient à dire que les éléments de toute séquence d'acquisition compatible cinématiquement sont deux à deux 20 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 25 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 30 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 17 303 9 72 8 graphe G, alors qu'il est de complexité exponentielle pour un graphe quelconque. Dans le présent mode de mise en oeuvre, le graphe G est représenté au moyen d'une matrice carrée M de booléens, dite matrice d'adjacence, dont 5 le terme situé à la ligne i et à la colonne j est : - Mil = 1 si Di et Dl sont compatibles cinématiquement, - = 0 sinon. Les conditions de compatibilité entre deux opportunités d'acquisition discrètes impliquent que M est une matrice triangulaire supérieure stricte. 10 Il est à noter que la mise en oeuvre 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 oeuvre est avantageuse dans la mesure où les tâches effectuées pendant les étapes 100 à 400, et plus particulièrement 15 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 oeuvre des étapes 100 à 400 permet avantageusement de 20 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 25 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. 30 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. 18 3039728 Pour ce faire, on procède en deux temps. Dans un premier temps, on détermine pour chaque opportunité d'acquisition discrète Di 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 5 cela, on introduit un ensemble P(Di) de toutes les opportunités d'acquisition discrètes précédant Di dans l'ensemble D, et compatibles cinématiquement avec D. Alors le poids objectif obj(D') est égal au poids de Di auquel s'ajoute le maximum des poids objectifs des opportunités d'acquisition discrètes de P(D'), c'est-à-dire 10 obj(D') = c[111 + maxsep(Di )(obj(s)). En faisant désormais l'hypothèse que l'ensemble D est de cardinal 0(D), et en introduisant une variable previous(D5 qui stocke l'opportunité d'acquisition discrète de l'ensemble P(Di) maximisant le poids de la séquence optimale menant à Di, la détermination du poids objectif de tous les éléments 15 de l'ensemble D se fait dynamiquement au moyen d'un premier algorithme s'écrivant : Pour i allant de 1 à 0(D) : obj(D') = c[d] Fin de Pour ; Pour i allant de 1 à 0(D) : previous(D') = -1 Fin de Pour ; 20 Pour i allant de 1 à 0(D) : Pour j allant de (i+1) à 0(D) : Si el = / : Si obj(D') + c[D1] > obj(D') : previous(D9 = D' ; 25 obje = obj(D') + c[d] ; Fin de Si ; Fin de Si ; Fin de Pour ; Fin de Pour ; 30 Il 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, 19 3039728 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 5 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 10 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 15 sens de la maximalité du poids. Dans un mode particulier de mise en oeuvre du procédé, l'étape 500 comporte pour chaque opportunité d'acquisition discrète Di de l'ensemble D une sélection partielle d'opportunités d'acquisition discrètes précédant Di dans l'ensemble D, et compatibles cinématiquement avec D', parmi toutes les 20 opportunités d'acquisition discrètes précédant Di dans l'ensemble D, et compatibles cinématiquement avec D. Une telle manière de procéder revient à restreindre l'ensemble P(Di) pour chaque opportunité d'acquisition discrète Di de l'ensemble D. Ainsi, ledit ensemble P(Di) peut être avantageusement remplacé par un ensemble Pc(D5 25 défini comme un sous-ensemble de ce dernier, et comportant les opportunités d'acquisition discrètes précédant Di dans l'ensemble D, compatibles cinématiquement avec Di ainsi que proches de D. Par proches de Di, on entend ici que les éléments de Pc(D5 occupent respectivement une position dans D de sorte que la différence entre cette position et celle de Di reste 30 bornée par une constante prédéfinie. Dès lors, le poids objectif obj(Di) s'obtient par la formule obj(Di) = c[Dil + maxsepaDi)(obj(s)). 20 3039728 A cet effet, on introduit une variable maxNextVertice(Di) qui contient pour chaque opportunité d'acquisition Di le nombre maximum de sommets du graphe G visités au cours de la recherche de la séquence optimale. Le scalaire maxNextVertice(Di) constitue ainsi un majorant du cardinal de l'ensemble 5 Pc(Di). Dès lors, le premier algorithme devient : Pour i allant de 1 à 0(D) : obj(D') = c[d] Fin de Pour ; Pour i allant de 1 à 0(D) : previous(D') = -1 Fin de Pour ; Pour i allant de 1 à 0(D) : 10 Pour j allant de (i+1) à maxNextVertice(D9 : Si Mid = 1 : Si obj(D9 + c[D1] > obj(D9 : previous(D') = D' ; obje = obj(D9 + c[d] ; 15 Fin de Si ; Fin de Si ; Fin de Pour ; Fin de Pour ; 20 Limiter l'ensemble P(Di) à l'ensemble Pc(D5 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 Pc(Di) contient uniquement des opportunités d'acquisition discrètes séparées de soixante secondes au plus de Di le long de l'orbite 25 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ù Pc(d) n'est pas utilisé. La figure 2 représente un mode préférentiel de mise en oeuvre 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 30 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é 21 303 9 72 8 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. 5 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. 10 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 15 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. 20 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 25 vérification est mise en oeuvre 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 30 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. 22 303 9 72 8 La figure 3 représente un mode particulier de mise en oeuvre 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. 5 Par agrégation de zones terrestres, on entend ici le fait de regrouper certaines zones terrestres Zi,...,ZN, 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. 10 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 15 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 oeuvre 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 20 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 25 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 contigües, 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 30 le temps de manoeuvre 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 23 3039728 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 oeuvre 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 5 variable pathMesh stockant les zones Zi 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(Di) la fonction qui renvoie la zone terrestre 10 à 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 à 0(D) : obj(D') = el Fin de Pour ; Pour i allant de 1 à 0(D) : previous(D') = -1 Fin de Pour ; 15 Pour i allant de 1 à 0(D) : Pour j = (i+1)...maxNextVertice(Di) : Si Mid = 1 : Si obj(Di) + c[D1 > obj(Di) : Si mesh(d) E pathMesh(d) ; 20 previouse = d ; obj(d) = obj(Di) + c[d] ; pathMesh(DI) = {pathMesh(Di) meshel ; Fin de Si ; Fin de Si ; 25 Fin de Si ; Fin de Pour ; Fin de Pour ; En stockant dans la variable pathMesh l'historique des zones 30 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 Zi 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 24 303 9 72 8 effet, dans l'hypothèse où une telle zone Zi 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 5 opportunités d'acquisition discrètes associées à Zi 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 10 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/tristéréo. C'est pourquoi, dans ledit mode préférentiel de mise en oeuvre illustré par la figure 4, le procédé comporte aussi une étape 450 de prise en compte 15 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 Zi devant être acquises une pluralité de fois sont mis à jour et dans laquelle une variable mustHit force 20 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 Zi devant être acquise une pluralité de fois appartiennent au même intervalle temporel d'accessibilité visuelle Ti et sont en outre distinctes les unes 25 des autres. Lesdites opportunités ont donc des dates de début différentes et permettent ainsi par exemple d'acquérir ladite zone Zi 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 Zi devant être acquises respectivement en 30 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 Zi devant être acquises plusieurs fois. 25 303 9 72 8 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 5 (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 à 10 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 : 15 - 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 20 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 oeuvre, ledit poids de chaînage dépend de la priorité de la requête, du critère qualité w[Zi] de la zone Zi associée à ladite requête ainsi que dudit critère de chaînage. Par exemple, la prise 25 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[Zi] 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[ah. Dès lors, le poids de chaînage de la séquence multi-acquisition est un 30 vecteur dont une seule composante est non nulle, - le graphe G est mis à jour selon un graphe Gmuiti. D'une part, le graphe Gmuiti comporte des sommets correspondant 26 3039728 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 Gmei comporte des sommets correspondant respectivement à chaque opportunité 5 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 10 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 Zi et s'écrit (ail:, al, cti3.), alors on a 15 mustHit(ah = al et mustHit(al) = 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 (ail:, al, agi, cti3.), où aj est une opportunité associée à une requête monovue d'une zone Zi différente de Z, alors on a 20 mustHit(ah = al , mustHit(al) = aj et mustHit(aJ) = ai3 Il est à noter que dans le cas où une requête correspondant à une demande d'acquisitions multiples est associée à plusieurs séquences multiacquisition possibles, le graphe Gmei comporte des sommets correspondant à des opportunités d'acquisition discrètes identiques en termes de date de début, 25 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 (ail:, al, al, ah de sorte que seules les séquences multiacquisition (ail:, ah et (ail ah sont cinématiquement compatibles et adaptées à 30 satisfaire ladite requête. Alors le sommet ai': va apparaître deux fois dans le graphe Gmei, à savoir une première fois avec un poids égal au poids de chaînage de la séquence (ail:, ah et de sorte que mustHit(al) = al , ainsi qu'une 27 303 9 72 8 deuxième fois avec un poids égal au poids de chaînage de la séquence (ail, ah et de sorte que mustHit(ail) = ai,. Les sommets al, al et ai, n'apparaissent quant à eux qu'une seule fois dans le graphe Gmei et ont un poids nul. 5 Le graphe Gmei est représenté au moyen d'une matrice Mmulti selon un mode de mise en oeuvre 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 Gmei tel que décrit ci-avant, que ladite matrice Mmulti a une taille 10 supérieure à la matrice M. En outre le graphe Gmuiti 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 15 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 20 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 25 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 30 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, 28 3039728 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 : c1. toutes les acquisitions programmées dans le plan courant doivent 5 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 10 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 15 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 Gme associé 20 à la matrice Mmulti à 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. REVENDICATIONS1. Procédé de planification d'acquisition d'images de zones terrestres Zi,...,ZN, 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 Zi étant associée à une requête 1=1; de priorité pi de sorte qu'à ladite zone Zi correspond un intervalle temporel d'accessibilité visuelle Ti le long de ladite orbite défilante, caractérisé en ce qu'il comporte les étapes successives suivantes : - une étape (100) de détermination pour chaque intervalle Ti d'opportunités d'acquisition discrètes de la zone Zi de sorte qu'à chaque dite opportunité d'acquisition discrète sont associées une date de début incluse dans ledit intervalle Tb une durée d'exécution, une contrainte locale cinématique ainsi qu'un poids, - une étape (200) de regroupement des opportunités d'acquisition discrètes des zones Zi dans un ensemble D, - une étape (300) 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 Di, - 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, - une étape (500) de détermination d'une séquence optimale d'opportunités d'acquisition discrètes ayant un poids maximal, et compatible cinématiquement. 2. Procédé selon la revendication 1, comportant 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 Zi 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. 3039728 29 3. Procédé selon l'une des revendications 1 à 2, dans lequel l'étape (500) de détermination d'une séquence optimale comporte la mise à jour dynamique d'une variable pathMesh stockant les zones Zi associées aux opportunités d'acquisition discrètes de ladite séquence optimale, de sorte à 5 empêcher la duplication d'opportunités d'acquisition discrètes de zones devant être acquises une seule fois. 4. Procédé selon l'une des revendications 1 à 3, dans lequel l'étape (500) comporte pour chaque opportunité d'acquisition discrète Di de l'ensemble D une sélection partielle d'opportunités d'acquisition discrètes précédant D' 10 dans D, et compatibles cinématiquement avec Di, parmi toutes les opportunités d'acquisition discrètes précédant Di dans D, et compatibles cinématiquement avec D'. 5. Procédé selon l'une des revendications 1 à 4, comportant une étape (600) de vérification de respect d'au moins une contrainte cumulative 15 imposée audit engin spatial, ladite étape (600) de vérification étant consécutive à l'étape (500) de détermination d'une séquence optimale. 6. Procédé selon la revendication 5, dans lequel ladite au moins une contrainte cumulative comporte une contrainte de taille mémoire. 7. Procédé selon l'une des revendications 5 à 6, dans lequel ladite 20 au moins une contrainte cumulative comporte une contrainte de puissance électrique consommée. 8. Procédé selon l'une des revendications 5 à 7, dans lequel ladite au moins une contrainte cumulative comporte une contrainte de temps maximum de fonctionnement d'instruments embarqués dans ledit engin. 25 9. Procédé selon l'une des revendications 5 à 8, dans lequel l'étape (600) 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. 30 10. Procédé selon l'une des revendications 1 à 9, dans lequel les dates de début des opportunités d'acquisition discrètes de chaque intervalle Ti sont échantillonnées temporellement selon un pas de temps constant. 11. Procédé selon la revendication 10, dans lequel les pas de temps 3039728 30 respectifs des intervalles T1,...,TN sont égaux entre eux. 12. Procédé selon la revendication 10, dans lequel les pas de temps respectifs d'au moins deux intervalles Ti et Ti sont différents. 13. Procédé selon la revendication 10, dans lequel le pas de temps 5 de chaque intervalle Ti est ajusté de sorte que ledit intervalle Ti est associé à 20 opportunités d'acquisition discrètes. 14. Procédé selon l'une des revendications 1 à 13, dans lequel 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 10 opportunité d'acquisition discrète. 15. Procédé selon l'une des revendications 1 à 14, dans lequel le poids de chaque opportunité d'acquisition discrète associée à une zone Zi dépend de la priorité pi de la requête Ri associée à ladite zone Z. 16. Procédé selon l'une des revendications 1 à 15, dans lequel 15 chaque zone Zi est associée à un critère qualité wi dépendant des caractéristiques géométriques de ladite zone Zi ainsi que des conditions météorologiques prévues au cours de l'intervalle temporel d'accessibilité visuelle Ti. 17. Procédé selon la revendication 16, dans lequel le poids de 20 chaque opportunité d'acquisition discrète associée à une zone Zi dépend de la priorité pi de la requête Ri de sorte à être un vecteur dont la seule composante non nulle est égale au critère qualité wi, la position de ladite composante non nulle dans ledit vecteur dépendant de la priorité de ladite requête R. 18. Procédé selon la revendication 17, dans lequel ladite composante 25 non nulle dépend du critère qualité wi et d'un critère d'intérêt d'acquisition associé à ladite opportunité d'acquisition discrète. 19. Procédé selon l'une des revendications 17 à 18, dans lequel l'étape (500) de détermination d'une séquence optimale comporte une comparaison des poids respectifs des opportunités d'acquisition discrètes 30 compatibles cinématiquement entre elles. 20. Procédé selon l'une des revendications 1 à 19, comportant une étape 50 d'agrégation de zones terrestres, antérieure à l'étape 100 de détermination d'opportunités d'acquisition discrètes. 3039728 31 21. Procédé selon l'une des revendications 1 à 20, dans lequel ledit engin spatial est un satellite d'observation fonctionnant selon un mode peigne.