gdal/wgs84toxy.py

11 lines
417 B
Python
Raw Permalink Normal View History

import transformer
import pandas as pd
if __name__=='__main__':
proj=transformer.Projector(transformer.EPSG_WGS_84,transformer.EPSG_CGCS_2000_108E)
bzx=pd.read_excel(r'd:\3040\Desktop\宝中线.xlsx')
lat=bzx['纬度']
lon=bzx['经度']
y,x=proj.transform(lat.to_numpy(),lon.to_numpy())
xy_df=pd.DataFrame({'y':y,'x':x})
print(xy_df)
xy_df.to_excel(r'd:\3040\Desktop\宝中线xy.xlsx')