df_data = pd.read_excel(os.path.join(r"J:\LUCC_suc35\PRO\f analysis\d cluster\a table", "city_station.xls"))
geoseries= [Point(xy) for xy in zip(df_data.lon, df_data.lat)]gdf = gpd.GeoDataFrame(df_data, crs=4326, geometry=geoseries) gdf.to_file(os.path.join(r"J:\LUCC_suc35\PRO\f analysis\d cluster\a table", "bb.shp"))