#encoding:utf-8 import pandas as pd import geopandas import pyproj df= geopandas.read_file(r'./static/input/pointtest.csv',encoding='utf-8') df[['lon','lat']] = df[['lon','lat']].apply(pd.to_numeric) gdf = geopandas.GeoDataFrame(df, geometry=geopandas.points_from_xy(df.lon, df.lat)) gdf.crs = pyproj.CRS.from_user_input('EPSG:4326') #给输出的shp增加投影 # gdf.rename(columns={'name':'ave_price'},inplace=True) #对原来的字段名进行更改 # gdf.rename(columns={'addrees':'area_ave_price'},inplace=True) gdf.to_file(r'./static/output/pointfeatures.shp', driver='ESRI Shapefile',encoding='utf-8')
利用geopandas将csv文件转为shp
最新推荐文章于 2024-11-05 11:05:00 发布