by
1 Solution
Accepted Solutions
So, assuming you have a table like this:
Display the locations
- right-click the table in the table of contents, "Display XY Data"
- Choose where you want to save the point feature class, choose the appropriate coordinate fields and your coordinate system, click OK.
Display the view polygons
This is a little more complicated.
You could probably do that using a combination of geoprocessing tools, but I didn't find an easy one.
This is a case where a small Python script gets the job done faster.
- Open the Python window
- Copy/Paste these functions and hit enter twice
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)
angles = range(start, end, 1)
if start > end:
angles = list(range(start, 360, 1)) + list(range(0, end, 1))
arc_geos = [center_geo.pointFromAngleAndDistance(angle, radius, "PLANAR") for angle in angles]
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, ["OID@", 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])
- Call the function
create_view_sheds("Observations", "memory/Viewsheds", "X", "Y", "BearingLeft", "BearingRight", "Distance", 25832)
So, assuming you have a table like this:
Display the locations
- right-click the table in the table of contents, "Display XY Data"
- Choose where you want to save the point feature class, choose the appropriate coordinate fields and your coordinate system, click OK.
Display the view polygons
This is a little more complicated.
You could probably do that using a combination of geoprocessing tools, but I didn't find an easy one.
This is a case where a small Python script gets the job done faster.
- Open the Python window
- Copy/Paste these functions and hit enter twice
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)
angles = range(start, end, 1)
if start > end:
angles = list(range(start, 360, 1)) + list(range(0, end, 1))
arc_geos = [center_geo.pointFromAngleAndDistance(angle, radius, "PLANAR") for angle in angles]
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, ["OID@", 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])
- Call the function
create_view_sheds("Observations", "memory/Viewsheds", "X", "Y", "BearingLeft", "BearingRight", "Distance", 25832)
by
by
by
by
It's like a never ending sudoku problem!
Ok, so I adjusted the title of my table and it now looks like this:
I now have an error that says:
Traceback (most recent call last):
File "<string>", line 1, in <module>
File "<string>", line 37, in create_view_sheds
RuntimeError: Cannot find field 'OID@'
If I try and break it down, is it saying that for this line:
for oid, x, y, bearing_left, bearing_right, distance in s_cursor:
It cannot find a field labelled 'OID@"?
I cannot say thank you enough to reflect my gratitude for the help on this. Hopefully one day my skills will be good enough to help someone else out!
Edit: I think the issue might be because I'm using a simple.csv file from Excel to upload the data. I'll have to figure out how to properly assign it an ObjectID.
