Latest Contributions by CCWeedcontrol

I have this polygon that I need to cut into pieces, 10 horizontal and 10 vertical blocks. My though was converting the polygon to polylines, (pic 1) then break the polyline into 4 parts,  then creating the 10 horizontal and 10 vertical lines, then converting those into polygons. The issue I am having is that when creating the segments, the west and south boundary have 20 segments. Not know why 20 segments were being created for the west and south boundary. I did a arcpy.management.FeatureVerticesToPoints and that it was creating two segments for the west and south polylines (pic 2). I am not sure how this can happen given that the polyline comes from the polygon.       import arcpy # Define your workspace and input feature class arcpy.env.workspace = r"C:\Temp\Sec.gdb" input_fc = "Sec_1N2W30" output_fc = "Polygons_Output" # This will store the polygons segment_fc = "Segmented_Polylines" # This will store the segmented polylines # Convert features to polygons arcpy.management.FeatureToPolygon(input_fc, output_fc) # Create an empty Polyline feature class for segmented polylines arcpy.CreateFeatureclass_management(arcpy.env.workspace, segment_fc, "POLYLINE") # Create a cursor to insert segmented polylines with arcpy.da.InsertCursor(segment_fc, ["SHAPE@"]) as insert_cursor: # Iterate over each polygon in the output_fc with arcpy.da.SearchCursor(output_fc, ["SHAPE@"]) as search_cursor: for row in search_cursor: polygon = row[0] # Get the exterior ring of the polygon exterior_ring = polygon.boundary() # Iterate over each part of the exterior ring for part in exterior_ring: if len(part) < 2: continue for i in range(len(part) - 1): p1 = part[i] p2 = part[i + 1] line = arcpy.Polyline(arcpy.Array([p1, p2])) length = line.length segment_length = length / 10 # Create points at each segment interval points = [p1] for j in range(1, 10): # 9 additional points to create 10 segments segment_point = line.positionAlongLine(j * segment_length) points.append(segment_point.firstPoint) points.append(p2) # Ensure no duplicate points by slightly adjusting them unique_points = [] for point in points: if not unique_points or unique_points[-1] != point: unique_points.append(point) # Create polylines from the unique points and insert them for k in range(len(unique_points) - 1): segment_line = arcpy.Polyline(arcpy.Array([unique_points[k], unique_points[k + 1]])) insert_cursor.insertRow([segment_line]) print("Segmented polylines created successfully in", segment_fc)       ... View more

I'm attempting to execute a straightforward select by location operation using "Have their center in" with Pro. I initially thought this was comparable to ArcMap's Select by Location with the condition "Have their centroid in the source layer feature." However, when I try to perform the same operation in ArcGIS Pro using "Have their center in," I encounter the following error: "Error 160306: The server was not found." My objective is to select features from Layer B whose centroid falls within the boundaries of polygons in Layer A. ... View more

I have a feature class and a table that I need to join. The table and feature class share a common field, The table will have duplicate Acc (ID's), and the Feature class will have just one Acc (ID's). Both are in the same file geodatabase, both have objectIds. The following doesn't crate the one to many join, what am I missing?     arcpy.env.workspace = "C:/Temp/test.gdb" fc = "C:/Temp/test.gdb/TaxPs" sourceTB = "C:/Temp/test.gdb/table" arcpy.MakeFeatureLayer_management(fc, "TaxPs") arcpy.AddJoin_management("TaxPs", "Acc", sourceTB, "Acc", "KEEP_COMMON") arcpy.management.CopyFeatures(fc, os.path.join(arcpy.env.workspace, "Taxps_Joined"))       ... View more

I need to copy features that are not Null. I need to use the arcpy.management.CopyFeatures because I am doing a join, that is one to many. The issue I am having is that my query to select "IS NOT NULL" doesn't work. I don't get an error, what is the correct syntax for this query?   Source = "Parcels" Lyr = "Par" fieldName = "TAXID" select = """{0} IS NOT NULL""".format(fieldName) arcpy.management.MakeFeatureLayer(TAXSource, Lyr) arcpy.SelectLayerByAttribute_management(Lyr, "NEW_SELECTION", select) arcpy.management.CopyFeatures(Lyr, "C:Temp/test.gdb/TestPar")   Lyr = "TaxPar" #fieldName = "TAXID" select = "TAXID IS NOT NULL" #"'TAXID = IS NOT NULL'" arcpy.management.MakeFeatureLayer(TAXSource, Lyr, select) arcpy.management.CopyFeatures(Lyr, "C:Temp/test.gdb/Testpar")   ... View more