Latest Contributions by GerryGabrisch

‎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.

‎12-28-2010

Having suffered enough frustrations with 10 bugs and crashes, I decided to finally reinstall 9.3.1 throughout our department. Since we have many mxds as 10 mxds I thought I could write a nice little script to convert them all to 9.3.1 mxds and get on with work using 9.3.1. Unfortunately, I too, get the �??this map created in a later version of ArcGIS�?? message. Too bad this function, like so much of 10, did nothing but waste my time. Maybe, if this method works sp2, and someone is still tired of 10, here is the code to convert all the v10 mxd to 9.3 mxds in a new directory. [HTML] #Convert all the v10 mxds in TheFolder, and any subdirectory in TheFolder, to #an earlier version of mxd to CopyToHere directory. try: import arcpy, os, sys TheFolder = r"Z:\GISpublic\GerryG\ATLAS\Copyof2010LummiAtlas" CopyToHere = r"Z:\GISpublic\GerryG\ATLAS\Copyof2010LummiAtlas\93mxds\\" for root, dirs, files in os.walk(TheFolder): for file in files: if file.endswith("mxd") or file.endswith("MXD"): print root print file mxd = arcpy.mapping.MapDocument(root + "\\"+ file) #Only make copies of mxd that have a dateSaved value (v10 only) if mxd.dateSaved: print mxd.dateSaved mxd.saveACopy(CopyToHere + file, "9.3") print "Done!" sys.exit() except: print arcpy.GetMessages(2) print ErrorDesc.message print "Error" sys.exit()[/HTML]

‎12-28-2010

I am not sure if this is your problem but when you are finished with a cursor be sure to use something like del row, rows By not deleting the row and rows (or whatever you name them) you can get schema locks when you try to run another cursor on the same data.

Online Status	Offline
Date Last Visited	4 weeks ago

My Ideas

Latest Contributions by GerryGabrisch

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

ArcPyMapping saveACopy() to previous version

Re: Nested Loop - Site Comparison

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...