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

Thiessen Polygons Not Generating Continuous Tessellations in 10.3

GerryGabrisch

I noticed that the Create Thiessen Polygon Tool in ArcGIS 10.3 is not producing a continuous irregular tessellation across the analysis area.

The resulting polygons have gaps and overlaps as confirmed by topology! I have not noticed this on other versions of ArcGIS.

Can anyone confirm that this is happening on their installation of ArcGIS 10.3? Is there a fix, or does anyone have a script to replace Create Thiessen Polygons that works correctly?

Find more posts tagged with

Accepted answers

DanPatterson_Retired

Clam Art

Well, you are going to have one H of a time trying to do that point distribution all at once.

In short, my toolbox didn't generate errors...don't know why yours did.

Shapefile in ... shapefile out as indicated in the tool help (for a litany of reasons)

The algorithm is Steve Fortune's original in case Ghislain is interested....I kept it old school without modern data structures, since it is a thing of beauty.

I just did one portion and had no issues except for the fact that

the distribution is hardly ideal for Voronoi or Delaunay triangulations
the coordinates are Geographic making the results essentially useless since the triangles etc are not geodesic

having said that, here are a few images and a package for you to examine. If you can play with the data and show me where you were getting the gaps, I can examine further (hope the mpk works... if not, I will zip and ship)

ADDENDUM

I am also guessing that is the coordinate system aka Geographic that is causing the problem

clam.mpk.zip

All comments

DanPatterson_Retired

Gerry My ArcScripts mirror... Triangulation tools, shapefiles as input and output, should work, let me know if it doesn't

GerryGabrisch

Thanks, but it threw an error.

DanPatterson_Retired

A duplicate point ... I will have to check version numbers since I downloaded it from the beta site and I had corrected that, so if you have duplicate points you will have to wait...sorry

GerryGabrisch

I see that there are spatially duplicate points but with different attributes! That might be a problem.

I have contacted the data creator to sort out the differences.

Create Thiessen Polygons did not care if there were spatially duplicate points. I guess if you are going to be doing some coding, and want to produce a better tool, let the user know there are spatially duplicate points as a guard against such a thing.

Thanks for throwing the error triangulation tools!

curtvprice

How large is the gap?

GerryGabrisch

Curtis, in this case the gaps and overlaps are a few feet over polygons that are 10 to 12 feet across.

I did get confirmation from ESRI that this is a bug and will be fixed in 10.4.

DanPatterson_Retired

Gerry ... coming soon...I am having to port from python 2 to python 3, so it affected anything using "cmp" which was abandoned in version 2 so it can be used in ArcGIS Pro as well.

The duplicate point issue has also been solved regardless of input feature geometry. If duplicates are removed, I have had no problems with all the test cases I have run it on. I am assuming that you are using points as input? If so, just remove duplicates. I am also making it more generic so that it doesn't matter what the input geometry is. All features are converted to points, then all duplicates removed.

So if you are in a pinch, remove your duplicate points. In any event, I will post my version long before 10.4 comes out.

gprince

Hey Gerry,

can you share the gpk of your original case that produces gaps?

Or open a ticket with tech support. We'd like to look at this case.

thanks

GerryGabrisch

Ghislain,

Thanks for looking into the issue.

I have attached the gpk to this response.

CreateThiessenPolygons.gpk.zip

DanPatterson_Retired

Clam Art

Well, you are going to have one H of a time trying to do that point distribution all at once.

In short, my toolbox didn't generate errors...don't know why yours did.

Shapefile in ... shapefile out as indicated in the tool help (for a litany of reasons)

The algorithm is Steve Fortune's original in case Ghislain is interested....I kept it old school without modern data structures, since it is a thing of beauty.

I just did one portion and had no issues except for the fact that

the distribution is hardly ideal for Voronoi or Delaunay triangulations
the coordinates are Geographic making the results essentially useless since the triangles etc are not geodesic

ADDENDUM

I am also guessing that is the coordinate system aka Geographic that is causing the problem

clam.mpk.zip

GerryGabrisch

Thanks for your help, Dan. Rest assured, I have a polygon layer that represents the bounding area of the survey that I will use to clip the Voronoi polygons with.

Projected - Geographic coordinate system? I had no idea this would make such a big difference in the resulting polygons. I have projected the data and re-created the Voronoi polygons and the gap and overlap problems went away. Best to read the fine print!

Problem Solved.

Join the Conversation

Signed in members can post, follow updates, and more. New here? Register a free account.

Quick Links

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