{"id":10320,"date":"2024-07-24T22:23:28","date_gmt":"2024-07-24T13:23:28","guid":{"rendered":"https:\/\/developers.agirobots.com\/jp\/?p=10320"},"modified":"2024-07-24T22:24:51","modified_gmt":"2024-07-24T13:24:51","slug":"point-in-polygon-python","status":"publish","type":"post","link":"https:\/\/developers.agirobots.com\/jp\/point-in-polygon-python\/","title":{"rendered":"Python\u3067\u70b9\u306e\u9818\u57df\u5185\u5916\u5224\u5b9a\u3092\u3059\u308b\u65b9\u6cd5"},"content":{"rendered":"\n

\u81ea\u52d5\u8d70\u884c\u30ed\u30dc\u30c3\u30c8\u306e\u8d70\u884c\u30d7\u30ed\u30b0\u30e9\u30e0\u3092\u4f5c\u6210\u3057\u3066\u3044\u308b\u3068\u3001\u4e8b\u524d\u306b\u6c7a\u3081\u305f\u9818\u57df\u306b\u3044\u308b\u306e\u304b\u3044\u306a\u3044\u306e\u304b\u3001\u3068\u3044\u3046\u3053\u3068\u3092\u5224\u5b9a\u3057\u305f\u304f\u306a\u308b\u3053\u3068\u304c\u7d50\u69cb\u3042\u308a\u307e\u3059\u3002\u3053\u3053\u3067\u306f\u3001\u305d\u306e\u65b9\u6cd5\u306b\u3064\u3044\u3066\u8abf\u3079\u305f\u5185\u5bb9\u306e\u5099\u5fd8\u9332\u3092\u4f5c\u6210\u3057\u307e\u3057\u305f\u3002<\/p>\n\n\n\n

\u7d50\u8ad6<\/h2>\n\n\n\n

shapely\u3068\u3044\u3046\u4fbf\u5229\u306a\u30e9\u30a4\u30d6\u30e9\u30ea\u304c\u516c\u958b\u3055\u308c\u3066\u3044\u308b\u306e\u3067\u3001\u305d\u308c\u3092\u4f7f\u7528\u3059\u308b\u3068\u7c21\u5358\u306b\u5b9f\u73fe\u3067\u304d\u307e\u3059\u3002<\/p>\n\n\n\n

shapely\u3092\u7528\u3044\u3066\u9818\u57df\u5185\u5916\u5224\u5b9a\u3092\u3059\u308b\u65b9\u6cd5<\/h2>\n\n\n\n

shapely\u3092\u7528\u3044\u305f\u9818\u57df\u5224\u5b9a\u306f\u3001\u591a\u89d2\u5f62\u3068\u3057\u3066\u9818\u57df\u3092\u5b9a\u7fa9\u3001\u5224\u5b9a\u3068\u3044\u3046\u9806\u756a\u3067\u884c\u3044\u307e\u3059\u3002<\/p>\n\n\n\n

\u4f8b\u3048\u3070\u3001\u4ee5\u4e0b\u306e\u3088\u3046\u306b\uff14\u3064\u306e\u70b9\u3067\u8868\u3055\u308c\u308b\u9818\u57df\u306e\u5185\u5916\u5224\u5b9a\u3092\u3059\u308b\u3053\u3068\u3092\u8003\u3048\u3066\u307f\u307e\u3059\u3002
\u70b9\u306f\u3001\u5de6\u304b\u3089\u9806\u756a\u306b\u3001(0.8, 4.3)\u3001(1.2, 1)\u3001(2.5, 5.7)\u3001 (5.1, 3.4)\u3068\u3057\u307e\u3059\u3002<\/p>\n\n\n

\n
\"\"<\/figure><\/div>\n\n\n

\u307e\u305a\u306f\u5fc5\u8981\u306a\u30af\u30e9\u30b9\u3092\u30a4\u30f3\u30dd\u30fc\u30c8\u3057\u307e\u3059\u3002\u4f7f\u7528\u3059\u308b\u306e\u306f\u3001Point\u3068\u3001Polygon\u3067\u3059\u3002<\/p>\n\n\n\n

from shapely.geometry import Point, Polygon<\/code><\/pre><\/div>\n\n\n\n
Polygon\u30af\u30e9\u30b9\u3092\u7528\u3044\u3066\u3001\u5148\u307b\u3069\u306e\uff14\u3064\u306e\u70b9\u3067\u8868\u3055\u308c\u308b\u30dd\u30ea\u30b4\u30f3\u306f\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u3057\u3066\u8a2d\u5b9a\u3059\u308b\u3053\u3068\u304c\u3067\u304d\u307e\u3059\u3002<\/p>\n\n\n\n
# 4\u3064\u306e\u70b9\u3092\u5b9a\u7fa9\np1 = (1.2, 1)\np2 = (0.8, 4.3)\np3 = (2.5, 5.7)\np4 = (5.1, 3.4)\n\n# \u591a\u89d2\u5f62\u3092\u5b9a\u7fa9\npolygon = Polygon([p1, p2, p3, p4])<\/code><\/pre><\/div>\n\n\n\n
\u3067\u306f\u3001\u4efb\u610f\u306e\u70b9\u304c\u3053\u306e\u30dd\u30ea\u30b4\u30f3\u306e\u5185\u5916\u306b\u3042\u308b\u306e\u304b\u3092\u5224\u5b9a\u3057\u3066\u307f\u307e\u3057\u3087\u3046\u3002<\/p>\n\n\n\n
\u4eca\u56de\u306f\u5177\u4f53\u4f8b\u3068\u3057\u3066\u3001(3, 3)\u3068(5, 2)\u306b\u95a2\u3057\u3066\u5224\u5b9a\u3057\u3066\u307f\u308b\u3053\u3068\u306b\u3057\u307e\u3059\u3002\u70b9\u306f\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u8868\u3057\u307e\u3059\u3002<\/p>\n\n\n\n
point33 = Point((3, 3))\npoint52 = Point((5, 2))<\/code><\/pre><\/div>\n\n\n\n
\u5224\u5b9a\u306f\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u3057\u3066\u5b9f\u884c\u3057\u307e\u3059\u3002<\/p>\n\n\n\n
is_inside33 = polygon.contains(point33)\nis_inside52 = polygon.contains(point52)\n\nprint(is_inside33)  # \u7d50\u679c\u306fTrue\nprint(is_inside52)  # \u7d50\u679c\u306fFalse<\/code><\/pre><\/div>\n\n\n\n
\u30e9\u30f3\u30c0\u30e0\u306a\u5927\u91cf\u306a\u70b9\u306b\u3064\u3044\u3066\u9818\u57df\u5185\u5916\u306e\u5224\u5b9a\u3092\u3057\u3066\u307f\u308b<\/h2>\n\n\n\n
xy\u8ef8\u305d\u308c\u305e\u308c\u306b\u3064\u3044\u3066\u30010~6\u306e\u7bc4\u56f2\u3067\u30e9\u30f3\u30c0\u30e0\u306b\u30b5\u30f3\u30d7\u30ea\u30f3\u30b0\u3057\u3066\u70b9\u3092\u5b9a\u7fa9\u3057\u3001\u305d\u308c\u304c\u9818\u57df\u306e\u5185\u5074\u304b\u5916\u5074\u304b\u3092\u63cf\u753b\u3057\u3066\u307f\u305f\u3068\u3053\u308d\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u306a\u308a\u307e\u3057\u305f\u3002\u5185\u5074\u306e\u70b9\u306f\u8d64\u70b9\u3001\u5916\u5074\u306e\u70b9\u306f\u9ec4\u8272\u70b9\u3067\u8868\u3057\u3066\u3044\u307e\u3059\u3002<\/p>\n\n\n
\n
\"\"<\/figure><\/div>\n\n\n
\u3053\u308c\u3092\u8a66\u3057\u306b\u5b9f\u884c\u3057\u3066\u307f\u305f\u65b9\u306f\u4ee5\u4e0b\u306e\u30b3\u30fc\u30c9\u3092\u5b9f\u884c\u3057\u3066\u307f\u3066\u304f\u3060\u3055\u3044\u3002<\/p>\n\n\n\n
import numpy as np\nimport matplotlib.pyplot as plt\nfrom shapely.geometry import Point, Polygon\n\n# 4\u3064\u306e\u70b9\u3092\u5b9a\u7fa9\np1 = (1.2, 1)\np2 = (0.8, 4.3)\np3 = (2.5, 5.7)\np4 = (5.1, 3.4)\n\n# \u591a\u89d2\u5f62\u3092\u5b9a\u7fa9\npolygon = Polygon([p1, p2, p3, p4])\n\n# \u30e9\u30f3\u30c0\u30e0\u306a\u70b9\u3092\u751f\u6210\nnum_points = 1000\nrandom_points = np.random.rand(num_points, 2) * 6  # 0\u304b\u30896\u306e\u7bc4\u56f2\u3067\u30e9\u30f3\u30c0\u30e0\u306a\u70b9\u3092\u751f\u6210\n\n# \u591a\u89d2\u5f62\u5185\u306b\u3042\u308b\u70b9\u3092\u8d64\u8272\u3001\u5916\u306b\u3042\u308b\u70b9\u3092\u7dd1\u8272\u3067\u5206\u985e\ninside_points = []\noutside_points = []\n\nfor point in random_points:\n    if polygon.contains(Point(point)):\n        inside_points.append(point)\n    else:\n        outside_points.append(point)\n\ninside_points = np.array(inside_points)\noutside_points = np.array(outside_points)\n\n# \u591a\u89d2\u5f62\u3092\u63cf\u753b\npolygon_x, polygon_y = polygon.exterior.xy\nplt.plot(polygon_x, polygon_y, 'b-', label='Polygon')\n\n# \u591a\u89d2\u5f62\u5185\u306e\u70b9\u3092\u8d64\u8272\u3001\u5916\u306e\u70b9\u3092\u7dd1\u8272\u3067\u63cf\u753b\nif inside_points.size > 0:\n    plt.plot(inside_points[:, 0], inside_points[:, 1], 'ro', markersize=2, label='Inside Points')\nif outside_points.size > 0:\n    plt.plot(outside_points[:, 0], outside_points[:, 1], 'yo', markersize=2, label='Outside Points')\n\n# # \u5404\u9802\u70b9\u3092\u63cf\u753b\nplt.plot(*zip(*[p1, p2, p3, p4]), 'bo', label='Vertices')\n\n# \u63cf\u753b\u306e\u8a2d\u5b9a\nplt.xlabel('X')\nplt.ylabel('Y')\nplt.title('Random Points in Polygon')\nplt.legend()\nplt.grid(True)\n\n# \u30b0\u30e9\u30d5\u3092\u8868\u793a\nplt.show()<\/code><\/pre><\/div>\n\n\n\n
<\/p>\n","protected":false},"excerpt":{"rendered":"
\u81ea\u52d5\u8d70\u884c\u30ed\u30dc\u30c3\u30c8\u306e\u8d70\u884c\u30d7\u30ed\u30b0\u30e9\u30e0\u3092\u4f5c\u6210\u3057\u3066\u3044\u308b\u3068\u3001\u4e8b\u524d\u306b\u6c7a\u3081\u305f\u9818\u57df\u306b\u3044\u308b\u306e\u304b\u3044\u306a\u3044\u306e\u304b\u3001\u3068\u3044\u3046\u3053\u3068\u3092\u5224\u5b9a\u3057\u305f\u304f\u306a\u308b\u3053\u3068\u304c\u7d50\u69cb\u3042\u308a\u307e\u3059\u3002\u3053\u3053\u3067\u306f\u3001\u305d\u306e\u65b9\u6cd5\u306b\u3064\u3044\u3066\u8abf\u3079\u305f\u5185\u5bb9\u306e\u5099\u5fd8\u9332\u3092\u4f5c\u6210\u3057\u307e\u3057\u305f\u3002 \u7d50\u8ad6 … <\/p>\n","protected":false},"author":1,"featured_media":10326,"comment_status":"closed","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"_vk_print_noindex":"","sitemap_hide":"","_veu_custom_css":"","veu_display_promotion_alert":"common","vkexunit_cta_each_option":"","_jetpack_memberships_contains_paid_content":false,"footnotes":""},"categories":[170,136],"tags":[],"class_list":["post-10320","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-blog","category-robot-research-make"],"veu_head_title_object":{"title":"","add_site_title":""},"jetpack_featured_media_url":"https:\/\/developers.agirobots.com\/jp\/wp-content\/uploads\/2024\/07\/image-2.png","jetpack-related-posts":[],"jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/developers.agirobots.com\/jp\/wp-json\/wp\/v2\/posts\/10320","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/developers.agirobots.com\/jp\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/developers.agirobots.com\/jp\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/developers.agirobots.com\/jp\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/developers.agirobots.com\/jp\/wp-json\/wp\/v2\/comments?post=10320"}],"version-history":[{"count":3,"href":"https:\/\/developers.agirobots.com\/jp\/wp-json\/wp\/v2\/posts\/10320\/revisions"}],"predecessor-version":[{"id":10327,"href":"https:\/\/developers.agirobots.com\/jp\/wp-json\/wp\/v2\/posts\/10320\/revisions\/10327"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/developers.agirobots.com\/jp\/wp-json\/wp\/v2\/media\/10326"}],"wp:attachment":[{"href":"https:\/\/developers.agirobots.com\/jp\/wp-json\/wp\/v2\/media?parent=10320"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/developers.agirobots.com\/jp\/wp-json\/wp\/v2\/categories?post=10320"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/developers.agirobots.com\/jp\/wp-json\/wp\/v2\/tags?post=10320"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}