ST_Polygonize/Shapely cannot polygonize when overlapping with shared node/point

Question

I'm trying to use both polygonize_full (Shapely, python) and ST_Polygonize (postgis) to get all polygons from a set of line strings.

There is one case where it fails: when there is an "overlapping polygon" with a shared "node".

Pictured here:

The result set is empty, but shapely provides a 'dangles' that provides the entire "failed" linestring,

Has anyone encountered this? Is there a way to gracefully solve this? I couldn't find anything in the docs.

I would expect these two be two separate polygons that overlap each other, despite their shared node.

Here's the WKT I'm using for example.

MULTILINESTRING((-122.43682090001732 37.767652910517015,-122.4366733785215 37.767652910517015),(-122.4366733785215 37.767652910517015,-122.4366733785215 37.76759990327466),(-122.4366733785215 37.76759990327466,-122.43682090001732 37.76759990327466),(-122.43682090001732 37.76759990327466,-122.43682090001732 37.767652910517015),(-122.43682090001732 37.767652910517015,-122.43675920920998 37.767717579352684),(-122.43675920920998 37.767717579352684,-122.43675920920998 37.767626406895836),(-122.43675920920998 37.767626406895836,-122.43682090001732 37.767652910517015))

Answer 1

You may try using shapely.ops.unary_union() before polygonize your multistring.

e.g.:

GoodMultistring= shapely.ops.unary_union(BadMultistring)
polygons = shapely.ops.polygonize(GoodMultistring)