Saving back ds after modifications to the same file format as provided by AVIRIS

[ohadshapira]

Hi HyperCoast team,

First of all, thank you!
Such an amazing project and it has been very useful in the past few weeks. Have been using a different tool and got a major upgrade when I made the transition to yours.

I have a question:
I'm trying to modify a hyperspectral image but couldn't figure out how can I save it back in the same format.
For example, How can I save the new ds in the same file format as provided by AVIRIS?

`
import hypercoast

filepath = "ang20210401t150456_rfl_v2z1"
ds = hypercoast.read_aviris(filepath)
ds = ds.coarsen(y=4, x=4, boundary='pad').mean() # downsample ds

How to save back ds?
`

Thank you,
Ohad

[giswqs]

Make sure the DataArray has crs and transform. Check out the link below on how to set crs and transform
https://corteva.github.io/rioxarray/html/getting_started/crs_management.html#API-Documentation

[giswqs]

See this example

HyperCoast/hypercoast/pace.py

Line 669 in f3a7863

data.rio.write_crs("EPSG:4326", inplace=True)

[ohadshapira]

Having an issue with the file formatting. let me write it in a different framing:

I'm trying to read, save and read new file:

filepath = "ang20210401t150456_rfl_v2z1"

ds = hypercoast.read_aviris(filepath)

hypercoast.aviris_to_image(ds,output="test",cellsize=12,crs="EPSG:3857")

ds = hypercoast.read_aviris("test") # trying to read saved file

I'm getting this error from the last command:

Traceback (most recent call last):
  File "C:\Users\Home\PycharmProjects\hsi\.venv\Lib\site-packages\xarray\core\dataset.py", line 1477, in _construct_dataarray
    variable = self._variables[name]
               ~~~~~~~~~~~~~~~^^^^^^
KeyError: 'wavelength'
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
  File "C:\Users\Home\PycharmProjects\hsi\.venv\Lib\site-packages\xarray\core\dataset.py", line 1576, in __getitem__
    return self._construct_dataarray(key)
           ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "C:\Users\Home\PycharmProjects\hsi\.venv\Lib\site-packages\xarray\core\dataset.py", line 1479, in _construct_dataarray
    _, name, variable = _get_virtual_variable(self._variables, name, self.sizes)
                        ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "C:\Users\Home\PycharmProjects\hsi\.venv\Lib\site-packages\xarray\core\dataset.py", line 213, in _get_virtual_variable
    raise KeyError(key)
KeyError: 'wavelength'
The above exception was the direct cause of the following exception:
Traceback (most recent call last):
  File "C:\Users\Home\PycharmProjects\hsi\.venv\Lib\site-packages\IPython\core\interactiveshell.py", line 3577, in run_code
    exec(code_obj, self.user_global_ns, self.user_ns)
  File "<ipython-input-6-7b1bb90b840a>", line 1, in <module>
    ds = hypercoast.read_aviris("test9.tif")
         ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "C:\Users\Home\PycharmProjects\hsi\.venv\Lib\site-packages\hypercoast\aviris.py", line 52, in read_aviris
    wavelength = ds["wavelength"].values.tolist()
                 ~~^^^^^^^^^^^^^^
  File "C:\Users\Home\PycharmProjects\hsi\.venv\Lib\site-packages\xarray\core\dataset.py", line 1578, in __getitem__
    raise KeyError(
KeyError: "No variable named 'wavelength'. Variables on the dataset include ['band', 'x', 'y', 'spatial_ref', 'band_data']"

How can I save an AVIRIS file to be compatible with the read_aviris function?

[giswqs]

You can save the DataArray as an image and use add_raster() method to load the image.

hypercoast.aviris_to_image(dataset, wavelengths=[400, 500, 600], output="aviris.tif")
m = hypercoast.Map()
m.add_raster("aviris.tif")
m