<div dir="ltr"><div dir="ltr">Sorry for the multiple messages, but one more data point.<div><br></div><div>If the VRT is composed of more than one image, everything seems to work.  I can literally do a "gdalbuildvrt out.vrt in.tif in.tif" and "gdalinfo out.vrt" on the multiple image VRT does not report the existence of overviews (it does for the single image VRT).  Extracting the nodata mask on out.vrt works (e.g. its not returning the first band of the image), although I'm not sure if its just doing a nearest neighbor resampling of the full resolution mask or if its implicitly cracking open the underlying geotiff's overview level.  I suppose it's the former or I wouldn't see a difference in behavior.</div><div><br></div><div>Is there some optimization in the VRT driver for when their is just a single image in the source?  Even's blog post (<a href="http://erouault.blogspot.com/2014/10/warping-overviews-and-warped-overviews.html">http://erouault.blogspot.com/2014/10/warping-overviews-and-warped-overviews.html</a>) on warped VRTs (see the last section of the post) makes me think its going straight to the overviews in the single image case, but this is turned off when more than one image composes the VRT.  Presumably the handling of masks in the former case is where the problem lies.</div><div><br></div><div>Thanks for any clues!</div><div>Patrick</div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, May 15, 2019 at 11:17 AM Patrick Young <<a href="mailto:patrick.mckendree.young@gmail.com">patrick.mckendree.young@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi all,<br><div><br></div><div>I'm experiencing some strange behavior when extracting a mask band using gdal_translate from a VRT vs the source GTiff at a reduced resolution.  I have a 3 band, 8 bit JPEG compressed (YCbCR) RGB geotiff with an internal no data mask and with internal overviews.  Making a VRT of this and running</div><div><br></div><blockquote style="margin:0px 0px 0px 40px;border:none;padding:0px"><div>gdal_translate -b mask in.vrt mask.tif</div></blockquote><div><br></div><div>Returns the expected no data mask.  </div><div><br></div><div>However, if I add in "-outsize 10% 0", I instead get the first band back as the mask tif rather than the mask band.  Any -outsize setting has the same effect.</div><div><br></div><div>Repeating this with the source geotiff rather than the vrt works as expected.</div><div><br></div><div>Interestingly, running</div><div><br></div><blockquote style="margin:0px 0px 0px 40px;border:none;padding:0px"><div>gdalinfo -oo OVERVIEW_LEVEL=1 in.vrt</div><div><br></div></blockquote>reports no mask flags set at all (they are reported without the -oo option), but running on the original geotiff rather than the vrt does report the mask flags.    <div><br></div><div>Anyone have thoughts on what the VRT driver is doing? </div><div><br></div><div>Thanks,</div><div>Patrick<br><div><br></div><div><br></div><div><br></div><div><br></div></div></div>
</blockquote></div>