rio-cloudmask
rio-cloudmask copied to clipboard
Rasterio plugin for identifying clouds in multi-spectral satellite imagery
Hello, I found a small typo in the algorithm. According to the paper, it is band 5 that is used for brightness test, i.e. swir1 and not nir as coded...
@perrygeo I prefer your clean framework here a million times over the prohibiting interfaces of the other implementations of FMask. 😄 But are you planning on adding some of the...
When the TIRS band 10 is all-zero, as is the case with e.g. [LC80460282016097LGN00](http://landsat-pds.s3.amazonaws.com/L8/046/028/LC80460282016097LGN00/index.html), 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",...
There are a few global operations (scene-wide quantiles, others?) that preclude a simplistic block iteration approach. But we should be able to trade a bit of speed, maybe doing two...
We should test each function in `equations.py` using hypothesis to suss out any edge cases.
Proximate failure: the `temp_land` function returns a single `nan` which prevents the `tlow, thigh = temp_land(pcps, water, tirs1)` from unpacking (expects two variables) But the bigger question is, what should...
The naive cloudshadow approach in #3 is better than nothing but barely acceptable. Let's study the full-blown cloud/shadow matching and implement it here.