filter_band process seems missing

[danielFlemstrom]

If i get it right, the following should be possible to write in the python client:
rgb_cube = cube.filter_bands(['oa04', 'oa06', 'oa09']) and seems be accepted in the openeo-python-client,as

However I cannot find any traces of that process in here. Is someone working on it, or should we proceed in another way?

All the best
/Daniel

[LukeWeidenwalker]

Afaik this process simply hasn't been implemented yet on our end! I don't have the exact historic context, but I presume that it just hasn't been required for any of the usecases we've implemented so far. To my knowledge, this isn't on our roadmap for the near future, but we're more than happy to accept any PRs that take a crack at this!

@ValentinaHutter please correct me if you have been working on this at any point!

edit: see this issue with missing processes that are on the roadmap

[danielFlemstrom]

I cannot see the roadmap though, maybe i do not have permissions ?

[LukeWeidenwalker]

Ah, apologies - this is indeed on a private openEOPlatform repo!

[clausmichele]

It wouldn't be difficult to implement it and it is indeed an useful process!
@danielFlemstrom, if you'd like to try, you can do something like this https://github.com/SARScripts/openeo_odc_driver/blob/ad6546feb086d39ae7ef62ad22fd6096336534a4/openeo_odc_driver.py#L898

[ValentinaHutter]

I now added this process to our development backend, so if you want to test it you can connect to that using

eodc_backend = "https://openeo-dev.eodc.eu/v1.0"
conn = openeo.connect(eodc_backend).authenticate_oidc(provider_id="egi")

a simple example would be

spatial_extent = {'west':10.454955,'east':10.573729,'south':46.102185,'north':46.213657,'crs':'EPSG:4326'}
data = conn.load_collection('boa_sentinel_2', spatial_extent=spatial_extent, bands=["B02", "B03", "B04", "B05", "B08"] , temporal_extent=["2019-05-01", "2019-06-01"] )
data_bands = data.filter_bands(bands=["B02", "B03"])
data_saved = data_bands.save_result(format="NetCDF")
job = data_saved.create_job()
job.start_job()

I want to update our live backend this week as well, so I will inform you once you can use it there, too.

[clausmichele]

@ValentinaHutter I've tried you code but I can't open the result files. Additionally the files do not have the .nc extensions.

results = job.get_results()

results.download_files('eodc_test')

[PosixPath('eodc_test/out_2019_05_02T12_00_00_EU010M_E047N014T1'),
 PosixPath('eodc_test/out_2019_05_02T12_00_00_EU010M_E048N014T1'),
 PosixPath('eodc_test/out_2019_05_09T12_00_00_EU010M_E047N014T1'),
 PosixPath('eodc_test/out_2019_05_09T12_00_00_EU010M_E048N014T1'),
 PosixPath('eodc_test/out_2019_05_14T12_00_00_EU010M_E047N014T1'),
 PosixPath('eodc_test/out_2019_05_14T12_00_00_EU010M_E048N014T1'),
 PosixPath('eodc_test/out_2019_05_24T12_00_00_EU010M_E047N014T1'),
 PosixPath('eodc_test/out_2019_05_24T12_00_00_EU010M_E048N014T1'),
 PosixPath('eodc_test/job-results.json')]

import xarray as xr

xr.open_dataset('eodc_test/eodc_test/out_2019_05_02T12_00_00_EU010M_E047N014T1')

---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
Input In [10], in <cell line: 2>()
      1 import xarray as xr
----> 2 xr.open_dataset('eodc_test/eodc_test/out_2019_05_02T12_00_00_EU010M_E047N014T1')

File ~/anaconda3/envs/eurac_env/lib/python3.8/site-packages/xarray/backends/api.py:479, in open_dataset(filename_or_obj, engine, chunks, cache, decode_cf, mask_and_scale, decode_times, decode_timedelta, use_cftime, concat_characters, decode_coords, drop_variables, backend_kwargs, *args, **kwargs)
    476     kwargs.update(backend_kwargs)
    478 if engine is None:
--> 479     engine = plugins.guess_engine(filename_or_obj)
    481 backend = plugins.get_backend(engine)
    483 decoders = _resolve_decoders_kwargs(
    484     decode_cf,
    485     open_backend_dataset_parameters=backend.open_dataset_parameters,
   (...)
    491     decode_coords=decode_coords,
    492 )

File ~/anaconda3/envs/eurac_env/lib/python3.8/site-packages/xarray/backends/plugins.py:148, in guess_engine(store_spec)
    140 else:
    141     error_msg = (
    142         "found the following matches with the input file in xarray's IO "
    143         f"backends: {compatible_engines}. But their dependencies may not be installed, see:\n"
    144         "https://docs.xarray.dev/en/stable/user-guide/io.html \n"
    145         "https://docs.xarray.dev/en/stable/getting-started-guide/installing.html"
    146     )
--> 148 raise ValueError(error_msg)

ValueError: did not find a match in any of xarray's currently installed IO backends ['netcdf4', 'h5netcdf', 'scipy', 'rasterio']. Consider explicitly selecting one of the installed engines via the ``engine`` parameter, or installing additional IO dependencies, see:
https://docs.xarray.dev/en/stable/getting-started-guide/installing.html
https://docs.xarray.dev/en/stable/user-guide/io.html

gdal also fails:

!gdalinfo eodc_test/eodc_test/out_2019_05_02T12_00_00_EU010M_E047N014T1

ERROR 4: eodc_test/eodc_test/out_2019_05_02T12_00_00_EU010M_E047N014T1: No such file or directory
gdalinfo failed - unable to open 'eodc_test/eodc_test/out_2019_05_02T12_00_00_EU010M_E047N014T1'.

[ValentinaHutter]

Are you sure you want to use this path as it is?

xr.open_dataset('eodc_test/eodc_test/out_2019_05_02T12_00_00_EU010M_E047N014T1')

to me it looks like the files are stored in this path:

'eodc_test/out_2019_05_02T12_00_00_EU010M_E047N014T1'

No such file or directory is an error that is normally thrown when the path is incorrect.

[clausmichele]

@ValentinaHutter sorry, such a stupid mistake! yes the files are readable, they are missing the file extension tough.