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

GRASS GIS trac at osgeo.org
Sat May 7 10:21:23 PDT 2016

#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
 Keywords:  i.atcorr  |        CPU:  x86-64
 Platform:  All       |
 I'm trying to add WorldView3 support but at time the needed helper script
 does not behave...:

 The issue (trivial?) is reproducible with landsat8 as well:

 cd imagery/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-
 line 477, in _evaluate
     out_of_bounds = self._check_bounds(x_new)
   File "/usr/lib64/python2.7/site-
 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.

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

More information about the grass-dev mailing list