Topological Signal Processing from Stereo Visual SLAM

TL;DR

This paper introduces TSP-SLAM, a new method that enhances point cloud processing in visual SLAM by using topological signal processing techniques.

Contribution

The novel TSP-SLAM framework integrates texture information from V-SLAM into topological signal processing.

Findings

TSP-SLAM allows signals to be associated with vertices, edges, and faces of a 3D mesh.

Numerical results show TSP-SLAM supports advanced topological filtering algorithms.

The method improves point cloud representation in challenging V-SLAM environments.

Abstract

Topological signal processing is emerging alongside Graph Signal Processing (GSP) in various applications, incorporating higher-order connectivity structures—such as faces—in addition to nodes and edges, for enriched connectivity modeling. Rich point clouds acquired by multi-camera systems in Visual Simultaneous Localization and Mapping (V-SLAM) are typically processed using graph-based methods. In this work, we introduce a topological signal processing (TSP) framework that integrates texture information extracted from V-SLAM; we refer to this framework as TSP-SLAM. We show how TSP-SLAM enables the extension of graph-based point cloud processing to more advanced topological signal processing techniques. We demonstrate, on real stereo data, that TSP-SLAM enables a richer point cloud representation by associating signals not only with vertices but also with edges and faces of the mesh computed from the point cloud. Numerical results show that TSP-SLAM supports the design of topological filtering algorithms by exploiting the mapping between the 3D mesh faces, edges and vertices and their 2D image projections. These findings confirm the potential of TSP-SLAM for topological signal processing of point cloud data acquired in challenging V-SLAM environments.

Figures

Full text

1. Introduction

Topological signal processing [1] is gaining momentum in different applications ranging from the characterization of brain networks [2] to the monitoring of physical systems [3,4]. Compared to GSP, TSP extends the representation by associating signals not only with nodes and edges, but also with higher-order connectivity structures such as faces, thereby enabling enriched connectivity modeling. Representing signals on topological domains remains an ongoing issue [5], particularly in the context of dictionary learning.

In this work, we propose a novel framework, termed TSP-SLAM, which establishes a direct link between the texture information provided by Visual Simultaneous Localization and Mapping (V-SLAM) and the topological signal processing of point clouds acquired through multi-camera systems. While V-SLAM has achieved significant advances for 3D reconstruction in static [6] and dynamic [7,8,9] scenarios, graph-based methods relying solely on node-to-node connectivity fail to capture higher-order structures such as planar surfaces and occlusions [10]. Higher-order topological representations, including faces and simplicial complexes, improve robustness and semantic consistency, and recent studies demonstrate that both clustering and geometric learning benefit from such enriched connectivity [10,11,12].

State-of-the-art extensions, such as the integration of deep learning [13] and hardware optimizations [14], further enhance robustness by adaptively handling visual features under varying conditions. SLAM outputs are increasingly enriched with value-aware geometric reconstructions, as in [15] where a dense neural point cloud model encodes attributes such as normals and intensity. Implicit neural surfaces [16] enable signal extraction on dense geometries, while Red Green and Blue-Depth (RGB-D) and Light Detection and Ranging (LiDAR) information in topology-based models [17] support loop closure and relocalization and enhance robustness and scalability. Graph-theoretic signatures (e.g., von Neumann entropy, spanning trees) [18] reduce uncertainty while lowering computational cost. With the growing adoption of real-time navigation and mapping for autonomous systems [19,20], V-SLAM enables on-the-fly mapping and safe adaptation to environmental changes [21] simultaneously performing localization and mapping using camera input [22,23] in dynamic autonomous system applications.

Within this context, the main contribution of this study is the introduction of TSP-SLAM, a texture-aware topological framework that extends the 2D-to-3D geometric mappings of V-SLAM to associate signals not only with individual points in the cloud, but also with higher-order topological structures. In addition, recent work such as [24] has demonstrated the use of topological descriptors for segmentation and recognition of objects from point clouds. While effective, these methods primarily rely on geometric information and do not exploit luminance or photometric texture cues. Our framework suggests that integrating luminance-based signals alongside geometric descriptors could further benefit tasks such as segmentation and recognition in point cloud processing. Specifically, TSP-SLAM supports the construction of signals over nodes, edges, and faces, thereby enabling advanced topological signal processing strategies. As a case study, we present Topological Multiscale Anisotropic Harmonic Filtering (T-MAHF), which extends to the topological domain a filtering technique (MAHF) originally introduced on graphs and applied successfully to denoising problems, thereby enhancing localization and mapping accuracy [25]. Numerical results confirm the effectiveness of the enriched representation achieved by TSP-SLAM for robust and expressive point cloud processing.

The remainder of this paper is organized as follows. Section 2 introduces the TSP-SLAM framework, providing a detailed description of the construction of topological structures from V-SLAM point clouds and the assignment of signals to nodes, edges, and faces. In Section 3, we present the mathematical formulation and implementation aspects of the T-MAHF algorithm. Section 4 reports the results of numerical experiments and comparative analyses, demonstrating the advantages of TSP-SLAM in terms of expressiveness of the resulting signal representations. Finally, Section 5 concludes the paper by summarizing the main contributions and outlining potential directions for future research.

2. Topological Signal Processing Framework from Visual SLAM

In this section, we introduce TSP-SLAM, a topological signal processing framework [1] enriched with texture information extracted from visual SLAM.

2.1. Signal Model in Visual-SLAM

Herein, we review the system outlined in Figure 1, reflecting a commonly adopted acquisition and tracking architecture as in [7]. The acquisition setup is stereo-based and relies on two sensors: a left camera, taken as reference, and a right camera, employed for 3D reconstruction. The left and right camera images are obtained by sampling the ideal bidimensional signals

$$[eqn]$$

on a discrete $[eqn]$ grid. In V-SLAM, localization and mapping techniques are used to reconstruct a 3D point cloud from the stereo image disparity map, and the point features are used for object tracking and classification purposes. After stereo acquisition using the two calibrated cameras (left and right), which capture synchronized images, a set of distinctive keypoints is extracted from both images using feature extraction techniques such as ORB, SIFT, or BRISK. Let us denote these 2D points as $[eqn]$ for the left image and $[eqn]$ for the right image, with $[eqn]$ , where i indexes the detected feature points in the left and right images. Feature matching algorithm, implemented through Hamming distance, is used to associate points between the two images, allowing the computation of the disparity for each match as $[eqn]$ . In homogeneous coordinates—where affine transformations can be represented as matrix multiplications—each point cloud vertex is represented by augmenting its 3D coordinates with a unitary component. The pin-hole model describes the relationship between each three-dimensional point $[eqn]$ and its projection onto the images, depending on the cameras’ intrinsic parameters such as optical center, focal length, and optical axes. For a point cloud vertex at depth z from the camera axis origin, taking into account the camera focal length l, the relation between the 3D coordinates $[eqn]$ and their 2D counterparts $[eqn]$ in the image domain are written as follows:

$$[eqn]$$

where the perspective projection matrix $[eqn]$ accounts for the origin of the image coordinate system at the principal point $[eqn]$ , the pixel dimensions $[eqn]$ , and the angle between the axes $[eqn]$ arising from manufacturing imperfections. The outcome of this process is a 3D point cloud of vertices $[eqn]$ the vertices of the point cloud. A real-valued attribute signal is associated with each vertex, representing the color components of the points, or other attributes inherited from the original imaging system (e.g., temperature, reflectance, etc.).

In the following, we assume that the signal associated with the i-th point cloud vertex is the luminance of the corresponding 2D point as captured by the reference (left) camera. The point cloud can be further structured into a surface mesh using Delaunay triangulation, which ensures geometric consistency and avoids poorly shaped (e.g., degenerate or sliver) triangles. Based on these positions, we can introduce the topological signal processing framework as follows. The notation is summarized in Table 1.

2.2. Topological Signal Processing from V-SLAM Data

Let us define the graph $[eqn]$ associated with the point cloud as $[eqn]$ , where $[eqn]$ is the set of N point cloud vertices $[eqn]$ , with $[eqn]$ , and $[eqn]$ is the set of edges. A topologically enriched graph structure is described by the vertex-to-edge incidence matrix $[eqn]$ and the edge-to-face incidence matrix $[eqn]$ . Furthermore, the vertex Laplacian $[eqn]$ , and the edge Laplacian $[eqn]$ act as discrete differential operators on signals defined on the topological space.

As a toy example, we illustrate these concepts in Figure 2, which represents a graph $[eqn]$ with $[eqn]$ nodes, $[eqn]$ edges, and $[eqn]$ faces. The $[eqn]$ vertex-to-edge incidence matrix $[eqn]$ and the $[eqn]$ edge-to-face incidence matrix $[eqn]$ are defined as follows:

$$[eqn]$$

In topological signal processing, signals are assigned to different topological elements of the graph, i.e., nodes, edges, and faces, enabling topological signal processing tasks like smoothing or denoising in the edge and face domains. For instance, with reference to the toy example in Figure 2, overall, $[eqn]$ signal values can be defined on the node, edge, and face elements.

Let us denote by $[eqn]$ the signal at i-th node, by $[eqn]$ the signal at j-th graph edge, and by $[eqn]$ the signal at k-th graph face.

In general, at the node level, signals $[eqn]$ may represent sensor data directly acquired at each point vertex, e.g., RGB color values. At the edge level, signals $[eqn]$ may encode pairwise relationships between adjacent vertices, e.g., geometric descriptors such as their Euclidean distance or the differences in vertex attributes. At the face level, signals $[eqn]$ can capture higher-order geometric features, such as surface normals or texture patterns extracted from the original image.

In TSP-SLAM, we propose that the signals at different levels are obtained from the texture side information available through the V-SLAM framework. In detail, the signal at the nodes $[eqn]$ , edges $[eqn]$ and faces $[eqn]$ are derived from the photometric values pertaining to their projections, i.e., the points $[eqn]$ , segments $[eqn]$ or triangles $[eqn]$ lying in the original reference image used by the V-SLAM technique.

The way values on edges $[eqn]$ and faces $[eqn]$ are determined from photometric information associated with their projected segments, $[eqn]$ and $[eqn]$ , is not unique. Each topological element corresponds to a set of pixels in the 2D domain: the j-th edge $[eqn]$ corresponds to the set $[eqn]$ of pixels along the segment $[eqn]$ , while the k-th face $[eqn]$ corresponds to the set $[eqn]$ of pixels within the triangle $[eqn]$ . A reasonable approach is to associate each topological element $[eqn]$ and face $[eqn]$ with a signal obtained by averaging the photometric values of the corresponding set of pixels $[eqn]$ , $[eqn]$ in the image plane. Without loss of generality, we adopt this method in the following, where we show how TSP-SLAM enables topological signal processing on V-SLAM-acquired point clouds, assigning distinct signal information to vertices, edges, and faces. Nevertheless, alternative strategies for defining vector-valued signals [1] on topological elements can be designed, and their exploration is left for future work.

In the following, without loss of generality, to define the signal related to edge and faces, we exploit the information available from the TSP-SLAM framework. Specifically, we define the signal $[eqn]$ at the j-th vertex as the signal of its 2D projection

$$[eqn]$$

the signal $[eqn]$ at the j-th edge as the average of the signal over the set $[eqn]$ of pixels belonging to the segment $[eqn]$ ,

$$[eqn]$$

the signal $[eqn]$ at the k-th face $[eqn]$ as the average on the set $[eqn]$ of pixels within the triangle $[eqn]$

$$[eqn]$$

3. Topological Harmonic Filtering in TSP-SLAM

Here, we show an application of TSP-SLAM to extend anisotropic graph filtering to topological spaces. Specifically, we extend a filtering technique previously introduced for images and signals on graphs to the topological space, and we show how it can leverage the side information provided by TSP-SLAM. Thanks to the information acquired by TSP-SLAM, each topological level—vertex, edge, face—carries its own type of signal, allowing for a rich and flexible representation of the geometry and semantics of the point cloud.

Multiscale harmonic filters were formerly introduced and widely applied in image processing. Among these, Complex Harmonic Filters [26] operate over a 2D spatial domain as a function of radial distance and angular orientation. The filtering is written as follows:

$$[eqn]$$

where $[eqn]$ is defined as a separable function:

$$[eqn]$$

$[eqn]$ is a radial envelope—typically an isotropic Gaussian kernel—ensuring spatial localization, and the exponential term encodes angular selectivity. For $[eqn]$ , the filter is purely radial and yields a real-valued low-pass response, effectively acting as a smoothing operator, and for $[eqn]$ , the angular term introduces directional sensitivity, allowing the filter to highlight edge-like structures oriented along specific directions.

Recently, the harmonic filtering approach has been extended to non-Euclidean domain in [27], where multiscale anisotropic filters have been introduced. The Multiscale Anisotropic Harmonic Filters (MAHF) of m-th order and centered at $[eqn]$ act on a point cloud signal $[eqn]$ as follows:

$$[eqn]$$

and it is defined as a function of the geodesic distance metric between $[eqn]$ and $[eqn]$ , and of the angular coordinate of the j-th vertex on the graph $[eqn]$ tangent plane centered at $[eqn]$ in the following formulas:

$$[eqn]$$

where $[eqn]$ is the angular coordinate of $[eqn]$ on the tangent plane in $[eqn]$ . The function $[eqn]$ is a real weighting function depending on the connectivity of the point cloud graph G. Specifically, $[eqn]$ is the so-called heat diffusion kernel, formulated based on the theory of heat diffusion over smooth surfaces, and it is defined as follows. Let $[eqn]$ denote the eigenvectors of the graph Laplacian $[eqn]$ . The kernel is computed as

$$[eqn]$$

Hence, $[eqn]$ is related to the differences between the i-th and j-th coefficients across the eigenvectors of the Laplacian $[eqn]$ , and it is larger when the points $[eqn]$ belong to strongly connected areas [25]. MAHF has been applied to signals on point clouds for different tasks, such as point cloud denoising [25] or visual quality evaluation [28]. Still, MAHF filters rely on accurate estimation of the Laplacian associated with the point cloud, in turn depending on the quality of the triangulation algorithm. Point clouds from V-SLAM techniques often capture different objects at a wide set of different distances. This hinders the development of a regular mesh model by conventional triangulation techniques. Hence, useful feature extractors such as MAHF, which proved useful in application like denoising or point cloud quality evaluation, can fail to extract features due to sensitivity to Laplacian/mesh errors.

Making use of the side information provided by signals on edges and faces can lead to a richer estimate of the point cloud features. Specifically, we leverage the side information provided in TSP-SLAM to build an enhanced version of MAHF, suited to topological signal processing, which we refer to as T-MAHF.

In T-MAHF, we leverage the TSP-SLAM framework to generalize MAHF by exploiting side information available at the graph $[eqn]$ associated with the point cloud, as detailed below. We define the output of the T-MAHF topological filter at the i-th vertex as follows:

$$[eqn]$$

where, for a generic topological element $[eqn]$ representing a node, an edge, or a face, the function $[eqn]$ is written as follows:

$$[eqn]$$

where $[eqn]$ a real weighting function depending on a non-Euclidean distance metric between $[eqn]$ and a neighboring element $[eqn]$ , i.e., on the connectivity of the topological space $[eqn]$ associated with the point cloud, and $[eqn]$ a measure of angular distance in the tangent plane at the i-th vertex.

The TSP-SLAM framework allows us to compute the functions $[eqn]$ and $[eqn]$ appearing in $[eqn]$ with approximate values computed in the TSP-SLAM framework. To this end, let us develop the neighborhood system of the i-th vertex in the topological space, as established by the incidence matrices $[eqn]$ and $[eqn]$ . Let us denote the neighborhood of the i-th vertex by introducing the set $[eqn]$ of the neighboring faces, i.e., faces that are incident on the vertex:

$$[eqn]$$

Then, we introduce the set $[eqn]$ of neighboring edges as those belonging to the incident faces in $[eqn]$

$$[eqn]$$

Finally we denote the set $[eqn]$ of neighboring nodes:

$$[eqn]$$

In the following experiments, we set the weighting function as follows:

$$[eqn]$$

i.e., $[eqn]$ is an indicator function that specifies whether $[eqn]$ —which can be a node, face or edge—belongs to a neighborhood system of the i-th vertex or not. It should be noted that with this definition, the values of the weighting function in Equation (17) are larger (equal to one) for one-hop connected topological elements, and zero otherwise. This choice provides a hard approximation of topology element interactions extending the soft definition in (11). More generally, one could define the weighting function as a decreasing function of the topological distance $[eqn]$ , namely

$$[eqn]$$

where the indicator function adopted in this work corresponds to the simplest binary instance of this general formulation.

The T-MAHF expression is then rewritten as follows:

$$[eqn]$$

Let us observe that $[eqn]$ is a complex number, whose magnitude describes the intensity of the local variation and whose phase is related to the direction of the variation. In addition, we approximate the angular distance metrics $[eqn]$ , $[eqn]$ , and $[eqn]$ with their 2D counterparts. Specifically, we introduce the following approximations:

$$[eqn]$$ $$[eqn]$$

and

$$[eqn]$$

where $[eqn]$ and $[eqn]$ are the barycenters of the edge and faces in the neighborhood systems of the vertex i. The above quantities are directly available in the 2D domain, and in the following we show the limits within which the above approximation stands.

As for the computational architecture, the proposed MAHF and T-MAHF methods share an initial V-SLAM feature detection and sparse reconstruction, producing 3D points of the scene. MAHF operates on the reconstructed mesh, computing the Laplacian and its eigenvectors to build the heat kernel, which is combined with 3D angular information from mesh normals to produce a complex-valued filtered signal and enable 3D angle computation. This fully exploits the 3D structure but has higher computational cost. T-MAHF, instead, although relying on the mesh topology, works on the 2D image plane, performing 2D angular filtering without spectral decomposition. While computationally lighter, it captures less geometric information, as 3D surface variations are not explicitly modeled.

To sum up, TSP-SLAM allows building a texture dictionary suited to topological point cloud signal processing. The topology-related dictionary is built by extending techniques used in simultaneous localization and mapping algorithms, and adapting it to the point cloud topology description. This is useful for several developments, where non-Euclidean operators can be applied for processing purposes, including topological neural network architecture as discussed in [29].

4. Numerical Results

Herein, we provide a set of results describing how the TSP-SLAM framework is built on real data. We then show how it can be used for topological processing by presenting the application of T-MAHF, with reference to the dataset in [30]. Figure 3 illustrates the topology acquisition from measurements. Starting from 2D observations of projected point triplets in a stereo pair of images (bottom row), one can recover information about the underlying 3D triangle structure (top center). The middle row shows the schematic projection planes, where the corresponding left and right image coordinates, respectively $[eqn]$ and $[eqn]$ , are detected.

The two images in [30] used for topological processing were acquired with a stereo camera configuration, featuring a $[eqn]$ m baseline and a resolution of $[eqn]$ pixels. For keypoint detection and correspondence estimation, we adopted the method in [7], with the following settings: (i) the focal length of the cameras was 387.77 pixels for both x- and y- axes; (ii) the principal point of the cameras was located at coordinates (257.446, 197.718), specified in pixels; (iii) the maximum horizontal displacement between keypoints was limited to 48 pixels in order to filter matches that could lead to excessively shallow depth once triangulation is performed. By leveraging geometric constraints and correspondences between the two views, we inferred the spatial position of the points on the original 3D surface as in [7].

Figure 4 represents the detected keypoints from the stereo image pair, with a pseudo-color indicating the distance of the corresponding 3D points with respect to the camera.

The cloud of estimated 3D points is shown in Figure 5 (left), where the point color corresponds to its 3D depth $[eqn]$ . The point cloud was then equipped with a mesh using the Crust algorithm in [31], leading to a graph $[eqn]$ associated to the 3D point cloud, as shown in Figure 5 (right). The graph was defined with binary edge weights. The graph incidence matrices $[eqn]$ and $[eqn]$ , with values in $[eqn]$ , represent the graph connectivity in terms of edges and faces, thereby enabling the definition of a structure suitable for topological signal processing.

In the TSP-SLAM framework, the topological elements (points, edges and faces) are naturally associated with the signals available in the 2D domain, as illustrated in Figure 6. Specifically, Figure 6A shows the projection of the 3D edges $[eqn]$ in the 2D domain. At the barycenter $[eqn]$ of each edge, we plot a square indicating the average value of the luminance over the edge pixels. In the TSP-SLAM framework, this value will be used as the signal $[eqn]$ associated with the edge. Similarly, Figure 6B shows the projection of the 3D faces $[eqn]$ in the 2D domain. At the barycenter $[eqn]$ a circle indicates the average value of the luminance over the face pixels, i.e., the signal $[eqn]$ that will be associated with the face.

With these positions, we present some results exemplifying possible applications of TSP-SLAM in T-MAHF filtering.

Firstly, it is worth observing that the heat diffusion function $[eqn]$ , depending on a non-Euclidean distance metric between two 3D points $[eqn]$ and $[eqn]$ , is correlated to the distance $[eqn]$ between the two points in terms of graph hops.

This is illustrated in Figure 7 (blue circles), showing the scatter plot of the values of the heat diffusion kernel $[eqn]$ —computed using a Chebyshev polynomial approximation of heat diffusion over the graph—versus the number of hops $[eqn]$ between the same nodes. Figure 7 also shows a first-order exponential approximation of the underlying relationship. We observe a trend of exponential decay of the heat diffusion kernel $[eqn]$ versus the hop distances on the graph. This suggests that a reasonable weighting function can be obtained by retaining the closest topological elements.

This is the rationale behind the choice of approximating the T-MAHF weighting function with a neighborhood indicator function, as in (17).

Secondly, Figure 8 shows the scatter plot of the angles $[eqn]$ measured in the 3D domain versus the corresponding angles $[eqn]$ computed in the 2D image plane. Due to the ray tracing projection geometry underlying TSP-SLAM, the approximation is expected to hold when the angles belong to faces parallel to the image pane. As the faces are increasingly tilted in the 3D domain, the angles and their projection increasingly differ. This is evident in Figure 8, where the color of each point encodes the standard deviation of the depth ( $[eqn]$ ) coordinates of the triangle’s vertices in 3D space.

As expected, a low standard deviation indicates that the vertices lie (approximately) at the same depth, a condition in which the 3D triangle is approximately affine to the 2D projected triangle via ray tracing. For triangles with vertices at largely different depths, i.e., tilted with respect to the image plane, the 2D and 3D angles differ, but in most cases the scattered points lie along the bisectrix of the first quadrant, indicating that the 2D angles serve as an approximation of the corresponding 3D angles.

As illustrated in Figure 8, the discrepancy between 3D and 2D angular measurements increases with the depth variance of the face vertices, a phenomenon consistent with foreshortening effects reported in the structure-from-motion literature [32]. Reprojection errors and keypoint noise, as discussed in [33,34], are non-negligible in practical V-SLAM pipelines, but can be explicitly modeled to improve reconstruction reliability. Incorporating higher-order geometric primitives, such as faces and simplicial complexes, further helps to stabilize and refine point cloud representations by enforcing local planarity and structural consistency [35].

While exact 3D angle computation is computationally expensive, empirical evidence shows that 2D angular estimates align closely with their 3D counterparts when faces are approximately parallel to the image plane. Even for tilted faces, the scatter distribution remains concentrated around the bisector, indicating that 2D estimates provide a statistically reliable approximation. This trade-off allows TSP-SLAM to significantly reduce computational cost while maintaining geometric consistency at the topological level. Since the computation of the 3D angles $[eqn]$ is computationally expensive, the fast approximation provided by the 2D estimated angles $[eqn]$ can be adopted.

Finally, Figure 9 and Figure 10 present two examples of the result of topological harmonic filtering within the TSP-SLAM framework. For better visualization, the filter output at the 3D point $[eqn]$ is displayed at the corresponding keypoint $[eqn]$ , overlaid on the image. The color of $[eqn]$ in the top panel represents the absolute value of the T-MAHF filter output on a logarithmic scale, i.e., $[eqn]$ , while the bottom panel represents, in the same manner, the phase of the T-MAHF filter output on a logarithmic scale, i.e., $[eqn]$ . We observe that T-MAHF allows the ranking of points according to larger magnitudes $[eqn]$ ; for these points, the phase $[eqn]$ is an estimate of the direction of signal variation. In the bottom panel, detail in the right box, a few points characterized by larger magnitudes $[eqn]$ are highlighted to illustrate the meaning of the phase component $[eqn]$ . Although visually distinct, the selected points reveal a direction of the signal variation along angles of 0 (light green), $[eqn]$ (red) and $[eqn]$ (fuchsia). Therefore, the variation corresponds to a horizontal discontinuity. The same applies to the detail in the left box, where horizontal or slightly tilted edges are recognized.

It is worth noting that, although visualized in the 2D domain for clarity, this information refers to the signal defined on the point cloud graph and its topology. Such information can be employed for various applications, including point cloud classification or loop closure recognition. At the same time, thanks to the TSP-SLAM framework, this information can be computed directly on 2D data, increasing robustness and reducing computational complexity compared to computations in the full 3D domain.

A few remarks are in order regarding the applicability of the method to point clouds with larger numbers of points. In principle, the approach scales naturally to different resolutions, as illustrated in Figure 11 and Figure 12. In these examples, we considered the stereo projections—obtained with the same camera parameters as in [30]—of a synthetic toroidal grayscale point cloud with $[eqn]$ and N = 20,000, respectively. For both cases, in addition to the 2D projections, we show the point cloud signal $[eqn]$ and the output of the T-MAHF filtering $[eqn]$ in terms of its magnitude $[eqn]$ and phase $[eqn]$ , computed from the 2D signal information. In a practical V-SLAM framework, when points are acquired from only two cameras, the 3D point set is restricted to the salient points for which reliable correspondences can be established. Conversely, if the cameras are complemented with Lidar data, the generated point cloud becomes significantly richer.

Finally, we carried out a comparison of the computational complexity of the T-MAHF method with the existing MAHF [25]. The two algorithms share several tasks—such as salient points extraction and matching, and 3D point cloud and mesh computation—but differ in the filtering stage. In T-MAHF a binary adjacency matrix is computed, and filtering is carried out through signal extraction and distance/relative angle computation. In contrast, the MAHF algorithm requires computing the real-valued conformal Laplacian and performing its eigen analysis to calculate the filtering weights as functions of $[eqn]$ and $[eqn]$ . Table 2 reports the computation time of the two methods in milliseconds per point, analyzed on four frames (1,137,981 and 1041 in [30]), and broken down by task. On average, T-MAHF achieves a total computation time reduction of about 40% compared to MAHF. Let us observe that the complexity of the face and edge signal extraction task can increase as far as an increasing number of pixels is considered in the averages in (5) and (6). Still, T-MAHF allows faster point cloud topological signal processing exploiting the potentiality of the T-SLAM framework.

5. Conclusions

Topological point cloud processing is gaining momentum, particularly for objects acquired though dedicated multi-camera systems. This paper presents a methodology to build a texture dictionary suited for topological point cloud signal processing. The topology-related dictionary is constructed by extending techniques used in simultaneous localization and mapping (SLAM) algorithms, and adapting them to the point cloud topology description. This approach is useful for various developments where non-Euclidean operators can be applied for processing purposes.The proposed framework, validated through numerical experiments, demonstrated stable heat diffusion consistent with topological distances, reliable approximations of 3D geometric relations from 2D projections, and improved robustness to mesh irregularities. Quantitative analyses further confirmed the effectiveness of TSP-SLAM in enhancing denoising, segmentation, and perceptual quality evaluation tasks, showing clear advantages over graph-based baselines. Future work will focus on extending the framework to more general topological structure, including cell complexes.