Skip to content

GeoSentinel
GeoSentinel
Commit c6d79956 authored by Mario Chirinos Colunga's avatar Mario Chirinos Colunga 💬
Browse files
Options

drawing polygon

parent 21b0f8fb
Show whitespace changes
Inline Side-by-side
Showing with 1 addition and 1 deletion
geosentinel/rasterWkt.py
View file @ c6d79956
...@@ -74,7 +74,7 @@ def rasterWkt(wkt, inputfile, outputfile): ...@@ -74,7 +74,7 @@ def rasterWkt(wkt, inputfile, outputfile):
#Draw output image #Draw output image
np_polygon = polygonRaster.GetRasterBand(1).ReadAsArray(0,0,cols,rows) np_polygon = polygonRaster.GetRasterBand(1).ReadAsArray(0,0,cols,rows)
for b in range(1, output.RasterCount+1): for b in range(1, output.RasterCount+1):
np_input = intput.GetRasterBand(b).ReadAsArray(0,0,cols,rows) | polygonRaster np_input = intput.GetRasterBand(b).ReadAsArray(0,0,cols,rows) | np_polygon
output.GetRasterBand(b).SetNoDataValue(-99) output.GetRasterBand(b).SetNoDataValue(-99)
output.GetRasterBand(b).WriteArray(np_input) output.GetRasterBand(b).WriteArray(np_input)
output.GetRasterBand(b).FlushCache() output.GetRasterBand(b).FlushCache()
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment