from shapely.geometry import Point,Polygon
import shapely.geometry as geometry
import shapely.vectorized
from shapely.geometry import MultiPolygon, Polygon
from shapely.geometry import shape,Point
from shapely.wkt import loads
df0 = pd.read_excel('boundary.xlsx')
road = df0.iloc[:,[1,2]].values
for i in range(np.shape(road)[0]):
Main_poly = loads(road[i,1]) # 导入028边界
lons = df_1['lon'].tolist()
lats = df_1['lat'].tolist()
flag_condition=shapely.vectorized.contains(Main_poly, lons, lats)
temp_array1 = np.zeros(len(df_1))
temp_array1[flag_condition] = 1
df_1[str(road[i,0])] = temp_array1
#data_pickup_sample20211210_plot = df_1[df_1['flag']==1]