我正在尝试根据来自另一个栅格的输入创建一个栅格。
如果栅格值等于 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 并使用它。另一种方法是遍历集合并为集合中的每个值创建一个光栅对象。循环完成后合并它们。

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

谢谢

最佳答案

我正在寻找类似问题的答案,并开发了一种使用“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)

关于python - ArcGIS 10 python - 检查栅格值是否在集合中的条件语句,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/10608809/

10-12 20:12