Latest Contributions by ZacharyHart

I've created a simple script to share via a web app so that folks without desktop can convert shapes to KML/KMZ. The script writes the output KMZ as expected in the scratch folder where it belongs but doesn't appear in the GP Results window and therefore is a null output when published as a GP Service. I have a similar script which takes a zip of a multiple shapefiles, reprojects them and zips them back up; this script tool is configured identically and the output zip shows in the GP results and also is returned in the output when published as a GP service. Script Tool: Code (not really sure this will help): import arcpy  import zipfile import sys import os from os import makedirs from os.path import join   inZip = arcpy.GetParameterAsText(0)     # Create a folder in the scratch directory to extract zip to  outFolder = os.path.join(arcpy.env.scratchFolder, "zipContents")  os.mkdir(outFolder) # Create a folder in the scratch directory for the KML KMLfolder = join(arcpy.env.scratchFolder, "KML") makedirs(KMLfolder)   # Extract the zip contents  zip2Extract = zipfile.ZipFile(inZip, 'r')  zip2Extract.extractall(outFolder)  zip2Extract.close() # Set the workspace arcpy.env.workspace = outFolder # Create a list of shapefiles featureclasses = arcpy.ListFeatureClasses() print featureclasses # Use ListFeatureClasses to generate a list of inputs for infc in arcpy.ListFeatureClasses():     # Determine if the input has a defined coordinate system, can't project if it does not     dsc = arcpy.Describe(infc)     if dsc.spatialReference.Name == "Unknown":         print('skipped this fc due to undefined coordinate system: ' + infc)             else:                # Make Feature Layer from FC         tempLayer = arcpy.MakeFeatureLayer_management(infc, "tempLayer")         # Define KML output         baseName=os.path.basename(infc).split('.')[0]         outKML = os.path.join(KMLfolder, baseName + '.kmz')                #Convert Layer to KML         output = arcpy.LayerToKML_conversion(tempLayer, outKML)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

I have a tool which creates a number of points grouped together in 'clusters'. Inevitably, the tool will generate points that fall outside the project boundary. Removing these points is a simple task using select by location. However, I also want to remove any cluster of points that has less than a required threshold value. I'm having some success with using 'collections' much in the same way that Summary Stats would work. However, I'm running into difficulty using an IF statement with the collection....it seems to print anytime the value (regardless of value) is 'collection > X' and prints nothing when ' collection < X' and fails to work when I try ' collection = X'. So it appears the collection is a 'Counter'. Am I on the wrong track here?I couldn't get 'GetCount' to summarize by each value in a field. EDIT: it looks like 'sorted' may be something I can use here, but having a hard time finding examples of this. I want to remove all points in the cluster if the particular cluster contains less than a threshold value. Here is my code so far: def CleanSamplePlots (inPoints, standLayer, outPoints):     #inPoints are the raw output of the desired point generation process     #standLayer is the project stand dataset     #outPoints are the cleaned, validated points     arcpy.MakeFeatureLayer_management(inPoints, "ptLayer")     arcpy.SelectLayerByLocation_management ("ptLayer", "intersect", standLayer,"", "", 'INVERT')     if int(arcpy.GetCount_management("ptLayer").getOutput(0)) > 0:         arcpy.DeleteFeatures_management("ptLayer")     with arcpy.da.SearchCursor("ptLayer", "CLUSTER_ID") as cursor:         count_of_items = collections.Counter(row[0] for row in cursor)         for item in sorted(count_of_items.items(), key=lambda x:x[1]):             if count_of_items < 8:                 print "Cluster:{0} Points{1}".format(item[0], item[1]) Cluster:4 Points4 Cluster:14 Points4 Cluster:6 Points7 Cluster:7 Points8 Cluster:12 Points8 Cluster:21 Points9 Cluster:1 Points9 Cluster:2 Points9 Cluster:20 Points10....etc ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

My brain is fried...if anyone can/wants to modify the title of the post please do so. I'll attempt to post this w/o giving the background of the tool...if that is necessary I'll happily provide it. Here is my code: import arcpy, os, sys, math #Environment settings outGDB = arcpy.env.scratchGDB arcpy.env.workspace = outGDB arcpy.env.overwriteOutput = True #Local Variables #startingPt = arcpy.GetParameterAsText(0) #sampleNumber = int(arcpy.GetParameterAsText(1)) #interDist = int(arcpy.GetParameterAsText(2)) startingPt = "startingPoints" interDist = 660 outName = "cruisePoints" #Create target point FC sr = arcpy.Describe(startingPt).spatialReference outPoint = arcpy.CreateFeatureclass_management(outGDB,outName,"POINT","","","",sr) arcpy.AddField_management(outPoint, "CLUSTER_ID", "TEXT", "", "", "", "", "NULLABLE", "NON_REQUIRED", "") arcpy.AddField_management(outPoint, "PLOT_ID", "TEXT", "", "", "", "", "NULLABLE", "NON_REQUIRED", "") #Open Insert Cursor insert = arcpy.da.InsertCursor(outPoint,["SHAPE@XY","CLUSTER_ID","PLOT_ID"]) fields = ["SHAPE@X", "SHAPE@Y", "OID@"] with arcpy.da.SearchCursor(startingPt, fields) as cursor:     for row in cursor:        #Set Start points        startX = row[0]        startY = row[1]        cluster = row[2]        print(startX,startY)        #Calculate quadrant points        northX = row[0]        northY = row[1] + interDist        print(northX, northY)        eastX = row[0] + interDist        eastY = row[1]        southX = row[0]        southY = row[1] - interDist        westX = row[0] - interDist        westY = startY        start = arcpy.Point(startX, startY)        print(start)        north = arcpy.Point(northX, northY)        east = arcpy.Point(eastX, eastY)        south = arcpy.Point(southX, southY)        west = arcpy.Point(westX, westY)        insert.insertRow([start,row,cluster])        insert.insertRow([north,row,cluster])        insert.insertRow([east,row,cluster])        insert.insertRow([south,row,cluster])        insert.insertRow([west,row,cluster])               del cursor‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ When pasting into the Python Console in Desktop I get the following error (note the third print as well) ... View more