创建Shapefile文件并写入数据的例子
基本思路
使用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中的几何数据。
fromosgeoimportogr fromosgeoimportosr importjson importos os.environ['SHAPE_ENCODING']="utf-8" withopen('China.json')asf: china=json.load(f) #创建DataSource driver=ogr.GetDriverByName('ESRIShapefile') 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 forfinchina['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) delfeature ds.FlushCache() delds
以上这篇创建Shapefile文件并写入数据的例子就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持毛票票。
声明:本文内容来源于网络,版权归原作者所有,内容由互联网用户自发贡献自行上传,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任。如果您发现有涉嫌版权的内容,欢迎发送邮件至:czq8825#qq.com(发邮件时,请将#更换为@)进行举报,并提供相关证据,一经查实,本站将立刻删除涉嫌侵权内容。