从栅格文件中提取数据

问题描述

我有3个栅格文件，值的范围是0到2。我试图提取所有3个文件中值都大于1的栅格。我想要的输出是一个栅格，其中包含值大于1的栅格和值。

如何在Python或ArcGIS中获得所需的输出？

解决方法

使用 R 您可以

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

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