#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> _______________________________________________ grass-dev mailing list grass-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-dev