## A note on the sum-rate-distortion function of some lossy source coding problems involving infinite-valued distortion functions

## Prakash Ishwar

Department of Electrical and Computer Engineering Boston University, Boston, MA 02215 USA e-mail: pi@bu.edu .

## Abstract

For a number of lossy source coding problems it is shown that even if the usual single-letter sum-rate-distortion expressions may become invalid for non-infinite distortion functions, they can be approached, to any desired accuracy, via the usual valid expressions for appropriately truncated finite versions of the distortion functions.

## Index Terms

Source coding, rate-distortion function, sum-rate, non-finite distortion function.

## I. INTRODUCTION

In a number lossy source coding problems, the minimum sum-rate needed to attain a target expected distortion of D or less - the operational sum-rate-distortion function - is characterized by an informationtheoretic rate-distortion function R D ( ) which is expressed in terms of an optimization problem of the following form:

$$R ( D ) \coloneqq \inf _ { p _ { X, Y, U, \widehat { X } } \in \mathcal { A } } \text{s.t. } E [ d ( X, Y, U, \widehat { X } ) ] \leq D } f ( p _ { X, Y, U, \widehat { X } } ),$$

where X ∈ X , Y ∈ Y , U ∈ U , and ̂ X ∈ X ̂ are, respectively, the source, side-information, auxiliary, and reconstruction random variables taking values in finite alphabets, p X,Y,U,X ̂ is their joint pmf, f is a finite linear combination of conditional mutual informations involving some or all the random variables, A is a finite set of marginal consistency 1 and Markov-chain constraints that the random variables need to satisfy, and d is a real-valued, nonnegative distortion function which does not depend on U . Examples include Shannon rate-distortion [1] and Gray-Leiner conditional rate-distortion [2]-[4] where auxiliary random variables are not needed, Wyner-Ziv rate-distortion [5], [6] which uses a single auxiliary random variable, and Kaspi's two-way rate-distortion [7] which uses multiple auxiliary random variables (collectively denoted by U here for convenience).

For Shannon's lossy source coding problem, the fact that the usual information-theoretic expression of the form (1) will continue to coincide with the operational rate-distortion function when the distortion function d can take the value ∞ , was established by Pinkston in [8] (also see [9][Ch 9,Historical Notes and References]). That this should also be true for the Gray-Leiner conditional rate-distortion function should be expected, since the encoder can group source samples that have the same side-information value into multiple conditional sources and code them separately like in Shannon's lossy source coding problem. We have, however, been unable to locate a reference which discusses the extension of the Gray-Leiner conditional rate-distortion function to non-finite distortion functions.

An example of an infinite-valued distortion function is the so-called erasure distortion function [10][Chapter 13, Problem 7]. Here, the distortion function equals zero if the source and reconstruction symbols agree, it

equals one if the reconstruction is a special erasure symbol (irrespective of the source symbol), and it equals infinity otherwise. The erasure distortion function has appeared in a number of recent works: in [11][Sec. III-B] in the context of the CEO problem, in [13] in the context of the multiple descriptions coding problem, and in [12] in the context of some examples.

The erasure distortion function was used in [14] and [15][Secs. VII.B,C] to construct the first example which shows that the usual information-theoretic two-way rate-distortion function with two messages can be strictly smaller than the usual one-message Wyner-Ziv information-theoretic rate-distortion function. Specifically, it was shown that for sufficiently correlated doubly-symmetric binary sources [5][Sec.II,Eqn.(20)] and the erasure distortion function, it is possible to make the ratio of the one-message rate to the two-message sum-rate arbitrarily large while simultaneously making the ratio of the backward rate to the forward rate in the two-message sum-rate arbitrarily small.

However, in [16] it was shown that, contrary to the claim made in footnote 8 of [15][Appendix E], the usual information-theoretic rate-distortion functions for the Wyner-Ziv and two-way source coding problems, which are used in [14] and [15][Secs. VII.B,C], are in fact strictly smaller than their operational counterparts for the erasure distortion function and doubly-symmetric binary sources which satisfy a positivity condition. This implies that the usual information-theoretic rate-distortion expressions are not always valid for non-finite distortion functions, such as the erasure distortion function, in a distributed source coding setting even though they are for the Shannon and Gray-Leiner settings.

For the Wyner-Ziv and two-way source coding problems, the correct information-theoretic rate-distortion expressions that coincide with their operational counterparts were characterized in [16] for the erasure distortion function when the source and side information satisfy a positivity condition. It was shown, in particular, that the two-message two-way rate-distortion function for this problem exactly coincides with the one-message Wyner-Ziv rate-distortion function casting into doubt some of the conclusions reached in [14] and [15][Secs. VII.B,C] about the benefit of interaction (two-way coding) for lossy source reproduction.

The aim of this article is to demonstrate that all the conclusions reached in [14] and [15][Secs. VII.B,C] regarding the benefit of interaction for lossy source reproduction are correct if one amends footnote 8 of [15][Appendix E] to the following: 'Although the usual information-theoretic expressions for the WynerZiv and two-message rate-distortion functions are not operationally attainable for the non-finite erasure distortion function, they can be operationally approached, as closely as desired, by replacing the ∞ value in the erasure distortion function with a sufficiently large, but finite, positive real number.'. The possibility of making such an approach work is contained in a suggestion of anonymous reviewers reported in Remark 1.1 of [16]. Such an approach is also taken in [11][Sec. III-B] where it is remarked that an infinite-valued distortion measure is unforgiving of decoding errors that have nonzero probability even if they are negligible. Decoding errors with vanishingly small but nonzero probability are unavoidable in most distributed source coding problems.

Formally speaking, our main result is that if d n is a sequence of bounded distortion functions that monotonically increases to a non-finite distortion function d ∞ , then the corresponding sequence of informationtheoretic rate-distortion functions R n ( D ) associated with d n also monotonically increases to the usual information-theoretic rate-distortion function R ∞ ( D ) associated with d ∞ . Thus by making d n approach d ∞ , it is possible to make the operational sum-rate-distortion function for d n as close as desired to R ∞ ( D ) even though R ∞ ( D ) itself may not be operationally attainable for d ∞ . Thus intuitions, examples, and broad qualitative conclusions (such as the benefit of interaction) that can be formed on the basis of examining R ∞ ( D ) will be essentially correct even from an operational perspective in the sense that while they may not hold true operationally for d ∞ itself, they will for d n , for all sufficiently large n .

## II. MAIN RESULT

We first establish a general result and then discuss its application to rate-distortion functions.

## A. Statement

Consider the following finite-dimensional optimization problem

$$\psi _ { n } ( D ) \coloneqq \inf _ { p \in \mathcal { C } _ { n } ( D ) \coloneqq \mathcal { A } \cap \mathcal { B } _ { n } ( D ) } f ( p )$$

where p is a probability vector (prob.vec.) in R k , k a finite positive integer, f is a real-valued continuous function of p , A is a fixed, nonempty, compact subset of probability vectors in R k , D is a finite nonnegative real number, and B n ( D ) := { p prob.vec. : 〈 p, d n 〉 ≤ D } , where 〈· , ·〉 denotes the usual Euclideanspace inner product and d n , n = 1 2 , , . . . , is a sequence of vectors in R k with finite, nonnegative, and nondecreasing components some of which (but not all) increase to ∞ while the rest monotonically increase to some finite nonnegative real numbers. We denote the limit of the d n 's by d ∞ . Thus, B n ( D ) (and therefore also C n ( D ) ) is a nested nonincreasing sequence of compact subsets of probability vectors in R k . We define

$$\mathcal { B } _ { \infty } ( D ) \coloneqq \cap _ { n = 1 } ^ { \infty } \mathcal { B } _ { n } ( D ) = \left \{ p \text{ prob.vec.} \, \colon \langle p, d _ { \infty } \rangle \leq D \right \}$$

with the convention 0 · ∞ := 0 in the inner product and observe that the support of any p in B ∞ ( D ) excludes components where d ∞ equals ∞ (since D is finite) and that B ∞ is a compact subset of probability vectors in R k . Let

$$\psi _ { \infty } ( D ) \coloneqq \inf _ { p \in \mathcal { C } _ { \infty } ( D ) \colon = \mathcal { A } \cap \mathcal { B } _ { \infty } ( D ) } f ( p ).$$

The main result is that ψ n ( D ) ↑ ψ ∞ ( D ) for all finite nonnegative D for which C ∞ ( D ) is nonempty.

## B. Proof

To prove this result, first observe that the constraint sets C n ( D ) are nested, nonincreasing, and contain C ∞ ( D ) . This implies that

$$\psi _ { n } ( D ) \uparrow \lim _ { n \to \infty } \psi _ { n } ( D ) \leq \psi _ { \infty } ( D ).$$

The main result is proved by establishing the reverse inequality. Towards this end, we note that the minimands in (2) and (4) are continuous functions and that the constraint sets are compact. This implies that there is a sequence of probability vectors p ( D ) n ∈ C n ( D ) and a probability vector p ( D ) ∞ ∈ C ∞ ( D ) such that ψ n ( D ) = f ( p ( D ) n ) for all n and ψ ∞ ( D ) = f ( p ( D ) ∞ ) . Since all these probability vectors belong to the compact set C 1 ( D ) , there is a subsequence p ( D ) n j , j = 1 2 , , . . . , converging to a probability vector q ( D ) ∞ . We will shortly show that q ( D ) ∞ is in C ∞ ( D ) . Then, since f is continuous, ψ n j ( d ) = f ( p ( d ) n j ) ↑ f ( q ( D ) ∞ ) ≥ ψ ∞ ( D ) establishing the reverse inequality. Finally, to see why q ( D ) ∞ is in C ∞ ( D ) , note that for components of d n j that increase to infinity, the corresponding components of p ( D ) n j must converge to zero since D is finite. The remaining components of both d n j and p ( D ) n j converge to finite nonnegative real numbers. Thus the subsequence of their (nonnegative) inner products, which is no more than D (finite and nonnegative), converges to the inner product of their limits.

## C. Application to rate-distortion functions

Comparing (1) and (2) it is apparent that they have the same form. To demonstrate that (1) is, in fact, a special case of (2), we only need to verify that the minimand, minimizing variable, and constraint sets of (1) satisfy all the assumptions in Sec. II-A that (2) is required to satisfy. First note that p in (2) corresponds to p X,Y,U,X ̂ in (1) with k = | X × Y × U × X | ̂ . The f in (1) is a finite linear combination of conditional mutual informations involving some or all the variables X,Y,U,X ̂ . This is a real-valued (in fact also nonnegative and bounded) continuous function of p X,Y,U,X ̂ since conditional mutual informations are continuous functions of the joint pmf of all the variables that appear in them [17][Ch.2,Sec.2.3], marginal pmfs are linear (therefore continuous) functions of the joint pmf, and the composition of a finite number of continuous functions is continuous. The set A in (1) is the set of pmfs which satisfy certain marginal

consistency and Markov-chain constraints associated with the random variables X,Y,U,X ̂ . This set is compact because a marginal consistency constraint is a linear equality constraint on the joint pmf (hence it defines a closed hyperplane within the bounded simplex of joint pmfs) and Markov-chain constraints can be expressed as the zero-level sets of appropriate conditional mutual information functions which, as we just discussed, are continuous. The distortion function d in (1) corresponds to d n in (2). They are both real-valued (finite) and nonnegative. The set B n ( D ) in (2) then corresponds to the expected distortion constraint E d X, Y, U, X [ ( ̂ )] ≤ D in (1). Finally, in typical scenarios of the source coding problems discussed in Sec. I (including [14] and [15][Secs. VII.B,C]), D is finite, and the feasible set { p X,Y,U,X ̂ ∈ A s.t. E d X, Y, U, X [ ( ̂ )] ≤ D } is nonempty.

## ACKNOWLEDGMENT

The author would like to thank Nan Ma, Sandeep Pradhan, and Bob Gray for their feedback on an earlier draft of this work.

## REFERENCES

- [1] C. E. Shannon, 'Coding theorems for a discrete source with a fidelity criterion,' in IRE Conv. Rec. , vol. 7, 1959, pp. 142-163. (Also in Information and Decision Processes , R. E. Machol, Ed. New York: McGraw-Hill, 1960, pp. 93-126, and in Claude Elwood Shannon: Collected Papers , N. J. A. Sloane and A. D. Wyner, Eds. Piscataway, NJ: IEEE Press, 1993, pp. 325-350.)
- [2] R. M. Gray, 'A new class of lower bounds to information rates of stationary sources via conditional rate-distortion functions,' IEEE Trans. Inform. Theory , vol. IT-19, pp. 480-489, 1973.
- [3] B. M. Leiner, 'Rate-distortion theory for sources with side information,' Ph. D. dissertation, Stanford Univ., Stanford, CA, Aug. 1973.
- [4] B. M. Leiner and R. M. Gray, 'Rate-distortion for ergodic sources with side information,' IEEE Trans. Inform. Theory , vol. IT-20, pp. 672-675, 1974.
- [5] A. D. Wyner and J. Ziv, 'The rate-distortion function for source coding with side-information at the receiver,' IEEE Trans. Inform. Theory , vol. IT-22, pp. 1-11, 1976.
- [6] A. D. Wyner, 'The rate-distortion function for source coding with side information at the decoder-II: General sources,' Inform. Contr. , vol. 38, pp. 60-80, 1978.
- [7] A. H. Kaspi, 'Two-way source coding with a fidelity criterion,' IEEE Trans. Inform. Theory , vol. IT-31, pp. 735-740, 1985.
- [8] J. T. Pinkston, 'Encoding Independent Sample Information Sources,' MIT Research Laboratory of Electronics, Technical Report 462 , 1967.
- [9] R. G. Gallager, Information Theory and Reliable Communication , John Wiles &amp; Sons, Inc., 1968.
- [10] T. M. Cover and J. A. Thomas, Elements of Information Theory , John Wiles &amp; Sons, Inc., 1991.
- [11] A. B. Wagner and V. Ananthram, 'An improved outer bound for multi terminal source coding,' IEEE Trans. Inform. Theory , vol. IT-54, pp. 1919-1937, 2008.
- [12] T. Courtade and R. Wesel, 'Multiterminal source coding with an entropy-based distortion measure,' in Proc. IEEE International Symposium on Information Theory , St. Petersburg, Russia, Jul. 31 - Aug. 5, 2011, pp. 2040-2044.
- [13] E. Ahmed and A. B. Wagner, 'Erasure multiple descriptions,' IEEE Trans. Inform. Theory , vol. IT-58, pp. 1328-1344, 2012.
- [14] N. Ma and P. Ishwar, 'Interaction strictly improves the Wyner-Ziv rate-distortion function,' in Proc. IEEE International Symposium on Information Theory , Austin, TX, USA, Jun. 13-18, 2010, pp. 61-65.
- [15] N. Ma and P. Ishwar, 'The infinite-message limit of two-terminal interactive source coding,' IEEE Trans. Inform. Theory , vol. IT-59, pp. 4071-4094, 2013.
- [16] Y-K. Chia and H. F. Chong, 'On lossy source coding with side information under the erasure distortion measure,' in Proc. IEEE International Symposium on Information Theory , Honolulu, HI, USA Jun. 29 - Jul. 4, 2014.
- [17] R. W. Yeung, A First Course in Information Theory , Kluwer Academic / Plenum Publishers, New York, NY, 2002.
Prakash Ishwar
A note on the sum-rate-distortion function of some lossy source coding problems involving infinite-valued distortion functions
For a number of lossy source coding problems it is shown that even if the
usual single-letter sum-rate-distortion expressions may become invalid for
non-infinite distortion functions, they can be approached, to any desired
accuracy, via the usual valid expressions for appropriately truncated finite
versions of the distortion functions.