Discussion :

[Juliendu92]

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

[Juliendu92]

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

Pour ceux que ça intéresse, j'ai trouvé comment faire :
J'ai utilisé data[i,j] qui me donne la valeur de mon pixel pour la ligne i et la colonne j.
J'ai fait une boucle, j'ai rempli ma liste que j'ai enfin transformé en tableau.