2

如何在不将其转换为 WKT 字符串的情况下迭代 GEOS 多边形中的点？目前我们正在做这样的事情：

GEOSGeometry *geom = GEOSGeomFromWKT("POLYGON ((1 1, 2 1, 2 2, 1 2, 11))");
char geomAsWKT[900] = GEOSGeomToWKT(geom);
/* Iterate over the geomAsWKT to get the points */

我尝试过的所有功能（GEOSGeomGetPointN、GEOSGeom_getCoordSeq和其他一些功能）仅适用于LinearRing.

4

1 回答 1

3

要遍历多边形的点，您必须LinearRing通过调用GEOSGetExteriorRing来获取，如下所示。

此示例适用于MultiPolygons 或Polygons。此示例打印构成多面体的 LinearRings 中点的所有 x,y 坐标。如果 inputGeom 是多边形，它也可以工作。

GEOSGeometry *inputGeom = GEOSGeomFromWKT("MULTIPOLYGON (((30 20, 45 40, 10 40, 30 20)), ((15 5, 40 10, 10 20, 5 10, 15 5)))");
const GEOSGeometry *linearRing;
const GEOSCoordSequence *coordSeq;
int numGeom = GEOSGetNumGeometries(inputGeom);
int n;
for (n=0; n < numGeom; n++) {
    linearRing = GEOSGetExteriorRing(GEOSGetGeometryN(inputGeom, n));
    printf("%s\n", GEOSGeomToWKT(linearRing));
    unsigned int numPoints, p;
    numPoints = GEOSGeomGetNumPoints(linearRing);
    for (p=0; p < numPoints; p++) {
        double xCoord, yCoord;
        coordSeq = GEOSGeom_getCoordSeq(linearRing);
        GEOSCoordSeq_getX(coordSeq, p, &xCoord);
        GEOSCoordSeq_getY(coordSeq, p, &yCoord);

        printf("%f, %f\n", xCoord, yCoord);
    }
}

就我的目的而言，多边形中从来没有孔，所以这就是我所需要的。

于 2014-06-03T18:55:52.540 回答