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