Latest Contributions by curtvprice

newList is a list of field names, there should be 72 items in it Apparently not, if you're getting that error. I suggest using a try-except block to find which line is causing the problem (see my code for an example). ... View more

What is your ArcGIS version and SP? (See Help > About ArcMap) If not on the most recent, do update! ... View more

I keep getting a �??list index out of range�?� error message when doing this. How many fields are there in the line? If you only have 10, lineList[36] will raise an out of range error. Here's another shot at this: # skip 5 rows skipRows = 5 for n in range(skipRows + 1):     line = inFile.readline() row = skipRows + 1 # read and clean up data while line:     print 'Row', row, 'line info=', line[:72]     # clean line and split into list     lineList = [f.strip().replace("/","_") for f in line.split(',')]     print 'line info =', lineList     latidx,lonidx = 36, 37 # lat and long are in field 36 and 37     if len(lineList) > lonidx:         raise Exception, "line {0} too short".format(row)     for k in latidx, lonidx:         if lineList == "":             lineList = 0     outFile.write(','.join(newList) + "\n")     line = inFile.readline()     row += 1   ... View more

If you construct a an intermediate filename using a built-in variables, the system won't modify them when the model is validated. (Note, %scratchfolder% is 10.1 only, %scratchworkspace% would be used in earlier versions). For example: %scratchfolder%\clip_%n%.tif will be interpreted at runtime as D:\work\scratch\clip_0.tif D:\work\scratch\clip_1.tif ... Depending on the iteration, the built-in iteration variable used could be either %n% or %i%. Hope this helps! ... View more

Neda, Please post your Python code in a code block. [thread=48475]Here's how.[/thread] 1. You may want to look into the Split Raster tool. It may be a lot easier to tile your data using that instead of making fishnets and using this scripting approach. 2. I did some once off code (not tested) to help you with working out how a script would work: # only import arcpy - don't use arcgisscripting/gp at 10.x import arcpy from arcpy import env arcpy.env.overwriteOutput = True ## 1 # in_fc is your polygon feature class in_fc = r"X:\mohammadi\finalmunichGIS\clip0ffishnet9000\shapefilesoffc\fc_parent1" # in_tif is your input in_raster = "X:\mohammadi\finalmunichGIS\Testdaten-original data\munich_subset.tif" # output location outGDB = r"X:\mohammadi\finalmunichGIS\munich.gdb" arcpy.MakeFeatureLayer_management(in_fc,lyrFC) # get list of object ids objIDs = [] objIDField = arcpy.Describe(lyrFC).OIDFieldName Rows = arcpy.SearchCursor(lyrFC) while Row in Rows:     objIDs.append(Row.getValue(objIDField)) del Row, Rows # temp polygons tmpFC = "in_memory/tmpPoly" # one by one clip raster to polygons arcpy.env.workspace= outGDB for objID in objIDs:     where = "{0} = {1}".format(objIDField,objID)     arcpy.SelectLayerByAttributes_management(lyrFC,where)     arcpy.CopyFeatures_management(lyrFC,tmpFC)     ext = arcpy.Describe(tmpFC).Extent     extString = "{0} {1} {2} {3}".format(ext.XMin, ext.YMin, ext.XMax, ext.YMax)     out_raster = "clip{0}".format(objID)     arcpy.Clip_management(in_raster, extString, out_raster, tmpFC)     print arcpy.GetMessages(0) arcpy.Delete_management(tmpFC) ... View more

I am trying to calculation the correlation coefficients for 26 environmental variables. When I put in some of my rasters, I get the following error message.  ERROR 010240: Could not save raster dataset to E:\..... with output format GRID Stack 7.x. Failed to execute (BandCollectionStats). 1. It could be a path name issue The output grid stack name must be short (<10 characters), must start with a letter, and must not contain any non alphanumeric characters except "_". This is a limitation of the Esri grid format. (The extra-short 10 characters is because a grid stack needs to create output grid names <gridname>_c1, etc.). It's possible that some older and little-used tools may still need to write to the grid format as an intermediate step even if the destination is a tiff, so these naming limits may apply even for tiff output in this case. There are bugs sometimes if the path itself contains spaces or special characters, but these issues are becoming more rare with each version of ArcGIS. Encapsulating paths within the raster object at 10.x helped a lot with this. 2. There may be a file size problem with the raster table COUNT for a value. If an integer raster exceeds the 2.1G limit of COUNT for a values in an integer raster table, the raster attribute table build will fail. This usually brings up a different error message, but it's another possibility. This is a hard limit and you may need to tile your data to get it to work. The way to check this is to process a smaller area of the raster (by setting a limited geoprocessing extent) and seeing if it works with a smaller output. This is always a good idea when working out a complex processing chain so you don't have to wait a long time for tools to run that are just going to crash because of some issue you must resolve. When it's all working, you can just reset the extent and run it "for real". Another useful tip - especially when processing large rasters, set the output and scratch workspace to the same location (preferably a folder so the tired-and-true grid format is used). When a tool completes, the output raster can then be renamed instead of copied from scratch to current workspace paths. ... View more

Here's my response to Jordan copied from my email earlier today, with some edits for clarity: Jordan, I totally agree that it's tricky to integrate the two environments. Haven't tried Canopy yet, but I understand Canopy is pretty much the same with a juiced up UI and improved setup tools. I really do believe that the best way to handle this is to install python from the ArcGIS distribution, and separately install Canopy. (You need EPD 32-bit to import 32-bit arcpy and EPD 64-bit to import 64-bit arcpy.) If you install EPD 32 and 64 and ArcGIS Desktop x64 background geoprocessing, this adds up to no less than four python distributions to deal with. Jason Pardy from Esri suggested I just try putting the site-packages folders in the python path. From ArcGIS I have been successful accessing EPD packages by adding a similar .pth file to the ArcGIS python install site-packages folder. I haven't run into any problems with "circular searches", I think Python is smart enough not to look in a site-packages folder twice. One way to modify the Python path is to copy the file Desktop.pth from \Python27\ArcGIS10.1\lib\site-packages to the file zzArcGIS.pth in the 32-bit Canopy site-packages folder (the paths are loaded in alpha order). This should allow you to import arcpy just fine on the Canopy side. Another approach is to just add the paths to the PYTHONPATH variable (see bat scripts below). bat scripts are especially useful if you find you need full control of the PATH and other environment variables to get things to work right. I've attached an example. The way it's done in UNIX -- much less ugly than registering DLLs and messing with globally-applied environment variables! (When you run the scripts, the environment tweaks only apply to the shell and its children, not your global environment.)  A few examples are attached: :: epd32arc.bat :: start EPD32 python prompt with arcpy @echo off set EPDPATH=E:\python27_epd32 set AGSPATH=D:\ArcGIS\Desktop10.1 :: save paths set PATHENV=%PATH% set PPATHENV=%PYTHONPATH% :: set paths EPD with access to arcpy set PATH=%EPDPATH%;%EPDPATH%\scripts;%PATH% set PYTHONPATH=%PYTHONPATH%;%AGSPATH%\bin;%AGSPATH%\arcpy;%AGSPATH%\ArcToolbox\Scripts cmd /c %EPDPATH%\python.exe :: restore paths set PATH=%PATHENV% set PYTHONPATH=%PPATHENV% :: delete variables set EPDPATH= set AGSPATH= set PATHENV= set PPATHENV= :: arcmap_epd.bat - start ArcMap with EPD32 libraries available @echo off :: set the paths below to your install locations set EPDPATH=E:\python27_epd32 set AGSPATH=D:\ArcGIS\Desktop10.1 :: save paths set PATHENV=%PATH% set PPATHENV=%PYTHONPATH% :: set paths for ArcMap w/ EPD set PATH=%EPDPATH%;%EPDPATH%\scripts;%PATH% set PYTHONPATH=%PYTHONPATH%;%EPDPATH%\lib\site-packages echo Starting ArcMap w/ Enthought Python Distribution modules... start /b /d %AGSPATH%\bin ArcMap.exe :: restore paths set PATH=%PATHENV% set PYTHONPATH=%PPATHENV% :: delete variables set EPDPATH= set AGSPATH= set PATHENV= set PPATHENV= ... View more

Help reference - Desktop Help 10.1: Adding an ASCII or text file table See the section: Overriding how text files are formatted ... View more

If I'm understanding this right, the Slice tool could be used to find the area value at a given percentile. However, this could be derailed by the fact that the distribution of cell drainage areas is dramatically skewed, to the point where calculating a area distribution curve may be problematic. Most cells in a basin will have very small drainage areas, and just a few (in the drainage) have very large drainage areas. ... View more

Is the legend property checkbox set to FixedFrame in the ArcMap interface? The help implies that that has to be the case for this property to be used. ... View more

Sorry, I didn't quite get the picture. Clearly the validation is behaving differently inside Model Builder. I don't know why that should be the case. Another thing that bugs me about Field Calculator vs Calculate Field is that in table view when you right click a field and choose Calculate... the interface is similar, but not the same as, the real Calculate Field tool (that it eventually runs). Would be less confusing if they were the same! Both these issues are definitely worth submitting in an incident or on ideas.arcgis.com. ... View more