I'm new using python and I'm trying to create a Virtual Raster mosaicing 3 raster images of Landsat in a time series. So I'm trying with gdal.BuildVRT but it doesn't work? I'll show the code:
from osgeo import gdal import os import glob #SET PATH path = 'D:\Imagery\L8 OLI TIRS' x = len(files))/3 #IN THE FILE I HAVE 3 TILES FOR EACH DATE """ example with only one date NDVI """ #subset route: in_subRoute = files+'-'+files[0+x]+'-'+files[0+x*2] #string, every 'files' equals to a directory of a tile. subRoute = in_subRoute.split('-') #to list Dirmosaic = "D:\Results\mosaic.vrt" gdal.BuildVRT(Dirmosaic, subRoute)
but doesn't happen anything
What am I doing wrong?