On Thu, Mar 13, 2014 at 5:28 PM, César Augusto Ramírez Franco <[email protected]> wrote: > Hello Markus, I'm using the terminal on Debian Testing/SID > > I'm trying to make a mask for the pixels of a landsat 8 image inside a basin > > mask1 is my basin, and mask2 is the cloudmask of the landsat image > > I calculated the cloudmask with > > r.mapcalc "cloudmask = if( (BQA ==1 || BQA >= 49152), null(), 1)" > > Maybe I'm using the wrong logic inside the if statement? How do I get > r.mapcalc to make an intersection of the two masks other than summing them?
I would use the bitpatterns of the QA file. Some example you can find here: http://courses.neteler.org/processing-landsat8-data-in-grass-gis-7/ --> 8. Applying the Landsat 8 Quality Assessment (QA) Band Best Markus _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
