Le portail francophone de la géomatique

Bonjour,

Votre requête n'est pas très claire :
- vous souhaitez itérer sur tous les polygones de votre shapefile si j'ai bien compris ?
- vous voulez des stats agrégées pour chaque polygone ou une liste des valeurs de pixel dans chacun ? Savez-vous qu'il y a une fonction, dans QGIS, de statistiques zonales qui permet de récupérer les valeurs agrégées d'un raster le long de polygones ?

- provider est un nom plutôt saugrenu pour une QgsVectorLayer. vlayer conviendrait mieux
- extent est un objet d'étendue qui n'est pas itérable, d'où le problème au niveau de "in extent"
- for devrait plutôt être associé aux polygones, donc qqch du type :
for feat in layer.getFeatures() :
(croisement feat*raster)

Bonjour,

Sur le principe, le code pourrait ressembler à ceci :

iface = qgis.utils.iface

layers = iface.legendInterface().layers()

# on repere la couche raster et la couche vecteur dans la legende (en considerant que seulement deux couches présentes dedans)
rlayer = filter(lambda x : isinstance(x, QgsRasterLayer), layers)[0]
vlayer = filter(lambda x : isinstance(x, QgsVectorLayer), layers)[0]

def createLayerWithFeature(feat) :
    
    memLayer = QgsVectorLayer("Polygon", "masque", "memory")
    memLayer.dataProvider().addFeatures([feat])
    
    return memLayer
    
def returnListOfRasterValuesByPolygonFeature(vlayer, rlayer) :

    l = list()
    for feat in vlayer.getFeatures() :
        
        # on croise le polygone avec le raster
        memLayer = createLayerWithFeature(feat)
        ret = processing.runalg("saga:clipgridwithpolygon", rlayer, memLayer, 0, None)
        rLayerMasked = processing.getObject(ret["OUTPUT"])
        
        # on lit ses valeurs
        rows = res.height()
        cols = res.width()

        block = res.dataProvider().block(1, res.extent(), cols, rows) # en assumant que c'est la premiere bande que l on souhaite extraire
        l2 = list()
        for i in range(rows) :
            for j in range(cols) :
                l2.append(block.value(i,j))
            
        l.append(l2)
        
    return l

# execution
returnListOfRasterValuesByPolygonFeature(vlayer, rlayer)

Bonjour,

je vous remercie de votre réponse.

1)
mais si j'ai pas bien compris la ligne de code suivante, je pense que l'algo clipgridwithpolygon prend en paramètre : le raster d'entrée, le polygone qui sert de mask et l'output ? Pourquoi le 0 et None ?

ret = processing.runalg("saga:clipgridwithpolygon", rlayer, memLayer, 0, None).

2)
Sachant que mon objectif est de récupérer les pixels qui intersectent avec le polygone et de leur appliquer une transparence, est ce que je dois nécessairement passer par un clip.?

Merci d'avance.

Bonjour,

"2) vous souhaitez appliquer une transparence à des pixels ? C'est différent de ce que vous demandiez au départ...En gros, vous souhaitez faire un trou (masque) dans le raster à l'endroit du polygone ?"

Je suis toujours dans la même logique afin je crois... je récupère les pixels qui intersectent avec mon shapefile (polygone) dans une liste ensuite je leur applique un transparent "pour faire un trou (masque) dans le raster à l'endroit du polygone" exactement comme vous le dites.

C'est exactement ça.

Là c'est une première étape... la suite consistera à partir d'un buffer sur des objets ponctuels et linéaires pour créer le trou.

Merci beaucoup de votre "réactivité"

Bonjour,

je crois que récupérer les pixels est plus pertinent dans mon cas.

Je cherche à créer un outil semblable au plugin Mapswipelayer de qgis. Le trou que je vais créer va me permettre de visualiser une partie de l'image qui est en arrière plan.

Bonjour,

D'accord.

Merci beaucoup