Suite

Sortie du modèle de correction atmosphérique pour l'image Landsat 8 ?


Voici ma procédure, observation et question/demande :

Source de données : Earth Explorer Landsat 8 Grass Version SIG : 7.0.2 RC1 Windows 7 32 bits OS

J'ai essayé d'effectuer une correction atmosphérique en utilisant la référence du wiki Grass

https://grasswiki.osgeo.org/wiki/Atmospheric_correction

  • Étape 01 :

    J'ai ignoré l'utilisation de l'outil Toar et utilisé le calculateur de carte raster pour générer :

    une).Réflectance = REFLECTANCE_MULT_BAND_X(0.00002)*DN + REFLECTANCE_ADD_BAND_X (-0.1)

    b).Réflectance avec correction de l'angle solaire = Réflectance/ Sin(Ɵ)

    Ici Sin(Ɵ) = Sin(49.52284282) = 0.760664796730206

  • Étape 02 :

    Branché à i.atcorr les bandes de réflectance générées à l'étape 01 pour les bandes 4, 3, 2 et 8 avec leur fichier de contrôle respectif.

  • Étape 03 : Bande IHS Pan Sharpened 4,3,2 avec bande 8 générée à l'étape 02

  • Étape 04 : Génération d'un composite RVB en utilisant Rouge = bande 4, Vert = bande 3 et Bleu = bande 2. (Toutes les bandes dérivées à l'étape 03)

Je joins l'image résultante pour comparaison. Visiblement, j'ai raté quelque chose. L'image ajustée à l'atmosphère aurait dû être plus claire. Cela ne ressemble même pas à une image RVB.

Joindre également le fichier de métadonnées, les sorties de commande (afin que vous puissiez suivre mes étapes, pour remarquer s'il y avait quelque chose d'inhabituel) et un fichier de contrôle pour référence. J'ai utilisé le modèle Tropical Atmospheric et le modèle Urban Aerosol. Les valeurs utilisées dans le fichier de contrôle ont été collectées à partir du wiki de correction atmosphérique.

Quelqu'un peut-il me dire comment obtenir une image correcte ajustée à l'atmosphère ?

De plus, j'ai remarqué qu'il n'y avait plus d'extension i.landsat.dehaze disponible. Pourriez-vous s'il vous plaît me dire comment le désembuage peut être possible en l'absence de l'extension?

Fichier de métadonnées ================================================ ================================ GROUPE = L1_METADATA_FILE GROUP = METADATA_FILE_INFO ORIGINE = "Image avec l'aimable autorisation de l'US Geological Survey" REQUEST_ID = "0501510272804_00036" LANDSAT_SCENE_ID = "LC81370442015300LGN00" FILE_DATE = 2015-10-27T08:52:24Z STATION_ID = "LGN" PROCESSING_SOFTWARE_VERSION = "LPGS_2.5.1" END_GROUPTY = METADATA_FICHIER_INFO GROUPE DE DONNÉES = "LGN" PROCESSING_SOFTWARE_VERSION = "LPGS_2.5.1" END_GROUPTY = METADATA_SORTIE_FILE_INFO GROUPE DE DONNÉES = PRODUCTL2000 "GEOTIFF" SPACECRAFT_ID = "LANDSAT_8" SENSOR_ID = "OLI_TIRS" WRS_PATH = 137 WRS_ROW = 44 NADIR_OFFNADIR = "NADIR" TARGET_WRS_PATH = 137 TARGET_WRS_ROW = 44 DATE_ACQUIRED = 2015 ENTERUL24_TIME_27 SCENE_1950" ENTERUL24_TIME_PRO = 2015 ENTERUL24_TIME_CENE_1950 : 150 CORNER_UL_LON_PRODUCT = 89.24071 CORNER_UR_LAT_PRODUCT = 24.18946 CORNER_UR_LON_PRODUCT = 91.51825 CORNER_LL_LAT_PRODUCT = 22.01936 CORNER_LL_LON_PRODUCT = 89.29959 CORNER_LR_LAT_PRODUCT = 22.05469 CORNER_LR_LON_PRODUCT_ UL52 = 91.541 ON_X_PRODUCT = 117900.000 CORNER_UL_PROJECTION_Y_PRODUCT = 2676000.000 CORNER_UR_PROJECTION_X_PRODUCT = 349500.000 CORNER_UR_PROJECTION_Y_PRODUCT = 2676000.000 CORNER_LL_PROJECTION_X_PRODUCT = 117900.000 CORNER_LL_PROJECTION_Y_PRODUCT = 2439600.000 CORNER_LR_PROJECTION_X_PRODUCT = 349500.000 CORNER_LR_PROJECTION_Y_PRODUCT = 2439600.000 PANCHROMATIC_LINES = 15761 PANCHROMATIC_SAMPLES = 15441 REFLECTIVE_LINES = 7881 REFLECTIVE_SAMPLES = 7721 THERMAL_LINES = 7881 THERMAL_SAMPLES = 7721 FILE_NAME_BAND_1 = "LC81370442015300LGN00_B1.TIF" FILE_NAME_BAND_2 = "LC81370442015300LGN00_B2.TIF" FILE_NAME_BAND_3 = "LC81370442015300LGN00_B3.TIF" FILE_NAME_BAND_4 = "LC81370442015300LGN00_B4.TIF" FILE_NAME_BAND_5 = "LC81370442015300LGN00_B5.TIF" FILE_NAME_BAND_6 = "LC81370442015300LGN00_B6.TIF" FILE_NAME_BAND_7 = "LC81370442015300LGN00_B7.TIF" FILE_NAME_BAND_8 = "LC81370442015300LGN00_B8.TIF" FILE_NAME_BAND_9 = "LC81370442015300LGN00_B9.TIF" FILE_NAME_BAND_10 = "LC81370442015300LGN00_B10.TIF" FILE_NAME_ BAND_11 = "LC81370442015300LGN00_B11.TIF" FILE_NAME_BAND_QUALITY = "LC81370442015300LGN00_BQA.TIF" METADATA_FILE_NAME = "LC81370442015300LGN00_MTL.txt" BPF_NAME_OLI = "LO8BPF20151027040956_20151027054850.01" BPF_NAME_TIRS = "LT8BPF20151027040602_20151027054456.01" CPF_NAME = "L8CPF20151001_20151231.01" RLUT_FILE_NAME = "L8RLUT20150303_20431231v11.h5" END_GROUP = PRODUCT_METADATA GROUP = IMAGE_ATTRIBUTES cloud_cover = 4,67 CLOUD_COVER_LAND = 4,75 IMAGE_QUALITY_OLI = 9 IMAGE_QUALITY_TIRS = 9 TIRS_SSM_POSITION_STATUS = "nominale" ROLL_ANGLE = -0,001 SUN_AZIMUTH = 150,12621205 SUN_ELEVATION = 49,52284282 EARTH_SUN_DISTANCE = 0,9939432 GROUND_CONTROL_POINTS_VERSION = 2 GROUND_CONTROL_POINTS_MODEL = 298 GEOMETRIC_RMSE_MODEL = 6,823 = 4,812 GEOMETRIC_RMSE_MODEL_Y GEOMETRIC_RMSE_MODEL_X = 4,837 GROUND_CONTROL_POINTS_VERIFY = 88 GEOMETRIC_RMSE_VERIFY = 4.846 END_GROUP = IMAGE_ATTRIBUTES GROUP = MIN_MAX_RADIANCE RADIANCE_MAXIMUM_BAND_1 = 769.35406 RADIANCE_MINIMUM_BAND_1 = -63.53351 RADIANCE_MAXI MUM_BAND_2 = 787,82794 RADIANCE_MINIMUM_BAND_2 = -65,05909 RADIANCE_MAXIMUM_BAND_3 = 725,97668 RADIANCE_MINIMUM_BAND_3 = -59,95139 RADIANCE_MAXIMUM_BAND_4 = 612,18445 RADIANCE_MINIMUM_BAND_4 = -50,55439 RADIANCE_MAXIMUM_BAND_5 = 374,62637 RADIANCE_MINIMUM_BAND_5 = -30,93677 RADIANCE_MAXIMUM_BAND_6 = 93,16617 RADIANCE_MINIMUM_BAND_6 = -7,69369 RADIANCE_MAXIMUM_BAND_7 = 31,40199 RADIANCE_MINIMUM_BAND_7 = -2,59319 RADIANCE_MAXIMUM_BAND_8 = 692,82434 RADIANCE_MINIMUM_BAND_8 = -57,21366 RADIANCE_MAXIMUM_BAND_9 = 146,41248 RADIANCE_MINIMUM_BAND_9 = -12,09079 RADIANCE_MAXIMUM_BAND_10 = 22,00180 RADIANCE_MINIMUM_BAND_10 = 0,10033 RADIANCE_MAXIMUM_BAND_11 = 22,00180 RADIANCE_MINIMUM_BAND_11 = 0,10033 END_GROUP = MIN_MAX_RADIANCE GROUP = MIN_MAX_REFLECTANCE REFLECTANCE_MAXIMUM_BAND_1 = 1,210700 REFLECTANCE_MINIMUM_BAND_1 = -0,099980 REFLECTANCE_MAXIMUM_BAND_2 = 1,210700 REFLECTANCE_MINIMUM_BAND_2 = -0,099980 REFLECTANCE_MAXIMUM_BAND_3 = 1,210700 REFLECTANCE_MINIMUM_BAND_3 = -0,099980 RÉFLECTANCE _MAXIMUM_BAND_4 = 1,210700 REFLECTANCE_MINIMUM_BAND_4 = -0,099980 REFLECTANCE_MAXIMUM_BAND_5 = 1,210700 REFLECTANCE_MINIMUM_BAND_5 = -0,099980 REFLECTANCE_MAXIMUM_BAND_6 = 1,210700 REFLECTANCE_MINIMUM_BAND_6 = -0,099980 REFLECTANCE_MAXIMUM_BAND_7 = 1,210700 REFLECTANCE_MINIMUM_BAND_7 = -0,099980 REFLECTANCE_MAXIMUM_BAND_8 = 1,210700 REFLECTANCE_MINIMUM_BAND_8 = -0,099980 REFLECTANCE_MAXIMUM_BAND_9 = 1,210700 REFLECTANCE_MINIMUM_BAND_9 = -0,099980 END_GROUP = MIN_MAX_REFLECTANCE GROUP = MIN_MAX_PIXEL_VALUE QUANTIZE_CAL_MAX_BAND_1 = 65535 QUANTIZE_CAL_MIN_BAND_1 = 1 QUANTIZE_CAL_MAX_BAND_2 = 65535 QUANTIZE_CAL_MIN_BAND_2 = 1 QUANTIZE_CAL_MAX_BAND_3 = 65535 QUANTIZE_CAL_MIN_BAND_3 = 1 QUANTIZE_CAL_MAX_BAND_4 = 65535 QUANTIZE_CAL_MIN_BAND_4 = 1 QUANTIZE_CAL_MAX_BAND_5 = 65535 QUANTIZE_CAL_MIN_BAND_5 = 1 QUANTIZE_CAL_MAX_BAND_6 = 65535 QUANTIZE_CAL_MIN_BAND_6 = 1 QUANTIZE_CAL_MAX_BAND_7 = 65535 QUANTIZE_CAL_MIN_BAND_7 = 1 QUANTIZE_CAL_MAX_BAND_8 = 65535 QUANTIZE_CAL_MIN_BAND_8 = 1 QUANTIZE_ CAL_MAX_BAND_9 = 65535 QUANTIZE_CAL_MIN_BAND_9 = 1 QUANTIZE_CAL_MAX_BAND_10 = 65535 QUANTIZE_CAL_MIN_BAND_10 = 1 QUANTIZE_CAL_MAX_BAND_11 = 65535 QUANTIZE_CAL_MIN_BAND_11 = 1 END_GROUP = MIN_MAX_PIXEL_VALUE GROUP = RADIOMETRIC_RESCALING RADIANCE_MULT_BAND_1 = 1.2709E-02 RADIANCE_MULT_BAND_2 = 1.3014E-02 RADIANCE_MULT_BAND_3 = 1.1993E-02 RADIANCE_MULT_BAND_4 = 1.0113E-02 = RADIANCE_MULT_BAND_5 6.1886E-03 RADIANCE_MULT_BAND_6 = 1.5390E-03 RADIANCE_MULT_BAND_7 = 5.1874E-04 RADIANCE_MULT_BAND_8 = 1.1445E-02 RADIANCE_MULT_BAND_9 = 2.4186E-03 RADIANCE_MULT_BAND_10 = 3.3420E-04 RADIANCE_MULT_BAND_BAND_8 = 1.1445E-02 RADIANCE_MULT_BAND_9 = 2.4186E-03 RADIANCE_MULT_BAND_10 = 3.3420E-04 RADIANCE_3.511_MULT_BAND_AD210D = 3.3420E-03 RADIANCE_MULT_BAND_AD_042 = 3.3420E-04 RADIANCE_MULT420_ADBAND-04D_9 = 2.4186E-03 -59,96339 RADIANCE_ADD_BAND_4 = -50,56450 RADIANCE_ADD_BAND_5 = -30,94296 RADIANCE_ADD_BAND_6 = -7,69523 RADIANCE_ADD_BAND_7 = -2,59370 RADIANCE_ADD_BAND_8 = -57,22510 RADIANCE_ADD_BAND_9 = -12,09321 RADIANCE_ADD_BAND_10 = 0,10000 RADIANCE_ADD_BAND_11 = 0,10000 REFLECTANCE_MULT_BAND_1 = 2.0000E-05 RE FLECTANCE_MULT_BAND_2 = 2.0000E-05 REFLECTANCE_MULT_BAND_3 = 2.0000E-05 REFLECTANCE_MULT_BAND_4 = 2.0000E-05 REFLECTANCE_MULT_BAND_5 = 2.0000E-05 REFLECTANCE_MULT_BAND_6 = 2.0000E-05 REFLECTANCE_MULT_BAND_BAND_7 = 2.0000E-05 REFLECTANCE_MULT_BAND_6 = 2.0000E-05 REFLECTANCE_MULT_MULT_BAND-058E_7 -0,100000 REFLECTANCE_ADD_BAND_2 = -0,100000 REFLECTANCE_ADD_BAND_3 = -0,100000 REFLECTANCE_ADD_BAND_4 = -0,100000 REFLECTANCE_ADD_BAND_5 = -0,100000 REFLECTANCE_ADD_BAND_6 = -0,100000 REFLECTANCE_ADD_BAND_7 = -0,100000 REFLECTANCE_ADD_BAND_8 = -0,100000 REFLECTANCE_ADD_BAND_9 = -0,100000 END_GROUP = RADIOMETRIC_RESCALING GROUP = TIRS_THERMAL_CONSTANTS K1_CONSTANT_BAND_10 = 774,8853 K1_CONSTANT_BAND_11 = 480,8883 K2_CONSTANT_BAND_10 = 1321,0789 K2_CONSTANT_BAND_11 = 1201.1442 END_GROUP = TIRS_THERMAL_CONSTANTS GROUP = PROJECTION_PARAMETERS MAP_PROJECTION = "UTM" DATUM = "WGS84" ELLIPSOID = "WGS84" UTM_ZONE = 46 GRID_CELL_SIZE_PANCHROMATIC = 15,00 GRID_CELL_SIZE_REFL = 30.00 GRID_CELL_SIZE_THERMAL = 30.00 ORIENTATION = "NORTH_UP" RESAMPLING_OPTION = "CUBIC_CONVOLUTION" END_GROUP = PROJECTION_PARAMETERS END_GROUP = L1_METADATA_FILE END Journal de sortie de commande =============================== ================================================== ===== Bande de correction atmosphérique 4 (ven 13 novembre 23:54:00 2015) i.atcorr -r --verbose [email protected] parameters=E:LandCover_ClassificationAeronetAtcorr_B4.txt output=ATM_B4 * * ******************************* 6s version 4.2b ******************* ************* * * identité des conditions géométriques * * ------------------------------ - * * Observation Landsat 8 * * * * mois : 10 jour : 27 * * angle zénithal solaire : 40,69° angle azimutal solaire : 150,19° * * angle zénithal de vue : 0,00° angle azimutal de vue : 0,00 deg * * angle de diffusion : 139,31 deg différence d'angle azimutal : 150,19 deg * * * * description du modèle atmosphérique * * ----------------------------- * * identité du modèle atmosphérique : * * tropical (uh2o=4.12g/cm2,uo3 =.247cm-atm ) * * * * identité type aérosols : * * modèle aérosols urbains * * * * identité optique condition : * * visibilité : 7.54 km opt. épais. 550nm : 0.5500 * * * * condition spectrale * * ------------------ * * Landsat 8 Red * * valeur de la fonction filtre : * * wl inf= 0.626 mic wl sup = 0,682 mic * * * * type de cible * * ----------- * * sol homogène * * réflectance constante sur les spectres 0,000 * * * * description de l'élévation de la cible * * ------- --------------------- * * pression au sol [mb] 1012,54 * * altitude au sol [km] 0,004 * * teneur en gaz au niveau cible : * * uh2o= 4,110 g/cm2 uo3= 0,248 cm-atm * * * * correction atmosphérique activée * * -------------------------------- * Correction atmosphérique… Correction atmosphérique terminée. (Ven 13 novembre 23:56:08 2015) Commande terminée (2 min 8 sec) Bande de correction atmosphérique 3 (Sam 14 novembre 00:04:55 2015) i.atcorr -r --verbose [email protected] parameters=E :LandCover_ClassificationAeronetAtcorr_B3.txt output=ATM_B3 * ****************************** 6s version 4.2b *** *************************** * * identité des conditions géométriques * * ---------------- --------------- * * Observation Landsat 8 * * * * mois : 10 jour : 27 * * angle zénithal solaire : 40,69 degrés angle azimutal solaire : 150,19 degrés * * angle zénithal de vue : 0.00 deg vue angle azimutal: 0.00 deg * * angle de diffusion: 139.31 deg différence d'angle azimutal: 150.19 deg * * * * description du modèle atmosphérique * * -------------------- --------- * * identité du modèle atmosphérique : * * tropical (uh2o=4.12g/cm2,uo3 =.247cm-atm) * * * * identité du type d'aérosols : * * Modèle d'aérosols urbains * * * * état optique identité : * * visibilité : 7.54 km opt. épais. 550nm : 0.5500 * * * * condition spectrale * * ------------------ * * Landsat 8 Green * * valeur de la fonction filtre : * * wl inf= 0.513 mic wl sup = 0,600 mic * * * * type de cible * * ----------- * * sol homogène * * réflectance constante sur les spectres 0,000 * * * * description de l'élévation de la cible * * ------- --------------------- * * pression au sol [mb] 1012,54 * * altitude au sol [km] 0,004 * * teneur en gaz au niveau cible : * * uh2o= 4,110 g/cm2 uo3= 0,248 cm-atm * * * * correction atmosphérique activée * * -------------------------------- * Correction atmosphérique… Correction atmosphérique terminée. (Sam 14 novembre 00:06:28 2015) Commande terminée (1 min 32 sec) Bande de correction atmosphérique 2 (Sam 14 novembre 00:12:59 2015) i.atcorr -r --verbose [email protected] parameters=E :LandCover_ClassificationAeronetAtcorr_B2.txt output=ATM_B2 * ****************************** 6s version 4.2b *** *************************** * * identité des conditions géométriques * * ---------------- --------------- * * Observation Landsat 8 * * * * mois : 10 jour : 27 * * angle zénithal solaire : 40,69 degrés angle azimutal solaire : 150,19 degrés * * angle zénithal de vue : 0.00 deg vue angle azimutal: 0.00 deg * * angle de diffusion: 139.31 deg différence d'angle azimutal: 150.19 deg * * * * description du modèle atmosphérique * * -------------------- --------- * * identité du modèle atmosphérique : * * tropical (uh2o=4.12g/cm2,uo3 =.247cm-atm) * * * * identité du type d'aérosols : * * Modèle d'aérosols urbains * * * * état optique identité : * * visibilité : 7.54 km opt. épais. 550nm : 0.5500 * * * * condition spectrale * * ------------------ * * Landsat 8 Blue * * valeur de la fonction filtre : * * wl inf= 0.436 mic wl sup = 0,527 mic * * * * type de cible * * ----------- * * sol homogène * * réflectance constante sur les spectres 0,000 * * * * description de l'élévation de la cible * * ------- --------------------- * * pression au sol [mb] 1012,54 * * altitude au sol [km] 0,004 * * teneur en gaz au niveau cible : * * uh2o= 4,110 g/cm2 uo3= 0,248 cm-atm * * * * correction atmosphérique activée * * -------------------------------- * Correction atmosphérique… Correction atmosphérique terminée. (Sam 14 novembre 00:14:32 2015) Commande terminée (1 min 32 sec) Bande de correction atmosphérique 8 (Sam 14 novembre 00:17:58 2015) i.atcorr -r --verbose [email protected] parameters=E :LandCover_ClassificationAeronetAtcorr_B8.txt output=ATM_B8 * ****************************** 6s version 4.2b *** *************************** * * identité des conditions géométriques * * ---------------- --------------- * * Observation Landsat 8 * * * * mois : 10 jour : 27 * * angle zénithal solaire : 40,69 degrés angle azimutal solaire : 150,19 degrés * * angle zénithal de vue : 0.00 deg vue angle azimutal: 0.00 deg * * angle de diffusion: 139.31 deg différence d'angle azimutal: 150.19 deg * * * * description du modèle atmosphérique * * -------------------- --------- * * identité du modèle atmosphérique : * * tropical (uh2o=4.12g/cm2,uo3 =.247cm-atm) * * * * identité du type d'aérosols : * * Modèle d'aérosols urbains * * * * état optique identité : * * visibilité : 7.54 km opt. épais. 550nm : 0.5500 * * * * condition spectrale * * ------------------ * * Landsat 8 Panchromatic * * valeur de la fonction filtre : * * wl inf= 0.488 mic wl sup = 0,692 mic * * * * type de cible * * ----------- * * sol homogène * * réflectance constante sur les spectres 0,000 * * * * description de l'élévation de la cible * * ------- --------------------- * * pression au sol [mb] 1012,54 * * altitude au sol [km] 0,004 * * teneur en gaz au niveau cible : * * uh2o= 4,110 g/cm2 uo3= 0,248 cm-atm * * * * correction atmosphérique activée * * -------------------------------- * Correction atmosphérique… Correction atmosphérique terminée. (Sam 14 nov 00:19:53 2015) Commande terminée (1 min 54 sec) Pan Sharpening (Sam 14 nov 00:23:58 2015) i.pansharpen --verbose [email protected] [email protected] blue= [email protected] [email protected] output=Méthode PAN_RGB=ihs Exécution d'une netteté panoramique avec une image panoramique haute résolution : 15.000000 Utilisation de l'algorithme IHS<->RGB Transformation en espace colorimétrique IHS… Mise en correspondance de l'histogramme… Retransformation en espace colorimétrique RVB et accentuation… Attribution tables de couleurs égalisées en gris pour les images de sortie… Les cartes de sortie pan-nettes suivantes ont été générées : PAN_RGB_red PAN_RGB_green PAN_RGB_blue Pour visualiser la sortie, exécutez : g.region -p raster=PAN_RGB_red d.rgb r=PAN_RGB_red g=PAN_RGB_green b=PAN_RGB_blue Si désiré , combinez les canaux en une seule carte RVB avec « r.composite ». Les couleurs des canaux peuvent être rééquilibrées à l'aide de i.colors.enhance. Grouper  n'existe pas encore. Création… Ajout de cartes raster au groupe … Ajout d'une carte raster  au groupe Ajout d'une carte raster  au groupe Ajout d'une carte raster  au groupe (Sam 14 nov 00:26:31 2015) Commande terminée (2 min 32 sec) Composite (Sam 14 nov 00:31:31 2015) r.composite --verbose [email protected] [email protected] blue [email protected] output=RGB_PAN.tif Création d'une table de couleurs pour la carte raster en sortie… Ecriture de la carte raster … r.composite complet. Carte raster  établi. (Sam 14 novembre 00:32:04 2015) Commande terminée (33 sec) Exemple de fichier de contrôle pour la bande 2 ============================ ================================================== === 18 10 27 4,41 90,373 23,113 1 3 0 0,55 -0,004 -1000 116

Je suggérerais de faire un pas à la fois.

Essayez d'abord le module i.landsat.toar et comparez les valeurs avec vos résultats à votre étape 1.

Deuxièmement, vous devez être prudent avec le module i.pansharpen car je pense qu'il n'accepte que les données 8 bits, donc les données flottantes produites à partir de la correction atmosphérique vont vous poser des problèmes.

Si vous souhaitez affiner les données corrigées atmosphériquement, essayez d'utiliser mapcalc (assurez-vous que les paramètres de votre région ont une résolution de 15 m x 15 m). Les formules sont écrites dans la doc i.pansharpen et disponibles ici : http://www.sensorsportal.com/HTML/DIGEST/august_2014/Vol_177/P_RP_0183.pdf

Vérifiez donc vos bandes individuelles avec une palette de couleurs en niveaux de gris, puis si elles ont toutes l'air bien, essayez le composite.


Voir la vidéo: Landsat 8 Image Preprocessing in QGIS. Atmospheric Correction. LULC and Change Detection Series (Octobre 2021).