[GRASS-dev] Strange Python error in i.atcorr's create_iwave.py script

Markus Neteler neteler at osgeo.org
Fri May 6 08:54:33 PDT 2016


I'm trying to add WorldView3 support but at time the needed helper
script does not behave...

Reproducible with landsat as well:

cd imageey/i.atcorr/

python create_iwave.py sensors_csv/landsat_8.csv
Getting sensor name from csv file: landsat_8
Traceback (most recent call last):
  File "create_iwave.py", line 250, in <module>
  File "create_iwave.py", line 237, in main
    write_cpp(bands, values, sensor, os.path.dirname(inputfile))
  File "create_iwave.py", line 156, in write_cpp
    fi, li = interpolate_band(values[:,[0,b+1]])
  File "create_iwave.py", line 94, in interpolate_band
    filter_f = f(np.arange(response[0,0], response[-1,0] + 2.5, 2.5))
  File "/usr/lib64/python2.7/site-packages/scipy/interpolate/polyint.py",
line 79, in __call__
    y = self._evaluate(x)
  File "/usr/lib64/python2.7/site-packages/scipy/interpolate/interpolate.py",
line 477, in _evaluate
    out_of_bounds = self._check_bounds(x_new)
  File "/usr/lib64/python2.7/site-packages/scipy/interpolate/interpolate.py",
line 507, in _check_bounds
    raise ValueError("A value in x_new is above the interpolation "
ValueError: A value in x_new is above the interpolation range.

The script is unchanged. Any ideas how to fix that?


More information about the grass-dev mailing list