利用geopandas进行投影变换（以WGS1984投影到WGS1984—UTM—48N为例）

import geopandas
import rasterio.warp
from shapely.geometry import shape

# 导入需要投影的点shapefile
file_name = "hydrologic_station.shp"
data = gpd.read_file(file_name)

# 事先利用arcgis，将点shapefile文件进行UTM48投影，读取投影后的文件，查询目标crs
file_name1 = "station.shp"
data_project = gpd.read_file(file_name1)
destination_crs = data_project.crs

geometry = rasterio.warp.transform_geom(
    src_crs=data.crs,
    dst_crs=destination_crs,
    geom=data.geometry.values,
)
data_project_check = data.set_geometry(
    [shape(geom) for geom in geometry],
    crs=destination_crs,
) #确认geopandas投影后的点shapefile是否与arcgis投影后的保持一致

#之后获得批量经纬度坐标，可以利用destination_crs直接投影