Universal Rate-Distortion-Perception Representations for Lossy Compression

1 Introduction . Unlike in lossless compression , the decoder in a lossy compression system has flexibility in how to reconstruct the source . Conventionally , some measure of distortion such as mean squared error , PSNR or SSIM/MS-SSIM [ 36 , 37 ] is used as a quality measure . Accordingly , lossy compression algorithms are analyzed through rate-distortion theory , wherein the objective is to minimize the amount of distortion for a given rate . However , it has been observed that low distortion is not necessarily synonymous with high perceptual quality ; indeed , deep learning based image compression has inspired works in which authors have noted that increased perceptual quality may come at the cost of increased distortion [ 1 , 4 ] . This culminated in the work of Blau & Michaeli [ 5 ] who propose the rate-distortion-perception theoretical framework . The main idea was to introduce a third perception axis which more closely mimics what humans would deem to be visually pleasing . Unlike distortion , judgement of perceptual quality is taken to 35th Conference on Neural Information Processing Systems ( NeurIPS 2021 ) . be inherently no-reference . The mathematical proxy for perceptual quality then comes in the form of a divergence between the source and the reconstruction distributions , motivated by the idea that perfect perceptual quality is achieved when they are identical . Leveraging generative adversarial networks [ 11 ] in the training procedure has made such a task possible for complex data-driven settings with efficacy even at very low rates [ 33 ] . Naturally , this induces a tradeoff between optimizing for perceptual quality and optimizing for distortion . But in designing a lossy compression system , one may wonder where exactly this tradeoff lies : is the objective tightly coupled with optimizing the representations generated by the encoder , or can most of this tradeoff be achieved by simply changing the decoding scheme ? Our contributions are as follows . We define the notion of universal representations which are generated by a fixed encoding scheme for the purpose of operating at multiple perception-distortion tradeoff points attained by varying the decoder . We then prove a coding theorem establishing the relationship between this operational definition and an information universal rate-distortionperception function . Under MSE distortion loss , we study this function for the special case of the Gaussian distribution and show that the penalty in fixing the representation map with fixed rate can be small in many interesting regimes . For general distributions , we characterize the achievable distortionperception region with respect to an arbitrary representation and establish a certain approximate universality property . We then turn to study how the operational tradeoffs achieved by machine learning models on image compression under a fixed encoder compared to varying encoders . Our results suggest that there is not much loss in reusing encoders trained for a specific point on the distortion-perception tradeoff across other points . The practical implication of this is to reduce the number of models to be trained within deep-learning enhanced compression systems . Building on [ 30 , 31 ] , one of the key steps in our techniques is the assumption of common randomness between the sender and receiver which will turn out to reduce the coding cost . Throughout this work , we focus on the scenario where a rate is fixed in advance . We address the scenario when the rate is changed in the supplementary . 2 Related Works . Image quality measures include full-reference metrics ( which require a ground truth as reference ) , or no-reference metrics ( which only use statistical features of inputs ) . Common full-reference metrics include MSE , SSIM/MS-SSIM [ 36 , 37 ] , PSNR or deep feature based distances [ 16 , 41 ] . No-reference metrics include BRISQUE/NIQE/PIQE [ 23 , 24 , 35 ] and Fréchet Inception Distance [ 14 ] . Roughly speaking , one can consider the former set to be distortion measures and the latter set to be perception measures in the rate-distortion-perception framework . Since GANs capable of synthesizing highly realistic samples have emerged , using trained discriminators as a proxy for perceptual quality in deep learning based systems has also been explored [ 17 ] . This idea is principled as various GAN objectives can be interpreted as estimating particular statistical distances [ 3 , 25 , 26 ] . Rate-distortion theory has long served as a theoretical foundation for lossy compression [ 7 ] . Within machine learning , variations of rate-distortion theory have been introduced to address representation learning [ 2 , 6 , 34 ] , wherein a central task is to extract useful information from data on some sort of budget , and also in the related field of generative modelling [ 15 ] . On the other hand , distributionpreserving lossy compression problems have also been studied in classical information theory literature [ 27 , 28 , 40 ] . More recently , in an effort to reduce blurriness and other artifacts , machine learning research in lossy compression has attempted to incorporate GAN regularization into compressive autoencoders [ 1 , 5 , 33 ] , which were traditionally optimized only for distortion loss [ 21 , 32 ] . This has led to highly successful data-driven models operating at very low rates , even for high-resolution images [ 22 ] . An earlier work of Blau & Michaeli [ 4 ] studied only the perception-distortion tradeoff within deep learning enhanced image restoration using GANs . This idea was then incorporated with distributionpreserving lossy compression [ 33 ] to study the rate-distortion-perception tradeoff in full generality [ 5 ] . The work most similar to ours is [ 39 ] , who observe that an optimal encoder for the “ classic ” ratedistortion function is also optimal for perfect perceptual compression at twice the distortion . Our work investigates the intermediate regime and also includes common randomness as a modelling assumption , which in principle allows us to achieve perfect perceptual quality at lower than twice the distortion . The concurrent work [ 10 ] also establishes the achievable distortion-perception region as in our Theorem 4 and provides a geometric interpretation of the optimal interpolator in Wasserstein space . 3 Rate-Distortion-Perception Representations . The backbone of rate-distortion theory characterizes an ( operational ) objective expressing what can be achieved by encoders and decoders with a quantization bottleneck in terms of an information function which is more convenient to analyze . Let X ∼ pX be an information source to be compressed through quantization . The quality of the compressed source is measured by a distortion function ∆ : X × X → R≥0 satisfying ∆ ( x , x̂ ) = 0 if and only if x = x̂ . We distinguish between the one-shot scenario in which we compress one symbol at a time , and the asymptotic scenario in which we encode n i.i.d . samples from X jointly and analyze the behaviour as n→ ∞ . The minimum rate needed to meet the distortion constraint D on average is denoted by R∗ ( D ) in the one-shot setting and by R ( ∞ ) ( D ) in the asymptotic setting . These are studied through the information rate-distortion function R ( D ) = inf pX̂|X I ( X ; X̂ ) s.t . E [ ∆ ( X , X̂ ) ] ≤ D , ( 1 ) where I ( X ; X̂ ) is the mutual information between a source X and reconstruction X̂ . The principal result of rate-distortion theory states that R ( ∞ ) ( D ) = R ( D ) [ 7 ] . Furthermore , it is also possible to characterize R∗ ( D ) using R ( D ) as we will soon see . In light of the discussion on perceptual quality , the flexibility in distortion function is not necessarily a good method to capture how realistic the output may be perceived . To resolve this , Blau & Michaeli [ 5 ] introduce an additional constraint to match the distributions of X and X̂ in the form of a nonnegative divergence between probability measures d ( · , · ) satisfying d ( p , q ) = 0 if and only if p = q . The one-shot rate-distortion-perception function R∗ ( D , P ) and asymptotic rate-distortion-perception function R ( ∞ ) ( D , P ) are defined in the same fashion as their rate-distortion counterparts , which we will later make precise . Definition 1 ( iRDPF ) . The information rate-distortion-perception function for a source X is defined as R ( D , P ) = inf pX̂|X I ( X ; X̂ ) s.t . E [ ∆ ( X , X̂ ) ] ≤ D , d ( pX , pX̂ ) ≤ P. The strong functional representation lemma [ 19 , 31 ] establishes relationships between the operational and information functions : R ( ∞ ) ( D , P ) = R ( D , P ) , ( 2 ) R ( D , P ) ≤ R∗ ( D , P ) ≤ R ( D , P ) + log ( R ( D , P ) + 1 ) + 5 . ( 3 ) These results hold also for R ( D ) = R ( D , ∞ ) . We make note that they were developed under a more general set of constraints for which stochastic encoders and decoders with a shared source of randomness were used . In practice , the sender and receiver agree on a random seed beforehand to emulate this behaviour . 3.1 Gaussian Case . We now present the closed form expression of R ( D , P ) for a Gaussian source under MSE distortion and squared Wasserstein-2 perception losses ( see also Figure 1 ( a ) and Figure 1 ( b ) ) . Recall that the squared Wasserstein-2 distance is defined as W 22 ( pX , pX̂ ) = inf E [ ∥X − X̂∥2 ] , ( 4 ) where the infimum is over all joint distributions of ( X , X̂ ) with marginals pX and pX̂ . Let µX = E [ X ] and σ2X = E [ ∥X − µX∥2 ] . Theorem 1 . For a scalar Gaussian source X ∼ N ( µX , σ2X ) , the information rate-distortionperception function under squared error distortion and squared Wasserstein-2 distance is attained by some X̂ jointly Gaussian with X and is given by R ( D , P ) = 1 2 log σ2X ( σX− √ P ) 2 σ2X ( σX− √ P ) 2− ( σ 2 X + ( σX− √ P ) 2−D 2 ) 2 if √ P ≤ σX − √ |σ2X −D| , max { 12 log σ2X D , 0 } if √ P > σX − √ |σ2X −D| . When √ P > σX − √ |σ2X −D| , the perception constraint is inactive and R ( D , P ) = R ( D ) . The choice of W 22 ( · , · ) perception loss turns out to not be essential ; we show in the supplementary that R ( D , P ) can also be expressed under the KL-divergence . 3.2 Universal Representations . Whereas the RDP function is regarded as the minimal rate for which we can vary an encoder-decoder pair to meet any distortion and perception constraints ( D , P ) , the universal RDP ( uRDP ) function generalizes this to the case where we fix an encoder and allow only the decoder to adapt in order to meet multiple constraints ( D , P ) ∈ Θ . For example , one case of interest is when Θ is the set of all ( D , P ) pairs associated with a given rate along the iRDP function ; how much additional rate is needed if this is to be achieved by a fixed encoder , rather than varying it across each objective ? The hope is that the rate to use some fixed encoder across this set is not much larger than the rate to achieve any single point . As we will see for the Gaussian distribution , this is in fact the case in the asymptotic setting , and also approximately true in the one-shot setting . Below , we define the one-shot universal rate-distortion-perception function and the information universal rate-distortion-perception function , then establish a relationship between the two . In these definitions we assume X is a random variable and Θ is an arbitrary non-empty set of ( D , P ) pairs . Definition 2 ( ouRDPF ) . A Θ-universal encoder of rate R is said to exist if we can find random variable U , encoding function fU : X → CU and decoding functions gU , D , P : CU → X̂ , ( D , P ) ∈ Θ such that E [ ℓ ( fU ( X ) ) ] ≤ R , E [ ∆ ( X , X̂D , P ) ] ≤ D , d ( pX , pX̂D , P ) ≤ P , where CU is a uniquely decodable binary code specified by U , X̂D , P = gU , D , P ( fU ( X ) ) , and ℓ ( fU ( X ) ) denotes the length of binary codeword fU ( X ) . The random variable U acts as a shared source of randomness . The infimum of such R is called the one-shot universal rate-distortionperception function ( ouRDPF ) and denoted by R∗ ( Θ ) . When Θ = { ( D , P ) } , this specializes to the one-shot rate-distortion-perception function R∗ ( D , P ) . Definition 3 ( iuRDPF ) . Let Z be a representation of X ( i.e . generated by some random transform pZ|X ) . Let PZ|X ( Θ ) be the set of transforms pZ|X such that for each ( D , P ) ∈ Θ , there exists pX̂D , P |Z for which E [ ∆ ( X , X̂D , P ) ] ≤ D and d ( pX , pX̂D , P ) ≤ P , where X ↔ Z ↔ X̂D , P are assumed to form a Markov chain . Define R ( Θ ) = inf pZ|X∈PZ|X ( Θ ) I ( X ; Z ) . ( 5 ) We refer to this as the information universal rate-distortion-perception function ( iuRDPF ) and say that the random variable Z is a representation which is Θ-universal with respect to X . The conditional distributions pX̂D , P |Z induce stochastic mappings transforming the representations to reconstructions X̂D , P in order to meet specific ( D , P ) constraints . Note that we assume a shared source of stochasticity within the ouRDPF as a tool to prove the achievability of the iuRDPF , but not within the definition of the iuRDPF itself . Moreover , source X , reconstruction X̂D , P , representation Z , and random seed U are all allowed to be multivariate random variables . Theorem 2 . R ( Θ ) ≤ R∗ ( Θ ) ≤ R ( Θ ) + log ( R ( Θ ) + 1 ) + 5 . In practice , the overhead log ( R ( Θ ) + 1 ) + 5 either makes the upper bound an overestimate of R∗ ( Θ ) or is negligible compared toR ( Θ ) . This overhead vanishes completely in the asymptotic setting as we will show in the supplementary . We can therefore interpret R ( Θ ) as the rate required to meet an entire set Θ of constraints with the encoder fixed . Within the set Θ , it is clear that sup ( D , P ) ∈ΘR ( D , P ) characterizes the rate required to meet the most demanding constraint . Now define A ( Θ ) = R ( Θ ) − sup ( D , P ) ∈Θ R ( D , P ) , ( 6 ) which is the rate penalty incurred by meeting all constraints in Θ with the encoder fixed . Let Ω ( R ) = { ( D , P ) : R ( D , P ) ≤ R } . It is ideal to have A ( Ω ( R ) ) = 0 for each R so that achieving the entire tradeoff with a single encoder is essentially no more expensive than to achieve any single point on the tradeoff , thereby alleviating the need to design a host of encoders for different distortionperception objectives with respect to the same rate . One can also take the following alternative perspective . The proof of Theorem 2 shows that every representation Z can be generated from source X using an encoder of rate I ( X ; Z ) + o ( I ( X ; Z ) ) , and based on Z , the decoder can produce reconstruction X̂D , P by leveraging random seed U to simulate conditional distribution pX̂D , P |Z . Therefore , the problem of designing an encoder boils down to identifying a suitable representation . Given a representation Z ofX , we define the achievable distortion-perception region Ω ( pZ|X ) as the set of all ( D , P ) pairs for which there exists pX̂D , P |Z such that E [ ∆ ( X , X̂D , P ) ] ≤ D and d ( pX , pX̂D , P ) ≤ P . Intuitively , Ω ( pZ|X ) is the set of all possible distortion-perception constraints that can be met based on representation Z . If Ω ( pZ|X ) = Ω ( R ) for some representation Z with I ( X ; Z ) = R , then Z has the maximal achievable distortion-perception region in the sense that Ω ( pZ′|X ) ⊆ Ω ( pZ|X ) for any Z ′ with I ( X ; Z ′ ) ≤ R. In the supplementary material we establish mild regularity conditions for which the existence of such Z is equivalent to the aforementioned desired property A ( Ω ( R ) ) = 0 . We shall show that that this ideal scenario actually arises in the Gaussian case and an approximate version can be found more broadly . Theorem 3 . LetX ∼ N ( µX , σ2X ) be a scalar Gaussian source and assume MSE andW 22 ( · , · ) losses . Let Θ be any non-empty set of ( D , P ) pairs . Then A ( Θ ) = 0 . ( 7 ) Moreover , for any representation Z jointly Gaussian with X such that I ( X ; Z ) = sup ( D , P ) ∈Θ R ( D , P ) , ( 8 ) we have Θ ⊆ Ω ( pZ|X ) = Ω ( I ( X ; Z ) ) . ( 9 ) Next we consider a general source X ∼ pX and characterize the achievable distortion-perception region for an arbitrary representation Z under MSE loss . We then provide some evidence indicating that every reconstruction X̂D , P achieving some point ( D , P ) on the distortion-perception tradeoff for a given R likely has the property Ω ( pX̂D , P |X ) ≈ Ω ( R ) . Theorem 4 ( Approximate universality for general sources ) . Assume MSE loss and any perception measure d ( · , · ) . Let Z be any arbitrary representation of X . Then Ω ( pZ|X ) ⊆ { ( D , P ) : D ≥ E [ ∥X − X̃∥2 ] + inf pX̂ : d ( pX , pX̂ ) ≤P W 22 ( pX̃ , pX̂ ) } ⊆ cl ( Ω ( pZ|X ) ) , where X̃ = E [ X|Z ] is the reconstruction minimizing squared error distortion with X under the representation Z and cl ( · ) denotes set closure . In particular , the two extreme points ( D ( a ) , P ( a ) ) = ( E [ ∥X − X̃∥2 ] , d ( pX , pX̃ ) ) and ( D ( b ) , P ( b ) ) = ( E [ ∥X − X̃∥2 ] +W 22 ( pX̃ , pX ) , 0 ) are contained in cl ( Ω ( pZ|X ) ) . which suggests that Ω ( pX̂D1 , P1 |X ) is not much smaller than Ω ( R ) . Moreover , in this case we have D ( b ) ≤ 2E [ ∥X − X̃∥2 ] ≤ 2D1 , ( 11 ) which implies that ( 2D1 , 0 ) is dominated by extreme point ( D ( b ) , 0 ) and consequently must be contained in cl ( Ω ( pX̂D1 , P1 |X ) ) . Therefore , the optimal representation in the conventional ratedistortion sense can be leveraged to meet any perception constraint with no more than a two-fold increase in distortion . As a corollary , one recovers Theorem 2 in Blau & Michaeli ( R ( 2D , 0 ) ≤ R ( D , ∞ ) ) . Hence , the numerical connection between R ( 2D , 0 ) and R ( D , ∞ ) is a manifestation of the existence of approximately Ω ( I ( X ; Z ) ) -universal representations Z . These analyses motivate the study of practical constructions for which we seek to achieve multiple ( D , P ) pairs with a single encoder .

Summary: The rate-distortion-perception (RDP) function[5,23] is a generalization of the classical rate-distortion (RD) trade-off in lossy compression to also measure realism, which establishes a theoretical footing for generative image/video compression, a topic of active research[1,11,17,25].  This paper studies the RDP function under the constraint of a shared encoder for multiple D-P points, both from the theoretical side and experimentally. This is of significant practical interest, as detailed below.