Cloud and cloud shadow masking for old LANDSAT and ASTER images in GRASS

Question

I am trying to remove cloud and cloud shadows for Landsat 4-5 ASTER images in GRASS GIS.

I have knowledge that detecting cloud and cloud shadow pixels which is >(255*0.2), but I cannot figure out the appropriate syntax for r.mapcalc?

markusN · Accepted Answer · 2016-10-20 18:34:38Z

1

Try something like this:

r.mapcalc "cloudmask = if(lsatband > 255*0.2, null(), lsatband)"
r.mask cloudmask

However, are you sure that the formula is that simple?

answered Oct 20 '16 at 18:34

markusN

11.3k2244

add a comment |

1 Answer 1