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