网络编程
位置:首页>> 网络编程>> Python编程>> Python 根据数据模板创建shapefile的实现

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

作者:staHuri  发布时间:2023-07-27 00:10:32 

标签:Python,创建,shapefile

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


#!/usr/bin/env python
# -*- coding: utf-8 -*-
# @File : copyShapefile.py
# @Author: huifer
# @Date : 2018-4-28
from os.path import exists

import gdal

from osgeo import ogr
from os import remove

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)
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 shapefile
lyr_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

来源:https://blog.csdn.net/staHuri/article/details/80115704

0
投稿

猜你喜欢

手机版 网络编程 asp之家 www.aspxhome.com