Wie finden Sie, ob die Koordinaten von Interesse innerhalb eines Polygons gefunden werden oder nicht?

Question

Ich habe einen Code geschrieben, um herauszufinden, ob die Koordinaten innerhalb oder außerhalb eines Rechtecks ​​liegen. Ich habe festgestellt, dass es eine Funktion "contain_points" gibt, die dem Zweck dient, äquivalent zu "inpolygon" in Matlab. Ich konnte kein Dokument über die Implementierung oder ein Beispiel für diese Funktion finden. Kann jemand vorschlagen, wie das funktioniert?

Answer 1

Werfen Sie einen Blick auf shapely Bibliothek. Es ist der De-facto-Standard beim Umgang mit 2D-Polygonen in Python (und vielen anderen Sprachen, da es ein Wrapping der GEOS C++ - Bibliothek ist).

Ein Beispiel von dem, was Sie wollen:

from shapely.geometry import Point, Polygon 
p = Point(0.0, 0.0) 
poly = Polygon([(1, 1), (-1, 1), (-1, -1), (1, -1)]) 
is_included = poly.contains(p) 

is_includedTrue jetzt ist.

Answer 2

@Eguaio, ich habe versucht, Ihren Code zu ändern. Aber ich bekomme das Ergebnis, "Punkte sind nicht im Polygon enthalten".

from shapely.geometry import Point,Polygon 
a = Point([(0.0,0.0)]) 
b = Polygon([(1,1),(-1,1),(-1,-1),(1,1)]) 
if b.contains(a): 
    print "the points are included in the polygon" 
else: 
    print "the points are not included in the polygon"