Python 根据数据模板创建shapefile的实现
废话不多说，我就直接上代码让大家看看吧！
#!/usr/bin/envpython
#-*-coding:utf-8-*-
#@File:copyShapefile.py
#@Author:huifer
#@Date:2018-4-28
fromos.pathimportexists
importgdal
fromosgeoimportogr
fromosimportremove
gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8","YES")#路径中文
gdal.SetConfigOption("SHAPE_ENCODING","GBK")#属性中文
in_shapefile="dataSample/wang_point.shp"#数据模板
out_shapefile="shapefileAa.shp"#输出数据集
in_ds=ogr.Open(in_shapefile)#读取模板数据
in_lyr=in_ds.GetLayerByIndex(0)
ifexists(out_shapefile):
remove(out_shapefile)
drv=ogr.GetDriverByName("ESRIShapefile")#指定数据驱动
out_ds=drv.CreateDataSource(out_shapefile)#创建数据源
proj=in_lyr.GetSpatialRef()#获取模板坐标系
out_lyr=out_ds.CreateLayer(out_shapefile.split(".")[0],proj,ogr.wkbPoint)
#copytheschemaoftheoriginalshapefiletothedestinationshapefile
lyr_def=in_lyr.GetLayerDefn()
foriinrange(lyr_def.GetFieldCount()):#获取字段长度
out_lyr.CreateField(lyr_def.GetFieldDefn(i))#创建字段
feature=ogr.Feature(lyr_def)
wkt="POINT(88615.73000075345.486000)"
point=ogr.CreateGeometryFromWkt(wkt)
feature.SetGeometry(point)
#添加点
out_lyr.CreateFeature(feature)
#关闭特征
feature=None
#关闭数据
data_source=None
以上这篇Python根据数据模板创建shapefile的实现就是小编分享给大家的全部内容了，希望能给大家一个参考，也希望大家多多支持毛票票。
声明：本文内容来源于网络，版权归原作者所有，内容由互联网用户自发贡献自行上传，本网站不拥有所有权，未作人工编辑处理，也不承担相关法律责任。如果您发现有涉嫌版权的内容，欢迎发送邮件至：czq8825#qq.com（发邮件时，请将#更换为@）进行举报，并提供相关证据，一经查实，本站将立刻删除涉嫌侵权内容。