Latest Contributions by 2Quiker

The best I came with is something like this,  import arcpy from collections import defaultdict, deque # ------------------------------------------------ # ENVIRONMENT # ------------------------------------------------ arcpy.env.overwriteOutput = True arcpy.env.addOutputsToMap = False # ------------------------------------------------ # INPUTS # ------------------------------------------------ input_points = "Centroids_2_10" output_fc = r"C:\Temp\ClusterProject.gdb\Groups_Of_15" search_distance = "1320 Feet" target_size = 15 group_field = "GROUP_ID" # ------------------------------------------------ # 1. COPY INPUT # ------------------------------------------------ print("1. Copying input points...") arcpy.management.CopyFeatures(input_points, output_fc) oid_field = arcpy.Describe(output_fc).OIDFieldName # ------------------------------------------------ # 2. GENERATE NEAR TABLE (ALL PAIRS) # ------------------------------------------------ print("2. Generating Near Table...") near_table = "in_memory\\near_table" arcpy.analysis.GenerateNearTable( in_features=output_fc, near_features=output_fc, out_table=near_table, search_radius=search_distance, location="NO_LOCATION", angle="NO_ANGLE", closest="ALL", closest_count="", method="GEODESIC" ) # ------------------------------------------------ # 3. BUILD ADJACENCY GRAPH # ------------------------------------------------ print("3. Building adjacency graph...") adj = defaultdict(set) with arcpy.da.SearchCursor(near_table, ["IN_FID", "NEAR_FID"]) as cur: for a, b in cur: if a != b: adj[a].add(b) adj[b].add(a) # ------------------------------------------------ # 4. FIND CONNECTED COMPONENTS # ------------------------------------------------ print("4. Finding connected components...") visited = set() components = [] all_oids = [row[0] for row in arcpy.da.SearchCursor(output_fc, [oid_field])] for oid in all_oids: if oid in visited: continue queue = deque([oid]) visited.add(oid) component = {oid} while queue: current = queue.popleft() for neighbor in adj.get(current, []): if neighbor not in visited: visited.add(neighbor) component.add(neighbor) queue.append(neighbor) components.append(component) print(f" → Found {len(components)} connected components") # ------------------------------------------------ # 5. SPLIT EACH COMPONENT INTO GROUPS OF 15 # ------------------------------------------------ print("5. Splitting components into groups of 15...") groups = {} group_id = 1 for comp in components: remaining = set(comp) while len(remaining) >= target_size: seed = remaining.pop() subgroup = [seed] queue = deque([seed]) while queue and len(subgroup) < target_size: current = queue.popleft() for n in adj.get(current, []): if n in remaining: remaining.remove(n) subgroup.append(n) queue.append(n) if len(subgroup) == target_size: break # Assign subgroup for oid in subgroup: groups[oid] = group_id group_id += 1 print(f" → Created {group_id - 1} groups of {target_size}") # ------------------------------------------------ # 6. WRITE GROUP FIELD # ------------------------------------------------ print("6. Writing GROUP_ID field...") field_names = [f.name for f in arcpy.ListFields(output_fc)] if group_field not in field_names: arcpy.management.AddField(output_fc, group_field, "LONG") with arcpy.da.UpdateCursor(output_fc, [oid_field, group_field]) as cur: for oid, _ in cur: cur.updateRow((oid, groups.get(oid))) # ------------------------------------------------ # 7. OUTPUT ONLY VALID GROUPS (OPTIONAL) # ------------------------------------------------ print("7. Creating output feature class (groups of 15 only)...") valid_fc = output_fc + "_Only15" arcpy.management.MakeFeatureLayer( output_fc, "groups15_lyr", f"{group_field} IS NOT NULL" ) arcpy.management.CopyFeatures("groups15_lyr", valid_fc) arcpy.management.Delete("groups15_lyr") # ------------------------------------------------ # CLEANUP # ------------------------------------------------ arcpy.management.Delete("in_memory") print("\nDONE") ... View more

Ya, I wasn't sure if there was or not. It does seem like this should be a capability somehow. ... View more

If they have to be within a certain distance of one another, how do you intend to filter out points that don't meet the min/and/or/max distance criteria? - Filter them out yes or some how mark them, for example a numeric field with 0 or 1 If points are too close to one another, do you intend to drop one or more? or take an average location? (eg aggregate) - Yes If points are too far from on another, do you intend to assign them to the nearest cluster perhaps breaking you number limit. - No, I am looking for exactly cluster of 15.   I have attached an image   I have the following code, which does mark the output feature class with 1, but it doesn't appear correct. The highlighted is what has 1 but you can see points that look like they should be part of the cluster import arcpy import pandas as pd from collections import Counter # --- INPUTS --- input_points = "Centroids_2_10" out_feature_class = r"C:\Temp\ClusterProject.gdb\PointClusters" clusters_exact15 = r"C:\Temp\\ClusterProject.gdb\Clusters_Exactly15" # --- CONSTANTS --- search_radius = 800 # feet target_cluster_size = 15 # Set overwrite output arcpy.env.overwriteOutput = True def find_clusters_of_exact_size(): """ Find clusters of exactly 15 points within 800 feet using graph-based method. This is the most reliable method for your specific requirement. """ print("="*70) print("FINDING CLUSTERS OF EXACTLY 15 POINTS WITHIN 800 FEET") print("="*70) # Step 1: Create near table print("\n1. Creating near table...") near_table = r"memory\near_table" arcpy.analysis.GenerateNearTable( in_features=input_points, near_features=input_points, out_table=near_table, search_radius=f"{search_radius} Feet", # Specify units method="PLANAR", closest="ALL", closest_count=0 # All points within radius ) print("2. Processing connections...") # Read near table near_data = [] with arcpy.da.SearchCursor(near_table, ['IN_FID', 'NEAR_FID', 'NEAR_DIST']) as cursor: for row in cursor: near_data.append(row) # Build adjacency dictionary adjacency = {} for in_fid, near_fid, dist in near_data: if in_fid != near_fid: # Skip self-connections adjacency.setdefault(in_fid, set()).add(near_fid) adjacency.setdefault(near_fid, set()).add(in_fid) print("3. Finding connected components (clusters)...") # Find connected components using BFS visited = set() all_clusters = [] for point_id in adjacency.keys(): if point_id not in visited: # Start new cluster cluster = set() stack = [point_id] while stack: current = stack.pop() if current not in visited: visited.add(current) cluster.add(current) # Add all neighbors for neighbor in adjacency.get(current, set()): if neighbor not in visited: stack.append(neighbor) all_clusters.append(cluster) # Find isolated points (not in adjacency at all) all_point_ids = set() with arcpy.da.SearchCursor(input_points, ["OID@"]) as cursor: for row in cursor: all_point_ids.add(row[0]) isolated_points = all_point_ids - visited for point_id in isolated_points: all_clusters.append({point_id}) print(f"\n4. Found {len(all_clusters)} total clusters") # Count cluster sizes cluster_sizes = [len(c) for c in all_clusters] size_counter = Counter(cluster_sizes) print("\nCluster size distribution:") for size, count in sorted(size_counter.items()): print(f" Size {size}: {count} clusters") # Find clusters with exactly 15 points exact_15_clusters = [c for c in all_clusters if len(c) == target_cluster_size] print(f"\n5. Found {len(exact_15_clusters)} clusters with exactly {target_cluster_size} points!") if not exact_15_clusters: print("No clusters found with exactly 15 points.") return # Step 5: Assign cluster IDs print("\n6. Assigning cluster IDs...") # Create mapping from point OID to cluster ID point_to_cluster = {} point_to_exact_cluster = {} # Assign IDs to all clusters for i, cluster in enumerate(all_clusters): cluster_id = i + 1 for point_id in cluster: point_to_cluster[point_id] = cluster_id # Assign special IDs to exact 15 clusters for i, cluster in enumerate(exact_15_clusters): exact_cluster_id = i + 1 for point_id in cluster: point_to_exact_cluster[point_id] = exact_cluster_id # Step 6: Create output feature classes print("\n7. Creating output feature classes...") # Create main output with all clusters arcpy.CopyFeatures_management(input_points, out_feature_class) # Add cluster fields arcpy.AddField_management(out_feature_class, "Cluster_ID", "LONG") arcpy.AddField_management(out_feature_class, "Exact15_ID", "LONG") arcpy.AddField_management(out_feature_class, "Is_Exact15", "SHORT") # Update with cluster assignments with arcpy.da.UpdateCursor(out_feature_class, ["OID@", "Cluster_ID", "Exact15_ID", "Is_Exact15"]) as cursor: for row in cursor: point_id = row[0] row[1] = point_to_cluster.get(point_id, -1) # -1 for isolated points row[2] = point_to_exact_cluster.get(point_id, -1) # -1 if not in exact 15 cluster row[3] = 1 if point_id in point_to_exact_cluster else 0 cursor.updateRow(row) print(f" Created: {out_feature_class}") # Create feature class with only exact 15 clusters where_clause = "Is_Exact15 = 1" temp_layer = "exact_clusters_layer" arcpy.MakeFeatureLayer_management(out_feature_class, temp_layer, where_clause) arcpy.CopyFeatures_management(temp_layer, clusters_exact15) arcpy.Delete_management(temp_layer) print(f" Created: {clusters_exact15}") # Step 7: Summary statistics print("\n" + "="*70) print("SUMMARY RESULTS") print("="*70) total_points = len(all_point_ids) isolated_count = len(isolated_points) clustered_points = total_points - isolated_count print(f"\nTotal points analyzed: {total_points}") print(f"Points in clusters: {clustered_points}") print(f"Isolated points (no neighbors within {search_radius} ft): {isolated_count}") print(f"Total clusters found: {len(all_clusters)}") print(f"Clusters with exactly {target_cluster_size} points: {len(exact_15_clusters)}") # Detailed breakdown of exact 15 clusters print(f"\nExact 15-point clusters (Cluster IDs):") for i, cluster in enumerate(exact_15_clusters, 1): print(f" Cluster {i}: {len(cluster)} points") # Suggestions for next steps print("\n" + "="*70) print("NEXT STEPS") print("="*70) print("\n1. View 'Clusters_Exactly15' to see the 15-point clusters") print("2. View 'PointClusters' to see all clusters (field: Cluster_ID)") print("3. Isolated points have Cluster_ID = -1") print("4. To adjust results, try:") print(f" - Increase search radius > {search_radius} ft for larger clusters") print(f" - Decrease search radius < {search_radius} ft for smaller clusters") print(f" - Target different cluster size (not 15)") # Clean up arcpy.Delete_management(near_table) return len(exact_15_clusters) def try_alternative_methods(): """Try alternative methods if graph-based doesn't work""" print("\n" + "="*70) print("ALTERNATIVE APPROACHES") print("="*70) print("\nIf the graph-based method doesn't work, try these:") print("\n1. SIMPLE BUFFER AND INTERSECT:") print(" - Buffer each point by 400 ft (half of 800)") print(" - Find overlapping buffers") print(" - Group points with overlapping buffers") print("\n2. SPATIAL JOIN:") print(" - Spatial join points to themselves within 800 ft") print(" - Create adjacency matrix from join results") print("\n3. MANUAL ITERATION:") print(" - For each point, find all neighbors within 800 ft") print(" - Build groups iteratively") print("\n4. USE EXISTING RESULTS:") print(" You already have 12 clusters of exactly 15 points!") print(" These are in 'Clusters_Exactly15' feature class") # Run the analysis if __name__ == "__main__": try: result = find_clusters_of_exact_size() if result == 0: print("\nNo exact 15-point clusters found. Trying alternative approaches...") try_alternative_methods() else: print(f"\nSUCCESS! Found {result} clusters of exactly 15 points.") print("\nThe graph-based method is complete and reliable.") except Exception as e: print(f"\nERROR: {str(e)}") print("\nTrying alternative approach...") try_alternative_methods()    ... View more

I forgot to mention that they have to be within a certain distance of each other. ... View more

I will have to try unchecking the feature. Thanks! ... View more

I have this issue as well, never been able to find answer, I have just used addDataFromPath() to re-add the layer. ... View more

account was a typeo, it's should be ParcelID. ... View more

Using what you posted, I didn't get an error. I have both the polygon and point layer added to my map, when start editing and create the point, the point is created but the ParcelID field is not populated in the point layer. The attribute rule is setup with the following settings, immediate calculation, Field - ParcelID, triggers - both insert & Update. ... View more

Thanks for the HTML tips, I made the adjustments you suggested, but nothing gets populated in the columns except for the Fire District column. I am still in Classic Map Viewer because we are still on an older version that doesn't have the "Source" button in the New Map Viewer.   Current HTML <table style="border-collapse: collapse; width: 100%;"> <tbody> <tr> <th style="border: 1px solid #ddd; padding: 8px;">Fire District</th> <th style="border: 1px solid #ddd; padding: 8px;">Contact Name</th> <th style="border: 1px solid #ddd; padding: 8px;">Phone</th> <th style="border: 1px solid #ddd; padding: 8px;">Email</th> <th style="border: 1px solid #ddd; padding: 8px;">Address</th> </tr> <tr> <td style="border: 1px solid #ddd; padding: 8px;">{FireDist}</td> <td style="border: 1px solid #ddd; padding: 8px;">{ContactName1}</td> <td style="border: 1px solid #ddd; padding: 8px;">{Phone1}</td> <td style="border: 1px solid #ddd; padding: 8px;">{Email1}</td> <td style="border: 1px solid #ddd; padding: 8px;">{Address1}</td> </tr> <tr> <td style="border: 1px solid #ddd; padding: 8px;"></td> <td style="border: 1px solid #ddd; padding: 8px;">{ContactName2}</td> <td style="border: 1px solid #ddd; padding: 8px;"></td> <td style="border: 1px solid #ddd; padding: 8px;">{Email2}</td> <td style="border: 1px solid #ddd; padding: 8px;">{Address2}</td> </tr> </tbody> </table>   Current Arcade code, // Debugging: Trim whitespace and standardize case for comparison var fireDistValue = Upper(Trim($feature.FireDist)); // Define contact information for Fire var Contact1 = { name: "xyz", phone: "xyzx", email: "xyz", address: "xyzxy" }; var Contact2 = { name: "xyz", phone: "xyzx", email: "xyzx", address: "xyzxy" }; // Function to format contact details function formatContact(contact) { return "Contact: " + contact.name + "\nPhone: " + contact.phone + "\nEmail: " + contact.email + "\nAddress: " + contact.address; } // Initialize result string var result = "Fire District: " + fireDistValue + "\n\n"; // Check the fire district and append the appropriate contact info if (fireDistValue == "xyzx xyzx FIRE") { result += formatContact(Contact1) + "\n\n" + formatContact(Contact2); } else { result += "No contact information available for this Fire District."; } // Return the result return result; ... View more

I clicked on the "Source" button and pasted the HTML code block, Click the source button once again. My expression name is Custom{expression/expr13}     table style="border-collapse: collapse; width: 100%;"> Fire District Contact Name Phone Email Address {FireDist} {ContactName1} {Phone1} {Email1} {Address1} {ContactName2} {Email2} {Address2}     I get , table style="border-collapse: collapse; width: 100%;"> Fire District Contact Name Phone Email Address xyzx xyzx FIRE   I can do <td rowspan="2" style="border: 1px solid #ddd; padding: 8px;">{expression/expr8}</td> <td style="border: 1px solid #ddd; padding: 8px;">{expression/expr13}</td>   But it puts all the info in both columns, it's not being separated out. ... View more