Latest Contributions by GerryGabrisch

Paul, Here is a link to a number of tool to create relative surface models.  There are tools to detrend by transects, Euclian planes, Voronoi polygons... I would give detrend by transect a try first.  White if you have troubles or questions. ftp://lnnr.lummi-nsn.gov/GIS_Scripts/DetrendingTools.zip ... View more

James, Assuming you have all the tiffs in one directory, and assuming the tiff file names are of the same length, you could sift through all the files in the folder, find which files were winter months, create a list of those file names, and then execute you code after the list of file names is complete.  This gets the file names, splits and slices out the month, and writes the file name to a new list. import  os indir = ""#YOUR DIRECTORY OF TIFFS #Create an empty list to store all the file names for winter months.... listofdatatoprocess = [] #Get the files in the directory.... for root, directories, files in os.walk(indir):     for filename in files:       filenamesplit = filename.split(".")         #extract the month...          month = filenamesplit[0][31:33]       if month == '11' or month == '12' or month == '01' or month == '02' or month == '03':             listofdatatoprocess.append(filename) for item in listofdatatoprocess:         #YOUR CODE HERE ... View more

But if you wanted to code a Boolean check box in a scripting tool using Python and ArcToolbox you would do something like this. myboolean = arcpy.GetParamterAsText(0) When you add the script to a ArcToolbox scripting tool you would define the parameter as type Boolean. Depending on whether the checkbox is checked or not, myboolean would have a value of "true' or 'false'. Then you could do something like if myboolean =='true':     your code here... ... View more

Dan, Thanks for the reply.  I am not trying to change the underlying data, just symbolize it in a different. In the grid layer properties I can see that the min and max values are between 0 and 1.  My question is, when I try to symbolize those data by classes I see a different range of values.  Why? ... View more

I had the same problem using 10.1 going from tables in a gdb to feature dataset feature classes.  Changing the output to something other than the name of the input table solved my two days of head scratching.  Thanks! ... View more

What do you mean by 'not regular'?  Are the cells not all perfect squares, is the distribution of the cells not rectangular, or are there cases where the 3 x 3 groupings will leave left-overs?  If the  cells are squares and the distribution of squares lends itself to 3 x 3 groupings across the surface then... Use the create fishnet tool to generate a feature class of cells that are 1500m x 1500m.  Then use a spatial join to populate the original cells with the OID of the 1500 x 1500 cells. ... View more

I have been able to successfully register folders using the manager for ArcGIS for Server 10.1.  Thanks to everyone that helped by responding with helpful posts. I tried to register a new folder today using the same methods and parameters that have worked in the past only today I get a message that the folder is being registered (and the spinning 'please wait bar') but nothing happens.  The spinning bar just spins! What would prevent me from registering a directory today when I have used the same methods in the past? ... View more

Thank you for your suggestions, but these usual ways don't work here on tidal flats (see attached) with low slopes. The contour lines make differences vizualising the river direction, but the method surface-slope not really. With this method I can't consider the direction. Also reclassification doesn't bring a benefit. Low slopes shouldn't make a difference.  Convert the contours to a surface model using the advice posted by others like topo to raster.  Generate a slope raster from the surface model.  Use focal statistics in the Spatial Analyst toolbox to return a new surface with the average slope over some user-defined  area.  Reclassify the focal statistics output. It also sounds like you want to view the direction of flow.  Is that right?  If you want to see which way the water would flow across your surface then generate a flow accumulation surface from your surface model  using the Spatial analyst /hydrology tools.  After you have created the flow accumulation surface, go into the flow accumulation symbology, change the gamma stretch to something like 20 and set the standard deviation to 0.5.  This will really bring out the flow directions. I hope this helps. ... View more

An alternative might be to create points along the line using this geoprocessing tool: ftp://lnnr.lummi-nsn.gov/GIS_Scripts/CreatePointsAlongALine/ Make the points with a 10 meter spacing.  The output features will be points on your line attributed with distance from the start node.   The tool even would with lines with bends. Next, this tool will populate the resulting point file with a new attribute with a modulo value.  ftp://lnnr.lummi-nsn.gov/GIS_Scripts/Modulo/ Given your example you calculate a modulo value of 20 on the distance values.  Using the new modulo values, select any point with a modulo other than zero and you will get a selection of the points at 10m, 30m, 50m from the start node. Of course, this works for one line of length 200m, or for all the lines in your line feature class with length of 200m.  If you have many lines of different lengths, and you have Python or Modelbuilder skills, you could iterate the lines, get the lengths, calculate the division disance, and run the process using the values from each iteration. Good luck! ... View more

Thanks, hua17! That is a more elegant solution than my hack as shown below.  infc =r'C:\MultiPartFeatureClass.shp' outfc = r'C:\FisrtPartFeatureClass.shp' newarray = []     for row in arcpy.da.SearchCursor(infc, ["OID@", "SHAPE@"]):         partnum = 0         for part in row[1]:             if partnum == 0:                 thispart = []                 for pnt in part:                     if pnt:                         thispart.append([pnt.X, pnt.Y])             partnum += 1         newarray.append(thispart)     features = []     for feature in newarray:         features.append(arcpy.Polyline(arcpy.Array([arcpy.Point(*coords) for coords in feature])))     arcpy.CopyFeatures_management(outfc, SingleParts) ... View more

Yes, unless you changed the default z factor parameter to something other than 1. ... View more

Try running this version with an error handler and see if you get an error message to post. # --------------------------------------------------------------------------- # NeighborScore.py # Created on: 2013-01-07 13:26:35.00000 # Updated on: 2014-02-12 # Description: Will create output tables based on the average and sum of neighboring polygon values # --------------------------------------------------------------------------- #B comment: I believe this is as the script was delivered on Feb 12 2014 #B comment: See most recent NeighborScore_Working version for most current version try:         import sys     import traceback     import scipy.ndimage     import numpy,arcpy     import string     import PIL.Image     arcpy.env.overwriteOutput = True     inpRaster=arcpy.GetParameterAsText(0)         desc = arcpy.Describe(inpRaster)     layersource = desc.catalogPath         imageread=scipy.ndimage.imread(layersource)     sobelimage=scipy.ndimage.sobel(imageread,0)         outimage=PIL.Image.fromarray(sobelimage)     Out=outimage.save('D:\Outputs\sobel.jpg')         except arcpy.ExecuteError:     print "error"     msgs = arcpy.GetMessages(2)     arcpy.AddError(msgs)      print msgs except:     print "error"     tb = sys.exc_info()[2]     tbinfo = traceback.format_tb(tb)[0]     pymsg = "PYTHON ERRORS:\nTraceback info:\n" + tbinfo + "\nError Info:\n" + str(sys.exc_info()[1])     msgs = "ArcPy ERRORS:\n" + arcpy.GetMessages(2) + "\n"     arcpy.AddError(pymsg)     arcpy.AddError(msgs)     print pymsg + "\n"     print msgs ... View more