Discussion :

[Juliendu92]

Bonsoir,
Je dois écrire un code python qui lit un fichier HDF.
J'ai trouver le module GDAL.

J'ai tenté ça :
from osgeo import gdal
gdal_dataset = gdal.Open ("/data/geospatial_10/ucfajlg/MOD12/MCD12Q1.A2005001.h17v03.005.2008310174635.hdf")
gdal_dataset.GetSubDatasets()

mais ça fonctionne pas, le message m'indique " AttributeError: 'NoneType' object has no attribute 'GetSubDatasets'"
Quelqu'un pourrait-il me débloquer ???

Julien

[VinsS]

Salut,

Le message d'erreur signifie que gdal.Open() a retourné None.

Normalement, il y a une erreur là derrière mais, d'après la doc [0], elle n'est pas accessible en python.

Alors première chose, vérifie ton fichier:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
 
from osgeo import gdal
import os
fname = "/data/geospatial_10/ucfajlg/MOD12/MCD12Q1.A2005001.h17v03.005.2008310174635.hdf"
if not os.path.isfile(fname):
    print("Le fichier n'existe pas !")
else:
    gdal_dataset = gdal.Open(fname, GA_ReadOnly)
    gdal_dataset.GetSubDatasets()

[0] http://www.gdal.org/gdal_tutorial.html

[Juliendu92]

Salut,

Le message d'erreur signifie que gdal.Open() a retourné None.

Normalement, il y a une erreur là derrière mais, d'après la doc [0], elle n'est pas accessible en python.

Alors première chose, vérifie ton fichier:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
 
from osgeo import gdal
import os
fname = "/data/geospatial_10/ucfajlg/MOD12/MCD12Q1.A2005001.h17v03.005.2008310174635.hdf"
if not os.path.isfile(fname):
    print("Le fichier n'existe pas !")
else:
    gdal_dataset = gdal.Open(fname, GA_ReadOnly)
    gdal_dataset.GetSubDatasets()

[0] http://www.gdal.org/gdal_tutorial.html

Salut, merci de prendre du temps pour me débloquer.

le fichier HDF est la bande du vert d'une image SPOT, ça ne fonctionne pas...
J'ai tapé le code ci-dessous et j'ai le message suivant : NameError: name 'GA_ReadOnly' is not defined
Il semblerait qu'il ne comprenne pas 'GA_ReadOnly'

from osgeo import gdal
import os

fname = "C:/0001_B0.HDF"
if not os.path.isfile(fname):
print("Le fichier n'existe pas !")
else:
gdal_dataset = gdal.Open(fname, GA_ReadOnly)
gdal_dataset.GetSubDatasets()

[VinsS]

Oui, il faut aussi l'importer

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
from gdalconst import GA_ReadOnly

Ce que tu aurais pu voir si tu avais jeté un oeil sur la doc que je t'ai indiqué.

[Juliendu92]

Oui, il faut aussi l'importer

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
from gdalconst import GA_ReadOnly

Ce que tu aurais pu voir si tu avais jeté un oeil sur la doc que je t'ai indiqué.

J'ai passé la soirée à la lire cette doc quand j'ai tapé GDAL sur Google,
j'ai essayé mais je n'y arrivais pas, je débute cette année avec la programmation, je suis un peu noyé par toutes les infos, tous les modules, les distributions de Python avec package... Merci de ton aide, ça fonctionne.