首先,我们找到了相关材料,即matplotlib的一个方法:
https://matplotlib.org/api/path_api.html
import numpy as np
import pandas as pd
import pickle
from matplotlib.path import Path
with open('kmeans_times_0_version_2.pkl', 'rb') as f: # 要判断的经纬度数据
plots = pickle.load(f)
# 例如我要判断该经纬度点在不在广州
boundary = pd.read_csv('guangzhou_wgs84.csv') # 广州外多边形的WGS84坐标经纬度点
path = Path(boundary, closed=True)
is_get_on_in_guangzhou = path.contains_points(plots_get_on[['GET_ON_LONGITUDE', 'GET_ON_LATITUDE']])
plots = plots[is_get_on_in_guangzhou]
plots