nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Rotate Polylines based on Existing Polyline Angles

AdamThompson3

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 Information Systems Stack… 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!

Encuentra más publicaciones etiquetadas con

Accepted answers

RandyBurton

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.

polygons

From here, you need to add additional parallel lines and clip them to the polygon.

<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">import</SPAN> math

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">draw_line</SPAN><SPAN class="punctuation token">(</SPAN>point<SPAN class="punctuation token">,</SPAN> distance<SPAN class="punctuation token">,</SPAN> bearing<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; angle <SPAN class="operator token">=</SPAN> <SPAN class="number token">90</SPAN> <SPAN class="operator token">-</SPAN> bearing
&nbsp;&nbsp;&nbsp; bearing <SPAN class="operator token">=</SPAN> math<SPAN class="punctuation token">.</SPAN>radians<SPAN class="punctuation token">(</SPAN>bearing<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; angle <SPAN class="operator token">=</SPAN> math<SPAN class="punctuation token">.</SPAN>radians<SPAN class="punctuation token">(</SPAN>angle<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; cosa <SPAN class="operator token">=</SPAN> math<SPAN class="punctuation token">.</SPAN>cos<SPAN class="punctuation token">(</SPAN>angle<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; cosb <SPAN class="operator token">=</SPAN> math<SPAN class="punctuation token">.</SPAN>cos<SPAN class="punctuation token">(</SPAN>bearing<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; x1<SPAN class="punctuation token">,</SPAN> y1 <SPAN class="operator token">=</SPAN> \
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="punctuation token">(</SPAN>point<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="punctuation token">(</SPAN>distance <SPAN class="operator token">*</SPAN> cosa<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> point<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="punctuation token">(</SPAN>distance <SPAN class="operator token">*</SPAN> cosb<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; x2<SPAN class="punctuation token">,</SPAN> y2 <SPAN class="operator token">=</SPAN> \
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="punctuation token">(</SPAN>point<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">-</SPAN> <SPAN class="punctuation token">(</SPAN>distance <SPAN class="operator token">*</SPAN> cosa<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> point<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">-</SPAN> <SPAN class="punctuation token">(</SPAN>distance <SPAN class="operator token">*</SPAN> cosb<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">return</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">[</SPAN>x1<SPAN class="punctuation token">,</SPAN> y1<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="punctuation token">[</SPAN>x2<SPAN class="punctuation token">,</SPAN> y2<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">]</SPAN>

polyFC <SPAN class="operator token">=</SPAN> <SPAN class="string token">"poly_line"</SPAN> <SPAN class="comment token"># the polygon feature layer</SPAN>
azimuth <SPAN class="operator token">=</SPAN> <SPAN class="number token">347</SPAN>

desc <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe<SPAN class="punctuation token">(</SPAN>polyFC<SPAN class="punctuation token">)</SPAN>
SR <SPAN class="operator token">=</SPAN> desc<SPAN class="punctuation token">.</SPAN>spatialReference
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>outputCoordinateSystem <SPAN class="operator token">=</SPAN> SR

points <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
lines <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>

<SPAN class="keyword token">with</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>polyFC<SPAN class="punctuation token">,</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="string token">'SHAPE@'</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">as</SPAN> cursor<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> row <SPAN class="keyword token">in</SPAN> cursor<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; centroid <SPAN class="operator token">=</SPAN> row<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">.</SPAN>centroid
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; points<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>PointGeometry<SPAN class="punctuation token">(</SPAN>centroid<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; dist <SPAN class="operator token">=</SPAN> <SPAN class="number token">0</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> part <SPAN class="keyword token">in</SPAN> row<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> pnt <SPAN class="keyword token">in</SPAN> part<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; cent_vert_dist <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>PointGeometry<SPAN class="punctuation token">(</SPAN>pnt<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>distanceTo<SPAN class="punctuation token">(</SPAN>centroid<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> cent_vert_dist <SPAN class="operator token">></SPAN> dist<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; dist <SPAN class="operator token">=</SPAN> cent_vert_dist
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># far_point = arcpy.PointGeometry(pnt)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># points.append(far_point)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; feature_info <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN> draw_line<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">(</SPAN>centroid<SPAN class="punctuation token">.</SPAN>X<SPAN class="punctuation token">,</SPAN> centroid<SPAN class="punctuation token">.</SPAN>Y<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> dist<SPAN class="punctuation token">,</SPAN> azimuth<SPAN class="punctuation token">)</SPAN> <SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> feature <SPAN class="keyword token">in</SPAN> feature_info<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lines<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>Polyline<SPAN class="punctuation token">(</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>Array<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN><SPAN class="operator token">*</SPAN>coords<SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">for</SPAN> coords <SPAN class="keyword token">in</SPAN> feature<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>

arcpy<SPAN class="punctuation token">.</SPAN>CopyFeatures_management<SPAN class="punctuation token">(</SPAN>points<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'in_memory\points'</SPAN><SPAN class="punctuation token">)</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>CopyFeatures_management<SPAN class="punctuation token">(</SPAN>lines<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'in_memory\lines'</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

All comments

RandyBurton

The tool you mention creates a hull rectangle for each polygon and makes the lines parallel to what is determined as the longest side of the hull rectangle. Are most or all of your polygons rectangles (and not irregular shapes)? If they are rectangles, would you always want the lines parallel to the most north/south sides?

AdamThompson3

Hey Randy,

Most are rectangles yes with a few irregular shapes. How would I go about changing this so it would always draw the lines parallel to the most north/south sides? I am not sure if that will work but its worth a shot to try and see how it looks visually.

RandyBurton

Here is some of the code from the stackexchange page (the part where it loops through the polygon feature):

<SPAN class="keyword token">for</SPAN> row <SPAN class="keyword token">in</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>polyFC<SPAN class="punctuation token">,</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="string token">"SHAPE@"</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN> spatial_reference<SPAN class="operator token">=</SPAN>SR<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># create inner buffer</SPAN>
&nbsp;&nbsp;&nbsp; polyBuff <SPAN class="operator token">=</SPAN> row<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">.</SPAN>buffer<SPAN class="punctuation token">(</SPAN>buffNum <SPAN class="operator token">*</SPAN> <SPAN class="operator token">-</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># create hull rectangle to establish a rotated area of interest</SPAN>
&nbsp;&nbsp;&nbsp; coordSplit <SPAN class="operator token">=</SPAN> row<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">.</SPAN>hullRectangle<SPAN class="punctuation token">.</SPAN>split<SPAN class="punctuation token">(</SPAN><SPAN class="string token">' '</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># collect corner coordinates</SPAN>
&nbsp;&nbsp;&nbsp; coordList <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Array<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>coordSplit<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN>coordSplit<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>coordSplit<SPAN class="punctuation token">[</SPAN><SPAN class="number token">2</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN>coordSplit<SPAN class="punctuation token">[</SPAN><SPAN class="number token">3</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>coordSplit<SPAN class="punctuation token">[</SPAN><SPAN class="number token">4</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN>coordSplit<SPAN class="punctuation token">[</SPAN><SPAN class="number token">5</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>coordSplit<SPAN class="punctuation token">[</SPAN><SPAN class="number token">6</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN>coordSplit<SPAN class="punctuation token">[</SPAN><SPAN class="number token">7</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>coordSplit<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN>coordSplit<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># create lines from hull rectangle</SPAN>
&nbsp;&nbsp;&nbsp; currentLines <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> pointNum <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">4</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>Array<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN>coordList<SPAN class="punctuation token">.</SPAN>getObject<SPAN class="punctuation token">(</SPAN>pointNum<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN>coordList<SPAN class="punctuation token">.</SPAN>getObject<SPAN class="punctuation token">(</SPAN>pointNum<SPAN class="operator token">+</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; hullRecLine <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Polyline<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Array<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN>coordList<SPAN class="punctuation token">.</SPAN>getObject<SPAN class="punctuation token">(</SPAN>pointNum<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN>coordList<SPAN class="punctuation token">.</SPAN>getObject<SPAN class="punctuation token">(</SPAN>pointNum<SPAN class="operator token">+</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; currentLines<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>hullRecLine<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># compare first and second line to determine if first line is short or long</SPAN>
&nbsp;&nbsp;&nbsp; firstLong <SPAN class="operator token">=</SPAN> <SPAN class="number token">0</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> currentLines<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">.</SPAN>length <SPAN class="operator token">></SPAN> currentLines<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">.</SPAN>length<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; firstLong <SPAN class="operator token">=</SPAN> <SPAN class="number token">1</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># calculate number of points needed along short axis</SPAN>
&nbsp;&nbsp;&nbsp; numPoints <SPAN class="operator token">=</SPAN> int<SPAN class="punctuation token">(</SPAN>math<SPAN class="punctuation token">.</SPAN>floor<SPAN class="punctuation token">(</SPAN>currentLines<SPAN class="punctuation token">[</SPAN>firstLong<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">.</SPAN>length<SPAN class="operator token">/</SPAN>lineSpaceNum<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># create and join points to create parallel lines</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> point <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">,</SPAN>numPoints<SPAN class="operator token">+</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; shortPoint1 <SPAN class="operator token">=</SPAN> currentLines<SPAN class="punctuation token">[</SPAN>firstLong<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">.</SPAN>positionAlongLine<SPAN class="punctuation token">(</SPAN>lineSpaceNum<SPAN class="operator token">*</SPAN>point<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; shortPoint2 <SPAN class="operator token">=</SPAN> currentLines<SPAN class="punctuation token">[</SPAN>firstLong <SPAN class="operator token">+</SPAN> <SPAN class="number token">2</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">.</SPAN>positionAlongLine<SPAN class="punctuation token">(</SPAN>currentLines<SPAN class="punctuation token">[</SPAN>firstLong <SPAN class="operator token">+</SPAN> <SPAN class="number token">2</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">.</SPAN>length <SPAN class="operator token">-</SPAN> <SPAN class="punctuation token">(</SPAN>lineSpaceNum<SPAN class="operator token">*</SPAN>point<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; parallel <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Polyline<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Array<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN>shortPoint1<SPAN class="punctuation token">.</SPAN>centroid<SPAN class="punctuation token">,</SPAN>shortPoint2<SPAN class="punctuation token">.</SPAN>centroid<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> SR<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># intersect parallel lines with buffer</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; parallelBuff <SPAN class="operator token">=</SPAN> parallel<SPAN class="punctuation token">.</SPAN>intersect<SPAN class="punctuation token">(</SPAN>polyBuff<SPAN class="punctuation token">,</SPAN><SPAN class="number token">2</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; parallels<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>parallelBuff<SPAN class="punctuation token">)</SPAN>
<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

In lines 6-17, the code is creating a hull rectangle and getting the lines that make up the rectangle. At lines 20-22, it is comparing the lines to see which side is the long side. The code then goes on to draw the lines.

If you were to compare the starting xy coordinates and ending coordinates of the first and second lines, you should be able to determine which is the most northerly of the lines. Then set the index ( firstLong ) to indicate which line to use.

AdamThompson3

hmm yeah I see what you mean, that could work! Also to note, I just received more data within this grid section where I have been given an ID for the section block and the average azimuth that current polylines that are drawn within it are. For example, BLOCK ID 1, Average Azimuth 347.57. How would I be able to introduce a portion where if the Block ID contains an Azimuth attribute, draw lines with the same spacing currently you set but with the angle that is specified within the table. Any help would be appreciated.

RandyBurton

It should be possible to use the polygon's centroid and draw lines using an azimuth. Several of Dan Patterson‌'s Py... blog‌ articles give code examples that might work into this project (or at least, to provide some ideas).

Geometry: Points in the field calculator

(see: Line direction or Azimuth to Compass Bearing and Convert Azimuth to Compass Bearing sections)

Origin, distances and bearings... geometry wanderings

Numpy Snippets # 3 ... Phish_Nyet ... creating sampling grids using numpy and arcpy

Hexagons, Rectangles and Triangles... Sampling Frameworks

AdamThompson3

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.

DanPatterson_Retired

Adam, that is the code that creates the polyline using the existing points etc.

You could... use angleAndDistanceTo (other, {method}) ...

from Polyline—ArcPy classes | ArcGIS Desktop

to find out what the current angle of the polyline is.

You could use pointFromAngleAndDistance (angle, distance, {method})

to replace your 2nd point.

The only downside I can think of is that you may need to decide whether you are going to rotate the existing polylines from the first point, the middle point or the last point on the existing polyline.

My Sampling Frameworks toolset allows you to create polygons... but in your case, you would want rotated rectangles, but you only want 2 of the side... you could at least pull out my code for a look-see to see if you can modify i... I haven't had any use cases for this sort of thing (ie parallel rotated lines, but I will put it on the 'todo' list

AdamThompson3

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

RandyBurton

polygons

From here, you need to add additional parallel lines and clip them to the polygon.

<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">import</SPAN> math

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">draw_line</SPAN><SPAN class="punctuation token">(</SPAN>point<SPAN class="punctuation token">,</SPAN> distance<SPAN class="punctuation token">,</SPAN> bearing<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; angle <SPAN class="operator token">=</SPAN> <SPAN class="number token">90</SPAN> <SPAN class="operator token">-</SPAN> bearing
&nbsp;&nbsp;&nbsp; bearing <SPAN class="operator token">=</SPAN> math<SPAN class="punctuation token">.</SPAN>radians<SPAN class="punctuation token">(</SPAN>bearing<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; angle <SPAN class="operator token">=</SPAN> math<SPAN class="punctuation token">.</SPAN>radians<SPAN class="punctuation token">(</SPAN>angle<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; cosa <SPAN class="operator token">=</SPAN> math<SPAN class="punctuation token">.</SPAN>cos<SPAN class="punctuation token">(</SPAN>angle<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; cosb <SPAN class="operator token">=</SPAN> math<SPAN class="punctuation token">.</SPAN>cos<SPAN class="punctuation token">(</SPAN>bearing<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; x1<SPAN class="punctuation token">,</SPAN> y1 <SPAN class="operator token">=</SPAN> \
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="punctuation token">(</SPAN>point<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="punctuation token">(</SPAN>distance <SPAN class="operator token">*</SPAN> cosa<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> point<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="punctuation token">(</SPAN>distance <SPAN class="operator token">*</SPAN> cosb<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; x2<SPAN class="punctuation token">,</SPAN> y2 <SPAN class="operator token">=</SPAN> \
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="punctuation token">(</SPAN>point<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">-</SPAN> <SPAN class="punctuation token">(</SPAN>distance <SPAN class="operator token">*</SPAN> cosa<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> point<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">-</SPAN> <SPAN class="punctuation token">(</SPAN>distance <SPAN class="operator token">*</SPAN> cosb<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">return</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">[</SPAN>x1<SPAN class="punctuation token">,</SPAN> y1<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="punctuation token">[</SPAN>x2<SPAN class="punctuation token">,</SPAN> y2<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">]</SPAN>

polyFC <SPAN class="operator token">=</SPAN> <SPAN class="string token">"poly_line"</SPAN> <SPAN class="comment token"># the polygon feature layer</SPAN>
azimuth <SPAN class="operator token">=</SPAN> <SPAN class="number token">347</SPAN>

desc <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe<SPAN class="punctuation token">(</SPAN>polyFC<SPAN class="punctuation token">)</SPAN>
SR <SPAN class="operator token">=</SPAN> desc<SPAN class="punctuation token">.</SPAN>spatialReference
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>outputCoordinateSystem <SPAN class="operator token">=</SPAN> SR

points <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
lines <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>

<SPAN class="keyword token">with</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>polyFC<SPAN class="punctuation token">,</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="string token">'SHAPE@'</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">as</SPAN> cursor<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> row <SPAN class="keyword token">in</SPAN> cursor<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; centroid <SPAN class="operator token">=</SPAN> row<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">.</SPAN>centroid
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; points<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>PointGeometry<SPAN class="punctuation token">(</SPAN>centroid<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; dist <SPAN class="operator token">=</SPAN> <SPAN class="number token">0</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> part <SPAN class="keyword token">in</SPAN> row<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> pnt <SPAN class="keyword token">in</SPAN> part<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; cent_vert_dist <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>PointGeometry<SPAN class="punctuation token">(</SPAN>pnt<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>distanceTo<SPAN class="punctuation token">(</SPAN>centroid<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> cent_vert_dist <SPAN class="operator token">></SPAN> dist<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; dist <SPAN class="operator token">=</SPAN> cent_vert_dist
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># far_point = arcpy.PointGeometry(pnt)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># points.append(far_point)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; feature_info <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN> draw_line<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">(</SPAN>centroid<SPAN class="punctuation token">.</SPAN>X<SPAN class="punctuation token">,</SPAN> centroid<SPAN class="punctuation token">.</SPAN>Y<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> dist<SPAN class="punctuation token">,</SPAN> azimuth<SPAN class="punctuation token">)</SPAN> <SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> feature <SPAN class="keyword token">in</SPAN> feature_info<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lines<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>Polyline<SPAN class="punctuation token">(</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>Array<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN><SPAN class="operator token">*</SPAN>coords<SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">for</SPAN> coords <SPAN class="keyword token">in</SPAN> feature<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>

arcpy<SPAN class="punctuation token">.</SPAN>CopyFeatures_management<SPAN class="punctuation token">(</SPAN>points<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'in_memory\points'</SPAN><SPAN class="punctuation token">)</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>CopyFeatures_management<SPAN class="punctuation token">(</SPAN>lines<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'in_memory\lines'</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

AdamThompson3

Hey Randy, Awesome script!! I am wondering now if this might be easier to implement then my previous one, or try and incorporate pieces of both. If I could keep the option parameter of spacing between lines somewhere that would be ideal. This is awesome though, I did not even think of drawing from a points centroid.

RandyBurton

Just some thoughts...

Distance is measured in the units of the spatial reference so you may need to do some conversion. Depending upon false northing value in the spatial reference, you may wish to adjust your azimuth setting.

For distance between lines, you would need to calculate an x value offset. With the desired distance being one side of a right triangle that is perpendicular to the line, you would use the hypotenuse value as an x offset. Adding/subtracting this value from the x value at the ends of the first line should provide parallel lines at the proper distance

If you have an odd number of lines, you can draw the first line through the centroid. If you want an even number of lines, use half the x offset value to move the first line to the right or left of the centroid and then add the rest of your lines.

Hope this helps.

AdamThompson3

Randy this helps a lot, awesome script! I am tweaking it just a bit where instead of one singular azimuth to go off of, I would like to pull my azimuth value from a field within the polygon shapefile. I have been attempting to use a search cursor to look on that field and then say row=azimuth, but have not been successful! What I get is this...

Here is the portion of code I tweaked;

Perhaps my syntax is a bit wonky somewhere. Thank you again for everything

Also This is what I am trying to pull as my azimuth value

each shape it iterates through should draw lines based on this value.

DanPatterson_Retired

interesting problem... if anyone needs to start lines from scratch, I documented an alternate line forming process to cover some of these ideas and some other thoughts

/blogs/dan_patterson/2019/01/17/transect-lines-parallel-lines-offset-lines

AdamThompson3

Hey Dan, wow this awesome!! Pretty sure exactly the last step I need, but I am struggling a bit with integrating it into my current script. So currently I am iterating through each grid shape, finding the center point and drawing a line based on a azimuth that is given to me. That is all good, now I would like to take that current line and replicate it X number of times and a X distance from it.

I see that your scipt indeed does this! I am not sure If I understand all the logic and code that is happening, there is a lot going on and I am not sure If I am to use all of it or only pieces. Given that I have already created a line, can I Incorporate some portions of your script where I am still iterating through each grid creating a line and replicating it. So it shows something like this (replicates created manually)

I am essentially trying to keep it all within one loop and trying to incorporate your code in this area of mine;

Thank you again for all the help, I am still looking at your code and trying to understand how I can implement it into mine. I am using the script Randy provided above with just a slight modification where I am grabbing the azimuth attribute from the table of the grid shapefile.

DanPatterson_Retired

Adam, I would stick with Randy's since it is designed to work with the azimuth within a table.

Mine is designed to work from scratch given a start location, known distance and azimuths.

It is in essence the line version of my...

Sampling Grids for ArcGIS Pro

I have a point version in my Point Tools toolbox, also on the code sharing site (search using Dan_Patterson should bring up a list

AdamThompson3

Thanks Dan, I will take a look and see what I can find, thank you!

AdamThompson3

Have been trying to replicate an inside buffer and intersect lines within that buffer while iterating through it all but keep running into this error... Not sure what is occurring.

Here is the code so far;

Once again have to thank you for all the help and cooperation with me!

RandyBurton

Here is the code I've been working with:

import arcpy
import math

def rotate(origin, point, angle):
    """
    Rotate a point counterclockwise by a given angle around a given origin.
    The angle should be given in radians.
    """
    ox, oy = origin
    px, py = point

    qx = ox + math.cos(angle) * (px - ox) - math.sin(angle) * (py - oy)
    qy = oy + math.sin(angle) * (px - ox) + math.cos(angle) * (py - oy)
    return qx, qy

def make_lines(origin, num_lines, length, spacing, bearing):
    '''
    origin = (x,y), num_lines = number of lines to draw
    length = 1/2 length of line, spacing = spacing between lines
    bearing = angle to rotate points
    '''
    angle = math.radians(360.0 - bearing)
    lines = []
    x, y = origin
    # set up line points before rotation
    xyz = x - (spacing*(num_lines-1))/2.0
    y1 = y + length*1.0
    y2 = y - length*1.0
    # rotate points around origin
    for lns in range(1,num_lines+1):
        qx1,qy1 = rotate((x,y),(xyz,y1),angle)
        qx2,qy2 = rotate((x,y),(xyz,y2),angle)
        lines.append([qx1,qy1,qx2,qy2])
        xyz += spacing
    return lines

# set input/output parameters
polyFC = arcpy.GetParameterAsText(0)      # input polygons : feature layer, input
angleField = arcpy.GetParameterAsText(1)  # field containing rotation: field, input (obtained from param0)
numLines = arcpy.GetParameterAsText(2)    # number of lines: long, input
lineSpacing = arcpy.GetParameterAsText(3) # line spacing : linear unit, input
buffDist = arcpy.GetParameterAsText(4)    # inner buffer distance : linear unit, input
outLines = arcpy.GetParameterAsText(5) # output parallel lines : feature class, output

desc = arcpy.Describe(polyFC)
SR = desc.spatialReference # get spatial reference of polygon feature
arcpy.env.overwriteOutput = True # output will be overwritten
arcpy.env.outputCoordinateSystem = SR # set default spatial reference

D2M = { # dictionary for converting distance to meters
    'CENTIMETERS': .01,
    'DECIMALDEGREES': 1.0, # no conversion
    'DECIMETERS': 0.1,
    'FEET': 0.3048, # 0.304800609601 Foot_US
    'INCHES': 0.0254,
    'KILOMETERS': 1000.0,
    'METERS': 1.0, # no conversion
    'MILES': 1609.344,
    'MILLIMETERS': 0.001,
    'NAUTICALMILES': 1852.0,
    'POINTS': 0.000352778,
    'UNKNOWN': 1.0, # no conversion, tool may select data frame default
    'YARDS': 0.9144,
    }

# parse numbers from parameters
numLines = int(numLines)

lineSpacing = lineSpacing.split(' ')
if lineSpacing[1].upper() not in D2M.keys():
    lineSpacing[1] = 'UNKNOWN'
spacing = float(lineSpacing[0]) * D2M[lineSpacing[1].upper()] * SR.metersPerUnit # convert spacing unit to meters, then to SR units

buffDist = buffDist.split(' ')
if buffDist[1].upper() not in D2M.keys():
    buffDist[1] = 'UNKNOWN'
buffNum = float(buffDist[0]) * D2M[buffDist[1].upper()] * SR.metersPerUnit # convert buffer unit to meters, then to SR units

# list to save lines
lines = []

with arcpy.da.SearchCursor(polyFC,['SHAPE@', angleField]) as cursor:
    for row in cursor:
        # buffer for trimming lines to polygon
        polyBuff = row[0].buffer(buffNum * -1)

        # get centroid of polygon
        centroid = row[0].centroid

        # find distance from centroid to farthest point, this will be used as 1/2 line length     
        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

        # pass centroid, dist for 1/2 line length, and other variables to make lines
        mkLns = make_lines(((centroid.X, centroid.Y)), numLines, dist, spacing, row[1])

        # make polyline and trim to buffer
        for ln in mkLns:
            pl = arcpy.Polyline(arcpy.Array([arcpy.Point(ln[0],ln[1]), arcpy.Point(ln[2],ln[3])]), SR)
            bufLn = pl.intersect(polyBuff,2)
            lines.append(bufLn)

# save lines to feature class
arcpy.CopyFeatures_management(lines,outLines)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

I added some code to better deal with the linear units; the original tool ignored the type of unit selection and used the default from the spatial reference. I create points for the ends of the lines and then rotate the points using the polygon's centroid as the rotation origin. The length for the lines (before buffering) is twice the distance from the centroid to the far point.

I've modified and added to the tool parameters:

Script tool dialog

And the results:

Line results

Hope this helps.

AdamThompson3

Hey Randy!

This looks almost exactly what I was trying to achieve, when I am in the office tomorrow I will give it a go! Thank you very much for this! I have been reading it over and researching a lot of the stuff you used so I can understand it fully! Outside of forum postings here, and if you're comfortable of course do you have any other form of communication I could possibly send questions should I have them? I am not sure of GeoNet's messaging or other communication means it has!

Thank you again.

Cheers

Adam

matthewgarner

I am just getting into the Python scripting world, so my understanding of all the logic is very basic at this point, but I have found the script you put together extremely helpful. The script works wonderful on solid polygons but I am having a issue when it comes to polygons with holes it is giving me this error message, if you have any idea on how I could correct this or lead me in the right direction I would really appreciate it.

-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

Thank you

POLYGON_WITH_HOLE_(FT).zip

RandyBurton

I haven't worked with polygons containing holes. From a quick search of Geonet and the web, I didn't find a solution for buffering a polygon with holes as you wish to do. I'm not sure how the intersect tool would deal with holes; perhaps it could remove portions of the lines where the hole is located (but it wouldn't include a buffer). Since this is an old question marked answered, I suggest you pose a new question with a link to this one.

Únete a la Conversación

Los miembros registrados pueden publicar, seguir actualizaciones y más. ¿Nuevo aquí? Regístrate gratis.

Enlaces Rápidos

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.