VRT ComputeStatistics() fails on VRTs that have rescaled values

[lpinner]

Possibly related to #5468? Is also apparent in GDAL 1.10

gdal_translate -of vrt -ot byte -scale -1 1 0 255 gdal/trunk/autotest/gcore/data/float32.tif test.vrt
gdalinfo -stats test.vrt
<snip>
Band 1 Block=20x20 Type=Byte, ColorInterp=Gray
ERROR 1: VRTSourcedRasterBand::IRasterIO() called recursively on the same band. It looks like the VRT is referencing itself.
ERROR 1: test.vrt, band 1: IReadBlock failed at X offset 0, Y offset 0
ERROR 1: test.vrt, band 1: Failed to compute statistics, no valid pixels found in sampling.

[lpinner]

I've tried this with a number of Float32 netcdfs and tifs and get the same errors. I originally came across the issue trying to generate stretched NDVI images (data range -1 to 1) for display.

[Even Rouault]

Yes, this is the same issue as #5468.

[lpinner]

Thanks Even, I've confirmed ComputeStatistics is working for rescaled VRTs using a trunk build from Tamas' GISInternals site.