Papoulis-Marks-Cheung Approach
The Papoulis-Marks-Cheung approach[1] is a theorem in multidimensional Shannon sampling theory that shows that the sampling density of a two-dimensional bandlimited function can be reduced to the support of the Fourier transform of the function. Applying a multidimensional generalization of a theorem by Athanasios Papoulis,[2] the approach was first proposed by Robert J. Marks II and Kwang Fai Cheung.[3][4] The approach has been called "elegant,"[5] "remarkably" closed,[6] and "interesting."[7]
The Theorem
[edit]The two-dimensional Fourier transform, or frequency spectrum, of a function is[8] where and are the spatial frequencies corresponding to and . When and are lengths, spatial frequency has units of cycles per unit length.
Prelee and Neuhoff[1] describe the Papoulis-Marks-Cheung approach as follows.
"Marks and Cheung focused on images with a given spectral support region and an initial base sampling lattice such that the induced spectral replicas of this support region do not overlap. They then showed that cosets of some sublattice could be removed from the base lattice until the sampling density was minimal … or approached minimal ... [This] allows the sampling rate to be reduced until it equals or approaches … [a] minimum." In this context, the limit is the area of the support of the spectrum."
In deriving their result, Marks and Cheung relied on Papoulis' generalized sampling expansion.[2][9]
Explanation
[edit]![](http://upload.wikimedia.org/wikipedia/commons/thumb/9/98/Papoulis-Marks-Cheung_Approach_Figure_1.jpg/220px-Papoulis-Marks-Cheung_Approach_Figure_1.jpg)
The Papoulis-Marks-Cheung approach is best explained by example. Consider from Figure 1 the half circle shown on the right half plane. A signal's spectrum, , is zero outside the half circle. Inside the circle, the spectrum' is arbitrary but is well behaved.[9] The half-circle, with unit radius, has an area of (cycles per unit length) squared.
According to the Papoulis-Marks-Cheung approach, the sampling density for the image can be reduced to samples per unit area. The Papoulis-Marks-Cheung approach informs how to do this.
To the right in Figure 1 is pictured a rectangular replication of the half circle which occurs when the two-dimensional function is sampled at spatial locations shown in Figure 2.
![](http://upload.wikimedia.org/wikipedia/commons/thumb/1/1c/Papoulis-Marks-Cheung_Approach_Figure_2.jpg/220px-Papoulis-Marks-Cheung_Approach_Figure_2.jpg)
This replication is a consequence of the multidimensional sampling theorem that shows that the sampling of a two-dimensional signal in the spatial domain results in spectrum replication in the Fourier domain. If the uniform sampling density were lower, the replications would overlap and an attempt at reconstruction of the original function would result in image aliasing. The sampling density to achieve this is equal to the area of the rectangular lattice cell of the spectrum replication. The corresponding area of the rectangle used in the replication is equal to (cycles per unit length) squared. As confirmed by Figure 2, the sampling density required to achieve the spectral replication is therefore samples per unit area. The Papoulis-Marks-Cheung approach says that this sampling density can be reduced to the area of the half circle, namely from to samples per unit area.
![](http://upload.wikimedia.org/wikipedia/commons/thumb/b/bb/Papoulis-Marks-Cheung_Approach_Figure_3.jpg/220px-Papoulis-Marks-Cheung_Approach_Figure_3.jpg)
To see how this reduction happens, consider Figure 3 where the rectangular lattice cell is divided into identical squares. Note that two of these squares lie totally in an area where the spectral replication is identically zero. These squares are shaded light green. Think of each of squares as spectra of different two-dimensional signals. All of the samples for the signals corresponding to the light green areas are zero and do not have to be considered. The area of the two green squares is . Since the samples corresponding to these squares do not have to be considered (they are all zero), the overall sampling density is reduced from samples per unit area to samples per unit area.
![](http://upload.wikimedia.org/wikipedia/commons/thumb/8/89/Papoulis-Marks-Cheung_Approach_Figure_4.jpg/220px-Papoulis-Marks-Cheung_Approach_Figure_4.jpg)
The corresponding reduction in sampling density is shown in Figure 4 where the red dots are locations where samples need not be taken. A single cell containing one red dot is shown shaded. The area of the cell is The corresponding reduction in sampling density is shown in Figure 4 where the red dots are locations where samples need not be taken. A single cell containing one red dot is shown shaded. The area of the cell is units. The sampling density is therefore, as also seen from the areas of two green squares in Figure 3, reduced by samples per unit area.
Extension
[edit]In the previous example, the squares in Figure 3 can be made arbitrarily small and increased in number so that, asymptotically, all of the area equal to zero can be covered. Thus, the sampling density can be reduced to the support of the spectrum, i.e., to the area where the spectrum is not identically zero.
The Papoulis-Marks-Cheung approach can straightforwardly be generalized to higher dimensions. Also, replication geometry need not be rectangular but can be any shape that will tile the entire plane such as parallelograms and hexagons.[10]
A more detailed mathematical description of the Papoulis-Marks-Cheung approach is available in the original paper by Marks and Cheung[4] and their derivative work.[11][9]
References
[edit]- ^ a b Prelee, Matthew A.; Neuhoff, David L. (May 2016). "Multidimensional Manhattan Sampling and Reconstruction". IEEE Transactions on Information Theory. 62 (5): 2772–2787. arXiv:1502.01694. doi:10.1109/TIT.2016.2542081. ISSN 0018-9448. S2CID 3199882.
- ^ a b Papoulis, A. (November 1977). "Generalized sampling expansion". IEEE Transactions on Circuits and Systems. 24 (11): 652–654. doi:10.1109/TCS.1977.1084284. ISSN 0098-4094.
- ^ Marks, Robert J. (1986-02-01). "Multidimensional-signal sample dependency at Nyquist densities". JOSA A. 3 (2): 268–273. Bibcode:1986JOSAA...3..268M. doi:10.1364/JOSAA.3.000268. ISSN 1520-8532.
- ^ a b Cheung, Kwan F.; Marks, Robert J. (1990-01-01). "Imaging sampling below the Nyquist density without aliasing". JOSA A. 7 (1): 92–105. Bibcode:1990JOSAA...7...92C. doi:10.1364/JOSAA.7.000092. ISSN 1520-8532.
- ^ Gardos, T.R.; Mersereau, R.M. (1991). "FIR filtering of images on a lattice with periodically deleted samples". [Proceedings] ICASSP 91: 1991 International Conference on Acoustics, Speech, and Signal Processing. pp. 2873–2876 vol.4. doi:10.1109/ICASSP.1991.151002. ISBN 0-7803-0003-3. S2CID 121408946.
- ^ Bresler, Y.; Feng, Ping (1996). "Spectrum-blind minimum-rate sampling and reconstruction of 2-D multiband signals". Proceedings of 3rd IEEE International Conference on Image Processing. Vol. 1. pp. 701–704 vol.1. doi:10.1109/ICIP.1996.559595. ISBN 0-7803-3259-8. S2CID 44441155.
- ^ Herley, Cormac; Wong, Ping Wah; Member, Senior (1999). "Minimum rate sampling and reconstruction of signals with arbitrary frequency support". IEEE Trans. Inform. Theory. 45 (5): 1555–1564. CiteSeerX 10.1.1.83.8638. doi:10.1109/18.771158.
- ^ Cohen, Leon (1995). Time-frequency analysis. Englewood Cliffs, N.J: Prentice Hall PTR. ISBN 0-13-594532-1. OCLC 31516509.
- ^ a b c Marks, Robert J. II (2009). Handbook of Fourier analysis & its applications. Oxford: Oxford University Press. ISBN 978-0-19-533592-7. OCLC 227191901.
- ^ Dudgeon, Dan E. (1984). Multidimensional digital signal processing. Russell M. Mersereau. Englewood Cliffs, NJ: Prentice-Hall. ISBN 0-13-604959-1. OCLC 9282699.
- ^ Cheung, Kwang F. (Dec 6, 2012). "A Multidimensional Extension of Papoulis' Generalized Sampling Expansion with Application in Minimum Sampling Density". In Marks, Robert J. II (ed.). Advanced Topics in Shannon Sampling and Interpolation Theory. Springer Science & Business Media. ISBN 978-1-4613-9757-1.