Computation and analysis of the N-Layer scalable rate-distortion function

Methods for determining and computing the rate-distortion (RD) bound for N-layer scalable source coding of a finite memoryless source are considered. Optimality conditions were previously derived for two layers in terms of the reproduction distributions q/sub y1/ and q/sub y2/|y/sub 1/. However, the ignored and seemingly insignificant boundary cases, where q/sub y1/=0 and q/sub y2/|y/sub 1/ is undefined, have major implications on the solution and its practical application. We demonstrate that, once the gap is filled and the result is extended to N-layers, it is, in general, impractical to validate a tentative solution, as one has to verify the conditions for all conceivable q/sub yi+1/,...,y/sub N/|y/sub 1/,...,y/sub i/ at each (y/sub 1/,...,y/sub i/) such that q/sub y1/,...,y/sub i/=0. As an alternative computational approach, we propose an iterative algorithm that converges to the optimal joint reproduction distribution q/sub y1/,...,y/sub N/, if initialized with q/sub y1/,...,y/sub N/>0 everywhere. For nonscalable coding (N=1), the algorithm specializes to the BlahutArimoto (1972) algorithm. The algorithm may be used to directly compute the RD bound, or as an optimality testing procedure by applying it to a perturbed tentative solution q. We address two additional difficulties due to the higher dimensionality of the RD surface in the scalable (N>1) case, namely, identifying the sufficient set of Lagrangian parameters to span the entire RD bound; and the problem of efficient navigation on the RD surface to compute a particular RD point.