Erreur lors de l'utilisation d'OGR pour déterminer si un centre de pixels est à l'intérieur d'un polygone

Question

J'essaye de développer (en utilisant C++ - MSVS 12.0) une fonction qui permet de déterminer quels pixels (d'une image raster) ont leur centre dans un polygone (précédemment peuplé en utilisant un fichier de formes). J'utilise GDAL 1.11.0 (juste installé, en utilisant devinstall) en construisant à partir de la source et en utilisant l'option INCLUDE_OGR_FRMTS = YES. Je peux utiliser GDAL et la plupart des fonctions OGR sans problème. Cependant, lorsque j'utilise le code suivant:

if (polygon->Contains(tmpPoint))

Je reçois le message d'erreur: ERREUR 6: GEOS supportent pas activé

Tout le monde sait comment résoudre ce problème?

J'utilise:

#include "ogrsf_frmts.h"

et ma fonction est déclarée:

void FindPixels(GDALDataset *image, OGRLayer *poLayer, OGRPolygon *polygon)

et une partie de mon code est:

OGRPoint *tmpPoint = NULL

OGRSpatialReference *spatialReference = NULL;

spatialReference = polygon->getSpatialReference();

tmpPoint = new OGRPoint();

tmpPoint->assignSpatialReference(spatialReference);

boucle

commencent:

tmpPoint->setX(imgTLX + (j * imgRes) + imgResHalf);

tmpPoint->setY(imgTLY - (i * imgRes) - imgResHalf);

if (polygon->Contains(tmpPoint))

Merci d'avance!

MB

Answer 1

Utilisez GDALRasterizeLayers pour graver l'image du polygone sur une trame. De cette façon, vous trouverez tous les pixels qui tombent dans le polygone, ou non. La valeur par défaut est de graver le pixel uniquement si le centre intersecte un polygone. Si la couche source comporte plusieurs polygones, vous devrez peut-être les distinguer en définissant un filtre d'attribut ou en utilisant des champs d'ID d'attribut gravés (bien que cela ne fonctionne pas si les polygones se chevauchent).