通过GDAL创建与写Shape文件时的wkt格式

通过GDAL创建与写Shape文件时的wkt格式

GDAL是一个符合OGC标准并且开源的,可以用于地理信息开发的空间数据转换库，利用它可以方便地对一些标准地理信息数据格式进行加载、修改及更新。

使用GDAL创建shp文件时，需要定义wkt字符串，各种基本要素类型的格式如下：

POINT(6 10)
LINESTRING(3 4,10 50,20 25)
POLYGON((1 1,5 1,5 5,1 5,1 1),(2 2,2 3,3 3,3 2,2 2))
MULTIPOINT(3.5 5.6, 4.8 10.5)
MULTILINESTRING((3 4,10 50,20 25),(-5 -8,-10 -8,-15 -4))
MULTIPOLYGON(((1 1,5 1,5 5,1 5,1 1),(2 2,2 3,3 3,3 2,2 2)),((6 3,9 2,9 4,6 3)))
GEOMETRYCOLLECTION(POINT(4 6),LINESTRING(4 6,7 10))
POINT ZM (1 1 5 60)
POINT M (1 1 80)

定义的生成要素的类型保持一致，最后是清空缓存，这里多说一句，字符串语法与postgis等开源gis一致，都遵循ogc国际标准。

wkt = ‘LINESTRING(3 4,10 50,20 25)’
line = ogr.CreateGeometryFromWkt(wkt)
feature.SetGeometry(line)
layer.CreateFeature(feature)

feature = None
data_source = None