Latest Contributions by JohannesLindner

Field names of joined tables have to be called using bracket notation. The resulting  Arcade expression should look like this: $feature["TableName.FieldName"]   ... View more

  def analyze_connections(area_fc, line_fc, area_id_field, line_id_field): """Analyzes how areas are connected by lines. area_fc: str, layer name or path of the areas (polygons) line_fc: str, layer name or path of the lines (polylines) area_id_field: str, uniqe identifier of the areas line_id_field: str, unique identifier of the lines Returns a dict of dicts, specifiying the areas and the count of direct connections to other areas """ areas = [row for row in arcpy.da.SearchCursor(area_fc, ["SHAPE@", area_id_field])] lines = [row for row in arcpy.da.SearchCursor(line_fc, ["SHAPE@", line_id_field])] areas_on_line = {line[1]: [area[1] for area in areas if not line[0].disjoint(area[0])] for line in lines} areas = [area[1] for area in areas] area_connections = dict() for area in areas: area_connections[area] = dict() for other_area in areas: connections = 0 for l, a in areas_on_line.items(): if area != other_area: if area in a and other_area in a: connections += 1 area_connections[area][other_area] = connections return area_connections def compare_connections(old, new): """Compares two connection dictionaries (eg output of analyze_connections). old: dict of dicts, the old connections new: dict of dicts, the new connections Returns a dict of dicts, specifying the difference between the connection counts of old and new connections. """ areas = sorted(list(old.keys())) diff = dict() for area in areas: diff[area] = dict() for other_area in areas: diff[area][other_area] = new[area][other_area] - old[area][other_area] return diff def print_adjacency_matrix(connections): """Takes a dict of dicts (eg output of analyze_connections) and prints a 2D matrix. """" areas = sorted(list(connections.keys())) rows = [[connections[area][other_area] for other_area in areas] for area in areas] for i, row in enumerate(rows): rows[i] = "\t".join([areas[i]] + [str(r) for r in row]) rows.insert(0, "\t".join([""] + areas)) print("\n".join(rows)) old_connections = analyze_connections("TestPolygons", "TestLines", "TextField", "TextField") print_adjacency_matrix(old_connections) # A B C D #A 0 1 0 1 #B 1 0 2 2 #C 0 2 0 1 #D 1 2 1 0 new_connections = analyze_connections("TestPolygons", "TestLines", "TextField", "TextField") print_adjacency_matrix(new_connections) # A B C D #A 0 0 1 0 #B 0 0 1 0 #C 1 1 0 1 #D 0 0 1 0 difference = compare_connections(old_connections, new_connections) print_adjacency_matrix(difference) # A B C D #A 0 -1 1 -1 #B -1 0 -1 -2 #C 1 -1 0 0 #D -1 -2 0 0   ... View more

arcpy.env.workspace = r"C:\gdblocation" fcList = arcpy.ListFeatureClasses("*") for fc in fcList: print(fc) fields = arcpy.ListFields(fc) field_names = [f.name for f in fields] null_counts = {fn: 0 for fn in field_names} with arcpy.da.SearchCursor(fc, field_names) as cursor: for row in cursor: for i, v in enumerate(row): if v is None: null_counts[field_names[i]] += 1 for field in fields: print("{0},{1},{2},{3}".format(field.name, field.type, field.length, field.domain, null_counts[field.name])) ... View more

var gasworks = FeaturesetByName(...) var gasworks_in_state = Intersects($feature, gasworks) var landfills = FeaturesetByName(...) var landfills_in_state = Intersects($feature, landfills) var waste = FeaturesetByName(...) var waste_in_state = Intersects($feature, waste) var point_counts = [ `gasworks: ${Count(gasworks_in_state)}`, `landfills: ${Count(landfills_in_state)}`, `waste: ${Count(waste_in_state)}`, ] return Concatenate(point_counts, TextFormatting.Newline) ... View more

You're almost there. The last step is to determine whether the date you calculated is before or after today: var startDate = $feature["EFFECTIVE_DT"]; var sixyears = DateAdd(startDate,6, 'years') return sixyears < Today()   This will return true for permits that passed the moratorium date, false for active permits.   JohannesLindner_0-1656578743535.png     Then you can change the legend labels: JohannesLindner_1-1656578785670.png   And that's it. JohannesLindner_2-1656578852149.png   ... View more

JohannesLindner_0-1656409948165.png   You need a Standard or Advanced license for that tool. I guess you have a Basic license?   To see your license level:   Open the Project settings JohannesLindner_2-1656410228268.png     Go to the Licensing tab. This might take a while to load. JohannesLindner_1-1656410098676.png   ... View more

Your syntax is OK. You don't need the elses here, because when you return, the rest of the code gets skipped and isn't executed. But that's a small, optional thing. Back() only works on arrays, not on feature sets, I just added an idea to change that. So to return the last feature, you need to sort the featureset. Assuming you have some key field (I use OBJECTID here), you could do it like that: var fsParcel = FeatureSetByName($datastore, "Parcel", ["OBJECTID", "APN", "SitusUnitNumber", "SitusStreetNumber"]) var fsParcelIntersect = OrderBy(Intersects(fsParcel, $feature), "OBJECTID") // loop through parcel features for (var Parcel in fsParcelIntersect) { // unit number is not empty and matches -> return APN if (Parcel.SitusUnitNumber != null && Parcel.SitusUnitNumber ==$feature.Unit){ return Parcel.APN } // unit number is empty, but street number matches -> return APN if (Parcel.SitusUnitNumber == null && Parcel.SitusStreetNumber == $feature.Add_Number){ return Parcel.APN } // both unit number and street number are empty/0 -> return last parcel's APN if (Parcel.SitusUnitNumber == null && Parcel.SitusStreetNumber == "0"){ // to get the last feature, reverse the featureset and call First() var LastParcel = First(OrderBy(fsParcelIntersect, "OBJECTID DESC")) return Iif(IsEmpty(LastParcel), Null, LastParcel.APN) } } // if we land here, none of the intersected parcels has the address's Unit // returns first intersect var Parcel = First(fsParcelIntersect) return Iif(IsEmpty(Parcel), Null, Parcel.APN)   I don't really understand your last check. If both SitusUnitNumber and SitusStreetNumber are empty (I guess that what "0" represents), why do you want to get the last parcel instead of just letting the loop get to the next parcel? ... View more