Landsat-8/OLI error raised on interpolate

Post Reply
dndat
Posts: 4
Joined: Sat May 30, 2020 5:26 am
company / institution: Space Technology Institute, VAST
Location: Hanoi, Vietnam

Landsat-8/OLI error raised on interpolate

Post by dndat »

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.
User avatar
fsteinmetz
Site Admin
Posts: 306
Joined: Fri Sep 07, 2018 1:34 pm
company / institution: Hygeos
Location: Lille, France
Contact:

Re: Landsat-8/OLI error raised on interpolate

Post by fsteinmetz »

Hello,
I could process the same Landsat8 file without an error with v4.13 (which I see you are using as well).
Do you get the same error on another L8 file ? Could you post your run script ?
dndat
Posts: 4
Joined: Sat May 30, 2020 5:26 am
company / institution: Space Technology Institute, VAST
Location: Hanoi, Vietnam

Re: Landsat-8/OLI error raised on interpolate

Post by dndat »

Hi, This is the script

from polymer.main import run_atm_corr, Level1, Level2
from polymer.level1_landsat8 import Level1_OLI
from polymer.level2_nc import Level2_NETCDF

run_atm_corr(
Level1_OLI('/media/dndat/ATC_Proc/LC08_L1TP_045011_20130813_20170309_01_T1/'),
Level2_NETCDF(outdir='/media/dndat/ATC_Proc/', ext='LC08_L1TP_045011_20130813_20170309_01_T1_polymer.nc', overwrite=True, tmpdir='/tmp/'), multiprocessing=-1
)

I did not get this error in tropical area. It happened for the scenes close to polar area
Cheers
dndat
Posts: 4
Joined: Sat May 30, 2020 5:26 am
company / institution: Space Technology Institute, VAST
Location: Hanoi, Vietnam

Re: Landsat-8/OLI error raised on interpolate

Post by dndat »

Hello
is it error come from l8_angles? I guess
cheers
User avatar
fsteinmetz
Site Admin
Posts: 306
Joined: Fri Sep 07, 2018 1:34 pm
company / institution: Hygeos
Location: Lille, France
Contact:

Re: Landsat-8/OLI error raised on interpolate

Post by fsteinmetz »

Possibly. I can not reproduce the error.
Did you run l8_angles with :
l8_angles LC08_L1TP_045011_20130813_20170309_01_T1_ANG.txt BOTH 1 -b 1 ?
dndat
Posts: 4
Joined: Sat May 30, 2020 5:26 am
company / institution: Space Technology Institute, VAST
Location: Hanoi, Vietnam

Re: Landsat-8/OLI error raised on interpolate

Post by dndat »

Yes, I did it
at other location, polymer did well

Reading coordinates from /media/dndat/ATC_Proc/LC08_L1TP_015043_20131201_20170307_01_T1/LC08_L1TP_015043_20131201_20170307_01_T1_B1.TIF
Image size is 7611x7761
Trying to download https://oceandata.sci.gsfc.nasa.gov/cgi ... MI_24h.hdf ...
2020-06-10 15:39:45 URL:https://oceandata.sci.gsfc.nasa.gov/ob/ ... MI_24h.hdf [201413/201413] -> "ANCILLARY/METEO/2013/335/N201333500_O3_AURAOMI_24h.hdf.tmp" [1]
Trying to download https://oceandata.sci.gsfc.nasa.gov/cgi ... 6h.hdf.bz2 ...
2020-06-10 15:39:48 URL:https://oceandata.sci.gsfc.nasa.gov/ob/ ... 6h.hdf.bz2 [830485/830485] -> "ANCILLARY/METEO/2013/335/N201333512_MET_NCEPR2_6h.hdf.bz2.tmp" [1]
Trying to download https://oceandata.sci.gsfc.nasa.gov/cgi ... 6h.hdf.bz2 ...
2020-06-10 15:39:52 URL:https://oceandata.sci.gsfc.nasa.gov/ob/ ... 6h.hdf.bz2 [829118/829118] -> "ANCILLARY/METEO/2013/335/N201333518_MET_NCEPR2_6h.hdf.bz2.tmp" [1]
Starting processing at 2020-06-10 15:39:58.583116
Initializing output file "/media/dndat/ATC_Proc/LC08_L1TP_015043_20131201_20170307_01_T1_polymer.nc"
Processing block: size (500, 400), offset (0, 0)
Processing block: size (500, 400), offset (0, 400)
Processing block: size (500, 400), offset (0, 800)
Processing block: size (500, 400), offset (0, 1200)
Processing block: size (500, 400), offset (0, 1600)
Processing block: size (500, 400), offset (0, 2000)
Processing block: size (500, 400), offset (0, 2400)
Processing block: size (500, 400), offset (0, 2800)
Processing block: size (500, 400), offset (0, 3200)
Processing block: size (500, 400), offset (0, 3600)
Processing block: size (500, 400), offset (0, 4000)
Processing block: size (500, 400), offset (0, 4400)
Processing block: size (500, 400), offset (0, 4800)
Processing block: size (500, 400), offset (0, 5200)
Processing block: size (500, 400), offset (0, 5600)
Processing block: size (500, 400), offset (0, 6000)
Processing block: size (500, 400), offset (0, 6400)
Processing block: size (500, 400), offset (0, 6800)
Processing block: size (500, 400), offset (0, 7200)
Processing block: size (500, 11), offset (0, 7600)
...
Post Reply