Hi there, Would be possible to write a 'VRT Pixel Function' [1] in Python? Probably not. But wouldn't that be nice? Ivan [1] - http://www.gdal.org/gdal_vrttut.html