Frank, folks, I found an error with gdalwarp. It appears that when the -rcs resampling option is used in combination with -srcnodata and -dstnodata, the nodata pixels still get interpolated. When another resampling option is used, the nodata pixels are (correctly) retained as such in the output image. Regards, Vincent.