I have an irregular polygon feature that I need to split into x number of equal area features. I need to be able to have the divisions be something other than north/south or east/west. Example: The black polygon in this case needs to be divided into 3 equal area polygons. I manually split these, but I need it to be more precise and the splits need to be aligned in parallel to each other. They need to be able to be calculated with an angle/rotation specification, and not just north/south or east/west. What I have tried: * This excellent code is the closest I have gotten to splitting an irregular polygon up into equal areas. Unfortunately, it only runs for north/south or east/west lines. I can't figure out how to rotate it and I don't have enough coding experience/time to efficiently edit it for my needs. * Fishnet Tool (as suggested by this, this, this, and this) and this has an angle of rotation. However, this creates a grid, but not equal areas. Also, it is for the feature extent not shape. * I cannot use the Parcel Fabrics tool (as suggested by this, this, this, and this) as I don't have the necessary licensing for it. It really blows my mind that this isn't a tool that exists/is readily available. Any advice/suggestions is welcomed. Or, if you really want to be a super hero, edit the code so I can specify an angle for rotation. Thank you very much!

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; } } }

Divide irregular polygon into equal areas by angle

pokateo_

I have an irregular polygon feature that I need to split into x number of equal area features. I need to be able to have the divisions be something other than north/south or east/west.

Example:

The black polygon in this case needs to be divided into 3 equal area polygons. I manually split these, but I need it to be more precise and the splits need to be aligned in parallel to each other. They need to be able to be calculated with an angle/rotation specification, and not just north/south or east/west.

Irregular polygon of equal areas

What I have tried:

This excellent code is the closest I have gotten to splitting an irregular polygon up into equal areas. Unfortunately, it only runs for north/south or east/west lines. I can't figure out how to rotate it and I don't have enough coding experience/time to efficiently edit it for my needs.
Fishnet Tool (as suggested by this, this, this, and this) and this has an angle of rotation. However, this creates a grid, but not equal areas. Also, it is for the feature extent not shape.
I cannot use the Parcel Fabrics tool (as suggested by this, this, this, and this) as I don't have the necessary licensing for it.

It really blows my mind that this isn't a tool that exists/is readily available.

Any advice/suggestions is welcomed. Or, if you really want to be a super hero, edit the code so I can specify an angle for rotation. Thank you very much!

Najděte více příspěvků označených jako

Accepted answers

DanPatterson_Retired

for the hero that wants to take it on, it is far easier to rotate polygon by the desired rotation angle, proceed with Alex's implementation, then rotate result back.

This will work with projected coordinates. Take the polygon, subtract the centroid value from the polygon points (ie a translation to the origin). Rotate by your desired angle (rotation). Do the work (rotate the reverse) and translate back to the original location.

All comments

DanPatterson_Retired

for the hero that wants to take it on, it is far easier to rotate polygon by the desired rotation angle, proceed with Alex's implementation, then rotate result back.

DarrenWiens2

I'd stick to editing Alex Tereshenko's code.

The main thing to edit is the template line coordinates used for cutting the polygons:

pnt_arr<SPAN class="punctuation token">.</SPAN>add<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>x_min<SPAN class="punctuation token">,</SPAN> y_max<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># top left</SPAN>
pnt_arr<SPAN class="punctuation token">.</SPAN>add<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>x_min<SPAN class="punctuation token">,</SPAN> y_min<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># bottom left‍‍‍‍</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Suppose you want a cut line 45 degrees from NW to SE. You'd change the coordinates to:

x_diff <SPAN class="operator token">=</SPAN> math<SPAN class="punctuation token">.</SPAN>tan<SPAN class="punctuation token">(</SPAN>math<SPAN class="punctuation token">.</SPAN>radians<SPAN class="punctuation token">(</SPAN><SPAN class="number token">45</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">*</SPAN> <SPAN class="punctuation token">(</SPAN>y_max <SPAN class="operator token">-</SPAN> y_min<SPAN class="punctuation token">)</SPAN>
pnt_arr<SPAN class="punctuation token">.</SPAN>add<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>x_min <SPAN class="operator token">-</SPAN> x_diff<SPAN class="punctuation token">,</SPAN> y_max<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
pnt_arr<SPAN class="punctuation token">.</SPAN>add<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>x_min<SPAN class="punctuation token">,</SPAN> y_min<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

... then keep trying cutting as you go right. Like picture below (maybe makes it more confusing...).

Unfortunately, I don't believe there's any way to avoid coding on this one. You'd also have to do some extra coding to accommodate all possible line angles.

pokateo_

Thank you so much, Dan!

I don't know why it didn't occur to me to manually rotate it, run the tool, and then rotate it back. That works perfectly. You're the hero I needed but not the one I deserved.

pokateo_

The picture does indeed make it better! Thanks for your thinking on this. I'll play around editing the code with this soon!

GISMallku

For anybody who runs across this thread and is looking for this functionality, please see Split irregular polygons into equal parts‌, add your vote and any specific user requirements in the comments.

This idea is Under Consideration for inclusion as a geoprocessing tool in ArcGIS Pro.

pokateo_

This is great to know! Thanks, Kory. Definitely looking forward to having this integrated.

SergeyTolstov

ArcGIS Pro now has this tool to divide polygons into equal areas:

Subdivide Polygon—Data Management toolbox | ArcGIS Desktop

Připojte se ke konverzaci

Přihlášení členové mohou přispívat, sledovat aktualizace a další. Jste tu noví? Zaregistrujte si bezplatný účet.

Rychlé odkazy

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