Frank, is it on purpose that gdal_merge.py does not set the nodata value of the resulting file? I 'fixed' this for myself by adding "t_band.SetNoDataValue(nodata)" in raster_copy_with_nodata (line 141). Is this appropriate? Vincent.