how to make a tool box from python script

MahdiMohammed — Mon, 22 Apr 2024 21:59:42 GMT

I have this code:

import arcpy

# Get the active map in the current project

active_map = arcpy.mp.ArcGISProject("CURRENT").activeMap

# Get the currently selected layer in tcrhe Contents pane

current_layer = active_map.listLayers()[0] # Assumes the first layer; adjust as needed

# Define the desired new name

new_layer_name = "HZ004"

# Check if the layer was found

if current_layer:

# Update the layer's name in the Contents pane

current_layer.name = new_layer_name

# Select and unselect the layer to trigger a refresh in the TOC

current_layer.visible = False

current_layer.visible = True

print(f"The layer name in the Contents pane has been successfully changed to {new_layer_name}.")

else:

print(f"No layer found in the project.")

# begin of required paramters

# name of the Lot layer that will be sampled

Polygon_Layer = "HZ004"

# sampling percentage

Sampling_Percent = 0.04

# the cell width and height for fishnet

# determine the size of each cell

cellSizeWidth = '25'

cellSizeHeight = '25'

#task Id and Lot Parameter

taskId = "NPA_CL_SY01_001"

lotNum = "LOT01"

Fishnet_Layer = "Fishnet004"

labels = 'NO_LABELS'

# end of required paramters

#begin of optional parameter - no need to change

# Number of rows and columns for fishnet if we don't determin width and height of cell

numRows = '0'

numColumns = '0'

# Create a point label feature class for the fishnet

#labels = 'NO_LABELS'

# name of fishnet layer

# Fishnet_Layer = "Fishnet004"

# end of optional parameter

import arcpy

# Get the current project

project = arcpy.mp.ArcGISProject("CURRENT")

# Get a list of all the layers in the project

layers = []

for map in project.listMaps():

layers.extend(map.listLayers())

# Find the feature class you want to use as the clip feature

clip_feature = None

for layer in layers:

if layer.isFeatureLayer:

if layer.name == Polygon_Layer :

clip_feature = layer

break

extent = project.activeView.camera.getExtent()

# fishnet options and creation

# Set the latitude and longitude coordinates for fishnet

latitude = extent.XMin

longitude = extent.YMin

maxLongitude = extent.YMax

maxLatitude = extent.XMax

originCoordinate = "%f %f" % (latitude,longitude)

yAxisCoordinate = "%f %f" % (latitude,maxLongitude)

conrnerCoordinate = "%f %f" % (maxLatitude,maxLongitude)

arcpy.CreateFishnet_management(Fishnet_Layer, originCoordinate,yAxisCoordinate , cellSizeWidth, cellSizeHeight, numRows, numColumns, conrnerCoordinate, labels, '', 'POLYGON')

arcpy.management.DefineProjection(Fishnet_Layer, arcpy.SpatialReference(arcpy.mp.ArcGISProject("CURRENT").activeMap.spatialReference.name))

# end of fishnet creation

arcpy.management.SelectLayerByLocation(Fishnet_Layer, "INTERSECT", Polygon_Layer, None, "NEW_SELECTION", "NOT_INVERT")

arcpy.management.AddField(Fishnet_Layer, "Box_Number", "TEXT", None, None, 15, "Box_Number", "NULLABLE", "NON_REQUIRED", '')

arcpy.management.CalculateField(Fishnet_Layer, "Box_Number", "SequentialNumber()", "PYTHON3", """# Calculates a sequential number

# More calculator examples at esriurl.com/CalculatorExamples

rec=0

def SequentialNumber():

global rec

pStart = 1

pInterval = 1

if (rec == 0):

rec = pStart

else:

rec = rec + pInterval

return rec""", "TEXT", "NO_ENFORCE_DOMAINS")

# Get the full path of the clip feature

clip_feature_path = clip_feature.dataSource

# Set the output feature class

output_fc = "Sampling_Polygon"

# Perform the clip

arcpy.analysis.Clip(Fishnet_Layer, clip_feature_path, output_fc)

import arcpy

import random

import string

# Set the workspace to the current workspace

arcpy.env.workspace = arcpy.mp.ArcGISProject("CURRENT").defaultGeodatabase

# Generate a random domain name starting with "Random_Status"

domain_name = "Sampling_Status_" + ''.join(random.choices(string.ascii_uppercase + string.digits, k=3))

# Check if the domain already exists and delete it if it does

if arcpy.Exists(domain_name):

arcpy.management.DeleteDomain(arcpy.env.workspace, domain_name)

# Create the domain

arcpy.management.CreateDomain(arcpy.env.workspace, domain_name, "Sampling Domain Status", "TEXT", "CODED")

arcpy.management.AddCodedValueToDomain(arcpy.env.workspace, domain_name, "Random Pending", "Random Pending")

arcpy.management.AddCodedValueToDomain(arcpy.env.workspace, domain_name, "Random Passed", "Random Passed")

arcpy.management.AddCodedValueToDomain(arcpy.env.workspace, domain_name, "Random Failed", "Random Failed")

arcpy.management.AddCodedValueToDomain(arcpy.env.workspace, domain_name, "Random Replaced", "Random Replaced")

arcpy.management.AddCodedValueToDomain(arcpy.env.workspace, domain_name, "Targeted Pending", "Targeted Pending")

arcpy.management.AddCodedValueToDomain(arcpy.env.workspace, domain_name, "Targeted Passed", "Targeted Passed")

arcpy.management.AddCodedValueToDomain(arcpy.env.workspace, domain_name, "Targeted Failed", "Targeted Failed")

# Get the "Sampling_Polygon" feature layer

sampling_polygon_layer = arcpy.mp.ArcGISProject("CURRENT").listMaps()[0].listLayers("Sampling_Polygon")[0]

# Check if the "Random_Status" field already exists and delete it if it does

if "Sampling_Status" in [f.name for f in arcpy.ListFields(sampling_polygon_layer)]:

arcpy.management.DeleteField(sampling_polygon_layer, "Sampling_Status")

# Add the "Random_Status" field to the layer

arcpy.management.AddField(sampling_polygon_layer, "Sampling_Status", "TEXT", field_length=50)

# Assign the random domain to the "Random_Status" field in the layer

arcpy.management.AssignDomainToField(sampling_polygon_layer, "Sampling_Status", domain_name)

arcpy.management.AddField("Sampling_Polygon", "X", "TEXT", None, None, None, "X", "NULLABLE", "NON_REQUIRED", '')

arcpy.management.AddField("Sampling_Polygon", "Y", "TEXT", None, None, None, "Y", "NULLABLE", "NON_REQUIRED", '')

arcpy.management.CalculateGeometryAttributes("Sampling_Polygon", "X INSIDE_X;Y INSIDE_Y", '', '', None, "SAME_AS_INPUT")

arcpy.management.AddField("Sampling_Polygon", "TaskID", "TEXT", None, None, None, "TaskID", "NULLABLE", "NON_REQUIRED", '')

# Use an update cursor to update the field

with arcpy.da.UpdateCursor("Sampling_Polygon", "TaskID") as cursor:

for row in cursor:

row[0] = taskId

cursor.updateRow(row)

# Clean up the cursor

del cursor

arcpy.management.AddField("Sampling_Polygon", "Lot_Number", "TEXT", None, None, None, "Lot_Number", "NULLABLE", "NON_REQUIRED", '')

# Use an update cursor to update the field

with arcpy.da.UpdateCursor("Sampling_Polygon", "Lot_Number") as cursor:

for row in cursor:

row[0] = lotNum

cursor.updateRow(row)

# Clean up the cursor

del cursor

import arcpy

import random

# Get the "Sampling_Polygon" feature layer

sampling_polygon_layer = arcpy.mp.ArcGISProject("CURRENT").listMaps()[0].listLayers("Sampling_Polygon")[0]

# Define the feature layer and field names

feature_layer = "Sampling_Polygon"

area_field = "SHAPE@AREA"

oid_field = arcpy.Describe(feature_layer).OIDFieldName

# Calculate the area threshold for selecting records

total_area = sum(row[0] for row in arcpy.da.SearchCursor(feature_layer, area_field))

area_threshold = total_area * Sampling_Percent # Sampling_Percent% of the total area

# Select a random box to start with

selected_boxes = []

selected_area = 0

with arcpy.da.SearchCursor(feature_layer, [area_field, oid_field, "SHAPE@"]) as cursor:

row = random.choice([row for row in cursor])

selected_boxes.append((row[1], row[2]))

selected_area += row[0]

# Determine the minimum distance between selected boxes

min_distance = max([box[1].length / 100 for box in selected_boxes])

# Exclude boxes that are within the minimum distance from the selected boxes

excluded_boxes = set()

while selected_area < area_threshold:

# Get the distance between the selected boxes and all other boxes

distances = {}

with arcpy.da.SearchCursor(feature_layer, [oid_field, "SHAPE@"]) as cursor:

for row in cursor:

if row[0] not in excluded_boxes:

distance = min([box[1].distanceTo(row[1]) for box in selected_boxes])

distances[row[0]] = distance

# Select the box with the maximum distance from the selected boxes

if not distances:

break

max_distance_oid = max(distances, key=distances.get)

if distances[max_distance_oid] >= min_distance:

excluded_boxes.add(max_distance_oid)

with arcpy.da.SearchCursor(feature_layer, [area_field, oid_field, "SHAPE@"], "{0} = {1}".format(oid_field, max_distance_oid)) as cursor:

row = next(cursor)

selected_boxes.append((row[1], row[2]))

selected_area += row[0]

# Update the minimum distance between selected boxes

min_distance = max([box[1].length / 100 for box in selected_boxes])

else:

# If no box can be selected, reduce the min_distance and try again

min_distance *= 0.9

# Select the randomly selected boxes by their ObjectIDs

selected_oids = [box[0] for box in selected_boxes]

where_clause = "{0} IN ({1})".format(oid_field, ','.join(str(oid) for oid in selected_oids))

arcpy.SelectLayerByAttribute_management(feature_layer, "NEW_SELECTION", where_clause)

arcpy.management.CalculateField("Sampling_Polygon", "Sampling_Status", "'Random Pending'", "PYTHON3", "", "TEXT", "NO_ENFORCE_DOMAINS")

# Set the workspace to the current workspace

arcpy.env.workspace = arcpy.mp.ArcGISProject("CURRENT").defaultGeodatabase

# Delete the "Fishnet_Sampling_Polygon" layer if it exists

if arcpy.Exists(Fishnet_Layer):

arcpy.management.Delete(Fishnet_Layer)

# Delete the "Exported_Polygon" layer if it exists

if arcpy.Exists(Polygon_Layer 😞

arcpy.management.Delete(Polygon_Layer )

import arcpy

# Get the active map in the current project

active_map = arcpy.mp.ArcGISProject("CURRENT").activeMap

# Get the currently selected layer in the Contents pane

current_layer = active_map.listLayers()[0] # Assumes the first layer; adjust as needed

# Check if a layer is selected

if current_layer:

# Get the layer name

polygon_layer_name = current_layer.name

# Set the output point feature class name

output_points_fc = f"{polygon_layer_name}_Vertices_Points"

# Use the Feature Vertices to Points tool

arcpy.management.FeatureVerticesToPoints(polygon_layer_name, output_points_fc, "ALL")

print(f"Points have been created from the vertices of the {polygon_layer_name} layer.")

else:

print("No layer found in the project.")

and its working when I run it in python window.
can any one help me to craete a tool box based on this code

topic Re: how to make a tool box from python script in Python Questions

how to make a tool box from python script

Re: how to make a tool box from python script

Re: how to make a tool box from python script