创建Shapefile文件并写入数据的例子


Posted in Python onNovember 26, 2019

基本思路

使用GDAL创建Shapefile数据的基本步骤如下:

使用osgeo.ogr.Driver的CreateDataSource()方法创建osgeo.ogr.DataSource矢量数据集

使用osgeo.ogr.DataSource的CreateLayer()方法创建一个图层

使用osgeo.ogr.FieldDefn()定义Shapefile文件的属性字段

创建osgeo.ogr.Feature对象,设置每个属性字段的值,使用Feature对象的SetGeometry()定义几何属性

创建Feature对象以后,使用osgeo.ogr.Layer的CreateFeature()添加Feature对象到当前图层

重复步骤4和5依次添加所有的Feature到当前图层即可

代码实现

下面的例子中,我们读取GeoJSON表示的中国省区数据,然后其转为Shapefile格式。

GeoJSON编码片段如下:

创建Shapefile文件并写入数据的例子

可以看到每个Feature都有一个properties字段和geometry字段,我们需要根据properties字段的信息创建Shapefile数据的属性表,根据geometry字段创建Shapefile中的几何数据。

from osgeo import ogr
from osgeo import osr
import json
import os
os.environ['SHAPE_ENCODING'] = "utf-8"


with open('China.json') as f:
 china = json.load(f)

# 创建DataSource
driver = ogr.GetDriverByName('ESRI Shapefile')
ds = driver.CreateDataSource('China.shp')

# 创建WGS84空间参考
srs = osr.SpatialReference()
srs.ImportFromEPSG(4326)

# 创建图层
layer = ds.CreateLayer('province', srs, ogr.wkbPolygon)
# 添加属性定义
fname = ogr.FieldDefn('Name', ogr.OFTString)
fname.SetWidth(24)
layer.CreateField(fname)
fcx = ogr.FieldDefn('CenterX', ogr.OFTReal)
layer.CreateField(fcx)
fcy = ogr.FieldDefn('CenterY', ogr.OFTReal)
layer.CreateField(fcy)

# 变量GeoJSON中的features
for f in china['features']:
 # 新建Feature并且给其属性赋值
 feature = ogr.Feature(layer.GetLayerDefn())
 feature.SetField('Name', f['properties']['name'])
 feature.SetField('CenterX', f['properties']['cp'][0])
 feature.SetField('CenterY', f['properties']['cp'][1])

 # 设置Feature的几何属性Geometry
 polygon = ogr.CreateGeometryFromJson(str(f['geometry']))
 feature.SetGeometry(polygon)
 # 创建Feature
 layer.CreateFeature(feature)
 del feature
ds.FlushCache()

del ds

以上这篇创建Shapefile文件并写入数据的例子就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
PHP魔术方法__ISSET、__UNSET使用实例
Nov 25 Python
Mac下Supervisor进程监控管理工具的安装与配置
Dec 16 Python
分析用Python脚本关闭文件操作的机制
Jun 28 Python
Django自定义分页效果
Jun 27 Python
浅谈使用Python变量时要避免的3个错误
Oct 30 Python
如何用python整理附件
May 13 Python
详解Numpy中的广播原则/机制
Sep 20 Python
Python 获取div标签中的文字实例
Dec 20 Python
Python元组常见操作示例
Feb 19 Python
Python提取特定时间段内数据的方法实例
Apr 01 Python
Python变量作用域LEGB用法解析
Feb 04 Python
Django单元测试中Fixtures用法详解
Feb 25 Python
python使用opencv在Windows下调用摄像头实现解析
Nov 26 #Python
使用Python实现 学生学籍管理系统
Nov 26 #Python
python redis 批量设置过期key过程解析
Nov 26 #Python
python3 tkinter实现添加图片和文本
Nov 26 #Python
使用Rasterio读取栅格数据的实例讲解
Nov 26 #Python
Python小程序之在图片上加入数字的代码
Nov 26 #Python
基于Python获取城市近7天天气预报
Nov 26 #Python
You might like
php xml-rpc远程调用
2008/12/19 PHP
Linux Apache PHP Oracle 安装配置(具体操作步骤)
2013/06/17 PHP
php获取新浪微博数据API实例
2013/11/12 PHP
php查看当前Session的ID实例
2015/03/16 PHP
PHP实现时间日期友好显示实现代码
2019/09/08 PHP
基于jquery跨浏览器显示的file上传控件
2011/10/24 Javascript
js实现特定位取反原理及示例
2014/06/30 Javascript
JS模仿编辑器实时改变文本框宽度和高度大小的方法
2015/08/17 Javascript
jquery.mousewheel实现整屏翻屏效果
2015/08/30 Javascript
jQuery EasyUI实现右键菜单变灰不可用效果
2015/09/24 Javascript
Angularjs全局变量被作用域监听的正确姿势
2016/02/06 Javascript
使用jQuery或者原生js实现鼠标滚动加载页面新数据
2016/03/06 Javascript
javascript检测移动设备横竖屏
2016/05/21 Javascript
用js实现简单算法的实例代码
2016/09/24 Javascript
NodeJS中的MongoDB快速入门详细教程
2016/11/11 NodeJs
JS实现的集合去重,交集,并集,差集功能示例
2018/03/13 Javascript
动态加载、移除js/css文件的示例代码
2018/03/20 Javascript
JS自定义滚动条效果
2020/03/13 Javascript
[48:32]VGJ.T vs Fnatic 2018国际邀请赛小组赛BO2 第一场 8.16
2018/08/17 DOTA
python通过zlib实现压缩与解压字符串的方法
2014/11/19 Python
Python面向对象class类属性及子类用法分析
2018/02/02 Python
使用python绘制3维正态分布图的方法
2018/12/29 Python
Python自动化导出zabbix数据并发邮件脚本
2019/08/16 Python
详解在Python中以绝对路径或者相对路径导入文件的方法
2019/08/30 Python
python要安装在哪个盘
2020/06/15 Python
如何基于Python代码实现高精度免费OCR工具
2020/06/18 Python
python3中celery异步框架简单使用+守护进程方式启动
2021/01/20 Python
H5新属性audio音频和video视频的控制详解(推荐)
2016/12/09 HTML / CSS
FLOS美国官网:意大利高级照明工艺的传奇
2018/08/07 全球购物
车间班组长的职责
2013/12/13 职场文书
园林设计专业毕业生求职信
2014/03/23 职场文书
说谎欺骗人检讨书300字
2014/11/18 职场文书
2014年师德师风工作总结
2014/11/25 职场文书
2015年车间安全管理工作总结
2015/05/13 职场文书
python如何正确使用yield
2021/05/21 Python
mysql数据插入覆盖和时间戳的问题及解决
2022/03/25 MySQL