Create VRT layer from wkt

[ustroetz]

Do we have a recipe that explains how to create a VRT layer?

If not, I have some ideas, but couldn't get it to work. Maybe somebody can help, then I'll add it as a recipe.

from osgeo import ogr

poly = ogr.CreateGeometryFromWkt('POLYGON((-0.95245901639344266 0.51475409836065578,-0.13606557377049189 0.4459016393442623,-0.45737704918032795 0.04918032786885251,-0.45737704918032795 0.04918032786885251,-0.95245901639344266 0.51475409836065578))')

# Create the output shapefile
vrtDriver = ogr.GetDriverByName("VRT")
outDataSource = vrtDriver.CreateDataSource('test.vrt')
outLayer = outDataSource.CreateLayer('test', geom_type=ogr.wkbPolygon )

# Create the feature and set values
featureDefn = outLayer.GetLayerDefn()
outFeature = ogr.Feature(featureDefn)
outFeature.SetGeometry(poly)
outLayer.CreateFeature(outFeature)