---
title: Persistent Homology Convolutions
url: https://www.emergentmind.com/topics/persistent-homology-convolutions
type: topic
---

# Persistent Homology Convolutions

Persistent homology convolutions designate several constructions that combine convolutional operators with persistent homology, but the locus of convolution differs across the literature. Depending on context, convolution is applied to the input signal before persistence is computed, to overlapping local windows through persistent homology-based feature extraction, to lattice-indexed invariants of multiparameter persistence modules through meet/join operations, to PH-derived real-valued functions on input or parameter space as a stabilizing average, or to sheaf-theoretic realizations of persistence objects as a thickening operator [2208.02107] [2507.14378] [2011.14057] [1512.01700] [1705.00955]. Across these formulations, the common objective is to make topological information compatible with locality, translation, order structure, or stability requirements that ordinary persistent homology or ordinary convolution alone do not directly address.

## 1. Scope and terminological landscape

The literature does not use a single standardized meaning for the phrase. A common misconception is that persistent homology convolutions always mean convolving persistence diagrams themselves. The principal lines of work summarized here do not adopt that interpretation. Instead, they treat convolution as motif detection before persistence, as a local topological receptive-field mechanism, as an order-theoretic operator on multiparameter persistence data, as a smoothing device on unstable PH-derived outputs, or as a categorical operation in sheaf-theoretic persistence [2208.02107] [2507.14378] [2011.14057] [1512.01700] [1705.00955].

| Construction | Role of convolution | Representative source |
|---|---|---|
| Convolutional persistence | Convolve data with filters before computing persistence | [2208.02107] |
| Persistent Homology Convolutions | Compute PH on local windows and combine via a kernel | [2507.14378] |
| Lattice-theoretic convolutions | Replace translations with meet/join on a finite lattice | [2011.14057] |
| Stabilizing unstable PH output | Convolve a PH-derived function on input/parameter space | [1512.01700] |
| Microlocal sheaf convolution | Use \(F\star G = Rs_!(F\boxtimes G)\) as a thickening operator | [1705.00955] |

Two recurrent themes organize these approaches. First, several papers argue that ordinary persistence can be too coarse when the relevant information is local or motif-dependent. Second, several papers argue that ordinary convolution may be algebraically mismatched to persistence data, either because topology should be extracted locally rather than globally, or because multiparameter persistence is indexed by a poset or lattice rather than by a translation group.

## 2. Convolutional persistence as motif detection

In "Convolutional Persistence Transforms" [2208.02107], convolution is applied to data defined over simplicial or cubical complexes before persistent homology is computed. The central idea is to interpret a filter as a local motif. For an image \(f:P\to \mathbb{R}\) on a rectangle \(P \subset \mathbb{Z}^d\), with filter \(g:B\to \mathbb{R}\) on another rectangle \(B\subset \mathbb{Z}^d\) and stride vector \(k\), the convolution is
\[
(f \ast g)(v) = \sum_{p \in B} g(p)f(p+v \odot k).
\]
For multi-channel data \(f:P\to \mathbb{R}^n\), \(g:B\to \mathbb{R}^n\), the convolution becomes
\[
(f \ast g)(v) = \sum_{p \in B} g(p)\cdot f(p+v \odot k).
\]
The scalar response field is then extended by the lower-\(*\) rule to a cubical complex on the output grid, and persistence is computed on the resulting filtration.

The simplicial analogue replaces image patches by a vertex signal on a connected simplicial complex. If \(X\in \mathbb{R}^{n\times d}\) encodes vertex features, \(A\) is an arbitrary fixed \(n\times n\) matrix, and \(W=w\) is a \(d\times 1\) weight vector, then the convolved signal is the scalar vertex function associated to \(AXw\). Persistent homology is computed on the lower-\(*\) extension of that scalar function over the original simplicial complex. This produces a transform viewpoint: for a family of filters, one obtains a family of persistence diagrams indexed by filters.

The paper’s conceptual claim is that the persistence diagram of the resulting convolution describes the way the motif is distributed across the simplicial complex. This is not merely smoothing. The convolved scalar field is a motif-response field, and persistent homology measures the topology of regions where the motif response is strong or weak. The trivial \(1\times 1\) filter recovers ordinary persistence, so convolutional persistence strictly generalizes standard persistence.

The theoretical contribution is unusually strong. The paper proves that, generically speaking, for any two labeled complexes one can find some filter for which they produce different persistence diagrams, so that the collection of all possible convolutional persistence diagrams is an injective invariant. This is established by showing convolutional persistence to be a special case of the Persistent Homology Transform. The paper also proves explicit stability bounds; for example, in the image setting,
\[
W_{\infty}(\operatorname{Diag}(\widehat{f_1 \ast g}),\operatorname{Diag}(\widehat{f_2 \ast g}))
\leq \|g\|_{1} \|f_1 - f_2\|_{\infty}.
\]
A further practical point is computational: with stride \(k=(2,\dots,2)\), persistence complexity becomes
\[
O\!\left(\left(\frac{p}{2}\right)^{d\omega}\right),
\]
and more generally, for \(M\) filters and stride-product \(\kappa=\Pi_i k_i\), the cost is
\[
O\!\left(M\left(\frac{|P|}{\kappa}\right)^{\omega}\right).
\]

The experiments use UCI digits, MNIST, Chinese handwritten digits, Devanagari characters, and a Kuramoto–Sivashinsky PDE parameter-estimation task. The reported empirical message is that convolutional persistence substantially outperforms ordinary persistence, even when one uses random filters and vectorizes the resulting diagrams by recording only their total persistences [2208.02107].

## 3. Local persistent homology convolutions for histopathology

"Classification of Histopathology Slides with Persistence Homology Convolutions" [2507.14378] introduces Persistent Homology Convolutions (PHCs) as a convolution-like operator that replaces local intensity-based feature extraction with local topological feature extraction. The motivation is domain-specific: in histopathology, topology is an important descriptor because disease-relevant abnormalities include cell and nucleus size variation, multinucleation, disorganization of tissue, and changes in the arrangement of cellular voids. The paper argues that previous methods used global topological summaries computed over an entire image, whereas PHC is designed to retain both locality and the convolutional notion of translation invariance of topological features.

The formalism begins with an \(N\times N\) grayscale image \(X\in\mathcal{G}_N\), where
\[
\mathcal{G}_N=[0,1]^{[[N]]\times[[N]]}, \qquad [[N]]=\{0,1,\dots,N\}.
\]
For \(M<N\), the map
\[
T_{(x,y)}:\mathcal{G}_N\to\mathcal{G}_M
\]
translates the \((x,y)\)-coordinate of \(X\) to the origin and restricts to the \(M\times M\) window \( [[M]]^2 \). PHC then slides a local window over the image, builds a filtration on each local subimage, computes persistent homology in that receptive field, vectorizes the resulting persistence diagram, and combines the local topological vectors via a kernel. In the implementation emphasized in the paper, \(\mathcal{F}(T_{(c i,c j)}(X))\) is the extended sublevel-set filtration of the height function on the \(M\times M\) subimage after thresholding; persistence is computed primarily in dimension \(p=1\); and the vectorization is a persistence image.

The local topological pipeline is highly specific. Images are converted from RGB to grayscale, resized from \(1024\times 1024\) to \(512\times 512\), thresholded, and then morphologically conditioned. The text contains an inconsistency: one part refers to “image erosion,” while another describes applying dilation once with a \(2\times 2\) kernel. The thresholded pixels are inserted into a simplex tree, along with edges between neighboring nonzero pixels, producing the adjacency complex of the image. The filtration is extended persistence of the height function \(f\), and the paper restricts attention to
\[
\mathrm{Ext}_1(f),
\]
explicitly excluding summaries generated when \(p=0\). Each local diagram is vectorized as a persistence image \(\mathbf{I}_{n\times n}(z)\), with \(n\in\{10,15,20,25\}\) and \(n=20\) fixed for reported experiments.

In experiments, \(M=32\) and \(c=32\), so windows are \(32\times 32\). For a \(512\times 512\) image, the output is
\[
(256,20,20),
\]
that is, 256 local persistence images, one per spatial locality. The persistent homology computation itself is not learned. The learnable part is the kernel \(k(x,y)\), described as “a linear operator (an arbitrary \(q\times q\) matrix) whose entries are optimized via backpropagation during model training.” The authors explicitly note that most of the PHC data generation is a separate preprocessing step applied on the image data outside of the training cycle except for optimizing the weights of \(k(x,y)\). PHC is therefore not end-to-end differentiable in the reported implementation.

The experimental setting uses the Osteosarcoma histopathology dataset from TCIA: 1144 RGB images, originally \(1024\times 1024\times 3\), with 3 classes—non-tumorous (47%), necrotic tumor (23%), and viable tumor (30%). After balancing by resampling and augmentation, the experimental dataset contains 1143 images. The study compares thresholded image data, global persistent homology, local persistent homology via PHC, images plus global PH, and images plus PHC. For the best accuracy-maximizing models, Height PHC attains accuracy \(0.9170\), precision \(0.9167\), sensitivity \(0.9389\), and specificity \(0.9585\), while Images + Height PHC attains accuracy \(0.9389\), precision \(0.9430\), sensitivity \(0.9651\), and specificity \(0.9716\). Height PHC also compares favorably in runtime: for a sample of 100 images, average per-image computation times were \(11.35\) sec for Height PHC and \(3496.0\) sec for Global Height PH. Across 1000 model initializations, Height PHC yields
\[
0.8480 \pm 0.0412 \text{ accuracy},
\]
compared with
\[
0.8051 \pm 0.0935 \text{ accuracy}
\]
for Image Data and
\[
0.8653 \pm 0.0561 \text{ accuracy}
\]
for Images + Height PHC. The paper interprets these results as evidence that PHC preserves where topological features occur and is less sensitive to hyperparameters [2507.14378].

## 4. Lattice-theoretic convolutions on multiparameter persistence modules

"Multidimensional Persistence Module Classification via Lattice-Theoretic Convolutions" [2011.14057] addresses a different object: not raw images, but multiparameter persistence modules themselves. The starting point is that persistent homology naturally produces data indexed by an ordered parameter set. In the one-parameter case, a filtration such as \(\mathrm{Rips}_r(\mathcal M)\) yields a persistence module
\[
\mathbb R \to \mathrm{Vect}_k,
\]
and the total order on \(\mathbb R\) permits interval decomposition and barcode summaries. In the multiparameter case, the paper considers a bifiltration
\[
\mathbb X_{r,t}=\mathrm{Rips}_r\{x\in \mathcal M\mid \rho(x)\le t\},
\]
which gives a \(2\)-parameter persistence module
\[
PH_i:\mathbb R_+\times \mathbb R\to \mathrm{Vect};\qquad (r,t)\mapsto H_i(X_{r,t}).
\]
Because the indexing category is now the poset \(\mathbb R_+\times \mathbb R\) with product order, there is “no complete compact barcode-like characterization” of such modules.

The paper therefore uses computable incomplete invariants. For a fixed homological degree \(i\), it encodes the module by the Hilbert function
\[
\mathrm{Hilb}:\mathbb R_+\times \mathbb R\to \mathbb Z_+;\qquad (r,t)\mapsto \dim(PH_i(r,t)),
\]
together with the multigraded Betti numbers
\[
\xi_j:\mathbb R_+\times \mathbb R\to \mathbb Z_+,\qquad j=0,1,2.
\]
This yields four channels over a discretized parameter grid: one Hilbert function plus \(\xi_0,\xi_1,\xi_2\). The conceptual claim is that ordinary convolutions on these arrays are not obviously the natural operators, because standard convolution assumes an abelian group such as \(\mathbb Z^2\), while the natural algebraic structure of the indexing set is a lattice.

The proposed alternative is lattice-theoretic signal processing. For signals \(f,g:L\to \mathbb R\) on a lattice \(L\), the paper defines
\[
(f *_\wedge g)(x)=\sum_{a\in L} f(x\wedge a)\, g(a), \qquad (f *_\vee g)(x)=\sum_{a\in L} f(x\vee a)\, g(a).
\]
For the finite lattice
\[
L=[m]\times [n], \qquad [n]=\{0,1,\dots,n\},
\]
meet and join are computed coordinatewise:
\[
(r,t)\wedge (r',t')=(\min(r,r'),\min(t,t')), \qquad
(r,t)\vee (r',t')=(\max(r,r'),\max(t,t')).
\]
The multichannel meet-convolution and join-convolution layers are
\[
\mathrm{MeetConv}(f)(x,y)^j = \sum_i \sum_{(a,b)\in [m]\times [n]} f_i(x\wedge a,\, y\wedge b)\, g^i_j(a,b),
\]
and
\[
\mathrm{JoinConv}(f)(x,y)^j = \sum_i \sum_{(a,b)\in [m]\times [n]} f_i(x\vee a,\, y\vee b)\, g^i_j(a,b).
\]
The actual layer used in experiments combines both directions as
\[
\alpha\cdot \mathrm{MeetConv}(x)+(1-\alpha)\cdot \mathrm{JoinConv}(x), \qquad \alpha\in [0,1],
\]
with
\[
\alpha=\tfrac12.
\]

A key technical issue is receptive-field design. For meet-convolution the appropriate neutral element is the maximum \((m,n)\), while for join-convolution it is the minimum \((0,0)\). The paper remarks that kernels supported only on a small contiguous block near one corner can produce degenerate receptive fields; for join-convolution, many neurons may see only themselves. The authors therefore hypothesize that appropriate lattice kernels should be supported on evenly spaced sublattices including both \((0,0)\) and \((m,n)\), and in experiments they use an evenly spaced \(4\times 4\) grid of lattice points.

The end-to-end pipeline uses a subset of Princeton ModelNet containing 10 object classes. For each CAD model, 3000 vertices are sampled to obtain a point cloud in \(\mathbb R^3\). A bifiltered simplicial complex is built using a geometric scale parameter and the codensity function
\[
\rho_{\mathrm{codense}}(x;k) = \left( \frac{1}{k}\sum_{y\in N_k(x)} d(x,y) \right)^{-1},
\]
with \(k=100\). Using RIVET, the authors compute degree-\(0\) multiparameter persistent homology on a discrete \(40\times 40\) grid, forming a \(40\times 40\times 4\) tensor. The lattice-convolutional CNN and a standard CNN baseline each have three convolutional layers followed by two fully connected layers, with hidden convolutional layers of 16 channels, a final convolutional layer of 8 channels, \(2\times 2\) max-pooling after the first two convolutional layers, a fully connected hidden layer of 32 features, a final softmax with cross-entropy loss, and the Adam optimizer with learning rate
\[
2\times 10^{-4}
\]
for 300 epochs. The main empirical result is modest: the lattice convolutional network “slightly underperforms” the standard convolutional classifier. The paper presents this not as superiority, but as a proof of concept for an order-aware convolutional formalism for multiparameter persistence modules [2011.14057].

## 5. Convolution as stabilization in input or parameter space

"Stabilizing the unstable output of persistent homology computations" [1512.01700] uses the word convolution in yet another sense. Here convolution does not act on images, on persistence modules, or on diagrams. Instead, some PH-derived quantity of interest is viewed as a real-valued function on the space of inputs or parameters, and convolution with a kernel is used to make that quantity Lipschitz-stable.

The formal setup is as follows. Let \(\mathbf{a}\in \mathbb{R}^n\) denote all real parameters entering a persistent homology pipeline \(C\), and let \(H:\mathbb{R}^n \to O\) denote the induced output map. One then chooses a real-valued summary \(p:O\to \mathbb{R}\), such as the persistence of the class born by a specific simplex, an indicator that a representative cycle intersects a chosen region, or the longest persistence after density thresholding. This defines
\[
h = p\circ H : \mathbb{R}^n \to \mathbb{R}.
\]
The function \(h\) may be discontinuous but measurable. The stabilized quantity is
\[
g = h*K,
\]
where
\[
(h*K)(t)=\int_{\mathbb{R}^d} h(s)K(t-s)\,ds
      =\int_{\mathbb{R}^d} h(t-s)K(s)\,ds.
\]
The paper discusses triangular, Epanechnikov, and Gaussian kernels.

The core theorems state that smoothing by convolution yields Lipschitz control under mild assumptions. If
\[
\|h\|_1=a
\]
and \(K\) is \(b\)-Lipschitz, then \(h*K\) is \(ab\)-Lipschitz. If \(\|h\|_\infty \le M\) locally and \(K\) is compactly supported and Lipschitz, then \(h*K\) is locally Lipschitz. If \(\|h\|_\infty \le M\) and \(K\) satisfies an \(L^1\)-translation bound, then \(h*K\) is globally Lipschitz. The conceptual interpretation is that exact localization of a generator or exact identification of a simplex may be unstable, while an average score under small perturbations is stable.

The method is computationally practical because convolution can be estimated by Monte Carlo:
\[
\frac1M \sum_{i=1}^M h(a-\epsilon_i) \to (h*K)(a)
\]
when the \(\epsilon_i\) are sampled i.i.d. from \(K\). The operational prescription is therefore to perturb the input many times, rerun persistent homology, and average the unstable output. The paper illustrates this with localization of representative cycles in a double-annulus example, with a brain-imaging example in which the stabilized probability for the generator of the \(28\)th longest bar to lie in a chosen ball is estimated as
\[
(h*K)(a)\approx 0.637,
\]
and with thresholding-based examples in which the persistence diagram itself is unstable with respect to preprocessing parameters. The paper is explicit that this framework complements stable diagram summaries by stabilizing information beyond the diagram, such as where a generator is located or which simplex created a feature [1512.01700].

## 6. Sheaf-theoretic convolution and derived persistence

"Persistent homology and microlocal sheaf theory" [1705.00955] develops a categorical notion of convolution on the derived category of sheaves on a finite-dimensional real vector space \(V\). This is not a convolution on persistence modules in the usual algebraic sense. Rather, it is a sheaf-theoretic convolution used to reinterpret persistence, stability, and higher-dimensional barcode-like structures.

The ambient category is
\[
D^b(\mathbf{k}_V),
\]
and the convolution functor is
\[
F\star G \eqdot Rs_!(F\boxtimes G),
\]
where \(s:V\times V \to V\) is addition. The category becomes a commutative monoidal category with unit \(\mathbf{k}_{\{0\}}\). A family of kernels \(K_a\) satisfies
\[
K_a\star K_b \simeq K_{a+b},
\]
so convolution with \(K_a\) behaves like a metric thickening operator. This is the basis for a pseudo-distance.

For \(F,G\in D^b(\mathbf{k}_V)\) and \(a\ge 0\), the paper says that \(F\) and \(G\) are \(a\)-isomorphic if there exist morphisms
\[
f:K_a\star F \to G,\qquad g:K_a\star G\to F
\]
satisfying compatibility conditions through \(K_{2a}\). The pseudo-distance is
\[
\dist(F,G)=\inf\Bigl(\{+\infty\}\cup\{a\in\mathbb R_{\ge0};\ \text{$F$ and $G$ are $a$-isomorphic}\}\Bigr).
\]
This is explicitly interleaving-like, but formulated in the derived sheaf category and expressed through convolution kernels \(K_a\) rather than shifts of persistence modules. The distance is symmetric, satisfies the triangle inequality, and can be finite even for sheaves concentrated in different cohomological degrees.

The main stability theorem states that for a locally compact space \(X\), continuous maps \(f_1,f_2:X\to V\), and any \(F\in D^b(\mathbf{k}_X)\),
\[
\dist(Rf_{1*}F,Rf_{2*}F)\le \|f_1-f_2\|, \qquad \dist(Rf_{1!}F,Rf_{2!}F)\le \|f_1-f_2\|.
\]
This is a sheaf-theoretic analogue of classical persistence stability. The broader microlocal framework introduces \(\gamma\)-sheaves, that is, sheaves on \(V\) whose microsupport lies in a cone determined by a closed convex proper cone \(\gamma\) with nonempty interior. The paper proves an equivalence between \(\gamma\)-sheaves and sheaves on the \(\gamma\)-topology \(V_\gamma\), approximation of constructible sheaves by piecewise-linear objects in the convolution pseudo-distance, and a higher-dimensional barcode analogue in which piecewise-linear \(\gamma\)-sheaves are constant on stratifications by \(\gamma\)-locally closed sets.

In dimension one, the framework recovers ordinary barcode theory. With \(V=\mathbb R\) and \(\gamma=\mathbb R_{\le 0}\), \(\gamma\)-locally closed sets are intervals \([a,b)\), and any \(\mathbb R\)-constructible \(\gamma\)-sheaf decomposes uniquely as a locally finite direct sum of interval sheaves. In higher dimensions, the paper does not obtain a full barcode classification. Instead, it establishes a richer categorical structure in which higher-dimensional “barcodes” embed fully faithfully but are not essentially surjective. This suggests that convolution, in this setting, functions as a conceptual higher-dimensional and categorical analogue of thickening and interleaving in persistence [1705.00955].

Source: https://www.emergentmind.com/topics/persistent-homology-convolutions