by
by
Oh gosh, I know so little! Yes, I just labelled the field in the excel sheet, not realizing what an OID was until I did some Googling.
I added a column to my excel sheet:
Which gives me this table in Arc GIS Pro:
I entered the code into the Python pane as such:
def circle_section(x, y, start, end, radius, wkid):
"""Creates a circle section.
x, y: float, coordinates of the center
start, end: int, start and end of the section in degrees, Nort=0°, East=90°
radius: float, radius of the section
wkid: int, well-known id of the coordinate system
returns an arcpy.Polygon in the chosen coordinate system
"""
start = int(start)
end = int(end)
center = arcpy.Point(x, y)
center_geo = arcpy.PointGeometry(center, wkid)
arc_geos = [center_geo.pointFromAngleAndDistance(angle, radius, "PLANAR") for angle in range(start, end, 1)]
points = [center] + [ag.firstPoint for ag in arc_geos] + [center]
return arcpy.Polygon(arcpy.Array(points), spatial_reference=wkid)
def create_view_sheds(in_table, out_features, x_field, y_field, bearing_left_field, bearing_right_field, distance_field, wkid):
"""Creates a polygon feature class with circle sections.
in_table: str, path to the input table or name of the table in the active map
out_features: str, path to the output feature class
*_field: str, names of the fields in the input table
wkid: int, well-known id of the coordinate system
"""
# create output feature class
from pathlib import Path
folder = str(Path(out_features).parent)
name = str(Path(out_features).name)
arcpy.management.CreateFeatureclass(folder, name, "POLYGON", spatial_reference=wkid)
arcpy.management.AddField(out_features, "FID", "LONG")
with arcpy.da.InsertCursor(out_features, ["SHAPE@", "FID"]) as i_cursor:
with arcpy.da.SearchCursor(in_table, ["ID", x_field, y_field, bearing_left_field, bearing_right_field, distance_field]) as s_cursor:
# for each table row
for oid, x, y, bearing_left, bearing_right, distance in s_cursor:
# create a circle section
poly = circle_section(x, y, bearing_left, bearing_right, distance, wkid)
# write the polygon and the table row's objectid into the feature class
i_cursor.insertRow([poly, oid])
create_view_sheds("Observations", "memory/Viewsheds", "X", "Y", "BearingLeft", "BearingRight", "Distance", 25832)The program seems to read it (there are no more error messages) and creates a Viewshed, but... it's blank.
I tried changing the order layer and style to see if it was hidden somehow, but it didn't do anything. Did I edit the code correctly?
by
