wkt raster

geosentinel/rasterWkt.py

@@ -116,10 +116,10 @@ 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 = inputImage.GetRasterBand(b).ReadAsArray(0,0,cols,rows)
+		np_input = np_polygon & inputImage.GetRasterBand(b).ReadAsArray(0,0,cols,rows)
 		output.GetRasterBand(b).SetNoDataValue(-99)
 		output.GetRasterBand(b).WriteArray(np_input)
 		output.GetRasterBand(b).FlushCache()