[WIP] Polgon Centerline

SWCD Tools.PolygonCenterline.pyt.xml

@@ -0,0 +1,2 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<metadata xml:lang="en"><Esri><CreaDate>20260622</CreaDate><CreaTime>17115200</CreaTime><ArcGISFormat>1.0</ArcGISFormat><SyncOnce>TRUE</SyncOnce></Esri></metadata>

SWCD Tools.pyt

@@ -1,6 +1,6 @@
 # -*- coding: utf-8 -*-
 
-from src import *
+from src import AgAssessment, BufferTools, FluvialGeomorphology, Hydrology, TerrainAnalysis, TerrainModification, TileDrainage, Utilities
 
 class Toolbox(object):
     def __init__(self):
@@ -10,47 +10,48 @@ class Toolbox(object):
         self.alias = "SWCD Tools"
 
         tools = [
-            SlopePolygon,
-            ContourPolygon,
-            BurnCulverts,
-            CollectRasters,
-            ExportLayouts,
-            StreamNetwork,
-            StreamElevation,
-            GenerateCrossSections,
-            LocalMinimums,
-            RelativeElevationModel,
-            TopographicPositionIndex,
-            PointPlots,
-            ShrubClusters,
-            BufferPotential,
-            StreamPowerIndex,
-            LeastAction,
-            RunoffCurveNumber,
-            CalculateEFH2,
-            WatershedSize,
-            StreambankDetection,
-            LandscapePosition,
-            VBET,
-            WatershedDelineation,
-            SubBasinDelineation,
-            RemoveUnused,
-            DecisionTree,
-            ImageDifferencingSetup,
-            ImageDifferencing,
-            ImageDifferencingClouds,
-            TopographicWetness,
-            PotentialWetlands,
-            BermAnalysis,
-            DamRemoval,
-            DefineParcels,
-            Agland,
-            NonAg,
-            GeocodeAddress,
-            Forest,
-            Process,
-            Export,
-            Restart
+            AgAssessment.DefineParcels,
+            AgAssessment.Agland,
+            AgAssessment.NonAg,
+            AgAssessment.Forest,
+            AgAssessment.Process,
+            AgAssessment.Export,
+            AgAssessment.Restart,
+            BufferTools.PointPlots,
+            BufferTools.ShrubClusters,
+            BufferTools.BufferPotential,
+            FluvialGeomorphology.StreamNetwork,
+            FluvialGeomorphology.LeastAction,
+            FluvialGeomorphology.StreambankDetection,
+            FluvialGeomorphology.StreamElevation,
+            FluvialGeomorphology.GenerateCrossSections,
+            FluvialGeomorphology.PolygonCenterline,
+            Hydrology.CalculateEFH2,
+            Hydrology.RunoffCurveNumber,
+            Hydrology.SubBasinDelineation,
+            Hydrology.WatershedDelineation,
+            Hydrology.WatershedSize,
+            TerrainModification.BermAnalysis,
+            TerrainModification.DamRemoval,
+            TerrainModification.BurnCulverts,
+            TerrainAnalysis.StreamPowerIndex,
+            TerrainAnalysis.LandscapePosition,
+            TerrainAnalysis.VBET,
+            TerrainAnalysis.TopographicWetness,
+            TerrainAnalysis.RelativeElevationModel,
+            TerrainAnalysis.PotentialWetlands,
+            TerrainAnalysis.TopographicPositionIndex,
+            TileDrainage.DecisionTree,
+            TileDrainage.ImageDifferencingSetup,
+            TileDrainage.ImageDifferencing,
+            TileDrainage.ImageDifferencingClouds,
+            Utilities.SlopePolygon,
+            Utilities.ContourPolygon,
+            Utilities.ExportLayouts,
+            Utilities.LocalMinimums,
+            Utilities.GeocodeAddress,
+            Utilities.RemoveUnused,
+            Utilities.CollectRasters,
         ]
 
         # List of tool classes associated with this toolbox

SWCD Tools.pyt.xml

@@ -1,2 +1,2 @@
 <?xml version="1.0"?>
-<metadata xml:lang="en"><Esri><CreaDate>20251119</CreaDate><CreaTime>09340100</CreaTime><ArcGISFormat>1.0</ArcGISFormat><SyncOnce>TRUE</SyncOnce><ModDate>20260625</ModDate><ModTime>181028</ModTime><scaleRange><minScale>150000000</minScale><maxScale>5000</maxScale></scaleRange><ArcGISProfile>ItemDescription</ArcGISProfile></Esri><toolbox name="SWCD Tools" alias="SWCD Tools"><arcToolboxHelpPath>c:\program files\arcgis\pro\Resources\Help\gp</arcToolboxHelpPath><toolsets/></toolbox><dataIdInfo><idCitation><resTitle>SWCD Tools</resTitle></idCitation><idPurp>A collection of ArcGIS tools related to conservation and natural resources</idPurp><searchKeys><keyword>Soil and Water Toolbox GIS ArcGIS</keyword></searchKeys></dataIdInfo><distInfo><distributor><distorFormat><formatName>ArcToolbox Toolbox</formatName></distorFormat></distributor></distInfo><mdHrLv><ScopeCd value="005"></ScopeCd></mdHrLv><mdDateSt Sync="TRUE">20251221</mdDateSt></metadata>
+<metadata xml:lang="en"><Esri><CreaDate>20251119</CreaDate><CreaTime>09340100</CreaTime><ArcGISFormat>1.0</ArcGISFormat><SyncOnce>TRUE</SyncOnce><ModDate>20260716</ModDate><ModTime>133353</ModTime><scaleRange><minScale>150000000</minScale><maxScale>5000</maxScale></scaleRange><ArcGISProfile>ItemDescription</ArcGISProfile></Esri><toolbox name="SWCD Tools" alias="SWCD Tools"><arcToolboxHelpPath>c:\program files\arcgis\pro\Resources\Help\gp</arcToolboxHelpPath><toolsets/></toolbox><dataIdInfo><idCitation><resTitle>SWCD Tools</resTitle></idCitation><idPurp>A collection of ArcGIS tools related to conservation and natural resources</idPurp><searchKeys><keyword>Soil and Water Toolbox GIS ArcGIS</keyword></searchKeys></dataIdInfo><distInfo><distributor><distorFormat><formatName>ArcToolbox Toolbox</formatName></distorFormat></distributor></distInfo><mdHrLv><ScopeCd value="005"></ScopeCd></mdHrLv><mdDateSt Sync="TRUE">20251221</mdDateSt></metadata>

src/AgAssessment/Agland.py

@@ -11,7 +11,7 @@
 import arcpy
 
 from .DefineParcels import AG_ASSESSMENT_GDB_NAME
-from ..helpers import license, reload_module, log, error
+from ..helpers import license, reload_module, log, warn, error
 from ..helpers import setup_environment as setup
 from ..helpers import validate_spatial_reference as validate
 
@@ -64,16 +64,16 @@ def execute(self, parameters, messages):
             m = None
             try:
                 m = project.listMaps(parcel)[0]
-            except:
-                log("unable to find map for {}, results may be incomplete".format(parcel))
+            except Exception:
+                warn("unable to find map for {}, results may be incomplete".format(parcel))
                 continue
 
             # get parcel layer or drop off of map
             parcel_lyr = None
             try:
                 parcel_lyr = m.listLayers("*_{}".format(parcel))[0]
-            except:
-                log("no appropriate parcel layer found for {}, results may be incomplete".format(parcel))
+            except Exception:
+                warn("no appropriate parcel layer found for {}, results may be incomplete".format(parcel))
                 continue
 
             # check how many pieces are selected

src/AgAssessment/DefineParcels.py

@@ -309,7 +309,7 @@ def execute(self, parameters, messages):
             # turn off parcel layer
             try:
                 parcel_layer.visible = False
-            except:
+            except Exception:
                 # parcel_layer is a shapefile
                 pass
 

src/AgAssessment/Forest.py

@@ -11,7 +11,7 @@
 import arcpy
 
 from .DefineParcels import AG_ASSESSMENT_GDB_NAME
-from ..helpers import license, reload_module, log, error
+from ..helpers import license, reload_module, log, warn, error
 from ..helpers import setup_environment as setup
 from ..helpers import validate_spatial_reference as validate
 
@@ -64,16 +64,16 @@ def execute(self, parameters, messages):
             m = None
             try:
                 m = project.listMaps(parcel)[0]
-            except:
-                log("unable to find map for {}, results may be incomplete".format(parcel))
+            except Exception:
+                warn("unable to find map for {}, results may be incomplete".format(parcel))
                 continue
 
             # get parcel layer or drop off of map
             parcel_lyr = None
             try:
                 parcel_lyr = m.listLayers("*_{}".format(parcel))[0]
-            except:
-                log("no appropriate parcel layer found for {}, results may be incomplete".format(parcel))
+            except Exception:
+                warn("no appropriate parcel layer found for {}, results may be incomplete".format(parcel))
                 continue
 
             # check how many pieces are selected

src/AgAssessment/NonAg.py

@@ -11,7 +11,7 @@
 import arcpy
 
 from .DefineParcels import AG_ASSESSMENT_GDB_NAME
-from ..helpers import license, reload_module, log, error
+from ..helpers import license, reload_module, log, warn, error
 from ..helpers import setup_environment as setup
 from ..helpers import validate_spatial_reference as validate
 
@@ -64,16 +64,16 @@ def execute(self, parameters, messages):
             m = None
             try:
                 m = project.listMaps(parcel)[0]
-            except:
-                log("unable to find map for {}, results may be incomplete".format(parcel))
+            except Exception:
+                warn("unable to find map for {}, results may be incomplete".format(parcel))
                 continue
 
             # get parcel layer or drop off of map
             parcel_lyr = None
             try:
                 parcel_lyr = m.listLayers("*_{}".format(parcel))[0]
-            except:
-                log("no appropriate parcel layer found for {}, results may be incomplete".format(parcel))
+            except Exception:
+                warn("no appropriate parcel layer found for {}, results may be incomplete".format(parcel))
                 continue
 
             # check how many pieces are selected

src/AgAssessment/Process.py

@@ -139,7 +139,7 @@ def execute(self, parameters, messages):
             m = None
             try:
                 m = project.listMaps(parcel)[0]
-            except:
+            except Exception:
                 warn("unable to find map for {}, results may be incomplete".format(parcel))
                 continue
 
@@ -151,7 +151,7 @@ def execute(self, parameters, messages):
             try:
                 lyt = project.listLayouts(parcel)[0]
                 layouts.append(lyt)
-            except:
+            except Exception:
                 warn("couldn't find layout for parcel {}, results may be incomplete".format(parcel))
                 continue
 

src/AgAssessment/Restart.py

@@ -12,7 +12,7 @@
 import arcpy
 
 from .DefineParcels import AG_ASSESSMENT_GDB_NAME
-from ..helpers import license, reload_module, log, empty_workspace
+from ..helpers import license, reload_module, log, warn, empty_workspace
 from ..helpers import setup_environment as setup
 
 class Restart(object):
@@ -81,8 +81,8 @@ def execute(self, parameters, messages):
                 with open(cache_file_path) as file:
                     cache = json.load(file)
                     parcels = cache["parcels"]
-            except:
-                log("Unable to find cache file and complete transaction. Please manually refresh the tool or manually clear out data.")
+            except Exception:
+                warn("Unable to find cache file and complete transaction. Please manually refresh the tool or manually clear out data.")
                 return
 
             # clear out maps, layouts, and feature classes
@@ -92,8 +92,8 @@ def execute(self, parameters, messages):
                 lyt = None
                 try:
                     lyt = project.listLayouts(parcel)[0]
-                except:
-                    log("couldn't find layout for parcel {}, results may be incomplete".format(parcel))
+                except Exception:
+                    warn("couldn't find layout for parcel {}, results may be incomplete".format(parcel))
                     continue
 
                 # delete map
@@ -103,8 +103,8 @@ def execute(self, parameters, messages):
                 m = None
                 try:
                     m = project.listMaps(parcel)[0]
-                except:
-                    log("unable to find map for parcel {}, results may be incomplete".format(parcel))
+                except Exception:
+                    warn("unable to find map for parcel {}, results may be incomplete".format(parcel))
                     continue
 
                 # delete map

src/BufferTools/BufferPotential.py

@@ -7,7 +7,7 @@
 # --------------------------------------------------------------------------------
 import arcpy
 
-from ..helpers import license, empty_workspace, reload_module, log, raster_and_layer
+from ..helpers import license, empty_workspace, reload_module, log, warn, raster_and_layer, ArealUnit, LinearUnit
 from ..helpers import setup_environment as setup
 from ..helpers import validate_spatial_reference as validate
 
@@ -162,9 +162,8 @@ def execute(self, parameters, messages):
 
         log("reading in parameters")
         stream = parameters[0].value
-        min_width = parameters[1].valueAsText
-        min_acres, min_acres_unit = parameters[2].valueAsText.split(" ")
-        min_acres = float(min_acres) * arcpy.ArealUnitConversionFactor(min_acres_unit, "AcresUS")
+        min_width = LinearUnit(parameters[1].valueAsText)
+        min_area = ArealUnit(parameters[2].valueAsText)
         extent = parameters[3].value
         output_file = parameters[4].valueAsText
         land_use_raster, _ = raster_and_layer(parameters[5].value)
@@ -188,7 +187,7 @@ def execute(self, parameters, messages):
 
         # pairwise buffer stream
         log("creating buffer polygon around stream")
-        arcpy.analysis.PairwiseBuffer(stream, scratch_stream_buffer, min_width, "ALL", "", "GEODESIC", "")
+        arcpy.analysis.PairwiseBuffer(stream, scratch_stream_buffer, str(min_width), "ALL", "", "GEODESIC", "")
 
         # clip land uses to buffer
         log("extracting land use data inside buffer area")
@@ -211,7 +210,7 @@ def execute(self, parameters, messages):
                     land_use_sql_query += " Or {} = '{}'".format(land_use_field, value)
             scratch_land_use = arcpy.sa.ExtractByAttributes(land_use_raster_clip, land_use_sql_query)
         else:
-            log("no valid land uses found in area, please try again with land uses found in analysis area")
+            warn("no valid land uses found in area, please try again with land uses found in analysis area")
             return
 
         # convert land usage output to polygon
@@ -241,11 +240,11 @@ def execute(self, parameters, messages):
 
         # calculate acreage
         log("calculating acreage of planting areas")
-        arcpy.management.AddField(scratch_dissolve, "Acres", "FLOAT", 2, 2)
-        arcpy.management.CalculateGeometryAttributes(scratch_dissolve, geometry_property=[["Acres", "AREA_GEODESIC"]], area_unit="ACRES_US")
+        arcpy.management.AddField(scratch_dissolve, min_area.unit, "FLOAT", 2, 2)
+        arcpy.management.CalculateGeometryAttributes(scratch_dissolve, geometry_property=[[min_area.unit, "AREA_GEODESIC"]], area_unit=min_area.full_unit())
 
         # drop acreage < threshold
-        sql_query = "Acres >= {}".format(min_acres)
+        sql_query = "{} >= {}".format(min_area.unit, min_area.area)
         arcpy.analysis.Select(scratch_dissolve, output_file, sql_query)
 
         # add output to map

src/BufferTools/PointPlots.py

@@ -18,7 +18,7 @@
 import arcpy
 import platform
 
-from ..helpers import license, empty_workspace, set_required_parameter, reload_module, log
+from ..helpers import license, empty_workspace, set_required_parameter, reload_module, log, warn
 from ..helpers import setup_environment as setup
 from ..helpers import validate_spatial_reference as validate
 
@@ -141,8 +141,8 @@ def execute(self, parameters, messages):
                 log("create sampling locations")
                 arcpy.management.CreateSpatialSamplingLocations(scratch_buffer, output_points, sampling_method="STRAT_POLY", strata_id_field=None, strata_count_method="PROP_AREA", num_samples=num, geometry_type="POINT", min_distance="{} Feet".format(radius*2))
 
-            except:
-                log("Failed to create {} point plots with a radius of {} feet. It is likely because the buffer is too narrow to fit all of the point plots.".format(num, radius))
+            except Exception:
+                warn("Failed to create {} point plots with a radius of {} feet. It is likely because the buffer is too narrow to fit all of the point plots.".format(num, radius))
 
                 radius = 11.8
                 num = int(math.ceil(acreage * 10))