rgb l2

geosentinel/rasterWkt.py

@@ -49,7 +49,23 @@ def rasterWkt(wkt, inputfile, outputfile):
 	sr.ImportFromWkt(wkt)
 	output.SetProjection(sr.ExportToWkt())
 
-	gdal.RasterizeLayer(output, [1], createLayer(wkt), burn_values=[0])
+
+	# shapefile
+	driver = ogr.GetDriverByName("ESRI Shapefile")
+	data_source = driver.CreateDataSource("myShape.shp")
+
+	srs = osr.SpatialReference()
+	srs.ImportFromEPSG(4326)
+	layer = data_source.CreateLayer(layerName, srs, ogr.wkbPolygon)
+	geometry = ogr.CreateGeometryFromWkt(wkt)
+	feature = ogr.Feature(layer.GetLayerDefn())
+	feature.SetGeometry(geometry)
+	layer.CreateFeature(feature)
+
+	gdal.RasterizeLayer(output, [1], layer, burn_values=[0])
+	
+	feature = None
+	data_source = None
 
 	output = None