wkt raster

geosentinel/rasterWkt.py

@@ -118,7 +118,7 @@ def rasterWkt(wkt, inputfile, outputfile):
 
 #	Draw output image
 	np_polygon = polygonRaster.GetRasterBand(1).ReadAsArray(0,0,cols,rows)
-#	np_polygon = getPolygonArray(inputImage, wkt)
+	np_polygon = getPolygonArray(inputImage, wkt)
 	print("Setup")
 	for b in range(1, output.RasterCount+1):
 		np_input = np_polygon & inputImage.GetRasterBand(b).ReadAsArray(0,0,cols,rows)
@@ -131,9 +131,9 @@ def rasterWkt(wkt, inputfile, outputfile):
 
 	print("end")
 	#Close files
-#	feature = None
-#	data_source = None
-#	polygonRaster = None
+	feature = None
+	data_source = None
+	polygonRaster = None
 	output = None
 	inputImage = None