Landsat-8/OLI error raised on interpolate
Posted: Sun Jun 07, 2020 2:38 pm
Could you help me to solve this problem, thank you very much
Reading coordinates from /media/dndat/ATC_Proc/LC08_L1TP_045011_20130813_20170309_01_T1/LC08_L1TP_045011_20130813_20170309_01_T1_B1.TIF
Image size is 8601x8641
Starting processing at 2020-06-07 21:31:59.069748
Initializing output file "/media/dndat/ATC_Proc/LC08_L1TP_045011_20130813_20170309_01_T1_polymer.nc"
Traceback (most recent call last):
File "dndatexample.py", line 15, in <module>
multiprocessing=-1)
File "/home/dndat/Software/polymer-v4.13/polymer/main.py", line 517, in run_atm_corr
for block in block_iter:
File "/home/dndat/anaconda3/lib/python3.6/multiprocessing/pool.py", line 735, in next
raise value
File "/home/dndat/anaconda3/lib/python3.6/multiprocessing/pool.py", line 119, in worker
result = (True, func(*args, **kwds))
File "/home/dndat/anaconda3/lib/python3.6/multiprocessing/pool.py", line 138, in _helper_reraises_exception
raise ex
File "/home/dndat/anaconda3/lib/python3.6/multiprocessing/pool.py", line 290, in _guarded_task_generation
for i, x in enumerate(iterable):
File "/home/dndat/Software/polymer-v4.13/polymer/main.py", line 434, in blockiterator
for block in level1.blocks(params.bands_read()):
File "/home/dndat/Software/polymer-v4.13/polymer/level1.py", line 175, in blocks
yield self.read_block(size, offset, bands_read)
File "/home/dndat/Software/polymer-v4.13/polymer/level1_landsat8.py", line 280, in read_block
block.ozone = self.ozone[block.latitude, block.longitude]
File "/home/dndat/Software/polymer-v4.13/polymer/ancillary.py", line 58, in __getitem__
return self.data[Idx(lat), Idx(lon)]
File "/home/dndat/Software/polymer-v4.13/polymer/luts.py", line 399, in __getitem__
keys = k.index(self.axes)
File "/home/dndat/Software/polymer-v4.13/polymer/luts.py", line 1112, in index
fill_value=fv)(self.value)
File "/home/dndat/anaconda3/lib/python3.6/site-packages/scipy/interpolate/polyint.py", line 79, in __call__
y = self._evaluate(x)
File "/home/dndat/anaconda3/lib/python3.6/site-packages/scipy/interpolate/interpolate.py", line 663, in _evaluate
below_bounds, above_bounds = self._check_bounds(x_new)
File "/home/dndat/anaconda3/lib/python3.6/site-packages/scipy/interpolate/interpolate.py", line 692, in _check_bounds
raise ValueError("A value in x_new is below the interpolation "
ValueError: A value in x_new is below the interpolation range.
Reading coordinates from /media/dndat/ATC_Proc/LC08_L1TP_045011_20130813_20170309_01_T1/LC08_L1TP_045011_20130813_20170309_01_T1_B1.TIF
Image size is 8601x8641
Starting processing at 2020-06-07 21:31:59.069748
Initializing output file "/media/dndat/ATC_Proc/LC08_L1TP_045011_20130813_20170309_01_T1_polymer.nc"
Traceback (most recent call last):
File "dndatexample.py", line 15, in <module>
multiprocessing=-1)
File "/home/dndat/Software/polymer-v4.13/polymer/main.py", line 517, in run_atm_corr
for block in block_iter:
File "/home/dndat/anaconda3/lib/python3.6/multiprocessing/pool.py", line 735, in next
raise value
File "/home/dndat/anaconda3/lib/python3.6/multiprocessing/pool.py", line 119, in worker
result = (True, func(*args, **kwds))
File "/home/dndat/anaconda3/lib/python3.6/multiprocessing/pool.py", line 138, in _helper_reraises_exception
raise ex
File "/home/dndat/anaconda3/lib/python3.6/multiprocessing/pool.py", line 290, in _guarded_task_generation
for i, x in enumerate(iterable):
File "/home/dndat/Software/polymer-v4.13/polymer/main.py", line 434, in blockiterator
for block in level1.blocks(params.bands_read()):
File "/home/dndat/Software/polymer-v4.13/polymer/level1.py", line 175, in blocks
yield self.read_block(size, offset, bands_read)
File "/home/dndat/Software/polymer-v4.13/polymer/level1_landsat8.py", line 280, in read_block
block.ozone = self.ozone[block.latitude, block.longitude]
File "/home/dndat/Software/polymer-v4.13/polymer/ancillary.py", line 58, in __getitem__
return self.data[Idx(lat), Idx(lon)]
File "/home/dndat/Software/polymer-v4.13/polymer/luts.py", line 399, in __getitem__
keys = k.index(self.axes)
File "/home/dndat/Software/polymer-v4.13/polymer/luts.py", line 1112, in index
fill_value=fv)(self.value)
File "/home/dndat/anaconda3/lib/python3.6/site-packages/scipy/interpolate/polyint.py", line 79, in __call__
y = self._evaluate(x)
File "/home/dndat/anaconda3/lib/python3.6/site-packages/scipy/interpolate/interpolate.py", line 663, in _evaluate
below_bounds, above_bounds = self._check_bounds(x_new)
File "/home/dndat/anaconda3/lib/python3.6/site-packages/scipy/interpolate/interpolate.py", line 692, in _check_bounds
raise ValueError("A value in x_new is below the interpolation "
ValueError: A value in x_new is below the interpolation range.