Télédétection Réflexion de surface Landsat et Albédo

Je suis assez verte en ce qui concerne RS et je fais ce projet d'école avec lequel j'ai besoin d'aide. Le projet consiste à mesurer les changements d'albédo pour la langue d'un glacier et j'utilise la programmation R pour le faire. J'ai donc obtenu mes données de l'explorateur de la Terre et j'ai opté pour le produit de réflexion de surface afin d'économiser quelques étapes. Cependant, les valeurs de ces cellules sont plus élevées que prévu, autour de 2000 selon la bande. Je dis plus haut parce que lorsque j'utilise une formule moyenne pondérée pour calculer l'albédo, que j'ai trouvée dans certains articles cités ci-dessous, le résultat me donne des valeurs de cellules d'environ 1000, pas le 0 à 1 que j'attendrais et l'albédo. Le document mentionne qu'il utilise également une réflectance de surface pour cette formule.

Ma question est la suivante: cette formule est-elle mauvaise? Albédo = 0,493Bande2 + 0,203 Bande + 0,150 Bande5 + 0,154Bande

Ou faut-il effectuer un autre traitement sur mon produit Earth Explorer Surface Reflectance. https://landsat.usgs.gov/landsat-surface-reflectance-high-level-data-products

J'utilise également des images Landsat 4-7.

sources:

Pimentel, Rafael et al. "Comparaison entre l'albédo des neiges obtenu à partir des images Landsat TM, ETM + et le produit d'albédo SPOT VEGETATION dans un site montagneux méditerranéen." Hydrologie 3.1 (2016): 10.

En outre, c'est vraiment similaire mais utilise TOA au lieu de la réflectance de surface: http://yceo.yale.edu/how-convert-landsat-dns-albedo

Les données de réflectance Landsat que vous avez téléchargées de l'USGS ont été mises à l'échelle en utilisant un facteur d'échelle de 0,0001. Donc, multiplier le nombre numérique par 0,0001 vous donnera une valeur comprise entre 0 et 1. La valeur de 2000 dont vous parlez, je crois, est la valeur saturée (qui devrait être 20000). Voir ce document pour plus de détails: https://landsat.usgs.gov/sites/default/files/documents/ledaps_product_guide.pdf

La conversion des données de réflectance de surface Landsat en albédo est en réalité beaucoup plus compliquée que l'application d'une seule équation. L'équation que vous avez énumérée n'est qu'une approximation rapide, qui devrait vous donner une idée générale de la façon dont l'albédo a changé. Donc, si vous utilisez simplement cela pour une comparaison simple, je dirais que c'est assez bon, tant que vous utilisez la réflectance de surface (au lieu de TOA). Si cela vous intéresse vraiment et que vous souhaitez en savoir plus, vous pouvez consulter ce document: http://www.sciencedirect.com/science/article/pii/S0034425711001477