topic Re: ERROR 999998 - tabulate area tool in Python script in Python Questions

ERROR 999998 - tabulate area tool in Python script

MollyMoore — Sat, 11 Dec 2021 13:04:54 GMT

I have a Python script that I would like to operate as a standalone script. When I copy and paste it into the Python window in a fresh ArcMap .mxd, it operates fine. However, when I attempt to run it as a standalone script within my Python IDE, it fails during the tabulate area tool that is written like this:

tab_area = TabulateArea(target_features,"unique_id",sdm,"Value",os.path.join(outTable,"tab_area_temp"),"30")

I am getting the ambiguous 999998 error message:

ExecuteError: ERROR 999998: Unexpected Error.

Does anyone have any ideas as to why this might be working fine when copied and pasted into the Python window in ArcMap, but failing as a standalone script?

Full Python script is below for reference:

import arcpy, os, datetime
from arcpy import env
from arcpy.sa import *
from itertools import groupby
from operator import itemgetter

print "start time: " + datetime.datetime.now().strftime("%Y-%m-%d %H:%M:%S")

arcpy.env.overwriteOutput = True

target_features = r'C:\Users\mmoore\Documents\ArcGIS\COA.gdb\PlanningUnit_Hex10acre'
sdm_folder = r'C:\Users\mmoore\Documents\ArcGIS\SDMs.gdb'
outTable = r'C:\Users\mmoore\Documents\ArcGIS\out_tables.gdb'
arcpy.env.workspace = sdm_folder
sdms = arcpy.ListRasters()

species_tables = []

for sdm in sdms:
    print sdm + " started at: " + datetime.datetime.now().strftime("%Y-%m-%d %H:%M:%S")
    sdm_path = os.path.join(sdm_folder,sdm)
    arcpy.CheckOutExtension("Spatial")
    tab_area = TabulateArea(target_features,"unique_id",sdm,"Value",os.path.join(outTable,"tab_area_temp"),"30")
    sdm_table = arcpy.TransposeFields_management(tab_area,"VALUE_1 Low;VALUE_2 Medium",os.path.join(outTable,sdm),"OccProb","PERCENTAGE","unique_id")
    arcpy.AddField_management(sdm_table,"ElSeason","TEXT","","",20,"ElSeason")
    with arcpy.da.UpdateCursor(sdm_table,["OccProb","PERCENTAGE","ElSeason"]) as cursor:
        for row in cursor:
            row[1] = (int(row[1])/40468.383184)*100
            cursor.updateRow(row)
            row[2] = sdm
            cursor.updateRow(row)
            if row[0] == "VALUE_1":
                row[0] = "Low"
            elif row[0] == "VALUE_2":
                row[0] = "Medium"
                cursor.updateRow(row)

    with arcpy.da.UpdateCursor(sdm_table, "PERCENTAGE") as cursor:
        for row in cursor:
            if row[0] > 10:
                pass
            else:
                cursor.deleteRow()

    case_fields = ["unique_id", "ElSeason"]
    max_field = "PERCENTAGE"
    sql_orderby = "ORDER BY {}, {} DESC".format(",".join(case_fields), max_field)

    # use groupby cursor to order by percent coverage
    with arcpy.da.UpdateCursor(sdm_table, "*", sql_clause=(None, sql_orderby)) as cursor:
        case_func = itemgetter(*(cursor.fields.index(fld) for fld in case_fields))
        for key, group in groupby(cursor, case_func):
            next(group)
            for extra in group:
                cursor.deleteRow()

    species_tables.append(sdm_table)
    arcpy.Delete_management(tab_area)
    print sdm + " finished at: " + datetime.datetime.now().strftime("%Y-%m-%d %H:%M:%S")

merge = arcpy.Merge_management(species_tables, os.path.join(outTable, "SGCN_SDMxPU"))
print "end time: " + datetime.datetime.now().strftime("%Y-%m-%d %H:%M:%S")‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

The full error message is as follows:

Traceback (most recent call last):
  File "H:\Scripts\COA\COA\COA_species_populator v.7_SDMraster.py", line 23, in <module>
    tab_area = TabulateArea(target_features,"unique_id",sdm,"Value",os.path.join(outTable,"tab_area_temp"),30)
  File "C:\Program Files (x86)\ArcGIS\Desktop10.5\ArcPy\arcpy\sa\Functions.py", line 6183, in TabulateArea
    processing_cell_size)
  File "C:\Program Files (x86)\ArcGIS\Desktop10.5\ArcPy\arcpy\sa\Utils.py", line 53, in swapper
    result = wrapper(*args, **kwargs)
  File "C:\Program Files (x86)\ArcGIS\Desktop10.5\ArcPy\arcpy\sa\Functions.py", line 6175, in Wrapper
    processing_cell_size)
  File "C:\Program Files (x86)\ArcGIS\Desktop10.5\ArcPy\arcpy\geoprocessing\_base.py", line 510, in <lambda>
    return lambda *args: val(*gp_fixargs(args, True))
ExecuteError: ERROR 999998: Unexpected Error.
Failed to execute (TabulateArea).‍‍‍‍‍‍‍‍‍‍‍‍‍

Re: ERROR 999998 - tabulate area tool in Python script

MitchHolley1 — Wed, 13 Dec 2017 20:53:17 GMT

Can you post the entire code here using the 'Syntax Highlighter' in the "More" menu. We need to know what the variables you're using refer to.

Re: ERROR 999998 - tabulate area tool in Python script

DanPatterson_Retired — Wed, 13 Dec 2017 20:53:25 GMT

too much is missing, but I can surmise that

either arcpy can't be found...
the paths to the files can't be found
arcpy.overwriteOutput hasn't been set to True or it isn't even in your script if you are not using a new file as output.
many other possible details, but they would be script and/or environment specific

Re: ERROR 999998 - tabulate area tool in Python script

MollyMoore — Wed, 13 Dec 2017 20:57:36 GMT

Just posted the whole script.

Re: ERROR 999998 - tabulate area tool in Python script

DanPatterson_Retired — Wed, 13 Dec 2017 21:08:39 GMT

One thing I did notice is that you checkout SA within the loop, but you try to import it at the top of the script. Within ArcMap, the rules are more accommodating. Perhaps try checking it out first, then import SA and I would throw in some print statements to make sure everything is going out as planned.

Re: ERROR 999998 - tabulate area tool in Python script

MitchHolley1 — Sat, 11 Dec 2021 13:04:56 GMT

So the code is erroring out on line 23? Could you post the entire error.

In the UpdateCursor loop starting at line 26, you do not need to update the row every time you specify a row value.

with arcpy.da.UpdateCursor(sdm_table,["OccProb","PERCENTAGE","ElSeason"]) as cursor:
    for row in cursor:
        row[1] = (int(row[1])/40468.383184)*100
        row[2] = sdm
        if row[0] == "VALUE_1":
            row[0] = "Low"
        elif row[0] == "VALUE_2":
            row[0] = "Medium"
        cursor.updateRow(row)

del cursor‍‍‍‍‍‍‍‍‍‍‍

Also, I would delete cursors after you're done with them.

Re: ERROR 999998 - tabulate area tool in Python script

MollyMoore — Wed, 13 Dec 2017 21:30:31 GMT

Thanks for the suggestions. I will make those updates.

I edited my post to include the full error message that I'm getting.

Re: ERROR 999998 - tabulate area tool in Python script

DanPatterson_Retired — Wed, 13 Dec 2017 21:45:00 GMT

If that is a layer or raster layer... a standalone script will fail unless it knows the path to the data ... which you appear to have.

that leaves the fact that you are trying to join a raster to a path... but I would still throw a print statement in

Re: ERROR 999998 - tabulate area tool in Python script

MitchHolley1 — Sat, 11 Dec 2021 13:04:59 GMT

Try specifying the processing cell size as an integer and not as a string.

Also, if you're looping over multiple rasters, you're going to need a specific output for each one. Right now, you're looping over the rasters and saving the TabulateArea output table as a single name. I would try something like this:

rastTable = r"{0}_{1}".format(os.path.join(outTable, "tab_area_temp"),
                              str(sdm))

tab_area = TabulateArea(target_features,
                        "unique_id",
                        sdm,
                        "Value",
                        rastTable,
                        30)‍‍‍‍‍‍‍‍‍‍‍‍‍

Re: ERROR 999998 - tabulate area tool in Python script

DanPatterson_Retired — Wed, 13 Dec 2017 22:12:22 GMT

good catch! But I though Molly said the script worked in ArcMap... strange.. if that is the case... string in one, but number in another? Hope it is that simple

Re: ERROR 999998 - tabulate area tool in Python script

MollyMoore — Thu, 14 Dec 2017 14:15:44 GMT

I tried the cell size as an integer and it is giving the same error.

I inserted some print statements and all the paths seem to be correct. I'm a bit baffled by this one.

Re: ERROR 999998 - tabulate area tool in Python script

DanPatterson_Retired — Thu, 14 Dec 2017 14:30:09 GMT

The only thing left to try is to move the origin files to a new folder location and specify a new destination location. Repeating, the same thing is getting you no where. If that doesn't work then there is something else fundamentally wrong.

Re: ERROR 999998 - tabulate area tool in Python script

MitchHolley1 — Sat, 11 Dec 2021 13:05:02 GMT

Could you try creating a layers out of the rasters before processing them. The below code creates layers for the target and input rasters.

target_features = r'C:\Users\mmoore\Documents\ArcGIS\COA.gdb\PlanningUnit_Hex10acre'

#target raster layer
arcpy.MakeRasterLayer_management(target_features, "targetLyr")

for sdm in sdms:
    print sdm + " started at: " + datetime.datetime.now().strftime("%Y-%m-%d %H:%M:%S")
    sdm_path = os.path.join(sdm_folder,sdm)
    arcpy.MakeRasterLayer_management(sdm_path,"{0}Lyr".format(sdm))#makes in-mem raster
    rastTable = r"{0}_{1}".format(os.path.join(outTable, "tab_area_temp"),
                                  str(sdm)
                                  )
    tab_area = TabulateArea("targetLyr",
                            "unique_id",
                            "{0}Lyr".format(sdm),
                            "Value",
                            rastTable,
                            30
                            )‍‍‍‍‍
    arcpy.Delete_management("{0}Lyr".format(sdm))‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Re: ERROR 999998 - tabulate area tool in Python script

MollyMoore — Thu, 14 Dec 2017 18:16:53 GMT

I have tried converting both inputs to layers and tried running the script with origin files in new locations. Both, to no avail - I am still getting the same error message.

Re: ERROR 999998 - tabulate area tool in Python script

MatthewDriscoll — Thu, 14 Dec 2017 19:13:07 GMT

I suspect the License is not getting checked out. Try catching an exceptions on that before running the Tabulate. See the Code Sample here: CheckOutExtension—Help | ArcGIS for Desktop

Re: ERROR 999998 - tabulate area tool in Python script

MollyMoore — Wed, 20 Dec 2017 14:50:37 GMT

I finally got this to work. But, I'm still wondering why it wouldn't work originally. Here are the different things I tried:

I was originally using a large (2.9 million features) feature class to delineate 'zones' within which to tabulate area of the raster. I tried clipping this to a smaller area and the stand-alone script then ran fine. It still would not work on the larger dataset.

Finally, I converted my zone feature class to a raster and the stand-alone script then ran on the whole dataset.