Latest Contributions by GerryGabrisch

Do you have an interpolated surface as a raster dataset?  If so, you can create the points using one of several python tools out there, then use add surface information tool in the 3D-ArcToolbox to populate the points with the interpolated value from the interpolated raster. ... View more

Or do you want to Locate Features Along Routes after you have created the points? ArcGIS Help 10.1 ... View more

Option 1: I want to create a new point along Hood every 50', so that the point created at the start of Hood is given the value of 566, and the next point is created 50' farther, with a value of 570 or so. So the 570 value would be added to the attribute table and the 570 values is the straight line distance between the newly created interval point and the point of origin of the line? ... View more

Can't you accomplish that task using the Spatial Join tool in ArcToolbox, then export the resulting table to a CSV file? ... View more

Here is an updated version of my Find and Delete Spatially Duplicate Records in a Shapefile.  This code works in ArcGIS 10.1+ and does not require and advance license, it works on any feature class, and it runs much faster relying on geometry and not selections. The geometry compare only uses X and y and not z values, sorry but the method comes that way from ESRI.  There is no consideration of the attribute table but the tool preserves the original FC. try:     print "Find And Delete Spatially Duplicate Records in a Feature Class..."     print "A Two-Bit Algorithms product\nby Gerry Gabrisch (geraldg@lummi-nsn.gov)"     print "Copy?Right! 2015...Free to hack with with attribution."         import sys, traceback, time     import arcpy     from arcpy import env         env.overwriteOutput = True     #the input FC with duplicate records....     inFC= r"C:\gTemp\Converted_Graphics.shp"     #The output FC that will have duplicate records purged....     outFC = r"C:\gTemp\Converted_GraphicsTEST.shp"         arcpy.CopyFeatures_management(inFC, outFC)         with arcpy.da.SearchCursor(inFC, ["OID@", "SHAPE@"]) as cursor:         for row in cursor:             with arcpy.da.UpdateCursor(outFC, ["OID@", "SHAPE@"]) as cursor2:                 dupcount = 0                 for row2 in cursor2:                     if cursor[1].equals(cursor2[1]):                         dupcount +=1                         if dupcount >=2:                             dupcount = 1                             cursor2.deleteRow()     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])     print pymsg + "\n" ... View more

Here is code that will sift though a directory (and all sub directories) and create a convex hull for all GIS data.  The convex hulls are written to a new feature class.  You can then add the feature class to an ArcMap project and see all the data that covers your area of interest.  The output feature class contains an attribute listing the path to the data. try:     print "Create Convex Hulls for All Vector Data..."     print "Created by Gerry Gabrisch (geraldg@lummi-nsn.gov)"     print "Copy?Right! Lummi Indian Business Council 10-1-2013"        import arcpy, math, os, sys, traceback, winsound     winsound.Beep(1000,500)     winsound.Beep(2500,500)         inDir = "Z:\\Data"     appendeddata = "C:\gtemp\DataExtentsZData.shp"        arcpy.env.overwriteOutput = True     datapaths = []        def GetBoundingGeometry(inFC,datapaths):         outFeatureClass = "C:\\gtemp\\convexhulltemp2.shp"         appendeddata = "C:\gtemp\DataExtentsZData.shp"         try:             arcpy.MinimumBoundingGeometry_management(inFC, outFeatureClass,"CONVEX_HULL", "ALL")             arcpy.Append_management(outFeatureClass, appendeddata, "NO_TEST","","")             desc = arcpy.Describe(inFC)             datapaths.append(desc.catalogPath)         except:             pass          for root, dirs, files in os.walk(inDir):         print "working in ", root         #Get all geodatabase feature classes and feature dataset feature classes...         if root.endswith(".gdb"):             arcpy.env.workspace = root             FeatureClasses = arcpy.ListFeatureClasses("*", "All")             for features in FeatureClasses:                 desc = arcpy.Describe(features)                 GetBoundingGeometry(features, datapaths)             datasets = arcpy.ListDatasets("*", "FEATURE")             for dataset in datasets:                 arcpy.env.workspace = root + "\\" + dataset                 FeatureClasses = arcpy.ListFeatureClasses("*", "All")                 for features in FeatureClasses:                     GetBoundingGeometry(features, datapaths)                else:             arcpy.env.workspace = root             FeatureClasses = arcpy.ListFeatureClasses("*", "All")             for features in FeatureClasses:                 desc = arcpy.Describe(features)                 if desc.dataType == "CoverageFeatureClass":                     pass                 else:                     GetBoundingGeometry(features, datapaths)     rows = arcpy.UpdateCursor(appendeddata)     counter = 0     for row in rows:         row.filepath = datapaths[counter]         counter += 1         rows.updateRow(row)               print "done"     winsound.Beep(2500,1000) except arcpy.ExecuteError:     msgs = arcpy.GetMessages(2)     arcpy.AddError(msgs)     print 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)     print pymsg + "\n"     print msgs ... View more

The output file name needs a '.shp' maybe? output_layer = "C:\\Users\\Fred\\Desktop\\GIS\\ModelBuilder Project\\TEST\\test" + str(thectr) +".shp" ... View more

Have you tried making a TIN from the points?  If you have shorelines, use those as hard break lines so the TIN does not "take short cuts. ... View more

Here is an ArcToolbox version of the script above that works in 10.1 and ought to work in any version of 10.1 or later. ftp://lnnr.lummi-nsn.gov/GIS_Scripts/MinimumColorMap.zip ftp://lnnr.lummi-nsn.gov/GIS_Scripts/ ... View more

Kim, you can symbolize any distribution of polygons with a four color map such that no two colors touch.  below is some code that will create a new attribute called fourcolor.  You can symbolized off of the resulting values.  You will have the change the path to the inFC dataset to point to your data. See here for more information....ESRI Mapping Centre - Maps: Los Angeles Ethnicity Here is an ESRI four-color tool but I do not know if it will work with your version of ArcGIS Esri Mapping Centre - ArcGIS Resources Good luck. try:     print"Minimum Color Map"     print "Copyright 2011, Gerry Gabrisch\n"     import arcpy, os, sys, string, traceback, operator             inFC = r"On-Reservation Clam Harvest Areas"         print "Creating a minimum color map for "+ inFC         print "Adding field 'FourColor...'"         listoffields = []     for field in arcpy.ListFields(inFC):         listoffields.append(field.name)     if 'GerryTemp' not in listoffields:         arcpy.AddField_management(inFC, "GerryTemp", "SHORT")     if 'FourColor' not in listoffields:         arcpy.AddField_management(inFC, "FourColor", "SHORT")     del listoffields         counter = 0     rows = arcpy.UpdateCursor(inFC)     for row in rows:         row.GerryTemp = counter         rows.updateRow(row)         counter +=1             listofconnections = []         g1 = arcpy.Geometry()     geometrylist1 = arcpy.CopyFeatures_management(inFC, g1)         g2 = arcpy.Geometry()     geometrylist2 = arcpy.CopyFeatures_management(inFC, g2)         print "Sorting...Please wait..."     counter1 = 0     for geometry1 in geometrylist1:         connection = []         connection.append(counter1)         counter  = 0         for geometry2 in geometrylist2:             #connection.append(counter)                         if (geometry1.touches(geometry2) or geometry1.overlaps(geometry2)) and counter != counter1:                                 connection.append(counter)             counter +=1         counter1 += 1                listofconnections.append(connection)            #sort the list based on the elements (more lists) by length     listofconnections = sorted(listofconnections, key = len, reverse = True)         iscolor1 = [];iscolor2 = [];iscolor3 =[];iscolor4 = []     cantbecolor1 =[];cantbecolor2 =[];cantbecolor3 =[]; cantbecolor4 =[]         counter = 0     for item in listofconnections:         #Set the polygon with the largest number of neighbors to color #1         if counter == 0:             iscolor1.append(item[0])             counter2 = 0             for item2 in item:                 if counter2 >  0:                     cantbecolor1.append(item2)                 counter2 += 1         #For any disconnected polygons, give them a value of 4                elif len(item) == 1:             iscolor4.append(item[0])             counter2 = 0             for item2 in item:                 if counter2 >  0:                     cantbecolor4.append(item2)                 counter2 += 1         #For any remaining polygons, sort out the colors...         else:             if item[0] not in cantbecolor1:                 iscolor1.append(item[0])                 counter2 = 0                 for item2 in item:                     if counter2 >  0:                         cantbecolor1.append(item2)                     counter2 += 1             elif item[0] not in cantbecolor2:                 iscolor2.append(item[0])                 counter2 = 0                 for item2 in item:                     if counter2 >  0:                         cantbecolor2.append(item2)                     counter2 += 1             elif item[0] not in cantbecolor3:                 iscolor3.append(item[0])                 counter2 = 0                 for item2 in item:                     if counter2 >  0:                         cantbecolor3.append(item2)                     counter2 += 1             else:                 iscolor4.append(item[0])                 counter2 = 0                 for item2 in item:                     if counter2 >  0:                         cantbecolor4.append(item2)                     counter2 += 1         counter += 1        print "Writing color codes to the attribute table..."        rows = arcpy.UpdateCursor(inFC)     for row in rows:            if row.GerryTemp in iscolor1:             row.FourColor = 1         if row.GerryTemp in iscolor2:             row.FourColor = 2            if row.GerryTemp in iscolor3:             row.FourColor = 3            if row.GerryTemp in iscolor4:             row.FourColor = 4         rows.updateRow(row)                     arcpy.DeleteField_management(inFC, "GerryTemp")       print "All done."     except arcpy.ExecuteError:     msgs = arcpy.GetMessages(2)     arcpy.AddError(msgs)     print 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)     print pymsg + "\n"     print msgs ... View more

And you can use HTML tags to format the item description for a nicer looking tool. ... View more

I wrote them. ... View more

But you need to have the tabular data converted to vector data firs,t and to do that use Make XY Event Layer in ArcToolbox. 1. Convert your tabular data to vector data. 2. Create your TIN or terrain as you see fit. 3.  Use the geoprocessing tools in ArcToolbox to generate a slope surface. ... View more

Oh, and about the elevation under stream bank problem.  Unless you know where these areas are and surveyed them and created a new surface model based on these data you are out-of-luck and stuck with the original surface model. ... View more

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