Problema con el cálculo de QGIS Zonal Statistics MEAN

Question

Necesito calcular algunas estadísticas zonales a partir de archivos raster para algunos polígonos en un shapefile. Estoy utilizando el siguiente código que funciona:

    eligibleLayer = QgsVectorLayer("/home/usr/Desktop/polygons.shp", "polygons", "ogr")
    for raster in findRasters(windMapsDirectory, '*.tif'):
        (infilepath, infilename) = os.path.split(raster) 
        windMapName = infilename
        list1 = re.split('\W+',windMapName)
        list2 = re.split("_", list1[0])
        list3 = list2[2]
        string1 = ''.join(list3)
        if string1 == 'K':
            zoneStat = QgsZonalStatistics(eligibleLayer, windMapName, 'K_', 1)
            zoneStat.calculateStatistics(None)

El problema es que sólo necesito la media de las estadísticas zonales (sin recuento ni suma). He intentado utilizar:

        if string1 == 'K':
            zoneStat = QgsZonalStatistics(eligibleLayer, windMapName, 'K_', 1, QgsZonalStatistics.Mean)

Sin embargo, obtengo el siguiente error: el objeto de tipo 'QgsZonalStatistics' no tiene el atributo 'Mean' . Estoy utilizando QGIS 2.8.6.

Answer 1

Por alguna razón su código no está produciendo un QgsZonalStatistics objeto. Imprima 'zoneStats' para corroborarlo. Se debe observar algo similar a:

<qgis._analysis.QgsZonalStatistics object at 0x9fc316ec>

He probado el siguiente código:

from qgis.analysis import QgsZonalStatistics

mapcanvas = iface.mapCanvas()

layers = mapcanvas.layers()
provider = layers[1].dataProvider()

path = provider.dataSourceUri()

zoneStat = QgsZonalStatistics(layers[0], path,"", 1, QgsZonalStatistics.Mean)

print zoneStat

zoneStat.calculateStatistics(None)

para esta situación:

y, tras ejecutar el script, se ha añadido sólo el valor medio; como se puede observar en la siguiente imagen:

A mí me funciona perfectamente.

Answer 2

Un poco tarde, pero acabo de tener el mismo problema usando QGIS 2.14 ( type object 'QgsZonalStatistics' has no attribute 'Mean' , con un QgsZonalStatistics objeto).

Resulta que funciona bien con QGIS 2.18, así que probablemente hubo algún cambio en la API de QGIS que no está documentado, o un error que fue corregido.