Intersection Properties of Finite Disk Collections

1. Introduction

One of the new techniques developed for the analysis of large clusters of information, known as Big Data, is Topological Data Analysis (TDA). In TDA, simplicial complexes associated with the data are constructed. These structures include the Vietoris-Rips complex, the Čech complex, and the piecewise linear lower star complex, among others. Of special interest to us is the generalized Čech complex structure. Although the standard Čech complex is formed by intersecting a collection of disks with a fixed radius, the generalized version allows varying radii. This flexibility enables us to highlight specific data points by assigning or weighting them with larger and/or more rapidly expanding balls to them, while de-emphasizing others by using smaller and/or slower growing balls. This approach proves valuable for handling noisy data sets, offering an alternative to discarding data that may not meet a specific significance threshold [1].

Understanding the patterns of intersections and the timing of intersections among a set of disks in ${\mathbb{R}}^{\mathrm{d}}$, each with potentially different radii, is a fundamental problem. This leads to the exploration of the generalized Čech complex structure, which captures the intersection information of these disks, regardless of their radii. Rescaling the radii by the same factor, we obtain a filtered generalized Čech complex, where the associated simplicial complexes evolve as the scale parameter varies. In particular, in [2], algorithms are provided to calculate the generalized Čech complex in ${\mathbb{R}}^2$, and [3] presents an algorithm to determine the Čech scale for a collection of disks in the plane.

To establish the necessary foundation for our study, Section 2 introduces crucial concepts and notation that will be used throughout the article and we focus on analyzing the intersection of a disk system in ${\mathbb{R}}^{\mathrm{d}}$. We start by investigating the intersection of two disks in Section 2.1 and then expand our analysis to a system of m disks in Section 2.2. By applying Helly’s Theorem, we prove that it is sufficient to examine the intersection of all subsystems consisting of $\mathrm{d}+1$ disks in order to determine if the system has a nonempty intersection.

In Section 3, we define Vietoris-Rips systems and Čech systems, together with presenting results regarding the Rips scale and the Čech scale, as well as their connections. In Section 3.1, we present an algorithm that can determine whether the intersection of the system is empty or non-empty. This is achieved by exclusively computing the poles of subsystems of disks (or spheres). Finally, in Section 3.2, we introduce the algorithm that computes an approximation to the Čech scale using the numerical bisection method.

Additionally, in Section 4 we incorporate the concept of an minimal axis-aligned bounding box (AABB) into our methodology. An AABB is a rectangular parallelepiped whose faces are perpendicular to the basis vectors. These bounding boxes frequently arise in spatial subdivision problems, such as ray tracing [4] and collision detection [5]. In this paper, we study AABBs to enclose the intersection of a finite collection of disks. This approach proves valuable for discerning whether the collection intersects at a singular point or not. In this section, we also provide an algorithm for constructing the AABB of a disk system.

2. Intersection properties of sphere systems

Throughout this work, we will refer to a d-disk system M, or simply a disk system, as a finite collection of closed disks in ${\mathbb{R}}^{\mathrm{d}}$ with positive and not necessarily equal radii, i.e.,

$$M=\{{\mathrm{D}}_{\mathrm{i}}({\mathrm{c}}_{\mathrm{i}};{\mathrm{r}}_{\mathrm{i}})\subset {\mathbb{R}}^{\mathrm{d}}\mid {\mathrm{c}}_{\mathrm{i}}\in {\mathbb{R}}^{\mathrm{d}},{\mathrm{r}}_{\mathrm{i}}>0,1\le \mathrm{i}\le \mathrm{m}<\infty \}.$$

Moreover, in order to study the intersection properties of a disk system M with the approach addressed in Section 3 and Section 4 of this work, we will conduct a study in this section of the intersection properties of the spheres corresponding to the boundaries of each disk in M, which we call a sphere system and denote by $\partial M$,

$$\partial M=\{\partial {\mathrm{D}}_{\mathrm{i}}\subset {\mathbb{R}}^{\mathrm{d}}\mid {\mathrm{D}}_{\mathrm{i}}\in M\},$$

where ∂ denotes the topological boundary operator.

Following the notation in [6], we introduce the following generalization of the sphere.

Definition 1.

An i-sphere in ${\mathbb{R}}^{\mathrm{d}}$ is the intersection of a sphere with an affine subspace of dimension i.

Of course, the notions of a sphere (as a $(\mathrm{d}-1)$-dimensional surface) and a d-sphere in ${\mathbb{R}}^{\mathrm{d}}$ agree. However, an i-sphere in ${\mathbb{R}}^{\mathrm{d}}$ can also be viewed as the intersection of d-spheres. For instance, the intersection of two spheres typically occurs in a hyperplane, forming a $(\mathrm{d}-1)$-sphere in ${\mathbb{R}}^{\mathrm{d}}$. When another d-sphere intersects this configuration, the result may be a $(\mathrm{d}-1)$-sphere, a $(\mathrm{d}-2)$-sphere, a 0-sphere (a single point), or it might even be empty, all within the same hyperplane. For a disk system $M=\left\{{\mathrm{D}}_{\mathrm{i}}({\mathrm{c}}_{\mathrm{i}};{\mathrm{r}}_{\mathrm{i}})\right\}$ composed of m disks, where $\{{\mathrm{c}}_1,\dots ,{\mathrm{c}}_{\mathrm{m}}\}$ is a set in general position in ${\mathbb{R}}^{\mathrm{d}}$, the maximum dimension of the affine subspace associated with the i-sphere, obtained from the intersection of all the spheres in $\partial M$, is at most $\mathrm{d}-\mathrm{m}+1$, or equivalently, $\mathrm{i}=\mathrm{d}-\mathrm{m}+1$. This conclusion is drawn from [6, Theorem 2.1] and the fact that the affine hull of $\{{\mathrm{c}}_1,\dots ,{\mathrm{c}}_{\mathrm{m}}\}$ is of dimension $\mathrm{m}-1$. Consequently, the following result is proven.

Lemma 1.

Let $M=\{{\mathrm{D}}_1({\mathrm{c}}_1;{\mathrm{r}}_1),\dots ,{\mathrm{D}}_{\mathrm{m}}({\mathrm{c}}_{\mathrm{m}};{\mathrm{r}}_{\mathrm{m}})\}$ be disk system such that $\{{\mathrm{c}}_1,\dots ,{\mathrm{c}}_{\mathrm{m}}\}$ is a set in general position in ${\mathbb{R}}^{\mathrm{d}}$. Then, the possibilities for the set ${\cap }_{{\mathrm{D}}_{\mathrm{i}}\in M}\partial {\mathrm{D}}_{\mathrm{i}}$ are:

• the empty set;
• a single point;
• a $(\mathrm{d}-\mathrm{m}+1)$-sphere.

Remarkable points in i-spheres that will play a key role in the rest of the article are the poles. Let ${\mathrm{\pi }}_{\mathrm{i}}:{\mathbb{R}}^{\mathrm{d}}⟶\mathbb{R}$ be the canonical projection on the i-th factor for $\mathrm{i}=1,...,\mathrm{d}$, and let $\{{\mathrm{e}}_1,{\mathrm{e}}_2,...,{\mathrm{e}}_{\mathrm{d}}\}$ be the standard basis of ${\mathbb{R}}^{\mathrm{d}}$.

Definition 2.

Let ${\mathrm{e}}_{\mathrm{q}}$ be the q-th vector of the canonical base of ${\mathbb{R}}^{\mathrm{d}}$. An ${\mathrm{e}}_{\mathrm{q}}$-north (south) pole of an i-sphere S in ${\mathbb{R}}^{\mathrm{d}}$ is a point on S whose projection on the q-th coordinate is maximum (minimum). In other words, $\mathrm{x}\in \mathrm{S}$ is the ${\mathrm{e}}_{\mathrm{q}}$-north pole if ${\mathrm{\pi }}_{\mathrm{q}}(\mathrm{y})\le {\mathrm{\pi }}_{\mathrm{q}}(\mathrm{x})$ for all $\mathrm{y}\in \mathrm{S}-\left\{\mathrm{x}\right\}$, where ${\mathrm{\pi }}_{\mathrm{q}}$ represents the projection onto the q-th coordinate.

We denote the ${\mathrm{e}}_{\mathrm{q}}$-north pole of S by ${\mathrm{s}}_{\mathrm{q}}^{+}$ and the ${\mathrm{e}}_{\mathrm{q}}$-south pole by ${\mathrm{s}}_{\mathrm{q}}^{-}$.

An i-sphere can have a single ${\mathrm{e}}_{\mathrm{q}}$-pole (north or south) or an infinite number of them, which occurs when a normal vector to the affine space containing the i-sphere is aligned with the vector ${\mathrm{e}}_{\mathrm{q}}$. We are interested in finding the ${\mathrm{e}}_{\mathrm{q}}$-poles of $(\mathrm{d}-\mathrm{m}+1)$-spheres originating from disk systems $M=\{{\mathrm{D}}_1({\mathrm{c}}_1;{\mathrm{r}}_1),\dots ,{\mathrm{D}}_{\mathrm{m}}({\mathrm{c}}_{\mathrm{m}};{\mathrm{r}}_{\mathrm{m}})\}$, by taking the intersection ${\cap }_{\mathrm{j}=1}^{\mathrm{m}}\partial {\mathrm{D}}_{\mathrm{j}}$. Such $(\mathrm{d}-\mathrm{m}+1)$-spheres will be denoted by ${\mathrm{S}}_M(\mathrm{c};\mathrm{r})$, to emphasize the disk system M, as well as its center and radius.

Lemma 2.

Let $M=\{{\mathrm{D}}_1,\dots ,{\mathrm{D}}_{\mathrm{m}}\}$ be a d-disk system such that ${\bigcap }_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{D}}_{\mathrm{j}}\ne \varnothing$, and let p be a point in ${\bigcap }_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{D}}_{\mathrm{j}}$ such that ${\mathrm{\pi }}_{\mathrm{q}}(\mathrm{p})\le {\mathrm{\pi }}_{\mathrm{q}}(\mathrm{x})$ (resp. ${\mathrm{\pi }}_{\mathrm{q}}(\mathrm{p})\ge {\mathrm{\pi }}_{\mathrm{q}}(\mathrm{x})$) for every x in ${\bigcap }_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{D}}_{\mathrm{j}}$. Then, there exists an i-sphere $\mathrm{S}=\partial {\mathrm{D}}_{{\mathrm{j}}_1}\cap \dots \cap \partial {\mathrm{D}}_{{\mathrm{j}}_{\mathrm{i}}}$ such that p is in S and p is the ${\mathrm{e}}_{\mathrm{q}}$-south pole (resp. ${\mathrm{e}}_{\mathrm{q}}$-north pole) of S.

Proof. 

Since ${\cap }_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{D}}_{\mathrm{j}}\ne \varnothing$, then $\partial ({\cap }_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{D}}_{\mathrm{j}})\ne \varnothing$, $\partial ({\cap }_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{D}}_{\mathrm{j}})\subset {\cap }_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{D}}_{\mathrm{j}}$ due to the closedness of the sets ${\mathrm{D}}_{\mathrm{j}}$, for $\mathrm{j}=1,\dots ,\mathrm{m}$, and $\mathrm{p}\in \partial ({\cap }_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{D}}_{\mathrm{j}})$.

On the other hand, since $\partial ({\cap }_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{D}}_{\mathrm{j}})\subset {\cup }_{\mathrm{i}=1}^{\mathrm{m}}\partial {\mathrm{D}}_{\mathrm{j}}$, there exist indices ${\mathrm{j}}_1,\dots ,{\mathrm{j}}_{\mathrm{i}}$ such that $\mathrm{p}\in \partial {\mathrm{D}}_{{\mathrm{j}}_{\mathrm{r}}}$ for any $\mathrm{r}=1,\dots ,\mathrm{i}$; let $\mathrm{\Lambda }(\mathrm{p})=\{{\mathrm{j}}_1,\dots ,{\mathrm{j}}_{\mathrm{i}}\}\subseteq \{1,\dots ,\mathrm{m}\}$ be a maximal subset of indices such that $\mathrm{p}\in {\mathrm{D}}_{\mathrm{j}}$ if and only if $\mathrm{j}\in \mathrm{\Lambda }(\mathrm{p})$. We claim that $\mathrm{p}\in {\cap }_{\mathrm{r}=1}^{\mathrm{i}}\partial {\mathrm{D}}_{{\mathrm{j}}_{\mathrm{r}}}$ is the ${\mathrm{e}}_{\mathrm{q}}$-south pole of $\mathrm{S}:={\cap }_{\mathrm{r}=1}^{\mathrm{i}}\partial {\mathrm{D}}_{{\mathrm{j}}_{\mathrm{r}}}$.

In effect, let ${\mathrm{V}}_{\mathrm{p}}\subset {\mathbb{R}}^{\mathrm{d}}$ be an open neighborhood of p sufficiently small such that:

• Every $\mathrm{x}\in {\mathrm{V}}_{\mathrm{p}}\cap \partial ({\cap }_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{D}}_{\mathrm{j}})$ has as maximal set of indices a proper subset of $\mathrm{\Lambda }(\mathrm{p})$,
• $\mathrm{S}\cap {\mathrm{V}}_{\mathrm{p}}\subset \partial ({\cap }_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{D}}_{\mathrm{j}})$.

The first condition can be guaranteed by the finiteness of the disk system M, and the second condition is a consequence of the maximality of the set $\mathrm{\Lambda }(\mathrm{p})$. Therefore, ${\mathrm{\pi }}_{\mathrm{q}}(\mathrm{p})\le {\mathrm{\pi }}_{\mathrm{q}}(\mathrm{x})$ for every $\mathrm{x}\in \mathrm{S}\cap {\mathrm{V}}_{\mathrm{p}}$, which is equivalent to the fact that ${\mathrm{\pi }}_{\mathrm{q}}(\mathrm{p})\le {\mathrm{\pi }}_{\mathrm{q}}(\mathrm{x})$ for every $\mathrm{x}\in \mathrm{S}$, in the case of i-spheres.

   □

2.1. Sphere systems with two spheres

In the following two lemmas we provide the computations to determine the center, radius and poles for a $(\mathrm{d}-1)$-sphere given by the intersection of two disks in ${\mathbb{R}}^{\mathrm{d}}$.

Lemma 3.

Let $M=\{{\mathrm{D}}_1({\mathrm{c}}_1;{\mathrm{r}}_1),{\mathrm{D}}_2({\mathrm{c}}_2;{\mathrm{r}}_2)\}$ be a disk system with two d-disks such that $\partial {\mathrm{D}}_1\cap \partial {\mathrm{D}}_2$ is a $(\mathrm{d}-1)$-sphere $\mathrm{S}={\mathrm{S}}_M(\mathrm{c};\mathrm{r})$ with center c and radius r. Then,

$$\mathrm{c}=\frac{1}{2}\left(1+\frac{{\mathrm{r}}_2^2-{\mathrm{r}}_1^2}{\parallel {\mathrm{c}}_2-{\mathrm{c}}_1{\parallel }^2}\right){\mathrm{c}}_1+\frac{1}{2}\left(1+\frac{{\mathrm{r}}_1^2-{\mathrm{r}}_2^2}{\parallel {\mathrm{c}}_2-{\mathrm{c}}_1{\parallel }^2}\right){\mathrm{c}}_2$$

and

$$\mathrm{r}=\frac{2\sqrt{\mathrm{s}(\mathrm{s}-\parallel {\mathrm{c}}_2-{\mathrm{c}}_1\parallel )(\mathrm{s}-{\mathrm{r}}_1)(\mathrm{s}-{\mathrm{r}}_2)}}{\parallel {\mathrm{c}}_2-{\mathrm{c}}_1\parallel }$$

where $\mathrm{s}=\frac{1}{2}(\parallel {\mathrm{c}}_2-{\mathrm{c}}_1\parallel +{\mathrm{r}}_1+{\mathrm{r}}_2)$.

Proof. 

Let $\mathrm{\Pi }$ be the hyperplane containing the $(\mathrm{d}-1)$-sphere S, which is defined by the equation:

$$\sum _{\mathrm{i}=1}^{\mathrm{d}}({\mathrm{k}}_{\mathrm{i}}-{\mathrm{h}}_{\mathrm{i}}){\mathrm{x}}_{\mathrm{i}}-\frac{1}{2}\sum _{\mathrm{i}=1}^{\mathrm{d}}({\mathrm{k}}_{\mathrm{i}}^2-{\mathrm{h}}_{\mathrm{i}}^2)=\frac{{\mathrm{r}}_1^2-{\mathrm{r}}_2^2}{2},$$

(1)

where ${\mathrm{c}}_1=({\mathrm{h}}_1,\dots ,{\mathrm{h}}_{\mathrm{d}})$ and ${\mathrm{c}}_2=({\mathrm{k}}_1,\dots ,{\mathrm{k}}_{\mathrm{d}})$. Then the normal vector of the hyperplane $\mathrm{\Pi }$ is given by $\mathrm{N}:={\mathrm{c}}_2-{\mathrm{c}}_1=({\mathrm{k}}_1-{\mathrm{h}}_1,\dots ,{\mathrm{k}}_{\mathrm{d}}-{\mathrm{h}}_{\mathrm{d}})$, and the center c of S is determined by the intersection point of the hyperplane $\mathrm{\Pi }$ with the perpendicular line that passes through the center ${\mathrm{c}}_1$ of ${\mathrm{D}}_1$. This line can be parameterized as $\mathrm{\gamma }:\mathrm{t}\mapsto {\mathrm{c}}_1+\mathrm{t}\mathrm{N}=({\mathrm{x}}_1(\mathrm{t}),\dots ,{\mathrm{x}}_{\mathrm{d}}(\mathrm{t}))$, such that $\mathrm{\gamma }(0)={\mathrm{c}}_1$ and $\mathrm{\gamma }(1)={\mathrm{c}}_2$. We can compute the intersection point $\mathrm{c}=\mathrm{\gamma }({\mathrm{t}}_{*})$ of $\mathrm{\Pi }$ and $\mathrm{\gamma }([0,1])$, for any ${\mathrm{t}}_{*}\in [0,1]$, by substituting it in (1),

$$\sum _{\mathrm{i}=1}^{\mathrm{d}}({\mathrm{k}}_{\mathrm{i}}-{\mathrm{h}}_{\mathrm{i}})({\mathrm{h}}_{\mathrm{i}}+{\mathrm{t}}_{*}({\mathrm{k}}_{\mathrm{i}}-{\mathrm{h}}_{\mathrm{i}}))=\frac{1}{2}\left(\sum _{\mathrm{i}=1}^{\mathrm{d}}({\mathrm{k}}_{\mathrm{i}}^2-{\mathrm{h}}_{\mathrm{i}}^2)+({\mathrm{r}}_1^2-{\mathrm{r}}_2^2)\right).$$

And solving for ${\mathrm{t}}_{*}$, we obtain that ${\mathrm{t}}_{*}=\frac{1}{2}+\frac{{\mathrm{r}}_1^2-{\mathrm{r}}_2^2}{2\parallel {\mathrm{c}}_2-{\mathrm{c}}_1{\parallel }^2}$. Hence, the center of S is given by:

$$\begin{array}{c}\hfill \mathrm{c}={\mathrm{c}}_1+{\mathrm{t}}_{*}\mathrm{N}=\frac{1}{2}\left(1+\frac{{\mathrm{r}}_2^2-{\mathrm{r}}_1^2}{\parallel {\mathrm{c}}_2-{\mathrm{c}}_1{\parallel }^2}\right){\mathrm{c}}_1+\frac{1}{2}\left(1+\frac{{\mathrm{r}}_1^2-{\mathrm{r}}_2^2}{\parallel {\mathrm{c}}_2-{\mathrm{c}}_1{\parallel }^2}\right){\mathrm{c}}_2\end{array}$$

Next, we will compute the radius r of S. This radius can be determined as the height r of the triangle with base $\parallel {\mathrm{c}}_2-{\mathrm{c}}_1\parallel$ formed by the points ${\mathrm{c}}_1$, ${\mathrm{c}}_2$, and a point on S. Thus, by the Heron’s formula we have that

$$\mathrm{r}=\frac{2\sqrt{\mathrm{s}(\mathrm{s}-\parallel {\mathrm{c}}_2-{\mathrm{c}}_1\parallel )(\mathrm{s}-{\mathrm{r}}_1)(\mathrm{s}-{\mathrm{r}}_2)}}{\parallel {\mathrm{c}}_2-{\mathrm{c}}_1\parallel },$$

where $\mathrm{s}=\frac{1}{2}(\parallel {\mathrm{c}}_2-{\mathrm{c}}_1\parallel +{\mathrm{r}}_1+{\mathrm{r}}_2)$ correspond to the semi-perimeter.    □

We can proceed now to compute the poles of the $(\mathrm{d}-1)$-sphere $\partial {\mathrm{D}}_1\cap \partial {\mathrm{D}}_2$.

Lemma 4.

Let ${\mathrm{D}}_1({\mathrm{c}}_1;{\mathrm{r}}_1)$ and ${\mathrm{D}}_2({\mathrm{c}}_2;{\mathrm{r}}_2)$ be two d-disks such that $\partial {\mathrm{D}}_1\cap \partial {\mathrm{D}}_2$ is a $(\mathrm{d}-1)$-sphere $\mathrm{S}=\mathrm{S}(\mathrm{c};\mathrm{r})$ with center c and radius r. Then, the ${\mathrm{e}}_{\mathrm{q}}$-poles of S are ${\mathrm{s}}_{\mathrm{q}}^{\pm }=\mathrm{c}\pm {\sum }_{\mathrm{i}}^{\mathrm{d}}{\mathrm{x}}_{\mathrm{i}}{\mathrm{e}}_{\mathrm{i}}$, where

$${\mathrm{x}}_{\mathrm{i}}=\left\{\begin{array}{cc}{\displaystyle \frac{\mathrm{r}|{\mathrm{\pi }}_{\mathrm{i}}({\mathrm{c}}_2-{\mathrm{c}}_1){\mathrm{\pi }}_{\mathrm{q}}({\mathrm{c}}_2-{\mathrm{c}}_1)|}{\parallel {\mathrm{c}}_2-{\mathrm{c}}_1\parallel \sqrt{\parallel {\mathrm{c}}_2-{\mathrm{c}}_1{\parallel }^2-{\mathrm{\pi }}_{\mathrm{q}}{({\mathrm{c}}_2-{\mathrm{c}}_1)}^2}}},\hfill & \mathrm{i}\ne \mathrm{q}\hfill \\ \\ {\displaystyle \frac{\mathrm{r}\sqrt{\parallel {\mathrm{c}}_2-{\mathrm{c}}_1{\parallel }^2-{\mathrm{\pi }}_{\mathrm{q}}{({\mathrm{c}}_2-{\mathrm{c}}_1)}^2}}{\parallel {\mathrm{c}}_2-{\mathrm{c}}_1\parallel }},\hfill & \mathrm{i}=\mathrm{q}.\hfill \end{array}\right.$$

Proof. 

For simplicity, we translate the hyperplane $\mathrm{\Pi }$, which contains the $(\mathrm{d}-1)$-sphere S, as well as the sphere itself, to the origin; in such case, the corresponding equations are given by,

$$\begin{array}{cc}\hfill \sum _{\mathrm{i}=1}^{\mathrm{d}}({\mathrm{k}}_{\mathrm{i}}-{\mathrm{h}}_{\mathrm{i}}){\mathrm{x}}_{\mathrm{i}}& =0,\hfill \\ \hfill \sum _{\mathrm{i}=1}^{\mathrm{d}}{\mathrm{x}}_{\mathrm{i}}^2& ={\mathrm{r}}^2,\hfill \end{array}$$

where ${\mathrm{h}}_{\mathrm{i}}:={\mathrm{\pi }}_{\mathrm{i}}({\mathrm{c}}_1)$ and ${\mathrm{k}}_{\mathrm{i}}:={\mathrm{\pi }}_{\mathrm{i}}({\mathrm{c}}_2)$ for $\mathrm{i}=1,\dots ,\mathrm{d}$. In the case that ${\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}}={\mathrm{\pi }}_{\mathrm{q}}({\mathrm{c}}_2-{\mathrm{c}}_1)=0$, the normal vector $\mathrm{N}={\mathrm{c}}_2-{\mathrm{c}}_1$ of the hyperplane $\mathrm{\Pi }$ is orthogonal to the basis vector ${\mathrm{e}}_{\mathrm{q}}$. Therefore, the ${\mathrm{e}}_{\mathrm{q}}$-poles of S are $\mathrm{c}\pm \mathrm{r}{\mathrm{e}}_{\mathrm{q}}$, which agree with the formulae of the lemma.

On the other hand, suppose that ${\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}}\ne 0$. To find the ${\mathrm{e}}_{\mathrm{q}}$-poles of S, we will use the Lagrange multiplier method. Consider the following function:

$${\mathrm{x}}_{\mathrm{q}}=\mathrm{f}({\mathrm{x}}_1,{\mathrm{x}}_2,...,\widehat{{\mathrm{x}}_{\mathrm{q}}},...,{\mathrm{x}}_{\mathrm{d}})=\frac{-{\sum }_{\mathrm{j}\ne \mathrm{q}}^{\mathrm{d}}({\mathrm{k}}_{\mathrm{j}}-{\mathrm{h}}_{\mathrm{j}}){\mathrm{x}}_{\mathrm{j}}}{{\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}}},$$

(2)

subject to the restriction:

$$\mathrm{g}({\mathrm{x}}_1,{\mathrm{x}}_2,...,\widehat{{\mathrm{x}}_{\mathrm{q}}},...,{\mathrm{x}}_{\mathrm{d}})=\sum _{\mathrm{j}\ne \mathrm{q}}^{\mathrm{d}}{\mathrm{x}}_{\mathrm{j}}^2+{\left(\frac{-{\sum }_{\mathrm{j}\ne \mathrm{q}}^{\mathrm{d}}({\mathrm{k}}_{\mathrm{j}}-{\mathrm{h}}_{\mathrm{j}}){\mathrm{x}}_{\mathrm{j}}}{{\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}}}\right)}^2-{\mathrm{r}}^2=0$$

Let $\mathrm{\lambda }$ be the Lagrange multiplier, we define

$$\mathrm{h}({\mathrm{x}}_1,...,\widehat{{\mathrm{x}}_{\mathrm{q}}},...,{\mathrm{x}}_{\mathrm{d}},\mathrm{\lambda })=\mathrm{f}({\mathrm{x}}_1,...,\widehat{{\mathrm{x}}_{\mathrm{q}}},...,{\mathrm{x}}_{\mathrm{d}})+\mathrm{\lambda }\mathrm{g}({\mathrm{x}}_1,...,\widehat{{\mathrm{x}}_{\mathrm{q}}},...,{\mathrm{x}}_{\mathrm{d}})$$

For any $\mathrm{i}\ne \mathrm{q}$, consider the following system of equations:

$$\frac{\partial \mathrm{h}}{\partial {\mathrm{x}}_{\mathrm{i}}}=-\frac{{\mathrm{k}}_{\mathrm{i}}-{\mathrm{h}}_{\mathrm{i}}}{{\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}}}+2\mathrm{\lambda }{\mathrm{x}}_{\mathrm{i}}+2\mathrm{\lambda }\left(\frac{-{\sum }_{\mathrm{j}\ne \mathrm{q}}^{\mathrm{d}}({\mathrm{k}}_{\mathrm{j}}-{\mathrm{h}}_{\mathrm{j}}){\mathrm{x}}_{\mathrm{j}}}{{\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}}}\right)\left(-\frac{{\mathrm{k}}_{\mathrm{i}}-{\mathrm{h}}_{\mathrm{i}}}{{\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}}}\right)=0.$$

Then

$$\begin{array}{cc}\hfill -\frac{{\mathrm{k}}_{\mathrm{i}}-{\mathrm{h}}_{\mathrm{i}}}{{\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}}}+2\mathrm{\lambda }\left({\mathrm{x}}_{\mathrm{i}}+\frac{{\mathrm{k}}_{\mathrm{i}}-{\mathrm{h}}_{\mathrm{i}}}{{({\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}})}^2}\sum _{\mathrm{j}\ne \mathrm{q}}^{\mathrm{d}}({\mathrm{k}}_{\mathrm{j}}-{\mathrm{h}}_{\mathrm{j}}){\mathrm{x}}_{\mathrm{j}}\right)& =0\hfill \\ \hfill -({\mathrm{k}}_{\mathrm{i}}-{\mathrm{h}}_{\mathrm{i}})({\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}})+2\mathrm{\lambda }\left({({\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}})}^2{\mathrm{x}}_{\mathrm{i}}+({\mathrm{k}}_{\mathrm{i}}-{\mathrm{h}}_{\mathrm{i}})\sum _{\mathrm{j}\ne \mathrm{q}}^{\mathrm{d}}({\mathrm{k}}_{\mathrm{j}}-{\mathrm{h}}_{\mathrm{j}}){\mathrm{x}}_{\mathrm{j}}\right)& =0\hfill \end{array}$$

Solving this system of equations for $\mathrm{\lambda }$, we obtain that,

$$\mathrm{\lambda }=\frac{({\mathrm{k}}_{\mathrm{i}}-{\mathrm{h}}_{\mathrm{i}})({\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}})}{2\left({({\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}})}^2{\mathrm{x}}_{\mathrm{i}}+({\mathrm{k}}_{\mathrm{i}}-{\mathrm{h}}_{\mathrm{i}}){\sum }_{\mathrm{j}\ne \mathrm{q}}^{\mathrm{d}}({\mathrm{k}}_{\mathrm{j}}-{\mathrm{h}}_{\mathrm{j}}){\mathrm{x}}_{\mathrm{j}}\right)}$$

Comparing the last expression for two indices $\mathrm{i}\ne \tilde{\mathrm{i}}$, we have that,

$$\begin{array}{cc}\hfill {\mathrm{x}}_{\mathrm{i}}^2& =\frac{{\mathrm{r}}^2{({\mathrm{k}}_{\mathrm{i}}-{\mathrm{h}}_{\mathrm{i}})}^2}{{({\mathrm{k}}_{\tilde{\mathrm{i}}}-{\mathrm{h}}_{\tilde{\mathrm{i}}})}^2+{\sum }_{\mathrm{j}\ne \tilde{\mathrm{i}},\mathrm{q}}{({\mathrm{k}}_{\mathrm{j}}-{\mathrm{h}}_{\mathrm{j}})}^2+\frac{1}{{({\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}})}^2}{\left({\sum }_{\mathrm{j}\ne \mathrm{q}}{({\mathrm{k}}_{\mathrm{j}}-{\mathrm{h}}_{\mathrm{j}})}^2\right)}^2}\hfill \\ \hfill & =\frac{{\mathrm{r}}^2{({\mathrm{k}}_{\mathrm{i}}-{\mathrm{h}}_{\mathrm{i}})}^2}{{\sum }_{\mathrm{j}\ne \mathrm{q}}{({\mathrm{k}}_{\mathrm{j}}-{\mathrm{h}}_{\mathrm{j}})}^2+\frac{1}{{({\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}})}^2}{\left({\sum }_{\mathrm{j}\ne \mathrm{q}}{({\mathrm{k}}_{\mathrm{j}}-{\mathrm{h}}_{\mathrm{j}})}^2\right)}^2}\hfill \\ \hfill & =\frac{{\mathrm{r}}^2{({\mathrm{k}}_{\mathrm{i}}-{\mathrm{h}}_{\mathrm{i}})}^2{({\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}})}^2}{\left({({\mathrm{k}}_{\mathrm{q}}-{\mathrm{h}}_{\mathrm{q}})}^2+{\sum }_{\mathrm{j}\ne \mathrm{q}}{({\mathrm{k}}_{\mathrm{j}}-{\mathrm{h}}_{\mathrm{j}})}^2\right)\left({\sum }_{\mathrm{j}\ne \mathrm{q}}{({\mathrm{k}}_{\mathrm{j}}-{\mathrm{h}}_{\mathrm{j}})}^2\right)}\hfill \\ \hfill & =\frac{{\mathrm{r}}^2{\mathrm{\pi }}_{\mathrm{i}}{({\mathrm{c}}_2-{\mathrm{c}}_1)}^2{\mathrm{\pi }}_{\mathrm{q}}{({\mathrm{c}}_2-{\mathrm{c}}_1)}^2}{\parallel {\mathrm{c}}_2-{\mathrm{c}}_1{\parallel }^2\left(\parallel {\mathrm{c}}_2-{\mathrm{c}}_1{\parallel }^2-{\mathrm{\pi }}_{\mathrm{q}}{({\mathrm{c}}_2-{\mathrm{c}}_1)}^2\right)}\hfill \end{array}$$

Finally, for $\mathrm{i}=\mathrm{q}$, we can use the last expression to substitute it in (2) and obtain the desired result.    □

2.2. Sphere systems with more than two spheres

Now, let us proceed with the explicit calculation of the coefficients for the center c of the $(\mathrm{d}-\mathrm{m}+1)$-sphere $\mathrm{S}={\cap }_{\mathrm{j}=1}^{\mathrm{m}}\partial {\mathrm{D}}_{\mathrm{j}}$. We can achieve this by considering the disk system translated to ${\mathrm{c}}_{\mathrm{m}}$, denoted as ${\left\{{\mathrm{D}}_{\mathrm{j}}({\mathrm{c}}_{\mathrm{j}}-{\mathrm{c}}_{\mathrm{m}};{\mathrm{r}}_{\mathrm{j}})\right\}}_{\mathrm{j}=1}^{\mathrm{m}}$, and by defining the $(\mathrm{d}-\mathrm{m}+1)$-sphere $\mathrm{S}-\left\{{\mathrm{c}}_{\mathrm{m}}\right\}={\cap }_{\mathrm{j}=1}^{\mathrm{m}}\partial {\mathrm{D}}_{\mathrm{j}}({\mathrm{c}}_{\mathrm{j}}-{\mathrm{c}}_{\mathrm{m}};{\mathrm{r}}_{\mathrm{j}})$. This sphere is positioned at the intersection of hyperplanes (for more details, refer to [6]).

$${({\mathrm{c}}_{\mathrm{k}}-{\mathrm{c}}_{\mathrm{m}})}^{\mathrm{T}}\mathrm{x}=\frac{1}{2}({\mathrm{r}}_{\mathrm{m}}^2+\parallel {\mathrm{c}}_{\mathrm{k}}-{\mathrm{c}}_{\mathrm{m}}{\parallel }^2-{\mathrm{r}}_{\mathrm{k}}^2)$$

(3)

for all $\mathrm{k}=1,...,\mathrm{m}-1$. Utilizing the information that the center of $\mathrm{S}-\left\{{\mathrm{c}}_{\mathrm{m}}\right\}$ can be expressed as a combination of the centers ${\mathrm{c}}_{\mathrm{k}}-{\mathrm{c}}_{\mathrm{m}}$ and substituting it into (3), we obtain a linear system of equations with dimensions $(\mathrm{m}-1)\times (\mathrm{m}-1)$:

$$\sum _{\mathrm{j}=1}^{\mathrm{m}-1}{\mathrm{\lambda }}_{\mathrm{j}}({\mathrm{c}}_{\mathrm{k}}-{\mathrm{c}}_{\mathrm{m}})·({\mathrm{c}}_{\mathrm{j}}-{\mathrm{c}}_{\mathrm{m}})=\frac{1}{2}({\mathrm{r}}_{\mathrm{m}}^2+\parallel {\mathrm{c}}_{\mathrm{k}}-{\mathrm{c}}_{\mathrm{m}}{\parallel }^2-{\mathrm{r}}_{\mathrm{k}}^2)$$

for $\mathrm{k}=1,...,\mathrm{m}-1$. Solving the system of equations for $\mathrm{\lambda }=({\mathrm{\lambda }}_1,...,{\mathrm{\lambda }}_{\mathrm{m}-1})$, we find the center of S as follows:

$$\mathrm{c}={\mathrm{\lambda }}_1({\mathrm{c}}_1-{\mathrm{c}}_{\mathrm{m}})+...+{\mathrm{\lambda }}_{\mathrm{m}-1}({\mathrm{c}}_{\mathrm{m}-1}-{\mathrm{c}}_{\mathrm{m}})+{\mathrm{c}}_{\mathrm{m}}$$

The radius of the sphere S can be computed using the equation:

$${\mathrm{r}}^2={\mathrm{r}}_{\mathrm{k}}^2-{\parallel \mathrm{c}-{\mathrm{c}}_{\mathrm{k}}\parallel }^2$$

for any $\mathrm{k}\in \{1,...,\mathrm{m}-1\}$.

Now that we have determined the center and radius of S, as well as the affine space that contains it, we can proceed to compute its ${\mathrm{e}}_{\mathrm{q}}$-poles for each $\mathrm{q}\in \{1,2,...,\mathrm{d}\}$. These poles reside in the affine space that contains S and within a set that we define below.

Let S be an i-sphere in ${\mathbb{R}}^{\mathrm{d}}$, and let ${\mathrm{n}}_1,...,{\mathrm{n}}_{\mathrm{d}-\mathrm{i}}$ be orthogonal vectors to the affine space L that contains S. Consider the space M generated by these vectors together with the vector ${\mathrm{e}}_{\mathrm{q}}$ from the canonical basis of ${\mathbb{R}}^{\mathrm{d}}$. Let us denote ${\mathrm{n}}_{\mathrm{j}}={\left({\mathrm{n}}_{\mathrm{l}}^{(\mathrm{j})}\right)}_{\mathrm{l}=1}^{\mathrm{d}}$ for each $\mathrm{j}=1,...,\mathrm{d}-\mathrm{i}$. Then, we can define ${\mathrm{L}}_0$, the set L translated to the origin, as follows:

$$\begin{array}{cc}\hfill {\mathrm{L}}_0& ={〈{\left\{{\mathrm{n}}_{\mathrm{j}}\right\}}_{\mathrm{j}=1}^{\mathrm{d}-\mathrm{i}}〉}^{\perp }\hfill \\ \hfill & =\left\{\mathrm{x}\in {\mathbb{R}}^{\mathrm{d}}\mid \mathrm{x}·{\mathrm{n}}_{\mathrm{j}}=0,\forall \mathrm{j}=1,...,\mathrm{d}-\mathrm{i}\right\}\hfill \end{array}$$

The set M is defined as:

$$\begin{array}{cc}\hfill M& =〈{\left\{{\mathrm{n}}_{\mathrm{j}}\right\}}_{\mathrm{j}=1}^{\mathrm{d}-\mathrm{i}}\cup \left\{{\mathrm{e}}_{\mathrm{q}}\right\}〉\hfill \\ \hfill & =\left\{\mathrm{x}=\sum _{\mathrm{j}=1}^{\mathrm{d}-\mathrm{i}}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_{\mathrm{j}}+{\mathrm{\lambda }}_{\mathrm{d}-\mathrm{i}+1}{\mathrm{e}}_{\mathrm{q}}\in {\mathbb{R}}^{\mathrm{d}}\mid {\mathrm{\lambda }}_{\mathrm{j}}\in \mathbb{R},\forall \mathrm{j}=1,...,\mathrm{d}-\mathrm{i}+1\right\}\hfill \\ \hfill & =\left\{\left(\sum _{\mathrm{j}=1}^{\mathrm{d}-\mathrm{i}}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_1^{(\mathrm{j})},...,\sum _{\mathrm{j}=1}^{\mathrm{d}-\mathrm{i}}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_{\mathrm{q}}^{(\mathrm{j})}+{\mathrm{\lambda }}_{\mathrm{d}-\mathrm{i}+1},...,\sum _{\mathrm{j}=1}^{\mathrm{d}-\mathrm{i}}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_{\mathrm{d}}^{(\mathrm{j})}\right)\mid {\mathrm{\lambda }}_{\mathrm{j}}\in \mathbb{R}\right\}\hfill \end{array}$$

Refer to Figure 1 for a visual representation of the subspaces L and $M+\left\{\mathrm{c}\right\}$.

As mentioned above, the ${\mathrm{e}}_{\mathrm{q}}$-poles of S lie at the intersection of L and $M+\left\{\mathrm{c}\right\}$, where c is the center of S. To simplify the calculations, we will utilize ${\mathrm{L}}_0$ and M, and then translate them into c. The intersection of M and ${\mathrm{L}}_0$ can be expressed as follows:

$$\begin{array}{cc}\hfill M\cap {\mathrm{L}}_0& =\left\{\mathrm{x}\in M\mid \mathrm{x}·{\mathrm{n}}_{\mathrm{k}}=0,\forall \mathrm{k}=1,...,\mathrm{d}-\mathrm{i}\right\}\hfill \\ \hfill & =\left\{\sum _{\mathrm{j}=1}^{\mathrm{d}-\mathrm{i}}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_{\mathrm{j}}+{\mathrm{\lambda }}_{\mathrm{d}-\mathrm{i}+1}{\mathrm{e}}_{\mathrm{q}}\in {\mathbb{R}}^{\mathrm{d}}\left|\left(\sum _{\mathrm{j}=1}^{\mathrm{d}-\mathrm{i}}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_{\mathrm{j}}+{\mathrm{\lambda }}_{\mathrm{d}-\mathrm{i}+1}{\mathrm{e}}_{\mathrm{q}}\right)·{\mathrm{n}}_{\mathrm{k}}=0,{\mathrm{\lambda }}_{\mathrm{j}}\in \mathbb{R}\forall \mathrm{k}=1,...,\mathrm{d}-\mathrm{i}\right.\right\}\hfill \\ \hfill & =\left\{\sum _{\mathrm{j}=1}^{\mathrm{d}-\mathrm{i}}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_{\mathrm{j}}+{\mathrm{\lambda }}_{\mathrm{d}-\mathrm{i}+1}{\mathrm{e}}_{\mathrm{q}}\in {\mathbb{R}}^{\mathrm{d}}\left|\sum _{\mathrm{j}=1}^{\mathrm{d}-\mathrm{i}}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_{\mathrm{j}}·{\mathrm{n}}_{\mathrm{k}}+{\mathrm{\lambda }}_{\mathrm{d}-\mathrm{i}+1}{\mathrm{n}}_{\mathrm{q}}^{\mathrm{k}}=0,{\mathrm{\lambda }}_{\mathrm{j}}\in \mathbb{R}\forall \mathrm{k}=1,...,\mathrm{d}-\mathrm{i}\right.\right\}\hfill \end{array}$$

Let us consider a disk system in ${\mathbb{R}}^{\mathrm{d}}$, denoted ${\left\{{\mathrm{D}}_{\mathrm{j}}({\mathrm{c}}_{\mathrm{j}};{\mathrm{r}}_{\mathrm{j}})\right\}}_{\mathrm{j}=1}^{\mathrm{m}}$, where $\mathrm{m}<\mathrm{d}$. The intersection of their boundaries forms a $(\mathrm{d}-\mathrm{m}+1)$-sphere S. In this case, the subspace M has dimension m, or $\mathrm{d}\mathrm{i}\mathrm{m}(M)=\mathrm{m}-1$ if ${\mathrm{e}}_{\mathrm{q}}\in M$. We choose the normal vectors for the affine space containing S as ${\mathrm{n}}_{\mathrm{j}}={\mathrm{c}}_{\mathrm{j}}-{\mathrm{c}}_{\mathrm{m}}$, where $\mathrm{j}=1,...,\mathrm{m}-1$. Then

$$M=\left\{\left(\sum _{\mathrm{j}=1}^{\mathrm{m}-1}{\mathrm{\lambda }}_{\mathrm{j}}\left({\mathrm{c}}_1^{(\mathrm{j})}-{\mathrm{c}}_1^{(\mathrm{m})}\right),...,\sum _{\mathrm{j}=1}^{\mathrm{m}-1}{\mathrm{\lambda }}_{\mathrm{j}}\left({\mathrm{c}}_{\mathrm{q}}^{(\mathrm{j})}-{\mathrm{c}}_{\mathrm{q}}^{(\mathrm{m})}\right)+{\mathrm{\lambda }}_{\mathrm{m}},...,\sum _{\mathrm{j}=1}^{\mathrm{m}-1}{\mathrm{\lambda }}_{\mathrm{j}}\left({\mathrm{c}}_{\mathrm{d}}^{(\mathrm{j})}-{\mathrm{c}}_{\mathrm{d}}^{(\mathrm{m})}\right)\right)\left|{\mathrm{\lambda }}_{\mathrm{j}}\in \mathbb{R}\right.\right\}$$

By rewriting, we have

$$\begin{array}{cc}\hfill M\cap {\mathrm{L}}_0& =\left\{\sum _{\mathrm{j}=1}^{\mathrm{m}-1}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_{\mathrm{j}}+{\mathrm{\lambda }}_{\mathrm{m}}{\mathrm{e}}_{\mathrm{q}}\in {\mathbb{R}}^{\mathrm{d}}\left|\sum _{\mathrm{j}=1}^{\mathrm{m}-1}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_{\mathrm{j}}·{\mathrm{n}}_{\mathrm{k}}+{\mathrm{\lambda }}_{\mathrm{m}}{\mathrm{n}}_{\mathrm{q}}^{\mathrm{k}}=0,{\mathrm{\lambda }}_{\mathrm{j}}\in \mathbb{R},\forall \mathrm{k}=1,...,\mathrm{m}-1\right.\right\}\hfill \end{array}$$

If $\mathrm{S}(\mathrm{c};\mathrm{r})={\cap }_{\mathrm{j}=1}^{\mathrm{m}}\partial {\mathrm{D}}_{\mathrm{j}}$ is the $(\mathrm{d}-\mathrm{m}+1)$-sphere with center in c and radius r, then the ${\mathrm{e}}_{\mathrm{q}}$-poles of S are the ${\mathrm{e}}_{\mathrm{q}}$-poles of $\mathrm{S}-\left\{{\mathrm{c}}_{\mathrm{m}}\right\}$ but translated by c. The poles of $\mathrm{S}-\left\{{\mathrm{c}}_{\mathrm{m}}\right\}$ are located in $M\cap {\mathrm{L}}_0$. If p is an ${\mathrm{e}}_{\mathrm{q}}$-pole of $\mathrm{S}-\left\{{\mathrm{c}}_{\mathrm{m}}\right\}$, then it can be expressed as

$$\mathrm{p}=\sum _{\mathrm{j}=1}^{\mathrm{m}-1}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_{\mathrm{j}}+{\mathrm{e}}_{\mathrm{q}}{\mathrm{\lambda }}_{\mathrm{m}}$$

for some ${\mathrm{\lambda }}_{\mathrm{j}}\in \mathbb{R}$, $\mathrm{j}=1,...,\mathrm{m}$ and the following conditions holds:

$$\sum _{\mathrm{j}=1}^{\mathrm{m}-1}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_{\mathrm{j}}·{\mathrm{n}}_{\mathrm{k}}+{\mathrm{\lambda }}_{\mathrm{m}}{\mathrm{n}}_{\mathrm{q}}^{\mathrm{k}}=0$$

for each $\mathrm{k}=1,...,\mathrm{m}-1$ and

$${\parallel \mathrm{p}\parallel }^2={\mathrm{r}}^2$$

Thus, if $\mathrm{p}={\sum }_{\mathrm{j}=1}^{\mathrm{m}-1}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_{\mathrm{j}}+{\mathrm{e}}_{\mathrm{q}}{\mathrm{\lambda }}_{\mathrm{m}}$ is an ${\mathrm{e}}_{\mathrm{q}}$-pole of $\mathrm{S}-\left\{{\mathrm{c}}_{\mathrm{m}}\right\}$, the following equations are satisfied for ${\mathrm{\lambda }}_1,{\mathrm{\lambda }}_2,...,{\mathrm{\lambda }}_{\mathrm{m}}\in \mathbb{R}$:

$$\sum _{\mathrm{j}=1}^{\mathrm{m}-1}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_{\mathrm{j}}·{\mathrm{n}}_{\mathrm{k}}+{\mathrm{\lambda }}_{\mathrm{m}}{\mathrm{n}}_{\mathrm{q}}^{\mathrm{k}}=0$$

(4)

$$\sum _{\mathrm{i}=1}^{\mathrm{d}}{\left(\sum _{\mathrm{j}=1}^{\mathrm{m}-1}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_{\mathrm{i}}^{\mathrm{j}}\right)}^2+2{\mathrm{\lambda }}_{\mathrm{m}}\sum _{\mathrm{j}=1}^{\mathrm{m}-1}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_{\mathrm{q}}^{\mathrm{j}}+{\mathrm{\lambda }}_{\mathrm{m}}^2={\mathrm{r}}^2$$

(5)

for all $\mathrm{k}=1,...,\mathrm{m}-1$, with r the radius of the $(\mathrm{d}-\mathrm{m}+1)$-sphere S. From (4) we have the system

$$\left(\begin{array}{cccc}{\mathrm{n}}_1·{\mathrm{n}}_1& {\mathrm{n}}_1·{\mathrm{n}}_2& ....& {\mathrm{n}}_1·{\mathrm{n}}_{\mathrm{m}-1}\\ {\mathrm{n}}_2·{\mathrm{n}}_1& {\mathrm{n}}_2·{\mathrm{n}}_2& ....& {\mathrm{n}}_2·{\mathrm{n}}_{\mathrm{m}-1}\\ ...\\ {\mathrm{n}}_{\mathrm{m}-1}·{\mathrm{n}}_1& {\mathrm{n}}_{\mathrm{m}-1}·{\mathrm{n}}_2& ....& {\mathrm{n}}_{\mathrm{m}-1}·{\mathrm{n}}_{\mathrm{m}-1}\end{array}\right)\left(\begin{array}{c}{\mathrm{\lambda }}_1\\ {\mathrm{\lambda }}_2\\ ...\\ {\mathrm{\lambda }}_{\mathrm{m}-1}\end{array}\right)+{\mathrm{\lambda }}_{\mathrm{m}}\left(\begin{array}{c}{\mathrm{n}}_{\mathrm{q}}^1\\ {\mathrm{n}}_{\mathrm{q}}^2\\ ...\\ {\mathrm{n}}_{\mathrm{q}}^{\mathrm{m}-1}\end{array}\right)=0.$$

Let us denote A as the matrix ${({\mathrm{n}}_{\mathrm{i}}·{\mathrm{n}}_{\mathrm{j}})}_{\mathrm{i},\mathrm{j}}$ and B as the vector ${(-{\mathrm{n}}_{\mathrm{q}}^{\mathrm{j}})}_{\mathrm{j}=1}^{\mathrm{m}-1}$. Then, we have $\mathrm{A}\mathrm{\lambda }={\mathrm{\lambda }}_{\mathrm{m}}\mathrm{B}$, where $\mathrm{\lambda }=({\mathrm{\lambda }}_1,...,{\mathrm{\lambda }}_{\mathrm{m}-1})$. Solving for $\mathrm{\lambda }$, we obtain ${\mathrm{\lambda }}_{\mathrm{j}}={\mathrm{\lambda }}_{\mathrm{m}}({\mathrm{A}}^{-1}\mathrm{B})\left[\mathrm{j}\right]$ for each $\mathrm{j}=1,...,\mathrm{m}-1$ (where $({\mathrm{A}}^{-1}\mathrm{B})\left[\mathrm{j}\right]$ denotes the entry j of the $(\mathrm{m}-1)\times 1$ vector $({\mathrm{A}}^{-1}\mathrm{B})$). By substituting the value of ${\mathrm{\lambda }}_{\mathrm{j}}$ into (5), we obtain the quadratic equation:

$${\mathrm{\lambda }}_{\mathrm{m}}^2\left[\sum _{\mathrm{i}=1}^{\mathrm{d}}{\left(\sum _{\mathrm{j}=1}^{\mathrm{m}-1}({\mathrm{A}}^{-1}\mathrm{B})\left[\mathrm{j}\right]{\mathrm{n}}_{\mathrm{i}}^{\mathrm{j}}\right)}^2+2\sum _{\mathrm{j}=1}^{\mathrm{m}-1}({\mathrm{A}}^{-1}\mathrm{B})\left[\mathrm{j}\right]{\mathrm{n}}_{\mathrm{q}}^{\mathrm{j}}+1\right]-{\mathrm{r}}^2=0$$

Let us define

$${\mathrm{\Gamma }}_{\mathrm{i}}=\sum _{\mathrm{j}=1}^{\mathrm{m}-1}({\mathrm{A}}^{-1}\mathrm{B})\left[\mathrm{j}\right]{\mathrm{n}}_{\mathrm{i}}^{\mathrm{j}}$$

for all $\mathrm{i}=1,...,\mathrm{d}$. Solving this equation, we find:

$${\mathrm{\lambda }}_{\mathrm{m}}=\frac{\pm \mathrm{r}}{\sqrt{{\sum }_{\mathrm{i}=1}^{\mathrm{d}}{\mathrm{\Gamma }}_{\mathrm{i}}^2+2{\mathrm{\Gamma }}_{\mathrm{q}}+1}}$$

$${\mathrm{\lambda }}_{\mathrm{j}}=\frac{\pm \mathrm{r}({\mathrm{A}}^{-1}\mathrm{B})\left[\mathrm{j}\right]}{\sqrt{{\sum }_{\mathrm{i}=1}^{\mathrm{d}}{\mathrm{\Gamma }}_{\mathrm{i}}^2+2{\mathrm{\Gamma }}_{\mathrm{q}}+1}}$$

for $\mathrm{j}=1,...,\mathrm{m}-1$. Therefore, the ${\mathrm{e}}_{\mathrm{q}}$-poles of S, for $\mathrm{q}=1,...,\mathrm{d}$, are:

$$\mathrm{p}=\sum _{\mathrm{j}=1}^{\mathrm{m}-1}{\mathrm{\lambda }}_{\mathrm{j}}{\mathrm{n}}_{\mathrm{j}}+{\mathrm{e}}_{\mathrm{q}}{\mathrm{\lambda }}_{\mathrm{m}}+\mathrm{c}=\frac{\pm \mathrm{r}}{\sqrt{{\sum }_{\mathrm{i}=1}^{\mathrm{d}}{\mathrm{\Gamma }}_{\mathrm{i}}^2+2{\mathrm{\Gamma }}_{\mathrm{q}}+1}}\left(\sum _{\mathrm{j}=1}^{\mathrm{m}-1}({\mathrm{A}}^{-1}\mathrm{B})\left[\mathrm{j}\right]{\mathrm{n}}_{\mathrm{j}}+{\mathrm{e}}_{\mathrm{q}}\right)+\mathrm{c}$$

3. Vietoris-Rips and Čech systems

Our goal in this section is to provide a comprehensive understanding of the disk system, the Vietoris-Rips system, and the Čech system. Additionally, we introduce some results that establish a certain connection between both disk systems. Investigating the features and qualities of data and spaces can provide us with useful knowledge about their geometric and topological characteristics.

Before we look into the definitions of Vietoris-Rips and Čech systems, let us give a brief overview. These systems are essential in the field of topological data analysis for recognizing and comprehending the geometric structure of point cloud data. Vietoris-Rips complex and Čech complex share the goal of capturing the topology of the underlying metric space, both provide different ways of recognizing connections and associations among data points. The Vietoris-Rips complex tends to be more efficient and scalable for large datasets, while the Čech complex can be more accurate but computationally more expensive. The choice between the two depends on the nature of the dataset and the specific goals of the topological analysis. Now, let us move on to defining these fundamental concepts.

Definition 3.

Let $M=\{{\mathrm{D}}_1,{\mathrm{D}}_2,\dots ,{\mathrm{D}}_{\mathrm{m}}\}$ be a d-disk system. We say M is a Vietoris-Rips system if ${\mathrm{D}}_{\mathrm{i}}\cap {\mathrm{D}}_{\mathrm{j}}\ne \varnothing$ for each pair $\mathrm{i},\mathrm{j}\in \{1,2,\dots ,\mathrm{m}\}$. Furthermore, if the d-disk system M has the nonempty intersection property ${\bigcap }_{{\mathrm{D}}_{\mathrm{i}}\in M}{\mathrm{D}}_{\mathrm{i}}\ne \varnothing$, then M is called a Čech system.

For each $\mathrm{\lambda }\ge 0$, we define a collection of d-disks $M_{\mathrm{\lambda }}$ with the same centers as those in the d-disk system M, but with radii rescaled by $\mathrm{\lambda }$. When $\mathrm{\lambda }>0$, $M_{\mathrm{\lambda }}$ is a d-disk system again. $M_1$ is equal to M, and $M_0$ is the set of the centers of the d-disks in M.

In the field of topological data analysis, the Rips scale and the Čech scale are essential parameters for determining the closeness and connectivity between data points. These two scales offer different perspectives on how we measure and comprehend geometric relationships within point-cloud data. To understand their importance in capturing the underlying topological structure, let us look at their definitions. The Vietoris-Rips scale ${\mathrm{\nu }}_M$, of a d-disk system M is the smallest $\mathrm{\lambda }\in \mathbb{R}$ such that $M_{\mathrm{\lambda }}$ is a Vietoris-Rips system. Similarly, the Čech scale ${\mathrm{\mu }}_M$, of M is the smallest $\mathrm{\lambda }\in \mathbb{R}$ such that $M_{\mathrm{\lambda }}$ is a Čech system. This is

$$\begin{array}{cc}\hfill {\mathrm{\nu }}_M& =inf\{\mathrm{\lambda }\in \mathbb{R}\mid M_{\mathrm{\lambda }}\mathrm{i}\mathrm{s}\mathrm{a}\mathrm{V}\mathrm{i}\mathrm{e}\mathrm{t}\mathrm{o}\mathrm{r}\mathrm{i}\mathrm{s}-\mathrm{R}\mathrm{i}\mathrm{p}\mathrm{s}\mathrm{s}\mathrm{y}\mathrm{s}\mathrm{t}\mathrm{e}\mathrm{m}\}\hfill \\ \hfill {\mathrm{\mu }}_M& =inf\{\mathrm{\lambda }\in \mathbb{R}\mid M_{\mathrm{\lambda }}\mathrm{i}\mathrm{s}\mathrm{a}\mathrm{Č}\mathrm{e}\mathrm{c}\mathrm{h}\mathrm{s}\mathrm{y}\mathrm{s}\mathrm{t}\mathrm{e}\mathrm{m}\}\hfill \end{array}$$

Next, we present some easily observable properties for both scales. It can be easily seen that M is a Vietoris-Rips system if and only if ${\mathrm{\nu }}_M\le 1$ (in particular ${\mathrm{\nu }}_{M_{{\mathrm{\nu }}_M}}=1$); similarly, M is a Čech system if and only if ${\mathrm{\mu }}_M\le 1$.

Note that for a given d-disk system $M=\{{\mathrm{D}}_1,{\mathrm{D}}_2,\dots ,{\mathrm{D}}_{\mathrm{m}}\}$ the Vietoris-Rips scale is

$${\mathrm{\nu }}_M=\underset{\mathrm{i}<\mathrm{j}}{max}\{\parallel {\mathrm{c}}_{\mathrm{i}}-{\mathrm{c}}_{\mathrm{j}}\parallel /({\mathrm{r}}_{\mathrm{i}}+{\mathrm{r}}_{\mathrm{j}})\}$$

where ${\mathrm{c}}_{\mathrm{i}}$ and ${\mathrm{r}}_{\mathrm{i}}$ are the center and radii of ${\mathrm{D}}_{\mathrm{i}}$. An additional observation is that, in cases where the disk system consists of either one or two disks, the Vietoris-Rips scale coincides with the Čech scale. It is evident that every Čech system is also a Vietoris-Rips system; however, the reverse assertion, in general, is not true.

Conversely, if the d-disk system contains at least three disks, determining the Čech scale becomes more complex. In the context of Čech scale, the following remark is important and play a key role in implementation (see [3] for details).

Remark 1.

If ${\mathrm{\mu }}_M$ is the Čech scale for M, then the ${\mathrm{\mu }}_M$-rescaled system $M_{{\mathrm{\mu }}_M}$, has only one point in the intersection ${\bigcap }_{{\mathrm{D}}_{\mathrm{i}}\in M}{\mathrm{D}}_{\mathrm{i}}({\mathrm{c}}_{\mathrm{i}};{\mathrm{\mu }}_M{\mathrm{r}}_{\mathrm{i}})$.

As we have mentioned, a Čech system is also a Vietoris-Rips system, but the converse is not true. What we can affirm is that if a system is a Vietoris-Rips system, then the system rescaled by the factor $\sqrt{2\mathrm{d}/(\mathrm{d}+1)}$ is also a Čech system. This is established by the following lemma, the proof of which can be found in [3].

Lemma 5.

Let $M=\left\{{\mathrm{D}}_{\mathrm{i}}({\mathrm{c}}_{\mathrm{i}};{\mathrm{r}}_{\mathrm{i}})\right\}$ be a d-disk system in euclidean space ${\mathbb{R}}^{\mathrm{d}}$. If ${\mathrm{D}}_{\mathrm{i}}({\mathrm{c}}_{\mathrm{i}};{\mathrm{r}}_{\mathrm{i}})\cap {\mathrm{D}}_{\mathrm{j}}({\mathrm{c}}_{\mathrm{j}};{\mathrm{r}}_{\mathrm{j}})\ne \varnothing$ for every pair of disks in M, then

$$\bigcap _{{\mathrm{D}}_{\mathrm{i}}\in M}{\mathrm{D}}_{\mathrm{i}}({\mathrm{c}}_{\mathrm{i}};\sqrt{2\mathrm{d}/(\mathrm{d}+1)}{\mathrm{r}}_{\mathrm{i}})\ne \varnothing .$$

One of the implications of the previous result is that, for any given disk system M, we can bound the Čech scale using the Vietoris-Rips scale ${\mathrm{\nu }}_M$. This is stated by the following corollary.

Corollary 1.

If M is an arbitrary d-disk system and ${\mathrm{\nu }}_M$ is its Vietoris-Rips scale, then its Čech scale satisfies ${\mathrm{\nu }}_M\le {\mathrm{\mu }}_M\le \sqrt{2\mathrm{d}/(\mathrm{d}+1)}{\mathrm{\nu }}_M$. Therefore, for every d-disk system M, the rescaled disk system $M_{\sqrt{2\mathrm{d}/(\mathrm{d}+1)}{\mathrm{\nu }}_M}$ is always a Čech system. In particular, if $\sqrt{2\mathrm{d}/(\mathrm{d}+1)}{\mathrm{\nu }}_M\le 1$ then $M_{{\mathrm{\nu }}_M}$ is a Čech system.

3.1. Algorithm for determining Čech system.

In the previous section, we have determined the ${\mathrm{e}}_{\mathrm{q}}$-poles for the intersection of any number of disks in ${\mathbb{R}}^{\mathrm{d}}$. If any of these poles is in all disks of the disk system, it indicates that the system conforms to the criteria of a Čech system. It is important to recognize that this result streamlines our calculation process, focusing on specific points to establish whether the system exhibits a non-empty intersection.

Given a system of m disks in ${\mathbb{R}}^{\mathrm{d}}$ where $\mathrm{m}>\mathrm{d}$, it is enough to verify if every subsystem of $\mathrm{d}+1$ disks qualifies as a Čech system to conclude that the entire system of disks has a non-empty intersection. This assertion is supported by the Helly’s Theorem.

Now, we introduce an algorithm that determines whether a disk system qualifies as a Čech system. In simpler terms, if the disk system exhibits a non-empty intersection, the algorithm outputs "TRUE"; otherwise, it outputs "FALSE". The algorithm operates by seeking poles within the intersections of the disk boundaries, which, as we have observed, correspond to i-spheres. It initiates the search for poles within individual disks and then progresses to the pairwise intersections of the disk boundaries ($(\mathrm{d}-2)-\mathrm{s}\mathrm{p}\mathrm{h}\mathrm{e}\mathrm{r}\mathrm{e}$), continuing the process iteratively. If a pole is found within the remaining disks, the system is classified as a Čech system.

Theorem 1.

Let $M=\{{\mathrm{D}}_1,\dots ,{\mathrm{D}}_{\mathrm{m}}\}$ be a d-disk system. Then, M is a Čech system if and only if Cech.system$(M)=$ TRUE.

Proof. 

If Cech.system$(M)=$ TRUE, then the Cech.system algorithm (Algorithm 1) found a pole contained in the intersection ${\cap }_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{D}}_{\mathrm{j}}$, therefore ${\cap }_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{D}}_{\mathrm{j}}\ne \varnothing$ and it follows that M is a Čech system.

On the other hand, if ${\bigcap }_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{D}}_{\mathrm{j}}\ne \varnothing$, let p be a point in ${\bigcap }_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{D}}_{\mathrm{j}}$ satisfying ${\mathrm{\pi }}_1(\mathrm{p})\le {\mathrm{\pi }}_1(\mathrm{x})$ for every x in ${\bigcap }_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{D}}_{\mathrm{j}}$. By Lemma 2, it follows that p belongs to an i-sphere and must be an ${\mathrm{e}}_1$-south pole. Therefore, by the exhaustive search of Algorithm 1 across all poles, its output is Cech.system$(M)=$ TRUE.    □

Algorithm 1: Cech.system

3.2. Algorithm to compute the Čech scale

Finding the minimum parameter for which the rescaled system of disks has a non-empty intersection is significant because it helps identify a critical threshold at which the disks come into contact. This parameter, known as the Čech scale, provides valuable information about the proximity or overlap of the disks, which can be crucial in various applications such as collision detection in computer graphics, spatial packing problems, and modeling physical phenomena. In this section, we introduce an algorithm to compute an approximation of the Čech scale for a system of m disks in ${\mathbb{R}}^{\mathrm{d}}$.  

Algorithm 2: Cech.scale

The given code presents an algorithm to compute an approximation of the Čech scale of a disk system in Euclidean space using Algorithm 1 and a precision parameter $\mathrm{\eta }>0$. It initializes the scale factor $\mathrm{\lambda }$ to the Rips scale ${\mathrm{\nu }}_M$. If this scale satisfies Cech.system$(M_{{\mathrm{\nu }}_M})=$ TRUE, it indicates that the Čech scale has been found. Otherwise, we initiate a cycle in which we compute Cech.system of the system rescaled by a factor $\mathrm{\lambda }$. The Čech scale is known to fall between the Rips scale and the value $\sqrt{2\mathrm{d}/(\mathrm{d}+1)}{\mathrm{\nu }}_M$ (Generalized Vietoris-Rips, Corollary 1). To approximate the Čech scale, we employ the bisection method as long as the interval enclosing the Čech scale has a length greater than $\mathrm{\eta }$. Finally, the algorithm returns an approximation of the Čech scale.

Utilizing the previously described algorithm, we can construct the filtered generalized Čech complex for a disk system M. Let $\mathcal{C}(M)$ denote the set of Čech subsystems, and ${\mathcal{C}}_{\mathrm{\lambda }}(M)$ the set of Čech subsystems for the rescaled disk system $M_{\mathrm{\lambda }}$. The Čech filtration of the M system forms a maximal chain of Čech complexes

$${\mathcal{C}}_{*}(M):{\mathcal{C}}_0(M)⊊{\mathcal{C}}_{{\mathrm{\lambda }}_1}(M)⊊{\mathcal{C}}_{{\mathrm{\lambda }}_2}(M)⊊...⊊{\mathcal{C}}_{{\mathrm{\mu }}_M}(M),$$

where each ${\mathrm{\lambda }}_{\mathrm{i}}$ represents the Čech scale of the system $M_{{\mathrm{\lambda }}_{\mathrm{i}}}$. Since the Čech scale of a disk system indicates the factor by which we must rescale the system to make it Čech, defining a level of the filtration, ${\mathcal{C}}_{\mathrm{\lambda }}(M)$, simply requires determining the Čech scale of the system $M_{\mathrm{\lambda }}$.

4. Minimal Axis-Aligned Bounding Box.

In this section, we introduce the concept of the minimal axis-aligned bounding box (AABB) for the intersection of d-disks and present methods for its computation. The AABB provides a simplified representation of the disk intersection, making it easier to obtain valuable information about the disks. This information could be useful for computing the Čech scale of a disk system.

Definition 4.

Let M be a disk system in ${\mathbb{R}}^{\mathrm{d}}$. The minimal axis-aligned bounding box of M, denoted as $\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M)$, is defined as the smallest axis-aligned bounding box that contains the intersection $\mathrm{D}={\cap }_{{\mathrm{D}}_{\mathrm{i}}\in M}{\mathrm{D}}_{\mathrm{i}}$, given by

$$\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M):=\bigcap _{\mathrm{D}\subset \tilde{\mathrm{B}}}\tilde{\mathrm{B}}$$

where $\tilde{\mathrm{B}}$ ranges over all axis-aligned bounding boxes that contain D.

Note that the AABB can be expressed as

$$\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M)=\prod _{\mathrm{k}=1}^{\mathrm{d}}[inf{\mathrm{\pi }}_{\mathrm{k}}(\partial \mathrm{D}),sup{\mathrm{\pi }}_{\mathrm{k}}(\partial \mathrm{D})]$$

where ${\mathrm{\pi }}_{\mathrm{k}}:{\mathbb{R}}^{\mathrm{d}}⟶\mathbb{R}$ is the canonical projection onto the k-th factor, and $\partial \mathrm{D}$ denotes the boundary of D. In other words, the AABB of a disk system M is given by the Cartesian product of intervals, where each interval is determined by the minimum and maximum values of the corresponding projection of the disk boundaries.

4.1. Minimal axis-aligned bounding box for two disks

Let’s consider the situation when the disk system M is composed of two disks ${\mathrm{D}}_1$ and ${\mathrm{D}}_2$ in ${\mathbb{R}}^{\mathrm{d}}$. If ${\mathrm{D}}_1\cap {\mathrm{D}}_2\ne \varnothing$ and ${\mathrm{D}}_1\ne {\mathrm{D}}_2$, the subset $\partial {\mathrm{D}}_1\cap \partial {\mathrm{D}}_2$ can take one of three forms: an empty set, a single common point (when the disks are tangent), or a $(\mathrm{d}-2)$-dimensional sphere. In the last case, we denote the $(\mathrm{d}-2)$-dimensional sphere or the $(\mathrm{d}-1)$-sphere $\partial {\mathrm{D}}_1\cap \partial {\mathrm{D}}_2$ by ${\mathrm{S}}_{1,2}$. To calculate $inf{\mathrm{\pi }}_{\mathrm{i}}(\partial ({\mathrm{D}}_1\cap {\mathrm{D}}_2))$ and $sup{\mathrm{\pi }}_{\mathrm{i}}(\partial ({\mathrm{D}}_1\cap {\mathrm{D}}_2))$ for each $\mathrm{i}\in \{1,2,...,\mathrm{d}\}$, we can use:

$$\begin{array}{c}\hfill inf{\mathrm{\pi }}_{\mathrm{i}}(\partial ({\mathrm{D}}_1\cap {\mathrm{D}}_2))=\left\{\begin{array}{cc}{\mathrm{\pi }}_{\mathrm{i}}({\mathrm{c}}_1-{\mathrm{r}}_1{\mathrm{e}}_{\mathrm{i}})\hfill & \mathrm{if}{\mathrm{c}}_1-{\mathrm{r}}_1{\mathrm{e}}_{\mathrm{i}}\in {\mathrm{D}}_2,\hfill \\ {\mathrm{\pi }}_{\mathrm{i}}({\mathrm{c}}_2-{\mathrm{r}}_2{\mathrm{e}}_{\mathrm{i}})\hfill & \mathrm{if}{\mathrm{c}}_2-{\mathrm{r}}_2{\mathrm{e}}_{\mathrm{i}}\in {\mathrm{D}}_1,\hfill \\ inf({\mathrm{\pi }}_{\mathrm{i}}(\partial {\mathrm{D}}_1\cap \partial {\mathrm{D}}_2))\hfill & \mathrm{otherwise},\hfill \end{array}\right.\end{array}$$

(6)

$$\begin{array}{c}\hfill sup{\mathrm{\pi }}_{\mathrm{i}}(\partial ({\mathrm{D}}_1\cap {\mathrm{D}}_2))=\left\{\begin{array}{cc}{\mathrm{\pi }}_{\mathrm{i}}({\mathrm{c}}_1+{\mathrm{r}}_1{\mathrm{e}}_{\mathrm{i}})\hfill & \mathrm{if}{\mathrm{c}}_1+{\mathrm{r}}_1{\mathrm{e}}_{\mathrm{i}}\in {\mathrm{D}}_2,\hfill \\ {\mathrm{\pi }}_{\mathrm{i}}({\mathrm{c}}_2+{\mathrm{r}}_2{\mathrm{e}}_{\mathrm{i}})\hfill & \mathrm{if}{\mathrm{c}}_2+{\mathrm{r}}_2{\mathrm{e}}_{\mathrm{i}}\in {\mathrm{D}}_1,\hfill \\ sup({\mathrm{\pi }}_{\mathrm{i}}(\partial {\mathrm{D}}_1\cap \partial {\mathrm{D}}_2))\hfill & \mathrm{otherwise}.\hfill \end{array}\right.\end{array}$$

Indeed, by Lemma 2, the extremes of the AABB are the projections of certain poles, either from the $(\mathrm{d}-1)$-sphere or some d-sphere. The d-spheres represent the boundaries of each disk, with poles given by ${\mathrm{c}}_{\mathrm{j}}\pm {\mathrm{r}}_{\mathrm{j}}{\mathrm{e}}_{\mathrm{i}}$, and the $(\mathrm{d}-1)$-sphere is ${\mathrm{S}}_{1,2}=\partial {\mathrm{D}}_1\cap \partial {\mathrm{D}}_2$, whose poles are computed using Lemma 4. It is worth noting that there are no further options for $(\mathrm{d}-\mathrm{m}+1)$-spheres in the case of a two-disk system.

To simplify the notation, we will use ${\mathrm{B}}_{\mathrm{i},\mathrm{j}}$ to denote the AABB of the intersection of disks ${\mathrm{D}}_{\mathrm{i}}$ and ${\mathrm{D}}_{\mathrm{j}}$. Figure 2 illustrates the AABB of the intersection of two disks in the plane.

According to (6) and Lemma 3, we have a method to calculate the axis-aligned boundary box (AABB) for systems of two disks.

Knowing how to compute the AABB of two disks is not sufficient to determine the AABB for a disk system with more than two disks in ${\mathbb{R}}^{\mathrm{d}}$. In the following examples, we demonstrate that the AABB of a disk system is not simply the intersection of all AABBs of two disks in ${\mathbb{R}}^3$.

Example 1. Let

$$M=\{{\mathrm{D}}_1((4,1,0);\sqrt{2}),{\mathrm{D}}_2((4,-1,0);\sqrt{2}),{\mathrm{D}}_3((0,0,0);3)\}\subset {\mathbb{R}}^3$$

be a Vietoris-Rips system in ${\mathbb{R}}^3$ with the following projection onto the $\mathrm{x}\mathrm{y}$-plane:

By computing the boxes ${\mathrm{B}}_{\mathrm{i},\mathrm{j}}$ for all $\mathrm{i}<\mathrm{j}$, we obtain:

$$\begin{array}{cc}\hfill {\mathrm{B}}_{1,2}& =[3,5]\times [-\sqrt{2},\sqrt{2}]\times [-1,1]\hfill \\ \hfill {\mathrm{B}}_{1,3}& =[4-\sqrt{2},3]\times [0,1.41]\times [-0.72,0.72]\hfill \\ \hfill {\mathrm{B}}_{2,3}& =[4-\sqrt{2},3]\times [-1.4,0]\times [-0.72,0.72]\hfill \end{array}$$

Therefore, ${\cap }_{\mathrm{i}<\mathrm{j}}{\mathrm{B}}_{\mathrm{i},\mathrm{j}}=[3,3]\times [0,0]\times [-0.72,0.72]$. However, the disk system intersects at the point $\mathrm{P}=(3,0,0)$, which means $\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M)=\left\{\mathrm{P}\right\}$. In other words, $\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M)$ is not equal to ${\cap }_{\mathrm{i}<\mathrm{j}}{\mathrm{B}}_{\mathrm{i},\mathrm{j}}$.

We know that if $\mathrm{D}\ne \varnothing$, then ${\cap }_{\mathrm{i}<\mathrm{j}}{\mathrm{B}}_{\mathrm{i},\mathrm{j}}\ne \varnothing$. However, the converse is not always true. The following example illustrates this fact.

Example 2. Let $\mathrm{N}=\{{\mathrm{D}}_1,{\mathrm{D}}_2,{\mathrm{D}}_3((0,1,0);\sqrt{10}),{\mathrm{D}}_4((3,0,1);0.9)\}\subset {\mathbb{R}}^3$ be a Vietoris-Rips system. The projection of disks ${\mathrm{D}}_1,{\mathrm{D}}_2,$ and ${\mathrm{D}}_3$ onto the $\mathrm{x}\mathrm{y}$-plane is illustrated in Figure 4.

We will now compute the intersections ${\mathrm{B}}_{\mathrm{i},\mathrm{j}}$ for different pairs of disks:

$$\begin{array}{cc}\hfill {\mathrm{B}}_{1,2}& =[3,5]\times [-\sqrt{2},\sqrt{2}]\times [-1,1]\hfill \\ \hfill {\mathrm{B}}_{1,3}& =[4-\sqrt{2},\sqrt{10}]\times [0,2]\times [-1,1]\hfill \\ \hfill {\mathrm{B}}_{1,4}& =[2.6,3.9]\times [-0.4,0.9]\times [0.100007,1.3]\hfill \\ \hfill {\mathrm{B}}_{2,3}& =[2.6,3]\times [-0.8,0]\times [-0.44,4.47]\hfill \\ \hfill {\mathrm{B}}_{2,4}& =[2.6,3.9]\times [-0.9,0.4]\times [0.100007,1.3]\hfill \\ \hfill {\mathrm{B}}_{3,4}& =[2.1,3.11]\times [-0.73,0.9]\times [0.1,1.7]\hfill \end{array}$$

The intersection of all pairwise intersections, ${\cap }_{\mathrm{i}<\mathrm{j}}{\mathrm{B}}_{\mathrm{i},\mathrm{j}}$, is given by $[3,3]\times [0,0]\times [0.1,1.7]$. However, the intersection of disks ${\mathrm{D}}_1,{\mathrm{D}}_2,$ and ${\mathrm{D}}_3$ is a single point P, which is not contained in ${\mathrm{D}}_4$ (by construction). Therefore, D is an empty set, but ${\cap }_{\mathrm{i}<\mathrm{j}}{\mathrm{B}}_{\mathrm{i},\mathrm{j}}$ is not.

These examples clearly illustrate that when dealing with the AABB of three disks or more, knowing the AABB for pairs of disks is insufficient. In Example Section 4.1, we observe that the intersection of $\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(\{{\mathrm{D}}_{\mathrm{i}},{\mathrm{D}}_{\mathrm{j}}\})$ is not equal to the AABB of the disk system M. Similarly, in Example Section 4.1, we find that the intersection of three disks is empty, yet the intersection of $\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(\{{\mathrm{D}}_{\mathrm{i}},{\mathrm{D}}_{\mathrm{j}}\})$ contains points.

Therefore, the next crucial step is to determine how to calculate the AABB of a disk system consisting of more than two disks in ${\mathbb{R}}^{\mathrm{d}}$.

4.2. Minimal axis-aligned bounding box for more than two disks

Given a system M consisting of m disks in ${\mathbb{R}}^{\mathrm{d}}$, we can compute the ${\mathrm{e}}_{\mathrm{q}}$-poles for any subcollection of disks. Using these ${\mathrm{e}}_{\mathrm{q}}$-poles, we can determine the axis-aligned bounding box (AABB) of M.

If $\mathrm{m}\le \mathrm{d}$, we calculate the ${\mathrm{e}}_{\mathrm{q}}$-poles of the $(\mathrm{d}-\mathrm{m}+1)$-sphere $\cap \partial {\mathrm{D}}_{\mathrm{i}}$, and with these poles, we define the AABB of M by taking $inf{\mathrm{\pi }}_{\mathrm{q}}\partial \mathrm{D}={\mathrm{\pi }}_{\mathrm{q}}(\mathrm{p})$, where p is the ${\mathrm{e}}_{\mathrm{q}}$-south pole of the $(\mathrm{d}-\mathrm{m}+1)$-sphere (similarly for $sup{\mathrm{\pi }}_{\mathrm{q}}\partial \mathrm{D}$).

In the case where $\mathrm{m}=\mathrm{d}+1$, we consider $\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M(1)),...,\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M(\mathrm{m}))$ as a collection of minimal axis-aligned boxes for the disk system $M(\mathrm{i})=M-\left\{{\mathrm{D}}_{\mathrm{i}}\right\}$ in ${\mathbb{R}}^{\mathrm{d}}$. If the intersection of any $\mathrm{d}+1$ of these sets is nonempty, then the intersection of the entire collection gives us the minimal axis-aligned box of the disk system M. This can be expressed as $\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M)=\cap \mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M(\mathrm{i}))$. The next theorem confirms this finding.

Theorem 2

(Helly’s theorem for minimal axis-aligned boxes). Let M be a Rips system with $\mathrm{d}+1$ disks in ${\mathbb{R}}^{\mathrm{d}}$. Then, the minimal axis-aligned bounding box for the intersection set $\mathrm{D}={\cap }_{\mathrm{j}=1}^{\mathrm{d}+1}{\mathrm{D}}_{\mathrm{j}}$ satisfies:

$$\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M)=\bigcap _{\mathrm{j}=1}^{\mathrm{d}+1}\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M(\mathrm{j}))$$

where $M(\mathrm{j})=M-\left\{{\mathrm{D}}_{\mathrm{j}}\right\}$.

Proof. 

We know that $\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M)\subseteq {\bigcap }_{\mathrm{j}=1}^{\mathrm{d}+1}\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M(\mathrm{j}))$. Now, our objective is to establish the reverse inclusion, that is, ${\bigcap }_{\mathrm{j}=1}^{\mathrm{d}+1}\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M(\mathrm{j}))\subseteq \mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M)$. In order to derive a contradiction, suppose that the reverse inclusion is not true. By definition, we have:

$$\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M)=\prod _{\mathrm{i}=1}^{\mathrm{d}}\left[inf{\mathrm{\pi }}_{\mathrm{i}}(\partial \mathrm{D}),sup{\mathrm{\pi }}_{\mathrm{i}}(\partial \mathrm{D})\right]$$

and

$$\bigcap _{\mathrm{j}=1}^{\mathrm{d}+1}\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M(\mathrm{j}))=\prod _{\mathrm{i}=1}^{\mathrm{d}}\left[\underset{\mathrm{k}}{max}\{inf{\mathrm{\pi }}_{\mathrm{i}}(\partial {\cap }_{\mathrm{j}\ne \mathrm{k}}{\mathrm{D}}_{\mathrm{j}})\},\underset{\mathrm{k}}{min}\{sup{\mathrm{\pi }}_{\mathrm{i}}(\partial {\cap }_{\mathrm{j}\ne \mathrm{k}}{\mathrm{D}}_{\mathrm{j}})\}\right]$$

Without loss of generality, let’s assume that:

$${\mathrm{\pi }}_1(\mathrm{p})>\underset{\mathrm{k}}{max}{\left\{inf{\mathrm{\pi }}_1({\cap }_{\mathrm{j}\ne \mathrm{k}}\partial {\mathrm{D}}_{\mathrm{j}})\right\}}_{\mathrm{k}=1}^{\mathrm{d}+1}$$

where $\mathrm{p}\in \mathrm{D}$ satisfies ${\mathrm{\pi }}_1(\mathrm{p})=inf{\mathrm{\pi }}_1(\partial \mathrm{D})$.

Let ${\mathrm{q}}_{\mathrm{j}}$ be the point in ${\cap }_{\mathrm{k}\ne \mathrm{j}}{\mathrm{D}}_{\mathrm{k}}$ such that ${\mathrm{\pi }}_1({\mathrm{q}}_{\mathrm{j}})=inf{\mathrm{\pi }}_1({\cap }_{\mathrm{k}\ne \mathrm{j}}{\mathrm{D}}_{\mathrm{k}})$ for each $\mathrm{j}=1,2,...,\mathrm{d}+1$. Note that ${\mathrm{q}}_{\mathrm{j}}\notin {\mathrm{D}}_{\mathrm{j}}$ because ${\mathrm{q}}_{\mathrm{j}}$ is not in D. Now, let ${\mathrm{\gamma }}_{\mathrm{j}}$ be the line segment that connects ${\mathrm{q}}_{\mathrm{j}}$ and p.

Choose $\mathrm{ϵ}>0$ small enough such that the hyperplane $\mathrm{P}:{\mathrm{x}}_1={\mathrm{\pi }}_1(\mathrm{p})-\mathrm{ϵ}$ does not contain any ${\mathrm{q}}_{\mathrm{j}}$, and ${\mathrm{\pi }}_1({\mathrm{q}}_{\mathrm{j}})<{\mathrm{\pi }}_1(\mathrm{p})-\mathrm{ϵ}$ for all $\mathrm{j}=1,2,...,\mathrm{d}+1$ (such hyperplane P exists because ${\mathrm{\pi }}_1(\mathrm{p})>{\mathrm{\pi }}_1({\mathrm{q}}_{\mathrm{j}})$ for all $\mathrm{j}=1,...,\mathrm{d}+1$). Since ${\mathrm{q}}_{\mathrm{j}}$ is in ${\cap }_{\mathrm{k}\ne \mathrm{j}}{\mathrm{D}}_{\mathrm{k}}$, the hyperplane P intersects every disk in $M(\mathrm{j})$ for each j, and ${\mathrm{\gamma }}_{\mathrm{j}}\subset {\cap }_{\mathrm{k}\ne \mathrm{j}}{\mathrm{D}}_{\mathrm{k}}$ intersects P at a point that is in ${\cap }_{\mathrm{k}\ne \mathrm{j}}{\mathrm{D}}_{\mathrm{k}}$ (see Figure 5). Furthermore, ${\mathrm{D}}_{\mathrm{k}}\cap \mathrm{P}$ is a $(\mathrm{d}-1)$-dimensional disk. Therefore, we have a collection $\mathcal{D}=\{{\mathrm{D}}_{\mathrm{j}}\cap \mathrm{P}|\mathrm{j}=1,2,...,\mathrm{d}+1\}$ of $\mathrm{d}+1$ disks, each of dimension $\mathrm{d}-1$, such that every subset A of $\mathcal{D}$ consisting of d disks has the non-empty intersection property. By Helly’s Theorem, the intersection of all $(\mathrm{d}-1)$-disks in $\mathcal{D}$ is not empty. Therefore, there exists a point $\mathrm{q}\in \mathrm{D}$ with ${\mathrm{\pi }}_1(\mathrm{q})<{\mathrm{\pi }}_1(\mathrm{p})$. However, this contradicts the fact that $\mathrm{p}\in \mathrm{D}$ is such that ${\mathrm{\pi }}_1(\mathrm{p})=inf{\mathrm{\pi }}_1(\partial \mathrm{D})$.

Therefore, we conclude that ${\bigcap }_{\mathrm{j}=1}^{\mathrm{d}+1}\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M(\mathrm{j}))=\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M)$.

   □

Lemma 6.

Let M be a Vietoris-Rips system of $\mathrm{d}+1$ disks in ${\mathbb{R}}^{\mathrm{d}}$. If $\mathrm{D}={\cap }_{\mathrm{i}=1}^{\mathrm{d}+1}{\mathrm{D}}_{\mathrm{i}}=\varnothing$ and $\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M(\mathrm{j}))\ne \varnothing$ for each $\mathrm{j}=1,...,\mathrm{d}+1$, then ${\cap }_{\mathrm{i}=1}^{\mathrm{d}+1}\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M(\mathrm{i}))$ consists only of intervals of the form $[\mathrm{a},\mathrm{b}]$ with $\mathrm{a}>\mathrm{b}$ (inverted intervals).

Proof. 

Suppose that ${\cap }_{\mathrm{i}=1}^{\mathrm{d}+1}\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M(\mathrm{i}))$ contains an interval. Without loss of generality, let’s assume it is the interval

$$[inf{\mathrm{\pi }}_1(\partial \mathrm{D}(\mathrm{k})),sup{\mathrm{\pi }}_1(\partial \mathrm{D}(\mathrm{l}))]$$

for some $\mathrm{k},\mathrm{l}\in \{1,2,...,\mathrm{d}+1\}$ where $\mathrm{D}(\mathrm{k})={\cap }_{\mathrm{i}\ne \mathrm{k}}{\mathrm{D}}_{\mathrm{i}}$. By definition of ${\cap }_{\mathrm{i}=1}^{\mathrm{d}+1}\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M(\mathrm{i}))$, we have $inf{\mathrm{\pi }}_1(\partial \mathrm{D}(\mathrm{k}))\ge inf{\mathrm{\pi }}_1(\partial \mathrm{D}(\mathrm{i}))$ for all $\mathrm{i}\ne \mathrm{k}$ and $sup{\mathrm{\pi }}_1(\partial \mathrm{D}(\mathrm{l}))\le sup{\mathrm{\pi }}_1(\partial \mathrm{D}(\mathrm{j}))$ for all $\mathrm{j}\ne \mathrm{l}$. Now, consider a hyperplane $\mathrm{P}:{\mathrm{x}}_1=\mathrm{p}$ with $inf{\mathrm{\pi }}_1(\partial \mathrm{D}(\mathrm{k}))\le \mathrm{p}\le sup{\mathrm{\pi }}_1(\partial \mathrm{D}(\mathrm{k}))$. The hyperplane P intersects $\mathrm{D}(\mathrm{j})$ for every $\mathrm{j}=1,...,\mathrm{d}+1$ (since P cuts through all the boxes $\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M(\mathrm{j}))$). We also know that $\mathrm{P}\cap {\mathrm{D}}_{\mathrm{i}}$ is a $\mathrm{d}-1$-dimensional disk for each i.

Therefore, we have a collection of d disks in a $\mathrm{d}-1$ dimensional space, and by Helly’s Theorem, this collection must have a non-empty intersection. This implies that $\mathrm{D}\ne \varnothing$, which contradicts our assumption that D is empty. Hence, all intervals in ${\cap }_{\mathrm{i}=1}^{\mathrm{d}+1}\mathrm{A}\mathrm{A}\mathrm{B}\mathrm{B}(M(\mathrm{i}))$ must be inverted intervals of the form $[\mathrm{a},\mathrm{b}]$ with $\mathrm{a}>\mathrm{b}$.    □

As an example, we provide the lower bounds of the AABB for a system of three disks in ${\mathbb{R}}^{\mathrm{d}}$, this is, we compute $inf{\mathrm{\pi }}_{\mathrm{i}}(\partial \mathrm{D})$ for each $\mathrm{i}=1,...,\mathrm{d}$ (analogous $sup{\mathrm{\pi }}_{\mathrm{i}}\partial \mathrm{D}$) with $\mathrm{D}={\cap }_{\mathrm{j}=1}^3{\mathrm{D}}_{\mathrm{j}}$. Let $M=\{{\mathrm{D}}_1,{\mathrm{D}}_2,{\mathrm{D}}_3\}$ be a Rips system in ${\mathbb{R}}^{\mathrm{d}}$, then, for each $\mathrm{i}=1,...,\mathrm{d}$, the computation of the AABB for the disk system M is given by:

$$\begin{array}{c}\hfill inf{\mathrm{\pi }}_{\mathrm{i}}(\partial \mathrm{D})=\left\{\begin{array}{cc}{\mathrm{\pi }}_{\mathrm{i}}({\mathrm{c}}_{\mathrm{k}}-{\mathrm{r}}_{\mathrm{k}}{\mathrm{e}}_{\mathrm{i}})\hfill & \mathrm{if}{\mathrm{c}}_{\mathrm{k}}-{\mathrm{r}}_{\mathrm{k}}{\mathrm{e}}_{\mathrm{i}}\in \mathrm{D}\hfill \\ \underset{\mathrm{j}<\mathrm{k}}{max}\{inf{\mathrm{\pi }}_{\mathrm{i}}(\partial {\mathrm{D}}_{\mathrm{j}}\cap \partial {\mathrm{D}}_{\mathrm{k}})\}\hfill & \mathrm{if}(*)\hfill \\ inf{\mathrm{\pi }}_{\mathrm{i}}({\cap }_{\mathrm{j}=1}^3\partial {\mathrm{D}}_{\mathrm{j}})\hfill & \mathrm{otherwise}\hfill \end{array}\right.\end{array}$$

where (*) denotes the case where $\mathrm{q}\in \mathrm{D}$ with $\mathrm{q}\in \partial {\mathrm{D}}_{\mathrm{j}}\cap \partial {\mathrm{D}}_{\mathrm{k}}$ such that

$${\mathrm{\pi }}_{\mathrm{i}}(\mathrm{q})=\underset{\mathrm{j}<\mathrm{k}}{max}\{inf{\mathrm{\pi }}_{\mathrm{i}}(\partial {\mathrm{D}}_{\mathrm{j}}\cap \partial {\mathrm{D}}_{\mathrm{k}})\}.$$

The preceding calculations are a result of Lemma 2. It is known that the extremes of the AABB lie in the projections of specific poles, either from a d-sphere, $(\mathrm{d}-1)$-sphere, or the $(\mathrm{d}-2)$-sphere.

Given a disk system M in ${\mathbb{R}}^{\mathrm{d}}$, we can compute the minimal axis-aligned bounding box (AABB) for the intersection of all disks in M. If the AABB is a point, then it represents the intersection of the disks. This property allows us to identify when the AABB is a point. If the AABB of M is not a point, we can rescale M by a scale factor $\mathrm{\lambda }$ such that the AABB of $M_{\mathrm{\lambda }}$ becomes a point. The value of $\mathrm{\lambda }$ is referred to as the Čech scale of the system M.

4.3. Minimal axis-aligned bounding box Algorithm

Now, we present an algorithm to calculate the minimal axis-aligned bounding box (AABB) of a system of m disks in ${\mathbb{R}}^{\mathrm{d}}$. As previously explained, the AABB’s extremes are defined by projecting the poles of specific i-spheres. Thus, in computing the AABB, we will determine the poles for each i-sphere within the disk system.

The algorithm starts by initializing the set of poles, P, as an empty set. In each iteration of the first loop, we determine the spheres formed by the intersection of the boundaries of the disk subcollections in the system M (for this, we require the center and radius, which are computed in Section 2.2). Next, we identify the poles of each sphere and if any of them are present in all the disks of M, we add them to the set P. If P is not empty, we proceed to calculate the extremes for each dimension of the AABB using the set of north poles for the upper bounds and the set of south poles for the lower bounds. In this case, the output is the product of the intervals defined by the computed extremes. If P is empty, it indicates that there is no intersection in the disk system.

Algorithm 3: AABB.minimal