I have a feature class which contains 300,000 polygon and multipolygon features in the same fc. The goal is to get a result layer of spatial join (must to be one to many join) between this feature class to itself and during the process I have to only use open source GIS software, better code than geoprocessing tools. I tried:
GeoPandas sjoin - doesn't have one to many option.
QGIS 1 - Tried to read the fc in QGIS then use join by location. The join crushed after 20%.
QGIS 2 - Tried to convert fc to shapefile then use multipart to singlepart than dissolve to get the first mix of polygon and multipolygon. Again the join crushed in the same time.
Mabye QGIS cannot handle this kind of geometries mix?