5 ms·
Finding a minimum polygon of a set of points
- TheLoneWolfling 12y agoI'm probably missing something, because I don't see what's so difficult about this. You can do this as a "rolling" algorithm. Start with three points, defining a triangle. Then repeatedly add the next point. To add a point, first check if it's inside the (convex) polygon found so far. If it is, you're done. Just continue to the next point to add. If it isn't, you walk along both directions along the (convex) polygon until you hit the first point that doesn't violate the convexity constraint. (Just compare slopes of the lines between the next point and the target point, and the target point and the point to be added). Then just splice it in in place of the chunk removed. Now, this isn't particularly fast - O(n^3) I think, but it's still faster than spinning up a LP solver. (Again, I'm probably missing something.)
- Trombone12 12y agoThe solution demonstrated on the page seems to ask for the n-gon that circumscribes n randomly placed points, which isn't the same as the smallest convex hull that contains them. I'm not sure why the polygon should have the same number sides as there are points to contain, especially since the original tweet asks for a convex hull, but there you go.
- TheLoneWolfling 12y agoOh, ok. That makes much more sense.
- Strilanc 12y agoFYI, the algorithm you were describing was to compute a convex hull. It can be done in O(n lg n) time by sorting the points by angle w.r.t. some starting point and then "gift wrapping". You can find where to add a new point in O(lg h) time with a binary search, but you may need to eject O(h) old points. Points can only be ejected once, so there might still be some nice aggregate bound in terms of n. http://en.wikipedia.org/wiki/Convex_hull_algorithms http://en.wikipedia.org/wiki/Convex_hull_algorithms
- laex 12y agoSome time back, I wrote a simple convex hull computing code in Go. If anyone's interested: https://github.com/hemantasapkota/go-convexhull https://github.com/hemantasapkota/go-convexhull
- jbattle 12y agoThe animation on wiki seems broken, but the gift wrapping algorithm is easy to understand and implement and gives good results (though different from what is presented here)... http://en.wikipedia.org/wiki/Gift_wrapping_algorithm http://en.wikipedia.org/wiki/Gift_wrapping_algorithm
- jgunsch 12y agoThe algorithm seems similar to the Graham scan [1], but in either case it doesn't calculate a _regular_ polygon, which makes the problem trickier. [1] http://en.wikipedia.org/wiki/Graham_scan http://en.wikipedia.org/wiki/Graham_scan
- stevenf 12y agoA simple, fast solution to this type of problem involves applying Welzl's algorithm for an expected O(n polylog n) solution.
- TheLoneWolfling 12y agoI thought Welzl's algorithm was for circles, not regular polygons.
- shasta 12y agoWell, there are sets of four points whose minimal containing square is not defined by any three of the points. That makes me suspicious that you're not going to get this plan to work. And playing with the program, I find minimal pentagons that seem to be defined by 5 points.
- mattdesl 12y agoRelated: https://github.com/mikolalysenko/alpha-shape https://github.com/mikolalysenko/alpha-shape
- nilkn 12y agoI played around with this a while ago: http://bl.ocks.org/zpconn/11387143 http://bl.ocks.org/zpconn/11387143
- nilkn 12y agoAs far as I can tell, this is actually finding the smallest regular n-gon which contains n points.
- deleted 12y ago[deleted]
- chaoxu 12y agoThe problem seems to be: find the smallest regular k-gon that covers all the points. (now, assume the total number of points is n, and k = O(n)) One can transform the problem to the following, and solve it in around O(n^7) time. http://dl.acm.org/citation.cfm?id=73853 http://dl.acm.org/citation.cfm?id=73853 There probably exist faster algorithm for this special case using smarter parametric search...