## LEARNED COMPRESSION OF POINT CLOUD GEOMETRY AND ATTRIBUTES IN A SINGLE MODEL THROUGH MULTIMODAL RATE-CONTROL

Michael Rudolph * 1,2 , Aron Riemenschneider †1 , and Amr Rizk ‡2 1 University of Duisburg-Essen, Germany 2 Leibniz University Hannover, Germany

## ABSTRACT

Point cloud compression is essential to experience volumetric multimedia as it drastically reduces the required streaming data rates. Point attributes, specifically colors, extend the challenge of lossy compression beyond geometric representation to achieving joint reconstruction of texture and geometry . State-of-the-art methods separate geometry and attributes to compress them individually. This comes at a computational cost, requiring an encoder and a decoder for each modality. Additionally, as attribute compression methods require the same geometry for encoding and decoding, the encoder emulates the decoder-side geometry reconstruction as an input step to project and compress the attributes.

In this work, we propose to learn joint compression of geometry and attributes using a single, adaptive autoencoder model, embedding both modalities into a unified latent space which is then entropy encoded. Key to the technique is to replace the search for trade-offs between rate, attribute quality and geometry quality, through conditioning the model on the desired qualities of both modalities, bypassing the need for training model ensembles. To differentiate important point cloud regions during encoding or to allow view-dependent compression for user-centered streaming, conditioning is pointwise, which allows for local quality and rate variation. Our evaluation shows comparable performance to state-of-the-art compression methods for geometry and attributes, while reducing complexity compared to related compression methods.

Keywords Virtual Reality, 6DoF, Point Cloud, Variable Rate, Geometry Compression, Attribute Compression

## 1 Introduction

Point clouds are a popular representation for immersive multimedia content [1, 2], benefiting from low processing overhead compared to other formats for three-dimensional data such as meshes. The idea of using point clouds as a data structure for streaming dynamical volumetric multimedia experiences over the Internet leads to significant bandwidth requirements, which establishes the need for efficient compression approaches . Standardization proposes solutions for point cloud compression through G-PCC [3], which uses octree data structures for the geometry and tailored approaches for the attributes like Region Adaptive Hierarchical Transform (RAHT) [4]. When considering point cloud sequences, V-PCC [5] utilizes video codecs to compress projections of the point cloud, thus allowing to exploit inter-frame correlations in 2D.

Inspired by and going beyond classical compression methods, Deep Learning has shown promising results for the compression of images [6-8], videos [9] and point clouds [10-16]. Similar to the traditional approach, learned point cloud compression mostly treats geometry and attribute compression separately . A geometry codec handles the

∗

michael.rudolph@ikt.uni-hannover.de

†

aron.riemenschneider@stud.uni-due.de

‡

amr.rizk@ikt.uni-hannover.de

Figure 1: Variable-Rate Coding using a joint architecture for geometry and attributes, allows to freely select attribute and geometric quality with only a single model . The figure shows one mapping (path) of attribute and geometry quality configurations q ( G/A ) that compose a quality (PCQM) Pareto front with the corresponding rates given in bit per point (bpp).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-c1e4d0a2-9689-496a-9258-66fd0280fffe.jpg)

compression of the points, while an attribute compression model handles the compression of colors given the knowledge of the undistorted geometry .

When combining the individual geometry and attribute codecs, whether learned or traditional, while relying on lossy compression, the reconstructed geometry at the decoder diverges from the geometry necessary as a prior for attribute decoding. This impacts the quality of the attribute reconstruction. To address this issue, the geometry is decoded and the attributes are projected on the reconstruction [3, 5, 17] before attribute encoding, ensuring they are compressed using the same geometric prior which would be available at the receiver side. We observe, however, that this introduces additional distortions and adds to the computational complexity of the encoding, requiring not only two individual encoding steps but additional geometry decoding and attribute projection, mostly done through costly nearest-neighbor searches, at the sender.

As depicted in Fig. 1, we propose a joint variable-rate compression method for point clouds in a single, adaptive autoencoder that embeds both modalities into a unified latent space. Figure 1 shows that using this joint model, we can individually gauge the quality of the geometry and the attributes. Using our single model dispenses with training an ensemble of models to fit a fixed number of individual tradeoffs between the bitrate and the qualities of the geometry and attributes. The contributions of this paper are:

- · We formulate the rate-distortion loss for multimodal compression of geometry and attributes to incorporate adaptive weighting of conditional loss components. This allows to train a single adaptive model covering a wide range of geometry and attribute qualities and thus rates.
- · We propose a joint compression model, handling the geometry and attributes in one architecture, avoiding the need for geometry decoding and reprojection on the sender side.
- · Utilizing conditional encoding, attribute and geometry quality can be modulated independently, allowing variable rate-encoding with a single model, as well as, region-of-interest encoding of point clouds.

## 2 Related Work

Attribute compression is essential for point cloud multimedia applications in contrast to solely compressing geometry, e.g., in LiDAR applications [18]. Merkuria et al. [19] developed a codec utilizing octree partitions and intra-frame prediction while projecting attributes to 2D grids and compressing them with legacy image codecs. MPEG standards [20]

distinguish between static point cloud compression using Geometry-based PCC (G-PCC) [3] and dynamic point cloud compression with Video-based PCC (V-PCC) [5]. While G-PCC relies on octree partitions for geometry compression and methods like Region Adaptive Hierarchical Transform (RAHT) [4] for attributes, V-PCC projects geometry and attribute information into video frames and leverages video codecs for compression. Driven by promising results in learned image compression [6-8], learning based solutions for point cloud geometry or attribute compression showed competitive results in both domains. However, most of these studies [10, 12-14, 16, 21] do not consider the case where both, geometry and attributes, are subjected to distortions, as it would be the case for a realistic transmission use case.

We note that early work explored learned compression of attributes and geometry in a single model [15], however, learned compression methods using separate models for both modalities [17] were shown to outperform the single model approach. The latter further boosts performance through wrapping RAHT [4] implemented in G-PCC [3] in the attribute branch and applying adaptive filters [22], which need to be optimized during each encoding, adding to the already large computational cost.

## 2.1 Learned Geometry Compression

Learning-based methods demonstrated competitive performance for point cloud compression, where early work [13, 16] proposed to use 3D convolutions on voxelized point clouds and leverage learned entropy models [6] to optimize a rate-distortion loss. The reconstruction is formulated as a binary classification problem using the (weighted) Binary Cross Entropy or the Focal Loss [23] to handle class imbalance due to point clouds sparsity in 3D space. The authors of [14] employed entropy models from [8], allowing to model spatial correlations combined with auto-regressive prediction based on already decoded coefficients. Utilizing the above reviewed contributions, the authors of [10] employ a Variational Autoencoder (VAE) architecture based on Voxception-ResNet [24], allowing to show on-par performance with V-PCC [5] in terms of geometry compression. We note, however, as all previously mentioned methods utilize 3D convolutions on sparse data represented in dense data structures, processing the point cloud in cubes is necessary to handle the memory requirements at the cost of limiting throughput and performance. To overcome this, sparse convolutions [25] have been employed in [11] to process large-scale point clouds in one pass and to exploit global context during hyper-encoding, while a multi-scale classification in each stage of the decoder prunes non-occupied voxels. Instead of using an encoder-decoder approach, [12] propose to encode the upsampling for each stage individually while relying on group-based decoding to leverage correlations between upsampled voxels. The model in [12] is then extended for inter-coding of dynamic point cloud sequences [21].

## 2.2 Learned Attribute Compression

Attribute compression is still dominated by traditional approaches, relying on Graph Transform [26] or RAHT [4], with the latter being part of G-PCC. In the field of learning-based approaches, models using sparse convolutions achieved comparable performance [12]. Combining traditional methods with learning-based approaches, the authors of [27] estimate the probabilities of transform coefficients from RAHT with deep entropy models, further improving compression efficiency. The attribute compression branch in YOGA [17] relies on a thumbnail point cloud compressed at low resolution through RAHT to then predict residuals of the attributes through a conditioned autoencoder similar to [12]. To further increase the quality of the reconstructed attributes, YOGA [17] employs adaptive filters [22], which are optimized during encoding and the parameters are added as side information to the bitstream to be used at the receiver.

## 2.3 Variable-Rate Control for Learned Compression

Numerous compression models for images [7, 8] and point clouds [10-16] based on the non-linear transform coding framework [6] achieve remarkable rate-distortion performance. A major drawback of these methods is that they are conditioned on a fixed tradeoff between compression and reconstruction at training time. As a result, an ensemble of models needs to be trained when multiple qualities are required, as it is commonly done in adaptive streaming [28]. Additionally, this further reduces flexibility as encoding configurations need to be fixed through hyper-parameters at training time.

Variable-rate control and progressive coding schemes are promising techniques to freely adapt the rate-distortion tradeoff using a single model for encoding and decoding. For example, for learned image compression, progressive coding is realized by iteratively refining the quantization grids [29]. Trit-planes [30] or dead-zone quantizers [31] extend this idea to better fit the distribution of quantization errors. Following a different approach, quantization or reconstruction results are compressed in cascading [32] or recurrent compression [33-35] models, allowing for progressive coding. In contrast to these architectural decisions, training with uniformly sampled tradeoff parameters λ while dropping a

corresponding number of feature channels in the entropy bottlenecks allows progressive coding through sorting feature channels by importance [36] and selecting a number of coefficients to decode.

Compared to progressive methods, adaptive rate-control targets conditioning a single model on a rate variable to allow controlling the rate-distortion of the compressed representation. Early work such as in [37] introduces a gain control layer, multiplying and shifting channels of the feature representation to control the quantization before encoding the transformed coefficients in the entropy bottleneck. By incorporating multiple gain control layers in the encoder and decoder, features are scaled hierarchically [38-40]. Additionally, invertible layers allow increasing the reconstruction quality [40]. YOGA [17] adopts these ideas for variable-rate compression with the layers proposed in [39] using sinusoidal embeddings, but requires an ensemble of entropy models as well as selecting appropriate rate parameters for the wrapped color compression using G-PCC.

Finally, we note that works such as [38] emphasizes regions of interest in images by conditioning compression models on 2D quality maps instead of a single scalar quality parameter. As the quality map is not available on the decoder side, it is inferred from the compressed representation for decoding. Extending this idea from image compression to point cloud compression is promising as it allows view-dependent encoding of three-dimensional content as done with traditional codecs [41, 42].

## 3 Challenge and Approach

## 3.1 Problem Statement

We consider the problem of compressing a point set P = ( g a , ) containing a set of point locations g , denoted geometry, and corresponding attributes a . Specifically, for the geometry, quantized coordinates g = { ( x , y , z i i i ) : x , y , z i i i ∈ N 0 , x , y i i , z i &lt; R } with maximum value R ∈ N are assumed 4 . Note that for real-valued coordinates, a voxelized representation could be easily derived by rounding and partitioning the full 3-dimensional space into blocks with size R 3 . We refer to attributes (and features) of a sparse tensor at the geometrical location u ∈ g as a u .

We observe that state-of-the-art compression methods [3, 5, 17] handle both modalities, i.e., geometry and attributes, separately. Specifically, methods such as [10, 12-14, 16, 21] or [12] solely compress geometry or attributes, respectively. Methods that compress both modalities such as [3, 5, 17] still resort to separating the individual compression procedures for geometry and attributes. In essence, they encode the geometry first, which is then decoded to project the attributes to the possibly distorted geometry after reconstruction. Finally, they compress the projected attributes on the distorted geometry. On the decoder side, the geometry is decoded and then used as side-information to allow attribute decoding. While this process allows compression of point cloud attributes under geometric distortion, it has a number of disadvantages: i) Encoding and decoding requires a sequence of non-parallelizable steps, resulting in high latency and ii) reprojecting the attributes on distorted geometries at encoding introduces additional distortions before compression.

Now, we aim at training a model for joint attribute and geometry compression. In a Lagrangian fashion, we seek to minimize the rate R under the quality constraints on geometry D ′ G and attributes D ′ A given the loss function

$$\mathcal { L } = \mathcal { R } + \lambda _ { A } \mathcal { D } _ { A } ^ { \prime } + \lambda _ { G } \mathcal { D } _ { G } ^ { \prime }$$

with corresponding deterministic, positive scalars λ A , λ G . This can also be understood as finding a tradeoff between geometry and attribute quality at a given rate. Compared to compression models that handle only a single modality in form of L = R + λ D , with fixed hyperparameter λ , the loss function above adds another dimension to the number of models required to train, tune and deploy, making ensembles of joint models inadequate for most use cases. Hence, a key to solving this challenge is providing a single, adaptive model. This stands in contrast to the image compression literature [6-8] that coined the rate-distortion tradeoff in a bijection sense between the rate R and the quality D . In this work, this bijection does not exist as the rate results from defining a tradeoff between geometry and attributes.

A seemingly promising approach is variable rate compression [38], in the sense

$$\mathcal { L } = \mathcal { R } + \mathcal { D } _ { A } + \mathcal { D } _ { G }$$

where we condition a single model for compression at multiple rates, thus only requiring a single training run and reducing the memory requirements at deployment. In contrast to (1), D A , D G now integrate locality weighting functions λ A , λ G instead of scalar quality parameters λ A , λ G allowing to define a location-specific tradeoff between the geometry, attributes and the rate. This is especially suited for multimedia point clouds, as it allows to differentiate important regions such as faces, or to offer view-dependent encoding. Additionally, instead of searching for tradeoffs between both modalities and the rate in (1) at training time, e.g. through grid search, attribute and geometric quality can be balanced post-training.

4 This allows using sparse data structures and thus memory and computational efficient convolutional neural networks.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-04a6e8ac-4e29-453a-9217-73fe784ba2f8.jpg)

(a) Transform Coding [7, 8]

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-90fbf5e9-924e-45cc-ad58-9ab654910031.jpg)

(b) Conditional extension [38]

Figure 2: Operational diagrams of the Mean-Scale Hyperprior model [7, 8] and the conditional addition presented by [38].

## 3.2 Intermezzo: From Non-Linear Transform Coding to Conditional Transform Coding

Next, we briefly review the non-linear transform coding framework from [6] in addition to its extension towards conditional variable rate coding [38]. The following review, also depicted as building blocks in Fig. 2, considers an arbitrary data vector x . In a later section, we show our building blocks where some are based on this framework.

Figure 2a shows the operational diagram of the mean-scale hyperprior [7, 8] where an input data vector x is non-linearly transformed through an encoder g a . The transform output, i.e. the coefficient vector y , is quantized as ˆ y = Q ( y ) and fed to the decoder g s producing the reconstruction ˆ x = g s ( ˆ y ) . Note that the quantized coefficients are usually entropy coded for transmission while the encoder and decoder non-linear transforms g a , and g s are learned.

Entropy coding the quantized coefficients ˆ y requires an estimate of their distribution. Given the observation that the coefficients vary locally in mean and scale, a hyperprior model allows to parameterize a local estimate in form of a Gaussian distribution [7, 8]. A hyperprior z = h s ( y ) is inferred from the coefficients y as depicted in Figure 2a. The quantized hyperprior ˆ z is then entropy coded using a factorized entropy model [7] and used to locally reconstruct the parameters of the Gaussian distributions necessary for encoding/decoding ˆ y using the hyperdecoder h s . Since backpropagation through a quantization operation results in gradients becoming zero almost everywhere, quantization is replaced through additive uniform noise [6] for training, i.e. ˜ y = ˜ ( Q y ) = y +∆ with ∆ being an independently and identically distributed uniform random variable, ∆ ∼ U [ -0 5 0 5] . , . . In the following, we refer to all variables affected by this quantization proxy through ˜ ) ( · .

Extending this approach to allow variable rate coding, the encoder g a in Fig. 2b is now conditioned on a variable c which is inferred through the conditional encoder p a from a quality variable q [38]. This quality variable allows to locally dictate the quality of corresponding elements in the data vector x . As the quality variable is not available for decoding, a proxy quality variable ˆ q is inferred from the hyperprior ˆ z . This allows to derive the condition variable ˆ c required in the synthesis transform g s .

## 3.3 Multimodal Conditioning for Geometry and Attributes

↦

Next, we show how to achieve variable rate compression of geometry and attributes in a single adaptive model through integrating the multimodal rate-distortion loss in (2) into the conditioned transform coding (cf. Sec. 3.2). To this end, we provide a per-point quality map q = [ q ( G ) , q ( A ) ] to dictate the target quality individually for both modalities, i.e. geometry and attributes, to the model. The quality map instantiates the quality variable q in Sec. 3.2 residing on the object geometry g . During training, each element of the quality map, which is set to q ( A ) u , q ( G ) u ∈ [0 , 1] , with u ∈ g , is mapped to a corresponding weighting map element λ ( A ) u , λ ( G ) u through a monotonically increasing function [38]. This allows to individually define the range of the weights in λ ( A ) , λ ( G ) and thus the target distortions for the respective modality. We use a quadratic weighting function , i.e. 5 T A : [0 , 1] → [ λ min , λ max ] with T A ( x ) = ( λ A max -λ A min ) x 2 + λ A min , such that λ ( A ) u = T A ( q ( A ) u ) , with u ∈ g . The same mapping holds for λ ( G ) u . The range [ λ A/G min , λ A/G max ] defines the minimum and maximum values in the form of distinct hyperparameters. At training time, quality maps q are either uniformly drawn or supplied as a gradient in a random direction, allowing the model to vary the compression ratio both globally and locally.

5 We note that exponential functions can be used as in [38], however, experiments show that a root function is more suitable.

Figure 3: Model Architecture, materializing the operating diagram in Fig. 2b. For Convolutions, b 3 indicates filters with kernel size b and dimension 3 . The arrows ↑ / ↓ symbolize upsampling and downsampling, respectively. A more detailed version of the architecture containing implementation details is given in the supplementary materials.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-53b36dca-aa03-4737-b3a1-8766cb612329.jpg)

## 3.3.1 Geometry Loss

For the geometry reconstruction, we use the Focal Loss [23] at every level of the decoder as proposed by [11] to allow iterative upsampling of the features and geometry in K steps 6 . In a nutshell, iterative upsampling in K steps is required for constructing the point cloud geometry from the downsampled representation of the encoder. More specifically, in each upsampling step k ∈ { 0 , . . . , K - } 1 , a set of possible point locations ˜ d ( k ) , is generated from the lower scale geometry ˜ g ( k +1) s.t. ˜ g ( k +1) ⊂ ˜ d ( k ) [43] . Finally, an occupancy probability 7 ˜ p ( k ) u at location u and scale k is inferred for the possible point locations ˜ d ( k ) . Incorporating the geometry loss map λ ( G,k ) u in each scale k , the distortion loss D G in (2) is

$$\mathcal { D } _ { G } = - \sum _ { k = 0 } ^ { K - 1 } \sum _ { \substack { u \in \tilde { d } ^ { ( k ) } \\ \lambda \dots } } \lambda _ { u } ^ { ( G, k ) } ( 1 - \dot { p } _ { u } ^ { ( k ) } ) ) ^ { \gamma } \log ( \dot { p } _ { u } ^ { ( k ) } )$$

$$\text{with} \, \dot { P } _ { u } ^ { ( k ) } = \begin{cases} \tilde { P } _ { u } ^ { ( k ) } & \text{if} \, u \in g ^ { ( k ) } \\ 1 - \tilde { P } _ { u } ^ { ( k ) } & \text{else} \end{cases},$$

where g ( k ) is the ground truth geometry at scale k and the hyperparameter γ known from the focal loss [23] reduces the gradients for well-classified samples, allowing the model to emphasize harder samples. As the loss map λ ( G,k ) u is only given on ground truth voxel geometry, i.e. u ∈ g , substitute values λ ( G,k ) u are derived through downsampling and upsampling using average pooling and transposed average pooling, respectively, resulting on the same base geometry as in ˜ d ( k ) .

## 3.3.2 Attribute Loss

The attribute reconstructions is typically cast as an element-wise regression problem [12, 15]. Using the L 2 -norm the attribute component D A in (2) becomes

$$\mathcal { D } _ { A } = \sum _ { u \in g \cap \tilde { g } } \lambda _ { u } ^ { ( A ) } \| \mathfrak { a } _ { u } - \tilde { \mathfrak { a } } _ { u } \| _ { 2 } ^ { 2 } \,,$$

where a u denotes the attributes at the geometrical location u ∈ g and ˜ a u , ˜ g are, as discussed in Sect. 3.2, affected by the quantization proxy. Recall that the attribute loss map λ ( A ) u contains individual, i.e. local, weighting to gauge the attributes to target distortion.

Note that the loss can only be computed on locations u ∈ g ∩ ˜ g , i.e. the intersection of ground truth geometry g and the reconstruction ˜ g , where we compare the ground truth attributes a u with the reconstruction ˜ a u .

6 corresponding to K resolutions also denoted as scales throughout this work.

7 Later, we omit the scale superscript for k = 0 , i.e. the original resolution ˜ g (0) := ˜ g g , (0) := g .

Figure 4: Conditioned feature extraction (CFE) block. A condition map c ( k ) = [ α ( k ) , β ( k ) ] is derived from the quality map q ( k ) at scale k for element-wise scaling and shifting of features according to (7). Consequently, the module learns a channel-wise weighting of features in t ( k ) according to the specified qualities q ( A ) and q ( G ) .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-ed7d824b-8c30-4318-8ed6-a8e74b0541da.jpg)

## 3.3.3 Rate Loss

The rate loss R is computed according to [7], aiming to minimize the cross-entropy as in

$$\mathcal { R } & = - \sum _ { u \in \mathfrak { g } ^ { ( K ) } } \sum _ { i } \log _ { 2 } \left ( p _ { \tilde { y } | \tilde { \mu }, \tilde { \sigma } } ( \tilde { \mathbf Y } _ { u, i } | \tilde { \mu } _ { u, i }, \tilde { \sigma } _ { u, i } ) \right ) \\ & \quad - \sum _ { u \in \mathfrak { g } ^ { ( H ) } } \sum _ { i } \log _ { 2 } \left ( p _ { \tilde { z } | \phi } ( \tilde { z } _ { u, i } | \phi ^ { ( i ) } ) \right ),$$

where Φ = { ϕ ( ) i } i are the learnable parameters of the factorized entropy model [7]. We use ˜ y u,i and ˜ z u,i to denote the transformed coefficient of ˜ y and ˜ z at the i th channel in the feature dimension at the geometry position u , respectively. In addition, ˜ µ u,i , σ ˜ u,i denote the local estimates of the means and scales that are obtained from the hyperprior ˜ z as in ( ˜ ˜ µ, σ ) = h s ( ˜ z ) . Hereby, the hyperprior ˜ z capturing the estimates resides on lower scale geometry g ( H ) , i.e. H &gt; K .

## 3.4 Model Architecture

We propose a sparse neural network architecture as depicted in Fig. 3, following the framework proposed by [38] (cf. 2b). Recall the notation conventions from Sec. 3.2, with ˆ ) ( · denoting variables affected by quantization (i.e. inference time), and their training time counterparts (quantization proxy) given as ˜ ) ( · .

The model receives a point cloud P = ( g a , ) . The encoder g a downsamples the point cloud P to the (resolution) scale k ∈ { 0 , . . . , K - } 1 , defined in Sect. 3.3.1, in combination with the quality map q residing on the same geometry g , to obtain a transformed representation ( ˆ g ( K ) , ˆ y ) .

The key to rate adaptive feature transformation at each scale k is the Conditional Feature Extraction (CFE) block (see Fig. 4) in the encoder that takes the quality map q as input. On the decoder side, the surrogate 8 quality map ˆ q ( K ) required for adaptive decoding at the lowest scale 9 is inferred from the hyperprior ˆ z and the transformed coefficients ˆ y . The second ingredient to our method is the Upsampling and Prunining (UP) block during decoding (see Fig. 5) which prunes the geometry at each scale k to maintain the sparsity of the input geometry, thus overcoming memory restrictions [11].

## 3.4.1 Conditional Feature Extraction on Geometry and Attribute Importance

We employ local conditioning of feature representations through a scale and shift operation, following common procedure to allow variable-rate control in image [37, 38] and point cloud compression [17]. Hereby, features t ( k ) residing on the geometry g ( k ) at a given scale k are transformed through the local scale and shift operation into t ′ ( k ) as

$$\mathfrak { t } _ { u } ^ { \prime ( k ) } = \alpha _ { u } ^ { ( k ) } \mathfrak { t } _ { u } ^ { ( k ) } + \beta _ { u } ^ { ( k ) } \quad \forall u \in \mathfrak { g } ^ { ( k ) }.$$

8 we call estimates at the decoder side surrogate .

9 in this work we set the number of resolution scales to K = 3 .

Figure 5: The Upsampling and Pruning (UP) block generates a locally dense geometry ˆ d ( k ) from the lower scale points ˆ g ( k +1) . Per-point occupancy probabilities p ( k ) are predicted from the upsampled and transformed features f ( k ) . Finally, the top-k most probable points are selected to prune the dense point cloud ( d ( k ) , f ( k ) ) , arriving at ( ˆ g ( k ) , ˆ f ( k ) ) .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-10711673-66d3-40f5-a987-a5de12bd2d58.jpg)

using the scales α ( k ) u and shifts β ( k ) u . In contrast to the related work [17, 38], which only predicts element-wise scales α ( k ) u and shifts β ( k ) u from a single quality map value q ( k ) u per geometry location u , joint adaptation of the geometry and attribute quality requires learning a mapping from geometry and attribute quality to the local transformation variables.

In essence, the feature extraction layer learns the importance of features for optimizing a given task by varying its scale and shift. This results in broader distributions at the bottleneck when aiming for higher rates, thus reducing quantization errors. Vice versa, dictating low qualities through q ( k ) results in more narrow distributions or in disabling certain feature channels completely through trivial distributions. Additionally, the additive shifting operation allows to encode the quality map q into the transformed coefficients y , embedding information required to infer the surrogate quality map from ˆ y .

## 3.4.2 Upsampling Sparse Features with Geometry Pruning

During encoding, lower scale geometry g ( k +1) is derived from g ( k ) through rounding in a dyadic scheme. As a result, when upsampling a representation ˆ g ( k +1) in the decoder, the correct voxels to be occupied in ˆ g ( k ) are unknown. Hence, to prevent information loss, voxels are generated at each possible location through a Generative Transposed Convolution [43], resulting in a dense geometry ˆ d ( k ) , which is a superset of g ( k ) . While the continuation of decoding with ˆ d ( k ) is conceptually possible, it leads to significant increase in memory and computational requirements. Geometry compression methods overcome this challenge by performing multiscale classification [11, 17], pruning the upsampled geometry ˆ d ( k ) at each scale k by inferring per-point occupancy probabilities ˆ p ( k ) u for all u ∈ ˆ d ( k ) . Note that attribute decoders, such as [12, 17], do not face this challenge by assuming the availability of undistorted geometry g ( k ) from a preceding step, which is in practice not true at the decoder.

Now, when combining both modalities, i.e. geometry and attributes, we need to predict the occupancy and prune the feature representation for the next scale. The challenge hereby is maintaining an expressive feature representation after pruning under potential geometry miss-predictions for attribute reconstruction. To increase robustness against these miss-predictions, we directly transform the dense, upsampled features, arriving at a dense feature representation, as shown in Fig. 5. Then, an occupancy prediction branch infers occupancy probabilities ˆ p k u from the transformed, dense feature representation f ( k ) . Finally, top-k selection is performed, keeping the k most probable voxel locations in a pruning step on the transformed, dense features.

## 4 Evaluation Results

## 4.1 Evaluation Setup

## 4.1.1 Implementation

We implement our approach in PyTorch, using MinkowskiEngine [25] to leverage sparse convolutions. For training, we collect a variety of human point clouds from the UVG dataset [44], select a number of frames for each point cloud sequence and crop them into cubes with edge length 128 , resulting in 18816 distinct samples. Opposed to projecting image textures on the samples during training [12], we solely subject the attributes to color jitter and randomly rotate each sample. For each sample, a random quality map ( g q , ) is derived according to Sec. 3.3. We train for 200 epochs

Figure 6: Pareto fronts for the PCQM-Rate tradeoff, optimized for the given point cloud.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-71232ce2-5adf-4dae-b186-74760546e435.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-c420fff8-05e3-4ab9-8fee-70c6eead177c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-f08eaf67-786c-4967-9978-87039cf6080d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-df1af472-9d74-484c-a9d1-0a6b0c469ff4.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-8a1da218-2dd6-49b7-9455-3ef029c99cff.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-8380ef16-5c8a-4ae4-8821-d2839897790f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-4492d38e-5d88-4598-9d40-a5099a3aad5c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-81fd1f5a-d58d-4a9f-918b-c3b5eb5d5d82.jpg)

using a batch size of 16 . Adam [45] is used for optimization with an initial learning rate of 10 -3 , which we decrease with factor 0 1 . every 50 epochs. To stabilize training, gradient norm clipping with threshold 1 0 . is integrated.

## 4.1.2 Baselines

We compare our method against YOGA [17], which we identify as the current state-of-the-art method for learned compression of geometry and attributes. Additionally, we use the tmc2v24 test model for V-PCC [5], using all-intra test conditions and the tmc13v23 implementation of G-PCC [3] employing octree compression for geometry and RAHT [4] for attributes. For both methods, we select the first 4 rate configurations specified in the respective test models, as we find that they result in a similar range of rates as our model.

## 4.1.3 Evaluation Metrics

We measure the bitrate in bits per point (bpp), using the cardinality | g | as divisor. Geometric quality is assessed using the D1-PSNR [46] metric, computing the mean distances between points of the tested geometry g and their nearest neighbor in the reference ˆ g . Utilizing the derived neighbor association, the Y-PSNR and YUV-PSNR metrics [47] are computed to assess the attribute quality, with the former only using the luminance channel and the later referring to the weighted average over all channels using weights ( 6 8 , 1 8 , 1 8 ) , following evaluation protocols in related literature [12, 17]. We consistently present symmetric metrics, i.e. computation is performed for swapped source and reference while reporting the stricter result. Finally, PCQM [48] is used as a joint metric, incorporating geometric and attribute features. Note that for all presented results, the reported bitrates correspond to the compression of geometry and attributes together.

We select 4 frames each from the 8iVFBv2 [49] and MVUB [50] dataset for testing, which are specified in the supplementary material. If not specified otherwise, uniform quality maps q ( A ) and q ( G ) are used in our model to allow fair comparison.

## 4.2 Perceptual Quality Optimized Pareto Front

Following the evaluation protocol specified in [17], we traverse the cross product of quality configurations for our approach with step size 0 05 . , resulting in 400 configurations for the point clouds in 8iVFBv2 [49]. Then, the Pareto front for the tradeoff between PCQM [48] and rate is determined and depicted in Fig. 6. Two examples of the chosen configuration pairs can be found in Fig. 7, which shows the pairs of attribute quality q ( A ) and geometric quality q ( G ) contributing to the Pareto front as a path over all possible combinations. The underlying contour in each plot shows the bitrate, attribute quality and geometric quality, resulting from all possible combinations. Similarly, Pareto fronts for

Figure 7: Mapping of attribute and geometry quality configurations q ( A ) , q ( G ) to the metrics bit per point (bpp), Y-PSNR and D1-PSNR for redandblack (top) and loot (bottom). The paths shows configuration pairs q ( A ) , q ( G ) that stem from the Pareto fronts in Fig. 6.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-8df7a569-ebb2-41e3-bcd3-02ef701b3e3f.jpg)

V-PCC [5] and G-PCC [3] are determined by traversing the full range of attribute QPs for each geometry setting. We also compare against the results reported in the original paper of YOGA [17] for the 8iVFBv2 [49] dataset.

From Fig. 6 we find that our model shows improved PCQM performance compared to all baselines for medium to high rates, indicating comparable or better perceived quality of the reconstructions at this compression ratio. However, when aiming for lower rates, the quality deteriorates faster compared to related methods. We attribute this to the fact that our approach is the only method solely relying on learned compression for the point clouds, while, e.g. YOGA [17] not only inherits low rate compression capabilities by compressing a thumbnail attribute representation using RAHT [4] but also employs adaptive filters [22] which are optimized during encoding on the specific content to counteract quality degradation.

We note that the shown results present optimized configurations for the respective point clouds. This is also apparent when comparing the selected configuration pairs for the point clouds redandblack and loot contributing to the Pareto fronts in Fig. 6, which are visualized as paths in Fig. 7. Depending on the individual attribute and geometric properties of the content, we notice that the paths and thus the chosen set of configurations to arrive at the Pareto fronts, differ significantly. Recall that this evaluation required traversing a large number of possible configurations to select the best tradeoff points to arrive at the pareto fronts and is thus not practical. Nevertheless, this evaluation provides insights into the capabilities of the tested compression methods.

## 4.3 Fixed-Configuration Encoding Performance

Now, we set the conditioning of our model to a fixed set of configurations, i.e. ( q ( G ) u , q ( A ) u ) ∈ { (0 05 . , 0 1) . , (0 1 . , 0 2) . , (0 2 . , 0 4) . , (0 4 . , 0 8) . } ∀ u ∈ g , independent of the input point cloud. This allows to compare against G-PCC [3] and V-PCC [5] using their respective standard configurations. We exclude YOGA [17], as such a set of configurations is not publicly available.

Selected rate-distortion curves are presented in Fig. 8 and corresponding Bjøntegaard-Deltas to G-PCC [3] are shown in Tab. 1. We provide further rate-distortion curves in the supplementary material. Note that our approach allows for significantly higher geometric quality according to D1-PSNR, thus small overlap of the resulting rate-distortion curves in the y-axis. This renders the Bjøntegaard-rate Delta ∆ r for the D1-PSNR metric less expressive on a subset of the tested point clouds.

Over all tested point clouds, our approach shows significantly higher geometric quality compared to G-PCC [3] and V-PCC [5]. The corresponding attribute rate-distortion curves according to the Y-PSNR metric show notable gains over

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-e22bfa6b-3305-4376-be97-81af6ed299bd.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-f587e353-6c95-4f3f-8a4e-799bd2496b35.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-5f6631a1-9992-4be0-9eaa-10912b144b63.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-d0203985-5bbf-46ad-b7b7-a79ba168193a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-26ca31da-bb2d-48c2-9a75-af09c1f3c727.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-a78d78ca-3974-4b81-bc52-dccd598676a7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-0d327abb-e098-459d-b119-95279f585f73.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-8047d866-fc0d-4a5a-9767-0d8b8f6a12ef.jpg)

Y-PSNR [dB]

Figure 8: Comparison of rate-distortion for geometry (D1-PSNR, top row) and attributes (Y-PSNR, bottom row) for selected point clouds from 8iVFBv2 [49] and MVUB [50].

Table 1: BD-Results for our approach compared to G-PCC [3].

|             | D1-PSNR   | D1-PSNR   | Y-PSNR   | Y-PSNR   | YUV-PSNR   | YUV-PSNR   | 1 - PCQM   | 1 - PCQM   |
|-------------|-----------|-----------|----------|----------|------------|------------|------------|------------|
| Point Cloud | ∆ r ↓     | ∆ D ↑     | ∆ r ↓    | ∆ D ↑    | ∆ r ↓      | ∆ D ↑      | ∆ r ↓      | ∆ D ↑      |
| longdress   | -82.5%    | 7.49 dB   | -3.6%    | 1.07 dB  | 22.6%      | 0.10 dB    | -13.7%     | 0.0016     |
| soldier     | -82.6%    | 6.78 dB   | -39.8%   | 2.94 dB  | -24.3%     | 1.70 dB    | -43.0%     | 0.0063     |
| loot        | -82.5%    | 7.49 dB   | -45.8%   | 3.28 dB  | -37.4%     | 1.89 dB    | -51.8%     | 0.0060     |
| redandblack | -60.4%    | 5.63 dB   | -19.5%   | 1.24 dB  | -3.1%      | 0.59 dB    | -17.6%     | 0.0023     |
| andrew9     | -72.8%    | 5.16 dB   | -16.3%   | 0.38 dB  | 0.4%       | -0.02 dB   | -19.3%     | 0.0030     |
| david9      | 79.0%     | 6.91 dB   | -48.2%   | 2.70 dB  | -42.6%     | 2.20 dB    | -52.3%     | 0.0069     |
| phil9       | -70.7%    | 5.05 dB   | -34.2%   | 1.63 dB  | -26.6%     | 1.05 dB    | -25.1%     | 0.0071     |
| sarah9      | -71.1%    | 6.68 dB   | -36.1%   | 2.90 dB  | -33.8%     | 2.36 dB    | -13.1%     | 0.0039     |

G-PCC, but underperform when compared against V-PCC. Considering attribute quality according to the YUV-PSNR metric, Tab. 1 shows similar results, with the exception of sequences with complex textures, i.e longdress , redandblack and andrew9 . For the Bjøntegaard-rate metrics ∆ r and ∆ D of the joint metric PCQM [48], we see high potential for rate-savings in the fixed configuration range compared to G-PCC [3].

## 4.4 Latency Analysis

Using a single, joint model in our approach allows to significantly simplify the encoding process, only requiring a single pass through the encoder. This stands in contrast to the standard procedure, which requires projection of the attributes to the reconstructed geometry during attribute encoding, as it is done in G-PCC [3], V-PCC [5] and YOGA [17]. Consequently, the encoding latency of our model ( 0 75 . s ) in Table 2 is more then one order of magnitude lower then the one for the implementations of traditional codecs [3, 5]. When comparing to the learned approach in YOGA [17], we also notice a lower encoding latency by one order of magnitude compared to the results reported in the paper. We note that the original paper utilizes different hardware, however, the encoding speed advantage of our approach, due to the single pass, is apparent.

Considering the decoding latency, we find that our method is competitive with respect to the related approaches, although all methods may require further optimization to allow application in practical use cases. Note, however, that the latency for all approaches significantly depends on the content (i.e. number of points) and the quality configuration. A more detailed, per-content and per-configuration analysis is given in the appendix.

Figure 9: Renders of point cloud soldier [49] at the same rate.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-9941d384-1388-4c23-a8c1-96b815622925.jpg)

Table 2: Comparison of complexity of related works compressing the test point clouds of the 8i dataset [49]. † indicates results from the original paper [17] on different hardware.

| Method   | Enc. time   | Dec. time   | Model size   |
|----------|-------------|-------------|--------------|
| G-PCC    | 11.24s      | 1.51s       | n/a          |
| V-PCC    | 57.29s      | 1.76s       | n/a          |
| YOGA     | † 8.15s     | † 3.23s     | † 169.5 MB   |
| Ours     | 0.75s       | 1.84s       | 120.1 MB     |

## 4.5 Showcase: View-dependent Compression

Recall that conditioning our model on a quality map ( g q , ) instead of a scalar quality parameter allows to perform region-of-interest rate control. To showcase this capability, we derive a view-dependent compression of a point cloud given some arbitrary user viewing angle. This is achieved by instantiating quality maps using a gradient along the viewing direction, assigning higher qualities to the front of the point cloud which decreases to the back. For comparison, we consistently apply the same procedure using a uniform quality map in Fig. 10a and a quality map according to the step function in Fig. 10c. For all compressed point clouds in Fig. 10, the geometry quality map q ( G ) is instantiated following an identical trend. The reconstructions are rendered from the given viewing-angle and standard image metrics, i.e. SSIM and the PSNR of the MSE in the YUV color space, are computed. Considering the image metrics of the renders, we find some potential for reducing the rate by locally dictating the quality q in both image metrics, i.e., PSNR and SSIM.

## 5 Conclusion and Outlook

We presented a method for training a single, joint compression model for point cloud geometry and attributes. Key to our approach is extending variable-rate compression techniques to multimodal compression, allowing to gauge the model on a tradeoff between the triplet of the quality of both modalities and the rate. Equipped with this, we proposed a joint compression architecture operating on a sparse tensor representation. Our technique achieves comparable compression performance to the state of the art in learned point cloud compression [17] which, however, has significantly higher encoding latency as it requires geometry decoding at the sender to encode the attributes, while outperforming G-PCC [3] and further closing the gap to V-PCC [5]. Our adaptive, variable-rate model illustrates the local degrees of freedom and flexibility in the joint compression of geometry and attributes.

However, compressing point clouds according to content specific configurations allows to fully utilize our model capabilities, but requires computationally expensive grid search in the space of possible configuration pairs. In some use cases, this may result in resorting to fixed configuration sets as evaluated. Tackling a similar issue in video compression, per-title encoding [51, 52] offers more efficient methods for optimizing a configuration set on specific content.

Deriving meaningful quality maps to locally condition the point cloud geometry and attribute quality using our method is a promising extension of this work. Our evaluation gives a perspective on the capabilities of view-dependent

Figure 10: Rendered views of the point cloud redandblack for different quality maps and corresponding image metrics.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-ca2de9c1-ff97-4d3a-b482-c0194bbd387a.jpg)

and region-of-interest based compression. Possible directions are incorporating per-point importance, e.g. through harnessing fixation maps of the content [53] or incorporating task-specific optimization [38].

## References

- [1] Sueyoon Lee, Irene Viola, Silvia Rossi, Zhirui Guo, Ignacio Reimat, Kinga Ławicka, Alina Striner, and Pablo Cesar. Designing and evaluating a vr lobby for a socially enriching remote opera watching experience. Transactions on Visualization and Computer Graphics , 2024.
- [2] Irene Viola, Jack Jansen, Shishir Subramanyam, Ignacio Reimat, and Pablo Cesar. Vr2gather: A collaborative, social virtual reality system for adaptive, multiparty real-time communication. IEEE MultiMedia , 30(2):48-59, 2023.
- [3] ISO/IEC JTC 1/SC 29. ISO/IEC 23090-9:2023, Information technology - Coded representation of immersive media - Part 9: Geometry-based point cloud compression . ISO/IEC, 2023.
- [4] Ricardo L De Queiroz and Philip A Chou. Compression of 3d point clouds using a region-adaptive hierarchical transform. IEEE Transactions on Image Processing , 25(8):3947-3956, 2016. doi:10.1109/TIP.2016.2575005.
- [5] ISO/IEC JTC 1/SC 29. ISO/IEC 23090-5:2021, Information technology - Coded representation of immersive media - Part 5: Visual volumetric video-based coding (V3C) and video-based point cloud compression (V-PCC) . ISO/IEC, 2021.
- [6] Johannes Ballé, Valero Laparra, and Eero P Simoncelli. End-to-end optimization of nonlinear transform codes for perceptual quality. In Picture Coding Symposium (PCS) , pages 1-5, 2016. doi:10.1109/PCS.2016.7906310.
- [7] Johannes Ballé, David Minnen, Saurabh Singh, Sung Jin Hwang, and Nick Johnston. Variational image compression with a scale hyperprior. arXiv preprint arXiv:1802.01436 , 2018.
- [8] David Minnen, Johannes Ballé, and George D Toderici. Joint autoregressive and hierarchical priors for learned image compression. Advances in Neural Information Processing Systems (NeurIPS) , 31, 2018.
- [9] Eirikur Agustsson, David Minnen, Nick Johnston, Johannes Balle, Sung Jin Hwang, and George Toderici. Scalespace flow for end-to-end optimized video compression. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR) , pages 8503-8512, 2020.

| [10]   | Jianqiang Wang, Hao Zhu, Haojie Liu, and Zhan Ma. Lossy point cloud geometry compression via end-to-end learning. IEEE Transactions on Circuits and Systems for Video Technology , 31(12):4909-4923, 2021.                                                                                                               |
|--------|--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [11]   | Jianqiang Wang, Dandan Ding, Zhu Li, and Zhan Ma. Multiscale point cloud geometry compression. In Data Compression Conference (DCC) , pages 73-82, 2021.                                                                                                                                                                 |
| [12]   | Jianqiang Wang, Dandan Ding, Zhu Li, Xiaoxing Feng, Chuntong Cao, and Zhan Ma. Sparse tensor-based multiscale representation for point cloud geometry compression. IEEE Transactions on Pattern Analysis and Machine Intelligence , 2022.                                                                                |
| [13]   | Maurice Quach, Giuseppe Valenzise, and Frederic Dufaux. Learning convolutional transforms for lossy point cloud geometry compression. In IEEE International Conference on Image Processing (ICIP) , pages 4320-4324, 2019.                                                                                               |
| [14]   | Maurice Quach, Giuseppe Valenzise, and Frederic Dufaux. Improved deep point cloud geometry compression. In IEEE International Workshop on Multimedia Signal Processing (MMSP) , pages 1-6, 2020.                                                                                                                         |
| [15]   | Evangelos Alexiou, Kuan Tung, and Touradj Ebrahimi. Towards neural network approaches for point cloud compression. In Applications of digital image processing XLIII , volume 11510, pages 18-37. SPIE, 2020.                                                                                                            |
| [16]   | André F.R. Guarda, Nuno M.M. Rodrigues, and Fernando Pereira. Point cloud coding: Adopting a deep learning-based approach. In Picture Coding Symposium (PCS) , pages 1-5, 2019.                                                                                                                                          |
| [17]   | Junteng Zhang, Tong Chen, Dandan Ding, and Zhan Ma. Yoga: Yet another geometry-based point cloud compressor. In Proceedings of the 31st ACM International Conference on Multimedia , pages 9070-9081, 2023. doi:10.1145/3581783.3613847.                                                                                 |
| [18]   | Yun He, Xinlin Ren, Danhang Tang, Yinda Zhang, Xiangyang Xue, and Yanwei Fu. Density-preserving deep point cloud compression. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR) , pages 2323-2332, 2022.                                                                       |
| [19]   | Rufael Mekuria, Kees Blom, and Pablo Cesar. Design, implementation, and evaluation of a point cloud codec for tele-immersive video. IEEE Transactions on Circuits and Systems for Video Technology , 27(4):828-842, 2016.                                                                                                |
| [20]   | Sebastian Schwarz, Marius Preda, Vittorio Baroncini, Madhukar Budagavi, Pablo Cesar, Philip A Chou, Robert A Cohen, Maja Krivoku´a, c Sébastien Lasserre, Zhu Li, et al. Emerging mpeg standards for point cloud compression. IEEE Journal on Emerging and Selected Topics in Circuits and Systems , 9(1):133-148, 2018. |
| [21]   | Jianqiang Wang, Dandan Ding, Hao Chen, and Zhan Ma. Dynamic point cloud geometry compression using multiscale inter conditional coding. arXiv preprint arXiv:2301.12165 , 2023.                                                                                                                                          |
| [22]   | Dandan Ding, Junjie Wang, Guangkun Zhen, Debargha Mukherjee, Urvang Joshi, and Zhan Ma. Neural adaptive loop filtering for video coding: Exploring multi-hypothesis sample refinement. IEEE Transactions on Circuits and Systems for Video Technology , 2023.                                                            |
| [23]   | Tsung-Yi Lin, Priya Goyal, Ross Girshick, Kaiming He, and Piotr Dollar. Focal loss for dense object detection. In Proceedings of the IEEE International Conference on Computer Vision (ICCV) , 2017.                                                                                                                     |
| [24]   | Andrew Brock, Theodore Lim, James MRitchie, and Nick Weston. Generative and discriminative voxel modeling with convolutional neural networks. arXiv preprint arXiv:1608.04236 , 2016.                                                                                                                                    |
| [25]   | Christopher Choy, JunYoung Gwak, and Silvio Savarese. 4d spatio-temporal convnets: Minkowski convolutional neural networks. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR) , 2019.                                                                                          |
| [26]   | Cha Zhang, Dinei Florencio, and Charles Loop. Point cloud attribute compression with graph transform. In IEEE International Conference on Image Processing (ICIP) , pages 2066-2070, 2014.                                                                                                                               |
| [27]   | Guangchi Fang, Qingyong Hu, Hanyun Wang, Yiling Xu, and Yulan Guo. 3dac: Learning attribute compression for point clouds. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR) , pages 14819-14828, 2022.                                                                         |
| [28]   | ISO/IEC JTC 1/SC 29. ISO/IEC 23009-1:2022, Information technology -Dynamic adaptive streaming over HTTP (DASH) -Part 1: Media presentation description and segment formats . ISO/IEC, 2022.                                                                                                                              |
| [29]   | Yadong Lu, Yinhao Zhu, Yang Yang, Amir Said, and Taco S Cohen. Progressive neural image compression with nested quantization and latent ordering. In IEEE International Conference on Image Processing (ICIP) , pages 539-543. IEEE, 2021.                                                                               |
| [30]   | Jae-Han Lee, Seungmin Jeon, Kwang Pyo Choi, Youngo Park, and Chang-Su Kim. Dpict: Deep progressive image compression using trit-planes. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR) , pages 16113-16122, 2022.                                                           |

| [31]   | Shaohui Li, Han Li, Wenrui Dai, Chenglin Li, Junni Zou, and Hongkai Xiong. Learned progressive image compression with dead-zone quantizers. IEEE Transactions on Circuits and Systems for Video Technology , 2022.                                                                                                                                                   |
|--------|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [32]   | Chunlei Cai, Li Chen, Xiaoyun Zhang, Guo Lu, and Zhiyong Gao. A novel deep progressive image compression framework. In Picture Coding Symposium (PCS) , pages 1-5, 2019.                                                                                                                                                                                             |
| [33]   | Enmao Diao, Jie Ding, and Vahid Tarokh. Drasic: Distributed recurrent autoencoder for scalable image compres- sion. In Data Compression Conference (DCC) , pages 3-12. IEEE, 2020.                                                                                                                                                                                   |
| [34]   | Khawar Islam, L Minh Dang, Sujin Lee, and Hyeonjoon Moon. Image compression with recurrent neural network and generalized divisive normalization. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR) Workshops , pages 1875-1879, 2021.                                                                                     |
| [35]   | Nick Johnston, Damien Vincent, David Minnen, Michele Covell, Saurabh Singh, Troy Chinen, Sung Jin Hwang, Joel Shor, and George Toderici. Improved lossy image compression with priming and spatially adaptive bit rates for recurrent networks. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR) , pages 4385-4393, 2018. |
| [36]   | Ali Hojjat, Janek Haberer, and Olaf Landsiedel. Progdtd: Progressive learned image compression with double-tail- drop training. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR) Workshops , pages 1130-1139, June 2023.                                                                                                  |
| [37]   | Ze Cui, Jing Wang, Shangyin Gao, Tiansheng Guo, Yihui Feng, and Bo Bai. Asymmetric gained deep image compression with continuous rate adaptation. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR) , pages 10532-10541, 2021.                                                                                             |
| [38]   | Myungseo Song, Jinyoung Choi, and Bohyung Han. Variable-rate deep image compression through spatially- adaptive feature transform. In Proceedings of the IEEE/CVF International Conference on Computer Vision (ICCV) pages 2380-2389, 2021.                                                                                                                          |
| [39]   | Zhihao Duan, Ming Lu, Jack Ma, Yuning Huang, Zhan Ma, and Fengqing Zhu. Qarv: Quantization-aware resnet vae for lossy image compression. IEEE Transactions on Pattern Analysis and Machine Intelligence , 2023.                                                                                                                                                      |
| [40]   | Shilv Cai, Zhijun Zhang, Liqun Chen, Luxin Yan, Sheng Zhong, and Xu Zou. High-fidelity variable-rate image compression via invertible activation transformation. arXiv preprint arXiv:2209.05054 , 2022.                                                                                                                                                             |
| [41]   | Shishir Subramanyam, Irene Viola, Alan Hanjalic, and Pablo Cesar. User centered adaptive streaming of dynamic point clouds with low complexity tiling. In Proceedings of the 28th ACM International Conference on Multimedia pages 3669-3677, 2020.                                                                                                                  |
| [42]   | Michael Rudolph and Amr Rizk. View-adaptive streaming of point cloud scenes through combined decompo- sition and video-based coding. In Proceedings of the 1st International Workshop on Advances in Point Cloud Compression, Processing and Analysis , pages 41-49, 2022.                                                                                           |
| [43]   | JunYoung Gwak, Christopher BChoy, and Silvio Savarese. Generative sparse detection networks for 3d single-shot object detection. In European Conference on Computer Vision (ECCV) , 2020.                                                                                                                                                                            |
| [44]   | Guillaume Gautier, Alexandre Mercat, Louis Fréneau, Mikko Pitkänen, and Jarno Vanne. Uvg-vpc: Voxelized point cloud dataset for visual volumetric video-based coding. In 15th International Conference on Quality of Multimedia Experience (QoMEX) , pages 244-247. IEEE, 2023.                                                                                      |
| [45]   | Diederik P Kingma and Jimmy Ba. Adam: A method for stochastic optimization. arXiv preprint arXiv:1412.6980 2014.                                                                                                                                                                                                                                                     |
| [46]   | Dong Tian, Hideaki Ochimizu, Chen Feng, Robert Cohen, and Anthony Vetro. Geometric distortion metrics for point cloud compression. In IEEE International Conference on Image Processing (ICIP) , pages 3460-3464, 2017.                                                                                                                                              |
| [47]   | Rufael Mekuria, Sebastien Laserre, and Christian Tulvan. Performance assessment of point cloud compression. In IEEE Visual Communications and Image Processing (VCIP) , pages 1-4, 2017.                                                                                                                                                                             |
| [48]   | Gabriel Meynet, Yana Nehmé, Julie Digne, and Guillaume Lavoué. Pcqm: A full-reference quality metric for colored 3d point clouds. In Twelfth International Conference on Quality of Multimedia Experience (QoMEX) pages 1-6. IEEE, 2020.                                                                                                                             |
| [49]   | Eugene d'Eon, Bob Harrison, Taos Myers, and Philip A Chou. 8i voxelized full bodies-a voxelized point cloud dataset. ISO/IEC JTC1/SC29 Joint WG11/WG1 (MPEG/JPEG) input document WG11M40059/WG1M74006 , 7 (8):11, 2017.                                                                                                                                              |
| [50]   | Charles Loop, Qin Cai, S Orts Escolano, and Philip A Chou. Microsoft voxelized upper bodies-a voxelized point cloud dataset. ISO/IEC JTC1/SC29 Joint WG11/WG1 (MPEG/JPEG) input document m38673 M , 72012:2016,                                                                                                                                                      |

- [51] Jan De Cock, Zhi Li, Megha Manohara, and Anne Aaron. Complexity-based consistent-quality encoding in the cloud. In IEEE International Conference on Image Processing (ICIP) , pages 1484-1488, 2016.
- [52] Hadi Amirpour, Christian Timmerer, and Mohammad Ghanbari. Pstr: Per-title encoding using spatio-temporal resolutions. In IEEE International Conference on Multimedia and Expo (ICME) , pages 1-6, 2021.
- [53] Xuemei Zhou, Irene Viola, Evangelos Alexiou, Jack Jansen, and Pablo Cesar. Qava-dpc: Eye-tracking based quality assessment and visual attention dataset for dynamic point cloud in 6 dof. In International Symposium on Mixed and Augmented Reality (ISMAR) , pages 69-78. IEEE, 2023.
- [54] Johannes Ballé, Nick Johnston, and David Minnen. Integer networks for data compression with latent-variable models. In International Conference on Learning Representations (ICLR) , 2018.

Figure 11: Implementation details for our architecture.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-82d128d5-3772-430d-ba2e-d692b22a4220.jpg)

Table 3: Tested point clouds

| Dataset      | Sequence                           | Frame              | Resolution          |
|--------------|------------------------------------|--------------------|---------------------|
| 8iVFBv2 [49] | longdress soldier loot redandblack | 1300 690 1200 1550 | 1024 1024 1024 1024 |
| MVUB [50]    | phil9 sarah9 david9 andrew9        | 0 0 0 0            | 512 512 512 512     |

## A Architecture

We extend the visualization in Fig. 3, presenting a more detailed version in Fig. 11 with implementation channel sizes for reproducability. All convolution layers, denoted Conv are three-dimensional sparse convolutions, for which we use MinkowskiEngine [25]. For layers with the extension ↑ , Generative Transposed Convolutions [43] are used, allowing to generate geometry locations and corresponding features during upsampling. Note that the hyperdecoder h s is prone to non-deterministic behaviour [54], which arises when using MinkowskiEngine [43]. To this end, we reduce the kernel size for generative convolutions in the hyperdecoder to 2 , as we found that this allows to counteract non-determinism in combination with sorting point locations before and after the convolution.

The geometry representation g ( 3 ) at the entropy bottleneck is encoded using G-PCC in lossless geometry mode, which is common practice [10, 11, 17]. The base geometry g ( 5 ) required to decode the hyperpriors z is derived through dyadic downsampling of g ( 3 ) .

In addition, the number of voxels to keep at each pruning stage in the UP block (cf. 3.4.2) is appended as side-information in the header of the bitstream.

## B Additions to the Evaluation

## B.1 Test Data

The point cloud frames select for testing are presented in Tab. 3, consisting of the 8iVFBv2 [49] and MVUB dataset [50]. While point clouds in the 8iVFBv2 dataset [49] consist of watertight surfaces and high quality geometry and textures at voxel resolution 1024 , point clouds in MVUB dataset [50] are frontal scans of human upper bodies containing visible (geometry) noise at voxel resolution 512 . Incorporating both datasets in the evaluation allowed us to assess the performance of our model over broad input conditions.

Table 4: Encoding latency measured on a AMD EPYC 7542 using a NVIDIA RTX4090 graphics card over all configurations on the point clouds of the 8i dataset [49] listed in the paper. † indicates that the results are obtained from the original paper [17], as the model is not open-sourced.

| Method   | Encoding Latency [s]   | Encoding Latency [s]   | Encoding Latency [s]   | Encoding Latency [s]   | Encoding Latency [s]   |
|----------|------------------------|------------------------|------------------------|------------------------|------------------------|
|          | Loot                   | Redandblack            | Soldier                | Longdress              | Average                |
| G-PCC r1 | 6 . 90 ± 0 . 005       | 6 . 56 ± 0 . 007       | 9 . 49 ± 0 . 008       | 7 . 30 ± 0 . 014       | 7 . 57                 |
| G-PCC r2 | 5 . 85 ± 0 . 006       | 5 . 52 ± 0 . 008       | 8 . 03 ± 0 . 008       | 6 . 21 ± 0 . 011       | 6 . 41                 |
| G-PCC r3 | 8 . 50 ± 0 . 008       | 8 . 03 ± 0 . 011       | 11 . 64 ± 0 . 011      | 7 . 34 ± 0 . 009       | 9 . 30                 |
| G-PCC r4 | 13 . 64 ± 0 . 014      | 12 . 87 ± 0 . 014      | 18 . 60 ± 0 . 024      | 9 . 03 ± 0 . 019       | 14 . 90                |
| V-PCC    | 55 . 08 ± 0 . 513      | 56 . 14 ± 0 . 498      | 65 . 60 ± 0 . 732      | 52 . 3 ± 0 . 772       | 57 . 28                |
| YOGA     | n/a                    | n/a                    | n/a                    | n/a                    | 8 . 15                 |
| Ours     | 0 . 66 ± 0 . 006       | 0 . 66 ± 0 . 008       | 0 . 95 ± 0 . 011       | 0 . 74 ± 0 . 007       | 0 . 75                 |

Table 5: Comparison of decoding latency, measured on a AMD EPYC 7542 using a NVIDIA RTX4090 graphics card over all configurations on the point clouds of the 8i dataset [49] listed in the paper. † indicates that the results are obtained from the original paper [17], as the model is not open-sourced.

| Method   | Decoding Latency [s]   | Decoding Latency [s]   | Decoding Latency [s]   | Decoding Latency [s]   | Decoding Latency [s]   |
|----------|------------------------|------------------------|------------------------|------------------------|------------------------|
|          | Loot                   | Redandblack            | Soldier                | Longdress              | Average                |
| G-PCC r1 | 0 . 08 ± 0 . 002       | 0 . 07 ± 0 . 002       | 0 . 10 ± 0 . 001       | 0 . 08 ± 0 . 001       | 0 . 08                 |
| G-PCC r2 | 0 . 28 ± 0 . 002       | 0 . 27 ± 0 . 002       | 0 . 38 ± 0 . 002       | 0 . 30 ± 0 . 003       | 0 . 31                 |
| G-PCC r3 | 1 . 06 ± 0 . 004       | 1 . 01 ± 0 . 004       | 1 . 45 ± 0 . 005       | 1 . 14 ± 0 . 005       | 1 . 16                 |
| G-PCC r4 | 2 . 29 ± 0 . 005       | 2 . 18 ± 0 . 007       | 3 . 14 ± 0 . 009       | 2 . 46 ± 0 . 008       | 2 . 52                 |
| V-PCC    | 1 . 60 ± 0 . 021       | 1 . 59 ± 0 . 018       | 2 . 12 ± 0 . 023       | 1 . 71 ± 0 . 022       | 1 . 76                 |
| YOGA     | n/a                    | n/a                    | n/a                    | n/a                    | 3 . 23                 |
| Ours     | 1 . 66 ± 0 . 009       | 1 . 71 ± 0 . 009       | 2 . 23 ± 0 . 009       | 1 . 81 ± 0 . 009       | 1 . 84                 |

## B.2 Extended Results for the Latency

Complementary to the results presented in Table 2 in Sec. 4.4, we showcase per content results of the encoding latency in Table 4 and decoding latency in Table 5. Since YOGA [17] is not open-sourced at the time of our evaluation, we can only report the averaged latencies from their results, which were measured on different hardware. For all other approaches, we use the configurations for the generation of the pareto-fronts in Sec. 4.2 and measure the times for each encoding and decoding. We notice that the encoding and decoding latency of G-PCC [3] heavily relies on the quality configuration for the geometry and thus split the results for each tested geometry setting. It becomes apparent that the latency highly depends on the content, more specifically on the number of points.

Overall, we observe similar behaviour as presented in Sec. 4.4: Our method allows significantly faster encoding times as other methods, which can be mainly attributed to the fact that the geometry does not need to be reconstructed for attribute re-projection during the encoding process. For decoding, we find our method performing comparable to related work.

## B.3 Complementary Rate Distortion Results for Fixed-Configuration Encoding

We showcase the rate-distortion results on all point clouds according to the evaluation procedure described in Sec. 4.3, as only selected curves were shown in the main document. Fig. 12 shows results for the 8iVFBv2 dataset [49] and Fig. 13 shows results for the MVUB dataset [50]. Overall, we find increased performance for the geometry metric D1-PSNR over V-PCC [5] and G-PCC [3]. For attribute quality our approach achieves a better rate-distortion tradeoff then G-PCC but falls short of the attribute qualities achieved by V-PCC. Especially the point cloud andrew9 , which shows a checkerboard texture pattern on the shirt, is hard to compress using our learned approach. Finally, for the joint metric PCQM [48], we find that our approach consistently achieves higher quality at the given rates compared to G-PCC [3]. For higher rate configurations, our method can deliver comparable and sometimes better quality as V-PCC [5], but falls of when increasing the compression ratio.

Additionally, we notice that the range of compressed bitrates for our model is small compared to the rate-distortion rates achievable by both traditional methods.

Figure 12: Rate distortion results on the 8iVFBv2 [49] dataset.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-4c3bf68b-d42e-410a-a4bc-ac88c40d5363.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-7ee0975f-9c03-452f-b64f-b763db51ef2e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-c0d8ba42-2f34-42f0-8697-409380e59bf9.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-89d1c6a3-7f1c-48c6-a532-453751edaa12.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-83ddbf23-968f-411b-b226-624a7cf7e134.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-0a6a3000-7742-409c-b22d-0d1db98b69b9.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-5bce410c-85eb-4a5e-b1cf-a44b0093e341.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-0f9290cf-84ed-43b0-b662-1982be24ba47.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-37fc587c-e44f-4921-a1d5-23aa8db82726.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-b7c89d4a-b253-43a5-a9ed-d2721ab254bc.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-0aeb99ed-1148-492b-a6ea-9c1a0fdd98c7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-b6a188de-3a3e-4d23-805d-d40459e27dc7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-0ded363b-61fe-4f50-be24-fc1d402e712a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-a2a4effc-9619-4c27-8c18-445e4b5961ee.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-72206abd-d88e-4a52-8e9b-7845a19493a1.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-4c35a65d-2e95-4e5f-9b83-1eee25957a46.jpg)

Figure 13: Rate distortion results on the MVUB [50] dataset.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-697adf82-eee3-4e53-b518-b3f09ffd09ce.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-55d02dc6-bd9a-44d5-bd17-7b3bc376eada.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-5c5b6d85-1043-4f25-bc8b-fd12f0bb4b75.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-ee8aea3a-86f9-4ee6-ab23-253ef328d1e4.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-bfa33df3-0da3-4649-a793-f04dc93e8888.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-21aa3637-a721-474c-b4d6-a9d3683e362c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-1b450ede-8a23-4db4-969d-365e5b98df0a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-ae269ff3-cc89-4e05-9460-4386e2aa0fef.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-912cbdd5-0ebd-4023-8a31-6cf9cfcbc8d4.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-a1b2bb6d-fb28-4c71-8ba1-b9a7b91e40a8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-ee422bde-2a4c-4db0-bd3e-fac5d15ad0f6.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-a3dd5656-9f91-4cb8-b765-5e17d8dadca7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-0132eb48-97cb-4db6-97e1-af4b24fb3917.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-5d41397d-421a-4dd7-b6f5-03e571cf958b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-4c3ebb96-9d02-4d86-97f8-e5990f99f358.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00599v1-figure-f74a44ec-56e6-4f3a-8041-bdff2b172e7e.jpg)
Michael Rudolph Aron Riemenschneider Amr Rizk
Learned Compression of Point Cloud Geometry and Attributes in a Single Model through Multimodal Rate-Control
Point cloud compression is essential to experience volumetric multimedia as
it drastically reduces the required streaming data rates. Point attributes,
specifically colors, extend the challenge of lossy compression beyond geometric
representation to achieving joint reconstruction of texture and geometry.
State-of-the-art methods separate geometry and attributes to compress them
individually. This comes at a computational cost, requiring an encoder and a
decoder for each modality. Additionally, as attribute compression methods
require the same geometry for encoding and decoding, the encoder emulates the
decoder-side geometry reconstruction as an input step to project and compress
the attributes. In this work, we propose to learn joint compression of geometry
and attributes using a single, adaptive autoencoder model, embedding both
modalities into a unified latent space which is then entropy encoded. Key to
the technique is to replace the search for trade-offs between rate, attribute
quality and geometry quality, through conditioning the model on the desired
qualities of both modalities, bypassing the need for training model ensembles.
To differentiate important point cloud regions during encoding or to allow
view-dependent compression for user-centered streaming, conditioning is
pointwise, which allows for local quality and rate variation. Our evaluation
shows comparable performance to state-of-the-art compression methods for
geometry and attributes, while reducing complexity compared to related
compression methods.