Latest Contributions by GerryGabrisch

Whether the TIN interpolation (straight lines with a constant slope between points) or some raster interpolation method is more accurate is unknown.  The density of your points, the  weights you give them, and the raster interpolation method will all influence the resulting surface model in different ways.  You could assess the quality of your surfaces using root mean square error (http://en.wikipedia.org/wiki/Root-mean-square_deviation).  In order to do that you would need a statistically sound number of surveyed point locations (not your GPS points)  with z values to calculate RMSE. ... View more

I have had good luck speeding up tools with 64-bit background geoprocessing. http://resources.arcgis.com/en/help/main/10.1/index.html#//002100000040000000 ... View more

If you have the bathymetric model of the lake you could use this tool to get a table of areas and volumes at different graduations. ftp://lnnr.lummi-nsn.gov/GIS_Scripts/MultiVolumes10/ Next, use the bathymetric model to generate contour lines at the same graduations.  Find the contour that matches your lake levels and use the multivolumes10 output table to get the volumes. ... View more

This code will populate an existing attribute with a random number.  Run the script.  Then, for example, select values less than 0.2 to select about 20% of your points. #Populates an existing attribute with a randomly generated value try:     import operator, os, random, traceback, sys, arcpy     arcpy.AddMessage("Adding random numbers to an existing attribute.")     arcpy.AddMessage("Copyright 2011, Gerry Gabrisch, gerry@gabrisch.us")     #Path to the feature class.     fc = arcpy.GetParameterAsText(0)     #The Attribute to populate...     attribute = arcpy.GetParameterAsText(1)       uc = arcpy.UpdateCursor(fc)     for row in uc:         x = random.random()         row.setValue(attribute, x)         uc.updateRow(row)     arcpy.AddMessage("Done")     except arcpy.ExecuteError:     msgs = arcpy.GetMessages(2)     arcpy.AddError(msgs)      arcpy.AddMessage(msgs) except:     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)     arcpy.AddMessage(pymsg + "\n")     arcpy.AddMessage(msgs) ... View more

Sorry for the problems.  I just tested the tool and it runs for me. From the link download the toolbox  and the DivideLine.py into a new folder on your computer.  Add the toolbox to ArcToolbox and run like any scripting tool. ... View more

The tool here creates a new point file without altering the original feature geometries. ftp://lnnr.lummi-nsn.gov/GIS_Scripts/CreatePointsAlongALine/ ... View more

Scripts and models can have the parameters coded into them which will give you that message.  Don't you see an OK button on the bottom, and if so, what happens when you click it. ... View more

If you developed a script, and you have hard-coded all the parameters like the path to the data, the attribute to populate....then you can click OK and the script will run. ... View more

A series of map at the same extent each showing a different feature in a feature class? Here is some code that will produce a new map for each water sample sites. Each output map is the same extent and only one sample site appears on each output map.  Each output map has the title changed based on an attribute and each output map is named after the attribute as well.  It would be easy to code this for pdf. The (ham-fisted) code would require some changes...a valid path to your mxd, "WaterSampleSiteMaster" (line 10) is the feature layer in my mxd and SiteID (line 13 and elsewhere) is an attribute in the Water sample site feature class that is used to label the map outputs and the map title. Good Luck.  try:     import arcpy, sys, traceback     print "Working..."         mxdPath = r"Z:\GISpublic\GerryG\WaterResources\WaterSampleSiteAutoMap\WaterSampleSiteAutoMap22.mxd"         mxd = arcpy.mapping.MapDocument(mxdPath)     df = arcpy.mapping.ListDataFrames(mxd)[0]         lyr = arcpy.mapping.ListLayers(mxd, "WaterSampleSiteMaster", df)[0]     print lyr.dataSource     theAttribute = "SiteID"     #GetAListOfEachAttribute     listOfAttributes = []     print "Start cursor"     rows = arcpy.SearchCursor(lyr.dataSource)     for row in rows:         listOfAttributes.append(row.getValue(theAttribute))     for elm in arcpy.mapping.ListLayoutElements(mxd, "TEXT_ELEMENT"):         if elm.name == 'gbglabeler':             elm = elm             break     for item in listOfAttributes:         GBGquery = '"SiteID" =' + '\''+ item + '\''         print  "Working on: ",item         print GBGquery         lyr.definitionQuery = GBGquery         elm.text = "Site: " + item                 outputjpg = "Z:\GISpublic\GerryG\WaterResources\WaterSampleSiteAutoMap\SiteMaps\\" + item + ".jpg"         arcpy.RefreshActiveView()         arcpy.mapping.ExportToJPEG (mxd, outputjpg)     print "done" except arcpy.ExecuteError:     msgs = arcpy.GetMessages(2)     arcpy.AddError(msgs)  #If the error is with the script, or with python, find error and report it to the screen... except:     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) ... View more

Brian, Thanks for the information.  Can you point me in the direction of some online help on the subject? ... View more

Hmmm, Must have been some hang up between Arcpy, Komodo, ArcGIS or Citrix.  I tried to save the original mxd, got an error that the file location folder was in use.  After shutting ArcMap and Komodo down the original code worked fine. ... View more

Thanks for the reply, Matt.  I hope you are doing well.  Yes, SiteID is an attribute storing text.  I tried the field delimeter but that just returned the same text string as the brute force way and nothing showed on the output jpeg.  Any other thoughts? ... View more

This tool will create a perpendicular line at a user defined length from the center ( or end nodes) of the line.  It is not exactly what you want but the code is a starting point to solving your problem. ftp://lnnr.lummi-nsn.gov/GIS_Scripts/createperpendicularlines.zip ... View more