FR 2952744 A1 20110520 FR 1053632 A 20100510 La présente invention appartient au domaine de l'exploitation d'images obtenues par au moins deux barrettes d'un instrument d'observation à défilement embarqué à bord d'un véhicule (avion, drone, hélicoptère, satellite, etc.) survolant un corps céleste, tel que la Terre. Une barrette d'un instrument à défilement comporte une pluralité de cellules sensibles (typiquement de 100 à 100000) généralement agencées en ligne, éventuellement sur plusieurs lignes. La ligne de visée de la barrette est dirigée vers la surface terrestre, et forme avec la ligne de cellules sensibles un plan, dit « plan de visée ». On peut définir la ligne de visée principale de l'instrument soit comme la ligne de visée du pixel central d'une des deux barrettes, soit comme la ligne de visée du pixel central d'une barrette fictive qui serait placée par exemple à égale distance des deux barrettes. Au cours d'une prise de vue, chaque barrette acquiert successivement, par défilement de son plan de visée sur la scène observée, des images lignes représentant différentes portions linéaires de ladite scène. Chaque cellule sensible de la barrette enregistre dans un pixel l'intensité lumineuse d'une zone de la portion linéaire visée, dans la longueur d'onde considérée. Le défilement du plan de visée de chaque barrette s'effectue par exemple par déplacement du véhicule. C'est le cas généralement pour un instrument embarqué dans un satellite placé sur une orbite basse. En général, l'acquisition des images lignes d'une même barrette est réalisée durant une période de temps fixe (durée d'acquisition), et avec une fréquence d'acquisition sensiblement constante. Pour une application typique dans le domaine spatial, un instrument d'observation est embarqué à bord d'un satellite en orbite à 500 à 800 kilomètres d'altitude, sa ligne de visée étant nominalement dirigée vers le centre de la Terre ; la focale de l'instrument est de 0.5 à 15 mètres, le champ de vue est de 1 à 20 degrés, chaque barrette comporte une ligne de 5000 à 50000 cellules sensibles, la vitesse de défilement est comprise entre 6.5 et 6.8 kilomètres par seconde, la fréquence d'acquisition est de 500 à 5000 Hertz ; la fauchée (longueur de la portion linéaire couverte par une image ligne à la surface terrestre) est de 10 à 100 kilomètres, la dimension d'une zone imagée par un pixel est de 0.5 à 20 mètres à la surface de la Terre. Les barrettes sont généralement séparées de 5 à 50 centimètres dans le plan focal ; l'écart temporel entre la prise d'image d'un point de la surface terrestre par la première barrette et sa prise d'image par la seconde barrette est de 0.05 à 5 secondes. Les images lignes successives acquises par une même barrette forment une matrice de pixels dite « image composée », représentant la scène observée durant la prise de vue. On comprend que, du fait que les images lignes sont acquises à des instants différents et du fait des mouvements de l'instrument au cours de la prise de vue, l'image composée répartie dans un plan suivant une grille régulière ne représente pas de façon exacte la scène telle qu'elle devrait être vue idéalement par un observateur surplombant ladite scène. Par exemple, du fait de ces mouvements, une route parfaitement droite pourra être curviligne dans l'image composée. Dans le cas le plus général, au cours d'une prise de vue, deux barrettes différentes fourniront deux images composées différentes dont les pixels acquis à un instant donné représentent des zones différentes de la scène observée. En général, la disposition des barrettes, par exemple parallèles dans le plan focal de l'instrument d'observation, est telle que ces images composées possèdent un recouvrement spatial important lorsque le véhicule défile au-dessus de la scène à observer. Les différences entre les deux images composées viennent du fait que les lignes de visée associées aux deux barrettes sont différentes. De ce fait, les images composées ne se recouvrent pas tout à fait au début et à la fin de la prise de vue. Mais les différences entre les deux images composées proviennent surtout des mouvements de l'instrument eux-mêmes, car ceux-ci vont affecter à tout instant des pixels différents dans chaque image composée, du fait précisément de l'écart entre les lignes de visée des deux barrettes. On connaît, du brevet FR 2 899 344, un procédé de restitution de la ligne de visée dans lequel on détermine, par corrélation, des zones caractéristiques homologues entre les images composées d'une même scène acquises par deux barrettes différentes, et on restitue les variations angulaires de la ligne de visée en fonction des distances géométriques entre les positions (c'est-à-dire les indices des pixels) des zones caractéristiques dans chacune desdites images composées. Cependant, ce procédé possède de nombreuses limitations. Par exemple, la capacité de pouvoir déterminer ou non des zones caractéristiques homologues dans les deux images composées dépend fortement de la nature de la scène observée. Par exemple, pour des scènes faiblement texturées, telles que des scènes de mer/océan et/ou montagne enneigée, il sera en pratique impossible de déterminer des zones caractéristiques homologues en nombre suffisant pour avoir de bonnes performances sur la restitution de la ligne de visée. De plus, le procédé du brevet FR 2 899 344 est sensible à des effets de parallaxe et/ou à la présence de nuages, qui peuvent distordre ou empêcher la corrélation des zones caractéristiques homologues. Dans le procédé du brevet FR 2 899 344, les paramètres estimés par corrélation entre zones caractéristiques homologues sont les variations angulaires relatives d'un instant d'échantillonnage à l'autre, ce qui nécessite une étape supplémentaire de filtrage et d'intégration pour restituer les variations angulaires absolues, étape qui peut s'avérer pénalisante en termes de complexité de calcul. De plus le procédé du brevet FR 2 899 344, qui s'applique très bien à une paire d'images composées, n'est pas facilement généralisable au traitement simultané de plus de deux images composées possédant des décalages temporels différents, voire des fréquences d'acquisition différentes. La présente invention a pour but d'apporter des solutions à toutes ces limitations, en proposant un procédé d'estimation du mouvement de l'instrument d'observation adapté à l'utilisation simultanée des images lignes acquises par toutes les barrettes de l'instrument d'observation, même lorsque lesdites barrettes sont en nombre égal ou supérieur à trois. La présente invention a également pour but de proposer un procédé d'estimation qui soit robuste même lorsque la scène observée comporte une ou des zones faiblement texturées. La présente invention a également pour but de proposer un procédé d'estimation adapté à tenir compte, outre des images composées fournies par les différentes barrettes, d'informations fournies notamment par les senseurs d'attitude embarqués dans le véhicule et/ou de tenir compte d'une connaissance a priori de certaines caractéristiques du mouvement de l'instrument d'observation. La présente invention concerne un procédé d'estimation du mouvement d'un instrument d'observation embarqué dans un véhicule survolant un corps céleste, l'instrument d'observation comportant au moins deux barrettes dont les lignes de visée sont différentes, chaque barrette réalisant l'acquisition successive, par défilement de sa ligne de visée sur le corps céleste, d'images lignes formant une image composée, des parties des images composées représentant sensiblement une même scène du corps céleste. Le mouvement estimé est déterminé par optimisation d'une fonction de ressemblance image associant, à une hypothèse de mouvement donnée, une ressemblance entre des valeurs de pixels mis en correspondance d'au moins une paire d'images composées, des pixels mis en correspondance étant des pixels représentant théoriquement le même point de la scène observée compte tenu de l'hypothèse de mouvement considérée. Le procédé comporte, pour chaque hypothèse de mouvement considérée au cours de l'optimisation de la fonction de ressemblance image, les étapes : - de détermination, pour l'au moins une paire d'images composées, d'au moins une transformation spatio-temporelle en fonction de l'hypothèse de mouvement considérée et d'un modèle géométrique de l'instrument d'observation, - de mise en correspondance des pixels de l'au moins une paire d'images composées au moyen de l'au moins une transformation spatio-temporelle. De préférence, la fonction de ressemblance image est optimisée itérativement, par itération des étapes suivantes jusqu'à ce qu'un critère d'arrêt soit vérifié : - détermination de l'au moins une transformation spatio-temporelle de l'au moins une paire d'images composées d'après l'hypothèse de mouvement de l'itération en cours, - mise en correspondance des pixels de l'au moins une paire d'images composées, - calcul, pour la fonction de ressemblance image, de données images en fonction desdits pixels mis en correspondance, - mise à jour de l'hypothèse de mouvement pour l'itération suivante en fonction des données images. De préférence, les données image comprennent au moins un gradient de la fonction de ressemblance image, représentatif de la variation, au voisinage de l'hypothèse de mouvement de l'itération en cours, de la ressemblance entre les valeurs des pixels mis en correspondance de l'au moins une paire d'images composées. De préférence, la fonction de ressemblance image est une fonction qui augmente lorsqu'une distance entre les valeurs de deux pixels mis en correspondance diminue, et on détermine le mouvement estimé par maximisation de la fonction de ressemblance image. De préférence, l'hypothèse de mouvement est décrite par l'ensemble des angles d'attitude de l'instrument d'observation à des instants échantillonnés de l'acquisition, la fonction de ressemblance image est calculée en fonction de ces angles, et le mouvement estimé est déterminé comme étant l'ensemble des angles d'attitude optimisant la fonction de ressemblance image. De préférence, le mouvement estimé est déterminé par optimisation d'une fonction de ressemblance externe associant, à une hypothèse de mouvement donnée, une ressemblance entre ladite hypothèse de mouvement et des mesures, représentatives de mouvements de l'instrument, effectuées par un ou des capteurs différents des barrettes de l'instrument d'observation. De préférence, la fonction de ressemblance externe est optimisée itérativement, par itération des étapes suivantes jusqu'à ce qu'un critère d'arrêt soit vérifié : - calcul, pour la fonction de ressemblance externe, de données externes en fonction des mesures du ou des capteurs et de l'hypothèse de mouvement de l'itération en cours, - mise à jour de l'hypothèse de mouvement pour l'itération suivante en fonction des données externes. De préférence, les données externes comprennent au moins un gradient de la fonction de ressemblance externe, représentatif de la variation, au voisinage de cette hypothèse de mouvement, de la ressemblance entre l'hypothèse de mouvement de l'itération en cours et les mesures du ou des capteurs. De préférence, la fonction de ressemblance image et la fonction de ressemblance externe sont optimisées conjointement, par optimisation d'une fonction globale cumulant ladite fonction de ressemblance image et ladite fonction de ressemblance externe. De préférence, une première estimation du mouvement est déterminée comme étant l'hypothèse de mouvement optimisant la fonction de ressemblance image, puis, partant de cette première estimation du mouvement, une seconde estimation du mouvement est déterminée comme étant l'hypothèse de mouvement optimisant la fonction de ressemblance externe. De préférence, au moins un capteur est un gyroscope ou un capteur stellaire ou un capteur angulaire d'un miroir orientable de l'instrument d'observation. De préférence, on détermine le mouvement estimé en fonction d'un modèle temporel prédéfini du mouvement de l'instrument d'observation. De préférence, le mouvement estimé est déterminé par optimisation d'une fonction de ressemblance modèle associant, à une hypothèse de mouvement donnée, une ressemblance entre ladite hypothèse de mouvement et le modèle temporel de mouvement. De préférence, le procédé comporte une étape de mise à jour du modèle temporel de mouvement en fonction d'une hypothèse de mouvement. De préférence, la fonction de ressemblance modèle est optimisée itérativement, par itération des étapes suivantes jusqu'à ce qu'un critère d'arrêt soit vérifié : - calcul, pour la fonction de ressemblance de modèle, de données modèle en fonction du modèle temporel de mouvement et de l'hypothèse de mouvement de l'itération en cours, - mise à jour de l'hypothèse de mouvement pour l'itération suivante en fonction des données modèle. De préférence, les données modèle comprennent au moins un gradient de la fonction de ressemblance de modèle, représentatif de la variation, au voisinage de cette hypothèse de mouvement, de la ressemblance entre l'hypothèse temporelle de l'itération en cours et le modèle temporel de mouvement. De préférence, la fonction de ressemblance image et la fonction de ressemblance modèle sont optimisées conjointement, par optimisation d'une fonction globale cumulant ladite fonction de ressemblance image et ladite fonction de ressemblance modèle. De préférence, le modèle temporel de mouvement est un modèle de mouvement paramétrique, de préférence un modèle dont les fonctions de base sont du type polynômes, ou ondelettes, ou sinusoïdes, ou splines. De préférence, le modèle temporel de mouvement est un modèle de mouvement probabiliste, de préférence un modèle de type autorégressif. De préférence, on réduit, dans la fonction de ressemblance image, les contributions de pixels mis en correspondance dont la ressemblance évolue, au cours de l'optimisation, de façon anormale par rapport à la ressemblance des autres pixels mis en correspondance. De préférence, la réduction des contributions de pixels mis en correspondance est effectuée au moyen d'une fonction d'estimation robuste. De préférence, la réduction des contributions de pixels mis en correspondance est effectuée au moyen d'au moins une matrice de pondération et, lorsque la fonction de ressemblance image est optimisée itérativement, le procédé comporte de préférence une étape de mise à jour de l'au moins une matrice de pondération en fonction des contributions des pixels mis en correspondance à ladite fonction de ressemblance image. De préférence, le procédé est itéré selon un schéma d'optimisation 30 multi-résolution. L'invention sera mieux comprise à la lecture de la description suivante, donnée à titre d'exemple nullement limitatif, et faite en se référant aux figures qui représentent : - Figure 1 : une représentation très schématique d'un instrument d'observation à défilement survolant la Terre, - Figure 2, 3, 4 et 5: des diagrammes illustrant des modes particuliers de mise en oeuvre d'un procédé selon l'invention. La présente invention concerne un procédé d'estimation du mouvement d'un instrument d'observation 10, embarqué dans un véhicule (avion, drone, hélicoptère, satellite, etc.) survolant un corps céleste, au cours de la prise de vue. Dans la suite de la description, on se place de manière non limitative dans le cas où le corps céleste survolé et observé est la Terre. La mise en oeuvre de ce procédé d'estimation du mouvement de l'instrument d'observation 10 est par exemple assurée par un dispositif de calcul (microcontrôleur, ordinateur muni d'au moins un microprocesseur, etc.). Le dispositif de calcul comporte des moyens de mémorisation (disque dur magnétique, mémoire flash, disque optique, etc.) dans lesquels est mémorisé un produit programme d'ordinateur, sous la forme d'un ensemble d'instructions de code de programme à exécuter pour mettre en oeuvre les étapes du procédé selon l'invention, à partir des images lignes acquises par l'instrument d'observation 10. Suivant certains modes de réalisation, le dispositif de calcul comporte également un ou des circuits intégrés, de type ASIC, FPGA, etc., adaptés à mettre en oeuvre tout ou partie des étapes du procédé selon l'invention. Le dispositif de calcul peut être embarqué avec l'instrument d'observation 10 dans le véhicule survolant la Terre dans le cas d'un traitement en temps réel. Dans le cas d'un traitement en temps différé, le dispositif de calcul est indifféremment embarqué avec l'instrument d'observation 10 dans le véhicule, ou déporté, c'est-à-dire localisé dans une station au sol ou dans un second véhicule survolant la Terre. A) Définitions et notations Des définitions et notations utilisées dans la suite de la description sont données ci-après. On comprend que le choix d'une convention particulière n'est pas limitatif de l'invention qui pourrait être décrite de manière équivalente en adoptant d'autres conventions sans que l'invention en soit modifiée. A.1) Instrument d'observation, imaqes lignes, imaqes composées On considère un instrument d'observation 10 comportant Nb barrettes b;, 1 i Nb. On suppose pour simplifier que chaque barrette comporte le même nombre K de cellules sensibles. On se place dans le cas non limitatif où les barrettes b; sont agencées de façon sensiblement parallèle dans le plan focal de l'instrument d'observation 10. Le défilement de la ligne de visée de chaque barrette b; permet d'acquérir successivement des images lignes correspondant à des portions linéaires différentes d'une scène S observée, c'est-à-dire des portions disjointes ou qui se recouvrent au plus partiellement. Dans la suite de la description, on se place dans le cas où le défilement de la ligne de visée d'une barrette s'effectue par déplacement du véhicule par rapport à la Terre (un tel instrument est connu sous le nom de « push broom » dans la littérature anglo-saxonne). Cet exemple n'est cependant pas limitatif de l'invention. Par exemple, le défilement de la ligne de visée d'une barrette peut être effectué par déplacement de ladite ligne de visée par rapport au véhicule, au moyen d'un mécanisme de balayage portant soit la barrette, soit un miroir de visée se trouvant sur le chemin optique de la barrette. Cela peut être le cas pour un satellite d'observation en orbite géostationnaire, stabilisé selon trois axes. A des fins de simplification des équations, et de manière nullement limitative de l'invention, on considère que la période d'acquisition d'images lignes, désignée par AT, est constante au cours du temps, et la même pour toutes les barrettes b;. On notera toutefois que dans le cas d'application le plus courant, notamment pour un système d'observation par satellite, on utilise une barrette fonctionnant dans une bande panchromatique ainsi que d'autres barrettes fonctionnant dans des bandes spectrales dédiées (par exemple vert, rouge, proche infrarouge), la fréquence d'acquisition de la barrette panchromatique étant généralement quatre fois plus élevée que celle des autres barrettes. Au cours de la prise de vue, N images lignes sont acquises par chacune des barrettes b; aux instants (n-1)v.T, 1 n N. Chaque barrette b; acquiert ainsi une image composée B; constituée par l'ensemble {B;(k,n), 1 k K, 1 n N}, où B;(k,n) est la valeur du pixel d'indice (k,n). L'indice k correspond à une dimension spatiale : la position de la cellule sensible le long de la barrette ; l'indice n représente une dimension temporelle : l'instant de la prise de vue. D'où l'appellation « grille spatio-temporelle » que l'on associera à l'image composée. Pour simplifier les notations, on indexera les instants d'acquisition de sorte que l'instant d'acquisition d'indice n est l'instant (n-1)..T. Au cours d'une prise de vue, deux barrettes différentes fourniront deux images composées différentes dont les pixels acquis à un même instant donné représentent des zones différentes de la scène observée. La disposition des barrettes dans le plan focal de l'instrument d'observation est telle que ces images composées possèdent un recouvrement spatial important lorsque le véhicule défile au-dessus de la scène à observer. Le procédé ne traitera que la partie des images composées qui correspond à ce recouvrement spatial. A.2) Repère instrument La figure 1 représente très schématiquement en perspective le survol d'une scène par un instrument d'observation 10 comportant au moins deux barrettes bl, b2 considérées parallèles dans le plan focal de l'instrument. On associe à l'instrument d'observation 10 un repère instrument RINSTR dont l'origine M représente la position de l'instrument d'observation 10 dans le véhicule. Le repère instrument RINSTR comporte trois axes XINSTR, YINSTR et ZINSTR, orthogonaux entre eux. L'axe XINSTR est parallèle aux barrettes b1 et b2, l'axe ZINSTR, correspondant à la direction de visée de l'instrument 10, est orienté vers la Terre, et l'axe YINSTR est orthogonal aux axes XINSTR et ZINSTR, et sensiblement dans la direction de déplacement de l'instrument d'observation 10, par déplacement du véhicule. A.3) Mouvement de l'instrument d'observation Dans le cas d'un défilement du plan de visée d'une barrette obtenu par le déplacement du véhicule, les mouvements de translation de l'instrument d'observation 10 proviennent du survol du véhicule au dessus de la scène. Ces mouvements de translation s'effectuent en général à des fréquences relativement basses et sont peu observables par le procédé, en particulier lorsque l'altitude du véhicule est haute. En outre, ils sont difficiles à distinguer des mouvements de rotation de l'instrument d'observation 10. C'est pourquoi on supposera pour simplifier que les mouvements de translation de l'instrument sont connus avec une précision suffisante, et que l'on cherche à restituer avant tout les mouvements de rotation de l'instrument d'observation 10. Il est cependant à noter, comme on le verra par la suite, que le procédé peut être facilement généralisé à l'estimation du mouvement complet de l'instrument, aussi bien en translation qu'en rotation, et peut s'étendre à l'estimation des mouvements de toute partie mobile de l'instrument, comme un miroir de visée. Les mouvements de rotation de l'instrument d'observation 10 (supposé fixe par rapport au véhicule) proviennent notamment des mouvements de rotation du véhicule porteur, dits « mouvements d'ensemble » de l'instrument d'observation 10. Ces mouvements sont en général estimés grâce à des capteurs dédiés embarqués à bord du véhicule, comme par exemple des gyroscopes, et des senseurs stellaires si le véhicule est un satellite. Cependant, cette estimation est entachée d'erreurs dues en particulier aux erreurs de mesures. D'autres mouvements de l'instrument, et plus particulièrement de sa ligne de visée, dits « mouvements parasites », proviennent de déplacements relatifs en translation et/ou en rotation entre différents éléments optiques de l'instrument (plan focal, miroirs, etc.), dus par exemple à des vibrations, ces déplacements affectant la ligne de visée de l'instrument. Les mouvements parasites sont beaucoup moins facilement mesurables, car il faudrait pour cela un grand nombre d'accéléromètres et de gyromètres répartis sur toute la structure de l'instrument d'observation. Les mouvements d'ensemble et les mouvements parasites se distinguent en général par leur amplitude et leur spectre fréquentiel. Les mouvements d'ensemble sont de l'ordre du centième de degré à quelques degrés, et sont relativement basses fréquences (de zéro à quelques Hertz). Les mouvements parasites sont un ordre de grandeur plus petits que les mouvements d'ensemble, et sont plus hautes fréquences (de quelques dizaines à quelques centaines de Hertz). Les mouvements d'ensemble de l'instrument d'observation 10 se traduisent par des rotations du repère instrument RINSTR par rapport à un repère de référence terrestre RT (par exemple le repère géographique local au point de référence de la scène observée, non représenté sur les figures) qui traduisent l'évolution temporelle de l'attitude de l'instrument par rapport audit repère de référence RT. Les mouvements parasites se traduisent par des petites rotations exprimées par rapport au repère instrument RINSTR, celles-ci se décomposant en rotations autour de l'axe YINSTR, dites « rotations en roulis », autour de l'axe XINSTR, dites « rotations en tangage », et autour de l'axe ZINSTR, dites « rotations en lacet ». Il est à noter que des erreurs d'estimation du mouvement d'ensemble de l'instrument d'observation 10 peuvent se formuler aux petits angles de la même façon que les mouvements parasites, en termes d'erreurs de roulis, tangage, lacet. Le mouvement de l'instrument d'observation 10 peut donc être décrit comme les valeurs prises au cours du temps d'angles d'attitude (angles d'Euler par exemple) par rapport à un repère de référence. D'autres formulations équivalentes à base d'autres angles, de quaternions ou de matrices de rotation peuvent être utilisées sans que l'invention en soit modifiée dans son principe. On désigne par OR(n), 1 n N, le mouvement réel de l'instrument d'observation 10 à l'instant d'acquisition d'indice n, par exemple exprimé dans le repère de référence terrestre RT. Comme nous l'avons vu, ce mouvement est la somme des mouvements d'ensemble et des mouvements parasites. Il est à noter que, dans les cas où le balayage de la ligne de visée s'effectue par rotation d'un miroir de visée présent sur le chemin optique de l'instrument, le ou les angles de rotation du miroir sont inclus dans le vecteur OR(n) à tout instant n,1nN. On désigne par OR l'ensemble de tous les OR(n), 1 n N, représentant le mouvement réel de l'instrument au cours de la prise de vue. On désigne par OEST(n), 1 n N, le mouvement estimé de l'instrument d'observation 10 à l'instant d'acquisition d'indice n exprimé dans le repère de référence RT. On désigne par C'EST l'ensemble de tous les OEST(n), 1 n N, représentant le mouvement estimé de l'instrument d'observation 10 au cours de la prise de vue. On désigne par O = {8(n), 1 n N} une hypothèse de mouvement de l'instrument 10 au cours de la prise de vue aux instants d'acquisition d'indice n, exprimée dans le repère de référence RT. A.4) Modèle qéométrique, transformation spatio-temporelle, et mise en correspondance des pixels des imaqes composées On suppose connu un modèle géométrique R;mod(k) de l'instrument qui associe à la cellule sensible d'indice k de la barrette b;, sa direction de visée exprimée dans le repère instrument RINSTR. Ce modèle géométrique peut également modéliser les distorsions optiques introduites par l'instrument d'observation 10, qui sont par exemple préalablement estimées par calibration. On considère un repère cartographique Rc = (O,Xc,Yc,Zc) (non représenté sur les figures), par exemple tel que (Xc,Yc) est le plan tangent à la surface terrestre en un point de référence O de la scène observée, dans lequel pourront se projeter les images composées acquises par les barrettes en vue de la reconstitution de l'image réelle de la scène survolée. Le plan (O,Xc,Yc) sera désigné par « référentiel image » dans la suite de la description. Pour simplifier les équations, on se place dans le cas non limitatif où les repères RT et Rc sont confondus. On note cependant que l'on pourrait considérer d'autres référentiels image, éventuellement non planaires, comme par exemple une projection sur l'ellipsoïde terrestre, sur un modèle connu de relief, ou encore l'une des images composées acquises par l'une des barrettes qui sont déformées par le mouvement de l'instrument d'observation 10. A partir d'une hypothèse de mouvement donnée O de l'instrument, chaque pixel d'indice q = (k,n) de l'image composée, acquis par la cellule sensible d'indice k de la barrette b; à l'instant n, peut être projeté dans le référentiel image en calculant l'intersection de la droite de visée du pixel d'indice q = (k,n) avec le plan (O,Xc,Yc). Pour cela, on calcule la direction de visée d;(k,n), d'une droite passant par la position M de l'instrument 10 et par la cellule sensible d'indice k à l'instant n, dans le repère cartographique Rc à partir de l'hypothèse de mouvement 8(n) de l'instrument 10 à l'instant n et du modèle géométrique de l'instrument : di(k,n) = Re(n) .Rimod(k) où Re(n) est une rotation qui transforme le repère instrument RINSTR dans le repère cartographique Rc. Le projeté p = (x,y) du pixel d'indice q = (k,n) dans le référentiel image (O,Xc,Yc) correspond à l'intersection de la droite de direction d;(k,n) avec le plan (O,Xc,Yc). On obtient ainsi une transformation V; qui, en fonction d'une hypothèse de mouvement O, associe à un pixel d'indice q = (k,n) de l'image composée, le point p = (x,y) du référentiel image, qui est le point effectivement représenté par le pixel B;(k,n) p = Vi(q,O) (1) Dans cette relation (1), les variables sont p et q, l'hypothèse de mouvement O étant considérée comme un ensemble de paramètres dont dépend la transformation V;, 1 i Nb. Le calcul de la projection peut être fait de manière analytique si le modèle géométrique de l'instrument d'observation 10 n'est pas trop compliqué. On peut également faire des tirages de rayons dans le cas d'un modèle géométrique d'instrument plus complexe. Un modèle de relief de la scène peut également être utilisé pour prendre en compte les décalages qu'il introduit. On désigne par p;k'n = (Xik,n ylk,n) la position du pixel d'indice (k,n) ainsi projeté dans le référentiel image (O,Xc,Yc), exprimée en mètres ou toute autre unité. On désigne par « grille spatio-temporelle » de l'image composée, l'ensemble des points p;k'n projetés des pixels de l'image composée dans le référentiel image. On désigne par « image rectifiée par le mouvement » l'image composée qui associe à la position p;k,n la valeur du pixel B;(k,n). Pour les besoins du procédé, il est nécessaire de calculer les transformations géométriques W;, inverses de V;, dites « transformations spatio-temporelles ». Etant donnée une hypothèse de mouvement O donnée, les transformations spatio-temporelles associent, à un point p = (x,y) du référentiel image dans la scène imagée par une barrette b; au cours de la prise de vue, la position q du pixel dans l'image composée B; acquise par la barrette b; dont le projeté dans le référentiel image est le point p = (x,y) compte tenu de l'hypothèse de mouvement O. Pour cela, il faut considérer que l'hypothèse de mouvement O = {0(n), 1 n N} est interpolée à des valeurs non entières du temps : O = {0(t), t réel dans l'intervalle [0,(n-1).AT]}, et que l'image composée B; rectifiée par le mouvement est également interpolée suivant tout procédé d'interpolation connu de l'homme de l'art, tel qu'une interpolation cubique, de Shannon tronquée, ou utilisant des splines, de sorte qu'on peut définir la valeur de pixels B;(a,t) pour des positions continues a d'une cellule sensible virtuelle le long de la barrette b;, et pour un instant d'acquisition t continu. Considérant le point p = (x,y) du référentiel image, on peut positionner ce point par rapport à la grille spatio-temporelle de l'image B; rectifiée par le mouvement, ce qui donne une première approximation entière (k',n') de la position du pixel dans l'image composée. On obtient des valeurs réelles (a,t) raffinant cette approximation par exemple en travaillant sur de petites variations et en linéarisant le problème. De façon plus générale, toute méthode connue de l'homme de l'art peut être utilisée pour inverser la transformation V. On appellera W; la transformation inverse de V; qui à tout point p = (x,y) du référentiel image associe la position q = (a,t) du pixel qui image ce point dans l'image composée B; compte tenu de l'hypothèse de mouvement O : q = Wi(p,O) (2) les variables p et q pouvant être égales à un nombre non entier par interpolation. Les transformations spatio-temporelles W;, 1 i Nb, permettent de mettre en correspondance les pixels des images composées d'après une hypothèse de mouvement O, c'est-à-dire de déterminer, pour tout point p considéré du référentiel image, les valeurs des pixels B;(W;(p,O)), 1 i Nb, qui représentent théoriquement ce point p dans les images composées. Par « représentent théoriquement ce point », on entend que les pixels mis en correspondance ne représentent effectivement le même point que lorsque l'hypothèse de mouvement O considérée est égale au mouvement réel OR. Lorsque l'hypothèse de mouvement O considérée sera très différente du mouvement réel OR, on comprend donc que les pixels mis en correspondance, bien que représentant théoriquement le même point, représenteront en pratique des points différents de la scène observée, et que les valeurs de ces pixels seront peu corrélées. Dans le cas où l'hypothèse de mouvement correspond à la réalité (O = 15 OR), on devrait avoir, aux erreurs d'interpolation près, pour tout point p dans le référentiel image et quels que soient i et j, 1 i,j Nb : Bi(Wi(p,OR )) = B;(W;(p,OR )) (3) car ces deux pixels mis en correspondance des images composées interpolées correspondent au même point p du référentiel image. 20 On comprend donc que l'on peut évaluer la cohérence d'une hypothèse de mouvement O donnée par comparaison des pixels mis en correspondance B;(W;(p,O)), 1 i Nb. Notons que la mise en correspondance des images composées peut s'effectuer à partir de la transformation spatio-temporelle : 25 Wii=WioWi-1=WioV; (4) qui fait passer directement d'une image composée B; à une autre image composée Bi. Pour calculer Wj;, on considère une paire de barrettes b;, bj, et une hypothèse de mouvement O. On considère un pixel q = (k,n) de l'image composée acquise par la barrette b;. On projette ce pixel dans le référentiel 30 image suivant l'hypothèse de mouvement O, ce qui donne le point pi = V;(q,O), puis on calcule la position qj = (a,t) de l'image composée Bj qui se projette en Pi: qi = Wi(Vi(q,O)) = Wii(q,O) (5) Dans tous les cas, on parlera d'une transformation spatio-temporelle mettant en correspondance les pixels des images composées B;, Bi obtenues par les barrettes b;, bi. Notons que cette transformation spatio-temporelle dépend également du référentiel image choisi, du modèle géométrique de l'instrument, et éventuellement d'un modèle de terrain que l'on introduit pour effectuer la mise en correspondance. B) Description du procédé d'estimation B.1) Principe qénéral du procédé d'estimation Selon l'invention, le mouvement estimé est déterminé par optimisation d'une fonction, dite « fonction F1 de ressemblance image » associant, à une hypothèse de mouvement O donnée, une ressemblance entre les valeurs de pixels mis en correspondance d'au moins une paire d'images composées. Tel qu'on l'a vu précédemment, des pixels mis en correspondance sont des pixels représentant théoriquement le même point de la scène observée compte tenu de l'hypothèse de mouvement considérée. Pour chaque hypothèse de mouvement O considérée au cours de l'optimisation de la fonction F1 de ressemblance image, le procédé comporte 20 les étapes : - 200 de détermination, pour l'au moins une paire d'images composées, d'au moins une transformation spatio-temporelle en fonction de l'hypothèse de mouvement O considérée et d'un modèle géométrique de l'instrument d'observation 10, 25 - 201 de mise en correspondance des pixels de l'au moins une paire d'images composées au moyen de l'au moins une transformation spatio-temporelle. Les opérations effectuées au cours des étapes 200 et 201 ont déjà été décrites au cours de la section A.4) précédente. 30 La fonction F1 de ressemblance image comprend par exemple une composante F12 prenant en compte les images composées B1, B2 acquises par les barrettes b1, b2. Les pixels de ces images composées sont mis en correspondance au moyen des transformations spatio-temporelles W1 et W2, ou de l'unique transformation spatio-temporelle W21 (ou W12). De préférence, on choisit comme composante F12 une fonction avantageusement convexe qui augmente lorsqu'une distance entre les matrices des valeurs des pixels mis en correspondance des deux images composées B1 et B2 diminue. Dans ce cas l'optimisation consiste à déterminer une hypothèse de mouvement qui tend à maximiser cette composante F12. Néanmoins, on pourrait de manière équivalente choisir une fonction concave qui diminue lorsqu'une distance entre les matrices des valeurs des pixels mis en correspondance des deux images composées B1 et B2 diminue. Dans ce cas l'optimisation consiste à déterminer l'hypothèse de mouvement qui tend à minimiser cette composante F12. La distance peut être quadratique, linéaire, sous forme de produit scalaire, avec ou sans pondération dépendant du pixel, etc. On note qu'il s'agit ici d'une distance entre les valeurs de pixels représentant théoriquement le même point, et non, tel que décrit dans le brevet FR 2 899 344, d'une distance géométrique entre les positions (indices) des pixels représentant effectivement le même point (zone caractéristique homologue) dans chacune des images composées. Par exemple, on peut prendre pour fonction de ressemblance image l'une des fonctions suivantes: i F12(0) = exp ù Lai; . hi(B,(W,(p, O)) ù B2(W2(p,O))) (6a) F12(0) = p=(i,1) / (6b) agi; . hi(BI(W,(p,O)) ù B2(W2(p, O))) p=(i,1) où p = (i,j) correspond à un quadrillage de tout ou partie de la zone du référentiel image qui est imagée par les deux barrettes b1 et b2, a;; sont des pondérations positives ou nulles, h1 est une p-norme. Les exemples donnés aux expressions (6a) et (6b) ne sont pas limitatifs. On note que l'optimisation de la composante F12 de l'expression (6a) est une maximisation, tandis que l'optimisation de la composante F12 de l'expression (6b) est une minimisation. Il est à noter que le quadrillage p = (i,j) peut être régulier, ou plus dense dans des zones pré-spécifiées de l'image composée dont les propriétés sont jugées plus favorables que d'autres pour des opérations de mise en correspondance et évaluation de ressemblance. Un autre exemple de fonction de ressemblance image peut être sous forme de produit scalaire : F12(0) _ agi; . h1(B1(W,(p,o)). B2(W2(p,o))) (7) p=(i,1) Tel qu'indiqué précédemment, la mise en correspondance pourrait se traiter directement d'image composée à image composée sans passer par la grille du référentiel image. Dans ce cas, la composante F12 de l'expression (6a) pourrait être réécrite comme suit : i F12(0)=exp û Lakn.h1(B1(p)ûB2(W21(p,o))) (8) p=(k,n) / où cette fois-ci, p = (k,n) correspond à tout ou partie des indices des pixels de l'image composée B1 dans la zone de recouvrement des images composées B1 et B2. Les pondérations a;; peuvent être réglées individuellement avant l'optimisation de la fonction de ressemblance image. En outre, ces pondérations peuvent être choisies toutes égales, par exemple égales à un. Il est à noter que la composante F12 peut se présenter sous d'autres formes. En particulier, des opérateurs peuvent être introduits sur les images composées pour tenir compte du fait qu'elles ont été acquises dans des domaines de longueurs d'ondes différents. De tels opérateurs sont considérés connus de l'homme de l'art, et ne sont pas considérés ici à des fins de simplification des expressions et à titre nullement limitatif de l'invention. Une composante F12 selon l'un quelconque des exemples ci-dessus est particulièrement avantageuse car elle ne nécessite pas une connaissance a priori de la scène observée. En considérant ce type de ressemblance image entre images composées, on exploite une connaissance a priori du système d'acquisition qui est celle que les images composées B;, 1 i Nb, représentent sensiblement la même scène observée et sont théoriquement superposables (une fois mises en correspondance) en l'absence de mouvements parasites de l'instrument d'observation 10 (mises à part les différences liées aux domaines de longueurs d'ondes d'acquisition dans le cas d'un instrument d'observation 10 multi-spectral, les erreurs d'interpolation, les bruits de mesure, les défauts de calibration, etc.). En outre, ce type de composante est insensible aux textures ambiguës telles les textures périodiques ou les textures à forte direction principale (des lignes par exemple). De manière plus générale, dans le cas où on prend en compte Nb barrettes, Nb étant égal ou supérieur à trois, la fonction FI de ressemblance image peut facilement prendre en compte tout ou partie des paires d'images composées disponibles, en cumulant les composantes F;i des différentes paires d'images composées considérées. Les composantes F;i sont avantageusement cumulées par produit dans le cas de composantes du type donné à l'expression (6a) : F,(0) = f ft . Fi; (0) (9a) La prise en compte simultanée des images lignes acquises par toutes les barrettes considérées permettra en général d'améliorer la robustesse du procédé par rapport à l'état de l'art, en présence de zones faiblement texturées, en particulier si les images lignes sont acquises dans des domaines de longueurs d'ondes différents. La figure 2 représente schématiquement un mode préféré de mise en oeuvre du procédé d'estimation, dans lequel l'optimisation de la fonction FI de ressemblance image est itérative, par itération des étapes suivantes jusqu'à ce qu'un critère d'arrêt soit vérifié : - 200 détermination des transformations spatio-temporelles W ou W;i, pour chacune des paires d'images composées considérées, d'après l'hypothèse de mouvement Ok de l'itération k en cours, - 201 mise en correspondance des pixels des paires d'images composées considérées au moyen des transformations spatio- temporelles déterminées à l'étape 200, - 202 calcul, pour la fonction FI de ressemblance image, de données images utiles pour effectuer l'optimisation, - 203 mise à jour de l'hypothèse de mouvement pour l'itération suivante, en fonction des données images. L'hypothèse de mouvement 01 considérée à la première itération correspond par exemple au mouvement idéal de l'instrument d'observation 10 tel qu'il devrait être si les ordres envoyés au système de contrôle d'attitude du véhicule étaient parfaitement respectés, et en l'absence de mouvements parasites au cours de l'acquisition. Suivant un autre exemple non limitatif, l'hypothèse de mouvement 01 est une première estimation du mouvement de l'instrument d'observation 10, par exemple déterminée en fonction de capteurs dédiés (senseurs d'attitude) différents de l'instrument d'observation 10. Les opérations effectuées au cours des étapes 200 et 201 ont déjà été décrites au cours de la section A.4) précédente. Au cours de l'étape 202 de calcul des données images, on détermine par exemple la valeur de la fonction FI de ressemblance image d'après l'hypothèse de mouvement Ok de l'itération k en cours. De préférence, on détermine, au cours de l'étape 202 de calcul des données images, un ou plusieurs gradients de la fonction FI de ressemblance image. Un tel gradient est par exemple déterminé en calculant la valeur de la fonction FI pour les hypothèses de mouvement Ok et Ok + 8Ok, où 8Ok est une variation de l'hypothèse de mouvement de faible norme. Ces gradients donnent une information sur la direction dans laquelle faire évoluer l'hypothèse de mouvement Ok pour optimiser la fonction FI de ressemblance image. Au cours de l'étape 203 de mise à jour, l'hypothèse de mouvement Ok+' est par exemple déterminée en fonction de l'hypothèse de mouvement Ok, à laquelle on ajoute un terme de mise à jour AOk qui est déterminé d'après les données images. A®k est par exemple un gradient calculé au cours de l'étape 202, permettant d'optimiser la fonction FI de ressemblance image. Si le critère d'arrêt est vérifié, le mouvement estimé C'EST correspond soit à l'hypothèse de mouvement Ok soit à l'hypothèse de mouvement Ok+'. Par exemple le critère d'arrêt est vérifié lorsque la norme du terme de mise à jour DOk est inférieure à un seuil prédéfini. On comprend que cet exemple n'est pas limitatif de l'invention, et que tout autre critère adapté, connu de l'homme de l'art, peut être mis en oeuvre dans le procédé d'estimation. B.2) Description de modes particuliers de mise en oeuvre du procédé Prise en compte de mesures externes Dans un mode particulier de mise en oeuvre du procédé d'estimation, compatible avec l'un quelconque des modes précédemment décrits, le mouvement estimé est également déterminé par optimisation d'une fonction F2, dite « de ressemblance externe ». La fonction F2 de ressemblance externe associe, à une hypothèse de mouvement O donnée, une ressemblance entre ladite hypothèse de mouvement et des mesures représentatives des mouvements de l'instrument d'observation 10, lesdites mesures étant effectuées par un ou des capteurs différents des barrettes dudit instrument d'observation. Ces capteurs sont désignés par « capteurs externes » dans la suite de la description. Du fait de la détermination du mouvement estimé en prenant en compte des mesures faites par l'instrument d'observation 10 (dans la fonction FI de ressemblance image) et des mesures faites par un ou des capteurs externes (fonction F2 de ressemblance externe), le procédé d'estimation devient plus robuste à la nature de la scène observée, dans la mesure où une la présence d'une zone faiblement texturée (mer/océan, nuage, montagne enneigée, etc.) ne serait défavorable, en termes de robustesse, qu'à l'optimisation par rapport à la fonction FI de ressemblance image. De préférence, au moins un capteur externe est un gyroscope ou un capteur stellaire (parfois appelé « star tracker » dans la littérature anglo-saxonne), embarqué dans le même véhicule que l'instrument d'observation 10. En complément ou en alternative, lorsque le défilement de la ligne de visée est obtenu au moyen de miroirs de visée, au moins un capteur externe est un capteur angulaire mesurant le débattement de ces miroirs à balayage, les mouvements de ces miroirs étant inclus dans le vecteur OR(n), pour 1 n N et faisant partie desdits mouvements de l'instrument à estimer. Suivant encore d'autres exemples, il est possible d'utiliser des données issues des images composées acquises par les barrettes, comme par exemple la position connue dans le repère cartographique Rc d'éléments particuliers et remarquables dans l'image composée, appelés communément amers ou points d'appui dans l'image composée. Dans la suite de la description, on se place dans le cas non limitatif de mesures du mouvement obtenues par un gyroscope et un capteur stellaire. Les mesures effectuées par un capteur stellaire sont des mesures absolues qui peuvent par exemple être exprimées sous la forme suivante : ysT(n) = 6R(n) + £sT(n), 1 L'expression (10) peut également s'exprimer sous la forme : MSTYST = MsT'OR + MsT'EST (11) où MsT est un opérateur de blanchissement du bruit de mesure EsT, de sorte que MsT.EsT possède sensiblement les caractéristiques d'un bruit blanc. Les mesures effectuées par un gyroscope sont des mesures différentielles qui peuvent par exemple être exprimées sous la forme suivante : yG(n) = OR(n) - OR(n-1 ) + EG(n), 2 n N (12) où EG désigne un bruit blanc de mesure. On désigne par YG le vecteur comprenant l'ensemble {yG(n), 2 n N}, par EG le vecteur comprenant l'ensemble {EG(n), 1 n N}. L'expression (12) peut également s'exprimer sous la forme : YG = MG'OR + EG(n) (13) où MG est un opérateur de convolution associé au filtre ayant pour réponse impulsionnelle [1, -1]. La fonction F2 de ressemblance externe, dans le cas d'un capteur stellaire et d'un gyroscope, s'exprime de préférence sous la forme suivante : F2(0)=exp(ùh2(Msr (Ysr -0))ùh3(YG ù MG .0)) (14) où h2 et h3 sont des p-normes. D'autres formes de la fonction F2 de ressemblance externe sont possibles, et l'on peut introduire par exemple des pondérations fréquentielles rendant compte de la bande passante des différents capteurs externes. De préférence, la fonction FI de ressemblance image et la fonction F2 de ressemblance externe sont optimisées conjointement. Par exemple, on optimise une fonction globale F cumulant ladite fonction de ressemblance image et ladite fonction de ressemblance externe. Rien n'exclut toutefois de mettre en oeuvre d'autres schémas d'optimisation qu'une optimisation conjointe des fonctions FI de ressemblance image et F2 de ressemblance externe. Suivant un autre exemple, on détermine l'hypothèse de mouvement O qui optimise la fonction FI, puis à partir de cette hypothèse de mouvement, on optimise la fonction F2. Cette optimisation alternée des fonctions FI et F2 peut en outre être itérée plusieurs fois. Suivant encore un autre exemple, au cours d'un processus itératif d'optimisation, par exemple par la méthode du gradient, on considère une fonction FI ou F2 différente d'une itération à l'autre, c'est à dire en alternant d'une itération à l'autre la fonction FI ou F2 dont on calcule le gradient. Une fonction globale F, pour une optimisation conjointe, prend avantageusement la forme d'une somme, éventuellement pondérée, des fonctions FI et F2. Suivant d'autres exemples, la fonction globale F pourra prendre la forme d'un produit des fonctions FI et F2. On note que, en prenant l'opposé d'un logarithme d'un produit de fonctions FI et F2, on obtient de manière équivalente pour l'invention une fonction globale sous la forme d'une somme des logarithmes des fonctions FI et F2. 24 Des exemples de fonctions globales F cumulant les fonctions FI et F2 peuvent s'exprimer comme suit : ( F(0)=exp -L Lai; h1(B1(W1(p,0))ùB2(W2(p,0))) x i Au cours de l'étape 204 de calcul des données externes, on détermine par exemple la valeur de la fonction F2 de ressemblance externe d'après l'hypothèse de mouvement C'k de l'itération k en cours et d'après les mesures effectuées par les capteurs externes. De préférence, on détermine, au cours de l'étape 204 de calcul des données externes, un ou plusieurs gradients de la fonction F2 de ressemblance externe, afin d'avoir une information sur la direction dans laquelle faire évoluer l'hypothèse de mouvement C'k de l'itération k en cours pour optimiser la fonction F2 de ressemblance externe. (15a) (15b) Au cours de l'étape 203 de mise à jour du mouvement estimé, on prend en compte les données images et les données externes, soit cumulativement (en tenant compte de l'ensemble des données images et externes à chaque itération k), soit en alternance (en tenant compte alternativement des données image à l'étape k du procédé, et des données externe à l'étape k+1). L'hypothèse de mouvement Ok+' pour l'itération (k+1) est par exemple déterminée en fonction de l'hypothèse de mouvement Ok, à laquelle on ajoute un terme de mise à jour DOk qui est déterminé d'après les données images et les données externes. De préférence, dans un mode cumulatif, le terme de mise à jour DOk est une somme pondérée des gradients calculés au cours de étapes 202 et 204, lesdits gradients étant pondérés par exemple par les valeurs des fonctions respectivement FI de ressemblance image et F2 de ressemblance externe. Dans un mode alterné, le terme de mise à jour sera alternativement le terme de gradient correspondant à la fonction FI de ressemblance image, et le terme de gradient correspondant à la fonction F2 de ressemblance externe. Prise en compte d'un modèle de mouvement Dans un mode particulier de mise en oeuvre, compatible avec l'un quelconque des modes de mise en oeuvre précédemment décrits, le mouvement de l'instrument d'observation 10 est supposé suivre un modèle, dit « modèle de mouvement ». Le mouvement de l'instrument d'observation 10 est continu au cours 25 du temps, de type généralement vibratoire, de sorte que les mouvements 8(n) sont corrélés d'un instant d'acquisition à l'autre. L'utilisation d'un modèle de mouvement permet de contraindre la recherche du mouvement estimé, de sorte que le nombre d'inconnues à déterminer est généralement inférieur au nombre total des paramètres 8(n), 1 30 n N, à estimer en l'absence de modèle de mouvement, égal à 3N (roulis, tangage, lacet). Un autre intérêt consiste à permettre ainsi de contraindre le contenu fréquentiel du mouvement estimé, notamment si l'on cherche à estimer plus particulièrement les mouvements parasites, dont on a vu que leur spectre était généralement principalement compris entre quelques dizaines de Hertz et quelques centaines de Hertz. Suivant un premier exemple, le modèle de mouvement est paramétrique, c'est-à-dire qu'il est décrit par un ensemble limité de paramètres. De préférence, le modèle de mouvement paramétrique utilise des combinaisons linéaires de fonctions de base, pouvant être de type polynômes, ondelettes, sinusoïdes, splines, etc. Un tel modèle de mouvement s'exprime par exemple, pour les angles de roulis cp(n) sous la forme suivante : M-1 9(n) = LÀm . pm(n) • wm(n), 1 n N (17) m=0 expression dans laquelle les paramètres du modèle de mouvement paramétrique sont les suivants : - M est l'ordre du modèle de mouvement paramétrique, - wm(n), 0 m M-1, sont les fonctions de base, - Mn), 0 m M-1, sont des coefficients adaptant localement le modèle de mouvement au cours des instants d'acquisition, - Àm sont des coefficients scalaires. Les modèles de mouvement pour les angles de tangage a(n) et de lacet ip(n) suivent le même format que celui donné à l'expression (17), les paramètres étant éventuellement différents. Suivant un autre exemple, le modèle de mouvement est probabiliste. De préférence, il s'agit d'un modèle de mouvement autorégressif, s'exprimant par exemple, sous la forme suivante : L 8(n) = L al . 8(n -1) + e(n), 1 n N (18) 1=1 expression dans laquelle les paramètres du modèle de mouvement autorégressif sont les suivants : - L est l'ordre du modèle de mouvement autorégressif, - a,, 1 I L, sont les coefficients du modèle de mouvement autorégressif, - e(n) est la réalisation d'un vecteur de bruits blancs gaussiens de 30 même dimension que 8(n). La prise en compte du modèle de mouvement dans le procédé d'estimation peut s'effectuer suivant différentes méthodes. Suivant un premier exemple de mise en oeuvre, on remplace dans la transformation spatio-temporelle W;(p,O), l'ensemble des paramètres 8(n) (angles de roulis cp(n), de tangage a(n) et de lacet ip(n)) par leur modèle paramétrique de mouvement du type donné dans l'expression (17). Le mouvement estimé est alors déterminé directement non pas en calculant l'ensemble des valeurs angulaires 8(n), mais en calculant les paramètres Àm du modèle de mouvement qui optimisent la fonction FI de ressemblance image et, le cas échéant, la fonction F2 de ressemblance externe. Dans une solution alternative, le mouvement estimé est déterminé par optimisation d'une fonction, dite « fonction F3 de ressemblance modèle ». La fonction F3 de ressemblance modèle associe, à une hypothèse de mouvement donnée, une ressemblance entre ladite hypothèse de mouvement donnée et le modèle de mouvement. Dans le cas d'un modèle de mouvement paramétrique du type donné à l'expression (17), la fonction F3 de ressemblance modèle pour le roulis cp(n) s'exprime par exemple sous la forme : ( M-1 \ F3(~P)=exp ûLh4 9(n)ûLÀm.pm(n).wm(n) (19) n \, m=0 où h4 est une p-norme. Les composantes identiques liées aux autres angles de tangage et de lacet sont cumulées par multiplication dans ce cas. Dans le cas d'un modèle de mouvement probabiliste du type donné à l'expression (18), la fonction F3 de ressemblance modèle s'écrira par exemple : ( L-1 F3 (0) = exp û L h4 0(n) û L al 0(n û I) (20) n 1=0 / / où h4 est une p-norme. L'optimisation de la fonction F3 de ressemblance modèle peut se faire conjointement avec l'optimisation de la fonction FI de ressemblance image, au moyen d'une fonction globale F cumulant ladite fonction de ressemblance image et ladite fonction de ressemblance modèle. De préférence, la fonction globale F cumule également la fonction F2 de ressemblance externe, et l'optimisation s'effectue conjointement pour les trois fonctions FI, F2 et F3. Suivant un autre exemple, on calcule l'hypothèse de mouvement O optimisant la fonction FI de ressemblance image, ou une fonction globale F cumulant la fonction FI de ressemblance image et la fonction F2 de ressemblance externe, ou une fonction globale F cumulant les trois fonctions FI, F2 et F3 (les paramètres du modèle étant fixés), puis, étant donnée cette hypothèse de mouvement, on calcule les paramètres du modèle (paramétrique ou probabiliste) optimisant la fonction F3 de ressemblance modèle pour cette hypothèse de mouvement. On peut itérer plusieurs fois cette optimisation alternée. La figure 4 représente un mode particulier de mise en oeuvre du procédé dans lequel le mouvement estimé C'EST est déterminé itérativement en tenant compte des images composées, des mesures effectuées par les capteurs externes et d'un modèle de mouvement probabiliste. Outre les étapes décrites précédemment en référence à la figure 3, le procédé d'estimation comporte également les étapes : - 205 de calcul de données, dites « données modèle », en fonction du modèle temporel de mouvement et de l'hypothèse de mouvement C'k de l'itération k en cours, - 203 mise à jour de l'hypothèse de mouvement pour l'itération suivante, en fonction des données modèle. Au cours de l'étape 205, on calcule des données modèle qui correspondent par exemple à la valeur de la fonction F3 de ressemblance modèle. De préférence, on détermine, au cours de l'étape 205 de calcul des données modèle, un ou plusieurs gradients de la fonction F3 de ressemblance modèle, afin d'avoir une information sur la direction dans laquelle faire évoluer l'hypothèse de mouvement C'k de l'itération k en cours pour optimiser la fonction F3 de ressemblance modèle. Les données modèle sont préférentiellement utilisées au cours de l'étape 203 de mise à jour de l'hypothèse de mouvement, tel que décrit précédemment pour la prise en compte des données externes. De préférence, le procédé comporte une 206 de mise à jour du modèle de mouvement en fonction de l'hypothèse de mouvement Ok+' pour l'itération (k+1) suivante. Au cours de cette étape 206, on met à jour le modèle de mouvement en fonction de l'hypothèse de mouvement Ok+' de l'itération (k+1). Cette mise à jour est une estimation d'un nouveau modèle de mouvement, par exemple une estimation au sens du maximum de vraisemblance, des moindres carrés, maximum a posteriori, etc. Masquaqe de pixels mis en correspondance Dans un mode particulier de mise en oeuvre du procédé, compatible avec l'un quelconque des modes précédemment décrits, on met en oeuvre un masquage des pixels mis en correspondance dont la ressemblance évolue, au cours de l'optimisation, de façon anormale par rapport à la ressemblance des autres pixels mis en correspondance. Par « masquage », on entend que les contributions desdits pixels dont la ressemblance évolue de façon anormale sont réduites, dans la fonction FI de ressemblance image, par rapport aux contributions des autres pixels. Les pixels mis en correspondance dont la ressemblance image évolue de façon anormale sont par exemple des pixels mis en correspondance dont les valeurs sont significativement différentes, en particulier dont les différences sont significativement plus importantes que celles statistiquement constatées, en moyenne, sur l'ensemble des pixels. En d'autres termes, les pixels mis en correspondance dont la ressemblance évolue de façon anormale correspondent à des mesures qui apparaissent comme aberrantes, compte tenu des mesures disponibles. De telles mesures aberrantes sont généralement liées à la présence de nuages, à des effets de parallaxe dus aux incertitudes sur le relief du fait que les lignes de visée des barrettes sont différentes, aux différences de domaines de longueurs d'ondes d'acquisition, à la prise en compte incorrecte du relief, etc. Le masquage des pixels mis en correspondance dont la ressemblance image évolue de façon anormale est par exemple effectué en considérant pour la fonction FI de ressemblance image une fonction hi d'estimation robuste. De telles fonctions d'estimation robuste, par exemple un M-estimateur, sont connues de l'homme de l'art. Suivant un autre exemple non limitatif, on détermine une matrice de pondération Mp des pixels des images composées. La matrice de pondération est préférentiellement de même dimension que les images composées. Certains éléments de la matrice de pondération peuvent être de valeur nulle (par exemple sur des zones hors recouvrement des images composées). Elle est par exemple mise en oeuvre comme suit : i F1 (0)=exp ûLLh,(MP(p) .(Bi(Wi(p,O))ûB;(p,Oe (21) i La figure 5 représente un mode particulier de mise en oeuvre du masquage par matrice de pondération, dans le cas où le mouvement estimé est déterminé itérativement, en tenant compte des images composées, des mesures effectuées par les capteurs externes et d'un modèle de mouvement probabiliste. Outre les étapes décrites en référence à la figure 4, le procédé d'estimation comporte également une étape 207 de mise à jour de la matrice de pondération Mkp pour l'itération k en cours, en fonction de l'hypothèse de mouvement Ok+1, étape 207 au cours de laquelle on détermine la matrice de pondération Mk+1p pour l'itération (k+1) suivante. De préférence, les éléments de la matrice de pondération Mk+1p sont calculés d'après les différences entre la valeur des pixels mis en correspondance d'après l'hypothèse de mouvement 0k+1, c'est-à-dire d'après les différences B1(W1(p,Ok+1)) û B2(W2(p,Ok+1)) Pour définir les coefficients de pondération, on peut par exemple calculer la répartition des différences de pixels mis en correspondance par intervalles de valeur, en déterminant pour chaque intervalle le nombre de couples de pixels dont la différence est dans cet intervalle, puis les pondérations de la matrice Mk+1p sont déterminées pour chaque pixel comme étant le rapport entre le nombre d'occurrences correspondant à l'intervalle dans lequel se situe le pixel et le nombre total de pixels dans les images. Il est à noter qu'on peut utiliser la même matrice de pondération pour toutes les paires d'images composées, ou déterminer une matrice de pondération pour chaque paire, en fonction des images composées de la paire. Schéma d'optimisation multi-résolution Dans encore un mode particulier de mise en oeuvre du procédé d'estimation, compatible avec l'un quelconque des modes précédemment décrits, ledit procédé d'estimation est itéré selon un schéma d'optimisation multi-résolution. En d'autres termes, lors d'une première itération du procédé d'estimation, on considère un premier sous-ensemble de pixels des images composées, et on détermine un premier mouvement estimé pour ce premier sous-ensemble. On considère dans une seconde itération du procédé un second sous-ensemble des pixels des images composées recalées avec le premier mouvement estimé, comportant davantage de pixels que le premier sous-ensemble, et on détermine un second mouvement estimé. On itère le procédé d'estimation en augmentant progressivement le nombre de pixels du sous-ensemble considéré, jusqu'à atteindre le sous-ensemble final, pouvant correspondre à l'intégralité des pixels des images composées, du moins dans la zone de recouvrement. En itérant le procédé d'estimation selon un schéma d'optimisation multi-résolution, on s'attend à améliorer la précision du mouvement estimé, et/ou à accélérer la convergence du procédé d'estimation vers un mouvement estimé qui satisfasse le critère d'arrêt. C) Remarques complémentaires, qénéralisations De manière plus générale, la portée de la présente invention ne se limite pas aux modes de mise en oeuvre décrits ci-dessus à titre d'exemples non limitatifs, mais s'étend au contraire à toutes les modifications à la portée de l'homme de l'art. Notamment, la détermination du mouvement estimé C'EST par optimisation de la fonction FI de ressemblance image et, le cas échéant, de la fonction F2 de ressemblance externe et/ou de la fonction F3 de ressemblance modèle, peut mettre en oeuvre tout procédé d'optimisation connu de l'homme de l'art. Le procédé peut s'appliquer à toute forme de barrettes, éventuellement autre que linéaire, par exemple linéaire discontinue, en V, ou tout autre forme géométrique, du moment qu'il existe un recouvrement entre les images composées d'au moins deux barrettes durant une prise de vue, et qu'il est possible d'exprimer un modèle géométrique de l'instrument d'observation 10 dans le repère instrument RINSTR. Il est à noter que les cellules sensibles des barrettes bi à bNb peuvent être agencées sur plusieurs lignes, le nombre de cellules sensibles par ligne étant cependant supérieur d'au moins un facteur dix au nombre de lignes. On comprend donc que les images lignes acquises par les différentes barrettes sont dans ce cas également constituées d'une pluralité de lignes de pixels. La présente invention n'en est pas modifiée, dans la mesure où les images composées sont formées à partir d'images lignes acquises à des instants différents, lesdites images lignes représentant des portions différentes de la scène observée à la surface de la Terre. En outre, rien n'exclut que les barrettes soient des lignes de cellules sensibles d'un même capteur se présentant sous la forme d'une matrice bidimensionnelle de cellules sensibles, à partir du moment où elles réalisent l'acquisition d'un même point de la scène observée à des instants différents. Rien n'exclut non plus de considérer des barrettes agencées non pas au plan focal d'un même instrument, mais au plan focal d'au moins deux instruments différents, dont les lignes de visée sont proches. On note que tout ce qui a été décrit reste valable si on désigne par « instrument », l'ensemble des au moins deux instruments, et par « ligne de visée de l'instrument », une ligne de visée éventuellement fictive par rapport à laquelle se déduisent les lignes de visée des au moins deux instruments à partir d'un modèle géométrique de l'ensemble des instruments. On note également que l'utilisation de toutes les images composées, et en particulier de toutes les paires possibles d'images composées, dans la fonction FI de ressemblance image, constitue un mode préféré de mise en oeuvre nullement limitatif de l'invention. Suivant un autre exemple, on utilise uniquement les images lignes acquises par les barrettes bl et bNb, c'est-à-dire les plus espacées dans le plan focal de l'instrument d'observation 10. S'il existe des instants d'acquisition pour lesquels on ne dispose pas d'image ligne acquise par l'une ou les deux barrettes bl et bNb, le mouvement estimé est par exemple interpolé pour lesdits instants d'acquisition, ou contraint par le modèle temporel de mouvement, qu'il soit paramétrique ou probabiliste. Suivant encore un autre exemple, on utilise toutes les images composées, mais on ne considère qu'un sous-ensemble de toutes les paires possibles d'images composées, de préférence un sous-ensemble dans lequel chaque image composée est utilisée au moins une fois. En plus des mouvements en rotation de l'instrument, on peut introduire les mouvements en translation de l'instrument d'observation 10 dans le mouvement à estimer. Il suffit d'augmenter le vecteur O d'hypothèse de mouvements de rotation de l'instrument en un vecteur (X,O) où X représente l'hypothèse de mouvement de translation au cours du temps, et exprimer les transformations spatio-temporelle et les fonctions de ressemblance FI, F2 et F3 en fonction de (X,O), puis à calculer l'hypothèse de mouvement (X,O) maximisant les fonctions de ressemblance. Dans le cas où l'on introduit les mouvements en translation parmi les variables à estimer, il est avantageux de considérer des mesures externes données par un capteur externe tel qu'un récepteur de navigation par satellites, une centrale inertielle, ou tout autre capteur permettant de fournir une estimation de la position du véhicule dans le référentiel terrestre RT, et de traiter ces mesures externes dans la fonction F2 de ressemblance externe comme décrit pour les mouvements en rotation. REVENDICATIONS1 - Procédé d'estimation du mouvement d'un instrument d'observation (10) embarqué dans un véhicule survolant un corps céleste, l'instrument d'observation (10) comportant au moins deux barrettes (bi,b2) dont les lignes de visée sont différentes, chaque barrette réalisant l'acquisition successive, par défilement de sa ligne de visée sur le corps céleste, d'images lignes formant une image composée, des parties des images composées représentant sensiblement une même scène du corps céleste, caractérisé en ce que le mouvement estimé est déterminé par optimisation d'une fonction de ressemblance image associant, à une hypothèse de mouvement donnée, une ressemblance entre des valeurs de pixels mis en correspondance d'au moins une paire d'images composées, des pixels mis en correspondance étant des pixels représentant théoriquement le même point de la scène observée compte tenu de l'hypothèse de mouvement considérée, et en ce que le procédé comporte, pour chaque hypothèse de mouvement considérée au cours de l'optimisation de la fonction de ressemblance image, les étapes : - (200) de détermination, pour l'au moins une paire d'images composées, d'au moins une transformation spatio-temporelle en fonction de l'hypothèse de mouvement considérée et d'un modèle géométrique de l'instrument d'observation (10), - (201) de mise en correspondance des pixels de l'au moins une paire d'images composées au moyen de l'au moins une transformation spatio-temporelle. 2 - Procédé selon la revendication 1, dans lequel la fonction de ressemblance image est optimisée itérativement, par itération des étapes suivantes jusqu'à ce qu'un critère d'arrêt soit vérifié : - (200) détermination de l'au moins une transformation spatio-temporelle de l'au moins une paire d'images composées d'après l'hypothèse de mouvement de l'itération en cours, - (201) mise en correspondance des pixels de l'au moins une paire d'images composées, - (202) calcul, pour la fonction de ressemblance image, de donnéesimages en fonction desdits pixels mis en correspondance, - (203) mise à jour de l'hypothèse de mouvement pour l'itération suivante en fonction des données images. 3 - Procédé selon la revendication 2, dans lequel les données image comprennent au moins un gradient de la fonction de ressemblance image, représentatif de la variation, au voisinage de l'hypothèse de mouvement de l'itération en cours, de la ressemblance entre les valeurs des pixels mis en correspondance de l'au moins une paire d'images composées. 4 - Procédé selon l'une des revendications précédentes, dans lequel le mouvement estimé est déterminé par optimisation d'une fonction de ressemblance externe associant, à une hypothèse de mouvement donnée, une ressemblance entre ladite hypothèse de mouvement et des mesures, représentatives de mouvements de l'instrument, effectuées par un ou des capteurs différents des barrettes de l'instrument d'observation (10). 5 - Procédé selon la revendication 4, dans lequel la fonction de ressemblance externe est optimisée itérativement, par itération des étapes suivantes jusqu'à ce qu'un critère d'arrêt soit vérifié : - (204) calcul, pour la fonction de ressemblance externe, de données externes en fonction des mesures du ou des capteurs et de l'hypothèse de mouvement de l'itération en cours, - (203) mise à jour de l'hypothèse de mouvement pour l'itération suivante en fonction des données externes. 6 - Procédé selon la revendication 5, dans lequel les données externes comprennent au moins un gradient de la fonction de ressemblance externe, représentatif de la variation, au voisinage de cette hypothèse de mouvement, de la ressemblance entre l'hypothèse de mouvement de l'itération en cours et les mesures du ou des capteurs. 7 - Procédé selon l'une des revendications 4 à 6, dans lequel la fonction de ressemblance image et la fonction de ressemblance externe sont optimisées conjointement, par optimisation d'une fonction globale cumulant ladite fonction de ressemblance image et ladite fonction de ressemblance externe. 8 - Procédé selon l'une des revendications 4 à 7, dans lequel au moins un capteur est un gyroscope ou un capteur stellaire ou un capteur angulaire d'un miroir orientable de l'instrument d'observation (10). 9 - Procédé selon l'une des revendications précédentes, dans lequel on détermine le mouvement estimé en fonction d'un modèle temporel prédéfini du mouvement de l'instrument d'observation (10). - Procédé selon la revendication 9, dans lequel le mouvement estimé est déterminé par optimisation d'une fonction de ressemblance modèle associant, à une hypothèse de mouvement donnée, une ressemblance 10 entre ladite hypothèse de mouvement et le modèle temporel de mouvement. 11 - Procédé selon la revendication 10, comportant une étape (206) de mise à jour du modèle temporel de mouvement en fonction d'une hypothèse de mouvement. 12 - Procédé selon l'une des revendications 10 à 11, dans lequel la fonction de ressemblance modèle est optimisée itérativement, par itération des étapes suivantes jusqu'à ce qu'un critère d'arrêt soit vérifié : - (205) calcul, pour la fonction de ressemblance de modèle, de données modèle en fonction du modèle temporel de mouvement et de l'hypothèse de mouvement de l'itération en cours, - (203) mise à jour de l'hypothèse de mouvement pour l'itération suivante en fonction des données modèle. 13 - Procédé selon la revendication 12, dans lequel les données modèle comprennent au moins un gradient de la fonction de ressemblance de modèle, représentatif de la variation, au voisinage de cette hypothèse de mouvement, de la ressemblance entre l'hypothèse temporelle de l'itération en cours et le modèle temporel de mouvement. 14 - Procédé selon l'une des revendications 10 à 13, dans lequel la fonction de ressemblance image et la fonction de ressemblance modèle sont optimisées conjointement, par optimisation d'une fonction globale cumulant ladite fonction de ressemblance image et ladite fonction de ressemblance modèle. 15 - Procédé selon l'une des revendications 9 à 14, dans lequel le modèletemporel de mouvement est un modèle de mouvement paramétrique, de préférence un modèle dont les fonctions de base sont du type polynômes, ou ondelettes, ou sinusoïdes, ou splines. 16 - Procédé selon l'une des revendications 9 à 14, dans lequel le modèle temporel de mouvement est un modèle de mouvement probabiliste, de préférence un modèle de type autorégressif. 17 - Procédé selon l'une des revendications précédentes, dans lequel on réduit, dans la fonction de ressemblance image, les contributions de pixels mis en correspondance dont la ressemblance évolue, au cours de l'optimisation, de façon anormale par rapport à la ressemblance des autres pixels mis en correspondance. 18 - Procédé selon la revendication 17, dans lequel la réduction des contributions de pixels mis en correspondance est effectuée au moyen d'une fonction d'estimation robuste. 19 - Procédé selon la revendication 17, dans lequel la réduction des contributions de pixels mis en correspondance est effectuée au moyen d'au moins une matrice de pondération. - Procédé selon la revendication 19, comportant, lorsque la fonction de ressemblance image est optimisée itérativement, une étape (207) de mise 20 à jour de l'au moins une matrice de pondération en fonction des contributions des pixels mis en correspondance à ladite fonction de ressemblance image.