---
title: Radial Persistence Transform in Graph & Manifold Analysis
url: https://www.emergentmind.com/topics/radial-persistence-transform
type: topic
---

# Radial Persistence Transform in Graph & Manifold Analysis

Searching arXiv for the cited papers to ground the article.
“Radial Persistence Transform” is not a formal term used in the cited papers, but the phrase accurately describes two closely related constructions in topological shape analysis. In the graph setting, it denotes a radial organization of directions and associated augmented persistence information around each vertex, together with a radial binary multi-search that reconstructs edges from a finite subset of augmented persistence diagrams [2212.13206]. In the manifold-with-boundary setting, the Extended Persistent Homology Transform (XPHT) is defined for height functions indexed by directions, and the paper explicitly states that “similar results could hold for other kinds of functions, such as radial functions,” which suggests a radial analogue built from extended persistence of distance-type filtrations [2208.14583]. Under this usage, the radial perspective is not a separate transform with a single canonical definition, but a structured way of indexing and exploiting persistence information by angular or center-based parameters.

## 1. Conceptual scope and relation to established transforms

The underlying framework is the Persistent Homology Transform (PHT) and its augmented and extended variants. For a simplicial complex \(K\) embedded in \(\mathbb{R}^d\), the PHT/APHT assigns to each direction \(s \in \mathbb{S}^{d-1}\) the persistence data of the lower-star filtration induced by the height function in that direction [2212.13206]. In the notation of the graph reconstruction paper, the augmented transform is
\[
X = \{\,(\dgm{}{s}, s)\,\}_{s\in \mathbb{S}^{d-1}}.
\]
A faithful discretization is defined there as “a finite subset of \(X\) from which all other elements of \(X\) can be deduced” [2212.13206].

The radial interpretation arises because the graph reconstruction algorithm does not merely sample arbitrary directions. It organizes directions by angular position in a fixed projection plane, radially around each vertex, and then uses those directions to localize adjacency information [2212.13206]. The paper states that the improvement in edge reconstruction comes from a “radial binary (multi-)search” exploiting the fact that graph edges can be ordered radially with respect to a reference plane [2212.13206]. This suggests a graph-focused radial version of APHT rather than a distinct transform with new invariance properties.

In the XPHT setting, the indexing parameter remains directional: each unit vector is assigned the extended persistence module of the corresponding height function [2208.14583]. However, the same paper states that similar results could hold for “other kinds of functions, such as radial functions” [2208.14583]. A plausible implication is that a radial persistence transform can be understood as a direct analogue of XPHT in which directions are replaced by centers or radial parameters, while the persistence object is extended persistence rather than ordinary persistence.

## 2. Persistent-homological foundations

For a simplicial complex \(K\) and a filter function \(f\), the sublevel sets
\[
K_t := \{ \sigma \in K \mid f(\sigma) \le t\}, \quad t \in \mathbb{R},
\]
form an increasing filtration, and persistent homology records the birth and death of homological features as \(t\) varies [2212.13206]. The resulting persistence diagram is a multiset of birth–death pairs \((b,d)\) [2212.13206].

The graph reconstruction paper uses the augmented persistence diagram (APD), in which every simplex is assigned to exactly one birth–death pair and some pairs may lie on the diagonal \(b=d\) [2212.13206]. For graphs, the augmentation has two consequences stated explicitly in the paper: every vertex appears in the diagram, and every edge corresponds either to a \(0\)-dimensional death or to a \(1\)-dimensional birth [2212.13206]. This simplex-level encoding is what allows the algorithm to count edges and compute indegrees from persistence data.

For a graph \((V,E)\) GP-immersed in \(\mathbb{R}^d\), and a direction \(s\in \mathbb{S}^{d-1}\), the lower-star filtration is defined by
\[
f_s(v) = s\cdot v
\]
for vertices \(v \in V\), and
\[
f_s([v_0,v_1]) = \max\{s\cdot v_0, s\cdot v_1\}
\]
for edges \([v_0,v_1] \in E\) [2212.13206]. The paper denotes the \(i\)-dimensional APD by \(\dgm{i}{s}\) and the union over all dimensions by
\[
\dgm{}{s} = \bigsqcup_i \dgm{i}{s}.
\]

The XPHT paper enlarges this directional paradigm from ordinary persistence to extended persistence. For a bounded function \(f : M \to \mathbb{R}\), it defines the extended persistence module over the parameter space
\[
\Theta = O \cup R,
\]
where
\[
O = \{(t,\mathrm{Ord}) : t \in \mathbb{R}\}, \qquad
R = \{(t,\mathrm{Rel}) : t \in \mathbb{R}\},
\]
and the order is
\[
\begin{aligned}
(s,\mathrm{Ord}) &< (t,\mathrm{Ord}) &&\text{iff } s < t,\\
(s,\mathrm{Rel}) &< (t,\mathrm{Rel}) &&\text{iff } s > t,\\
(s,\mathrm{Ord}) &< (t,\mathrm{Rel}) &&\text{for all } s,t.
\end{aligned}
\]
The module assigns
\[
V_{(t,\mathrm{Ord})} = H_k(M_t,\emptyset), \qquad
V_{(t,\mathrm{Rel})} = H_k(M, M^t),
\]
with
\[
M_t = f^{-1}(-\infty,t], \qquad M^t = f^{-1}[t,\infty)
\]
[2208.14583]. Extended persistence differs from ordinary persistence in that essential classes become finite intervals spanning the ordinary and relative parts of the filtration, which yields finite Wasserstein distances even when two shapes have different Betti numbers [2208.14583].

## 3. Radial organization in graph reconstruction

The graph reconstruction paper improves a prior finite-discretization result for APHT on graphs by replacing a linear scan over candidate edges with a radial binary multi-search [2212.13206]. Previous work by Belton et al. reconstructed a GP-immersed graph with \(n=|V|\) using \(n^2-n+d+1\) augmented diagrams in
\[
O\big(dn^{d+1} + n^4 + (d+n^2)\,T_{\text{oracle}}\big)
\]
time, with edge reconstruction as the bottleneck because each pair of vertices is tested for adjacency [2212.13206]. The new method reduces edge reconstruction to \(O(m\log n)\) diagrams, where \(m=|E|\), rather than \(O(n^2)\), and the time improves accordingly from the previous \(O(n^4)\) bottleneck [2212.13206].

The radial step begins by fixing the orthogonal projection \(\pi:\mathbb{R}^d\to\mathbb{R}^2\) onto the \((e_1,e_2)\)-plane and ordering projected vertices around a chosen center vertex \(v\) by angle around \(\pi(v)\) [2212.13206]. The paper’s general position assumption requires: every set of \(d+1\) points is affinely independent; no three points are colinear after projection into the \((e_1,e_2)\)-plane; and every point has a unique height with respect to \(e_2\) [2212.13206]. These conditions guarantee that each vertex has a well-defined cyclic order of all other vertices around it in the projection plane and that indegree computations are unambiguous.

The central counting quantity is the indegree of a vertex in a direction. The paper defines
\[
\indeg{v}{s}
\]
as the number of edges incident to \(v\) with height \(s\cdot v\) [2212.13206]. It then proves an edge-counting lemma: for a filter function \(f\), the edges with function value \(c\) are in one-to-one correspondence with
\[
\big\{(b,d)\in \dgm{1}{f} : b=c\big\}\ \cup\ \big\{(b,d)\in \dgm{0}{f} : d=c\big\},
\]
and if \(c=s\cdot v\) is unique among vertices, the cardinality of this multiset is exactly \(\indeg{v}{s}\) [2212.13206]. The paper further states that if \(v\) has unique height in direction \(s\), then \(\indeg{v}{s}\) can be computed from \(\dgm{}{s}\) in \(\zeroindegTimeFull\) time, using \(\zeroindegDGMComp\) diagrams [2212.13206].

This is the mechanism by which radial localization becomes possible. By choosing a direction perpendicular to a selected angular boundary, one obtains an indegree count at the height of \(v\); by subtracting contributions from already identified edges outside the target sector, one recovers the number of edges from \(v\) into that angular region [2212.13206].

## 4. Edge arcs, radial binary multi-search, and sweep reconstruction

To formalize angular sectors, the paper defines an edge arc \(\eA\) centered at \(v\) with start and stop angles \((\alpha_1,\alpha_2)\), a radially ordered vertex array \(\eA.\vV\), and an edge count \(\eA.\eC\) equal to the number of edges between \(v\) and vertices in \(\eA.\vV\) [2212.13206]. Geometrically, \(\eA\) represents a sector in the upper half-plane above \(v\) between angles \(\alpha_1\) and \(\alpha_2\) [2212.13206].

For such an arc, the paper proves an arc-count lemma. Let \(s\in \mathbb{S}^{d-1}\) be the direction perpendicular to \(\alpha_2\) so that the arc is entirely below \(s\cdot v\), and let \(E_*\) be the set of edges with height \(s\cdot v\) that are not in \(\eA\). If no other vertex has height \(s\cdot v\), then
\[
\eA.\eC = \indeg{v}{s} - |E_*|
\]
[2212.13206]. This is the key relation that converts persistence information into an exact edge count for an angular sector.

The binary-search step is implemented by Algorithm 3.1 in the paper, \(\splitArc{\eA}{\varbigedges}{\theta}\), which splits an edge arc into two subarcs of approximately equal cardinality in the candidate vertex list [2212.13206]. The split is chosen by selecting an angle between consecutive rays, using the minimal angular separation \(\theta\), and setting a perpendicular direction
\[
s := e^{i(\alpha - \pi/2)}
\]
in the \((e_1,e_2)\)-plane [2212.13206]. Theorem 3.2 states that the algorithm uses \(\splitDGMComp\) diagrams and time \(\splitTimeFull\), while producing left and right subarcs whose vertex lists partition the original list in clockwise order and satisfy
\[
|\eA_\ell.\vV| = \lceil |\eA.\vV|/2\rceil, \qquad
|\eA_r.\vV| = \lfloor |\eA.\vV|/2\rfloor
\]
[2212.13206].

Outgoing edges of a vertex \(v\) are then found by Algorithm 3.2, \(\findUpEdges{v}{V_v}{in_v}{\theta}{\adgm}\) [2212.13206]. The paper describes the procedure as follows. One computes
\[
\varindeg := \indeg{v}{-e_2},
\]
the number of edges from vertices above \(v\), initializes a stack with a single edge arc covering the entire upper half-plane, and repeatedly applies three rules: discard arcs with \(\eA.\eC=0\); if \(\eA.\eC = |\eA.\vV|\), declare every candidate vertex in the arc adjacent to \(v\); otherwise split the arc and continue [2212.13206]. Theorem 3.3 states that this finds all outgoing edges of \(v\), in clockwise order, using \(\upDGMCompFull\) diagrams and \(\upTimeFull\) time [2212.13206]. Because each split halves the candidate set up to rounding, each recovered outgoing edge requires \(O(\log n)\) angular bisections [2212.13206].

Full reconstruction is obtained by sweeping vertices in increasing \(e_2\)-height. Algorithm 3.3, \(\findEdges{V}\), first obtains the APD in direction \(-e_2\), precomputes clockwise orderings of vertices above each vertex, computes the minimal angular separation \(\theta\), and then processes vertices bottom-to-top, using already reconstructed edges from below as the incoming set \(in_v\) when recovering outgoing edges [2212.13206]. Theorem 3.4 states that \(\findEdges\) reconstructs \(E\) using \(\edgeDGMCompFull\) diagrams and \(\edgeTimeFull\) time [2212.13206]. The paper summarizes the improvement qualitatively as fewer diagrams, output sensitivity with dependence on the actual number of edges \(m\), and dimension-friendly behavior because the radial structure lives in a two-dimensional projection [2212.13206].

## 5. Extended persistence and radial-function analogues

The XPHT paper defines the Extended Persistent Homology Transform by assigning to each direction \(v \in S^{n-1}\) the tuple
\[
XPHT(M)(v) = \big( XPH_0(M,h_v^M), XPH_1(M,h_v^M), \ldots, XPH_{n-1}(M,h_v^M) \big)
\]
[2208.14583]. It then defines a distance between shapes by integrating Wasserstein distances between their extended persistence modules:
\[
d^{XPHT}_p(M_1,M_2)^p = \int_{v\in S^{n-1}} \sum_{k=0}^{n-1} W_p\big( XPH_k(M_1,h_v), XPH_k(M_2,h_v) \big)^p \, dv
\]
[2208.14583]. Because extended persistence turns essential classes into finite intervals, this distance remains finite even when shapes have different Betti numbers [2208.14583].

The paper’s principal theoretical contribution for manifolds with boundary is that the extended persistence of a manifold \(A\) for a height function can be deduced from the extended persistence of the boundary \(X=\partial A\), together with the sign labels of boundary critical points [2208.14583]. Theorem 3.2 states that if
\[
XPH_k(X,h_v) = \bigoplus_{[b_i,d_i)\in S_X} \mathcal{I}_{[b_i,d_i)},
\]
then
\[
XPH_k(A,h_v) = \bigoplus_{[b_i,d_i)\in J_A^k} \mathcal{I}_{[b_i,d_i)},
\]
where \(J_A^k\) selects intervals according to whether their births correspond to positive or negative critical points of \(h_v^A\) with the appropriate Morse indices [2208.14583]. The paper further gives explicit formulas for essential classes in terms of minima and maxima of the height function on connected boundary components [2208.14583].

These results are directional rather than radial, but the paper explicitly states that “it is reasonable to expect that similar results could hold for other kinds of functions, such as radial functions” [2208.14583]. It also presents a consistent radial substitution in notation:
\[
r_c(x) = \|x-c\| \quad\text{or}\quad r_c(x)=\|x-c\|^2
\]
as the candidate radial function centered at \(c\) [2208.14583]. This suggests a radial extended persistent homology transform in which the parameter space is a family of centers \(c\) rather than a sphere of directions, and the persistence modules are those of radial filtrations instead of linear height filtrations. The paper even writes the corresponding transform and distance formulas as a natural extension:
\[
RPT(M)(c) = \big( XPH_0(M,r_c),\dots,XPH_{n-1}(M,r_c) \big),
\]
and
\[
d^{REXPHT}_p(M_1,M_2)^p = \int_{\mathcal{C}} \sum_{k=0}^{n-1} W_p\big(XPH_k(M_1,r_c),XPH_k(M_2,r_c)\big)^p\,dc
\]
[2208.14583]. Since the paper does not formalize these as theorems, they are best read as a technically grounded extension rather than an established definition.

## 6. Assumptions, computational practice, and limitations

Both lines of work rely on strong regularity assumptions. In the graph case, correctness depends on the general position conditions in Assumption 2.1: affine independence, no three projected vertices colinear, and unique heights in the reference direction \(e_2\) [2212.13206]. The paper notes that a version of the unique-height condition can be enforced algorithmically by basis tilting [2212.13206]. It also assumes exact augmented persistence diagrams returned by an oracle, and its implementation model uses efficient lookup in the diagram by birth and death values [2212.13206]. The method is tailored to graphs and one-skeleta; the paper explicitly states that “Radially ordering higher dimensional simplices is not well-defined, and this issue prevents the methods presented here from being immediately transferrable” [2212.13206].

The XPHT framework likewise assumes Morse-theoretic regularity for the function under study, either in smooth or PL form [2208.14583]. For binary images, the paper builds a PL manifold with boundary from 8-adjacent foreground and 4-adjacent background pixels, proves that the resulting boundary is a disjoint union of simple closed PL curves, and then computes \(PH_0\) of boundary filtrations by union–find [2208.14583]. It handles non-generic height ties by identifying \(0\)-critical vertices and flat critical segments, and classifies positivity of critical points via explicit geometric tests involving determinants and rotated edge vectors [2208.14583]. Directions are sampled as
\[
v_i = (\cos\frac{2\pi i}{K}, \sin\frac{2\pi i}{K}), \quad i=0,\dots,K-1,
\]
and antipodal duality reduces the number of required height filtrations by a factor of two when \(K\) is even [2208.14583].

The principal limitation of the graph-based radial method is scope: it reconstructs one-skeleta, not higher-dimensional simplicial structure [2212.13206]. The principal limitation of the radial extension suggested by XPHT is that it remains a proposal rather than a completed theory in the cited paper [2208.14583]. The paper states that the Morse-theoretic and boundary-based decomposition principles “are not specific to linear functions” and that radial functions are a natural candidate [2208.14583]. This suggests that a fully developed radial persistence transform would require proofs of injectivity, stability, and computational reductions comparable to those already established for PHT, APHT, and XPHT.

Taken together, the two papers locate the radial persistence idea at the intersection of directional topological transforms and structured parameterization. In one case, the radial structure is explicit and algorithmic: directions are organized by angle around vertices to enable output-sensitive graph reconstruction from APHT data [2212.13206]. In the other, the radial structure is prospective but theoretically motivated: extended persistence for radial functions is presented as a natural generalization of the height-based XPHT framework [2208.14583].

Source: https://www.emergentmind.com/topics/radial-persistence-transform