I have 99 LS8 scenes (mountainous terrain) and before processing them any further I would like to mask out cloudy and snow pixel since they might mess up the topographic correction I have to apply. My ultimate goal is a NDVI/EVI time series, so I can treat each scene for itself. I tried the USGS surface product since it contains already a cloud mask but it contains too many artefacts for my study area and also maps out landslides, which are my subject of study. I also tried FMask, providing really bad results.
Does someone have any recommendation how to handle this problem? I mean, in the end it might be quicker to build a separate mask for every scene instead of struggeling but an automatic approach would be ideal (like a DT).