Convex hull area in Python?
I have a set of points A. I get the convex hull CH_A of A.
Then, I have extra points, point set B. I add B into A and get a bigger point set. I obtain the convex hull CH_AB of this bigger set containing both A and B.
I want to quantify how much I have to pay to add B into set A. I am thinking about using an additional area to quantify this cost.
Say CH_A has an area of Area_A, then CH_AB has an area of Area_AB. Then, I want to calculate the marginal cost as
(Area_AB - Area_A) / Area_A
How may I get the area of the convex hull in Python?
Answer
You could just use the ConvexHull class from scipy.spatial. It will not only give you the hull's area, but it will compute the hull for you as well. But if you do use it, BEWARE! In 2D, the attribute you want to use is not area, it is volume, because the former will actually give you the hull's perimeter.
That's because the attributes are named after their values in 3D, where area will indeed be the hull's area, and volume, well, its volume. For 2D hulls, the names are the same, but what they actually contain is not quite what it says on the tin. Worse, the documentation does not warn you about this.
(You can easily check this with trivial examples such as an isosceles right triangle of length 1 on its legs: the perimeter should be 2+sqrt(2), or about 3.414213562, and the area should be 0.5.)