The output backscatter applies two nodata values: 0 (black area) and nan (white area):
0 is not defined as nodata value in the metadata.
Temporary workaround: set a consistent nodata value in post:
ds = gdal.Open("gt.tif", GA_Update)
ar = ds.ReadAsArray()
ar[~np.isfinite(ar)] = -9999
ar[ar == 0] = -9999
band = ds.GetRasterBand(1)
band.SetNoDataValue(-9999)
band.WriteArray(ar)
band.FlushCache()
band = None
ds = None
The output backscatter applies two nodata values: 0 (black area) and nan (white area):
0 is not defined as nodata value in the metadata.
Temporary workaround: set a consistent nodata value in post: