netCDF GDAL NoData Value=32769 after linear_scale_range(0,65335,0,65335)

[VincentVerelst]

When applying the linear_scale_range process to convert dtype to UInt16 the resulting netCDF has a NoData Value of 32769 instead of the expected 65335, using gdalinfo.

An example job is CDSE j-25010912043543edbb291a8946dbb1e1

This issue can be reproduced with the following code:

import openeo

c = openeo.connect('openeo.dataspace.copernicus.eu').authenticate_oidc()

bbox = {'north': 52.5, 'south': 52.45, 'east': 13.4, 'west': 13.35}

s2 = c.load_collection(
    collection_id="SENTINEL2_L2A",
    spatial_extent=bbox,
    temporal_extent=["2020-01-01", "2020-01-10"],
    bands = ["B04"]
)
s2 = s2.linear_scale_range(0, 65335, 0, 65335)

job = s2.execute_batch(
    out_format="netcdf",
    title="test title",
)