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
投稿猜你喜欢
- jQuery由美国人John Resig创建，至今已吸引了来自世界各地的众多javascript高手加入其team，包括来自德国的J&
前言YOLOV5模型从发布到现在都是炙手可热的目标检测模型，被广泛运用于各大场景之中。因此，我们不光要知道如何进行yolov5模型的训练，而- -- 建立测试表和数据create table webuser (username varchar(10));insert into web
- 1、安装scikit-learn1.1Scikit-learn 依赖Python (>= 2.6 or >= 3.3),NumP
- 1. 获取时间函数功能返回格式CURDATE(), CURRENT_DATE()返回当前日期，只包含年、月、日2022-03-08CURTI
- 题目：有四个数字：1、2、3、4，能组成多少个互不相同且无重复数字的三位数？各是多少？程序分析：可填在百位、十位、个位的数字都是1、2、3、
- 1.直方图直方图: (1) 图像中不同像素等级出现的次数 (2) 图像中具有不同等级的像素关于总像素数目的比值。我们使用cv2.calcHi
- 在python3.6版本中去掉了os.path.walk()函数os.walk()函数声明：walk(top,topdown=True,on
- 字典的创建名称={”键名“：”键值“}举例：person={
- 前言今天也是偏实战的内容，作为一个并发复习课，很简单，我们来看看怎样实现一个并发任务 manager。在微服务的场景下，我们有很多任务的执行
- 1.使用render方法return render(request,'index.html') 返回的页面内容是index.
- 使用一个遵循buffer protocol的对象就可以和numpy交互了.这个buffer_protocol要有哪些东西呢? 要有如下接口:
- jinjia和vue.js默认的模板转义符都是{{}}目前的解决办法是修改vue.js的转义符,将原来的{{}}替换为其他标签,我改为{[]
- 如果有空格就用％20代替，如果有其它字符就用％ASCII代替，如果有汉字等四个字节的字符，就用两个％ASCII来代替。不过有时候我们也需要将
- 1. *表示匹配任意多个字符 \d*表示匹配任意多个数字字符import retext = "
- 创建项目scrapy startproject zhaoping创建爬虫cd zhaopingscrapy genspider hr zha
- 曾有位网友遇到这样一个问题：产品名称为“A&T Plastic”，在产品列表中需要做这样的超链接，<a href="
- 视频观看视频函数的参数定义函数时，我们把参数的名字和位置确定下来，函数的接口定义就完成了。参数在函数名后的括号内指定。您可以根据需要添加任意
- 在使用 SQLAlchemy 时，那些看似很小的选择可能对这种对象关系映射工具包的性能产生重要影响。对象关系映射Object-relatio
- 在pycharm中，可以通过venv来建立工程，运行等等。但是一旦把这个工程文件夹拿到其他地方运行，而且不是在venv环境中运行，就有可能遇
