Latest Contributions by RandyBurton

I think a sort is required between line 7 and 8.  Interesting code, James.  Thanks for sharing. # line 7 df = df.sort(['area_id','length'],ascending=[True,False]) # line 8‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

I am still getting the folder instead of the feature with 10.5.0.  I added a few lines to check the working directory for Python and found if I change the cwd to the arcpy workspace then describe finds the feature.  Perhaps there's a problem with the Python path/environment? print "Arcpy: {}".format(arcpy.env.workspace) print "Python: {}".format(os.getcwd()) for fc in arcpy.ListFeatureClasses("*"):         desc = arcpy.Describe(fc)         print "{} - {}".format(fc, desc.dataType) os.chdir(arcpy.env.workspace) print "\nNow: {}".format(os.getcwd()) for fc in arcpy.ListFeatureClasses("*"):         desc = arcpy.Describe(fc)         print "{} - {}".format(fc, desc.dataType) # Output (using IDLE): Arcpy: C:\Temp\describe\Default.gdb Python: C:\Temp\describe point_fc - Folder line_fc - FeatureClass polygon_fc - FeatureClass Now: C:\Temp\describe\Default.gdb point_fc - FeatureClass line_fc - FeatureClass polygon_fc - FeatureClass‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Try this one: import arcpy arcpy.env.overwriteOutput=True testVector = r"C:\Users\adam\Desktop\testarcpy\cont_int.shp" result = arcpy.da.SearchCursor(testVector, ["glac_id","length"],                                where_clause = "glac_id = 'G012039E78878N'",                                sql_clause = (None, 'ORDER BY length DESC')                                ).next() print result[0] # glac_id print result[1] # length ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Try this for the where clause (double quotes around entire where clause, single quotes only around the id code): where_clause = "glac_id = 'G012039E78878N'",‍‍‍ ... View more

In your original question, you were wanting to "group by" area_id.  My thoughts were to use the 3 area_id values to feed  the where_clause in 3 SearchCursors.  You would use only the first row of each search. ... View more

I did try this on a computer running ver 10.2.1, and describe gave the correct results.  I will try again with the 10.5.0 version later today. ... View more

In my latest test I was using "point_fc" for a folder name, with the same name used for a feature inside the geodatabase.  I first noticed it when I moved a folder using the same name as a polygon feature that was part of a dataset into the directory with the file geodatabase.  I was using describe to look at a feature's shapeType property - which is not a property when describing a folder! ... View more

Version is ArcMap 10.5 ... View more

Perhaps something like (to get only first row): OID = arcpy.da.SearchCursor(myFile, ["OBJECTID","length"],                               where_clause = "area_id = 'some_code'",                               sql_clause = (None, 'ORDER BY length DESC')                               ).next()[0] # use [0] to return just OBJECTID, remove for tuple‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Are you using a file geodatabase? Or an SQL server?  Some commands do not work with a file geodatabases. TOP is only supported by SQL Server and MS Access databases. See: SearchCursor ... View more

There are several ways to work around this issue.  Renaming the folder is one option.  Or this adjustment to the code will also work: for fc in arcpy.ListFeatureClasses("*"):         desc = arcpy.Describe(gdb+'\\'+fc) # include gdb/workspace ‍‍ I was expecting Describe to use the declared workspace in a stand-alone script.  But, apparently, it starts with the folder that contains the geodatabase.  This is an example where scripts work slightly different between ArcMap's Python window and a stand-alone script. Using hasattr to check for properties is a good suggestion. ... View more

Have you considered the Add Geometry Attributes tool (Data Management > Features)?  It will add a Point_X and Point_Y to an attribute table and populate it with the geometry values.  You will need to select GCS_WGS_1984 as the coordinate system, and the tool will populate the fields with longitude and latitude. ... View more

Using ListFeatureClasses to loop through a geodatabase: import arcpy import os gdb = r'C:\Path\To\geodatabase.gdb' arcpy.env.workspace = gdb sfDir = 'C:/Temp/' # directory to put shapefiles, subdirectories will be created for fc in arcpy.ListFeatureClasses("*","Polygon"):     print "Creating shapefile: {}".format(fc)     if not os.path.exists('{}{}'.format(sfDir, fc)):         os.makedirs('{}{}'.format(sfDir, fc))     arcpy.FeatureClassToShapefile_conversion(fc, '{}{}'.format(sfDir, fc))     arcpy.AddField_management('{}{}/{}.shp'.format(sfDir,fc,fc),                               field_name = "LABEL",                               field_type = "TEXT",                               field_length = "25" )     arcpy.CalculateField_management('{}{}/{}.shp'.format(sfDir,fc,fc),                                     field = "LABEL",                                     expression = '"{}"'.format(fc) )‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ This puts shape files in a subdirectories to keep them grouped. It uses the feature's name for the subdirectory, the shape file name and to populate a "LABEL" field. This assumes that the feature name will be valid for such purposes. It has not been tested with datasets, but a version using ListDatasets as an outer loop could be created if this doesn't pick up all the desired features. ... View more