1. 导入相关库
import pandas as pd
from osgeo import gdal
from shapely import wkt
import geopandas as gpd
2. 类型转换
- Geopandas字符串类型的几何数据转为geometry类型
geometry_csv=pd.read_csv('data/merged.csv')
geometry_csv['geometry'] = geometry_csv['geometry'].apply(wkt.loads)
geometry = gp.GeoDataFrame(geometry_csv, geometry='geometry')
geometry.tail(2)
3. csv文件转为shp文件完整代码
data=pd.read_csv(r'data/geo.csv')
data['lon']=data['lon'].astype('str')
data['lat']=data['lat'].astype('str')
data['point']="point ("+data['lon']+" "+data['lat']+")"
data['point'] = data['point'].apply(wkt.loads)
out_shp = gpd.GeoDataFrame(data)
out_shp.set_geometry('point',inplace=True)
out_shp.to_file('data/geo.shp',crs='EPSG:4326',encoding='utf-8')