Combine Tool Fails in Arcpy but not in Model Builder

Discussion created by jxyu on Jul 2, 2013
Latest reply on Jul 2, 2013 by rzufelt

I have a model builder tool that uses the Combine tool in one of the steps, and have increased the maximum number of unique values to render in the ArcCatalog and ArcMap's Raster Dataset option to accommodate the number of unique value. The tool works. However, when I trying to call the combine tool with the same inputs in arcpy, it is not recognizing the setting, and would give the following error message: "ExecuteError: ERROR 999999: Error executing function. ("esriDataSourcesRaster.RasterCalcUniqueValues") Too many unique values".

I am using 10.1, are there modules and environmental setting that I should be calling to resolve the issue?

Thank you all!

The codes runs right up to the combine step.

import arcpy, os, sys
from arcpy import env
from arcpy.sa import *

env.overwriteOutput = "True"

##RegionFile= arcpy.GetParameterAsText(0)
RegionFile = r'C:\Temp\RegionFile.csv'

#constant input
inBFE = r'C:\Work\SDE_Connections_us\SV2_FEMA.sde\FEMA.FEMA.Polyline_01_16_13\FEMA.FEMA.BFE_01_16_13'
NED30M = r'\\rsarcsv3\Model_08\ArcSDE\SDE_Connections\SV2_NED_30m.sde\NED_30m.NED_30M.NED_30M'

rows = arcpy.SearchCursor(RegionFile, "", "", "Region", "")
currentRegion = ""
for row in rows:
        currentRegion = row.Region
        print currentRegion
        MF = r'C:\Work\Flowlines'
        fGDB = arcpy.CreateFileGDB_management(r'C:\Temp', "R" + currentRegion, "CURRENT")
        env.workspace = fGDB
##        fGDB = r'C:\Temp\R.gdb'
        outClipBFE = str(fGDB) + "\clipBFE" 
        outIntersetBFE_MF = str(fGDB) + "\BFE_MF_Intersect"
        outpt2Ras = str(fGDB) + "\BFE_MF_Pt2Ras"
        outint100yr = str(fGDB) + "\int100yr"
        outintNED30 = str(fGDB) + "\intNED20"
        outCombine = str(fGDB) + "\combineBFE_NED"

        twentyoneLyrFolder = r'C:\Work\21Layers'
        airFldMap = twentyoneLyrFolder + "\Region" + currentRegion + "\wl_R" + currentRegion + "_tif_09.tif"
        regionBND = "C:\Work\BOUNDARY" + currentRegion + "_BuffClip"

        env.cellSize = FldMap
        env.extent = FldMap
        env.snapRaster = NED30M
        clipBFE = arcpy.Clip_analysis(inBFE, regionBND, outClipBFE)
        intersectBFE_MF = arcpy.Intersect_analysis([clipBFE, MF], outIntersetBFE_MF, "ALL", "", "POINT")
        pt2Ras = arcpy.PointToRaster_conversion(outIntersetBFE_MF, "OBJECTID", outpt2Ras, "MOST_FREQUENT", "NONE", "30")

        #multiplying raster by 1000
        int100yr = (SetNull(airFldMap <= -9999, FldMap))* 1000
        intNED30 = Raster(NED30M) * 1000
        #combine Rasters
        combineBFE_NED= Combine([outpt2Ras, outintNED30, outint100yr])