ользую Python, и я определил широты и долготы (в градусах) многоугольника на карте. Моя цель - проверить, является ли общая точкаP координатx,y попадает в такой многоугольник. Поэтому я хотел бы иметь функцию, которая позволяет мне проверять такое состояние и возвращатьTrue или жеFalse если точка находится внутри или снаружи многоугольника.

В этом примере точка находится снаружи, поэтому результат будетFalse

Вопрос: Есть ли библиотека / пакет, который позволяет достичь моей цели? если да, какой из них вы рекомендуете? Вы могли бы привести небольшой пример того, как его использовать?

Вот код, который я написал до сих пор:

import numpy as np

# Define vertices of polygon (lat/lon)
v0 = [7.5, -2.5] 
v1 = [2, 3.5]
v2 = [-2, 4]
v3 = [-5.5, -4]
v4 = [0, -10]
lats_vect = np.array([v0[0],v1[0],v2[0],v3[0],v4[0]])
lons_vect = np.array([v0[1],v1[1],v2[1],v3[1],v4[1]])

# Point of interest P
x, y = -6, 5 # x = Lat, y = Lon

## START MODIFYING FROM HERE; DO NOT MODIFY POLYGON VERTICES AND DATA TYPE
# Check if point of interest falls within polygon boundaries
# If yes, return True
# If no, return falls

Чтобы нарисовать многоугольник и интересную точку, я использовал cartopy и написал следующие строки кода:

import cartopy.crs as ccrs
import matplotlib.pyplot as plt
ax = plt.axes(projection=ccrs.PlateCarree())
ax.stock_img() 

# Append first vertex to end of vector to close polygon when plotting
lats_vect = np.append(lats_vect, lats_vect[0])
lons_vect = np.append(lons_vect, lons_vect[0])
plt.plot([lons_vect[0:-1], lons_vect[1:]], [lats_vect[0:-1], lats_vect[1:]],
         color='black', linewidth=1, 
         transform=ccrs.Geodetic(),
         )   

plt.plot(y, x, 
        '*',          # marker shape
        color='blue',  # marker colour
        markersize=8  # marker size
        )  

plt.show()  

Заметка:

точки связаны друг с другом Великими кругами!Я исследовал в промежутке, и я закончил тем, что нашел некоторые подобные вопросыкак этот но у меня не было успеха, так как все они используют.shp файлы, которых у меня нет.