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数据库的一个表。
关注公众号
获取免费资源
Copyright © Since 2014.
开源地理空间基金会中文分会
吉ICP备05002032号
Powered by TorCMS