drawing polygon

0bd35fc9 · Mario Chirinos Colunga · 6aefb8fe · 0bd35fc9
Commit 0bd35fc9 authored Aug 20, 2018 by Mario Chirinos Colunga 💬
Hide whitespace changes
Inline Side-by-side

Showing with 1 addition and 1 deletion

rasterWkt.py geosentinel/rasterWkt.py +1 -1

No files found.
--- a/geosentinel/rasterWkt.py
+++ b/geosentinel/rasterWkt.py
@@ -76,7 +76,7 @@ def rasterWkt(wkt, inputfile, outputfile):
 	#Draw output image
 	np_polygon = polygonRaster.GetRasterBand(1).ReadAsArray(0,0,cols,rows)
 	for b in range(1, output.RasterCount+1):
-		np_input = np_polygon | intput.GetRasterBand(b).ReadAsArray(0,0,cols,rows)
+		np_input = np_polygon & intput.GetRasterBand(b).ReadAsArray(0,0,cols,rows)
 		output.GetRasterBand(b).SetNoDataValue(-99)
 		output.GetRasterBand(b).WriteArray(np_input)
 		output.GetRasterBand(b).FlushCache()