Latest Contributions by JohannesLindner

To format code: JohannesLindner_0-1670479234964.png JohannesLindner_1-1670479257402.png   I have added another layer area intersect, it appears as a field but it isn't calculating the intersect area That's because you put the intersection with the WHA features into a separate loop. They get calculated and then overwritten without being stored in the output dict.   I would also like to group the phases Fortunately, that is quite simple in this case. I used When() to get a default value if the phase is in neither of those groups (eg when it's null).     // blablabla // loop over the treatment area features for(var ta in fsTa) { // get the intersecting sugar area var intSugar = Intersects(fsSugar, ta) var sugarArea = 0 for(var s in intSugar) { sugarArea += AreaGeodetic(Intersection(s, ta), "hectares") } // get the intersecting WHA area var intWHA = Intersects(fsWHA, ta) var WHAArea = 0 for(var w in intWHA) { WHAArea += AreaGeodetic(Intersection(w, ta), "hectares") } // calculate the status value var status = When(Includes([1, 2, 3, 4], ta.Phase), "Status 1", ta.Phase == 5, "Status 2", "Default") // append to the output dict var f = {attributes: {TA_Status: status, TA_Ha: AreaGeodetic(ta, "hectares"), TA_Sugar: sugarArea, TA_WHA: WHAArea}} Push(combinedDict.features, f) } // convert to Featureset and return return Featureset(Text(combinedDict))     ... View more

AFAIK, there is no out of the box tool for that. This thread lists a few algorithms the get an incircle: algorithm - Largest circle inside a non-convex polygon - Stack Overflow   I took one of the approaches: tesselate the polygon with a Voronoi diagram (Thiessen polygons in ArcGIS) find the node with the greates distance to the polygon's border, that's the center of the incircle def create_incircle(polygon): add = arcpy.env.addOutputsToMap arcpy.env.addOutputsToMap = False # densify polygon to get rid of curves, which only have 2 vertices polygon = polygon.densify("DISTANCE", 10, 1) # extract vertices into a point fc vertices = arcpy.management.CreateFeatureclass("memory", "vertices", "POINT", spatial_reference=polygon.spatialReference) with arcpy.da.InsertCursor(vertices, ["SHAPE@"]) as cursor: for part in polygon: for v in part: cursor.insertRow([v]) # create Thiessen polygons thiessen = arcpy.analysis.CreateThiessenPolygons(vertices, "memory/thiessen") # intersect these with themselves to get candidate points test_points = arcpy.analysis.Intersect([thiessen], "memory/TestPoints", output_type="POINT") # get the point inside the polygon with the greatest distance to the polygon's boundary candidates = [] boundary = polygon.boundary() with arcpy.da.SearchCursor(test_points, ["SHAPE@"]) as cursor: for p, in cursor: if p.disjoint(polygon): continue candidates.append([p, p.distanceTo(boundary)]) center = sorted(candidates, key=lambda c: c[1])[-1] # buffer that point ith the distance to the boundary, return arcpy.env.addOutputsToMap = add return center[0].buffer(center[1]) # create the ouput fc incircles = arcpy.management.CreateFeatureclass("memory", "Incircles") arcpy.management.AddField(incircles, "FID", "LONG") # create the incircles with arcpy.da.InsertCursor(incircles, ["SHAPE@", "FID"]) as i_cursor: with arcpy.da.SearchCursor("TestPolygons", ["SHAPE@", "OID@"]) as s_cursor: for shp, oid in s_cursor: ic = create_incircle(shp) i_cursor.insertRow([ic, oid])   JohannesLindner_0-1670429585805.png   This is quite slow (1 - 2 seconds per feature), probably because the densification in line 5 adds lots of vertices. This can probably be optimized much, but it works as a quick-and-dirty approach. Sadly, it requires an Advanced license for the Thiessen polygons. ... View more

// Access data layers from portal var port = Portal( 'http://www.argis.com') var fsTa = FeatureSetByPortalItem(port, 'itemid', 0) var fsSugar = FeatureSetByPortalItem(port,'itemid', 0) //Create empty dictionary var sugarDict = { fields:[ {name: "TA_Name", type: "esriFieldTypeString"}, {name: "TA_Ha", type: "esriFieldTypeDouble"}, {name: "TA_Sugar", type: "esriFieldTypeDouble"}, ], geometryType: "", features: [], } // loop over the treatment area features for(var ta in fsTa) { // get the intersecting sugar areas var intSugar = Intersects(fsSugar, ta) // loop over those features and get the sum of the intersections var sugarArea = 0 for(var s in intSugar) { sugarArea += Area(Intersection(s, ta), "hectares") } // append to the output dict var f = {attributes: {TA_Name: ta.Name, TA_Ha: Area(ta, "hectares"), TA_Sugar: sugarArea}} Push(sugarDict.features, f) } // convert to Featureset and return return Featureset(Text(sugarDict)) ... View more

You can get the symbology of a layer. Form there, you'll have to follow the rabbit hole to the single symbology items. This script will loop through a list of layers in the active map, analyze their symbology and add a new field to each layer with the symbology's label value: layer_names = ["TestPolygons"] aprx = arcpy.mp.ArcGISProject("current") for name in layer_names: # get the layer layer = aprx.activeMap.listLayers(name)[0] # get the symbology fields and items renderer = layer.symbology.renderer fields = renderer.fields items = renderer.groups[0].items # create a dict {values: label} item_dict = {tuple(i.values[0]): i.label for i in items} # add a new field arcpy.management.AddField(layer, "Unit", "TEXT") # insert the label values with arcpy.da.UpdateCursor(layer, ["Unit"] + fields) as cursor: for row in cursor: key = tuple([str(x) for x in row[1:]]) # item.values stores values as string try: row[0] = item_dict[key] cursor.updateRow(row) except KeyError: # no symbology found pass   JohannesLindner_0-1670407418636.png   ... View more

Your first problem (Having to halve the area): No idea. I can't reproduce this, all areas are the same as when i measure them manually in Pro. The most obvious reason could be that you have doubled soil polygons. Also, compare the soil area sum against the parcel area. They should be the same.   Your second problem (combine areas of the same type): Use a dictionary with the soil type as keys. If you encounter a type the first time, store its area. If you encounter it again, add its area to the stored one.   //list soils in a parcel by type and acreage. //get the features from the soil layer is overlaped by the parcel layer var gensoils = Intersects(FeatureSetByName($map,"gensoils"),$feature) // create a dictionary that stores the areas for each soil type var soil_areas = Dictionary() // create a variable that stores the sum of all soil areas var soil_area_sum = 0 // loop over the soil polygons for(var s in gensoils) { // get the area of the intersection var a = Area(Intersection($feature, s), "acre") // add it to the sum soil_area_sum += a // store it in the dict var type = s.muname if(HasKey(soil_areas, type)) { soil_areas[type] += a } else { soil_areas[type] = a } } // get some general result info var result = [ `The total area of the parcel is ${Round(Area($feature, "acres"))} Acre(s).`, `The total area of soil in this parcel is ${Round(soil_area_sum)} Acre(s).`, `There are ${Count(gensoils)} type(s) of General Soils.` ] // loop over the soil types and append their area and percentage for(var type in soil_areas) { var a = soil_areas[type] var p = a / soil_area_sum * 100 Push(result, `${type} ---- ${Round(p, 2)}% or ${Round(a, 2)} ac(s)`) } return Concatenate(result, TextFormatting.NewLine) JohannesLindner_0-1670402793022.png JohannesLindner_1-1670402808026.png   ... View more

var distances = [] // get the polygon parts var rings = Geometry($feature).rings // loop over the parts for(var r in rings) { var ring = rings[r] // loop over the vertices (skip first) for(var v = 1; v < Count(ring); v++) { Push(distances, Distance(ring[v - 1], ring[v])) } } return Max(distances)   This will get the longest side of any polygon, no matter the part count or the form. You could simplify it by only extracting the first ring and checking only the first three vertices (giving you the two sides), but it would not make any noticeable time difference. ... View more

Starting from Dan's suggestion of convexHull(), which seemed good: with arcpy.da.UpdateCursor("TestPolygons", ["SHAPE@", "IntegerField"]) as cursor: for shp, n in cursor: n = len(shp[0]) - len(shp.convexHull()[0]) # only for singlepart polygons cursor.updateRow([shp, n])   This works well for simple geometries: JohannesLindner_0-1670325926043.png   But it starts to break down for more complex geometries, where convex vertices are not always on the convex hull: JohannesLindner_1-1670326259996.png   I was ready to start fussing with angles between polylines when I realized that you can quite simply say whether a vertex is concave or convex: cretae a triangle of the vertex and its neighbours. If that triangle is inside the polygon, the vertex is convex, else it's concave. JohannesLindner_2-1670329512390.png   def get_concave_vertex_count(polygon): n = 0 for part in polygon: for i, vertex in enumerate(part): prev_vertex = part[i-1] try: next_vertex = part[i+1] except IndexError: next_vertex = part[0] triangle = arcpy.Polygon(arcpy.Array([prev_vertex, vertex, next_vertex])) n += triangle.touches(polygon) # touch: only boundaries intersect return n with arcpy.da.UpdateCursor("PolygonLayer", ["SHAPE@", "ConcaveVertices"]) as cursor: for shp, n in cursor: n = get_concave_vertex_count(shp) cursor.updateRow([shp, n]) ... View more

It is absolutely possible. I changed your example table to reflect that the counts are not taken daily. I used the same values for depots A and B: JohannesLindner_4-1670317581961.png     For the production, I just assumed that depot A fills 10 sandbags per day, while depot B manages to fill 20. JohannesLindner_3-1670317519975.png       // load the count fs var count_fs = Featureset(Text({ fields: [ {name: "DateOfCount", type: "esriFieldTypeDate"}, {name: "Depot", type: "esriFieldTypeString"}, {name: "Count", type: "esriFieldTypeInteger"}, ], features: [ {attributes: {DateOfCount: Number(Date(2022,11,1)), Depot: "A", Count: 100}}, {attributes: {DateOfCount: Number(Date(2022,11,6)), Depot: "A", Count: 50}}, {attributes: {DateOfCount: Number(Date(2022,11,22)), Depot: "A", Count: 150}}, {attributes: {DateOfCount: Number(Date(2022,11,1)), Depot: "B", Count: 100}}, {attributes: {DateOfCount: Number(Date(2022,11,6)), Depot: "B", Count: 50}}, {attributes: {DateOfCount: Number(Date(2022,11,22)), Depot: "B", Count: 150}}, ], geometryType: "" })) //return count_fs // load the production fs var production_fs = { fields: [ {name: "DateOfProduction", type: "esriFieldTypeDate"}, {name: "Depot", type: "esriFieldTypeString"}, {name: "Count", type: "esriFieldTypeInteger"}, ], features: [], geometryType: "" } for(var i = 0; i < 30; i++) { Push(production_fs.features, {attributes: {DateOfProduction: Number(DateAdd(Date(2022,11,1), i, "days")), Depot: "A", Count: 10}}) Push(production_fs.features, {attributes: {DateOfProduction: Number(DateAdd(Date(2022,11,1), i, "days")), Depot: "B", Count: 20}}) } production_fs = Featureset(Text(production_fs)) //return production_fs // create the output dict var out = { fields: [ {name: "DateOfCount", type: "esriFieldTypeDate"}, {name: "Depot", type: "esriFieldTypeString"}, {name: "Count", type: "esriFieldTypeInteger"}, {name: "Taken", type: "esriFieldTypeInteger"}, {name: "Produced", type: "esriFieldTypeInteger"}, ], features: [], geometryType: "" } // loop over depots var depots = ["A", "B"] for(var d in depots) { // find all counts of the depot and order by date var depot = depots[d] var depot_count_fs = Filter(count_fs, "Depot = @depot") depot_count_fs = OrderBy(depot_count_fs, "DateOfCount") // loop over depot_count_fs var previous_depot_count = First(depot_count_fs) for(var depot_count in depot_count_fs) { // find all entries in the production_fs between now and the previous count date var start_date = previous_depot_count.DateOfCount var end_date = depot_count.DateOfCount var production_between_dates = Filter(production_fs, "Depot = @depot AND DateOfProduction > @start_date AND DateOfProduction <= @end_date") // calculate the values var current_count = depot_count.Count var produced = Sum(production_between_dates, "Count") var taken = previous_depot_count.Count + produced - current_count // append new feature Push(out.features, {attributes: {DateOfCount: Number(depot_count.DateOfCount), Depot: depot, Count: current_count, Taken: taken, Produced: produced}}) // set previous previous_depot_count = depot_count } } return Featureset(Text(out))   JohannesLindner_5-1670317869128.png The "Taken" and "Produced" columns list the sandbags that were taken and produced between the current and the previous count.     Lines 1-35 define the sample data. They should of course be replaced with code loading your actual data sets. Depending on your database and date format, you might have to change the sql clause in line 64. ... View more

The ArcGIS Pro label engine has a huge amount of options for placing and fitting labels and for conflict resolution. Take a look at the docs: Label with the Maplex Label Engine—ArcGIS Pro | Documentation   Some solutions for your problem: Use straight label placement instead of horizontal: JohannesLindner_0-1669973204836.png JohannesLindner_1-1669973220423.png   Don't try to place the labels at fixed positions: JohannesLindner_2-1669973511849.png   Play with the Overrun and Reduce Size parameters JohannesLindner_3-1669973607286.png JohannesLindner_4-1669973633452.png Just force all labels to appear JohannesLindner_5-1669973720374.png   JohannesLindner_6-1669973744568.png   ... View more