Minimum enclosing circle/ellipse 2

· Source: Yet Another Math Programming Consultant · Field: Science & Research — Mathematics & Computational Sciences · Depth: Advanced, medium

Summary

This content addresses the problem of finding the minimum enclosing circle for a given set of $n$ circles of varying sizes. It begins by illustrating a dataset of 15 circles with their x, y coordinates and radii, and then introduces a preprocessing step to identify and remove circles entirely contained within another, reducing the dataset. The article presents two Nonlinear Programming (NLP) formulations for this optimization problem, including one that uses a small constant $\varepsilon = 0.0001$ to prevent gradient issues and another that removes the square root by squaring both sides, adding a lower bound constraint for the radius. Finally, it reformulates the problem as a Second-Order Conic Program (SOCP), which introduces additional variables but eliminates the need for an initial point or explicit lower bound on the radius. All three models yield identical solutions, with an optimal enclosing circle having a center at (14.296, 12.600) and a radius of 16.875.

Key takeaway

For a Research Scientist or Software Engineer working on geometric optimization problems, understanding the equivalence and nuances of NLP and SOCP formulations for the minimum enclosing circle problem is crucial. You should consider preprocessing your data to simplify the problem and evaluate both NLP (e.g., CONOPT) and SOCP (e.g., CPLEX) solvers, potentially using one's solution as a warm start for the other to achieve optimal precision.

Key insights

The smallest enclosing circle for a set of circles can be found using NLP or SOCP models.

Principles

Method

The method involves preprocessing to remove contained circles, then solving an NLP or SOCP model to minimize the enclosing circle's radius, subject to constraints ensuring all remaining circles are contained.

In practice

Topics

Best for: Research Scientist, Software Engineer

Related on AIssential

Open in AIssential →

Editorial summary, takeaway, and curation by AIssential. Original article published by Yet Another Math Programming Consultant.