python将wkt解析为shp数据
时间: 2024-03-18 11:44:41 浏览: 134
解析wkt数据
要将WKT解析为SHP数据,可以使用Python中的GDAL库。以下是一个简单的示例代码:
```python
from osgeo import ogr
# WKT数据
wkt = 'POLYGON ((30 10, 40 40, 20 40, 10 20, 30 10))'
# 创建空的矢量数据源
driver = ogr.GetDriverByName('ESRI Shapefile')
data_source = driver.CreateDataSource('output.shp')
# 创建矢量图层
layer = data_source.CreateLayer('polygon', geom_type=ogr.wkbPolygon)
# 添加字段
field_defn = ogr.FieldDefn('id', ogr.OFTInteger)
layer.CreateField(field_defn)
# 创建要素
feature_defn = layer.GetLayerDefn()
feature = ogr.Feature(feature_defn)
feature.SetGeometry(ogr.CreateGeometryFromWkt(wkt))
feature.SetField('id', 1)
# 将要素添加到图层中
layer.CreateFeature(feature)
# 释放资源
feature = None
data_source = None
```
在上面的代码中,我们首先定义了一个WKT字符串。然后,我们使用GDAL库创建了一个空的矢量数据源,然后创建了一个矢量图层,并添加了一个整数类型的字段。接下来,我们使用WKT字符串创建一个要素,并将其添加到图层中。最后,我们释放了资源并保存了SHP文件。
阅读全文