'GIS-like' operations about geopandas HOT 5 CLOSED

This would be powerful. I might call the parameter something other than by_geom; perhaps something like all_geoms? We could implement this API immediately (albeit inefficiently) by using unary_union on the other geometry, couldn't we?

from geopandas.

@kjordahl yes and no: The standard GIS 'layer-style' way of working is based on something like this:

for each geom in layer1:
    for each geom in layer2:
        difference layer1 geom by layer2 geom
        update layer1 geom

which is equivalent to using unary_union on the other geometry first. However, when it comes to things like intersect or union, things are slightly different. Now, additional geometries might be created, because parts of layer1 geom might intersect with several layer2 geoms, leading to a larger number of geoms in the output:

for each geom in layer1:
    for each geom in layer2:
        intersect layer1 geom with layer2 geom
        add intersect to output

Having said all that, since we aren't necessarily working in the realm of layers with GeoPandas, perhaps we should drop this conceptual model and adopt an API that is equivalent to a unary_union on other first (with more efficient use of spatial indices down the line). That way instead of saying "intersect this layer with this other layer", we are saying "intersect each of these geometries with these other ones", which I think fits the pandas way of doing things better. Alternatively, if we wanted to stick with a layer concept, then we could always just return a geometry collection for things like intersect and union? The only problem with this is that fewer operations support this, and it would definitely make i/o more difficult for the user.
Long story short: yes, we can do this now if we are willing to do things differently than a GIS (which I think we should be willing to do). Also, having the control afforded by all_geoms means we should probably still be able to do the GIS 'layer-approach' with some clever index alignments..?

from geopandas.

If we decide which way to go, I'll try to submit a pull request for some of this over the weekend.

from geopandas.

@cfarmer Do the spatial join and overlay functionality cover everything in this issue? Or are there additional "GIS-like" operations/concepts that we should look at?

from geopandas.

Closing for now. If something's missing from sjoin and overlay, let's bring up a new issue.

from geopandas.