encode shapes

[vincentsarago]

ref #1

[vincentsarago]

FYI, I tried to use Mapbox's wagyu to create valid geometry but it didn't really change the results

from wagyu import enums as wEnums
from wagyu.linear_ring import LinearRing as wLinearRing
from wagyu.point import Point as wPoint
from wagyu.polygon import Polygon as wPolygon
from wagyu.wagyu import Wagyu

        w = Wagyu()
        wext = wPolygon(
            [
                wLinearRing(
                    [wPoint(x, y) for x, y in exterior.coords]
                )
            ]
        )
        w.add_polygon(wext, wEnums.PolygonKind.SUBJECT)
        for interior in list(interiors):
            wint = wPolygon(
                [
                    wLinearRing(
                        [wPoint(x, y) for x, y in interior.coords]
                    )
                ]
            )
            w.add_polygon(wint, wEnums.PolygonKind.SUBJECT)

        result = w.execute(
            wEnums.OperationKind.UNION,
            wEnums.FillKind.EVEN_ODD,
            wEnums.FillKind.EVEN_ODD
        )
        for p in result.polygons:
            for l in p.linear_rings:
                feature.add_ring(len(l))
                for p in l.points:
                    feature.set_point(int(p.x), int(p.y))

[vincentsarago]

Here we use shapely to split the geometry between interior/exterior

[vincentsarago]

seems that applying a 0 buffer creates a valid geometry (for GEOS)

[RichardScottOZ]

Read that the other day too, when I had a large polygonize output

[vincentsarago]

so this is still WIP. (see #1)

I've deployed a tiler using https://github.com/developmentseed/titiler-mvt but I'm getting a lot of warnings because the geometry are still invalid.

I may try to use pygeos instead of shapely (and maybe try wagyu again)