The most isolated feature can be obtained by searching for a point, which represents a feature in the Euclidian space, with the highest isolation among the other points that surround it. Isolation criteria of a point can be reformulated as the distance to another nearest point from the set, i.e. the smallest distance among the neighbours. We establish the neighbourhood relationships by using the Delaunay graph that is the dual of the Voronoi graph, which avoids considering points that are not in the proximity.
The complexity of Delaunay computation is limited by O(nlog n) in the worst case. The isolation search is limited by O(nk), where k is the largest number of neighbours in the Voronoi diagram, and, since it can be considered a constant for a large and evenly distributed point set in space, the complexity is instead O(n). Taking into account both operations, we stay within the O(nlog n) limits. An improvement to the constant factor in the isolation search could have been made by caching the distance between any two neighbours, improving the performance of this operation roughly twice-fold.
Python (version 3) language was used for the implementation. It was tempting to use C++, but the overhead of the build setup and the dependencies like the Boost Polygon library was not worth it in the end. We use the Delauney module of the SciPy python library that is in turn based on the Qhull library. NumPy library was used for the algebraic operations and as a dependency. SciPy for Python can be installed by executing:
pip3 install scipy