ArcGIS 10 python - 条件语句检查以查看栅格值是否在集合中

Question

我试图根据另一个栅格的输入创建一个栅格。 如果栅格值等于set（）中包含的数字，那么我希望它为1，否则为0

我尝试过以下方法：

 ConfusedRaster = arcpy.Con(inraster in repeatSet, 1, 0)

和

 ConfusedRaster = arcpy.Con(inraster, 1, 0, "inraster in repeatSet")

这些都不奏效。 我相信它们不起作用，因为where子句只接受Map Algebra表达式： ArcGIS Help

我还可以考虑另外两种方法来做这件事。 一个是将它转换为NumPyArray并使用它。 另一个是循环遍历集合并为集合中的每个值创建栅格对象。 循环完成后合并它们。

有没有人对如何解决这个问题有任何建议或意见？

谢谢

Answer 1

我正在寻找类似问题的答案，并使用'ExtractByAttributes'中的SQL子句开发了一种方法。

repeatList = list(repeatSet)
ras1 = arcpy.sa.ExtractByAttributes(inraster, 'VALUE IN (' + str(repeatList).strip('[]') + ')') 
ConfusedRaster = arcpy.sa.Con(arcpy.sa.IsNull(ras1) == 0, 1, 0)