Python 根据数据模板创建shapefile的实现

Python  /  管理员 发布于 7年前   307

废话不多说，我就直接上代码让大家看看吧！

#!/usr/bin/env python# -*- coding: utf-8 -*-# @File : copyShapefile.py# @Author: huifer# @Date : 2018-4-28from os.path import existsimport gdalfrom osgeo import ogrfrom os import removegdal.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)if exists(out_shapefile):  remove(out_shapefile)drv = ogr.GetDriverByName("ESRI Shapefile") # 指定数据驱动out_ds = drv.CreateDataSource(out_shapefile) # 创建数据源proj = in_lyr.GetSpatialRef() # 获取模板坐标系out_lyr = out_ds.CreateLayer(out_shapefile.split(".")[0], proj, ogr.wkbPoint)# copy the schema of the original shapefile to the destination shapefilelyr_def = in_lyr.GetLayerDefn()for i in range(lyr_def.GetFieldCount()): # 获取字段长度  out_lyr.CreateField(lyr_def.GetFieldDefn(i)) # 创建字段  feature = ogr.Feature(lyr_def)  wkt = "POINT(88615.730000 75345.486000)"  point = ogr.CreateGeometryFromWkt(wkt)  feature.SetGeometry(point)  # 添加点  out_lyr.CreateFeature(feature)  # 关闭 特征  feature = None  # 关闭数据data_source = None

以上这篇Python 根据数据模板创建shapefile的实现就是小编分享给大家的全部内容了，希望能给大家一个参考，也希望大家多多支持。