条件函数 (空间分析)

需要 Spatial Analyst 许可。

获得 Image Analyst 许可后可用。

摘要

可对输入栅格的每个输入像元执行 if/else 条件评估。

了解有关使用条件函数工具执行条件评估的详细信息

插图

条件函数图示
OutRas = Con(InRas1, 40, 30, "Value >= 2")

使用情况

  • 如果真栅格数据或可选假栅格数据为浮点型,则输出栅格数据也将为浮点型。 如果真表达式和可选假栅格数据均为整型,则输出栅格数据也将为整型。

  • 如果输入条件栅格(Python 中的 in_conditional_raster)是单波段栅格,并且输入条件为真时所取的栅格数据或常量值(Python 中的 in_true_raster_or_constant)栅格或可选的输入条件为假时所取的栅格数据或常量值(Python 中的 in_false_raster_or_constant)栅格是常量,输出将是单波段栅格。

  • 如果所有输入都是多波段栅格,则输出将是多波段栅格。 如果真输入或可选的假输入为常量,则输出栅格也将是多波段的。 每个多波段输入中的波段数必须相同。

  • 该工具将使用来自其他输入的相应波段在条件栅格的每个波段上执行操作。 如果条件输入是多波段栅格,并且真或假栅格输入是常量,则该工具将使用多波段输入中每个波段的常量值执行操作。

  • 如果表达式的评估结果非零,则将被视为 true。

  • 如果未指定输入条件为假时所取的栅格数据或常量值,则将为表达式结果不为 true 的那些像元分配 NoData。

  • 如果 NoData 不满足表达式,则像元不会接收输入条件为假时所取的栅格数据值;像元值仍是 NoData。

  • 表达式使用 SQL 查询。 有关创建查询的详细信息,请参阅以下主题:

  • 要在 Python 中使用 {where_clause},应以引号括起来。 例如,"Value > 5000”

    有关使用 Python 指定查询的详细信息,请参阅帮助。

  • 在 Python 中,您可避免使用 {where_clause},其通过将地图代数表达式用作 in_conditional_raster 来指定 Value 字段。

    例如,以下表达式:

    • Con("elev", 0, 1, "value > 1000")

    可被重新编写,如下所示:

    • Con(Raster("elev") > 1000, 0, 1)

    有关详细信息,请参阅下列代码示例或在地图代数中查看构建复杂语句

  • 逻辑表达式的最大长度为 4,096 个字符。

  • 如果其中至少一个输入是有相同变量数的多维栅格数据,则该工具将对具有相同维度值的所有剖切片执行运算。 输出将是 CRF 格式的多维栅格。 输入中的变量必须至少具有一个常见维度和一个常见维度值以供此工具进行处理,否则将发生错误。

    如果任意两个输入为多维栅格并共享一个变量(但名称不同),则请在运行工具之前取消选中匹配多维变量地理处理环境(在 Python 中设置 arcpy.env.matchMultidimensionalVariable = False)。

    如果输入条件栅格值是多维栅格,而输入条件为真时所取的栅格数据或常量值输入条件为假时所取的栅格数据或常量值参数设置为常量值,则工具将使用常量值针对所有剖切片的所有变量执行运算,而输出将是多维栅格。

参数

标注说明数据类型
输入条件栅格数据

此输入栅格表示所需条件的 true 或 false 结果。

可为整型或浮点型。

Raster Layer
输入条件为真时所取的栅格数据或常量值

条件为 true 时,值将作为输出像元值的输入。

它可以是整型栅格、浮点型栅格或常量值。

Raster Layer; Constant
输入条件为假时所取的栅格数据或常量值
(可选)

条件为 false 时,值将作为输出像元值的输入。

它可以是整型栅格、浮点型栅格或常量值。

Raster Layer; Constant
表达式
(可选)

此逻辑表达式用于确定为 true 或 false 的输入像元。

Where 子句遵循 SQL 表达式的一般格式。 如果您单击编辑 SQL 模式按钮 SQL 查询,则可以直接输入,例如 VALUE > 100。 如果处于编辑子句模式 编辑子句 中,则可以通过单击添加子句模式按钮来开始构建表达式。

SQL Expression

返回值

标注说明数据类型
输出栅格

输出栅格。

Raster

Con(in_conditional_raster, in_true_raster_or_constant, {in_false_raster_or_constant}, {where_clause})
名称说明数据类型
in_conditional_raster

此输入栅格表示所需条件的 true 或 false 结果。

可为整型或浮点型。

Raster Layer
in_true_raster_or_constant

条件为 true 时,值将作为输出像元值的输入。

它可以是整型栅格、浮点型栅格或常量值。

Raster Layer; Constant
in_false_raster_or_constant
(可选)

条件为 false 时,值将作为输出像元值的输入。

它可以是整型栅格、浮点型栅格或常量值。

Raster Layer; Constant
where_clause
(可选)

此逻辑表达式用于确定为 true 或 false 的输入像元。

表达式遵循 SQL 表达式的一般格式。 where_clause 的示例为 "VALUE > 100"

SQL Expression

返回值

名称说明数据类型
out_raster

输出栅格。

Raster

代码示例

Con 示例 1(Python 窗口)

在本例中,当输入条件栅格值大于 2000 时,原始值将保留在输出中;否则,该值将为 NoData。

import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outCon = Con("elevation", "elevation", "", "VALUE > 2000")
outCon.save("C:/sapyexamples/output/outcon.img")

# Execute Con using a map algebra expression instead of a where clause
outCon2 = Con(Raster("elevation") > 2000, "elevation")
outCon2.save("C:/sapyexamples/output/outcon2")
Con 示例 2(Python 窗口)

在此例中,原始值将被保留在除 NoData 外的输出中,其将被用 0 代替。

import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outCon = Con(IsNull("elevation"),0, "elevation")
outCon.save("C:/sapyexamples/output/outcon")
Con 示例 3(Python 窗口)

在此例中,两个不同的栅格用于创建条件栅格。

import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
inRaster1 = Raster("landuse")
inRaster2 = Raster("landuse2")
outCon = Con(((inRaster1 == 1) & (inRaster2 == 5)), inRaster1 + inRaster2, 99)
outCon.save("C:/sapyexamples/output/outcon")
Con 示例 4(Python 窗口)

在本例中,Con 中将使用多个 Con 工具。

import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
inRas1 = Raster("inRaster")
outCon = Con(inRas1 < 45,1, Con((inRas1 >= 45) & (inRas1 < 47),2, Con((inRas1 >= 47) & (inRas1 < 49),3, Con(inRas1 >= 49,4))))
outCon.save("C:/sapyexamples/output/outcon")
Con 示例 5(独立脚本)

在本例中,如果输入条件栅格数据值大于或等于 1500,则输出值将为 1;如果输入条件栅格数据值小于 1500,则输出值将为 0。

# Name: Con_Ex_02.py
# Description: Performs a conditional if/else evaluation 
#              on each cell of an input raster.
# Requirements: Spatial Analyst Extension

# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *

# Set environment settings
env.workspace = "C:/sapyexamples/data"

# Set local variables
inRaster = Raster("elevation")
inTrueRaster = 1
inFalseConstant = 0
whereClause = "VALUE >= 1500"

# Execute Con
outCon = Con(inRaster, inTrueRaster, inFalseConstant, whereClause)

# Execute Con using a map algebra expression instead of a where clause
outCon2 = Con(inRaster >= 1500, inTrueRaster, inFalseConstant)

# Save the outputs 
outCon.save("C:/sapyexamples/output/outcon")
outCon2.save("C:/sapyexamples/output/outcon2")