Re: bug in regionprops underlying library qhull?
Hi Stefan, I do not know well the algorithm that qhul implements, and it is lucky to be more complicated than what is needed in this case since it seems to work in higher dimensions than 2D. Nevertheless, I do not see any reason why the points should be pathological: here my modifications:  a/skimage/morphology/convex_hull.py +++ b/skimage/morphology/convex_hull.py @@ 51,9 +51,15 @@ def convex_hull_image(image): 'scipy >= 0.9.') # Find the convex hull  chull = Delaunay(coords).convex_hull  v = coords[np.unique(chull)] + #chull = Delaunay(coords).convex_hull + #v1 = coords[np.unique(chull)] + import vigra + v=vigra.geometry.convexHull(np.copy(coords).astype(np.float32)).view(np.ndarray) + #print v1 + #print v2 + #raise As you can see I just exchanged the computation from scipy with that of another library to which I personally had contributed implementing the standard monotone chain convex hull algorithm. Everything seems fine now. Il giorno venerdì 24 agosto 2012 13:43:43 UTC+2, Stefan van der Walt ha scritto:
Hi Luca
On Wed, Aug 22, 2012 at 1:41 AM, Luca Fiaschi
javascript:> wrote: I confirm the bug to be really in qhull. I substituted the computation of the convex hull in scipy (qhull) with another library (vigranumpy) and thins fix the problem.
Strange, do you know if the points are pathological for some reason? Perhaps we can prefilter them to prevent any crashes. Where did you get these coordinates from?
Stéfan
participants (1)

Luca Fiaschi