nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Raster to Polygon fails with empty input

JosephChampagne

When I attempt to run the Raster to Polygon tool with a NoData raster (801 rows and 1400 columns), a corrupted shapefile is output and the following error message appears:

ERROR 010067: Error in executing grid expression.
Failed to execute (RasterToPolygon).

This is an issue because I'm trying to batch process a large number of watershed rasters (some of which are empty) using a Python script (832 tiles), and the script is interrupted when it encounters the first empty raster. Is there any known workaround for this issue? Alternately, is there any way to programmatically check for these invalid inputs?

Finde mehr Beiträge mit dem Tag

Accepted answers

All comments

SteveLynch

look at the GetRasterProperties gp tool

-Steve

NeilAyres

Here is an example of getting the % of nulls in a raster.


# calculate % no data cells in a raster.
import sys, os, arcpy
from arcpy import env
from arcpy.sa import *


if arcpy.CheckExtension("Spatial"):
    arcpy.CheckOutExtension("Spatial")
else:
    print "No SA licence"
    exit


HomeDir = "c:/Data/Swala/BurkinaFaso/IPSurvey/"
fgdb = "IPRasterData.gdb"
InData = os.path.join(HomeDir, fgdb)


env.workspace = InData
env.overwriteOutput = True


ListRast = arcpy.ListRasters()
# ListTable = arcpy.ListTables()


for r in ListRast:
    print "Processing", r
    Ras = arcpy.Raster(r)
    NullRas = IsNull(Ras)
    NCol = NullRas.width
    NRow = NullRas.height
    NCell = NCol * NRow
    Mean = NullRas.mean
    NNull = Mean * NCell
    print "{} % Null {}".format(r, NNull / NCell * 100)

BenLeslie1

Hi - I'm just trying to work round this error in my script....

"To test if a Raster has all NoData run the IsNull tool. If the output is all 1 then it only has NoData"

How do you test whether the whole raster is 1?

SteveLynch

In 10SP1 we fixed it so that Raster To Polygon with all NoData will throw
ERROR 010151: No features found

To test if a Raster has all NoData run the IsNull tool. If the output is all 1 then it only has NoData

Steve

automate_fire.png

DanielBennett

Hi,

Same problem as above. What is the fastest (python script) test that tells me if a raster is all NoData?

I'm thinking SearchCursor:

rows = arcpy.SearchCursor("C:/Data/raster", " ??WHERE NODATA??", "", "", "")
if is_empty(rows):
empty = True
else:
empty = False

Nimm an der Unterhaltung teil

Angemeldete Mitglieder können Beiträge verfassen, Updates folgen und mehr. Neu hier? Registriere ein kostenloses Konto.

Schnellzugriffe

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.