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

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

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

image

About this issue

  • Original URL
  • State: open
  • Created 6 years ago
  • Comments: 20 (20 by maintainers)

Most upvoted comments

@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.

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.