Intersection properties of finite disk collections

Jesús F. Espinoza [email protected] and Cynthia G. Esquer-Pérez [email protected] Departamento de Matemáticas, Universidad de Sonora, México

Abstract.

In this work we study the intersection properties of a finite disk system in the euclidean space. We accomplish this by utilizing subsets of spheres with varying dimensions and analyze specific points within them, referred to as poles. Additionally, we introduce two applications: estimating the common scale factor for the radii that makes the re-scaled disks intersects in a single point, this is the Čech scale, and constructing the minimal Axis-Aligned Bounding Box (AABB) that encloses the intersection of all disks in the system.

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}^{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 [4], 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}^{d}$ . We start by investigating the intersection of two disks in Subsection 2.1 and then expand our analysis to a system of $m$ disks in Subsection 2.2. By applying Helly’s Theorem, we prove that it is sufficient to examine the intersection of all subsystems consisting of $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 Subsection 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 Subsection 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 [5] and collision detection [2]. 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}^{d}$ with positive and not necessarily equal radii, i.e.,

M=\{D_{i}(c_{i};r_{i})\subset\mathbb{R}^{d}\mid c_{i}\in\mathbb{R}^{d},r_{i}>0,1\leq i\leq m<\infty\}.

Moreover, in order to study the intersection properties of a disk system $M$ with the approach addressed in Sections 4 and 3 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 D_{i}\subset\mathbb{R}^{d}\mid D_{i}\in M\},

where $\partial$ 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}^{d}$ is the intersection of a sphere with an affine subspace of dimension $i$ .

Of course, the notions of a sphere (as a $(d-1)$ -dimensional surface) and a $d$ -sphere in $\mathbb{R}^{d}$ agree. However, an $i$ -sphere in $\mathbb{R}^{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 $(d-1)$ -sphere in $\mathbb{R}^{d}$ . When another $d$ -sphere intersects this configuration, the result may be a $(d-1)$ -sphere, a $(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=\{D_{i}(c_{i};r_{i})\}$ composed of $m$ disks, where $\{c_{1},\dots,c_{m}\}$ is a set in general position in $\mathbb{R}^{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 $d-m+1$ , or equivalently, $i=d-m+1$ . This conclusion is drawn from [6, Theorem 2.1] and the fact that the affine hull of $\{c_{1},\dots,c_{m}\}$ is of dimension $m-1$ . Consequently, the following result is proven.

Lemma 2.

Let $M=\{D_{1}(c_{1};r_{1}),\ldots,D_{m}(c_{m};r_{m})\}$ be disk system such that $\{c_{1},\dots,c_{m}\}$ is a set in general position in $\mathbb{R}^{d}$ . Then, the possibilities for the set $\cap_{D_{i}\in M}\partial D_{i}$ are:

(1)

the empty set;
(2)

a single point;
(3)

a $(d-m+1)$ -sphere.

Remarkable points in $i$ -spheres that will play a key role in the rest of the article are the poles. Let $\pi_{i}:\mathbb{R}^{d}\longrightarrow\mathbb{R}$ be the canonical projection on the $i$ -th factor for $i=1,...,d$ , and let $\{e_{1},e_{2},...,e_{d}\}$ be the standard basis of $\mathbb{R}^{d}$ .

Definition 3.

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

We denote the $e_{q}$ -north pole of $S$ by $s_{q}^{+}$ and the $e_{q}$ -south pole by $s_{q}^{-}$ .

An $i$ -sphere can have a single $e_{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 $e_{q}$ . We are interested in finding the $e_{q}$ -poles of $(d-m+1)$ -spheres originating from disk systems $M=\{D_{1}(c_{1};r_{1}),\ldots,D_{m}(c_{m};r_{m})\}$ , by taking the intersection $\cap_{j=1}^{m}\partial D_{j}$ . Such $(d-m+1)$ -spheres will be denoted by $S_{M}(c;r)$ , to emphasize the disk system $M$ , as well as its center and radius.

Lemma 4.

Let $M=\{D_{1},\ldots,D_{m}\}$ be a $d$ -disk system such that $\bigcap_{j=1}^{m}D_{j}\neq\emptyset$ , and let $p$ be a point in $\bigcap_{j=1}^{m}D_{j}$ such that $\pi_{q}(p)\leq\pi_{q}(x)$ (resp. $\pi_{q}(p)\geq\pi_{q}(x)$ ) for every $x$ in $\bigcap_{j=1}^{m}D_{j}$ . Then, there exists an $i$ -sphere $S=\partial D_{j_{1}}\cap\cdots\cap\partial D_{j_{i}}$ such that $p$ is in $S$ and $p$ is the $e_{q}$ -south pole (resp. $e_{q}$ -north pole) of $S$ .

Proof.

Since $\cap_{j=1}^{m}D_{j}\neq\emptyset$ , then $\partial(\cap_{j=1}^{m}D_{j})\neq\emptyset$ , $\partial(\cap_{j=1}^{m}D_{j})\subset\cap_{j=1}^{m}D_{j}$ due to the closedness of the sets $D_{j}$ , for $j=1,\ldots,m$ , and $p\in\partial(\cap_{j=1}^{m}D_{j})$ .

On the other hand, since $\partial(\cap_{j=1}^{m}D_{j})\subset\cup_{i=1}^{m}\partial D_{j}$ , there exist indices $j_{1},\ldots,j_{i}$ such that $p\in\partial D_{j_{r}}$ for any $r=1,\ldots,i$ ; let $\Lambda(p)=\{j_{1},\ldots,j_{i}\}\subseteq\{1,\ldots,m\}$ be a maximal subset of indices such that $p\in D_{j}$ if and only if $j\in\Lambda(p)$ . We claim that $p\in\cap_{r=1}^{i}\partial D_{j_{r}}$ is the $e_{q}$ -south pole of $S:=\cap_{r=1}^{i}\partial D_{j_{r}}$ .

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

(1)

Every $x\in V_{p}\cap\partial(\cap_{j=1}^{m}D_{j})$ has as maximal set of indices a proper subset of $\Lambda(p)$ ,
(2)

$S\cap V_{p}\subset\partial(\cap_{j=1}^{m}D_{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 $\Lambda(p)$ . Therefore, $\pi_{q}(p)\leq\pi_{q}(x)$ for every $x\in S\cap V_{p}$ , which is equivalent to the fact that $\pi_{q}(p)\leq\pi_{q}(x)$ for every $x\in 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 $(d-1)$ -sphere given by the intersection of two disks in $\mathbb{R}^{d}$ .

Lemma 5.

Let $M=\{D_{1}(c_{1};r_{1}),D_{2}(c_{2};r_{2})\}$ be a disk system with two $d$ -disks such that $\partial D_{1}\cap\partial D_{2}$ is a $(d-1)$ -sphere $S=S_{M}(c;r)$ with center $c$ and radius $r$ . Then,

c=\frac{1}{2}\left(1+\frac{r_{2}^{2}-r_{1}^{2}}{\|c_{2}-c_{1}\|^{2}}\right)c_{1}+\frac{1}{2}\left(1+\frac{r_{1}^{2}-r_{2}^{2}}{\|c_{2}-c_{1}\|^{2}}\right)c_{2}

and

r=\frac{2\sqrt{s(s-\|c_{2}-c_{1}\|)(s-r_{1})(s-r_{2})}}{\|c_{2}-c_{1}\|}

where $s=\frac{1}{2}(\|c_{2}-c_{1}\|+r_{1}+r_{2})$ .

Proof.

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

(1)

\sum_{i=1}^{d}(k_{i}-h_{i})x_{i}-\frac{1}{2}\sum_{i=1}^{d}(k_{i}^{2}-h_{i}^{2})=\frac{r_{1}^{2}-r_{2}^{2}}{2},

where $c_{1}=(h_{1},\ldots,h_{d})$ and $c_{2}=(k_{1},\ldots,k_{d})$ . Then the normal vector of the hyperplane $\Pi$ is given by $N:=c_{2}-c_{1}=(k_{1}-h_{1},\ldots,k_{d}-h_{d})$ , and the center $c$ of $S$ is determined by the intersection point of the hyperplane $\Pi$ with the perpendicular line that passes through the center $c_{1}$ of $D_{1}$ . This line can be parameterized as $\gamma:t\mapsto c_{1}+tN=(x_{1}(t),\ldots,x_{d}(t))$ , such that $\gamma(0)=c_{1}$ and $\gamma(1)=c_{2}$ . We can compute the intersection point $c=\gamma(t_{*})$ of $\Pi$ and $\gamma([0,1])$ , for any $t_{*}\in[0,1]$ , by substituting it in (1),

\sum_{i=1}^{d}(k_{i}-h_{i})(h_{i}+t_{*}(k_{i}-h_{i}))=\frac{1}{2}\left(\sum_{i=1}^{d}(k_{i}^{2}-h_{i}^{2})+(r_{1}^{2}-r_{2}^{2})\right).

And solving for $t_{*}$ , we obtain that $t_{*}=\frac{1}{2}+\frac{r_{1}^{2}-r_{2}^{2}}{2\|c_{2}-c_{1}\|^{2}}$ . Hence, the center of $S$ is given by:

\displaystyle c=c_{1}+t_{*}N=\frac{1}{2}\left(1+\frac{r_{2}^{2}-r_{1}^{2}}{\|c_{2}-c_{1}\|^{2}}\right)c_{1}+\frac{1}{2}\left(1+\frac{r_{1}^{2}-r_{2}^{2}}{\|c_{2}-c_{1}\|^{2}}\right)c_{2}

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

r=\frac{2\sqrt{s(s-\|c_{2}-c_{1}\|)(s-r_{1})(s-r_{2})}}{\|c_{2}-c_{1}\|},

where $s=\frac{1}{2}(\|c_{2}-c_{1}\|+r_{1}+r_{2})$ correspond to the semi-perimeter. ∎

We can proceed now to compute the poles of the $(d-1)$ -sphere $\partial D_{1}\cap\partial D_{2}$ .

Lemma 6.

Let $D_{1}(c_{1};r_{1})$ and $D_{2}(c_{2};r_{2})$ be two $d$ -disks such that $\partial D_{1}\cap\partial D_{2}$ is a $(d-1)$ -sphere $S=S(c;r)$ with center $c$ and radius $r$ . Then, the $e_{q}$ -poles of $S$ are $s_{q}^{\pm}=c\pm\sum_{i}^{d}x_{i}e_{i}$ , where

x_{i}=\begin{cases}\dfrac{r|\pi_{i}(c_{2}-c_{1})\pi_{q}(c_{2}-c_{1})|}{\|c_{2}-c_{1}\|\sqrt{\|c_{2}-c_{1}\|^{2}-\pi_{q}(c_{2}-c_{1})^{2}}},&i\neq q\\ \\ \dfrac{r{\sqrt{\|c_{2}-c_{1}\|^{2}-\pi_{q}(c_{2}-c_{1})^{2}}}}{\|c_{2}-c_{1}\|},&i=q.\end{cases}

Proof.

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

	$\displaystyle\sum_{i=1}^{d}(k_{i}-h_{i})x_{i}$	$\displaystyle=0,$
	$\displaystyle\sum_{i=1}^{d}x_{i}^{2}$	$\displaystyle=r^{2},$

where $h_{i}:=\pi_{i}(c_{1})$ and $k_{i}:=\pi_{i}(c_{2})$ for $i=1,\ldots,d$ . In the case that $k_{q}-h_{q}=\pi_{q}(c_{2}-c_{1})=0$ , the normal vector $N=c_{2}-c_{1}$ of the hyperplane $\Pi$ is orthogonal to the basis vector $e_{q}$ . Therefore, the $e_{q}$ -poles of $S$ are $c\pm re_{q}$ , which agree with the formulae of the lemma.

On the other hand, suppose that $k_{q}-h_{q}\neq 0$ . To find the $e_{q}$ -poles of $S$ , we will use the Lagrange multiplier method. Consider the following function:

(2)

x_{q}=f(x_{1},x_{2},...,\widehat{x_{q}},...,x_{d})=\frac{-\sum_{j\neq q}^{d}(k_{j}-h_{j})x_{j}}{k_{q}-h_{q}},

subject to the restriction:

g(x_{1},x_{2},...,\widehat{x_{q}},...,x_{d})=\sum_{j\neq q}^{d}x_{j}^{2}+\left(\frac{-\sum_{j\neq q}^{d}(k_{j}-h_{j})x_{j}}{k_{q}-h_{q}}\right)^{2}-r^{2}=0

Let $\lambda$ be the Lagrange multiplier, we define

h(x_{1},...,\widehat{x_{q}},...,x_{d},\lambda)=f(x_{1},...,\widehat{x_{q}},...,x_{d})+\lambda g(x_{1},...,\widehat{x_{q}},...,x_{d})

For any $i\neq q$ , consider the following system of equations:

\frac{\partial h}{\partial x_{i}}=-\frac{k_{i}-h_{i}}{k_{q}-h_{q}}+2\lambda x_{i}+2\lambda\left(\frac{-\sum_{j\neq q}^{d}(k_{j}-h_{j})x_{j}}{k_{q}-h_{q}}\right)\left(-\frac{k_{i}-h_{i}}{k_{q}-h_{q}}\right)=0.

Then

	$\displaystyle-\frac{k_{i}-h_{i}}{k_{q}-h_{q}}+2\lambda\left(x_{i}+\frac{k_{i}-h_{i}}{(k_{q}-h_{q})^{2}}{\sum_{j\neq q}^{d}(k_{j}-h_{j})x_{j}}\right)$	$\displaystyle=0$
	$\displaystyle-(k_{i}-h_{i})(k_{q}-h_{q})+2\lambda\left((k_{q}-h_{q})^{2}x_{i}+(k_{i}-h_{i}){\sum_{j\neq q}^{d}(k_{j}-h_{j})x_{j}}\right)$	$\displaystyle=0$

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

\lambda=\frac{(k_{i}-h_{i})(k_{q}-h_{q})}{2\left((k_{q}-h_{q})^{2}x_{i}+(k_{i}-h_{i}){\sum_{j\neq q}^{d}(k_{j}-h_{j})x_{j}}\right)}

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

	$\displaystyle x_{i}^{2}$	$\displaystyle=\frac{r^{2}(k_{i}-h_{i})^{2}}{(k_{\tilde{i}}-h_{\tilde{i}})^{2}+\sum_{j\neq\tilde{i},q}{(k_{j}-h_{j})}^{2}+\frac{1}{(k_{q}-h_{q})^{2}}\left(\sum_{j\neq q}(k_{j}-h_{j})^{2}\right)^{2}}$
		$\displaystyle=\frac{r^{2}(k_{i}-h_{i})^{2}}{\sum_{j\neq q}{(k_{j}-h_{j})}^{2}+\frac{1}{(k_{q}-h_{q})^{2}}\left(\sum_{j\neq q}(k_{j}-h_{j})^{2}\right)^{2}}$
		$\displaystyle=\frac{r^{2}(k_{i}-h_{i})^{2}(k_{q}-h_{q})^{2}}{\left((k_{q}-h_{q})^{2}+\sum_{j\neq q}(k_{j}-h_{j})^{2}\right)\left(\sum_{j\neq q}{(k_{j}-h_{j})}^{2}\right)}$
		$\displaystyle=\frac{r^{2}\pi_{i}(c_{2}-c_{1})^{2}\pi_{q}(c_{2}-c_{1})^{2}}{\\|c_{2}-c_{1}\\|^{2}\left(\\|c_{2}-c_{1}\\|^{2}-\pi_{q}(c_{2}-c_{1})^{2}\right)}$

Finally, for $i=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 $(d-m+1)$ -sphere $S=\cap_{j=1}^{m}\partial D_{j}$ . We can achieve this by considering the disk system translated to $c_{m}$ , denoted as $\{D_{j}(c_{j}-c_{m};r_{j})\}_{j=1}^{m}$ , and by defining the $(d-m+1)$ -sphere $S-\{c_{m}\}=\cap_{j=1}^{m}\partial D_{j}(c_{j}-c_{m};r_{j})$ . This sphere is positioned at the intersection of hyperplanes (for more details, refer to [6]).

(3)

(c_{k}-c_{m})^{T}x=\frac{1}{2}(r_{m}^{2}+\|c_{k}-c_{m}\|^{2}-r_{k}^{2})

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

\sum_{j=1}^{m-1}\lambda_{j}(c_{k}-c_{m})\cdot(c_{j}-c_{m})=\frac{1}{2}(r_{m}^{2}+\|c_{k}-c_{m}\|^{2}-r_{k}^{2})

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

c=\lambda_{1}(c_{1}-c_{m})+...+\lambda_{m-1}(c_{m-1}-c_{m})+c_{m}

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

r^{2}=r_{k}^{2}-\|c-c_{k}\|^{2}

for any $k\in\{1,...,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 $e_{q}$ -poles for each $q\in\{1,2,...,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}^{d}$ , and let $n_{1},...,n_{d-i}$ be orthogonal vectors to the affine space $L$ that contains $S$ . Consider the space $M$ generated by these vectors together with the vector $e_{q}$ from the canonical basis of $\mathbb{R}^{d}$ . Let us denote $n_{j}=\left(n_{l}^{(j)}\right)_{l=1}^{d}$ for each $j=1,...,d-i$ . Then, we can define $L_{0}$ , the set $L$ translated to the origin, as follows:

	$\displaystyle L_{0}$	$\displaystyle=\left<\{n_{j}\}_{j=1}^{d-i}\right>^{\perp}$
		$\displaystyle=\left\{x\in\mathbb{R}^{d}\mid x\cdot n_{j}=0,\hskip 5.69046pt\forall j=1,...,d-i\right\}$

The set $M$ is defined as:

	$\displaystyle M$	$\displaystyle=\left<\{n_{j}\}_{j=1}^{d-i}\cup\{e_{q}\}\right>$
		$\displaystyle=\left\{x=\sum_{j=1}^{d-i}\lambda_{j}n_{j}+\lambda_{d-i+1}e_{q}\in\mathbb{R}^{d}\mid\lambda_{j}\in\mathbb{R},\hskip 5.69046pt\forall j=1,...,d-i+1\right\}$
		$\displaystyle=\left\{\left(\sum_{j=1}^{d-i}\lambda_{j}n_{1}^{(j)},...,\sum_{j=1}^{d-i}\lambda_{j}n_{q}^{(j)}+\lambda_{d-i+1},...,\sum_{j=1}^{d-i}\lambda_{j}n_{d}^{(j)}\right)\mid\lambda_{j}\in\mathbb{R}\right\}$

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

Refer to caption — Figure 1. Visualization of the subspaces $L$ and $M+\{c\}$ .

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

	$\displaystyle M\cap L_{0}$	$\displaystyle=\left\{x\in M\mid x\cdot n_{k}=0,\hskip 5.69046pt\forall k=1,...,d-i\right\}$
		$\displaystyle=\left\{\sum_{j=1}^{d-i}\lambda_{j}n_{j}+\lambda_{d-i+1}e_{q}\in\mathbb{R}^{d}\left\|\left(\sum_{j=1}^{d-i}\lambda_{j}n_{j}+\lambda_{d-i+1}e_{q}\right)\cdot n_{k}=0,\hskip 5.69046pt\lambda_{j}\in\mathbb{R}\hskip 5.69046pt\forall k=1,...,d-i\right.\right\}$
		$\displaystyle=\left\{\sum_{j=1}^{d-i}\lambda_{j}n_{j}+\lambda_{d-i+1}e_{q}\in\mathbb{R}^{d}\left\|\sum_{j=1}^{d-i}\lambda_{j}n_{j}\cdot n_{k}+\lambda_{d-i+1}n_{q}^{k}=0,\hskip 5.69046pt\lambda_{j}\in\mathbb{R}\hskip 5.69046pt\forall k=1,...,d-i\right.\right\}$

Let us consider a disk system in $\mathbb{R}^{d}$ , denoted $\{D_{j}(c_{j};r_{j})\}_{j=1}^{m}$ , where $m<d$ . The intersection of their boundaries forms a $(d-m+1)$ -sphere $S$ . In this case, the subspace $M$ has dimension $m$ , or $dim(M)=m-1$ if $e_{q}\in M$ . We choose the normal vectors for the affine space containing $S$ as $n_{j}=c_{j}-c_{m}$ , where $j=1,...,m-1$ . Then

M=\left\{\left(\sum_{j=1}^{m-1}\lambda_{j}\left(c_{1}^{(j)}-c_{1}^{(m)}\right),...,\sum_{j=1}^{m-1}\lambda_{j}\left(c_{q}^{(j)}-c_{q}^{(m)}\right)+\lambda_{m},...,\sum_{j=1}^{m-1}\lambda_{j}\left(c_{d}^{(j)}-c_{d}^{(m)}\right)\right)\left|\lambda_{j}\in\mathbb{R}\right.\right\}

By rewriting, we have

\displaystyle M\cap L_{0}

\displaystyle=\left\{\sum_{j=1}^{m-1}\lambda_{j}n_{j}+\lambda_{m}e_{q}\in\mathbb{R}^{d}\left|\sum_{j=1}^{m-1}\lambda_{j}n_{j}\cdot n_{k}+\lambda_{m}n_{q}^{k}=0,\hskip 5.69046pt\lambda_{j}\in\mathbb{R},\hskip 5.69046pt\forall k=1,...,m-1\right.\right\}

If $S(c;r)=\cap_{j=1}^{m}\partial D_{j}$ is the $(d-m+1)$ -sphere with center in $c$ and radius $r$ , then the $e_{q}$ -poles of $S$ are the $e_{q}$ -poles of $S-\{c_{m}\}$ but translated by $c$ . The poles of $S-\{c_{m}\}$ are located in $M\cap L_{0}$ . If $p$ is an $e_{q}$ -pole of $S-\{c_{m}\}$ , then it can be expressed as

p=\sum_{j=1}^{m-1}\lambda_{j}n_{j}+e_{q}\lambda_{m}

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

\sum_{j=1}^{m-1}\lambda_{j}n_{j}\cdot n_{k}+\lambda_{m}n_{q}^{k}=0

for each $k=1,...,m-1$ and

\|p\|^{2}=r^{2}

Thus, if $p=\sum_{j=1}^{m-1}\lambda_{j}n_{j}+e_{q}\lambda_{m}$ is an $e_{q}$ -pole of $S-\{c_{m}\}$ , the following equations are satisfied for $\lambda_{1},\lambda_{2},...,\lambda_{m}\in\mathbb{R}$ :

(4)

\sum_{j=1}^{m-1}\lambda_{j}n_{j}\cdot n_{k}+\lambda_{m}n_{q}^{k}=0

(5)

\sum_{i=1}^{d}\left(\sum_{j=1}^{m-1}\lambda_{j}n_{i}^{j}\right)^{2}+2\lambda_{m}\sum_{j=1}^{m-1}\lambda_{j}n_{q}^{j}+\lambda_{m}^{2}=r^{2}

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

\begin{pmatrix}n_{1}\cdot n_{1}&n_{1}\cdot n_{2}&....&n_{1}\cdot n_{m-1}\\ n_{2}\cdot n_{1}&n_{2}\cdot n_{2}&....&n_{2}\cdot n_{m-1}\\ ...\\ n_{m-1}\cdot n_{1}&n_{m-1}\cdot n_{2}&....&n_{m-1}\cdot n_{m-1}\\ \end{pmatrix}\begin{pmatrix}\lambda_{1}\\ \lambda_{2}\\ ...\\ \lambda_{m-1}\end{pmatrix}+\lambda_{m}\begin{pmatrix}n_{q}^{1}\\ n_{q}^{2}\\ ...\\ n_{q}^{m-1}\end{pmatrix}=0.

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

\lambda_{m}^{2}\left[\sum_{i=1}^{d}\left(\sum_{j=1}^{m-1}(A^{-1}B)[j]n_{i}^{j}\right)^{2}+2\sum_{j=1}^{m-1}(A^{-1}B)[j]n_{q}^{j}+1\right]-r^{2}=0

Let us define

\Gamma_{i}=\sum_{j=1}^{m-1}(A^{-1}B)[j]n_{i}^{j}

for all $i=1,...,d$ . Solving this equation, we find:

\lambda_{m}=\frac{\pm r}{\sqrt{\sum_{i=1}^{d}\Gamma_{i}^{2}+2\Gamma_{q}+1}}

\lambda_{j}=\frac{\pm r(A^{-1}B)[j]}{\sqrt{\sum_{i=1}^{d}\Gamma_{i}^{2}+2\Gamma_{q}+1}}

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

p=\sum_{j=1}^{m-1}\lambda_{j}n_{j}+e_{q}\lambda_{m}+c=\frac{\pm r}{\sqrt{\sum_{i=1}^{d}\Gamma_{i}^{2}+2\Gamma_{q}+1}}\left(\sum_{j=1}^{m-1}(A^{-1}B)[j]n_{j}+e_{q}\right)+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 7.

Let $M=\{D_{1},D_{2},\ldots,D_{m}\}$ be a $d$ -disk system. We say $M$ is a Vietoris-Rips system if $D_{i}\cap D_{j}\neq\emptyset$ for each pair $i,j\in\{1,2,\ldots,m\}$ . Furthermore, if the $d$ -disk system $M$ has the nonempty intersection property $\bigcap_{D_{i}\in M}D_{i}\neq\emptyset$ , then $M$ is called a Čech system.

For each $\lambda\geq 0$ , we define a collection of $d$ -disks $M_{\lambda}$ with the same centers as those in the $d$ -disk system $M$ , but with radii rescaled by $\lambda$ . When $\lambda>0$ , $M_{\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 $\nu_{M}$ , of a $d$ -disk system $M$ is the smallest $\lambda\in\mathbb{R}$ such that $M_{\lambda}$ is a Vietoris-Rips system. Similarly, the Čech scale $\mu_{M}$ , of $M$ is the smallest $\lambda\in\mathbb{R}$ such that $M_{\lambda}$ is a Čech system. This is

	$\displaystyle\nu_{M}$	$\displaystyle=\inf\{\lambda\in\mathbb{R}\mid M_{\lambda}\mbox{ is a Vietoris-Rips system}\}$
	$\displaystyle\mu_{M}$	$\displaystyle=\inf\{\lambda\in\mathbb{R}\mid M_{\lambda}\mbox{ is a \v{C}ech system}\}$

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 $\nu_{M}\leq 1$ (in particular $\nu_{M_{\nu_{M}}}=1$ ); similarly, $M$ is a Čech system if and only if $\mu_{M}\leq 1$ .

Note that for a given $d$ -disk system $M=\{D_{1},D_{2},\ldots,D_{m}\}$ the Vietoris-Rips scale is

\nu_{M}=\max_{i<j}\{\|c_{i}-c_{j}\|/(r_{i}+r_{j})\}

where $c_{i}$ and $r_{i}$ are the center and radii of $D_{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 8.

If $\mu_{M}$ is the Čech scale for $M$ , then the $\mu_{M}$ -rescaled system $M_{\mu_{M}}$ , has only one point in the intersection $\bigcap_{D_{i}\in M}D_{i}(c_{i};\mu_{M}r_{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{2d/(d+1)}$ is also a Čech system. This is established by the following lemma, the proof of which can be found in [3].

Lemma 9.

Let $M=\{D_{i}(c_{i};r_{i})\}$ be a $d$ -disk system in euclidean space $\mathbb{R}^{d}$ . If $D_{i}(c_{i};r_{i})\cap D_{j}(c_{j};r_{j})\neq\emptyset$ for every pair of disks in $M$ , then

\bigcap_{D_{i}\in M}D_{i}(c_{i};\sqrt{2d/(d+1)}\,r_{i})\neq\emptyset.

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 $\nu_{M}$ . This is stated by the following corollary.

Corollary 10.

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

3.1. Algorithm for determining Čech system.

In the previous section, we have determined the $e_{q}$ -poles for the intersection of any number of disks in $\mathbb{R}^{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}^{d}$ where $m>d$ , it is enough to verify if every subsystem of $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 ( $(d-2)-sphere$ ), continuing the process iteratively. If a pole is found within the remaining disks, the system is classified as a Čech system.

Input : A

d

-disk system

M=\{D_{j}\}_{j=1}^{m}

Output : A logical TRUE/FALSE to indicate if

M

is a Čech system

2 Initialize:

\texttt{Is\_Cech\_System}\leftarrow\texttt{FALSE}

3 for $k\leftarrow 1$ to $m$ do

4 Let

\mathcal{S}

be the set of

(d-k+1)

-spheres of

\partial M

5 for $S$ in $\mathcal{S}$ do

6 for $q\leftarrow 1$ to $d$ do

7 Compute the set

\{s_{q}^{\pm}\}

e_{q}

-poles of

S

8 for $s\leftarrow\{s_{q}^{\pm}\}$ do

9 if $s\in\cap_{j=1}^{d}D_{j}$ then

\texttt{Is\_Cech\_System}\leftarrow\texttt{TRUE}

11 Go to line 16

12 end if

14 end for

16 end for

18 end for

20 end for

return (Is_Cech_System)

Algorithm 1 Cech.system

Theorem 11.

Let $M=\{D_{1},\ldots,D_{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_{j=1}^{m}D_{j}$ , therefore $\cap_{j=1}^{m}D_{j}\neq\emptyset$ and it follows that $M$ is a Čech system.

On the other hand, if $\bigcap_{j=1}^{m}D_{j}\neq\emptyset$ , let $p$ be a point in $\bigcap_{j=1}^{m}D_{j}$ satisfying $\pi_{1}(p)\leq\pi_{1}(x)$ for every $x$ in $\bigcap_{j=1}^{m}D_{j}$ . By Lemma 4, it follows that $p$ belongs to an $i$ -sphere and must be an $e_{1}$ -south pole. Therefore, by the exhaustive search of Algorithm 1 across all poles, its output is Cech.system $(M)=$ TRUE. ∎

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}^{d}$ .

Input : A

d

-disk system

M

\mathbb{R}^{d}

and precision parameter

\eta>0

Output :

\mu_{M}

, a Čech scale approximation

2 Compute the Vietoris-Rips scale of

M

\nu_{M}

3 if Cech.system $(M_{\nu_{M}})=$ TRUE then

\mu_{M}\leftarrow\nu_{M}

5 Go to line 16

7else

8 Initialize:

\mu_{M}^{*}\leftarrow\nu_{M}

\mu_{M}\leftarrow\sqrt{2d/(d+1)}\nu_{M}

9 while $\mu_{M}-\mu_{M}^{*}>\eta$ do

10 Compute:

\lambda\leftarrow\dfrac{\mu_{M}^{*}+\mu_{M}}{2}

11 if Cech.system $(M_{\lambda})=$ TRUE then

12 Update:

\mu_{M}\leftarrow\lambda

13 else

14 Update:

\mu_{M}^{*}\leftarrow\lambda

15 end if

17 end while

19 end if

return (

\mu_{M}

)

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 $\eta>0$ . It initializes the scale factor $\lambda$ to the Rips scale $\nu_{M}$ . If this scale satisfies Cech.system $(M_{\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 $\lambda$ . The Čech scale is known to fall between the Rips scale and the value $\sqrt{{2d}/{(d+1)}}\nu_{M}$ (Generalized Vietoris-Rips, Corollary 10). To approximate the Čech scale, we employ the bisection method as long as the interval enclosing the Čech scale has a length greater than $\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 $\mathscr{C}(M)$ denote the set of Čech subsystems, and $\mathscr{C}_{\lambda}(M)$ the set of Čech subsystems for the rescaled disk system $M_{\lambda}$ . The Čech filtration of the $M$ system forms a maximal chain of Čech complexes

\mathscr{C}_{*}(M):\mathscr{C}_{0}(M)\subsetneq\mathscr{C}_{\lambda_{1}}(M)\subsetneq\mathscr{C}_{\lambda_{2}}(M)\subsetneq...\subsetneq\mathscr{C}_{\mu_{M}}(M),

where each $\lambda_{i}$ represents the Čech scale of the system $M_{\lambda_{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, $\mathscr{C}_{\lambda}(M)$ , simply requires determining the Čech scale of the system $M_{\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 12.

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

AABB(M):=\bigcap_{D\subset\tilde{B}}\tilde{B}

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

Note that the AABB can be expressed as

AABB(M)=\prod_{k=1}^{d}[\inf\pi_{k}(\partial D),\sup\pi_{k}(\partial D)]

where $\pi_{k}:\mathbb{R}^{d}\longrightarrow\mathbb{R}$ is the canonical projection onto the $k$ -th factor, and $\partial 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 $D_{1}$ and $D_{2}$ in $\mathbb{R}^{d}$ . If $D_{1}\cap D_{2}\neq\emptyset$ and $D_{1}\neq D_{2}$ , the subset $\partial D_{1}\cap\partial D_{2}$ can take one of three forms: an empty set, a single common point (when the disks are tangent), or a $(d-2)$ -dimensional sphere. In the last case, we denote the $(d-2)$ -dimensional sphere or the $(d-1)$ -sphere $\partial D_{1}\cap\partial D_{2}$ by $S_{1,2}$ . To calculate $\inf\pi_{i}(\partial(D_{1}\cap D_{2}))$ and $\sup\pi_{i}(\partial(D_{1}\cap D_{2}))$ for each $i\in\{1,2,...,d\}$ , we can use:

(6)

\displaystyle\inf\pi_{i}(\partial(D_{1}\cap D_{2}))=\begin{cases}\pi_{i}(c_{1}-r_{1}e_{i})&\text{if }c_{1}-r_{1}e_{i}\in D_{2},\\ \pi_{i}(c_{2}-r_{2}e_{i})&\text{if }c_{2}-r_{2}e_{i}\in D_{1},\\ \inf(\pi_{i}(\partial D_{1}\cap\partial D_{2}))&\text{otherwise},\end{cases}

\displaystyle\sup\pi_{i}(\partial(D_{1}\cap D_{2}))=\begin{cases}\pi_{i}(c_{1}+r_{1}e_{i})&\text{if }c_{1}+r_{1}e_{i}\in D_{2},\\ \pi_{i}(c_{2}+r_{2}e_{i})&\text{if }c_{2}+r_{2}e_{i}\in D_{1},\\ \sup(\pi_{i}(\partial D_{1}\cap\partial D_{2}))&\text{otherwise}.\end{cases}

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

To simplify the notation, we will use $B_{i,j}$ to denote the AABB of the intersection of disks $D_{i}$ and $D_{j}$ . Figure 2 illustrates the AABB of the intersection of two disks in the plane.

According to (6) and Lemma 5, 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}^{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 13.

Let

M=\{D_{1}((4,1,0);\sqrt{2}),D_{2}((4,-1,0);\sqrt{2}),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 $xy$ -plane:

By computing the boxes $B_{i,j}$ for all $i<j$ , we obtain:

	$\displaystyle B_{1,2}$	$\displaystyle=[3,5]\times[-\sqrt{2},\sqrt{2}]\times[-1,1]$
	$\displaystyle B_{1,3}$	$\displaystyle=[4-\sqrt{2},3]\times[0,1.41]\times[-0.72,0.72]$
	$\displaystyle B_{2,3}$	$\displaystyle=[4-\sqrt{2},3]\times[-1.4,0]\times[-0.72,0.72]$

Therefore, $\cap_{i<j}B_{i,j}=[3,3]\times[0,0]\times[-0.72,0.72]$ . However, the disk system intersects at the point $P=(3,0,0)$ , which means $AABB(M)=\{P\}$ . In other words, $AABB(M)$ is not equal to $\cap_{i<j}B_{i,j}$ .

We know that if $D\neq\emptyset$ , then $\cap_{i<j}B_{i,j}\neq\emptyset$ . However, the converse is not always true. The following example illustrates this fact.

Example 14.

Let $N=\{D_{1},D_{2},D_{3}((0,1,0);\sqrt{10}),D_{4}((3,0,1);0.9)\}\subset\mathbb{R}^{3}$ be a Vietoris-Rips system. The projection of disks $D_{1},D_{2},$ and $D_{3}$ onto the $xy$ -plane is illustrated in Figure 4.

We will now compute the intersections $B_{i,j}$ for different pairs of disks:

	$\displaystyle B_{1,2}$	$\displaystyle=[3,5]\times[-\sqrt{2},\sqrt{2}]\times[-1,1]$
	$\displaystyle B_{1,3}$	$\displaystyle=[4-\sqrt{2},\sqrt{10}]\times[0,2]\times[-1,1]$
	$\displaystyle B_{1,4}$	$\displaystyle=[2.6,3.9]\times[-0.4,0.9]\times[0.100007,1.3]$
	$\displaystyle B_{2,3}$	$\displaystyle=[2.6,3]\times[-0.8,0]\times[-0.44,4.47]$
	$\displaystyle B_{2,4}$	$\displaystyle=[2.6,3.9]\times[-0.9,0.4]\times[0.100007,1.3]$
	$\displaystyle B_{3,4}$	$\displaystyle=[2.1,3.11]\times[-0.73,0.9]\times[0.1,1.7]$

The intersection of all pairwise intersections, $\cap_{i<j}B_{i,j}$ , is given by $[3,3]\times[0,0]\times[0.1,1.7]$ . However, the intersection of disks $D_{1},D_{2},$ and $D_{3}$ is a single point $P$ , which is not contained in $D_{4}$ (by construction). Therefore, $D$ is an empty set, but $\cap_{i<j}B_{i,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 13, we observe that the intersection of $AABB(\{D_{i},D_{j}\})$ is not equal to the AABB of the disk system $M$ . Similarly, in Example 14, we find that the intersection of three disks is empty, yet the intersection of $AABB(\{D_{i},D_{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}^{d}$ .

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

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

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

In the case where $m=d+1$ , we consider $AABB(M(1)),...,AABB(M(m))$ as a collection of minimal axis-aligned boxes for the disk system $M(i)=M-\{D_{i}\}$ in $\mathbb{R}^{d}$ . If the intersection of any $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 $AABB(M)=\cap AABB(M(i))$ . The next theorem confirms this finding.

Theorem 15 (Helly’s theorem for minimal axis-aligned boxes).

Let $M$ be a Rips system with $d+1$ disks in $\mathbb{R}^{d}$ . Then, the minimal axis-aligned bounding box for the intersection set $D=\cap_{j=1}^{d+1}D_{j}$ satisfies:

AABB(M)=\bigcap_{j=1}^{d+1}AABB(M(j))

where $M(j)=M-\{D_{j}\}$ .

Proof.

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

AABB(M)=\prod_{i=1}^{d}\left[\inf\pi_{i}(\partial D),\sup\pi_{i}(\partial D)\right]

and

\bigcap_{j=1}^{d+1}AABB(M(j))=\prod_{i=1}^{d}\left[\max_{k}\{\inf\pi_{i}(\partial\cap_{j\neq k}D_{j})\},\min_{k}\{\sup\pi_{i}(\partial\cap_{j\neq k}D_{j})\}\right]

Without loss of generality, let’s assume that:

\pi_{1}(p)>\max_{k}\left\{\inf\pi_{1}(\cap_{j\neq k}\partial D_{j})\right\}_{k=1}^{d+1}

where $p\in D$ satisfies $\pi_{1}(p)=\inf\pi_{1}(\partial D)$ .

Let $q_{j}$ be the point in $\cap_{k\neq j}D_{k}$ such that $\pi_{1}(q_{j})=\inf\pi_{1}(\cap_{k\neq j}D_{k})$ for each $j=1,2,...,d+1$ . Note that $q_{j}\notin D_{j}$ because $q_{j}$ is not in $D$ . Now, let $\gamma_{j}$ be the line segment that connects $q_{j}$ and $p$ .

Choose $\epsilon>0$ small enough such that the hyperplane $P:x_{1}=\pi_{1}(p)-\epsilon$ does not contain any $q_{j}$ , and $\pi_{1}(q_{j})<\pi_{1}(p)-\epsilon$ for all $j=1,2,...,d+1$ (such hyperplane $P$ exists because $\pi_{1}(p)>\pi_{1}(q_{j})$ for all $j=1,...,d+1$ ). Since $q_{j}$ is in $\cap_{k\neq j}D_{k}$ , the hyperplane $P$ intersects every disk in $M(j)$ for each $j$ , and $\gamma_{j}\subset\cap_{k\neq j}D_{k}$ intersects $P$ at a point that is in $\cap_{k\neq j}D_{k}$ (see Figure 5). Furthermore, $D_{k}\cap P$ is a $(d-1)$ -dimensional disk. Therefore, we have a collection $\mathcal{D}=\{D_{j}\cap P|j=1,2,...,d+1\}$ of $d+1$ disks, each of dimension $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 $(d-1)$ -disks in $\mathcal{D}$ is not empty. Therefore, there exists a point $q\in D$ with $\pi_{1}(q)<\pi_{1}(p)$ . However, this contradicts the fact that $p\in D$ is such that $\pi_{1}(p)=\inf\pi_{1}(\partial D)$ .

Therefore, we conclude that $\bigcap_{j=1}^{d+1}AABB(M(j))=AABB(M)$ .

∎

Lemma 16.

Let $M$ be a Vietoris-Rips system of $d+1$ disks in $\mathbb{R}^{d}$ . If $D=\cap_{i=1}^{d+1}D_{i}=\emptyset$ and $AABB(M(j))\neq\emptyset$ for each $j=1,...,d+1$ , then $\cap_{i=1}^{d+1}AABB(M(i))$ consists only of intervals of the form $[a,b]$ with $a>b$ (inverted intervals).

Proof.

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

[\inf\pi_{1}(\partial D(k)),\sup\pi_{1}(\partial D(l))]

for some $k,l\in\{1,2,...,d+1\}$ where $D(k)=\cap_{i\neq k}D_{i}$ . By definition of $\cap_{i=1}^{d+1}AABB(M(i))$ , we have $\inf\pi_{1}(\partial D(k))\geq\inf\pi_{1}(\partial D(i))$ for all $i\neq k$ and $\sup\pi_{1}(\partial D(l))\leq\sup\pi_{1}(\partial D(j))$ for all $j\neq l$ . Now, consider a hyperplane $P:x_{1}=p$ with $\inf\pi_{1}(\partial D(k))\leq p\leq\sup\pi_{1}(\partial D(k))$ . The hyperplane $P$ intersects $D(j)$ for every $j=1,...,d+1$ (since $P$ cuts through all the boxes $AABB(M(j))$ ). We also know that $P\cap D_{i}$ is a $d-1$ -dimensional disk for each $i$ .

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

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

\displaystyle\inf\pi_{i}(\partial D)=\begin{cases}\pi_{i}(c_{k}-r_{k}e_{i})&\text{if $c_{k}-r_{k}e_{i}\in D$}\\ \max_{j<k}\{\inf\pi_{i}(\partial D_{j}\cap\partial D_{k})\}&\text{if (*)}\\ \inf\pi_{i}(\cap_{j=1}^{3}\partial D_{j})&\text{otherwise}\end{cases}

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

\pi_{i}(q)=\max_{j<k}\{\inf\pi_{i}(\partial D_{j}\cap\partial D_{k})\}.

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

Given a disk system $M$ in $\mathbb{R}^{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 $\lambda$ such that the AABB of $M_{\lambda}$ becomes a point. The value of $\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}^{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.

Input : A

d

-disk system

M=\{D_{j}\}_{j=1}^{m}

Output : The minimal AABB for the system

M

2 Initialize:

P=\emptyset

3 for $k\leftarrow 1$ to $m$ do

4 Let

\mathcal{S}

be the set of

(d-k+1)

-spheres of

\partial M

5 for $S$ in $\mathcal{S}$ do

6 for $q\leftarrow 1$ to $d$ do

7 Compute the set

\{s_{q}^{\pm}\}

e_{q}

-poles of

S

8 for $s\leftarrow\{s_{q}^{\pm}\}$ do

9 if $s\in\cap_{j=1}^{d}D_{j}$ then

10 Add:

P\leftarrow s

12 end if

14 end for

16 end for

18 end for

20 end for

21if $P\neq\emptyset$ then

22 for $q\leftarrow 1$ to $d$ do

a_{q}=min_{P}\{\pi_{q}(s_{q}^{-})\}

b_{q}=max_{P}\{\pi_{q}(s_{q}^{+})\}

26 end for

27 return (

\Pi_{q=1}^{d}[a_{q},b_{q}]

)

28else

29 return (The disk system does not intersect)

30 end if

Algorithm 3 AABB.minimal

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 Subsection 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.

5. Acknowledgements

C.G.E.P acknowledges CONACYT for the financial support provided through a National Fellowship (CVU-638165).

References

[1] G. Bell, A. Lawson, J. Martin, J. Rudzinski, and C. Smyth, Weighted persistent homology, Involve, a Journal of Mathematics, 12 (2019), pp. 823–837.
[2] P. Cai, Y. Cai, I. Chandrasekaran, and J. Zheng, Collision detection using axis aligned bounding boxes, Simulation, Serious Games and Their Applications, (2013), pp. 1–14.
[3] J. F. Espinoza, R. Hernández-Amador, H. A. Hernández-Hernández, and B. Ramonetti-Valencia, A numerical approach for the filtered generalized čech complex, Algorithms, 13 (2019), p. 11.
[4] N. K. Le, P. Martins, L. Decreusefond, and A. Vergne, Construction of the generalized cech complex, 2015.
[5] J. Mahovsky and B. Wyvill, Fast ray-axis aligned bounding box overlap tests with plucker coordinates, Journal of Graphics Tools, 9 (2004), pp. 35–46.
[6] D. S. Maioli, C. Lavor, and D. S. Gonçalves, A note on computing the intersection of spheres in $\mathbb{R}^{n}$ , The ANZIAM Journal, 59 (2017), p. 271–279.