problem avec GetrasterBand sous python 3

Version imprimable

Bonjour,

je sais pas si quelqu'un peut m'aider sur ce sujet mais j'essaie de travailler sur des raster de type MODIS format TIFF

mon problème est que je n'arrive pas à utiliser bien le package GDAL sur Python (peut-être)

j'essaie de lire la première bande dans une liste de raster TIFF mais je n'y arrive pas

voila mon code et l'erreur :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
import os
    import glob
    from osgeo import gdal
    import numpy as np
    import matplotlib.pyplot as plt
    import scipy.ndimage
    import pandas as pd
    import datetime as dt
 
 
    inDir = r'D:\Python\MODIS\Doy2000/'
    os.chdir(inDir)  
 
    outDir = os.path.normpath(os.path.split(inDir)[0] + os.sep + 'output') + '\\' 
    if not os.path.exists(outDir): os.makedirs(outDir)
 
    # Create list of MODIS rasters for MODIS EVI
 
    modis_EVI = glob.glob('MOD13Q1.006__250m_16_days_EVI_**.tif')
 
    # Read file 
    EVI = gdal.Open(modis_EVI [0])
 
EVI_meta = EVI.GetMetadata()  
 
    geotransform = EVI.GetGeoTransform()
    proj= EVI.GetProjection() 
 
    #Read the band 
    EVIBand = EVI.GetRasterBand(1)   
 
    # Import band as an array           
    EVIData = EVIBand.ReadAsArray().astype('float')

l'erreur :

AttributeError: 'NoneType' object has no attribute 'GetRasterBand'
 'NoneType' object has no attribute 'GetGeoTransform' 
'NoneType' object has no attribute 'GetProjection'

je sais pas quoi faire

merci

02/01/2019, 08h19
BufferBob
salut,

~~j'ai regardé rapidement, au minimum il y a un espace en trop ligne 22 :~~ (edit: c'est très faux, wiztricks donne la solution plus bas)
Code:

1 2 EVI = gdal.Open(modis_EVI [0]) # bad EVI = gdal.Open(modis_EVI[0]) # good
02/01/2019, 09h34
wiztricks

Salut,

Citation:

Envoyé par loula05

l'error :

AttributeError: 'NoneType' object has no attribute 'GetRasterBand'

Le NoneType' object correspond à la variable EVI assignée par gdal.Open(modis_EVI [0]). Et si vous lisez un peu la documentation de la bibliothèque osgeo, None est ce que retourne gdal.open en cas de problème pour ouvrir ou traiter le fichier.

- W
02/01/2019, 20h10
loula05

Citation:

Envoyé par wiztricks

Salut,

Le NoneType' object correspond à la variable EVI assignée par gdal.Open(modis_EVI [0]). Et si vous lisez un peu la documentation de la bibliothèque osgeo, None est ce que retourne gdal.open en cas de problème pour ouvrir ou traiter le fichier.

- W

Bonjour,

merci je comprends sinon je vais chercher de documentation
02/01/2019, 20h12
loula05
Citation:
Envoyé par BufferBob

salut,

j'ai regardé rapidement, au minimum il y a un espace en trop ligne 22 :

Code:

1 2 EVI = gdal.Open(modis_EVI [0]) # bad EVI = gdal.Open(modis_EVI[0]) # good
Bonjour,

est ce que l'espace peu provoquer un problème et avoir ce type des erreurs ?

merci
02/01/2019, 20h29
wiztricks
Citation:

Envoyé par loula05

est ce que l'espace peu provoquer un problème et avoir ce type des erreurs ?

Non. Pour vous en convaincre, vous lancez la console Python et vous essayez:
Code:

1 2 3 4 liste = [ 'toto.dat' ] open(liste [0], 'w') open(liste[0], 'w') open(liste [ 0 ], 'w')
- W
03/01/2019, 15h13
loula05

Bonjour,

moi j'ai un problème au niveau de getrasterband c'est pas au niveau de gdal.open

je comprend pas c'est quoi le probleme

voilà le problème :

Code:

1
2
3
4
5
6
7
8
9
evi = gdal.Open(modis_EVI[0])
 
Band = evi.GetRasterBand(1)
Traceback (most recent call last):
 
  File "<ipython-input-38-09828cf19778>", line 1, in <module>
    Band = evi.GetRasterBand(1)
 
AttributeError: 'NoneType' object has no attribute 'GetRasterBand'

03/01/2019, 15h32
BufferBob

Citation:

Envoyé par wiztricks

Non. Pour vous en convaincre, vous lancez la console Python et (...)

oops, effectivement, merci d'avoir corrigé

Citation:

Envoyé par wiztricks

Le NoneType' object correspond à la variable EVI assignée par gdal.Open(modis_EVI [0]). (...) None est ce que retourne gdal.open en cas de problème pour ouvrir ou traiter le fichier.

ça semble assez clair, c'est de gdal.Open() que provient le None, est-ce que modis_EVI[0] contient bien un nom de fichier ? est-ce que ce fichier existe bel et bien ? est-ce qu'on a le droit de lire dedans etc.

Citation:

Envoyé par BufferBob

oops, effectivement, merci d'avoir corrigé

ça semble assez clair, c'est de gdal.Open() que provient le None, est-ce que modis_EVI[0] contient bien un nom de fichier ? est-ce que ce fichier existe bel et bien ? est-ce qu'on a le droit de lire dedans etc.

quand j'ecris :
Code:

1 2 3 4 5 print(modis_EVI) ['MOD13Q1.006__250m_16_days_EVI_doy2000209_aid0001.tif', 'MOD13Q1.006__250m_16_days_EVI_doy2000225_aid0001.tif', 'MOD13Q1.006__250m_16_days_EVI_doy2000241_aid0001.tif', 'MOD13Q1.006__250m_16_days_EVI_doy2000257_aid0001.tif', 'MOD13Q1.006__250m_16_days_EVI_doy2000273_aid0001.tif', 'MOD13Q1.006__250m_16_days_EVI_doy2000289_aid0001.tif', 'MOD13Q1.006__250m_16_days_EVI_doy2000305_aid0001.tif', 'MOD13Q1.006__250m_16_days_EVI_doy2000321_aid0001.tif', 'MOD13Q1.006__250m_16_days_EVI_doy2000337_aid0001.tif', 'MOD13Q1.006__250m_16_days_EVI_doy2000353_aid0001.tif'] print(modis_EVI[0]) MOD13Q1.006__250m_16_days_EVI_doy2000209_aid0001.tif
il me donne le nom du fichier , mais quand j essai de lire la bande avec getrasterband s affiche l erreur

03/01/2019, 15h44
loula05
aussi quand je mis print(evi) me donne None
Code:

1 2 print(evi) None
03/01/2019, 16h03
wiztricks
Citation:

Envoyé par loula05

moi j'ai un problème au niveau de getrasterband c'est pas au niveau de gdal.open

Vous avez ce problème:
Code:

1 2 3 4 5 6 >>> evi = None >>> evi.getrasterband Traceback (most recent call last): File "<stdin>", line 1, in <module> AttributeError: 'NoneType' object has no attribute 'getrasterband' >>>
à la différence que dans votre code, c'est gdal.open qui retourne le None assigné à la variable "evi".

- W
03/01/2019, 16h12
loula05

Oui veuillez m'excuser c'est ça je viens de comprendre le problème c'est gdal.open quand j'écris prin(evi) me donne None, je sais pas comment je peux contrôler ce problème afin de continuer le travail je bloque.

:(
03/01/2019, 16h17
wiztricks

Citation:

Envoyé par loula05

je sais pas comment je peux controler ce problème afin de continuer le travail je bloque .

Vous ouvrez la documentation, et vous y cherchez comment lui afficher un message d'erreur plus clair (que simplement None) pour avoir une piste sur ce qui ne va pas.
Vous pouvez aussi essayer de remplacer gdal.open par un open "standard" pour vous assurer déjà que le fichier est accessible.
Impossible de trouver une solution sans comprendre la cause du problème et comme ces fichiers sont sur vos disques durs...

- W
03/01/2019, 16h33
loula05

Citation:

Envoyé par wiztricks

Vous ouvrez la documentation, et vous y cherchez comment lui afficher un message d'erreur plus clair (que simplement None) pour avoir une piste sur ce qui ne va pas.
Vous pouvez aussi essayer de remplacer gdal.open par un open "standard" pour vous assurer déjà que le fichier est accessible.
Impossible de trouver une solution sans comprendre la cause du problème et comme ces fichiers sont sur vos disques durs...

- W

mes fichiers sont des raster format TIF j'ai deja travailler sur ces données avec autre programme et ils posent aucun problème mais pour python je suis débutante j 'essai de travailler mais j'arrive pas a comprendre ce problème meme si j'ai lu quelque cours !

merci
03/01/2019, 17h32
loula05
j'ai remarqué et j'ai lu sur plusieurs site qu'il faut ajouter ces lignes apres gdal.Open !

est ce que c'est pour fixer l erreur
Code:

1 2 3 4 src_ds = gdal.Open( "INPUT.tif" ) if src_ds is None: print 'Unable to open INPUT.tif' sys.exit(1)
merci d'avance
03/01/2019, 17h54
wiztricks

Salut,

Citation:

Envoyé par loula05

est ce que c'est pour fixer l erreur

Cela vous donnera juste un message d'erreur plus explicite qu'AttributeError... (ce qui n'est pas si mal).

- W
03/01/2019, 18h56
loula05

Citation:

Envoyé par wiztricks

Salut,

Cela vous donnera juste un message d'erreur plus explicite qu'AttributeError... (ce qui n'est pas si mal).

- W

est ce que si j'ajoute le chemin complet de ma liste des rasters ici va eliminer l'erreur ?

Code:

evi = gdal.Open('c:/Doy2000/modis_EVI[0]')

sinon je comprends pas vraiment comment résoudre le problème
03/01/2019, 19h45
wiztricks

Citation:

Envoyé par loula05

est ce que si j'ajoute le chemin complet de ma liste des rasters ici va eliminer l'erreur ?

Code:

evi = gdal.Open('c:/Doy2000/modis_EVI[0]')

Si c'est un nom de chemin complet, il devrait se terminer par ".tif".

Citation:

sinon je comprends pas vraiment comment résoudre le problème

Si vous voulez comprendre, il va falloir prendre le temps d'apprendre à programmer, et pas juste recopier des lignes de code "pour voir". Et çà, on ne peut pas le faire à votre place.

- W
03/01/2019, 20h05
loula05

Citation:

Envoyé par wiztricks

Si c'est un nom de chemin complet, il devrait se terminer par ".tif".

Si vous voulez comprendre, il va falloir prendre le temps d'apprendre à programmer, et pas juste recopier des lignes de code "pour voir". Et çà, on ne peut pas le faire à votre place.

- W

je pose des questions ici c'est normal et je pense pas que quelqu'un sait programmer dès le début , en tout cas je suis dans un forum pour apprendre , si vous voulez pas aidez pas besoin de dire au gens vous copiez des codes .

merci
07/01/2019, 22h33
loula05

Bonjour ,

je vous remercie pour vos réponses , problème résolu c'etait un problème au niveau format de données ( extension TIF) j'ai du renomer ou re enregistrer mes données en format TIF encore 2em fois afin qu'il puisse gdal les lire et trouver

merci à tout