problem avec GetrasterBand sous python 3

**loula05** · 02/01/2019, 01h08

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

**BufferBob** · 02/01/2019, 08h19

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

**wiztricks** · 02/01/2019, 09h34

Salut,

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

**loula05** · 02/01/2019, 20h10

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

**loula05** · 02/01/2019, 20h12

Envoyé par BufferBob

salut,

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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

**wiztricks** · 02/01/2019, 20h29

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
liste = [ 'toto.dat' ]
open(liste [0], 'w')
open(liste[0], 'w')
open(liste [ 0 ], 'w')

- W

**loula05** · 03/01/2019, 15h13

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

**loula05** · 03/01/2019, 15h24

voilà le problème :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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'

**BufferBob** · 03/01/2019, 15h32

Envoyé par wiztricks

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

oops, effectivement, merci d'avoir corrigé

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.

**loula05** · 03/01/2019, 15h38

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

**loula05** · 03/01/2019, 15h44

aussi quand je mis print(evi) me donne None

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
print(evi)
None

**wiztricks** · 03/01/2019, 16h03

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

**loula05** · 03/01/2019, 16h12

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.

**wiztricks** · 03/01/2019, 16h17

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

**loula05** · 03/01/2019, 16h33

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

**loula05** · 03/01/2019, 17h32

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

**wiztricks** · 03/01/2019, 17h54

Salut,

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

**loula05** · 03/01/2019, 18h56

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

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

**wiztricks** · 03/01/2019, 19h45

Envoyé par loula05

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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

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

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

**loula05** · 03/01/2019, 20h05

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

problem avec GetrasterBand sous python 3 [Python 3.X]

Python

Discussions similaires

Partager

Partager