AxisError

[skimprem]

If you run the demcompare with example data, you will get an error:

demcompare sample_config.json --loglevel DEBUG

2024-10-15/15:32:30 :: ERROR ::  Demcompare Traceback (most recent call last):
  File "/home/roman/gitrepo/dems/.venv/lib/python3.12/site-packages/demcompare/demcompare.py", line 85, in main
    demcompare.run(args.config, loglevel=args.loglevel)
  File "/home/roman/gitrepo/dems/.venv/lib/python3.12/site-packages/demcompare/__init__.py", line 177, in run
    input_stats_ref = compute_dem_slope(input_stats_ref)
                      ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/home/roman/gitrepo/dems/.venv/lib/python3.12/site-packages/demcompare/dem_tools.py", line 1011, in compute_dem_slope
    lonr = np.roll(lon, 1, 1)
           ^^^^^^^^^^^^^^^^^^
  File "/home/roman/gitrepo/dems/.venv/lib/python3.12/site-packages/numpy/_core/numeric.py", line 1268, in roll
    axis = normalize_axis_tuple(axis, a.ndim, allow_duplicate=True)
           ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/home/roman/gitrepo/dems/.venv/lib/python3.12/site-packages/numpy/_core/numeric.py", line 1439, in normalize_axis_tuple
    axis = tuple([normalize_axis_index(ax, ndim, argname) for ax in axis])
                  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
numpy.exceptions.AxisError: axis 1 is out of bounds for array of dimension 1

[duboise-cnes]

Hello @skimprem, Thank you for the issue. We haven't tested demcompare in python 3.12 and it appears that rasterio is upgraded and some rasterio functions have not the same output, which breaks some internal slope computation function in demcompare. The quickest fix is to downgrade rasterio to 1.3.11:

pip uninstall rasterio pip install rasterio==1.3.11 We shall produce a fix through an upcoming commit, but no release is foreseen yet.

Thanks again for the bug anyway, hope the quick fix answers your need.