AssertionFailedException: found two shells in EdgeRing list

[hmy930303]

Please note

Expected behavior and actual behavior.

Expected unary_union successfully

Steps to reproduce the problem.

it happens when i use shapely union operation (unary_union()), all of all geoms is_valid is True. But sometime when i use set operation such as result = list(set(geoms_list)) , seems to make random order , the unary_union(result) runs successfully.

Code

import shapely
from shapely.ops import unary_union

def test_unary_union():
    # copy all datas from upload text file to here
    datas_string = "copy to here"
    geoms_data = [shapely.wkt.loads(geom_wkt) for geom_wkt in datas_string.split("\t")]
    # all geoms are valid
    assert all([geom.is_valid for geom in geoms_data])
    # this operation will raise error, i don't know why
    union = unary_union(geoms_data)
    # after set operation, the unary_union runs successfully (maybe random the order of geoms, the count of geoms is the same)
    new_geoms = list(set(geoms_data))
    union = unary_union(new_geoms)

Data

data.txt

Operating system

Ubuntu 22.04

Shapely version and provenance

2.0.2 with pip install GEOS ersion: 3.11.1

[theroggy]

I cannot reproduce the problem using geos 3.12.1, so most likely it has been solved.

You could try with shapely 2.0.3, which includes geos 3.11.3. geos 3.11.3 contains most of the same bugfixes as geos 3.12.1.