Comment trouver les points d'intersection avec un polygone dans geopandas?

J'ai essayé d'utiliser la "coupe" sur une geodataframe, pour voir les points qui se trouvent à l'intérieur d'un polygone. Toutefois, seule la première fonction dans le cadre de sera de retour en tant que vrai. Ce que je fais mal?

from geopandas.geoseries import *

p1 = Point(.5,.5)
p2 = Point(.5,1)
p3 = Point(1,1)

g1 = GeoSeries([p1,p2,p3])
g2 = GeoSeries([p2,p3])

g = GeoSeries([Polygon([(0,0), (0,2), (2,2), (2,0)])])

g1.intersects(g) # Flags the first point as inside, even though all are.
g2.intersects(g) # The second point gets picked up as inside (but not 3rd)

OriginalL'auteur Thomas Pingel | 2015-05-22

  1. 12

    Selon la la documentation:

    Opérations binaires peut être appliqué entre deux GeoSeries, auquel cas
    l'opération est réalisée elementwise. Les deux séries seront
    aligné par correspondance des indices.

    Vos exemples ne sont pas censés travailler. Donc, si vous voulez tester pour chaque point d'être en un seul polygone, vous aurez à faire:

    poly = GeoSeries(Polygon([(0,0), (0,2), (2,2), (2,0)]))
g1.intersects(poly.ix[0])

    Sorties:

        0    True
    1    True
    2    True
    dtype: bool

    Ou si vous voulez tester pour toutes les géométries dans un GeoSeries:

    points.intersects(poly.unary_union)

    Geopandas s'appuie sur Galbées géométriques de travail. Il est parfois utile (et plus facile à lire) pour l'utiliser directement. Le code suivant fonctionne également comme annoncé:

    from shapely.geometry import *

p1 = Point(.5,.5)
p2 = Point(.5,1)
p3 = Point(1,1)

poly = Polygon([(0,0), (0,2), (2,2), (2,0)])

for p in [p1, p2, p3]:
    print(poly.intersects(p))

    Vous pouvez également jeter un oeil à
    Comment gérer les erreurs d'arrondi dans Galbées pour les problèmes qui peuvent survenir avec des points sur les limites.

    Cela a aidé beaucoup dans l'explication de ce qui se passait. Nous sommes en fait à l'aide d'un GeoDataFrame, et j'ai donc été coulée le problème de manière un peu plus simple. Ton exemple fonctionne très bien, mais c'est plus pratique pour nous de garder les polygones (et leurs attributs) dans un GeoDataFrame pour d'autres raisons, sinon la galbées façon fonctionne très bien.
    Bon! J'ai ajouté votre exemple de ma réponse, que je trouve plus complet (avec la ref de la doc). Vous pourriez envisager de marquer la question comme résolue et/ou upvoting.

    OriginalL'auteur Fabzi

  2. 3

    Un moyen de contourner cela semble être d'obtenir soit une entrée particulière (qui ne fonctionne pas pour mon application, mais peuvent travailler pour quelqu'un d'autre:

    from geopandas.geoseries import *

p1 = Point(.5,.5)
p2 = Point(.5,1)
p3 = Point(1,1)

points = GeoSeries([p1,p2,p3])

poly = GeoSeries([Polygon([(0,0), (0,2), (2,2), (2,0)])])

points.intersects(poly.ix[0])

    D'une autre façon (plus utile à ma demande) est à l'intersection avec unaire de l'union des fonctionnalités pour la deuxième couche:

    points.intersects(poly.unary_union)

    OriginalL'auteur Thomas Pingel

  3. 0

    Vous pouvez facilement vérifier les points qui se trouve à l'intérieur d'un polygone avec cette simple fonction ci-dessous:

    import geopandas
from shapely.geometry import *

p1 = Point(.5,.5)
p2 = Point(.5,1)
p3 = Point(1,1)

g = Polygon([(0,0), (0,2), (2,2), (2,0)])

def point_inside_shape(point, shape):
    #point of type Point
    #shape of type Polygon
    pnt = geopandas.GeoDataFrame(geometry=[point], index=['A'])
    return(pnt.within(shape).iloc[0])

for p in [p1, p2, p3]:
    print(point_inside_shape(p, g))

    OriginalL'auteur Ioannis Nasios