Bonjour,

J'ai ouvert mon fichier .hdf avec GDAL, je tente d'extraire les données de mon fichier dans une zone géographique délimitée par la ligne I à Ifin ainsi que par la colonne J à J fin

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
               if not os.path.isfile(fichier):
                    print("Le fichier n'existe pas !")
               else:
                    gdal_dataset = gdal.Open(fichier, GA_ReadOnly)
                    gdal_dataset.GetSubDatasets()
                    driver = gdal.GetDriverByName('HFA')
                    driver.Register()
                    cols , rows = gdal_dataset.RasterXSize, gdal_dataset.RasterYSize
                    print cols, rows
                    for rows in range (I,Ifin):
                        for cols in range (J,Jfin):
                            data = gdal_dataset.ReadAsArray(0, 0, cols, rows)
Je cherche en fait l'équivalent très simple de Matlab qui lit mon fichier et extrait les données entre les lignes I et Ifin et les lignes J et Jfin :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
datatout=hdfread(fichier,'PIXEL DATA');
data=datatout(I:Ifin, J:Jfin);