Brevet FR3014588A1 — publié le 20150612

PROCEDE DE RESTITUTION DU MOUVEMENT HAUTE FREQUENCE DE LA LIGNE DE VISEE D'UN IMAGEUR Le domaine de l'invention est celui de la restitution de micro-vibrations d'un imageur monté sur un porteur tel qu'un satellite soumis à des perturbations d'attitude lors de son déplacement. Avec l'amélioration de la résolution des instruments optiques, les micro-vibrations peuvent devenir problématiques. Pour restituer les images acquises par un imageur monté sur un satellite, il est nécessaire de connaître avec précision l'orientation et le mouvement de sa ligne de visée. lo Dans le domaine satellitaire, le SCAO acronyme de l'expression Système de Commande d'Attitude et d'Orbite, permet de restituer l'attitude du satellite, dans le domaine des basses fréquences (0-4hz). Les perturbations de la ligne de visée de plus hautes fréquences aussi désignées micro-vibrations, dégrade la qualité des images sous forme notamment de distorsion 15 géométrique, et ce d'autant plus que ces images sont à haute résolution ; ces perturbations sont difficilement modélisables et donc difficilement prévisibles. Dans le cas où un imageur en vol se trouve confronté à des perturbateurs hautes fréquences, la solution consiste en une correction a 20 posteriori au cours d'une étape de traitement des images acquises. Des méthodes d'estimation de ces micro-vibrations à partir des images acquises par l'imageur ont été développées. Ces méthodes exploitent la particularité de l'imageur à observer la même scène à plusieurs instants différents ; l'imageur fonctionne en mode « pousse-balai » (ou 25 « push-broom » en anglais) et comporte différentes barrettes de détecteurs telles que les images (image 1, image 2) respectivement issues de ces barrettes (barrette 1, barrette 2) se recouvrent totalement ou partiellement, comme illustré figure 1. Cette caractéristique est présente 1. lors de l'utilisation d'un plan focal « à plat ». Dans ce type de plan 30 focal, l'ensemble des barrettes de détecteurs nécessaire pour couvrir toute la fauchée, c'est-à-dire toute la largeur de la scène, doivent être positionnées en quinconce comme sur la figure 1, pour des questions d'aménagements et de reconstruction de la ligne image. Ainsi les images issues de deux barrettes proches comportent une zone de recouvrement appelé « zone inter-barrette »(ZIB). Dans ces ZIBs un même point au sol est vu par une barrette paire et une barrette impaire à deux instants différents. 2. Un autre exemple de configuration est le cas des barrettes XS telles qu'installées sur la plateforme PLEIADE HAUTE RESOLUTION ou PHR. Les micro-vibrations du satellite et donc de l'imageur engendrent des distorsions géométriques dans les images qui sont mesurées sous forme 10 de décalages d'alignement entre les points homologues des différentes barrettes ; ces décalages peuvent être mesurés par des méthodes de type corrélation ou flot optique. Ces décalages sont aussi appelés disparités et l'ensemble des disparités sur une zone est désignée carte de disparités. La première étape de traitement d'image montrée figure 2, 15 exploite les zones de recouvrement (ZIBs dans l'exemple) afin d'en calculer les cartes de disparités engendrées par les micro-vibrations pvib(t). Cette étape d'estimation de la carte de disparités s'appuie sur une méthode classique de type corrélation dense entre images acquises par les barrettes. A partir de ces décalages ou plus précisément de cette carte de 20 disparités, l'étape suivante d'intégration temporelle permet de remonter à l'information de micro-vibration pvib(t) ; elle est réalisée par résolution d'un problème inverse de type déconvolution. Une solution de type filtre de Wiener comme décrit dans la demande de brevet EP 1 843 295, mène à des performances peu probantes comme illustré sur la figure 3a qui montre une 25 variation d'une micro-vibration haute fréquence ainsi estimée par rapport à la micro-vibration réelle, cette variation pouvant atteindre environ 25%. En conséquence, il demeure à ce jour un besoin pour un procédé de traitement de ces images permettant d'améliorer ces performances. 30 L'invention qui s'inscrit dans le cadre de l'étape d'intégration temporelle, est basée sur une méthode de type déconvolution comportant une étape de rotation dans le plan complexe de la carte de disparités. Plus précisément l'invention a pour objet, un procédé de 35 traitement d'un signal acquis par un imageur monté sur un porteur et soumis à des micro vibrations lors de son mouvement, l'imageur comportant des barrettes de photo détecteurs décalées entre elles dans le plan focal dans le sens du mouvement, et présentant deux à deux des zones de recouvrement dans le sens du mouvement, l'imageur étant apte à acquérir des images selon un mode « push-broom », le procédé comportant : une étape d'estimation d'une carte de disparités provoquées par les micro vibrations, à partir de l'acquisition d'au moins deux images et de mesures de décalage entre des zones de recouvrement de ces images, et une étape d'estimation des micro vibrations à partir de la carte de disparités. Il est principalement caractérisé en ce que l'étape d'estimation des micro vibrations comprend : un découpage du domaine fréquentiel des micro vibrations en 15 portions g) contenant chacune la fréquence d'une unique harmonique de micro vibration, pour chaque portion 0 du domaine fréquentiel : o calcul d'un déphasage optimal à partir d'un critère statistique sur la transformée de Fourier de la carte de 20 disparités réduite à cette portion S2 de domaine fréquentiel, o application de ce déphasage optimal à la carte de disparité correspondante, o passage dans l'espace de Fourier pour obtenir le spectre de cette carte de disparités déphasée, et sélection de la portion 25 de spectre correspondant à la portion S2 de domaine fréquentiel, o déconvolution de la portion de spectre sélectionnée au moyen d'un filtre de déconvolution, o retour dans l'espace direct par transformée de Fourier 30 inverse, et application du déphasage optimal inverse pour obtenir la micro vibration temporelle correspondant à la portion S2 de domaine fréquentiel considérée, - une concaténation des micro vibrations obtenues pour chaque portion de domaine fréquentiel. 35 Cela permet une amélioration notable de la qualité de restitution des perturbations d'attitude (ou micro-vibrations). Le déphasage optimal est avantageusement obtenu en maximisant l'écart de la distribution de la TF de la carte de disparités à une loi gaussienne (distribution caractéristique d'un bruit gaussien). Cet écart peut être calculé à partir du kurtosis ou d'une fonction des moments d'ordre supérieur. Selon une caractéristique de l'invention, les barrettes comportent 10 au moins une barrette panchromatique et une barrette multi-spectrale. Les barrettes peuvent être décalées entre elles en quinconce. Les perturbations d'attitude (ou micro-vibrations) présentent typiquement une fréquence supérieure à 10 Hz. L'invention a aussi pour objet un produit programme d'ordinateur, 15 ledit programme d'ordinateur comprenant des instructions de code permettant d'effectuer les étapes du procédé tel que décrit précédemment, lorsque ledit programme est exécuté sur un ordinateur. D'autres caractéristiques et avantages de l'invention apparaîtront 20 à la lecture de la description détaillée qui suit, faite à titre d'exemple non limitatif et en référence aux dessins annexés dans lesquels : la figure 1 déjà décrite représente schématiquement un exemple de configuration des barrettes d'un imageur, la figure 2 déjà décrite illustre les principales étapes d'un procédé 25 de restitution de micro-vibrations selon l'état de la technique, les figures 3 illustrent les performances d'un procédé de restitution d'une micro-vibration haute fréquence pvib(t) selon l'état de la technique (fig 3a), et celles d'un exemple de procédé de restitution selon l'invention (fig 3b), avec pour chaque figure la micro-vibration réelle (en trait continu) et la micro- 30 vibration restituée (points restitués sous forme de ronds), la figure 4 illustre les principales sous-étapes de l'étape d'intégration temporelle selon l'invention, la figure 5 montre schématiquement un exemple de module du spectre de carte de disparité comportant deux harmoniques et donc deux 35 portions de domaine fréquentiel qui doivent être traités séparément, la figure 6 montre schématiquement un exemple d'impact du déphasage sur l'histogramme de la transformée de Fourier de la carte de disparité (en pointillés : contribution du bruit, en trait plein : contribution d'une microvibration), la figure 7 montre un exemple de mesure de kurtosis (écart à distribution gaussienne) de la distribution de la transformée de Fourier de la carte de disparité en fonction du déphasage appliqué, permettant la détermination du déphasage optimal. D'une figure à l'autre, les mêmes éléments sont repérés par les 10 mêmes références. L'étape d'estimation de la carte de disparités s'appuie sur une méthode classique connue de l'homme du métier de type recherche de points homologues entre des images acquises par des barrettes différentes, 15 ou plus précisément dans les zones de recouvrement de ces images. Habituellement, on dispose dans le plan focal différents types de barrettes, par exemple des barrettes panchromatiques (PAN) et des barrettes multi-spectrales (XS). Elles sont généralement disposées en quinconce comme montré sur la figure 1 mais d'autres configurations 20 peuvent être envisagées : elles pourraient aussi être alignées et/ou inclinées l'une par rapport à l'autre, etc. Cette estimation est basée sur les mesures de disparités échantillonnées à haute ou très haute cadence (> 10 kHz). La disparité entre un couple d'image s'écrit : 25 D(k*Te) = pvib(k*Te) - pvib(t-AT(k*Te)), où AT(t) représente l'écart de temps entre la visibilité d'un même point au sol vu par les deux barrettes, Te représente le temps d'échantillonnage d'une ligne image, k un entier non nul. L'effet de l'écartement temporel entre les barrettes est donc celui 30 d'un filtre coupe-bande de la forme [1-6(t- AT(t))], où t = k*Te.(Te étant le temps d'échantillonnage d'une ligne image). Cet écartement temporel AT(t) qui dépend de l'écartement spatial dans le sens du mouvement de l'imageur, est généralement supposé constant, mais le procédé selon l'invention s'applique également dans le cas d'un écartement temporel variable. La carte de disparités peut par exemple être obtenue en appliquant les étapes de « Mise en correspondance des images » et d' « Estimation des variations angulaires » décrites dans la demande de brevet EP 1 843 295. L'invention qui s'inscrit dans le cadre de l'étape d'intégration temporelle, est basée sur une méthode d'intégration temporelle de type déconvolution du filtre coupe-bande rappelé ci-dessus, qui comporte une étape de rotation dans le plan complexe des données associées à la carte de disparités. Les performances insuffisantes de la méthode de l'état de la technique proviennent notamment de ce que le filtre de convolution comporte de nombreux zéros dans le cas considéré, c'est-à-dire aux fréquences de la forme K/AT, K étant un entier. Avant d'effectuer la dé-convolution, une étape de rotation dans le plan complexe est effectuée, comme montré figure 4 : A1) Le domaine fréquentiel des micro vibrations est découpé en portions n contenant chacune la fréquence d'une unique harmonique. Ce découpage est fait sur la base d'une détection des harmoniques des micro vibrations par exemple par seuillage ; selon une variante les portions sont prédéterminées, connues a priori. Dans l'exemple de la figure, on détecte deux harmoniques à environ 180 Hz et 280 Hz ; le domaine fréquentiel sera donc découpé en deux portions par exemple n1 =20 Hz-250Hz et Q2=250 Hz-750 Hz. Les étapes suivantes sont effectuées de façon indépendante pour 25 chacune de ces portions. A2) Pour chaque portion S2 de domaine fréquentiel, la représentation du spectre de la carte de disparités est la plus avantageuse lorsque le déphasage appliqué maximise l'écart entre la distribution du spectre et une gaussienne, comme illustré figure 6. Le bruit supposé gaussien a une 30 distribution gaussienne indépendamment du déphasage. En effet, on modélise les disparités selon la formule : D = Dpvib + bruit où « Dpvib » représente la disparité due aux micro vibrations et « bruit » représente le bruit de mesure des disparités. Après les étapes de 35 transformée de Fourier et de déphasage cette équation devient : TF(Deph(D,9)) = TF(Deph(Dpvib, (p)) + TF(Deph(bruit, (p)) où « TF » représente l'opération de transformée de Fourier, et Deph(D, (p) : D(k*Te) exp(i*k*(p)*D(k*Te) représente l'opération de déphasage d'un angle (p. Le troisième terme TF(Deph(bruit, (p)) a une distribution gaussienne quel que soit (p. Le second terme TF(Deph(Dpvib, (p)) a une distribution non gaussienne qui dépend de (p. Ainsi lorsque la distribution du premier terme TF(Deph(D,(p)) s'éloigne d'une gaussienne, c'est que la contribution du second terme augmente. Lorsque l'écart à une gaussienne est maximal, le déphasage (popt est optimal dans le sens où la micro vibration est la plus visible. Classiquement l'écart d'une distribution à une loi gaussienne peut être mesuré par les moments d'ordre supérieur, et par exemple par les moments d'ordre quatre normalisés appelés « kurtosis » selon : (popt = max(p( kurtosis( TFo(Deph(D,(p)) ) où TF0(f) représente la TF de f à la portion de spectre considéré. Un exemple de mesure de kurtosis en fonction du déphasage d'un cas typique est montré figure 7. Dans cet exemple, le déphasage optimal vaut -2.8rad. Un grand nombre de mesures issues de l'analyse en composantes indépendantes peuvent également être utilisées, notamment les mesures développées par Pierre Common (P. Comon. « Independent component analysis, a newconcept? » Signal Processing, 36(3):287-314, 1994). A l'issue de cette étape on a obtenu plusieurs déphasages optimaux : un déphasage optimal pour chaque portion de spectre correspondant au domaine fréquentiel SZ considéré. B1) Pour chaque portion S2 de domaine fréquentiel, on applique le déphasage optimal obtenu à la carte de disparité correspondante. B2) On passe dans l'espace de Fourier pour obtenir le spectre de cette carte de disparités déphasée. B3) On sélectionne la portion correspondante de ce spectre. C) La dé-convolution est appliquée dans l'espace déphasé S dec =TF(Deph(D, (p_opt)) * Filtre_deconvolution ; on peut utiliser comme filtre de déconvolution le filtre de Wiener comme dans notre exemple, ou le filtre inverse régularisé ou non calculé à partir de la réponse impulsionnelle du filtre initial coupe bande. D1) Le résultat de cette déconvolution est ramené dans l'espace direct en 5 appliquant successivement la transformée de Fourier inverse puis le déphasage optimal inverse. On obtient alors pvib(t) selon : pvibn(k*Te) = exp(-i*k*(popt)*TF-1(S_dec)(k*Te) pour la portion g) de domaine fréquentiel considérée. D2) Enfin les résultats de la dé-convolution de chaque portion 52 sont 10 concaténés sur l'ensemble des portions û considérées : pvib(k*Te) = Somme{Û} [ pvibn(k*Te) ]. En appliquant le procédé selon l'invention tel que décrit, on a obtenu les performances illustrées sur la figure 3b qui montre une variation 15 de la micro-vibration haute fréquence ainsi estimée par rapport à la micro- vibration réelle, inférieure à 1 %. L'imageur peut bien sûr être embarqué sur un porteur autre qu'un satellite, tel qu'un autre aéronef, un véhicule terrestre ou maritime, etc. 20 Ce procédé de traitement d'image peut notamment s'implémenter à partir d'un produit programme d'ordinateur, ce programme d'ordinateur comprenant des instructions de code permettant d'effectuer les étapes du procédé de reconstruction. Il est enregistré sur un support lisible par 25 ordinateur. Le support peut être électronique, magnétique, optique, électromagnétique ou être un support de diffusion de type infrarouge. De tels supports sont par exemple, des mémoires à semi-conducteur (Random Access Memory RAM, Read-Only Memory ROM), des bandes, des disquettes ou disques magnétiques ou optiques (Compact Disk - Read Only 30 Memory (CD-ROM), Compact Disk - Read/Write (CD-R/W) et DVD). REVENDICATIONS1. Procédé de traitement d'un signal acquis par un imageur monté sur un porteur et soumis à des micro vibrations lors de son mouvement, l'imageur comportant des barrettes de photo détecteurs décalées entre elles dans le plan focal dans le sens du mouvement, et présentant deux à deux des zones de recouvrement dans le sens du mouvement, l'imageur étant apte à acquérir des images selon un mode « push-broom », le procédé comportant : - une étape d'estimation d'une carte de disparités provoquées par les micro vibrations, à partir de l'acquisition d'au moins deux images et de mesures de décalage entre des zones de recouvrement de ces images, et une étape d'estimation des micro vibrations à partir de la carte de disparités, caractérisé en ce que l'étape d'estimation des micro vibrations comprend : un découpage du domaine fréquentiel des micro vibrations en portions Q contenant chacune la fréquence d'une unique harmonique (étape A1), pour chaque portion S2 du domaine fréquentiel : o calcul d'un déphasage optimal à partir d'un critère statistique sur la distribution de la transformée de Fourier de la carte de disparités réduite à cette portion ) de domaine fréquentiel (étape A2), o application de ce déphasage optimal à la carte de disparité correspondante (étape B1), o passage dans l'espace de Fourier (étape B2), et sélection de la portion de spectre de la carte de disparité déphasée correspondant à la portion s) de domaine fréquentiel (étape B3), o déconvolution de la portion de spectre sélectionnée au moyen d'un filtre de déconvolution (étape C),o retour dans l'espace direct par transformée de Fourier inverse, et application du déphasage optimal inverse pour obtenir la micro vibration temporelle correspondant à la portion Q de domaine fréquentiel considérée (étape Dl), une concaténation des micro vibrations obtenues pour chaque portion Q de domaine fréquentiel (étape D2). 2. Procédé de traitement d'un signal selon la revendication précédente, caractérisé en ce que le déphasage optimal est obtenu en maximisant un écart entre la distribution de la TF de la carte de disparités et une gaussienne. 3. Procédé de traitement d'un signal selon la revendication précédente, caractérisé en ce que cet écart est obtenu par une combinaison linéaire de moments d'ordre supérieur normalisée. 4. Procédé de traitement d'un signal selon la revendication précédente, caractérisé en ce que cet écart est obtenu par le kurtosis. 5. Procédé de traitement d'un signal selon l'une des revendications précédentes, caractérisé en ce que les barrettes comportent au moins une barrette panchromatique et une barrette multi-spectrale. 6. Procédé de traitement d'un signal selon l'une des revendications précédentes, caractérisé en ce que les barrettes sont décalées entre elles en quinconce. 7. Procédé de traitement d'un signal selon l'une des revendications précédentes, caractérisé en ce que les micro vibrations présentent une fréquence supérieure à 10 Hz. 8. Procédé de traitement d'un signal selon l'une des revendications précédentes, caractérisé en ce que le porteur est un satellite en orbite. 9. Un produit programme d'ordinateur, ledit programme d'ordinateur comprenant des instructions de code permettant d'effectuer les étapesdu procédé selon l'une quelconque des revendications précédentes, lorsque ledit programme est exécuté sur un ordinateur.5