I'm trying to read ROIs from Landsat-8 overview using GDAL /vsicurl/ , but the image is ungeoreferenced. I have geotransform from full size image, which I can copy. Is there a way to pass geotransform to gdalwarp to handle ungeoreferenced image as georeferenced?
Here is the command I want to execute from Python:
ds_source = 'http://landsat-pds.s3.amazonaws.com/L8/044/034/LC80440342015224LGN00/LC80440342015224LGN00_B3.TIF.ovr' # here I use -te_srs,-te to crop ROI; data should be taken from 2nd OVR level warp_options = ' -of MEM -r cubic -te_srs "+proj=longlat +datum=WGS84 +no_defs " -te -123.75 31.952162238 -112.5 40.9798980696 -dstalpha -ovr 2' output_dataset = gdal.Warp('', '/vsicurl/%s' % ds_source, options=warp_options)
Lines above produce error:
ERROR 1: -te_srs ignored since none of -t_srs, -s_srs is specified or the input dataset has no projection.