Issue: ElevationAboveStreamEuclidean Produces Empty TIFF

[Heljor]

Issue: ElevationAboveStreamEuclidean Produces Empty TIFF

When executing the elevation_above_stream_euclidean function, the resulting DEM outputs an empty .tif file.

Steps to Reproduce:

# Modify no data value
wbt.modify_no_data_value(
    i="dem.tif",
    new_value="0.0"


)

# Breach single cell pits
wbt.breach_single_cell_pits(
    dem="dem.tif",
    output="breach_single_cell_pits.tif"
)

# Fill depressions using Wang and Liu's method
wbt.fill_depressions_wang_and_liu(
    dem="breach_single_cell_pits.tif",
    output="filldepress.tif",
    fix_flats=True,
    flat_increment=None
)

# D8 pointer
wbt.d8_pointer(
    dem="filldepress.tif",
    output="D8Pointer.tif",
    esri_pntr=True
)

# D8 flow accumulation
wbt.d8_flow_accumulation(
    i="D8Pointer.tif",
    output="D8FlowAccumulation.tif",
    out_type="specific contributing area",
    log=False,
    clip=False,
    pntr=True,
    esri_pntr=True
)

# Extract stream
wbt.extract_streams(
    flow_accum="D8FlowAccumulation.tif",
    output="extractstream.tif",
    threshold="1000",
    zero_background=False
)

# Elevation above stream Euclidean
wbt.elevation_above_stream_euclidean(
    dem="breach_single_cell_pits.tif",
    streams="extractstream.tif",
    output="elevation_above_stream_euclidean.tif"
)

[Atreyagaurav]

I don't know why this has been ignored for so long, this was due to the code masking the output value from the streams raster's no_data instead of input DEM's.

Here is my original issue with a patch: #245

But you'd have to compile it. And looks like the code has been changed, so you can't apply that patch to the new version automatically. Hopefully it'll be fixed.