banniere

Le portail francophone de la géomatique


Toujours pas inscrit ? Mot de passe oublié ?
Nom d'utilisateur    Mot de passe              Toujours pas inscrit ?   Mot de passe oublié ?

Annonce

Rencontres QGIS 2025

L'appel à participation est ouvert jusqu'au 19 janvier 2025!

Pages: 1

 

#1 Wed 25 July 2012 09:34

phil008
Participant occasionnel
Date d'inscription: 10 Sep 2011
Messages: 23

Traiter une image MODIS (sous GDAL) pour faire une composition colorée

Bonjour,

Je n'arrive pas à traiter des images MODIS, pour obtenir une image en couleurs naturelles ('true colors'). D'après ce que j'ai compris en fouillant sur Google il faut télécharger une scene MOD02HKM, puis combiner les canaux 4, 3 et 1. Une fois telechargée je regarde sous GDAL les infos du fichier pour savoir quels canaux extraire:

 

Code:

gdalinfo MOD02HKM.A2011205.0905.005.2011205192823.hdf

[...]

SUBDATASET_1_NAME=HDF4_EOS:EOS_SWATH:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":MODIS_SWATH_Type_L1B:EV_500_RefSB
  SUBDATASET_1_DESC=[5x4060x2708] EV_500_RefSB MODIS_SWATH_Type_L1B (16-bit unsigned integer)
  SUBDATASET_2_NAME=HDF4_EOS:EOS_SWATH:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":MODIS_SWATH_Type_L1B:EV_500_RefSB_Uncert_Indexes
  SUBDATASET_2_DESC=[5x4060x2708] EV_500_RefSB_Uncert_Indexes MODIS_SWATH_Type_L1B (8-bit unsigned integer)
  SUBDATASET_3_NAME=HDF4_EOS:EOS_SWATH:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":MODIS_SWATH_Type_L1B:EV_250_Aggr500_RefSB
  SUBDATASET_3_DESC=[2x4060x2708] EV_250_Aggr500_RefSB MODIS_SWATH_Type_L1B (16-bit unsigned integer)
  SUBDATASET_4_NAME=HDF4_EOS:EOS_SWATH:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":MODIS_SWATH_Type_L1B:EV_250_Aggr500_RefSB_Uncert_Indexes
  SUBDATASET_4_DESC=[2x4060x2708] EV_250_Aggr500_RefSB_Uncert_Indexes MODIS_SWATH_Type_L1B (8-bit unsigned integer)
  SUBDATASET_5_NAME=HDF4_EOS:EOS_SWATH:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":MODIS_SWATH_Type_L1B:EV_250_Aggr500_RefSB_Samples_Used
  SUBDATASET_5_DESC=[2x4060x2708] EV_250_Aggr500_RefSB_Samples_Used MODIS_SWATH_Type_L1B (8-bit integer)
  SUBDATASET_6_NAME=HDF4_SDS:UNKNOWN:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":0
  SUBDATASET_6_DESC=[2030x1354] Latitude (32-bit floating-point)
  SUBDATASET_7_NAME=HDF4_SDS:UNKNOWN:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":1
  SUBDATASET_7_DESC=[2030x1354] Longitude (32-bit floating-point)
  SUBDATASET_8_NAME=HDF4_SDS:UNKNOWN:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":2
  SUBDATASET_8_DESC=[5x4060x2708] EV_500_RefSB (16-bit unsigned integer)
  SUBDATASET_9_NAME=HDF4_SDS:UNKNOWN:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":3
  SUBDATASET_9_DESC=[5x4060x2708] EV_500_RefSB_Uncert_Indexes (8-bit unsigned integer)
  SUBDATASET_10_NAME=HDF4_SDS:UNKNOWN:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":4
  SUBDATASET_10_DESC=[2x4060x2708] EV_250_Aggr500_RefSB (16-bit unsigned integer)
  SUBDATASET_11_NAME=HDF4_SDS:UNKNOWN:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":5
  SUBDATASET_11_DESC=[2x4060x2708] EV_250_Aggr500_RefSB_Uncert_Indexes (8-bit unsigned integer)
  SUBDATASET_12_NAME=HDF4_SDS:UNKNOWN:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":6
  SUBDATASET_12_DESC=[2x4060x2708] EV_250_Aggr500_RefSB_Samples_Used (8-bit integer)
  SUBDATASET_13_NAME=HDF4_SDS:UNKNOWN:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":7
  SUBDATASET_13_DESC=[16x10] Noise in Thermal Detectors (8-bit unsigned integer)
  SUBDATASET_14_NAME=HDF4_SDS:UNKNOWN:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":8
  SUBDATASET_14_DESC=[16x10] Change in relative responses of thermal detectors (8-bit unsigned integer)
  SUBDATASET_15_NAME=HDF4_SDS:UNKNOWN:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":9
  SUBDATASET_15_DESC=[203x16x10] DC Restore Change for Thermal Bands (8-bit integer)
  SUBDATASET_16_NAME=HDF4_SDS:UNKNOWN:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":10
  SUBDATASET_16_DESC=[203x2x40] DC Restore Change for Reflective 250m Bands (8-bit integer)
  SUBDATASET_17_NAME=HDF4_SDS:UNKNOWN:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":11
  SUBDATASET_17_DESC=[203x5x20] DC Restore Change for Reflective 500m Bands (8-bit integer)
  SUBDATASET_18_NAME=HDF4_SDS:UNKNOWN:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":12
  SUBDATASET_18_DESC=[203x15x10] DC Restore Change for Reflective 1km Bands (8-bit integer)

Quand j'extrais et combine sous GDAL les canaux 4, 3 et 1, (HDF4_SDS:UNKNOWN:"MOD02HKM.A2011205.0905.005.2011205192823.hdf":4 etc...) je n'obtiens pas ma composition colorée mais un truc avec des parasites. Même quand je visionne les canaux separement sous QGIS, c'est plein de parasite. Même resultat lorsque j'essaye avec d'autres images. Est ce que je me serais trompé de canaux. Je suis perdu, c'est vraiment pas simple de traiter des images MODIS. Comment faites vous pour obtenir une image en vraie couleurs?

Dernière modification par phil008 (Wed 25 July 2012 09:40)

Hors ligne

 

 

Pages: 1

Pied de page des forums

Powered by FluxBB