Brevet FR3025683A1 — publié le 20160311

Procédé et dispositif d'estimation de paramètres d'égalisation d'image numérique La présente invention concerne un procédé d'estimation de paramètres d'égalisation d'image numérique. Elle concerne également un dispositif d'estimation de paramètres d'égalisation d'image numérique, un procédé d'évaluation de bruit d'égalisation et un dispositif d'évaluation de bruit d'égalisation associé. La présente invention se situe dans le domaine de l'amélioration de la qualité d'images numériques, en particulier des images présentant du bruit structuré directionnel composé de rayures, par exemple horizontales et/ou verticales. De tels artefacts sont dus à l'acquisition d'images, et se présentent en particulier pour des images acquises par des capteurs d'images à détecteurs multiples. Des capteurs d'images à détecteurs multiples sont utilisés notamment dans le domaine de la télédétection par satellite. En particulier, un mode d'acquisition d'images connu dans le domaine de la télédétection consiste à utiliser un système d'acquisition d'images comportant une ou plusieurs barrettes de détecteurs fixé à un véhicule spatial ou aérien, perpendiculairement à la direction de déplacement du véhicule. Ainsi, le système d'acquisition d'images défile au regard du paysage observé, chaque ligne de l'image numérique acquise correspondant à l'acquisition effectuée à un instant donné. L'image numérique obtenue est composée d'une matrice bidimensionnelle d'échantillons appelés pixels, chaque échantillon ayant une valeur de radiométrie associée. Dans ce mode d'acquisition, chaque colonne de l'image numérique obtenue est associée physiquement à un détecteur de la matrice monodimensionnelle de détecteurs. Ce mode d'acquisition par déplacement rectiligne des barrettes de détecteurs, perpendiculaire à la direction de déplacement du porteur, est connu sous le nom d'acquisition par balayage, ou « push-broom » en anglais. Dans un cas de figure idéal, l'acquisition d'un paysage uniforme devrait fournir une image numérique uniforme. Cependant, en pratique, les détecteurs présentent des différences de sensibilité qui se traduisent par des artefacts sur une image numérique acquise, qui sont des rayures verticales dans le cas de l'acquisition d'images en mode « push-broom » décrite ci-dessus. Il existe diverses méthodes de correction des artefacts dus aux différences de sensibilité entre détecteurs d'un capteur d'image à multiples détecteurs, ces différences de sensibilité ayant pour effet l'obtention d'une image comportant du bruit d'acquisition sous forme de rayures notamment horizontales/verticales, ou selon d'autres directions. On connaît notamment une méthode de correction de ces artefacts, dite d'égalisation d'image numérique, utilisant un modèle mathématique de correction, représentatif de la dispersion des réponses de l'ensemble des détecteurs, et le calcul des 3025683 2 paramètres d'une transformation permettant de déduire, à partir d'une valeur de radiométrie réellement observée, une valeur de radiométrie corrigée correspondant à une détection homogène pour laquelle l'image d'un paysage uniforme est uniforme, sans bruit d'acquisition. 5 Selon les techniques classiques, une méthode d'égalisation numérique d'image comporte une phase d'étalonnage permettant d'obtenir des valeurs de paramètres de transformation et une phase effective de mise en oeuvre de la correction par application de la transformation choisie avec les paramètres obtenus par étalonnage. Typiquement, pour l'acquisition d'images par télédétection, la phase d'étalonnage nécessite une 10 acquisition d'images uniformes sur la totalité des détecteurs, par exemple des images de nuit ou de paysages de neige, en conditions d'observation sans nuage, ce qui est difficile à obtenir, car de nombreuses prises de vue sont nécessaires. De plus, une phase de validation des valeurs de paramètres obtenus par un expert est également nécessaire. L'étalonnage classique en vue de l'égalisation d'image numérique est donc long, et 15 nécessite une intervention d'expert. Dans le domaine de la télédétection par satellite, une méthode d'étalonnage améliorée a été proposée dans l'article « AMETHIST : A Method for Equalization Thanks to HISTograms », de P. Kubik et V. Pascal, publié dans SPIE Proceedings vol. 5570, pages 256-267, novembre 2004. 20 Cette méthode propose l'acquisition d'une image d'étalonnage pour l'égalisation, par rotation du système d'acquisition d'images porté par le satellite, de manière à aligner sensiblement la ou les barrettes de détecteurs à la direction de déplacement de celui-ci. Ainsi, chaque détecteur fournit successivement une valeur d'échantillon d'image correspondant à une même zone de paysage, faisant partie d'un paysage quelconque 25 survolé par le satellite porteur du système d'acquisition. Un ensemble de paramètres d'égalisation, correspondant à un modèle d'égalisation de la réponse des détecteurs, est calculé à partir de l'image d'étalonnage acquise, après un redressement de celle-ci. La méthode proposée met en oeuvre une mise en correspondance d'histogrammes pour le calcul des paramètres d'égalisation d'image. Cette méthode permet d'obtenir un 30 ensemble de paramètres d'égalisation à partir d'une acquisition d'une seule image d'étalonnage d'égalisation, ne nécessitant pas l'acquisition d'une pluralité d'images uniformes. Cependant, les paramètres d'égalisation ainsi calculés sont adaptés à l'image d'étalonnage d'égalisation acquise et ne permettent pas d'apporter une amélioration de 35 qualité équivalente pour tout type d'images, en particulier lorsque les détecteurs 3025683 3 présentent une variabilité en fonction du spectre électromagnétique du rayonnement optique détecté. Un des objectifs de l'invention est d'améliorer l'estimation des valeurs des paramètres d'égalisation d'un modèle de correction lorsque les détecteurs présentent une 5 variabilité telle que décrite ci-dessus, et d'améliorer par conséquent la correction d'images numériques présentant des artefacts linéaires acquises par déplacement d'un système d'acquisition d'images à détecteurs multiples. Un autre objectif de la présente invention est de fournir une évaluation de bruit d'égalisation améliorée, permettant d'évaluer la performance d'un système d'acquisition à différents niveaux de radiométrie. 10 A cet effet, l'invention propose, selon un premier aspect, un procédé d'estimation de paramètres d'égalisation d'image numérique acquise par balayage d'une zone de paysage par un porteur d'un système d'acquisition d'images, le système d'acquisition d'images comprenant au moins une barrette de détecteurs, la ou les barrettes de détecteurs étant positionnées sensiblement perpendiculairement à la direction de 15 déplacement du porteur du système d'acquisition d'images, chaque détecteur étant apte à fournir une valeur de radiométrie représentative d'un rayonnement optique d'une dite zone de paysage. Le système d'acquisition d'images est pilotable de manière à ce qu'au moins un sous-ensemble de détecteurs voisins soit apte à fournir successivement une valeur d'échantillon d'image correspondant à une même zone de paysage. 20 Le procédé comporte les étapes suivantes : - acquisition d'au moins une image d'égalisation pour une catégorie de paysage choisie parmi un ensemble prédéterminé de catégories de paysage, chaque catégorie de paysage ayant un spectre électromagnétique associé, chaque image d'égalisation comprenant au moins une matrice d'échantillons d'image, chaque échantillon d'image 25 ayant une valeur de radiométrie associée correspondant à une même zone de paysage de la catégorie de paysage choisie, - pour au moins une catégorie de paysage choisie, calcul d'un ensemble de paramètres d'égalisation d'un modèle de correction paramétré prédéterminé pour chacun des détecteurs à partir des échantillons d'image d'égalisation acquise pour ladite 30 catégorie de paysage. Avantageusement, le procédé d'estimation de paramètres d'égalisation d'image numérique de l'invention permet d'estimer des paramètres d'égalisation adaptés pour un ensemble de catégories de paysage couvrant un spectre de longueurs d'ondes d'acquisition donné. 3025683 4 Le procédé selon l'invention peut présenter une ou plusieurs des caractéristiques ci-dessous, prise indépendamment ou en combinaison selon toutes combinaisons techniquement possibles. Le procédé comporte une étape de formation d'une portion d'image d'égalisation 5 redressée pour ladite catégorie de paysage choisie à partir des images d'égalisation acquises, chaque échantillon d'une ligne ou d'une colonne de ladite image d'égalisation redressée ayant une valeur de radiométrie obtenue par un desdits détecteurs pour une même zone de paysage balayé, chaque échantillon d'une colonne ou d'une ligne de ladite image d'égalisation redressée ayant une valeur de radiométrie obtenue par un même 10 détecteur, et une étape de calcul d'un ensemble de paramètres d'égalisation pour chacun des détecteurs à partir de ladite image d'égalisation redressée. Le procédé comprend une étape de formation d'image d'égalisation redressée par concaténation de portions d'image d'égalisation redressées pour chaque catégorie de paysage, chaque portion d'image d'égalisation redressée par catégorie de paysage ayant 15 un même nombre de colonnes. Le procédé comprend, pour chaque catégorie de paysage, le calcul d'un ensemble de paramètres d'égalisation par catégorie de paysage à partir de la portion d'image d'égalisation redressée obtenue pour ladite catégorie de paysage. En outre, le procédé comprend une étape de combinaison des ensembles de 20 paramètres d'égalisation par catégorie de paysage pour obtenir un ensemble final de paramètres d'égalisation applicable pour tout type de paysage. La formation d'une portion d'image d'égalisation redressée pour ladite catégorie de paysage choisie comporte une étape de redressement de chacune des images d'égalisation acquises par application d'un filtre de redressement. 25 Le filtre de redressement est un filtre associant, à un échantillon d'image courant de la portion d'image redressée, la valeur de radiométrie d'un échantillon d'image d'égalisation acquise sélectionné comme plus proche voisin dudit échantillon d'image courant selon un critère géométrique de voisinage prédéterminé. Le calcul de paramètres d'égalisation comprend des étapes de : 30 -calcul d'un histogramme cumulé normalisé à partir d'un sous-ensemble d'échantillons de l'image d'égalisation redressée correspondant à un détecteur du système d'acquisition d'images ; -calcul et mémorisation des quantiles dudit histogramme cumulé normalisé pour un ensemble prédéterminé de valeurs de fréquences associées. 3025683 5 Le procédé comprend une étape de calcul d'un ensemble de quantiles moyens pour chacune des valeurs de fréquence à partir des quantiles préalablement calculés pour chaque détecteur. Il comprend également, une étape de calcul, pour chaque détecteur, des 5 paramètres d'égalisation dudit modèle de correction paramétré prédéterminé, permettant de minimiser une erreur quadratique moyennes entre chacun des quantiles moyens et un quantile correspondant à la même fréquence que ledit quantile moyen, obtenu pour les valeurs de radiométrie corrigées en appliquant ledit modèle de correction. Le modèle de correction paramétré est un modèle affine par morceaux. 10 Selon un mode de réalisation, le pilotage du système d'acquisition d'images comporte une rotation de ladite au moins une barrette de détecteurs pour l'aligner sensiblement à ladite direction de déplacement du porteur du système d'acquisition d'images, de manière à ce que chaque détecteur soit apte à fournir successivement une valeur d'échantillon d'image correspondant à une même zone de paysage. 15 L'invention propose, selon un deuxième aspect, un procédé d'évaluation de bruit d'égalisation d'image numérique acquise par balayage d'une zone de paysage par un porteur d'un système d'acquisition d'images, le système d'acquisition d'images comprenant au moins une barrette de détecteurs, la ou les barrettes de détecteurs étant positionnées sensiblement perpendiculairement à la direction de déplacement du porteur 20 du système d'acquisition d'images, chaque détecteur étant apte à fournir une valeur de radiométrie représentative d'un rayonnement optique d'une dite zone de paysage. Ce procédé comporte les étapes de : -acquisition d'au moins une image d'égalisation pour une catégorie de paysage choisie parmi un ensemble prédéterminé de catégories de paysage, chaque catégorie de 25 paysage ayant un spectre électromagnétique associé, chaque image d'égalisation comprenant au moins une matrice d'échantillons d'image, chaque échantillon d'image ayant une valeur de radiométrie associée correspondant à une même zone de paysage de la catégorie de paysage choisie, - calcul d'un vecteur moyen d'échantillons, associé à un ensemble de sous- 30 ensembles de l'image d'égalisation acquise, chaque sous-ensemble correspondant à une catégorie de paysage et à un niveau de radiométrie moyen, chaque échantillon du vecteur moyen d'échantillons correspondant à un détecteur du système d'acquisition, - calcul d'un critère d'évaluation de niveau de bruit d'égalisation utilisant les vecteurs moyens calculés pour chaque dit sous-ensemble. 3025683 6 Le procédé d'évaluation de bruit d'égalisation d'image numérique selon l'invention peut présenter une ou plusieurs des caractéristiques ci-dessous, prise indépendamment ou en combinaison selon toutes combinaisons techniquement possibles. Le calcul d'un critère d'évaluation de niveau de bruit d'égalisation comporte une 5 étape d'application d'un ensemble de coefficients d'égalisation associé à un modèle d'égalisation prédéterminé. L'étape de calcul d'un vecteur moyen d'échantillons comporte un redressement de l'image d'égalisation acquise pour chaque catégorie de paysage, chaque échantillon d'une ligne ou d'une colonne de ladite image d'égalisation redressée ayant une valeur de 10 radiométrie obtenue par un desdits détecteurs pour une même zone de paysage balayé, chaque échantillon d'une colonne ou d'une ligne de ladite image d'égalisation redressée ayant une valeur de radiométrie obtenue par un même détecteur du système d'acquisition d'images. L'étape de calcul d'un vecteur moyen d'échantillons comporte des étapes de : 15 - calcul d'un histogramme cumulé normalisé à partir d'un sous-ensemble d'échantillons de l'image d'égalisation redressée correspondant à un détecteur du système d'acquisition d'images ; -calcul et mémorisation des quantiles dudit histogramme cumulé normalisé pour un ensemble prédéterminé de valeurs de fréquences associées, 20 - formation d'un vecteur moyen d'échantillons à partir des quantiles associés à une valeur de fréquence donnée pour chaque détecteur du système d'acquisition d'images. L'invention propose, selon un troisième aspect, un dispositif d'estimation de paramètres d'égalisation d'image numérique acquise par balayage d'une zone de paysage par un porteur d'un système d'acquisition d'images, le système d'acquisition 25 d'images comprenant au moins une barrette de détecteurs, la ou les barrettes de détecteurs étant positionnées sensiblement perpendiculairement à la direction de déplacement du porteur du système d'acquisition d'images, chaque détecteur étant apte à fournir une valeur de radiométrie représentative d'un rayonnement optique d'une dite zone de paysage, le système d'acquisition d'images étant pilotable de manière à ce qu'au 30 moins un sous-ensemble de détecteurs voisins soit apte à fournir successivement une valeur d'échantillon d'image correspondant à une même zone de paysage. Le dispositif d'estimation de paramètres d'égalisation d'image numérique comporte : - un module d'obtention d'au moins une image d'égalisation pour une catégorie de 35 paysage choisie parmi un ensemble prédéterminé de catégories de paysage, chaque catégorie de paysage ayant un spectre électromagnétique associé, chaque image 3025683 7 d'égalisation comprenant au moins une matrice d'échantillons d'image, chaque échantillon d'image ayant une valeur de radiométrie associée correspondant à une même zone de paysage de la catégorie de paysage choisie, - pour au moins une catégorie de paysage choisie, un module de calcul d'un 5 ensemble de paramètres d'égalisation d'un modèle de correction paramétré prédéterminé pour chacun des détecteurs à partir des échantillons d'image d'égalisation acquise pour ladite catégorie de paysage. Selon un quatrième aspect, l'invention concerne un programme d'ordinateur pour l'estimation de paramètres d'égalisation d'image numérique acquise par balayage d'une 10 zone de paysage par un porteur d'un système d'acquisition d'images, le système d'acquisition d'images comprenant au moins une barrette de détecteurs, la ou les barrettes de détecteurs étant positionnées sensiblement perpendiculairement à la direction de déplacement du porteur du système d'acquisition d'images, chaque détecteur étant apte à fournir une valeur de radiométrie représentative d'un rayonnement optique 15 d'une dite zone de paysage, le système d'acquisition d'images étant pilotable de manière à ce qu'au moins un sous-ensemble de détecteurs voisins soit apte à fournir successivement une valeur d'échantillon d'image correspondant à une même zone de paysage, ledit programme comportant des instructions pour mettre en oeuvre les étapes de : 20 - obtention d'au moins une image d'égalisation pour une catégorie de paysage choisie parmi un ensemble prédéterminé de catégories de paysage, chaque catégorie de paysage ayant un spectre électromagnétique associé, chaque image d'égalisation comprenant au moins une matrice d'échantillons d'image, chaque échantillon d'image ayant une valeur de radiométrie associée correspondant à une même zone de paysage 25 de la catégorie de paysage choisie, - pour au moins une catégorie de paysage choisie, calcul d'un ensemble de paramètres d'égalisation d'un modèle de correction paramétré prédéterminé pour chacun des détecteurs à partir des échantillons d'image d'égalisation acquise pour ladite catégorie de paysage. 30 L'invention propose, selon un cinquième aspect, un dispositif d'évaluation de bruit d'égalisation d'image numérique acquise par balayage d'une zone de paysage par un porteur d'un système d'acquisition d'images, le système d'acquisition d'images comprenant au moins une barrette de détecteurs, la ou les barrettes de détecteurs étant positionnées sensiblement perpendiculairement à la direction de déplacement du porteur 35 du système d'acquisition d'images, chaque détecteur étant apte à fournir une valeur de radiométrie représentative d'un rayonnement optique d'une dite zone de paysage. 3025683 8 Le dispositif d'évaluation de bruit d'égalisation d'image numérique comporte : -un module d'obtention d'au moins une image d'égalisation pour une catégorie de paysage choisie parmi un ensemble prédéterminé de catégories de paysage, chaque catégorie de paysage ayant un spectre électromagnétique associé, chaque image 5 d'égalisation comprenant au moins une matrice d'échantillons d'image, chaque échantillon d'image ayant une valeur de radiométrie associée correspondant à une même zone de paysage de la catégorie de paysage choisie, - un module calcul d'un vecteur moyen d'échantillons, associé à un ensemble de sous-ensembles de l'image d'égalisation acquise, chaque sous-ensemble correspondant 10 à une catégorie de paysage et à un niveau de radiométrie moyen, chaque élément du vecteur correspondant à un détecteur du système d'acquisition, - un module de calcul d'un critère d'évaluation de niveau de bruit d'égalisation utilisant les vecteurs moyens calculés pour chaque dit sous-ensemble. Selon un sixième aspect, l'invention concerne un programme d'ordinateur pour 15 l'évaluation de bruit d'égalisation d'image numérique acquise par balayage d'une zone de paysage par un porteur d'un système d'acquisition d'images, le système d'acquisition d'images comprenant au moins une barrette de détecteurs, la ou les barrettes de détecteurs) étant positionnées sensiblement perpendiculairement à la direction de déplacement du porteur du système d'acquisition d'images, chaque détecteur étant apte à 20 fournir une valeur de radiométrie représentative d'un rayonnement optique d'une dite zone de paysage, le système d'acquisition d'images étant pilotable de manière à ce qu'au moins un sous-ensemble de détecteurs voisins soit apte à fournir successivement une valeur d'échantillon d'image correspondant à une même zone de paysage, ledit programme comportant des instructions pour mettre en oeuvre les étapes de : 25 - obtention d'au moins une image d'égalisation pour une catégorie de paysage choisie parmi un ensemble prédéterminé de catégories de paysage, chaque catégorie de paysage ayant un spectre électromagnétique associé, chaque image d'égalisation comprenant au moins une matrice d'échantillons d'image, chaque échantillon d'image ayant une valeur de radiométrie associée correspondant à une même zone de paysage 30 de la catégorie de paysage choisie, - calcul d'un vecteur moyen d'échantillons, associé à un ensemble de sous-ensembles de l'image d'égalisation acquise, chaque sous-ensemble correspondant à une catégorie de paysage et à un niveau de radiométrie moyen, chaque élément du vecteur correspondant à un détecteur du système d'acquisition, 35 - calcul d'un critère d'évaluation de niveau de bruit d'égalisation utilisant les vecteurs moyens calculés pour chaque dit sous-ensemble. 3025683 9 D'autres caractéristiques et avantages de l'invention ressortiront de la description qui en est donnée ci-dessous, à titre indicatif et nullement limitatif, en référence aux figures annexées, parmi lesquelles : 5 -la figure 1 illustre schématiquement le principe d'acquisition d'image par un capteur à multiples détecteurs dans le cas de l'imagerie satellite ; -la figure 2 illustre schématiquement un profil d'acquisition d'image numérique selon une première orientation du système d'acquisition ; - les figures 3 et 4 illustrent schématiquement des profils de réponse spectrale de 10 détecteur en fonction de la longueur d'onde ; -la figure 5 illustre schématiquement un profil d'acquisition d'image numérique selon une deuxième orientation du système d'acquisition ; -la figure 6 est un schéma représentant les blocs fonctionnels d'un dispositif programmable apte à mettre en oeuvre l'invention ; 15 -la figure 7 est un synoptique d'un procédé d'estimation de paramètres d'égalisation d'image numérique selon un mode de réalisation de l'invention ; -la figure 8 illustre une portion d'image d'égalisation redressée ; - la figure 9 illustre une image d'égalisation redressée complète ; -la figure 10 illustre un modèle de correction d'image ; 20 -la figure 11 est un synoptique d'un procédé de calcul de coefficients d'égalisation selon un mode de réalisation de l'invention ; - la figure 12 est un synoptique d'un procédé d'évaluation de bruit d'égalisation selon un mode de réalisation de l'invention. 25 L'invention sera décrite ci-après plus particulièrement dans le domaine de l'acquisition d'image par satellite, par des barrettes de détecteurs multiples, dans le mode d'acquisition par le mode dit «push-broom », dans lequel, en mode opérationnel, chaque colonne de l'image numérique obtenue est associée physiquement à un détecteur de la barrette de détecteurs du système d'acquisition. 30 Le contexte de l'acquisition d'images numériques par satellite est illustré schématiquement dans la figure 1. Un porteur de système d'acquisition d'images, qui est un satellite dans cet exemple, ayant une trajectoire de déplacement T comporte un système d'acquisition d'images 2, comprenant une ou plusieurs matrices monodimensionnelles ou barrettes 4 35 de détecteurs 6, et un système optique 8. Chaque détecteur 6 réalise l'acquisition d'une portion Pu du paysage P observé, également appelée zone du paysage observé, par 3025683 10 intégration du rayonnement optique pendant un temps d'intégration pour fournir une valeur de radiométrie associée à la zone du paysage. En mode d'acquisition d'images opérationnel, on enregistre une matrice d'image numérique I telle qu'illustrée à la figure 2, comprenant des échantillons d'image 10 ou 5 pixels, chaque pixel ayant une valeur de radiométrie. Les valeurs de radiométrie sont codées sur un nombre donné de bits, par exemple 8, 12, 16 ou 32. La matrice d'image numérique I est composée de lignes de pixels Lig, et de colonne de pixels Col,, chaque colonne de pixels étant associée à un détecteur 6 d'une barrette de détecteurs, chaque ligne correspondant à l'acquisition d'une zone du paysage 10 à un instant donné. En pratique, sur le satellite Pléiades HR (marque déposée), le système d'acquisition 2 comprend cinq matrices de détecteurs CCD (« charge coupled device »), localisées dans le plan focal du système optique 8. En mode opérationnel, le système d'acquisition permet d'acquérir, en conditions d'acquisition nadir, un paysage P de 20 km 15 de largeur, chaque pixel de l'image numérique acquise correspondant à une largeur de terrain Pu de l'ordre de 70 cm. Etant donné que les détecteurs 6 peuvent avoir, en pratique, des sensibilités différentes, même lorsque le paysage sous-jacent est uniforme, des différences de valeurs de radiométrie sont observées, ce qui donne lieu à des artefacts visibles, qui sont 20 des rayures verticales dans ce cas, particulièrement visibles lorsque le paysage sous- jacent observé est uniforme. De plus, les détecteurs ont des profils légèrement différents en fonction de la longueur d'onde du rayonnement optique. Les figures 3 et 4 illustrent schématiquement le profil de réponse spectrale de deux détecteurs D1 et D2 en fonction de la longueur d'onde 25 du signal optique acquis. Il a été constaté qu'en pratique les profils de réponse spectrale sont légèrement décalés. En outre, les paysages observés ont des profils de luminance spectrale différents, et couvrent un spectre électromagnétique variable, comme illustré sur les figures 3 et 4, qui illustrent respectivement un profil spectral d'un paysage P1 d'une première catégorie de paysage et un profil spectral d'un paysage P2 d'une deuxième 30 catégorie de paysage. Ainsi, la variabilité des détecteurs combinée avec la variabilité des paysages observés peut induire de grandes différences quant aux artefacts à corriger sur les images acquises. L'égalisation d'image a pour objectif de réduire voire supprimer ces artefacts directionnels dus à la non-uniformité des détecteurs. 3025683 11 Un objectif de l'invention est de proposer un procédé d'estimation de valeurs de paramètres d'égalisation, pour un modèle d'égalisation d'image donné, permettant de remédier à la non-uniformité spectrale et radiométrique de la sensibilité des détecteurs. A cet effet, il est proposé d'acquérir une ou plusieurs images d'étalonnage pour 5 l'égalisation, appelées images d'égalisation ci-après, à utiliser dans une phase préalable d'estimation de paramètres d'égalisation. Une image d'égalisation pour l'étalonnage est acquise, dans ce mode de réalisation, par rotation de la ou des barrettes 4 de détecteurs 6, de manière à l'orienter sensiblement parallèlement à la direction T de déplacement du satellite, comme illustré 10 schématiquement à la figure 5. Une image bidimensionnelle est formée par juxtaposition des lignes acquises à chaque temps d'échantillonnage. Lorsque le système d'acquisition d'images orienté de cette manière défile sensiblement au droit d'un paysage P observé, chacun des détecteurs capte successivement un rayonnement optique correspondant sensiblement à une même zone 15 du paysage P. Ainsi, dans l'image d'égalisation IE obtenue, le pixel 12 a une valeur de radiométrie correspondant à une zone de paysage Pu (hachurée sur la figure 1) obtenue par le détecteur D1, le pixel 14 a une valeur de radiométrie correspondant à la même zone de paysage Pu obtenue par le détecteur D2 et ainsi de suite. 20 Il est à noter que le cas illustré à la figure 5 est schématique et correspond à une acquisition «idéale ». En pratique, le guidage n'est pas idéal pour tous les détecteurs, mais le modèle géométrique d'acquisition d'image pour chaque détecteur est connu et mémorisé, et permet ainsi de récupérer les pixels correspondant sensiblement à une même zone du paysage vue par chacun des détecteurs. 25 Selon l'invention, l'acquisition d'image d'égalisation IE est effectuée sur un paysage P appartenant à une catégorie de paysage donnée parmi un ensemble de catégories de paysage. Chaque catégorie de paysage comprend des paysages relativement homogènes sur une zone observée, les profils de luminance spectrale des différents paysages d'une 30 même catégorie étant sensiblement proches entre eux, en termes de profil spectral et de gamme de radiométrie couverte. A titre d'exemple non limitatif, l'ensemble suivant de huit catégories de paysage est considéré : {Sable, Forêt, Mer, Végétation, Neige, Bord de côte, Paysage urbain, Nuage}. . 35 Les procédés de traitement d'images d'égalisation pour le calcul de coefficients d'égalisation et d'évaluation de bruit d'égalisation de l'invention sont mis en oeuvre par un 3025683 12 dispositif programmable de type ordinateur, station de travail, dont les principaux blocs fonctionnels sont illustrés à la figure 6. Un dispositif programmable 20 apte à mettre en oeuvre les étapes de traitement pour d'estimation de paramètres d'égalisation ou pour l'évaluation de bruit d'égalisation 5 d'image numérique selon l'invention comprend une unité centrale de traitement 28, ou CPU, apte à exécuter des instructions de programme informatique lorsque le dispositif 20 est mis sous tension. Dans un mode de réalisation, une unité centrale de traitement multiprocesseurs est utilisée, permettant d'effectuer des calculs parallèles. Le dispositif 20 comporte également des moyens de stockage d'informations 30, par exemple des 10 registres, aptes à stocker des instructions de code exécutable permettant la mise en oeuvre de programmes comportant des instructions de code aptes à mettre en oeuvre les procédés selon l'invention. Le dispositif 30 comporte des moyens de commande 24 permettant de mettre à jour des paramètres et de recevoir des commandes d'un opérateur. Lorsque le dispositif 15 programmable 30 est un dispositif embarqué, les moyens de commande 24 comprennent un dispositif de télécommunications permettant de recevoir des commandes et des valeurs de paramètres à distance. Alternativement et de manière optionnelle, les moyens de commande 24 sont des moyens de saisie de commandes d'un opérateur, par exemple un clavier. 20 De manière optionnelle, le dispositif programmable 30 comprend un écran 22 et un moyen supplémentaire de pointage 26, tel une souris. Les divers blocs fonctionnels du dispositif 30 décrits ci-dessus sont connectés via un bus de communication 32. Un système d'estimation de paramètres d'égalisation d'image numérique selon 25 l'invention comporte un système d'acquisition d'images tel que décrit ci-dessus et un dispositif programmable comportant des modules aptes à mettre en oeuvre les étapes d'un procédé d'estimation de paramètres d'égalisation d'image numérique selon l'invention. Un système d'évaluation de bruit d'égalisation d'image numérique selon l'invention 30 comporte un système d'acquisition d'images tel que décrit ci-dessus et un dispositif programmable comportant des modules aptes à mettre en oeuvre les étapes d'un procédé d'évaluation de bruit d'égalisation d'image numérique selon l'invention. Les principales étapes d'un procédé d'estimation de paramètres d'égalisation d'image numérique selon un mode de réalisation de l'invention sont illustrées à la figure 7. 35 Lors d'une première étape 40, une image d'égalisation lEk est acquise pour un paysage Pk appartenant à une catégorie de paysage parmi un ensemble donné de 3025683 13 catégories de paysage prédéterminé. L'ensemble de catégories de paysage est sélectionné en fonction d'un objectif de qualité pour des applications utilisant les images acquises par le satellite par la suite, en mode opérationnel. L'ensemble de catégories de paysage est sélectionné pour couvrir un paysage moyen, défini en fonction des 5 applications visées, à la fois en termes de gamme de radiométrie et en spectre électromagnétique couvert. La taille de l'image d'égalisation lEk acquise est telle que le nombre de colonnes est égal au nombre de détecteurs noté Nd, et le nombre de lignes est paramétrable, de préférence compris entre Nd+1000 et (1+K)xNd, avec K un entier choisi, par exemple K=5. 10 Selon une variante, plusieurs images d'égalisation {lEk} pour des paysages {Pk} d'une même catégorie de paysage sont acquises à des moments et à des endroits différents, afin de diminuer l'influence des conditions météorologiques sur les données d'image acquises, ainsi que l'influence de l'incertitude sur la connaissance a priori du paysage réellement observé, en particulier son niveau de radiométrie, son homogénéité, 15 l'adéquation avec la catégorie visée. Dans cette variante, les images d'égalisation acquises sont concaténées pour former une image d'égalisation lEk pour la catégorie de paysage considérée. Ensuite, lors d'une étape de pré-traitement et de redressement 42, l'image d'égalisation lEk acquise est redressée pour former une portion IR' d'image d'égalisation 20 redressée. La figure 8 illustre schématiquement une d'image redressée IR, obtenue à partie de l'image IE illustrée à la figure 5. Dans ce mode de réalisation, l'image redressée IR contient, dans chaque ligne, la valeur de radiométrie obtenue pour une même zone de paysage Pu d'un paysage P par 25 chacun des détecteurs D1, D2 et ainsi de suite. De préférence, le nombre de colonnes de l'image d'égalisation redressée IR est égal au nombre total de détecteurs Nd et le nombre total de lignes de l'image IR, correspondant à un paysage d'une catégorie de paysage donnée, est de préférence compris entre 1000 et 4xNd. 30 En variante, on calcule une image d'égalisation redressée IR contenant, dans chaque ligne, la valeur de radiométrie obtenue pour une même zone de paysage Pu d'un paysage P par chacun des détecteurs. Le redressement appliqué pour obtenir une image d'égalisation redressée IR à partir d'une image d'égalisation IE est effectué par application d'un filtrage de 35 redressement, permettant de calculer la valeur de radiométrie d'un pixel Wu de l'image IR à partir de valeurs de pixels Uf(,),f(J) de l'image d'égalisation IE. 3025683 14 La fonction f de correspondance géométrique est déterminée en fonction des caractéristiques de déplacement du système d'acquisition par rapport au paysage et de l'orientation des détecteurs, ces caractéristiques étant obtenues à partir des consignes de guidage du satellite et de mesures a posteriori à partir d'équipements à bord du satellite, 5 et permet de calculer une correspondance entre grilles d'échantillonnage des images respectives IE et IR. Optionnellement, un modèle numérique de terrain fournissant une connaissance a priori de l'altitude du terrain correspondant au paysage est également utilisé. Les valeurs f(i) et f(j) de coordonnées dans la grille d'échantillonnage de l'image 1E 10 correspondant aux coordonnées entières i,j dans la grille d'échantillonnage de l'image IR sont, dans le cas général, des valeurs réelles. Dans le mode de réalisation préféré, le filtrage de redressement appliqué est un filtrage de type plus proche voisin : Wl = Un), jo , le critère de voisinage géométrique étant la sélection du plus proche voisin dans la même colonne, défini par: io = i et 15 jo = Int(f (j)) où Int(x) est la fonction donnant la valeur entière la plus proche du nombre réel x . Ainsi, le filtrage de redressement sélectionne la valeur de radiométrie du pixel le plus proche des coordonnées réelles io, f(j) pour obtenir la valeur de radiométrie W. De retour à la figure 7, les étapes 40 et 42 sont répétées pour chaque catégorie de 20 paysage de l'ensemble des catégories de paysage considérées. Dans un cas particulier, l'ensemble de catégories de paysage ne contient qu'une seule catégorie de paysage, par exemple une seule des catégories : {Sable, Forêt, Mer, Végétation, Neige, Bord de côte, Paysage urbain}. Dans ce cas, les paramètres d'égalisation obtenus seront particulièrement adaptés pour la catégorie de paysage 25 sélectionnée, mais non adaptés pour réaliser une égalisation de paysages d'une autre catégorie. Après itération des étapes 40, 42 sur l'ensemble de catégories de paysage considérées, on dispose d'un ensemble de portions d'image d'égalisation redressées IR', qui sont mémorisées. 30 Lors de l'étape d'obtention 44 d'image d'égalisation redressée, une image d'égalisation redressée complète IR est obtenue à partir de l'ensemble des portions d'images d'égalisation redressées. Dans un mode de réalisation, l'image d'égalisation redressée complète IR est obtenue par concaténation de portions d'image d'égalisation redressée IR' correspondant 35 à chaque catégorie de paysage Pk. 3025683 15 La figure 9 illustre une image d'égalisation redressée complète 50, formée de portions 52, 54, 56, 58 correspondant respectivement aux paysages P1 à P4, chacun des paysages P' à P4 faisant partie d'une catégorie de paysage d'un ensemble de catégories de paysage données. 5 De préférence, les portions 52 à 58 sont extraites des images IR' obtenues précédemment, en conservant uniquement des lignes pour lesquelles il n'y a pas de saturation, donc pour lesquelles les valeurs de radiométrie n'ont pas été tronquées, de manière à améliorer la fiabilité de l'estimation ultérieure des paramètres d'égalisation. De préférence, les portions respectives 52 à 58 présentent chacune un même 10 nombre de lignes, par exemple égal à 10000. Enfin, l'image d'égalisation redressée complète IR est utilisée pour la détermination de paramètres d'égalisation lors de l'étape 46, les paramètres d'égalisation étant mémorisés après calcul pour une utilisation ultérieure, notamment dans un procédé d'égalisation d'images acquises en mode opérationnel. 15 Les paramètres d'égalisation dépendent d'un modèle mathématique de correction, représentatif d'une détection homogène pour l'ensemble des détecteurs. Un modèle linéaire, affine par morceaux ou fonctionnel peut être envisagé. Le modèle mathématique de correction peut être choisi en fonction de caractéristiques connues du comportement des détecteurs du système d'acquisition d'images. 20 De manière générale, un modèle mathématique de correction paramétré F(Ici la) par un nombre donné de paramètres ou coefficients {c1}, permet d'obtenir des valeurs de radiométrie Y corrigées à partir de valeurs de radiométrie X observées. On note : Y -F(Icil,X) (Eq1) 25 Pour chaque détecteur D,, un ensemble de paramètres fc'd du modèle de correction est déterminé. Par exemple, comme illustré à la figure 10, un modèle mathématique de correction comportant deux morceaux affines, appelé également bilinéaire, est utilisé. La figure 10 représente le modèle bilinéaire de correction dans un repère 30 orthonormé comprenant en abscisse les valeurs de radiométrie X obtenues par le système d'acquisition d'images et en ordonnée les valeurs de radiométrie Y corrigées par égalisation. Pour chaque détecteur D,, le modèle bilinéaire illustré est paramétré par les paramètres de pente cli=a, et c'2=a2, où al est une première pente entre une valeur de 3025683 16 radiométrie Xmin et une valeur seuil Xs, et a2 une deuxième pente entre la valeur de radiométrie X, et la valeur de radiométrie maximale X'x. La fonction de correction appliquée est la suivante : Y = 0 pour X compris entre 0 et Xmin (Eq2) 5 Y= ai - (X - X Win) pour X compris entre X.. et X, (Eq3) Y= a2 - X +(a1 - a2) - Z, - Xffill, pour X compris entre X, et X. (Eq4) La valeur seuil X, est prédéterminée ou à estimer en même temps que les valeurs des paramètres al et a2, auquel cas cette valeur est un troisième paramètre c'3 du modèle. 10 Selon une première variante, une seule valeur de seuil X, est utilisée. Selon une deuxième variante, une valeur de seuil X', par détecteur D, est utilisée. Les valeurs des paramètres {c'1, c'2, c'3} sont déterminées pour chacun des détecteurs D, du système d'acquisition d'images. Plusieurs modes de réalisation de l'étape de détermination des paramètres 15 d'égalisation sont envisagés. Les principales étapes de la détermination des paramètres d'égalisation selon un premier mode de réalisation sont illustrées à la figure 11. Dans le mode de réalisation considéré, l'image d'égalisation redressée complète IR comporte un nombre L de lignes, et tous les pixels d'une même ligne correspondent à une 20 même zone de paysage, la valeur de radiométrie associée étant obtenue par un détecteur différent. Chaque colonne de l'image d'égalisation redressée complète comprend des valeurs de radiométrie obtenues sensiblement par un même détecteur. Les paramètres du modèle mathématique de correction à appliquer pour éliminer les artefacts dus à la variabilité des détecteurs sont calculés à partir de l'image 25 d'égalisation redressée complète IR. Lors d'une première étape d'initialisation 80, F valeurs notées ch à qF, dites valeurs de fréquences, comprises entre 0 et 1, sont déterminées. Le nombre n de valeurs est un paramètre du procédé, et peut être adapté à la catégorie de paysage. Pour une catégorie de paysage de basse radiométrie, une valeur de n faible est choisie, par 30 exemple n=4. De manière plus générale, n est de préférence compris entre 1 et 100, par exemple n=35 pour une dynamique de radiométrie codée sur 12 bits. Par exemple, des valeurs l {a ..1---1 a ,rn- - -,qF} régulièrement espacées entre 0,1 et 0,9 sont sélectionnées. En variante, les valeurs de fréquence qn sont calculées à partir d'un histogramme 35 moyen des colonnes de l'image d'égalisation redressée complète IR. 3025683 17 Ensuite, pour chacun des détecteurs D, parmi tous les détecteurs du système d'acquisition d'images (étape 82), un histogramme cumulé normalisé, noté HC,, est calculé et mémorisé à l'étape de calcul 84. En pratique, un histogramme cumulé HC,() des radiométries de la colonne de 5 l'image d'égalisation redressée cumulée IR associée au détecteur D, est calculé. De manière connue, la valeur au point V d'un histogramme cumulé est égale au nombre de pixels du signal d'image considéré, donc de la colonne de l'image IR considérée dans le cas présent, qui ont une valeur de radiométrie inférieure ou égale à V. L'histogramme cumulé est normalisé par le nombre total de pixels de la colonne 10 d'image. Par conséquent, les valeurs HC,(V) sont comprises entre 0 et 1. Ensuite, pour le détecteur D, courant, on considère chaque valeur de fréquence qn parmi l'ensemble de valeurs de fréquence {q1,...,qF} à l'étape 86, et on détermine le quantile Xig" correspondant à la fréquence qn dans l'histogramme cumulé normalisé HC, (): 15 HC,(X q'') qn (Eq5) La valeur de quantile déterminée est mémorisée, par exemple dans un tableau. Lorsque les valeurs de quantile pour chacune des valeurs de fréquence qn pour le détecteur courant D, sont obtenues, le procédé retourne à l'étape 82 de sélection d'un détecteur suivant Di+, 20 Ainsi, les étapes 82 à 88 sont répétées pour tous les détecteurs du système d'acquisition d'images 2. Lorsque tous les détecteurs ont été traités, une étape 90 de calcul de quantiles moyens est mise en oeuvre. Lors de cette étape, pour chaque valeur de fréquence qn, on calcule un quantile 25 moyen, représentatif de la moyenne des quantiles déterminés détecteur par détecteur : 1 Nd Zq" Xig" . (Eq 6) où Nd est le nombre de détecteurs Nd i=1 Ensuite, lors de l'étape 92, chaque détecteur D, est considéré, et un ensemble de paramètres d'égalisation {QI} est calculé par rapport à l'ensemble des quantiles moyens préalablement calculés et mémorisés, par une méthode d'optimisation de type moindres 30 carrés dans ce mode de réalisation. Dans ce mode de réalisation, les paramètres d'égalisation {QI} sont choisis comme étant les paramètres qui minimisent l'erreur quadratique entre les quantiles calculés des valeurs de fréquence qn après application de la fonction mathématique de correction et les valeurs moyennes de quantiles ou quantiles moyens calculés. 3025683 18 Si on note, pour le détecteur D,, Y=F({ca,X) le signal de radiométrie corrigé, et Y " le quantile calculé après correction pour la fréquence qn, les paramètres {ca sont calculés pour minimiser la somme S(c') : F2 S({C; (Eq7) n=1 \, 5 Où Zq' est la valeur moyenne des quantiles pour la fréquence qn calculée par la formule (Eq6), F ({c;}, X tq' ) et un est la moyenne de l'écart-type des radiométries d'une colonne pour la valeur radiométrique moyenne Zq' pondérée par le nombre de lignes de IR. La valeur un est déterminée préalablement à la mise en oeuvre du procédé d'estimation, c'est un paramètre caractéristique du bruit d'acquisition de chaque 10 détecteur. Selon une variante, la valeur de un est calculée en utilisant un modèle paramétré en fonction de la valeur radiométrique Zq" . Dans ce mode de réalisation, un ensemble de paramètres d'égalisation {QI} définissant la fonction de correction de radiométrie est obtenu pour l'image d'égalisation redressée complète IR, cet ensemble de paramètres d'égalisation étant optimisé pour 15 l'ensemble des catégories de paysages considéré représenté dans IR. Selon un deuxième mode de réalisation, un ensemble de paramètres d'égalisation {cil par détecteur D, et par catégorie de paysage Pk est calculé. Ces paramètres sont destinés à être utilisés pour un paysage de la catégorie de paysage correspondante. Dans ce deuxième mode de réalisation, les étapes 80 à 92 sont effectuées pour 20 chaque détecteur D, et chaque catégorie de paysage Pk, donc il n'y a pas de formation d'image redressée complète par concaténation de portions d'image d'égalisation redressées. L'image utilisée en entrée est l'image lEk obtenue pour la catégorie de paysage Pk. Selon un troisième mode de réalisation, les étapes 80 à 90 sont effectuées pour 25 chaque détecteur et chaque paysage observé, pour obtenir un ensemble de paramètres {ci r}pour chaque détecteur D, et chaque paysage Pr. Les ensembles de paramètres sont ensuite fournis en entrée de l'étape 92. En sortie de l'étape 92, un ensemble de paramètres global pour chaque détecteur est alors obtenu, par exemple en appliquant une moyenne pondérée de paramètres de même rang obtenus pour un détecteur D, 30 donné, pour chaque paysage Pr traité. Selon un autre aspect, l'invention a pour objet un procédé d'évaluation de bruit d'égalisation sur chaque catégorie de paysage et à différents niveaux de radiométrie. 3025683 19 Avantageusement, un bilan précis de performance du système d'acquisition d'images peut être effectué en vol, donc lors d'une phase opérationnelle, pour différents niveaux de radiométrie et différentes catégories de paysage. Les méthodes connues d'évaluation de bruit d'égalisation ne permettent pas une caractérisation des performances d'un système 5 d'acquisition aussi fine. La figure 12 illustre un procédé d'utilisation des paramètres d'égalisation estimés pour l'évaluation de bruit d'égalisation pour plusieurs niveaux moyens de radiométrie et plusieurs catégories de paysage. Le procédé est composé de trois étapes principales, une étape 100 d'acquisition 10 d'une ou plusieurs images d'égalisation par catégorie de paysage considérée, analogue à l'étape d'acquisition d'image d'égalisation 40 décrite ci-dessus, une étape 102 de calcul d'un ensemble de sous-ensembles constitués de lignes moyennes ou vecteurs moyens de paysage, chaque ligne moyenne de paysage correspondant à un niveau radiométrique moyen, vu successivement par l'ensemble de détecteurs du système d'acquisition, et une 15 étape 104 d'évaluation de bruit d'égalisation par niveau radiométrique moyen et par catégorie de paysage. L'étape 102 est composée de deux étapes dans le mode de réalisation de la figure 12, une étape 106 de redressement de la ou des images d'égalisation acquises lors de l'étape d'acquisition 100, et une étape 108 de calcul des quantiles Xiq' de chaque 20 catégorie de paysage correspondant à la fréquence qn, chaque quantile Xig" correspondant à un quantile moyen Zq' indiquant un niveau de radiométrie associé au quantile Xiq". Chaque ligne moyenne de paysage avant égalisation est un vecteur moyen monodimensionnel, comprenant autant d'éléments que de détecteurs, chaque élément étant une valeur de quantile Xiq" correspondant à un quantile moyen Zq" 25 De préférence, les étapes 80 à 88 détaillée ci-dessus en référence à la figure 11 sont appliquées pour le calcul des quantiles Xiq". L'étape 104 d'évaluation de bruit d'égalisation par niveau radiométrique moyen et par type de paysage comprend une première étape optionnelle 110 d'application de paramètres d'égalisation associés au modèle de correction Y -F(cii,X) pour égaliser 30 chaque ligne moyenne de paysage correspondant à un niveau de radiométrie Zq' et à une catégorie de paysage. Selon un premier mode de réalisation les paramètres d'égalisation {c1} sont estimés par le procédé d'estimation décrit ci-dessus. 3025683 20 Selon un deuxième mode de réalisation, des paramètres d'égalisation prédéfinis, obtenus au préalable par une méthode d'estimation, sont appliqués. Ainsi, le procédé permet d'évaluer la performance du système d'acquisition combiné à un jeu de paramètres d'égalisation donné. 5 En variante, l'égalisation est effectuée au préalable avec un ensemble de paramètres d'égalisation prédéfini, et l'étape 108 est directement suivie de l'étape 112 décrite ci-après. L'étape 110, lorsqu'elle est effectuée, est suivie d'une étape 112 d'application d'un critère d'évaluation du niveau de bruit. 10 Dans un mode de réalisation, l'évaluation de niveau de bruit de la ligne moyenne égalisée se fait par une formule mathématique d'évaluation de la dispersion après égalisation de la ligne moyenne de paysage considérée. De préférence, une ligne moyenne égalisée est calculée pour chaque quantile Yiq" , associé un niveau de radiométrie moyen Zq' pour une catégorie de paysage 15 donnée. Dans un mode de réalisation, le niveau de bruit d'égalisation pour un signal d'image Y est calculé par : C(X) = Max stat ( Am(i) RMS _col(i)s (Eq 8) La fonction Max stat donne la valeur maximale statistiquement à P pour cent des 20 valeurs d'un vecteur, avec P=99,7% par exemple. Am(i) = max (Yiq")- min (Yiq" ) (Eq 9) 1\1, est un nombre de colonnes (typiquement Nc=50) et RMS col(k) est une valeur connue par ailleurs, caractéristique de la moyenne, sur les colonnes i à i+Nc-1, des écarts-types des radiométries d'une colonne pour la valeur radiométrique moyenne Zq" . 25 Bien évidemment, le critère C(-)) peut être remplacé par tout autre critère d'évaluation analogue. On peut donc calculer une évaluation Cn (-) pour chaque quantile, associé un niveau de radiométrie moyen Zq' pour une catégorie de paysage donnée. Ces valeurs d'évaluation peuvent ensuite être exploitées par exemple pour 30 cartographier les performances de l'égalisation d'image, par rapport aux valeurs de radiométrie et aux catégories de paysage. Ainsi, il est possible d'évaluer finement la 3025683 21 performance des paramètres d'égalisation calculés, sur l'ensemble du domaine de radiométries et de paysages d'intérêt. REVENDICATIONS1.- Procédé d'évaluation de bruit d'égalisation d'image numérique acquise par balayage d'une zone de paysage par un porteur d'un système d'acquisition d'images (2), le système d'acquisition d'images (2) comprenant au moins une barrette (4) de détecteurs (6), la ou les barrettes (4) de détecteurs (6) étant positionnées sensiblement perpendiculairement à la direction de déplacement du porteur du système d'acquisition d'images (2), chaque détecteur (6) étant apte à fournir une valeur de radiométrie représentative d'un rayonnement optique d'une dite zone de paysage, caractérisé en ce qu'il comporte les étapes de : -acquisition (100) d'au moins une image d'égalisation pour une catégorie de paysage choisie parmi un ensemble prédéterminé de catégories de paysage, chaque catégorie de paysage ayant un spectre électromagnétique associé, chaque image d'égalisation comprenant au moins une matrice d'échantillons d'image, chaque échantillon d'image ayant une valeur de radiométrie associée correspondant à une même zone de paysage de la catégorie de paysage choisie, - calcul (102) d'un vecteur moyen d'échantillons, associé à un ensemble de sous-ensembles de l'image d'égalisation acquise, chaque sous-ensemble correspondant à une catégorie de paysage et à un niveau de radiométrie moyen, chaque échantillon du vecteur moyen d'échantillons correspondant à un détecteur du système d'acquisition, - calcul (104) d'un critère d'évaluation de niveau de bruit d'égalisation utilisant les vecteurs moyens calculés pour chaque dit sous-ensemble. 2.- Procédé d'évaluation de bruit d'égalisation d'image numérique selon la revendication 1, caractérisé en ce que le calcul (104) d'un critère d'évaluation de niveau de bruit d'égalisation comporte une étape d'application d'un ensemble de coefficients d'égalisation associé à un modèle d'égalisation prédéterminé. 3.- Procédé d'évaluation de bruit d'égalisation d'image numérique selon l'une des revendications 1 ou 2, caractérisé en ce que l'étape de calcul (102) d'un vecteur moyen d'échantillons comporte un redressement (106) de l'image d'égalisation acquise pour chaque catégorie de paysage, chaque échantillon d'une ligne ou d'une colonne de ladite image d'égalisation redressée (IR') ayant une valeur de radiométrie obtenue par un desdits détecteurs pour une même zone de paysage balayé, chaque échantillon d'une colonne ou d'une ligne de ladite image d'égalisation redressée (IR') ayant une valeur de radiométrie obtenue par un même détecteur du système d'acquisition d'images. 3025683 23 4.- Procédé d'évaluation de bruit d'égalisation d'image numérique selon la revendication 3, caractérisé en ce que l'étape de calcul (102) d'un vecteur moyen d'échantillons comporte des étapes de : 5 - calcul d'un histogramme cumulé normalisé à partir d'un sous-ensemble d'échantillons de l'image d'égalisation redressée correspondant à un détecteur du système d'acquisition d'images ; -calcul et mémorisation des quantiles dudit histogramme cumulé normalisé pour un ensemble prédéterminé de valeurs de fréquences associées, 10 - formation d'un vecteur moyen d'échantillons à partir des quantiles associés à une valeur de fréquence donnée pour chaque détecteur du système d'acquisition d'images. 5.- Dispositif d'évaluation de bruit d'égalisation d'image numérique acquise par balayage d'une zone de paysage par un porteur d'un système d'acquisition d'images (2), 15 le système d'acquisition d'images (2) comprenant au moins une barrette (4) de détecteurs (6), la ou les barrettes (4) de détecteurs (6) étant positionnées sensiblement perpendiculairement à la direction de déplacement du porteur du système d'acquisition d'images (2), chaque détecteur (6) étant apte à fournir une valeur de radiométrie représentative d'un rayonnement optique d'une dite zone de paysage, 20 caractérisé en ce qu'il comporte : -un module d'obtention d'au moins une image d'égalisation pour une catégorie de paysage choisie parmi un ensemble prédéterminé de catégories de paysage, chaque catégorie de paysage ayant un spectre électromagnétique associé, chaque image d'égalisation comprenant au moins une matrice d'échantillons 25 d'image, chaque échantillon d'image ayant une valeur de radiométrie associée correspondant à une même zone de paysage de la catégorie de paysage choisie, - un module calcul d'un vecteur moyen d'échantillons, associé à un ensemble de sous-ensembles de l'image d'égalisation acquise, chaque sous-ensemble correspondant à une catégorie de paysage et à un niveau de radiométrie moyen, chaque élément du 30 vecteur correspondant à un détecteur du système d'acquisition, - un module de calcul d'un critère d'évaluation de niveau de bruit d'égalisation utilisant les vecteurs moyens calculés pour chaque dit sous-ensemble. 6.- Programme d'ordinateur comportant des instructions de code de programme 35 pour l'exécution des étapes du procédé d'évaluation de bruit d'égalisation d'image 3025683 24 numérique selon les revendications 1 à 4 lorsque ledit programme est exécuté par un dispositif programmable. 5