平面上の多角形中に点があるかどうか判定する

7878 ワード
shapely Python3 データサイエンス Python3 テキストリンク

下記の図のように、(5,5), (10,5), (10,10), (5,10)の青い点で囲まれた矩形があったとする。

ここに、以下のように(7,8)のオレンジの点と(12,10)の緑の点があったときに、矩形中にあるかどうかを判定する方法を紹介する。

Shapelyというライブラリを利用する。
これは、GEOSをベースとしたpythonライブラリで、平面上の幾何学図形の分析と操作のために使われる。

参考: https://pypi.org/project/Shapely/

$ pip install shapely

矩形の中に点があるかどうか

以下のコードで確認できる。

from shapely.geometry import Point
from shapely.geometry.polygon import Polygon

point = Point(7,8)
polygon = Polygon([(5,5), (10,5), (10,10), (5,10)])
print(polygon.contains(point)) # True

point = Point(12,10)
polygon = Polygon([(5,5), (10,5), (10,10), (5,10)])
print(polygon.contains(point)) # False

境界上の点は、矩形内の点ではないと判定される。

point = Point(5, 5)
polygon = Polygon([(5,5), (10,5), (10,10),(5,10)])
print(polygon.contains(point)) # False

5点以上の場合

以下の図のように、5点ある場合はどうするのか？（ (7.5, 12.5)を追加した。）

point = Point(7.5, 7.5)
polygon = Polygon([(5,5), (10,5), (10,10),(5,10), (7.5,12.5)]) # ここに追加するだけ。
print(polygon.contains(point)) # True

対象の図形を凸包として判定したい場合は、以下で確認可能だ。

from shapely.geometry import Point
from shapely.geometry.polygon import Polygon
from shapely.geometry import MultiPoint


point = Point(7.5, 7.5)
polygon = MultiPoint([(5,5), (10,5), (10,10),(5,10), (7.5,12.5)]).convex_hull
print(polygon.contains(point)) # True

この場合も、境界上の点はFalseと判定される。