Hi,
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>
main()
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?
thanks
Markus
_______________________________________________
grass-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-dev
