基本思路
使用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编码片段如下：
可以看到每个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文件并写入数据的例子就是小编分享给大家的全部内容了，希望能给大家一个参考，也希望大家多多支持。
Shapefile,写入数据
免责声明：本站文章均来自网站采集或用户投稿，网站不提供任何软件下载或自行开发的软件！ 如有用户或公司发现本站内容信息存在侵权行为，请邮件告知！ 858582#qq.com
P70系列延期，华为新旗舰将在下月发布
3月20日消息，近期博主@数码闲聊站 透露，原定三月份发布的华为新旗舰P70系列延期发布，预计4月份上市。
而博主@定焦数码 爆料，华为的P70系列在定位上已经超过了Mate60，成为了重要的旗舰系列之一。它肩负着重返影像领域顶尖的使命。那么这次P70会带来哪些令人惊艳的创新呢？
根据目前爆料的消息来看，华为P70系列将推出三个版本，其中P70和P70 Pro采用了三角形的摄像头模组设计，而P70 Art则采用了与上一代P60 Art相似的不规则形状设计。这样的外观是否好看见仁见智，但辨识度绝对拉满。