Crear condicional para verificar si el ráster tiene valor

Question

Tengo un problema con este código:

 ac = ExtractByMask(raster, "area.shp")
arcpy.RasterToPoint_conversion(ac, "point.shp", "VALUE")
 

El problema es que a veces el ráster no tiene valor. (Es una acumulación de flujo y la precisión de mi ráster no es suficiente para lo que estoy preguntando, pero ese no es el punto). Entonces, RasterToPoint no puede funcionar.

Entonces mi pregunta es:
¿Cómo puedo hacer una condición como: "if my raster has no value, don't do the RasterToPoint"

Answer 1

Hay varias formas de incorporar propiedades de trama en una instrucción condicional. Aquí hay tres métodos:

 import arcpy, os, numpy

ras = r"C:\path\to\raster.tif"

# Method 1 Using Spatial Analyst
if arcpy.sa.Raster(ras).maximum > 0:
    arcpy.RasterToPoint_conversion(ras, r"C:\temp\out.shp", "VALUE")

# Method 2 using the raster properties
if arcpy.GetRasterProperties_management(ras, "MAXIMUM").getOutput(0) > 0:
    arcpy.RasterToPoint_conversion(ras, r"C:\temp\out.shp", "VALUE")

# Method 3 Describing the raster via a numpy array
array = arcpy.RasterToNumPyArray(ras)

if numpy.max(array) > 0:
    arcpy.RasterToPoint_conversion(ras, r"C:\temp\out.shp", "VALUE")