I read about your Landsat 8 cloud segmentation approach here:
… and here:
I like the simplicity of your approach and implemented it in the google earth engine:
… but ran in to obvious problems:
(1) the floor command must be a mistake as it reduces the term to zero (term3_1)
(2) the output does not result in a cloud mask (term 1 results in high numbers > 1.5 , which is making the comparison with the SWIR2 band pointless)
I double-checked the formula (as provided in the resources). It is implemented correctly as it is stated.
I am wondering if there is anything I overlooked.
Any help to resolve this issue would be very much appreciated.