wkt raster

280726a1 · Mario Chirinos Colunga · 4e1631d2 · 280726a1
Commit 280726a1 authored Oct 09, 2018 by Mario Chirinos Colunga 💬
Hide whitespace changes
Inline Side-by-side

Showing with 6 additions and 6 deletions

rasterWkt.py geosentinel/rasterWkt.py +6 -6

No files found.
--- a/geosentinel/rasterWkt.py
+++ b/geosentinel/rasterWkt.py
@@ -106,14 +106,14 @@ def rasterWkt(wkt, inputfile, outputfile):
 	srs = osr.SpatialReference()
 	srs.ImportFromEPSG(4326)
 	layer = data_source.GetLayer()
-#	layer = data_source.CreateLayer("wkt", srs, ogr.wkbPolygon)
-#	geometry = ogr.CreateGeometryFromWkt(wkt)
-#	feature = ogr.Feature(layer.GetLayerDefn())
+	layer = data_source.CreateLayer("wkt", srs, ogr.wkbPolygon)
+	geometry = ogr.CreateGeometryFromWkt(wkt)
+	feature = ogr.Feature(layer.GetLayerDefn())

 	print("features")
-#	feature.SetStyleString("PEN(c:#FF0000,w:5px);")
-#	feature.SetGeometry(geometry)
-#	layer.CreateFeature(feature)
+	feature.SetStyleString("PEN(c:#FF0000,w:5px);")
+	feature.SetGeometry(geometry)
+	layer.CreateFeature(feature)

 	gdal.RasterizeLayer(polygonRaster, [1], layer, burn_values=[255])
 	polygonRaster.GetRasterBand(1).FlushCache()