Ouvrir un fichier HDF avec Python

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

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:

---

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

Citation:

---

Envoyé par 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:

---

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()

Oui, il faut aussi l'importer

Code:

---

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é.

Citation:

---

Envoyé par VinsS

Oui, il faut aussi l'importer

Code:

---

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.