Skip to content

Instantly share code, notes, and snippets.

@nmmmnu

nmmmnu/point_in_polygon.py

Forked from inside-code-yt/point_in_polygon.py
Created July 31, 2023 07:38
Show Gist options

  • Select an option

    Learn more about clone URLs
  • Save nmmmnu/5017c1bed9976139137f6b544e6e1052 to your computer and use it in GitHub Desktop.

Select an option

Learn more about clone URLs
Save nmmmnu/5017c1bed9976139137f6b544e6e1052 to your computer and use it in GitHub Desktop.
Download ZIP

Revisions

  1. @inside-code-yt inside-code-yt created this gist Feb 26, 2023.
    34 changes: 34 additions & 0 deletions point_in_polygon.py
    View file
    Original file line number Diff line number Diff line change
    @@ -0,0 +1,34 @@
    from polygenerator import random_polygon
    import matplotlib.pyplot as plt


    def is_inside(edges, xp, yp):
    cnt = 0
    for edge in edges:
    (x1, y1), (x2, y2) = edge
    if (yp < y1) != (yp < y2) and xp < x1 + ((yp-y1)/(y2-y1))*(x2-x1):
    cnt += 1
    return cnt%2 == 1


    def onclick(event):
    xp, yp = event.xdata, event.ydata
    if is_inside(edges, xp, yp):
    print("inside")
    plt.plot(xp, yp, "go", markersize=5)
    else:
    print("outside")
    plt.plot(xp, yp, "ro", markersize=5)
    plt.gcf().canvas.draw()


    polygon = random_polygon(num_points=20)
    polygon.append(polygon[0])
    edges = list(zip(polygon, polygon[1:] + polygon[:1]))
    plt.figure(figsize=(10, 10))
    plt.gca().set_aspect("equal")
    xs, ys = zip(*polygon)
    plt.gcf().canvas.mpl_connect('button_press_event', onclick)
    plt.plot(xs, ys, "b-", linewidth=0.8)
    plt.show()