Hey Everyone,
Currently what I am trying to do is rotate polylines in certain grid sections to the same orientation of existing polylines in that grid. This seems to be quite complicated as I have looked quite extensively for something like this, but to no success. Currently right now I have implemented a tool that was posted here arcpy - Drawing parallel lines inside polygons (Well Paths) using ArcGIS Desktop? - Geographic Infor... now this works for me and does what it should... but here is my problem posted below.
The grid lines that were drawn automatically (yellow) are horizontal where as I would like them be vertical. If there is a way through arcpy or other tools where I can rotate these lines to be the same direction of polylines that are currently there. This is just one example, this happens many times in different areas.
Thank you!
Accepted Solutions
Here's some code that I have been experimenting with that may give you some ideas. It takes a polygon feature class and creates a layer of points marking each polygon's centroid. It also creates a line through the centroid at a given azimuth.
From here, you need to add additional parallel lines and clip them to the polygon.
import arcpy
import math
def draw_line(point, distance, bearing):
angle = 90 - bearing
bearing = math.radians(bearing)
angle = math.radians(angle)
cosa = math.cos(angle)
cosb = math.cos(bearing)
x1, y1 = \
(point[0] + (distance * cosa), point[1] + (distance * cosb))
x2, y2 = \
(point[0] - (distance * cosa), point[1] - (distance * cosb))
return [[x1, y1], [x2, y2]]
polyFC = "poly_line" # the polygon feature layer
azimuth = 347
desc = arcpy.Describe(polyFC)
SR = desc.spatialReference
arcpy.env.overwriteOutput = True
arcpy.env.outputCoordinateSystem = SR
points = []
lines = []
with arcpy.da.SearchCursor(polyFC,['SHAPE@']) as cursor:
for row in cursor:
centroid = row[0].centroid
points.append(arcpy.PointGeometry(centroid))
dist = 0
for part in row[0]:
for pnt in part:
cent_vert_dist = arcpy.PointGeometry(pnt).distanceTo(centroid)
if cent_vert_dist > dist:
dist = cent_vert_dist
# far_point = arcpy.PointGeometry(pnt)
# points.append(far_point)
feature_info = [ draw_line((centroid.X, centroid.Y), dist, azimuth) ]
for feature in feature_info:
lines.append(
arcpy.Polyline(
arcpy.Array([arcpy.Point(*coords) for coords in feature])))
arcpy.CopyFeatures_management(points,'in_memory\points')
arcpy.CopyFeatures_management(lines,'in_memory\lines')Hey Randy, those are helpful thank you for that! I have kind of narrowed it down to where I need to change the code.
if I could somehow change it where, If I know the azimuth is 347 for example for lines drawn in that section, then I know the lines should be oriented N-S with that degree on them. This has proved to be extremely challenging, at least for my level of python knowledge.
Hey Dan thanks for the suggestions, I will play around with those two functions and see what I can do. I am thinking of not even rotating them, If I could just draw them based on azimuth direction from the start that would be great. I have been provided a table that outlines the grid, each section in that grid has its own ID and the average azimuth within each section. For example
If I could draw the lines from the start with these azimuth bearings that would be helpful. I have been exploring the other python modules Shapely, and even was able to create bearing lines within QGIS. But I am not sure how to automate it further. I appreciate all the help none the less and am always happy to hear suggestions and ideas. Thanks again
