Extract data from Raster files

Question

I have 3 raster files, the value ranges from 0 to 2. I am trying to extract grids where the value is greater than 1 in all 3 files. My desired output is a raster which contains the grids and values where value is more than 1.

How can I get the desired output in Python or ArcGIS?

Answer 1

With R you can do

library(raster)
s <- stack(system.file("external/rlogo.grd", package="raster")) / 150

x <- all(s > 1)
plot(x)