Discussion :

[loula05]

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]

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

[loula05]

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]

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]

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]

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'

[wiztricks]

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

[wiztricks]

Salut,

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]

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