Issues with overlay - oversimplifies in some cases, not sure what's wrong

Two geodataframes with boundaries (census tracts and proprietary boundaries):

after splits = gp.overlay(nhds, cts, how='intersection') result is weired:

[ozak]

@Casyfill I think that is what the geodataframe function is_valid does. Usually something like df.geometry = df.geometry.buffer(0) solves simple issues that cause geometries to not be valid.

[jorisvandenbossche]

well, I will dig into the gdal issues then, I guess

You could try to create a new clean environment? (that’s what also solved it in the issue you linked to). I suppose the problem is that you updates some libraries with conda-forge, and there is some conflict between other libraries from the default channel.