Skip to content

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

wkt raster

parent 1a97be40
Hide whitespace changes
Inline Side-by-side
Showing with 1 addition and 1 deletion
geosentinel/rasterWkt.py
View file @ da6e1968
...@@ -118,7 +118,7 @@ def rasterWkt(wkt, inputfile, outputfile): ...@@ -118,7 +118,7 @@ def rasterWkt(wkt, inputfile, outputfile):
np_polygon = polygonRaster.GetRasterBand(1).ReadAsArray(0,0,cols,rows) np_polygon = polygonRaster.GetRasterBand(1).ReadAsArray(0,0,cols,rows)
# np_polygon = getPolygonArray(intput, wkt) # np_polygon = getPolygonArray(intput, wkt)
for b in range(1, output.RasterCount+1): for b in range(1, output.RasterCount+1):
np_input = np_polygon & intput.GetRasterBand(b).ReadAsArray(0,0,cols,rows) np_input = intput.GetRasterBand(b).ReadAsArray(0,0,cols,rows)
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