Latest Contributions by JohnDye

‎10-16-2013

They get into this a bit on stack overflow, but I think the issue is you are using triple quotes interchangeably with comments, where a triple quote is really a doc string. Shazam Lucas! I didn't know that triple quotes were docstrings! I always just assumed they were just multi-line comments. After reading the StackOverflow thread and then looking at the actual python documentation referenced within it makes a little more sense. It seems that if you want to use the docstring in this manner, it needs to be within the code block. I'm not sure why using docstrings interchangable with comments would cause an issue though. It doesn't seem to cause any issues within any of my other scripts.

‎10-08-2013

There's no simple function to simply return a list of layers with a selection on them, but you could probably manage this by looping through the layers in the dataframe, getting the layers extent, then getting the layer's selected extent and comparing the two, and if the extents differ, then there is a selection applied to the layer and you can then append it to a list. I haven't tested this and I'm just writing this off the top of my head, but it should give you an idea of where to go with this: mxd = arcpy.mapping.MapDocument("Current") df = arcpy.mapping.ListDataFrames(mxd)[0] TOCLayerList = arcpy.mapping.ListLayers(mxd, "*", df) SelectedLayerList = [] for lyr in TOCLayerList: lyrExtent = lyr.getExtent() lyrSelectedExtent = lyr.getSelectedExtent() if not lyrExtent == lyrSelectedExtent: SelectedLayerList.append(lyr) print SelectedLayerList Again, I haven't tested this or played with it at all, but conceptually this is the best path I can think of to return only those layers with a selection applied. I'm not sure what the Layer method 'getSelectedExtent()' would return if there was no selection applied to a given layer, so you should toy with that as well and figure out what the return is and make a consideration for that in the logic so that you don't get any false positives. You should also note that if you plan to do this on an MXD with a large amount of layers in it, it could take a fair amount of time to complete. Hope this helps.

‎09-27-2013

'params' should be a list object containing your individual parameters. So change the parameter from 'params' to param0, create a list object named params, populating 'param0' in the list and return the list. def getParameterInfo(self): """Define parameter definitions""" param0 = arcpy.Parameter( displayName="Input Features", name="in_features", datatype="DEDatasetType", parameterType="Required", direction="Input") params = [param0] return params I would have thought you could do it the way you did, but apparently if your not getting anything returned, you can't. Easy fix though, just a bit more code.

‎09-26-2013

I reached out to Jason Pardy from ESRI to get the answer and am posting it here as a courtesy in case anyone else finds themself searching for a similar solution. The answer was to create a Polygon from the extent returned by the rectangle_geometry function and to use that polygon in the SBL Intersect. import pythonaddins mxd = arcpy.mapping.MapDocument("Current") df = arcpy.mapping.ListDataFrames(mxd)[0] class ShowMe_This(object): """Implementation for ShowMe.This (ComboBox)""" def __init__(self): self.items = ["Item1", "Item2", "Item3"] self.editable = False self.enabled = True self.dropdownWidth = 'XXXXXXXXXXXXXXXXXXXXXXXXXX' self.width = 'XXXXXXXXXXXXXXXXXXXXXXXXXX' def onSelChange(self, selection): global FeatureGeometry global SelFC if selection == "Item1": FeatureGeometry = "Point" SelFC = r"C:\ArcGIS\MyGeodatabase.gdb\Feature1" elif selection == "Item2": FeatureGeometry = "Polygon" SelFC = r"C:\ArcGIS\MyGeodatabase.gdb\Feature2" elif selection == "Item3": FeatureGeometry = "Polygon" SelFC = r"C:\ArcGIS\MyGeodatabase.gdb\Feature3" else: FeatureGeometry = "" SelFC = "" def refresh(self): self.refresh() class ShowMe_LikeThat(object): """Implementation for ShowMe.LikeThat (ComboBox)""" def __init__(self): self.items = ["Red", "Blue", "Orange", "Green", "Purple", "Pink", "Yellow"] self.editable = False self.enabled = True self.dropdownWidth = 'XXXXXXXXXXXXXXXXXXXXXXXXXX' self.width = 'XXXXXXXXXXXXXXXXXXXXXXXXXX' def onSelChange(self, selection): global Color Color = selection def refresh(self): self.refresh() class ShowMe_ForThis(object): """Implementation for ShowMe.ForThis (Tool)""" def __init__(self): self.enabled = True self.shape = "Rectangle" self.cursor = 3 def onRectangle(self, rectangle_geometry): SelExtent = rectangle_geometry # Get the Extent Extent = rectangle_geometry # Create a Polygon from the Extent Object ExtentPoly = arcpy.Polygon(arcpy.Array([arcpy.Point(Extent.XMin, Extent.YMin), arcpy.Point(Extent.XMin, Extent.YMax), arcpy.Point(Extent.XMax, Extent.YMax), arcpy.Point(Extent.XMax, Extent.YMin)])) # Select the Feature that intersects the Extent Polygon arcpy.management.SelectLayerByLocation(SelFC, "INTERSECT", ExtentPoly, "", "NEW_SELECTION") class ShowMe_Execute(object): """Implementation for ShowMe.Execute (Button)""" def __init__(self): self.enabled = True self.checked = False def onClick(self): # Do whatever you want to do. Technically, you could completely remove this class and just put it all in the ShowMe_ForThis class, so that your logic occurs as soon as the selections are complete.

‎09-24-2013

Hi guys, I'm trying to create a Python Addin to allow a user to dynamically select a feature in the Map Display Window in ArcMap, grab the ID of that feature and populate it in a variable, then select from a second feature class the features within the second feature class which have that ID, exporting those features to a new FC and then symbolizing with a layer file and adding the result to the data frame. Here's what I have so far: import arcpy import pythonaddins mxd = arcpy.mapping.MapDocument("Current") df = arcpy.mapping.ListDataFrames(mxd)[0] class ShowMe_This(object): """Implementation for ShowMe.This (ComboBox)""" def __init__(self): self.items = ["Predictive Trade Area", "Existing Trade Area", "ReDomiciled Customers"] self.editable = False self.enabled = True self.dropdownWidth = 'XXXXXXXXXXXXXXXXXXXXXXXXXX' self.width = 'XXXXXXXXXXXXXXXXXXXXXXXXXX' def onSelChange(self, selection): global FeatureGeometry global SelFC if selection == "ReDomiciled Customers": FeatureGeometry = "Point" SelFC = r"C:\ArcGIS\SYSTEM\COF\Data\MyGeodatabase.gdb\Customers" elif selection == "Predictive Trade Area": FeatureGeometry = "Polygon" SelFC = r"C:\ArcGIS\SYSTEM\COF\Data\MyGeodatabase.gdb\PTAs" elif selection == "Existing Trade Area": FeatureGeometry = "Polygon" SelFC = r"C:\ArcGIS\SYSTEM\COF\Data\MyGeodatabase.gdb\ETAs" else: FeatureGeometry = "" SelFC = "" def onEditChange(self, text): pass def onFocus(self, focused): pass def onEnter(self): pass def refresh(self): pass class ShowMe_LikeThat(object): """Implementation for ShowMe.LikeThat (ComboBox)""" def __init__(self): self.items = ["Red", "Blue", "Orange", "Green", "Purple", "Pink", "Yellow"] self.editable = False self.enabled = True self.dropdownWidth = 'XXXXXXXXXXXXXXXXXXXXXXXXXX' self.width = 'XXXXXXXXXXXXXXXXXXXXXXXXXX' def onSelChange(self, selection): global Color Color = selection def onEditChange(self, text): pass def onFocus(self, focused): pass def onEnter(self): pass def refresh(self): pass class ShowMe_ForThis(object): """Implementation for ShowMe.ForThis (Tool)""" def __init__(self): self.enabled = True self.shape = "Rectangle" self.cursor = 3 def onRectangle(self, rectangle_geometry): global SelExtent SelExtent = rectangle_geometry class ShowMe_Execute(object): """Implementation for ShowMe.Execute (Button)""" def __init__(self): self.enabled = True self.checked = False def onClick(self): # Establish a reference to the appropriate Layer File based on the users selections LYRfile = r"C:\ArcGIS\SYSTEM\COF\Data\Layer Files\ShowMe" + "\\" + FeatureGeometry + "\\" + Color + ".lyr" print LYRfile if FeatureGeometry == "": pythonaddins.MessageBox("You must select a Feature from the 'Show Me This' drop-down menu.", "ERROR", 0) elif Color == "": pythonaddins.MessageBox("You must select a color to symbolize the feature from the 'Like That' drop-down menu.", "ERROR", 0) elif FeatureGeometry == "" and Color == "": pythonaddins.MessageBox("You must select both a Feature from the 'Show Me This' drop-down menu and a color from the 'Like This' drop-down menu.", "ERROR", 0) else: pass # How do I select the feature from 'SelFC' that falls within the extent provided by 'SelExtent'??? Problem is, I can't figure out how to use the extent generated by 'SelExtent', which is really just the 'rectangle_geometry' produced by the tool, to select the particular feature from 'SelFC' that falls within the extent. I'm thinking maybe I should use onCircle instead as that produces a polygon object I might be able to do a Select By Location with for an Intersect, but thought I would see if any of the big brains here knew better... Thanks in advance

‎09-23-2013

Because Python Toolboxes can have multiple tools, all contained within the .pyt file, and each tool can have it's own documentation xml contained within the tool's xml file, it would be a nice improvement to see the My Toolboxes category in ArcToolbox support organizing those tools within Folders. Example: I have 2 Python Toolboxes, both loaded into the My Toolboxes directory on my system so that they will be visible in the My Toolboxes section in ArcCatalog. As you can see, the Location Analytics toolbox is hosted at the root of the My Toolboxes directory, while the 'GO' toolbox is hosted within a subfolder. However, because the GO toolbox is hosted within a subfolder, the My Toolboxes category in ArcToolbox does not recognize that it is present and as a result does not populate the 'GO' toolbox in ArcToolbox. I would like to be able to organize custom tools into folders in the My Toolbox directory, because all of my tools have accompanying documentation and many of my toolboxes contain multiple tools, resulting in a kind of jumble of files in the My Toolbox directory which could be better organized through the implementation of folders, instead of ending up with this: Alternatively, support for Shortcut links would also suffice as I could organize and store the tools elsewhere, and just place shortcut links in the My Toolbox directory which link back to the actual files so they can be organized into a folder.

‎06-20-2013

I have a Table and a Feature Class. I'm trying to create a Selection on the FC based on the Selection in the Table, and export the resulting Selected Features in the FC to a new FC, but I keep getting an invalid expression error. #Create the Variables BlockGroups = r"C:\ArcGIS\Business Analyst\US\Demographic Data\BA_2011.gdb\esri_bg" BG2Branch_ODCM_MDT = r"C:\ArcGIS\PTA_Test.gdb\CBSA_12420_ProForma_ODCM_MDT" #Make a Feature Layer from the Block Groups arcpy.MakeFeatureLayer_management(BlockGroups, "BlockGroups_lyr") # Create unique set of StoreIDs StoreSet = set([r[0] for r in arcpy.da.SearchCursor(BG2Branch_ODCM_MDT, ["DID"])]) # Create a for loop to iterate through the unique set of StoreIDs for Store in StoreSet: # Make a TableView of the Table using a Where Query to make a Table View only of the records where the 'DestID' is equal to the # current Store TableView = arcpy.MakeTableView_management(BG2Branch_ODCM_MDT, "BG2Branch_ODCM_MDT_tView", '"DestID" = ' + "'" + str(Store) + "'") TableView = TableView.getOutput(0) # Select from the BlockGroups Feature Layer, the records which have an 'ID' that matches the values held in the 'OriginID' field # in the current selection of the 'BG2Branch_ODCM_MDT' table. arcpy.SelectLayerByAttribute_management("BlockGroups_lyr", "NEW_SELECTION", '"ID" IN (SELECT "OriginID" FROM ' + str(TableView) + ")") # Export the selected Block Groups to a new Feature Class. arcpy.FeatureClassToFeatureClass_conversion("BlockGroups_lyr", r"C:\ArcGIS\SYSTEM\COF\Data\CapitalOne.gdb", "TradeArea_" + str(Store)) I keep receiving the error Runtime error Traceback (most recent call last): File "<string>", line 19, in <module> File "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\management.py", line 6435, in SelectLayerByAttribute raise e ExecuteError: ERROR 000358: Invalid expression Failed to execute (SelectLayerByAttribute).

‎06-14-2013

If executing the 'arcpy.AddField_management()' function against a File Geodatabase, the 'Scale' and 'Precision' parameters are ignored. This means that in order to set these properties, you have to create an empty feature class, and then use 'arcpy.FeatureClassToFeatureClass_conversion' and create a series of Field Maps in order to set the field properties. Not only is that inefficient (I have to create two feature classes, then go back and delete the first) but Scale and Precision are very basic properties of a Geodatabase field. FeatureClassToFeatureClass_conversion should not be required just to set such basic properties, I should not need to type like 10 lines of code or possibly more, depending on how many fields I need to add, just to be able to set 'Scale' and 'Precision' for a File Geodatabase field. You need to fix this like...yesterday, ESRI.

‎06-13-2013

Rhett, Yes, I tried executing Table to Table directly against the result layer ('outNALayer') and it fails. Here's why: The result layer object ('outNALayer') from the OD Cost Matrix is a Network Analyst Composite Layer which is made up of multiple sublayers, each of which references a particular Network Analyst object (Origins, Destinations, Point Barriers, Line Barriers, Polygon Barriers & Solved Route/Lines) used in the Network Analysis Solve function. See this article: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//004700000034000000 [INDENT][/INDENT]"The network analysis layer has its own properties, even though it is a composite layer made up of subordinate feature layers that in turn have their own sets of properties. Some of the network analysis layer properties are crucial to defining the analysis you want to perform. Properties include the impedance to be used and the restrictions to be obeyed, as well as the U-turn policy, output shape type, accumulation attributes, and properties for finding network locations. Additionally, each kind of analysis has specific properties of its own." The Solved Route ('Lines') is what I was after once the OD Cost Matrix completes because I wanted the distance from the origin to the destination. See this article: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//00470000004r000000#GUID-823C3FF5-FE11-4D30-B50E-4B3CBBC05933 [INDENT][/INDENT]"The Lines network analysis class stores information about the resulting paths from the origins to the destinations. As with other feature layers, its symbology can be accessed and altered from its Layer Properties dialog box. The Lines class is an output-only class; it is empty until the analysis is complete. Since OD cost matrix analyses tend to result in a multitude of lines, the line objects are not displayed in the Network Analyst window; rather, only a count of the lines is displayed there. You can inspect the individual lines in the attribute table of the Lines feature class. When you choose to display lines in the OD cost matrix analysis layer, they appear as straight lines. They are generalized. However, their cost attributes always report the least-cost network path." The issue is that there is no way to access the Composite Layer ('outNALayer')'s sublayers (Origins, Destinations, Point Barriers, Line Barriers, Polygon Barriers & Solved Route/Lines) outside of a dataframe. This is why the only solution was to load the Composite Layer ('outNALayer') into an MXD, reference the dataframe, and extract from there. You cannot access the sublayers through any arcpy function that I've been able to find. You have to go through the dataframe to get to those objects/sublayers.

‎06-12-2013

This is commonly done with the Summary Statistics tool Darren, That's exactly what I was looking for! 700+ tools makes it hard to find something when you don't know the name! Thanks!!

‎06-12-2013

This is commonly done with the Summary Statistics tool Darren, That's exactly what I was looking for! Thanks!!

‎06-12-2013

Hey Rhett, I pulled the Python snippet through arcmap and this was it returned: # Replace a layer/table view name with a path to a dataset (which can be a layer file) or create the layer/table view within the script # The following inputs are layers or table views: "BG2Branches\Lines" arcpy.TableToTable_conversion("BG2Branches/Lines","C:/Users/jdk588/Documents/ArcGIS/Default.gdb","ODCM_Test","#","""Name "Name" true true true 128 Text 0 0 ,First,#,/ODLines,Name,-1,-1;OriginID "OriginID" true true true 4 Long 0 0 ,First,#,/ODLines,OriginID,-1,-1;DestinationID "DestinationID" true true true 4 Long 0 0 ,First,#,/ODLines,DestinationID,-1,-1;DestinationRank "DestinationRank" true true true 4 Long 0 0 ,First,#,/ODLines,DestinationRank,-1,-1;Total_Length "Total_Length" true true true 8 Double 0 0 ,First,#,/ODLines,Total_Length,-1,-1""","#") As you can see, it still references the 'BG2Branches\Lines' which is accessible in ArcMap but not ArcCatalog. The issue is that the the BG2Branches is some sort of composite NA Layer, and is not written to disk anywhere. It's like a Group Layer of some sort, and there's no real method to get the sublayers without an instantiation in a dataframe. So I have figured out a solution that works without having to run the tool from ArcMap. Basically, I just save the result to a layer file, then use the mapping module to leverage an ArcGIS Template MXD and load the layer into the Template MXD's Dataframe. One that's done, I can extract the table from the Template MXD's Dataframe without ever having to actually open the MXD. After the 'TableToTable_conversion' is complete and I have a copy of the table on disk, I remove the layer from the Template MXD and then delete the actual layer file. It takes about 1-2 minutes for all of that to occur, so not too shabby. arcpy.na.Solve(outNALayer) arcpy.SaveToLayerFile_management(outNALayer, workspace + "\\" + outLayerFile, "RELATIVE") mxd = arcpy.mapping.MapDocument(r" C:\Program Files (x86)\ArcGIS\Desktop10.1\MapTemplates\Traditional Layouts\LetterPortrait.mxd") df = arcpy.mapping.ListDataFrames(mxd)[0] addlyr = arcpy.mapping.Layer(r???C:\Temp\BG2Branches.lyr") arcpy.mapping.AddLayer(df, addlyr) for lyr in df: if lyr.name == "Lines": arcpy.TableToTable_conversion(lyr, addlyr, "ODCM_Table") arcpy.mapping.RemoveLayer(df, lyr) arcpy.Delete_management(workspace + ???\\??? + str(addlyr) + ???.lyr???) else: arcpy.AddError (???ODCM Lines Not Found.???)

‎06-12-2013

I could be wrong, but I think your answer is here. "When a process ends prematurely or crashes, .lock files may be temporarily left behind in the geodatabase folder. ArcGIS eventually removes these files in future sessions as new locks are taken. In the meantime, such files do not continue to lock data, and as they take up no disk space, removing them provides no benefit. " In my experience, I've found this isn't neccesarily true. Especially when working with a Geodatabase that is on a Network Share. Numerous times I have tried to access an FC that someone else processed against, only to find that they are locking the file even though they don't even have ArcMap open (Note: If you look at the Geodatabase in Windows Explorer, the name of the Lock file reveals the computer the lock came from.) As a result, the first and last thing all of my scripts do is compact the Geodatabase in order to ensure no that no residual locks prevent the script from running and that no locks remain in the GDB upon exit. "The geoprocessing Compact tool and the copy and paste operations, available in the Catalog tree, also delete unused .lock files. The geoprocessing Compact tool compacts the geodatabase, removing all inactive .lock files in the process. Copying and pasting a file geodatabase removes all inactive .lock files from the source geodatabase before copying data to the new geodatabase." Hope that helps.

‎06-11-2013

I don't see where you are assigning ODCM_LYR to anything, is that just missing from the posted code. Hey Rhett, Thanks for following up. I appreciate the efforts. Sorry, I should have included the reference. In that particular test, 'ODCM_LYR' was referenced to the result object using: ODCM_LYR = arcpy.na.Solve(outNALayer) ODCM_LYR = ODCM_LYR.getOutput(0) also, if you are trying to convert from the numpyarray to FC, wouldn't that be NumPyArrayToFeatureClass instead of FeatureClassToNumPyArray ? In this test, I wasn't trying to convert from a numpy array to FC, I was trying to get the Lines NA Class into a Numpy Array and then follow that up with a NumpPyArraytoTable to get the table. I never got it into the array though =( Currently, I'm working around all of this by running the tool out of Arcmap and extracting the table directly using: arcpy.TableToTable_conversion("BG2Branches\Lines", workspace, str(cbsaFD) + "_ODCM") and it works just fine...but I dont think I should be forced to run this out of ArcMap only. It should be executable through Catalog or as a Standalone Script as well.

‎06-11-2013

Is there a BG2Branches Feature Class (not a .lyr)? Nope. That's the point. This is supposed to produce that Feature Layer which I can then convert to an FC.

Online Status	Offline
Date Last Visited	‎11-11-2020 02:23 AM

My Ideas

Latest Contributions by JohnDye

Re: Is it me or the IDE??

Re: Determining Selected Layers in the Table of Contents with ArcPy

Re: Tool has no parameter - Can't create python toolbox parameters ?

Re: Dynamically Select Feature in Map with Python Addin

Dynamically Select Feature in Map with Python Addin

Support Folders inside of User's 'My Toolboxes' directory

Exporting Subselection based on Selection in different table

Support 'Scale' and 'Precision' parameters 'arcpy.AddField_management' in File Geodatabase

Re: Convert Result Object to GDB Table

Re: Calculate Field Sum for Unique Set

Re: Calculate Field Sum for Unique Set

Re: Convert Result Object to GDB Table

Re: Mystery locks on feature classes

Re: Convert Result Object to GDB Table

Re: Convert Result Object to GDB Table

Re: Authenticating to Portal

Getting 498: Invalid Token from AGS Server after p...

Re: ArcGIS Data Store Under the Hood

Re: ArcGIS Data Store Under the Hood

Re: ArcGIS Data Store Under the Hood