Commit cab6a89b authored by Manuel Grizonnet's avatar Manuel Grizonnet

ENH: replace build vrt with gdal api

parent d1acf41d
......@@ -235,8 +235,9 @@ class snow_detector :
#build vrt
print "building bands vrt"
self.img=op.join(self.path_tmp, "lis.vrt")
call_subprocess(["gdalbuildvrt","-separate", self.img, sb_path_resampled, rb_path_resampled, gb_path_resampled])
gdal.BuildVRT(self.img,[sb_path_resampled, rb_path_resampled, gb_path_resampled], separate=True)
#Set bands parameters
self.nGreen=3
self.nRed=2
......
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