rio-cloudmask icon indicating copy to clipboard operation
rio-cloudmask copied to clipboard

Handle absence of TIRS Band 10

Open perrygeo opened this issue 9 years ago • 0 comments

When the TIRS band 10 is all-zero, as is the case with e.g. LC80460282016097LGN00, the equations 7&8 are filled with nans

    /Users/mperry/env/mapbox35/lib/python3.5/site-packages/numpy/lib/nanfunctions.py:1001: RuntimeWarning: All-NaN slice encountered
      warnings.warn("All-NaN slice encountered", RuntimeWarning)

Which leads to the masking of the entire scene, as if there were 100% clouds.

The fmask algorithm is capable of dealing with the absence of a thermal band. We just need to

  • [ ] figure out why some seemingly arbitrary scenes have blank thermal bands (is this a known TIRS issue @celoyd?)
  • [ ] adjust the implementation to fallback to a non-thermal mode.

perrygeo avatar Oct 04 '16 14:10 perrygeo