[GRASS-dev] [GRASS GIS] #3027: i.atcorr's create_iwave.py broken

[GRASS GIS]

#3027: i.atcorr's create_iwave.py broken
----------------------+---------------------------------
  Reporter:  neteler  |      Owner:  grass-dev@…
      Type:  defect   |     Status:  new
  Priority:  normal   |  Milestone:  7.0.5
 Component:  Imagery  |    Version:  svn-releasebranch70
Resolution:           |   Keywords:  i.atcorr
       CPU:  x86-64   |   Platform:  All
----------------------+---------------------------------

Comment (by neteler):

 Replying to [comment:5 nikosa]:
 > Patch updated to handle, now hopefully, all cases -- whether the *input*
 step in the wavelength slots is less, equal or greater than the
 *hardcoded* step `2.5` for the interpolation function.

 Thanks for that. Unfortunately it still fails:

 {{{
 python create_iwave.py sensors_csv/landsat_8.csv
 Getting sensor name from csv file: landsat_8
 Number of bands found: 9
 Traceback (most recent call last):
   File "create_iwave.py", line 267, in <module>
     main()
   File "create_iwave.py", line 254, in main
     write_cpp(bands, values, sensor, os.path.dirname(inputfile))
   File "create_iwave.py", line 173, in write_cpp
     fi, li = interpolate_band(values[:,[0,b+1]])
   File "create_iwave.py", line 108, in interpolate_band
     filter_f = f(np.arange(start, stop, step))
   File "/usr/lib/python2.7/dist-packages/scipy/interpolate/polyint.py",
 line 54, in __call__
     y = self._evaluate(x)
   File "/usr/lib/python2.7/dist-
 packages/scipy/interpolate/interpolate.py", line 448, in _evaluate
     out_of_bounds = self._check_bounds(x_new)
   File "/usr/lib/python2.7/dist-
 packages/scipy/interpolate/interpolate.py", line 478, 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.
 }}}

--
Ticket URL: <https://trac.osgeo.org/grass/ticket/3027#comment:7>
GRASS GIS <https://grass.osgeo.org>