<div dir="ltr"><div>I'm using GDAL v2.1.3</div><div><br></div>I have a tiff file with this info:<div><div><font face="monospace, monospace">Driver: GTiff/GeoTIFF</font></div><div><font face="monospace, monospace">Files: Chlorofyl-index.tif</font></div><div><font face="monospace, monospace"> Chlorofyl-index.tif.aux.xml</font></div><div><font face="monospace, monospace">Size is 5140, 6065</font></div><div><font face="monospace, monospace">Coordinate System is:</font></div><div><font face="monospace, monospace">PROJCS["WGS 84 / UTM zone 32N",</font></div><div><font face="monospace, monospace"> GEOGCS["WGS 84",</font></div><div><font face="monospace, monospace"> DATUM["WGS_1984",</font></div><div><font face="monospace, monospace"> SPHEROID["WGS 84",6378137,298.257223563,</font></div><div><font face="monospace, monospace"> AUTHORITY["EPSG","7030"]],</font></div><div><font face="monospace, monospace"> AUTHORITY["EPSG","6326"]],</font></div><div><font face="monospace, monospace"> PRIMEM["Greenwich",0,</font></div><div><font face="monospace, monospace"> AUTHORITY["EPSG","8901"]],</font></div><div><font face="monospace, monospace"> UNIT["degree",0.0174532925199433,</font></div><div><font face="monospace, monospace"> AUTHORITY["EPSG","9122"]],</font></div><div><font face="monospace, monospace"> AUTHORITY["EPSG","4326"]],</font></div><div><font face="monospace, monospace"> PROJECTION["Transverse_Mercator"],</font></div><div><font face="monospace, monospace"> PARAMETER["latitude_of_origin",0],</font></div><div><font face="monospace, monospace"> PARAMETER["central_meridian",9],</font></div><div><font face="monospace, monospace"> PARAMETER["scale_factor",0.9996],</font></div><div><font face="monospace, monospace"> PARAMETER["false_easting",500000],</font></div><div><font face="monospace, monospace"> PARAMETER["false_northing",0],</font></div><div><font face="monospace, monospace"> UNIT["metre",1,</font></div><div><font face="monospace, monospace"> AUTHORITY["EPSG","9001"]],</font></div><div><font face="monospace, monospace"> AXIS["Easting",EAST],</font></div><div><font face="monospace, monospace"> AXIS["Northing",NORTH],</font></div><div><font face="monospace, monospace"> AUTHORITY["EPSG","32632"]]</font></div><div><font face="monospace, monospace">Origin = (360810.191610000040000,5860662.620930000200000)</font></div><div><font face="monospace, monospace">Pixel Size = (0.125010000000000,-0.125010000000000)</font></div><div><font face="monospace, monospace">Metadata:</font></div><div><font face="monospace, monospace"> AREA_OR_POINT=Area</font></div><div><span style="font-family:monospace,monospace">Image Structure Metadata:</span><br></div><div><font face="monospace, monospace"> COMPRESSION=LZW</font></div><div><font face="monospace, monospace"> INTERLEAVE=BAND</font></div><div><font face="monospace, monospace">Corner Coordinates:</font></div><div><font face="monospace, monospace">Upper Left ( 360810.192, 5860662.621) ( 6d55'54.04"E, 52d52'39.50"N)</font></div><div><font face="monospace, monospace">Lower Left ( 360810.192, 5859904.435) ( 6d55'55.21"E, 52d52'14.98"N)</font></div><div><font face="monospace, monospace">Upper Right ( 361452.743, 5860662.621) ( 6d56'28.39"E, 52d52'40.10"N)</font></div><div><font face="monospace, monospace">Lower Right ( 361452.743, 5859904.435) ( 6d56'29.55"E, 52d52'15.58"N)</font></div><div><font face="monospace, monospace">Center ( 361131.467, 5860283.528) ( 6d56'11.80"E, 52d52'27.54"N)</font></div><div><font face="monospace, monospace">Band 1 Block=256x256 Type=Float32, ColorInterp=Gray</font></div><div><font face="monospace, monospace"> Min=-0.177 Max=33.644 </font></div><div><font face="monospace, monospace"> Minimum=-0.177, Maximum=33.644, Mean=18.605, StdDev=7.894</font></div><div><font face="monospace, monospace"> NoData Value=-10000</font></div><div><font face="monospace, monospace"> Metadata:</font></div><div><font face="monospace, monospace"> STATISTICS_MAXIMUM=33.644180297852</font></div><div><font face="monospace, monospace"> STATISTICS_MEAN=18.604999564386</font></div><div><font face="monospace, monospace"> STATISTICS_MINIMUM=-0.17713606357574</font></div><div><font face="monospace, monospace"> STATISTICS_STDDEV=7.8939827089249</font></div></div><div><br></div><div>I cut this file with a border shapefile and save it to a VRT file using this command:</div><div><font face="monospace, monospace">gdalwarp -of VRT -cutline border.shp -crop_to_cutline -dstnodata 0 -overwrite Chlorofyl-index.tif Chlorofyl-index-clipped.vrt</font><br></div><div><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>I open this vrt-file in MapWindow5 and in QGis and it looks as expected.</div><div><br></div><div>I now add external overviews to the tiff file using:</div><div><font face="monospace, monospace">gdaladdo -ro -r nearest Chlorofyl-index.tif 2 4 8 16</font><br></div><div dir="ltr"><br>
I reopen the vrt file and both MapWindow5 as QGis show a complete different part of the tiff.</div><div dir="ltr">In QGis it is much smaller and outside the original border and with MW5 it is bigger and less high than the border. Either way they are wrong.</div><div dir="ltr"><br></div><div dir="ltr">When I delete the .ovr file and reopen the vrt-file all is well again.</div><div dir="ltr">I already tried running my GDALWarp command after creating the overviews but that doesn't matter.</div><div dir="ltr"><br></div><div>Creating a tiff-file instead of a vrt-file results in a correct result.</div><div><br></div><div>Is this a known and expected behavior of the VRT-format?</div><div>I really want to keep using the vrt-file since it is only 700kb instead of the tiff-file which is 63MB and I will be doing more calculations later on, which also result in new files.</div><div><br></div><div>If needed I can send my input files.</div><div><br></div><div>Thanks,</div><div dir="ltr"><br>
Paul<br>
</div></div></div></div></div></div>
</div></div>