Latest Contributions by GerryGabrisch

‎10-14-2011

Here is a portion of the code that returns the error. Again, it works fine as a stand alone where it is taking the inFC as text just like the arcpy.GetParameterAsText(0) try: print "go" import arcpy, os, sys, string, traceback, operator inFC = arcpy.GetParameterAsText(0) g1 = arcpy.Geometry() geometrylist1 = arcpy.CopyFeatures_management(inFC, g1) g2 = arcpy.Geometry() geometrylist2 = arcpy.CopyFeatures_management(inFC, g2) for geometry1 in geometrylist1: for geometry2 in geometrylist2: if (geometry1.touches(geometry2) or geometry1.overlaps(geometry2)): print "Do more stuff here..." print "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

‎10-14-2011

I made a four color map algorithm for Python. When I execute the script in PythonWin, Komodo, or Idle, with the path to the feature class hard-coded into the script, things run fine. I know the syntax is correct. Once I add the script to ArcToolbox it throws an error (error message posted below). The error occurs at this line. if (geometry1.touches(geometry2) or geometry1.overlaps(geometry2)) and counter != counter1: File "E:\SCRIPT~1\MINIMU~1\minimumpolygoncolor.py", line 46, in <module> if (geometry1.touches(geometry2) or geometry1.overlaps(geometry2)) and counter != counter1: Error Info: 'unicode' object has no attribute 'touches' ArcPy ERRORS: I also developed a tool to remap data paths in an mxd. Again, things run fine if I post the code into the Python window. All the broken links are fixed. When I convert this tool to a Toolbox script, the paths get updated...and then removed from the mxd. Does anyone know why a script would act differently once they are brought into ArcToolbox?

‎10-07-2011

I have a USA nationwide dataset of wetland polygons and a local polygon layer of parcel. The CONUS wetlands use NAD_1983_Albers, and the parcels use Wa SP N NAD83. Both feature classes use single part geometries. A GIS user (on an older desktop running XP with 1 GB RAM) tried to clip the wetlands to the parcels and got an error message that no polygons were retured. I tried the same tool on my Windows 7 8GB computer (thinking this was a memory issue) and the tool returned a copy of the CONUS wetlands. ?! I manually selected out the CONUS wetlands for my study area and created a new shapefile from that selection, re-ran the clip tool, and everything works fine. No projections were altered in this process. Any thoughts on why this would happen are appreciated.

‎10-04-2011

Great, that worked thanks! What if there are more than one feature feature datasets in the gdb and each feature dataset has a feature class with the same name?

‎10-03-2011

We have migrated a bulk of our GIS data from shapefiles into feature classes within feature datasets (ArcGIS-file geodatabase-feature datasets). I am trying to develop a Python script to automatically replace the most popular datasets in an mxd with the paths to the new data. My problem is that both the data model (shapefile to feature dataset feature class) has changed, and the feature class name has changed. The replaceDataSource() and replaceWorkspacePaths() methods do not have parameter for connecting to feature datasets. See http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/MapDocument/00s30000000n000000/ The code posted below works if I move my feature class out of the feature dataset within the file geodatabase (meaning the syntax is good). Still, how to I set the replaceDataSource to a feature data set? try: import arcpy.mapping, os, sys, string, traceback mxd = arcpy.mapping.MapDocument("CURRENT") for dataframe in arcpy.mapping.ListDataFrames(mxd): BrokenLinkList = arcpy.mapping.ListBrokenDataSources(dataframe) for item in BrokenLinkList: print item.datasetName if item.datasetName == r"LummiReservationFullBoundary": print "droppin' in" #The path to the new directory... newWorkspace = r"Z:\Data\Boundaries\Administrative\Lummi.gdb\TribalBoundaries" #The new dataset name... datasetName = r"LummiReservation" #This should have a "FEATUREDATASET" options to replace FILEGDB_WORKSPACE. item.replaceDataSource (newWorkspace, "FILEGDB_WORKSPACE", datasetName, False) arcpy.RefreshTOC() #mxd.save() del mxd print "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

‎07-22-2011

Here is the code to populate an attribute in the attribute table with a 0 or 1 for 5% of your records. If you want to select 5% of those records that have some other attribute value, then just insert a select by attribute where I have it marked, then only 5% of the selected records will be assigned a 1. '''Selects a random number of records from a feature class and populates a user-defined field with either a 1 or a 0 based on the input parameters. For example, if percent_to_include is set to 5, then a random ~5 percent of the records will be asigned a value of 1, the rest a value of 0. ''' try: import operator, os, random, traceback, sys, arcpy ##Path to the feature class. fc = r"Z:\GISpublic\GerryG\RestorationandMitigationConcept\20110425IvyPlotMonitoring\transects.shp" ##The attribute to populate with 1 or 0... attribute = r"issample" ##The percent threshold to include as a random selection. percent_to_include = 5 ##INSERT YOUR SELECT BY ATTRIBIUTE HERE uc = arcpy.UpdateCursor(fc) for row in uc: x = random.random() * 100 #arcpy.AddMessage(x) if x <= percent_to_include: row.setValue(attribute, 1) else: row.setValue(attribute, 0) uc.updateRow(row) print "Done" 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]) print pymsg + "\n"

‎07-07-2011

No, you would think! The centroid returns the midpoint of the line if and only if that point falls on the line, otherwise it returns the coordinates of the label point (which seems to be defined by magic:)). I have a script that I tested on a line shaped like three sides of a square, and the centroid is not close to the midpoint of that line. Using the midpoint function from calculate geometry does return the true midpoint coordinates of that line, but there is no Arcpy/Python property that returns that same coordinate value.

‎07-06-2011

I would like to use Python to get the midpoint of a line. There are a number of geometry properties that will return the centroid, the label point�?� http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Polyline/000v000000n2000000/ �?�but what I am looking for is the coordinate values of the position which is exactly halfway along the line�??s length. It is possible to calculate the midpoint using the calculate geometry functions by right clicking a new field in the attribute table and selecting calculate geometry, but I can find no way to get that calculation into a Python code block. Any ideas?

‎06-10-2011

ESRI has posted a script that will do what you want. http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Writing_geometries/002z0000001v000000/ Will a little work you can get this script up and running, your polygons created, and learn a lot about scripting in Python.

‎06-07-2011

Why not just clip the Thiessen polygons to the county border?

‎04-08-2011

This code worked to remove a \n escape sequence from a string, maybe you can alter the code to work for a \u x = somestringwithanescapesequence" x = repr(x) x = x.replace(r'\n', "")

‎03-18-2011

Try this tool http://resources.arcgis.com/gallery/file/map-templates/details?entryID=BC061D5E-1422-2418-7FCE-CC0FBB0C0825

‎02-23-2011

try this... def tvol(a, b, area): if (a >=2): return int(round(a*b/100))

‎02-22-2011

Try using this... OCCAT__tif = "D:/Documentation/draft_12/Corridas/4/ScriptInputs/OCCAT_"+ instep + ".tif" OC_01_tif = "D:/Documentation/draft_12/Corridas/4/ScriptInputs/OC_"+ instep + ".tif" not this OCCAT__tif = "D://Documentation//draft_12//Corridas//4//ScriptInputs//OCCAT_"+ instep + ".tif" OC_01_tif = "D://Documentation//draft_12//Corridas//4//ScriptInputs//OC_"+ instep + ".tif"

‎12-29-2010

You can't. You need to have ArcGIS installed on the computer to take advantage of the geoprocessing functions in arcpy.

Online Status	Offline
Date Last Visited	‎05-15-2026 09:46 AM

My Ideas

Latest Contributions by GerryGabrisch

Re: ArcToolbox scripts behave differently then the hard-coded versions.

ArcToolbox scripts behave differently then the hard-coded versions.

Interesting results from the analysis -clip tool

Re: Replace Data Sources from shapefile to feature dataset feature class in an mxd us

Replace Data Sources from shapefile to feature dataset feature class in an mxd using

Re: Random Polygon Selector by Attribute/Area

Re: Midpoint of a line?

Midpoint of a line?

Re: Corner Points to Polygons...

Re: Run Thiessen polygons to a specific boundary shapefile

Re: Calculate Value tool fails, unicode error from paths with "\U" (rawunicodeescape)

Re: converting gpx file with track data to shapefile for ArcGIS v 10

Re: Rounding decimals up using the field calculator

Re: Problema con Script en Arcgis 10

Re: how can I import arcpy module in python if the application of arcgis not install

Please stop sending me direct mail

User cannot see a directory in Catalog but it is v...

Re: Cannot select or edit elements in Layout view

Re: Creating KML files with Attribute Data

Hey ESRI, your hillshade python code example is wr...