FR 3066301 A1 20181116 FR 1754234 A 20170515 1 DOMAINE TECHNIQUE La présente invention appartient au domaine des images optiques acquises par un instrument d'observation optique embarqué à bord d'un engin spatial, tel qu'un satellite, ou d'un engin aérien évoluant à haute altitude, et concerne plus particulièrement un procédé de traitement d'une image en vue de corriger les effets liés à la présence d'aérosols et/ou de nuages, lors de l'acquisition de ladite image, dans l'atmosphère entre l'engin spatial ou aérien et la scène observée représentée par ladite image. ÉTAT DE LA TECHNIQUE En effet, il est connu que la présence d'aérosols (brouillard, sable, etc.) ou de nuages peu épais au-dessus d'une scène observée depuis un satellite se traduit, dans des images mesurant le rayonnement reçu de la scène dans des bandes de longueurs d'ondes notamment optiques, par la présence d'un voile qui modifie la perception de ladite scène. La présence d'aérosols ou de nuages entraîne également un effet dit d'environnement (« adjacency effect » dans la littérature anglo-saxonne) qui se traduit par une image floutée du fait de la diffusion du rayonnement sur les particules présentes dans l'atmosphère. Il est donc souhaitable de pouvoir traiter de telles images afin d'essayer de supprimer en tout ou partie ledit voile et l'effet d'environnement. Il va de soi que cela est tout particulièrement intéressant dans les régions de la Terre qui souffrent de conditions atmosphériques difficiles, comme par exemple des régions tropicales qui, lorsqu'elles ne sont pas entièrement recouvertes par des nuages lors de la saison des pluies, sont fréquemment affectées par des brumes lors de la saison sèche. Ce type de traitement est applicable de manière plus générale dans toute région de la Terre, notamment lorsqu'il s'agit de produire des images mosaïques de zones de grandes dimensions à partir de plusieurs images, pour pouvoir exploiter des images de qualité moindre et ainsi accélérer la finalisation de l'image mosaïque. De tels traitements visant à supprimer les effets liés à la présence de d'aérosols et/ou de nuages sont utiles non seulement pour améliorer le rendu visuel d'une image, mais également pour fournir une représentation plus 3066301 2 robuste et normalisée qui permet in fine d'améliorer l'efficacité d'éventuels traitements ultérieurs, comme par exemple une classification de la scène, etc. En outre, il est souhaitable de pouvoir effectuer ces traitements de manière autonome, c'est-à-dire sans avoir besoin d'informations externes sur 5 les conditions atmosphériques dans lesquelles une image a été acquise et sans avoir à faire appel à une autre image représentant la même scène. Afin de corriger les effets de l'atmosphère sur une image acquise par un satellite, il est connu notamment de rechercher les pixels les plus sombres de l'image, et de supposer que ces pixels représentent des objets noirs dans la 10 scène (approche dite « dark prior » dans la littérature anglo-saxonne). Ces pixels sont alors utilisés pour caractériser le voile présent dans l'image et le supprimer. Une telle approche est par exemple décrite dans l'article scientifique : « Single Image Haze Removal Using Dark Channel Prior », Kaiming He, Jian Sun, Xiaoou Tang, CVPR 2009. 15 Toutefois, la détection automatique de pixels noirs est parfois compliquée et, surtout, de tels pixels n'existent pas toujours dans une image. Une alternative consiste à utiliser des objets de référence dans la scène, dont la réflectance est connue a priori, et de les utiliser pour caractériser le voile présent dans l'image. Cependant, il n'est pas toujours possible d'avoir de tels 20 objets de référence dans une scène, et leur détection peut en outre être rendue difficile par la présence d'aérosols et/ou de nuages dans l'atmosphère. Il est également connu de la demande internationale WO 2014/197021 de combiner l'approche précédente avec un modèle atmosphérique modélisant la propagation dans l'atmosphère en présence d'aérosols. Le modèle 25 atmosphérique utilisé est de type LOWTRAN. Le procédé décrit recherche automatiquement les pixels sombres dans l'image et estime de manière itérative l'épaisseur optique des aérosols (« aerosol optical depth » ou AOD dans la littérature anglo-saxonne). En supposant que les pixels sombres correspondent à des objets noirs, le procédé estime l'épaisseur optique des 30 aérosols en tenant compte de l'effet d'environnement. Toutefois, l'hypothèse de présence d'un objet noir n'est pas toujours vérifiée, de sorte que le traitement appliqué n'est pas toujours efficace. En outre, le procédé décrit ne permet pas de distinguer les aérosols des nuages, 3066301 3 alors que ceux-ci produisent des effets différents. D'une part, les effets des aérosols sont très marqués dans les longueurs d'ondes bleues alors que les effets des nuages sont relativement achromatiques. En outre, les aérosols et les nuages produisent des effets d'environnement différents du fait qu'ils se 5 situent à des altitudes différentes et ont des comportements de diffusion différents. EXPOSÉ DE L'INVENTION La présente invention a pour objectif de remédier à tout ou partie des limitations des solutions de l'art antérieur, notamment celles exposées ci-avant, 10 en proposant une solution qui permette de supprimer en tout ou partie les effets liés à la présence d'aérosols et/ou de nuages (non opaques) dans une image acquise par un instrument d'observation optique embarqué à bord d'un engin spatial, tel qu'un satellite, ou d'un engin aérien évoluant à haute altitude. En outre, la présente invention vise à proposer une solution qui 15 permette d'effectuer de tels traitements de manière autonome, c'est-à-dire sans avoir à recevoir d'informations a priori sur les conditions atmosphériques dans lesquelles l'image a été acquise et sans avoir à faire appel à une autre image représentant la même scène acquise à une date antérieure. A cet effet, et selon un premier aspect, l'invention concerne un 20 procédé de traitement d'une image optique acquise par un engin spatial ou aérien, ladite image représentant un rayonnement dans au moins une bande de longueurs d'ondes reçu d'une scène à la surface de la Terre, comportant : - une obtention d'un modèle de réflectance de la scène et d'un modèle atmosphérique, ledit modèle atmosphérique comportant au 25 moins un modèle d'aérosols et un modèle de nuages modélisant les effets de l'atmosphère dans des couches respectives d'altitudes respectives différentes, - une détermination de valeurs optimales de paramètres du modèle atmosphérique permettant d'optimiser ledit modèle atmosphérique 30 par rapport à l'image, - une détermination de valeurs de paramètres de correction pour chaque couche du modèle atmosphérique en fonction des valeurs optimales des paramètres du modèle atmosphérique, 3066301 4 - une conversion de l'image en une image au sol avec effets d'environnement en fonction des valeurs des paramètres de correction de chaque couche du modèle atmosphérique, - une conversion de l'image au sol avec effets d'environnement en 5 image au sol sans effets d'environnement en fonction des valeurs des paramètres de correction de chaque couche du modèle atmosphérique. Ainsi, le procédé de traitement d'image utilise à la fois un modèle de réflectance au sol qui modélise la réflectance intrinsèque de la scène observée, 10 et ne nécessite pas de détecter de pixels sombres dans l'image et de supposer que ceux-ci représentent des objets noirs de la scène. En outre, le modèle atmosphérique utilisé est un modèle multicouche comportant au moins deux couches d'altitudes respectives différentes, à savoir une couche modélisée par un modèle d'aérosols qui restitue les effets des 15 aérosols et, au-dessus, une couche modélisée par un modèle de nuages qui restitue les effets des nuages. Par conséquent, les effets des aérosols d'une part et les effets des nuages d'autre part sont modélisés par deux modèles différents, en prenant en compte les altitudes respectives auxquelles se trouvent les aérosols (par 20 exemple entre 0 et environ 2 kilomètres) et les nuages (par exemple entre environ 2 kilomètres et environ 4 kilomètres) dans l'atmosphère. Il en résulte que les effets respectifs des aérosols et des nuages peuvent être mieux estimés et corrigés dans l'image acquise par l'engin spatial ou aérien. Le procédé de traitement comporte alors une conversion de l'image en 25 image au sol avec effets d'environnement, c'est-à-dire en image au sol dans laquelle les effets d'environnement n'ont pas encore été corrigés, et une conversion de ladite image au sol avec effets d'environnement en image au sol sans effets d'environnement, c'est-à-dire en image au sol dans laquelle les effets d'environnement ont été corrigés. 30 Dans des modes particuliers de mise en oeuvre, le procédé de traitement d'image peut comporter en outre l'une ou plusieurs des caractéristiques suivantes, prises isolément ou selon toutes les combinaisons techniquement possibles. 3066301 5 Dans des modes particuliers de mise en oeuvre, les paramètres de correction comportent des réflectances intrinsèques et des atténuations intrinsèques de chaque couche du modèle atmosphérique, et la conversion en image au sol avec effets d'environnement comporte une correction de ce 5 réflectances intrinsèques et atténuations intrinsèques des différentes couches du modèle atmosphérique. Dans des modes particuliers de mise en oeuvre, les paramètres de correction comportent des proportions d'environnement transmis pour chaque couche du modèle atmosphérique, et la conversion en image au sol sans effets 10 d'environnement comporte une correction de contributions d'environnement introduites par les différentes couches en fonction des proportions d'arrière-plan transmis. Dans des modes particuliers de mise en oeuvre, la correction d'une contribution d'environnement en un pixel pour une couche comporte : 15 - un calcul d'une réflectance reçue dans un cône de vue dont le sommet correspond au pixel considéré, - une suppression de ladite réflectance reçue dans le cône de vue, pondérée par la proportion d'effet d'environnement transmis pour la couche considérée et le pixel considéré. 20 Dans des modes particuliers de mise en oeuvre, le procédé de traitement d'image comportant une identification de pixels de l'image pour lesquels un critère d'opacité de la couche du modèle de nuages est vérifié, dits « pixels masqués », et la valeur de chaque pixel masqué est forcée, dans l'image au sol sans effets d'environnement, à une valeur prédéterminée. 25 En effet, en cas de nuage opaque, il n'est pas possible de restituer la scène sous le nuage, de sorte qu'il n'est pas nécessaire de corriger les valeurs de pixels représentant la scène masquée sous le nuage. Dans des modes particuliers de mise en oeuvre, la valeur prédéterminée correspond à la valeur du pixel masqué dans l'image après 30 conversion en valeur en sommet de nuage (« top of cloud » dans la littérature anglo-saxonne) ou en valeur en sommet d'atmosphère (« top of atmosphere » dans la littérature anglo-saxonne). Dans des modes particuliers de mise en oeuvre, la détermination de 3066301 6 valeurs optimales de paramètres du modèle atmosphérique comporte une détermination, pour chaque paramètre, d'une carte de valeurs optimales dudit paramètre associées respectivement à différents pixels de l'image, et un filtrage spatial de ladite carte de valeurs optimales dudit paramètre. 5 Dans des modes particuliers de mise en oeuvre, pour chaque paramètre, la détermination d'une carte de valeurs optimales dudit paramètre comporte en outre une détermination d'une carte de valeurs de confiance associées aux valeurs optimales, et le filtrage spatial d'une carte de valeurs optimales est effectué en fonction de la carte de valeurs de confiance associée. 10 Dans des modes particuliers de mise en oeuvre, le filtrage spatial d'une carte de valeurs optimales comporte un calcul de cartes de valeurs optimales à résolution dégradée pour différentes valeurs de dégradation de résolution, en fonction de la carte de valeurs optimales et de la carte de valeurs de confiance associée, et un calcul d'une moyenne pondérée de la carte de 15 valeurs optimale et des cartes de valeurs optimales à résolution dégradée. Dans des modes particuliers de mise en oeuvre, l'image correspond à une image brute dont les valeurs des pixels ont été converties en valeurs en sommet d'atmosphère. Dans des modes particuliers de mise en oeuvre, l'image correspond à 20 une image haute résolution dont la résolution a été dégradée, comportant une détermination de facteurs de correction à appliquer à l'image pour obtenir l'image au sol sans effets d'environnement, et une conversion de l'image haute résolution en une image haute résolution au sol sans effets d'environnement en fonction desdits facteurs de correction. 25 Dans des modes particuliers de mise en oeuvre, le modèle atmosphérique comporte en outre un modèle de ciel modélisant les effets de l'atmosphère dans une couche d'altitude supérieure à la couche du modèle de nuages. Dans des modes particuliers de mise en oeuvre, le procédé de 30 traitement comporte une étape préalable de détermination si l'image comporte un nuage, un modèle atmosphérique sans modèle de nuages étant utilisé lorsqu'il est déterminé que l'image ne comporte pas de nuage. Selon un second aspect, l'invention concerne un produit programme 3066301 7 d'ordinateur comportant un ensemble d'instructions de code de programme qui, lorsqu'elles sont exécutées par un processeur, configurent ledit processeur pour mettre en oeuvre un procédé de traitement d'image selon l'un quelconque des modes de mise en oeuvre de l'invention. 5 Selon un troisième aspect, l'invention concerne un dispositif de traitement d'image acquise par un engin spatial ou aérien comportant des moyens configurés pour mettre en oeuvre un procédé de traitement d'image selon l'un quelconque des modes de mise en oeuvre de l'invention. PRÉSENTATION DES FIGURES 10 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 schématique d'un satellite d'observation en orbite terrestre, 15 - Figure 2: un diagramme représentant les principales étapes d'un procédé de traitement d'image, - Figure 3 : un diagramme représentant un mode préféré de mise en oeuvre de procédé de traitement d'image, - Figure 4: une représentation schématique des différentes couches 20 d'un exemple de modèle atmosphérique, - Figure 5: un diagramme représentant le fonctionnement pyramidal d'un filtrage spatial d'une carte de valeurs optimales d'un paramètre du modèle atmosphérique, - Figure 6: un diagramme représentant un exemple de calcul d'une 25 réflectance reçue dans un cône de vue, - Figure 7: un diagramme représentant une variante de mise en oeuvre du procédé de traitement d'image de la figure 3. Dans ces figures, des références identiques d'une figure à une autre désignent des éléments identiques ou analogues. Pour des raisons de clarté, 30 les éléments représentés ne sont pas à l'échelle, sauf mention contraire. DESCRIPTION DÉTAILLÉE DE MODES DE RÉALISATION Dans la suite de la description, on considère de manière non limitative le cas d'un traitement d'images acquises depuis un engin spatial de type 3066301 8 satellite. Il convient toutefois de préciser que l'invention est également applicable à des images acquises par un instrument d'observation optique embarqué à bord d'un engin aérien (avion, ballon, drone, etc.) évoluant à une altitude suffisamment élevée pour que ladite image soit exposée aux effets de 5 voile et d'environnement de la colonne atmosphérique sous-jacente. La figure 1 représente schématiquement un satellite 10 en orbite terrestre, configuré pour observer une scène à la surface de la Terre. L'orbite terrestre est par exemple une orbite défilante, de sorte que la scène observée par le satellite varie au cours du temps, du fait au moins du déplacement dudit 10 satellite 10 par rapport à la Terre. De manière conventionnelle, le satellite 10 embarque au moins un instrument d'observation 11 optique qui permet d'acquérir une image correspondant au rayonnement reçu dans une bande de longueurs d'ondes. La bande de longueurs d'ondes d'acquisition de l'image se trouve dans le 15 domaine des longueurs d'ondes optiques qui correspond au domaine des longueurs d'ondes visibles et infrarouges, notamment proche infrarouge (« Near Infrared » ou NIR) ou infrarouge à courte longueur d'onde (« Short Wavelength Infrared » ou SWIR). Généralement, le ou les instruments d'observation 11 embarqués par 20 le satellite 10 comportent un ensemble de détecteurs sensibles dans des bandes de longueurs d'ondes respectives différentes. L'image acquise d'une scène est alors multi-spectrale ou hyper-spectrale (si le nombre de bandes de longueurs d'ondes est important) et comporte alors plusieurs images élémentaires de la scène acquises dans des bandes de longueurs d'ondes 25 respectives différentes. Par exemple, il est possible d'avoir à la fois : - une image élémentaire dans la bande des longueurs d'ondes bleues (par exemple [450-520] nanomètres), - une image élémentaire dans la bande des longueurs d'ondes vertes (par exemple [530-590] nanomètres), 30 - une image élémentaire dans la bande des longueurs d'ondes rouges (par exemple [625-695] nanomètres). - une image élémentaire dans la bande des longueurs d'ondes proche infrarouge (par exemple [760-890] nanomètres). 3066301 9 L'instrument d'observation 11 peut également, suivant d'autres exemples, être un instrument panchromatique et acquérir directement une image dans une bande de longueurs d'ondes large, qui couvre par exemple directement toutes les longueurs d'ondes visibles. 5 La présente invention est applicable à une image au sens large, c'est- à-dire qu'elle est applicable sur une image (élémentaire) unique acquise dans une seule bande de longueurs d'ondes large ou étroite. Toutefois, la présente invention trouve une application particulièrement avantageuse dans le cas d'une image qui est multi-spectrale ou hyper-spectrale, c'est-à-dire dans le cas 10 d'une image d'une scène constituée d'une pluralité d'images élémentaires de la scène acquises dans des bandes de longueurs d'ondes respectives différentes. Une fois une image acquise par le ou les instruments d'observation 11 du satellite 10, ladite image est mémorisée pour être transmise à une station sol 20, lorsque le satellite 10 survole ladite station sol 20. Une fois transmise à 15 une station sol 20, l'image est ensuite fournie à un dispositif de traitement (non représenté sur les figures) qui peut effectuer en temps différé les traitements visant à corriger dans l'image les effets de la présence d'aérosols et/ou de nuages dans l'atmosphère lors de l'acquisition de l'image, en mettant en oeuvre un procédé 50 de traitement d'image. 20 Le dispositif de traitement comporte par exemple un circuit de traitement comportant un ou plusieurs processeurs et des moyens de mémorisation (disque dur magnétique, mémoire électronique, disque optique, etc.) dans lesquels sont mémorisés des données et un produit programme d'ordinateur, sous la forme d'un ensemble d'instructions de code de 25 programme à exécuter pour mettre en oeuvre tout ou partie des étapes du procédé 50 de traitement d'image. Alternativement ou en complément, le circuit de traitement comporte un ou des circuits logiques programmables (FPGA, PLD, etc.), et/ou un ou des circuits intégrés spécialisés (ASIC), et/ou un ensemble de composants électroniques discrets, etc. adaptés à mettre en 30 oeuvre tout ou partie des étapes du procédé 50 de traitement d'image. En d'autres termes, le circuit de traitement correspond à un ensemble de moyens configurés de façon logicielle (produit programme d'ordinateur spécifique) et/ou matérielle (FPGA, PLD, ASIC, etc.) pour mettre en oeuvre les 3066301 10 différentes étapes du procédé 50 de traitement d'image. La figure 2 représente schématiquement les principales étapes d'un procédé 50 de traitement d'image visant à corriger les effets de la présence d'aérosols et/ou de nuages dans l'image acquise d'une scène. Tel qu'illustré 5 par la figure 2, le procédé 50 de traitement d'image comporte des étapes de : - 51 obtention d'un modèle de réflectance au sol et d'un modèle atmosphérique modélisant les effets de l'atmosphère entre la scène et le satellite 10, ledit modèle atmosphérique comportant au moins un modèle d'aérosols et un modèle de nuages modélisant 10 les effets de l'atmosphère dans des couches respectives d'altitudes respectives différentes, - 52 détermination de valeurs optimales de paramètres du modèle atmosphérique permettant d'optimiser ledit modèle atmosphérique par rapport à l'image, 15 - 53 détermination de valeurs de paramètres de correction pour chaque couche du modèle atmosphérique en fonction des valeurs optimales des paramètres du modèle atmosphérique, - 54 conversion de l'image en une image au sol avec effets d'environnement en fonction des valeurs des paramètres de 20 correction de chaque couche du modèle atmosphérique, - 55 conversion de l'image au sol avec effets d'environnement en image au sol sans effets d'environnement en fonction des valeurs des paramètres de correction de chaque couche du modèle atmosphérique. 25 Des exemples de mise en oeuvre des étapes ci-dessus sont décrits plus en détail ci-après, en référence notamment à la figure 3 qui représente un mode préféré de mise en oeuvre du procédé 50 de traitement d'image. Dans l'exemple non limitatif illustré par la figure 3, le procédé 50 de traitement d'image comporte en outre une étape 56 de conversion des valeurs 30 de pixels de l'image en valeurs en sommet d'atmosphère, dites valeurs de réflectance TOA (« Top of Atmosphere » dans la littérature anglo-saxonne). Une telle conversion vise donc à obtenir une image d'une valeur physique de réflectance telle qu'on peut la mesurer dans l'espace à l'extérieur de 3066301 11 l'enveloppe atmosphérique de la Terre. A) Obtention de modèle de réflectance au sol et de modèle atmosphérique Tel qu'indiqué ci-dessus, la procédé 50 de traitement d'image utilise 5 un modèle de scène qui comporte principalement deux types de modèles : - un modèle de réflectance au sol qui modélise la réflectance intrinsèque de la scène observée, et - un modèle atmosphérique qui modélise la fonction de transfert de l'atmosphère entre la scène et le satellite 10, et plus 10 particulièrement entre la scène et le sommet de l'atmosphère. A.1) Modèle de réflectance au sol Le modèle de réflectance au sol est par exemple basé sur un modèle SAIL/PROSPECT, qui modélise notamment la réflectance de scènes 15 correspondant à des couverts de végétation, la végétation étant le contributeur principal dans la majorité des scènes observées depuis un satellite 10. Le modèle SAIL/PROSPECT, également connu sous le nom de modèle PROSAIL dans la littérature scientifique, prend en compte, de manière conventionnelle, la direction d'illumination de la scène par le Soleil ainsi que la 20 direction de visée de l'instrument (modèle BRDF pour « Bidirectional Reflectance Distribution Function »). Le modèle SAIL/PROSPECT a été validé de manière exhaustive et est couramment utilisé par la communauté scientifique. On peut citer à titre d'exemple la publication scientifique : « PROSPECT + SAIL Models : A Review 25 of Use for Vegetation Characterization », Jacquemoud S., Verhoef W., Baret W., Bacour C., Zarco-Tejada P.J., Asner G.P., François C., Ustin, S.L., Remote Sensing of Environment 113, S56-S66. Il est également possible d'enrichir le modèle de réflectance au sol au moyen par exemple : 30 - d'une connaissance a priori du couvert observé et d'une gamme prédéterminée de densité de végétation (allant par exemple de forêt tropicale à région semi-aride), - d'une composante d'eau, qui met par exemple en oeuvre un 3066301 12 modèle de transfert radiatif utilisant le même formalisme que le modèle SAIL et les propriétés optiques de l'eau telles que caractérisées dans des eaux potentiellement peu profondes et turbides, dites de type CASE II (voir par exemple « Variations in the 5 Light Absorption Coefficients of Phytoplankton, Nonalgal Particles, and Dissolved Organic Matter in Coastal Waters Around Europe », Babin et al., Journal of Geophysical Research, 108, 3211) s'il est connu a priori qu'une forte présence d'eau est possible dans la scène observée (par exemple zone côtière, lacs, etc.), 10 - d'une réponse spectrale du sol prédéterminée, sélectionnée en fonction d'une connaissance a priori du type sol observé dans la scène lorsque la contribution du sol est susceptible d'être importante par rapport à la contribution du couvert végétal (montagne, désert, etc.), 15 - d'une modélisation d'autres caractéristiques susceptibles d'influencer la réflectance de la scène observée, comme par exemple une présence importante de zones brulées, de zones enneigées, de surfaces artificielles présentant une signature spectrale prédéterminée, etc. 20 A.2) Modèle atmosphérique Le modèle atmosphérique est par exemple basé sur un modèle LOWTRAN (voir par exemple « Users Guide to LOWTRAN 7 », Kneisys F.X. et al., 1988, Air Force Geophysics Lab Hanscom AFB MA). 25 Tel qu'indiqué ci-dessus, le modèle atmosphérique est un modèle multicouche comportant au moins deux couches d'altitudes respectives différentes, à savoir une couche modélisée par un modèle d'aérosols qui modélise les effets des aérosols et, au-dessus, une couche modélisée par un modèle de nuages qui modélise les effets des nuages. 30 Dans la suite de la description, on se place de manière non limitative dans le cas où le modèle atmosphérique comporte en outre un modèle de ciel modélisant les effets de l'atmosphère dans une couche d'altitude supérieure à la couche du modèle de nuages. 3066301 13 La répartition des différentes couches du modèle atmosphérique est représentée sur la figure 4. Tel qu'illustré par la figure 4, l'atmosphère est modélisée suivant trois couches principales : - une couche basse 30, qui correspond par exemple aux altitudes 5 comprises entre 0 et environ 2 kilomètres auxquelles se situent les aérosols, et qui est modélisée par le modèle d'aérosols, - une couche intermédiaire 31, qui correspond par exemple aux altitudes comprises entre environ 2 et environ 4 kilomètres auxquelles peuvent se situer des nuages (à l'exception de certains 10 types de nuages, comme par exemple des cirrus), et qui est modélisée par le modèle de nuages, - une couche haute 32, au-dessus de la couche intermédiaire 31 et jusqu'au sommet de l'atmosphère, qui introduit principalement une diffusion Rayleigh modélisée par le modèle de ciel. 15 Le modèle d'aérosols et le modèle de nuages sont par exemple modélisés en tant que couches de milieu trouble (« turbid medium » dans la littérature anglo-saxonne) avec une fonction de phase Henvey-Greenstein et un paramètre d'asymétrie adaptés aux comportements respectifs des aérosols et des nuages. Pour le modèle d'aérosols et le modèle de nuages, les fonctions 20 de transfert sont par exemple exprimées selon le formalisme 4-flux tel que développé dans le modèle SAIL. Les propriétés optiques de nuages sont bien connues et décrites dans la littérature, et peuvent être utilisées pour paramétrer le modèle de nuages, et pour établir une relation entre l'épaisseur optique d'un nuage et la réflectance 25 apparente dudit nuage, par exemple pour utiliser l'épaisseur optique de nuage comme paramètre du modèle de nuages. Pour un modèle d'aérosols de type LOWTRAN, le paramètre directeur est généralement la distance de visibilité, en relation avec la charge d'aérosols. Les propriétés optiques des aérosols peuvent être déduites de l'appel du 30 modèle LOWTRAN, par comparaison des résultats fournis par ledit modèle LOWTRAN en considérant d'une part une absence d'aérosols (visibilité maximale) et, d'autre part, un type particulier d'aérosols et une valeur de référence de la distance de visibilité. Ainsi, il est possible d'établir des relations 3066301 14 (« look-up tables ») entre le paramètre distance de visibilité du modèle LOWTRAN et l'épaisseur optique de la couche aérosols, et d'utiliser ladite distance de visibilité comme paramètre dudit modèle d'aérosols. Il est également possible d'enrichir le modèle d'aérosols et/ou le 5 modèle de nuages au moyen par exemple d'un modèle de sable pour les zones géographiques dans lesquelles des particules de sable en quantité non négligeable peuvent être retrouvées dans la couche basse 30 et/ou dans la couche intermédiaire 31. Le modèle de ciel modélise, de manière conventionnelle, la diffusion 10 Rayleigh introduite par la couche haute 32. L'intensité de la diffusion Rayleigh, en relation avec l'épaisseur optique de la colonne atmosphérique moléculaire, peut être contrôlée de manière conventionnelle par un paramètre équivalent d'altitude. Ce paramètre dépend du relief de la scène, varie modérément selon les conditions de pression et température, et correspond au-dessus des 15 nuages opaques à l'altitude apparente du plafond nuageux. Ce paramètre est désigné ci-après par « altitude plancher de la colonne atmosphérique ». A.3) Sélection du modèle de scène Afin de pouvoir traiter des images de scènes situés à des endroits 20 différents à la surface de la Terre, il est possible, dans des modes particuliers de mise en oeuvre, de mémoriser dans une base de données une bibliothèque de modèles de scène. Chaque modèle de scène mémorisé dans cette base de données correspond à une combinaison particulière d'un modèle de réflectance au sol et d'un modèle atmosphérique adaptée à une combinaison 25 particulière de type de paysage et de conditions climatiques. En classifiant différentes zones à la surface de la Terre en fonction de leurs types de paysage et de leurs conditions climatiques, il est possible d'établir une ou plusieurs cartes géographiques permettant de sélectionner, pour chaque zone à la surface de la Terre, le modèle de scène le plus adapté 30 au type de paysage et aux conditions climatiques rencontrés dans cette zone. Ainsi, la sélection revient principalement, dans le cas où une bibliothèque globale de modèles de scène a préalablement été constituée dans une base de données, à identifier la zone dans laquelle se trouve la scène 3066301 15 observée et à obtenir dans la base de données le modèle de scène associé à ladite zone. B) Détermination de valeurs optimales de paramètres du modèle de scène 5 Tel qu'indiqué ci-dessus, le procédé 50 de traitement d'image comporte une étape 52 de détermination de valeurs optimales de paramètres du modèle atmosphérique permettant d'optimiser ledit modèle atmosphérique par rapport à l'image. Ladite étape 52 de détermination de valeurs optimales de paramètres 10 du modèle atmosphérique s'inscrit dans un processus plus général d'inversion du modèle de scène qui conduit à déterminer également des valeurs optimales de paramètres du modèle de réflectance au sol permettant d'optimiser le modèle de réflectance au sol par rapport à l'image. De préférence, le modèle de réflectance au sol et le modèle atmosphérique du modèle de scène sont 15 optimisés conjointement et simultanément. En effet, un jeu particulier de valeurs des paramètres du modèle de scène permet d'obtenir une image théorique de ladite scène au sommet de l'atmosphère, qui peut être comparée à l'image acquise par le satellite 10. Par exemple, l'optimisation vise à maximiser la ressemblance entre 20 l'image théorique et l'image acquise par le satellite 10, c'est-à-dire qu'elle vise à déterminer les valeurs dites optimales desdits paramètres qui permettent d'obtenir une ressemblance maximale entre l'image théorique et ladite image acquise par le satellite 10. Toutefois, d'autres types d'optimisation peuvent être considérés et un processus d'optimisation comporte généralement la définition 25 préalable d'une fonction de coût à optimiser, c'est-à-dire à minimiser ou à maximiser suivant le type de fonction de coût. Le choix d'une fonction de coût particulière ne constitue qu'une variante d'implémentation de l'invention. Par exemple, les paramètres dont les valeurs optimales doivent être déterminées du modèle atmosphérique correspondent à: 30 - la distance de visibilité pour le modèle d'aérosols, - l'épaisseur optique de nuage pour le modèle de nuages, et - l'altitude plancher de la colonne atmosphérique pour le modèle de ciel. 3066301 16 Toutefois, l'invention est également applicable en considérant d'autres types de paramètres, et le choix d'un type particulier de paramètre ne constitue qu'une variante d'implémentation de l'invention. A titre d'exemple, il est également possible de considérer comme paramètre du modèle d'aérosols 5 l'épaisseur optique d'aérosols (AOD). Il est à noter que plus le nombre de mesures dans le domaine spectral est important, et plus les valeurs optimales déterminées desdits paramètres seront précises et proches de la réalité. A ce titre, l'utilisation d'une image multi-spectrale ou hyper-spectrale, constituée de plusieurs images 10 élémentaires représentant le rayonnement reçu d'une même scène dans des bandes de longueurs d'ondes respectives différentes permet d'augmenter le nombre de mesures et donc d'améliorer les performances de l'inversion du modèle atmosphérique (et, de manière plus générale, du modèle de scène). A l'issue de l'étape 52 de détermination de valeurs optimales du 15 modèle de scène, on dispose donc par exemple, en chaque pixel de l'image (et, dans le cas d'une image multi-spectrale ou hyper-spectrale, pour chaque image élémentaire) : - d'une valeur optimale pour la distance de visibilité, - d'une valeur optimale pour l'épaisseur optique de nuages, 20 - d'une valeur optimale pour l'altitude plancher de la colonne atmosphérique. En outre, les méthodes connues d'optimisation permettent également de fournir une valeur confiance associée à chaque valeur optimale déterminée d'un paramètre, représentative de la vraisemblance de ladite valeur optimale 25 compte tenu par exemple des valeurs respectives du pixel correspondant dans l'image théorique et dans l'image acquise par le satellite 10. Ainsi, à l'issue de l'étape 52 de détermination de valeurs optimales du modèle de scène, on peut également disposer, en chaque pixel de l'image (et, dans le cas d'une image multi-spectrale ou hyper-spectrale, pour chaque image élémentaire) : 30 - d'une valeur de confiance associée à la valeur optimale de la distance de visibilité, - d'une valeur de confiance associée à la valeur optimale de l'épaisseur optique de nuages, 3066301 17 - d'une valeur de confiance associée à l'altitude plancher de la colonne atmosphérique. Il est à noter que, pour la couche haute 32, la valeur optimale de l'altitude basse de la colonne atmosphérique est en principe sensiblement la 5 même pour tous les pixels de l'image. Par conséquent, il est possible, pour le modèle de ciel de la couche haute 32, de ne déterminer qu'une seule valeur optimale de l'altitude plancher de la colonne atmosphérique, ou de ne déterminer que quelques valeurs optimales pour quelques pixels de ladite l'image (et, éventuellement, les valeurs de confiance associées). 10 Il est également à noter que, de manière analogue, des valeurs optimales (et éventuellement des valeurs de confiance associées) peuvent être obtenues pour les paramètres du modèle de réflectance au sol. Dans des modes particuliers de mise en oeuvre, les valeurs de confiance associées aux valeurs optimales des paramètres du modèle 15 atmosphérique, calculées au cours de l'inversion du modèle de scène, peuvent être pondérées afin de tenir compte de différents critères. Suivant un premier exemple, il est possible de tenir compte de la valeur optimale de l'épaisseur optique de nuage pour pondérer la valeur de confiance de la valeur optimale de la distance de visibilité de la couche basse 20 30 du même pixel. En effet, dans le cas d'un nuage épais, la charge d'aérosols sous le nuage ne peut pas être déterminée correctement. Il est donc avantageux, dans le cas où la valeur optimale de l'épaisseur optique de nuage est importante, de réduire la valeur de confiance de la valeur optimale de la distance de visibilité de la couche basse 30 déterminée pour le même pixel. 25 Suivant un autre exemple, il est possible de tenir compte de valeurs optimales des paramètres du modèle de réflectance au sol pour pondérer les valeurs de confiance des valeurs optimales de la distance de visibilité de la couche basse 30. En effet, la détermination de la charge d'aérosols (et donc de la distance de visibilité) est plus précise au-dessus de zones sombres telles 30 que des forêts. Il est donc avantageux d'augmenter la valeur de confiance de la valeur optimale de la distance de visibilité de la couche basse 30 pour un pixel représentant une zone de végétation dense et sombre de la scène. A l'issue de l'étape 52 de détermination de valeurs optimales du 3066301 18 modèle de scène, on dispose donc de préférence, pour chaque paramètre et pour chaque couche du modèle atmosphérique (et, le cas échéant, pour chaque image élémentaire) : - d'une carte de valeurs optimales associées respectivement à 5 différents pixels, - d'une carte de valeurs de confiance associées auxdites valeurs optimales. Tel qu'indiqué précédemment, pour la couche haute 32, la valeur optimale de l'altitude plancher de la colonne atmosphérique est en principe 10 sensiblement la même pour tous les pixels de l'image. Par conséquent, il est possible de déterminer des cartes de valeurs optimales et des cartes valeurs de confiance associées uniquement pour la couche basse 30 (modèle d'aérosols) et pour la couche intermédiaire 31 (modèle de nuages). Dans des modes préférés de mise en oeuvre, et tel qu'illustré par la 15 figure 3, le procédé 50 de traitement d'image comporte en outre une étape 520 de filtrage spatial des cartes de valeurs optimales des paramètres. Avantageusement, le filtrage spatial d'une carte de valeurs optimales est effectué en fonction de la carte de valeurs de confiance associée. La figure 5 représente schématiquement un mode préféré de mise en 20 oeuvre, dans lequel l'étape 520 de filtrage spatial d'une carte de valeurs optimales comporte un calcul de cartes de valeurs optimales à résolution dégradée pour différentes valeurs de dégradation de résolution, en fonction de la carte de valeurs optimales et de la carte de valeurs de confiance associée, et un calcul d'une moyenne pondérée de la carte de valeurs optimale et des 25 cartes de valeurs optimales à résolution dégradée. En effet, dans l'exemple non limitatif illustré par la figure 5, diverses cartes de valeurs optimales (et de cartes de valeurs de confiance associées) à résolution dégradée sont calculées. Plus particulièrement, l'exemple illustré par la figure 5 représente 30 quatre dégradations de résolution successives par un facteur 2 (étapes Si, S2, S3 et S4 sur la figure 5), de sorte à obtenir des cartes de valeurs optimales (et de cartes de valeurs de confiance associées) à résolution dégradée par respectivement un facteur 2, 4, 8 et 16. Il est à noter qu'il est possible de 3066301 19 considérer un nombre de dégradations de résolutions successives par un facteur 2 différent de quatre, et l'exemple illustré par la figure 5 n'est pas limité à un nombre particulier de dégradations de résolution successives. Par exemple, il est possible d'effectuer des dégradations de résolution successives 5 par un facteur 2 jusqu'à obtenir une carte de valeurs optimales à résolution dégradée ne comportant qu'un seul pixel suivant une dimension. Chaque dégradation de résolution est de préférence réalisée au moyen de filtres de convolution rapides qui calculent une moyenne pondérée des valeurs optimales selon les valeurs de confiance associées à des sous- 10 résolutions spatiales successives qui correspondent à une dégradation de la résolution initiale par des facteurs 2 à 16. La carte de valeurs optimales (et la carte de valeurs de confiance associées) obtenue à l'issue de l'étape S4 est ensuite soumise à une augmentation de résolution par un facteur 2 (étape S5). Une moyenne 15 pondérée de la carte de valeurs optimales obtenue à l'issue de l'étape S3 et de la carte de valeurs optimales obtenue à l'issue de l'étape S5 est ensuite calculée (étape S6), au moyen de coefficients de pondération respectifs k et (1k), k étant par exemple choisi égal à 0.5. Les cartes de valeurs de confiance associées sont combinées selon la même moyenne pondérée. 20 La carte de valeurs optimales (et la carte de valeurs de confiance associées) obtenue à l'issue de l'étape S6 est ensuite soumise à une augmentation de résolution par un facteur 2 (étape S7). Une moyenne pondérée de la carte de valeurs optimales obtenue à l'issue de l'étape S2 et de la carte de valeurs optimales obtenue à l'issue de l'étape S7 est ensuite 25 calculée (étape S8), au moyen desdits coefficients de pondération k et (1-k). Les cartes de valeurs de confiance associées sont combinées selon la même moyenne pondérée. La carte de valeurs optimales (et la carte de valeurs de confiance associées) obtenue à l'issue de l'étape S8 est ensuite soumise à une 30 augmentation de résolution par un facteur 2 (étape S9). Une moyenne pondérée de la carte de valeurs optimales obtenue à l'issue de l'étape Si et de la carte de valeurs optimales obtenue à l'issue de l'étape S9 est ensuite calculée (étape S10), au moyen desdits coefficients de pondération k et (1-k). 3066301 20 Les cartes de valeurs de confiance associées sont combinées selon la même moyenne pondérée. La carte de valeurs optimales (et la carte de valeurs de confiance associées) obtenue à l'issue de l'étape S10 est ensuite soumise à une 5 augmentation de résolution par un facteur 2 (étape S11). Une moyenne pondérée de la carte de valeurs optimales initiale et de la carte de valeurs optimales obtenue à l'issue de l'étape S11 est ensuite calculée (étape S12), au moyen desdits coefficients de pondération k et (1-k). Les cartes de valeurs de confiance associées sont combinées selon la même moyenne pondérée. 10 La carte de valeurs optimales et la carte de valeurs de confiance associées obtenues à l'issue de l'étape 512 correspondent au résultat du filtrage spatial dans l'exemple illustré par la figure 5. Dans des modes préférés de mise en oeuvre, le filtrage spatial de la carte de valeurs optimales de l'épaisseur optique de nuage ne concerne pas 15 les pixels pour lesquels un critère prédéterminé d'opacité est vérifié. Par exemple, le critère d'opacité est vérifié si la valeur optimale de l'épaisseur optique de nuage est supérieure à une valeur seuil prédéterminée représentative de la présence d'un nuage opaque. En effet, le filtrage spatial, qui est un filtre passe-bas, a tendance à 20 adoucir les transitions rapides de la carte de valeurs optimales et permet donc de supprimer des artefacts liées par exemple à du bruit. Toutefois, un tel filtrage spatial contribuerait également à adoucir une délimitation d'un nuage épais, ce qui contribuerait à étaler cette délimitation sur des pixels voisins du nuage épais. C'est pourquoi il est avantageux de ne pas effectuer de filtrage 25 spatial sur les valeurs optimales correspondant à un nuage épais. Par exemple, il est possible d'identifier les pixels, dits « pixels masqués », pour lesquels le critère d'opacité est vérifié. Pour ces pixels masqués, la carte de valeurs de confiance initiale est par exemple modifiée pour forcer les valeurs de confiance associées aux pixels masqués à zéro. La carte de valeurs de 30 confiance ainsi modifiée est utilisée au cours de l'étape 520 de filtrage spatial. A l'issue de l'étape 520 de filtrage spatial, les valeurs optimales des pixels masqués et les valeurs de confiance associées sont par exemple remplacées par leurs valeurs initiales, avant filtrage spatial. 3066301 21 C) Détermination de paramètres de correction pour chaque couche Tel qu'indiqué ci-dessus, le procédé 50 de traitement d'image comporte une étape 53 de détermination de valeurs de paramètres de 5 correction pour chaque couche du modèle atmosphérique en fonction des valeurs optimales des paramètres du modèle atmosphérique. De manière générale, les paramètres de correction d'une couche du modèle atmosphérique correspondent à des paramètres décrivant la fonction de transfert de cette couche, et qui peuvent par conséquent être utilisés pour 10 inverser cette fonction de transfert et corriger les effets de la couche correspondante du modèle de transfert dans l'image. Dans des modes préférés de mise en oeuvre, les paramètres de correction considérés sont par exemple, pour chaque couche (couche basse 30, couche intermédiaire 31 et couche haute 32) : 15 - une réflectance intrinsèque pi de la couche considérée, - une atténuation intrinsèque ri de la couche considérée, qui correspond à l'atténuation introduite par la couche considérée sur le rayonnement reçu d'une couche immédiatement inférieure ou, dans le cas de la couche basse 30, reçu de la scène, 20 - une proportion d'effet environnement xenvi dans la transmission de la couche considérée. La proportion d'effet d'environnement xenvi correspond à la contribution de l'environnement du point observé dans le rayonnement transmis à la couche supérieure, et peut s'exprimer de manière conventionnelle 25 selon l'expression suivante : Tdo xenvi = Too + Tdo expression dans laquelle : - TOO représente la transmission directe de rayonnement par la couche considérée, c'est-à-dire le coefficient de transmission du rayonnement qui provient effectivement du point observé au sol 30 selon la ligne de visée considérée, par traversée (sans interception) de la couche considérée, 3066301 22 - -cd, représente la transmission indirecte de rayonnement par la couche considérée, c'est-à-dire le coefficient de transmission du rayonnement collecté selon la ligne de visée considérée qui provient de l'ensemble du voisinage du point observé au sol, par 5 diffusion au passage de la couche considérée. Rien n'exclut cependant de considérer d'autres types de paramètres de correction et le choix d'un type particulier de paramètre de correction ne constitue qu'une variante d'implémentation de l'invention. Les valeurs de tels paramètres de correction peuvent être calculées, 10 pour chaque couche, directement à partir du modèle atmosphérique en considérant les valeurs optimales déterminées pour les différents paramètres du modèle atmosphérique. L'étape 53 de détermination des valeurs des paramètres de correction et l'étape 52 de détermination des valeurs optimales peuvent donc être exécutées simultanément. 15 Ainsi, dans l'exemple précédent, l'étape 53 de détermination de valeurs de paramètres de correction doit permettre de déterminer, pour chaque pixel (et, le cas échéant, pour chaque image élémentaire, c'est-à-dire pour chaque bande de longueurs d'ondes considérée) : - trois valeurs de réflectance intrinsèque pi (une valeur par couche 20 du modèle atmosphérique), - trois valeurs d'atténuation intrinsèque Ti (une valeur par couche du modèle atmosphérique), - trois valeurs de proportion d'effet d'environnement xenvi (une valeur par couche du modèle atmosphérique). 25 Tel qu'indiqué précédemment, la valeur optimale de l'altitude plancher de la colonne atmosphérique du modèle de ciel est sensiblement la même pour tous les pixels de l'image, de sorte que les valeurs des paramètres de correction pour la couche haute 32 sont en principe également sensiblement les mêmes pour tous les pixels de l'image. 30 Afin de limiter la quantité de données à mémoriser, il est avantageux de mémoriser des tables de correction associées respectivement à chaque paramètre du modèle atmosphérique et à chaque bande de longueurs d'ondes considérée. Par exemple, cela revient à déterminer, pour chaque bande de 3066301 23 longueurs d'ondes considérée : - une table de correction selon la distance de visibilité, - une table de correction selon l'épaisseur optique de nuage, - une table de correction selon l'altitude plancher de la colonne 5 atmosphérique. Une telle table de correction comporte, pour une ou plusieurs valeurs typiques du paramètre de modèle considéré (distance de visibilité, épaisseur optique de nuage, altitude plancher de la colonne atmosphérique), les valeurs correspondantes des paramètres de correction (réflectance intrinsèque, 10 atténuation intrinsèque, proportion d'effet d'environnement). La table ci-dessous donne un exemple non limitatif d'une telle table de correction. Valeurs typiques du 3 7 10 L.- .1 200 paramètre de modèle Réflectance L.- .1 L.- .1 L.- .1 L.- .1 L.- .1 intrinsèque Pi Atténuation L...1 L...1 L...1 L...1 L...1 intrinsèque Ti Proportion d'effet [...] [...] [...] L.- .1 L.- .1 d'environnement xenvi Tel qu'indiqué précédemment, la valeur optimale de l'altitude plancher de la colonne atmosphérique est sensiblement la même pour tous les pixels de l'image, de sorte que la table de correction correspondante peut ne comporter 15 qu'une seule valeur typique égale à la valeur optimale. Ainsi, pour récupérer les valeurs des paramètres de correction à utiliser pour un pixel de l'image, par exemple pour la couche basse 30 du modèle atmosphérique, il suffit d'obtenir la valeur optimale de la distance de visibilité (pour la bande de longueurs d'ondes de l'image considérée). Cette 20 valeur optimale permet par exemple d'identifier, dans la table de correction associée à la distance de visibilité (et associée à la bande de longueurs d'ondes considérée), les valeurs typiques les plus proches de ladite valeur optimale déterminée ainsi que les valeurs des paramètres de correction associées, et de calculer les valeurs finales des paramètres de correction par 25 interpolation desdites valeurs des paramètres de correction fournies par la table de correction associée à la distance de visibilité. 3066301 24 D) Conversion de l'image en image au sol avec effets d'environnement Le procédé 50 de traitement d'image comporte également une étape 54 de conversion de l'image en image au sol avec effets d'environnement en 5 fonction de valeurs des paramètres de correction. L'étape 54 de conversion en image au sol avec effets d'environnement vise à corriger les effets liés à la présence d'aérosols et/ou de nuages, à l'exception des effets d'environnement. Plus particulièrement, l'étape 54 de conversion en image au sol avec effets d'environnement corrige successivement les effets de chaque couche du 10 modèle atmosphérique, en considérant lesdites couches successivement par altitudes décroissantes. En d'autres termes, les effets de la couche haute 32 sont corrigés en premier, les effets de la couche intermédiaire 31 sont corrigés en second, et les effets de la couche basse sont corrigés en dernier. Ainsi, pour un pixel donné, et en se plaçant après obtention de valeurs 15 de réflectance en sommet d'atmosphère, désignée ci-après par 0 , TOA, on obtient successivement : - une valeur de réflectance pi qui correspond à la réflectance au sommet de la couche intermédiaire 31, - une valeur de réflectance 102 qui correspond à la réflectance au 20 sommet de la couche basse 30, - une valeur de réflectance p3 qui correspond à la réflectance au sol avec effets d'environnement. Pour chaque couche du modèle atmosphérique, la correction s'effectue en fonction des valeurs des paramètres de correction, et vise 25 notamment à supprimer la réflectance intrinsèque de la couche considérée et à compenser l'atténuation intrinsèque de ladite couche considérée, par exemple selon l'expression suivante : x (PK - PL) expression dans laquelle OK 2 et Yi est un gain qui compense l'atténuation intrinsèque de la couche considérée (Yi = 1/T1). 30 Tel qu'indiqué précédemment, les valeurs des paramètres de correction pi et ri dépendent du paramètre de modèle considéré (distance de visibilité, épaisseur optique de nuage) et de la valeur optimale dudit paramètre 3066301 25 de modèle considéré au niveau du pixel considéré (et sont fonction de la bande de longueurs d'ondes considérée). Pour l'altitude plancher de la colonne atmosphérique, la valeur optimale est essentiellement la même quel que soit le pixel considéré, et la correction est essentiellement la même pour tous les 5 pixels de l'image (et fonction de la bande de longueurs d'ondes considérée). E) Conversion en image au sol sans effets d'environnement Le procédé 50 de traitement d'image comporte également une étape 55 de conversion de l'image au sol avec effets d'environnement en image au 10 sol sans effets d'environnement en fonction des valeurs des paramètres de correction. L'étape 55 de conversion en image au sol sans effets d'environnement vise donc à corriger les effets d'environnement. Plus particulièrement, l'étape 55 de conversion en image au sol sans effets d'environnement corrige successivement les effets de chaque couche du 15 modèle atmosphérique, en considérant lesdites couches successivement par altitudes décroissantes. En d'autres termes, les effets de la couche haute 32 sont corrigés en premier, les effets de la couche intermédiaire 31 sont corrigés en second, et les effets de la couche basse sont corrigés en dernier. Ainsi, pour un pixel donné, et en se plaçant après obtention de valeurs 20 de réflectance au sol avec effets d'environnement, désignée ci-après par p3,0 (et qui correspond à la réflectance p3 calculée ci-dessus), on obtient successivement : - une valeur de réflectance p3,1 qui correspond à la réflectance après correction de l'effet d'environnement de la couche haute 32, 25 - une valeur de réflectance 103,2 qui correspond à la réflectance après correction de l'effet d'environnement de la couche haute 32 et de la couche intermédiaire 31, - une valeur de réflectance 103,3 qui correspond à la réflectance au sol sans effets d'environnement. 30 Pour chaque couche du modèle atmosphérique, la correction s'effectue en fonction des valeurs des paramètres de correction, et vise à supprimer en tout ou partie les contributions de l'environnement, par exemple 3066301 26 selon l'expression suivante : fBackground ( , p3,K - xenvi x V--3,K ) P3,K+1 = 1 - xenvi expression dans laquelle 0 REVENDICATIONS1 - Procédé (50) de traitement d'une image optique acquise par un engin (10) spatial ou aérien, ladite image représentant un rayonnement dans au moins une bande de longueurs d'ondes reçu d'une scène à la surface de la Terre, caractérisé en ce qu'il comporte : - une obtention (51) d'un modèle de réflectance au sol et d'un modèle atmosphérique, ledit modèle atmosphérique comportant au moins un modèle d'aérosols et un modèle de nuages modélisant les effets de l'atmosphère dans des couches respectives d'altitudes respectives différentes, - une détermination (52) de valeurs optimales de paramètres du modèle atmosphérique permettant d'optimiser ledit modèle atmosphérique par rapport à l'image, - une détermination (53) de valeurs de paramètres de correction pour chaque couche du modèle atmosphérique en fonction des valeurs optimales des paramètres du modèle atmosphérique, - une conversion (54) de l'image en une image au sol avec effets d'environnement en fonction des valeurs des paramètres de correction de chaque couche du modèle atmosphérique, - une conversion (55) de l'image au sol avec effets d'environnement en image au sol sans effets d'environnement en fonction des valeurs des paramètres de correction de chaque couche du modèle atmosphérique. 2 - Procédé (50) selon la revendication 1, dans lequel les paramètres de correction comportent des réflectances intrinsèques et des atténuations intrinsèques pour chaque couche du modèle atmosphérique, et la conversion en image au sol avec effets d'environnement comporte une correction des réflectances intrinsèques et des atténuations intrinsèques des différentes couches du modèle atmosphérique. 3 - Procédé (50) selon l'une des revendications 1 à 2, dans lequel les paramètres de correction comportent des proportions d'effet d'environnement pour chaque couche du modèle atmosphérique, et la conversion en image au sol sans effets d'environnement comporte une 3066301 31 correction de contributions d'environnement introduites par les différentes couches en fonction des proportions d'effet d'environnement. 4 - Procédé (50) selon la revendication 3, dans lequel la correction d'une contribution d'environnement en un pixel pour une couche comporte : - un calcul d'une réflectance reçue dans un cône de vue dont le sommet correspond au pixel considéré, - une suppression de ladite réflectance reçue dans le cône de vue, pondérée par la proportion d'effet d'environnement pour la couche considérée et le pixel considéré. Procédé (50) selon l'une des revendications précédentes, comportant une identification de pixels de l'image pour lesquels un critère d'opacité de la couche du modèle de nuages est vérifié, dits « pixels masqués », et dans lequel la valeur de chaque pixel masqué est forcée (57), dans l'image au sol sans effets d'environnement, à une valeur prédéterminée. Procédé (50) selon la revendication 5, dans lequel la valeur prédéterminée correspond à la valeur du pixel masqué dans l'image après conversion en valeur en sommet de nuage ou en valeur en sommet d'atmosphère. Procédé (50) selon l'une des revendications précédentes, dans lequel la détermination (52) de valeurs optimales de paramètres du modèle atmosphérique comporte une détermination, pour chaque paramètre, d'une carte de valeurs optimales dudit paramètre associées respectivement à différents pixels de l'image, et un filtrage spatial (520) de ladite carte de valeurs optimales dudit paramètre. Procédé (50) selon la revendication 7, dans lequel, pour chaque paramètre, la détermination d'une carte de valeurs optimales dudit paramètre comporte en outre une détermination d'une carte de valeurs de confiance associées auxdites valeurs optimales, et le filtrage spatial (520) d'une carte de valeurs optimales est effectué en fonction de la carte de valeurs de confiance associée. 30 9 - Procédé (50) selon la revendication 8, dans lequel le filtrage spatial (520) d'une carte de valeurs optimales comporte un calcul de cartes de valeurs optimales à résolution dégradée pour différentes valeurs de dégradation de résolution, en fonction de la carte de valeurs optimales et de la carte de 5 10 5 - 156 - 7 - 20 8 - 25 3066301 32 valeurs de confiance associée, et un calcul d'une moyenne pondérée de la carte de valeurs optimale et des cartes de valeurs optimales à résolution dégradée. 10- Procédé (50) selon l'une des revendications précédentes, dans lequel 5 l'image correspond à une image brute dont les valeurs des pixels ont été converties en valeurs en sommet d'atmosphère. 11 - Procédé (50) selon l'une des revendications précédentes, dans lequel l'image correspond à une image haute résolution dont la résolution a été dégradée, comportant une détermination (581) de facteurs de correction à 10 appliquer à l'image pour obtenir l'image au sol sans effets d'environnement, et une conversion (582) de l'image haute résolution en une image haute résolution au sol sans effets d'environnement en fonction desdits facteurs de correction. 12 - Procédé (50) selon l'une des revendications précédentes, dans lequel le 15 modèle atmosphérique comporte en outre un modèle de ciel modélisant les effets de l'atmosphère dans une couche d'altitude supérieure à la couche du modèle de nuages. 13- Procédé (50) selon l'une des revendications précédentes, comportant une étape préalable de détermination si l'image comporte un nuage, un 20 modèle atmosphérique sans modèle de nuages étant utilisé lorsqu'il est déterminé que l'image ne comporte pas de nuage. 14- Produit programme d'ordinateur caractérisé en ce qu'il comporte un ensemble d'instructions de code de programme qui, lorsqu'elles sont exécutées par un processeur, configurent ledit processeur pour mettre en 25 oeuvre un procédé (50) de traitement d'image selon l'une des revendications précédentes. 15- Dispositif de traitement d'image acquise par un engin (10) spatial ou aérien, caractérisé en ce qu'il comporte des moyens configurés pour mettre en oeuvre un procédé (50) de traitement d'image selon l'une des 30 revendications 1 à 13.