Python Concave Hull Polygon of a set of lines

僤鯓⒐⒋嵵緔 提交于 2020-06-12 07:41:29

问题


I'm looking for a python implementation for the Concave Hull problem. My problem is a bit different since I don't have a set of points, but a set of lines, where the result Concave-Hull will roughly bound along the lines (as in the left drawing).

I understand that there is no single "correct answer". But some approximation will be enough for my needs. One possible solution is to take each line and interpolate it to a range of let's say 20 points and find the concave hull of all the created points. Not sure about it.

Edit:

I think the lines add some value making the hull clearer and easier to find.

A good python implementation for the problem, even if not using the lines (just finding a concave hull from a list of points) will also be helpful


回答1:


This is an answer for your subquestion:

A good python implementation for the problem, even if not using the lines (just finding a concave hull from a list of points) will also be helpful

You could use alphashape. The tricky part is to choose an alpha that fits your needs. Alphashape comes with a function to find the optimum alpha value. Basically it starts with 0 (= convex hull) and increases alpha until it starts loosing points. From this optimum value we take 95 %, which is, of course, a rather arbitrary solution, but it'll give you a good approximation in many cases.

import alphashape
import matplotlib.pyplot as plt
from descartes import PolygonPatch

points = [(17, 158),(15, 135),(38, 183),(43, 19),(93, 88),(96, 140),(149, 163),(128, 248),(216, 265),(248, 210),(223, 167),(256, 151),(331, 214),(340, 187),(316, 53),(298, 35),(182, 0),(121, 42)]

alpha = 0.95 * alphashape.optimizealpha(points)
hull = alphashape.alphashape(points, alpha)
hull_pts = hull.exterior.coords.xy

fig, ax = plt.subplots()
ax.scatter(hull_pts[0], hull_pts[1], color='red')
ax.add_patch(PolygonPatch(hull, fill=False, color='green'))

One possible solution is to take each line and interpolate it to a range of let's say 20 points and find the concave hull of all the created points.

This will not give you the desired output as the concave hull will follow these additional (fake) points and it becomes more concave than it can be with the original points.
I think the best solution for the whole problem is to start with the concave hull of the points for the optimum alpha obtained from optimizealpha and then decrease it until your hull doesn't intersect any of your lines as suggested by @sgillen. This can be done similarly to finding the optimum alpha by using a bisection loop with testing any([polygon.crosses(line) for line in lines]).




回答2:


Here is a fine blog post (and associated github repo) on finding the concave hull for a set of points using python.

My recommendation to you is the following. Create a set of points using the endpoints of each line. Then use the linked to code to generate a concave hull for these points, with some guess for the value of alpha. Once this is done you can check if the generated hull intersects any of your lines, and if it does modify alpha. You can make the check for intersection and adjustment automated if you like.

You could also try adding the midpoints of your lines to your set of points, that may decrease the number of alphas you need to try.




回答3:


Although this question might already be answered, here is also my approach:

Like the other ones also said you first have to convert the endpoints of the lines to a list of points.

After that you could need this tidy function from the SciPy Library: scipy.spatial.ConvexHull(). Basically you just pass an numpy-array with the vertices to the function (created with numpy.array()) and it returns a hull-Object

Here is the Documentation

With the .points-attribute you can either get all the Points, or with the .vertices-attribute you can get the indices of the input list which form the hull. Also you can get things like Area or Volume (for 3D shapes) if you are interested in that.

~Okaghana



来源:https://stackoverflow.com/questions/57260352/python-concave-hull-polygon-of-a-set-of-lines

易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!