python 矢量图形库_Python与开源GIS：使用OGR定义矢量数据的字段

GIS

数据除了图形要素之外，更重要的就是属性数据了。关于

GIS

数据的属性，完全可

以用关系型数据库的观点来看待。

只不过在这个关系数据库中，

专门一列索引值，

并数据库

中的记录与图形对象对应起来，而且，这种对应还是“一一对应”。

我们可以给那个矩形添加属性数据。

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

# -*- coding:utf-8 -*-

from osgeo import ogr

import os,math

driver = ogr.GetDriverByName("ESRI Shapefile")

extfile = 'rect_field_demo.shp'

if os.access( extfile, os.F_OK ):

driver.DeleteDataSource( extfile )

extent = [400, 1100, 300, 600]

newds = driver.CreateDataSource(extfile)

layernew = newds.CreateLayer('rect',None,ogr.wkbPolygon)

fieldcnstr = ogr.FieldDefn("fd",ogr.OFTString)

fieldcnstr.SetWidth(32)

layernew.CreateField(fieldcnstr)

fieldf = ogr.FieldDefn("f",ogr.OFTReal)

layernew.CreateField(fieldf)

wkt = 'POLYGON ((%f %f,%f %f,%f %f,%f %f,%f %f))' % (extent[0],extent[3],

extent[1],extent[3], extent[1],extent[2],