创建数据表的WKT

创建数据表的WKT
发布日期: 2016-10-06 更新日期: 1970-01-01 编辑:yubiao 浏览次数: 3325

 import ogr, osr

database = 'test'
usr = 'postgres'
pw = ''
table = 'test'

wkt = "POINT (1120351.5712494177 741921.4223245403)"
point = ogr.CreateGeometryFromWkt(wkt)

connectionString = "PG:dbname='%s' user='%s' password='%s'" % (database,usr,pw)
ogrds = ogr.Open(connectionString)

srs = osr.SpatialReference()
srs.ImportFromEPSG(4326)

layer = ogrds.CreateLayer(table, srs, ogr.wkbPoint, ['OVERWRITE=YES'] )

layerDefn = layer.GetLayerDefn()

feature = ogr.Feature(layerDefn)
feature.SetGeometry(point)

layer.StartTransaction()
layer.CreateFeature(feature)
layer.CommitTransaction()
说明:

这个方法创建现有PostGIS数据库的一个表。