|
|
Compression of digital hologram for three-dimensional object using Wavelet-Bandelets transform |
Optics Express, Vol. 19, Issue 9, pp. 8019-8031 (2011)
http://dx.doi.org/10.1364/OE.19.008019
Acrobat PDF (1578 KB)
Abstract
In the transformation based compression algorithms of digital hologram for three-dimensional object, the balance between compression ratio and normalized root mean square (NRMS) error is always the core of algorithm development. The Wavelet transform method is efficient to achieve high compression ratio but NRMS error is also high. In order to solve this issue, we propose a hologram compression method using Wavelet-Bandelets transform. Our simulation and experimental results show that the Wavelet-Bandelets method has a higher compression ratio than Wavelet methods and all the other methods investigated in this paper, while it still maintains low NRMS error.
© 2011 OSA
1. Introduction
I. Yamaguchi and T. Zhang, “Phase-shifting digital holography,” Opt. Lett. 22(16), 1268–1270 (1997). [CrossRef] [PubMed]
M. L. Piao, N. Kim, and J. H. Park, “Phase contrast projection display using photopolymer,” J. Opt. Soc. Korea 12(4), 319–325 (2008). [CrossRef]
T. J. Naughton, Y. Frauel, B. Javidi, and E. Tajahuerce, “Compression of digital holograms for three-dimensional object reconstruction and recognition,” Appl. Opt. 41(20), 4124–4132 (2002). [CrossRef] [PubMed]
A. Shortt, T. J. Naughton, and B. Javidi, “Compression of digital holograms of three-dimensional objects using wavelets,” Opt. Express 14(7), 2625–2630 (2006). [CrossRef] [PubMed]
A. E. Shortt, T. J. Naughton, and B. Javidi, “Histogram approaches for lossy compression of digital holograms of three-dimensional objects,” IEEE Trans. Image Process. 16(6), 1548–1556 (2007). [CrossRef] [PubMed]
H. Yoshikawa and J. Tamai, “Holographic image compression by motion picture coding,” Proc. SPIE 2652, 2–9 (1996). [CrossRef]
Y. H. Seo, H. J. Choi, J. S. Yoo, G. S. Lee, C. H. Kim, S. H. Lee, S. H. Lee, and D. W. Kim, “Digital hologram compression technique by eliminating spatial correlations based on MCTF,” Opt. Commun. 283(21), 4261–4270 (2010). [CrossRef]
E. Darakis and T. J. Naughton, “Compression of digital hologram sequences using MPEG-4,” Proc. SPIE 7358, 735811 (2009). [CrossRef]
H. Yoshikawa and J. Tamai, “Holographic image compression by motion picture coding,” Proc. SPIE 2652, 2–9 (1996). [CrossRef]
Y. H. Seo, H. J. Choi, J. S. Yoo, G. S. Lee, C. H. Kim, S. H. Lee, S. H. Lee, and D. W. Kim, “Digital hologram compression technique by eliminating spatial correlations based on MCTF,” Opt. Commun. 283(21), 4261–4270 (2010). [CrossRef]
T. W. Ng and K. T. Ang, “Fourier-transform method of data compression and temporal fringe pattern analysis,” Appl. Opt. 44(33), 7043–7049 (2005). [CrossRef] [PubMed]
E. Le Pennec and S. Mallat, “Sparse geometric image representation with Bandelets,” IEEE Trans. Image Process. 14(4), 423–438 (2005). [CrossRef] [PubMed]
2. Wavelet transform and Wavelet-Bandelets transform
2.1. Wavelet transform
2.2. Bandelets transform
E. Le Pennec and S. Mallat, “Sparse geometric image representation with Bandelets,” IEEE Trans. Image Process. 14(4), 423–438 (2005). [CrossRef] [PubMed]
E. Le Pennec and S. Mallat, “Bandelets image approximation and compression,” Multiscale Model. Simul. 4(3), 992–1039 (2005). [CrossRef]
- • Implementation of Wavelet transform with data from a digital image.
- • A warped Wavelet transform with a sub-band filtering along the flow lines [14].
- • Computation of Bandelets coefficients from the warped Wavelet coefficients along the flow lines.
E. Le Pennec and S. Mallat, “Bandelets image approximation and compression,” Multiscale Model. Simul. 4(3), 992–1039 (2005). [CrossRef]
3. Proposed method
- First step: Digital hologram of a 3-D object is recoded using an optical system shown in Fig. 4 . Four holograms with phase shifts of 0, , π, and of reference beam are recorded to yield amplitude and phase of the object optical field based on the principle of phase-shifting digital holography [1,2
I. Yamaguchi and T. Zhang, “Phase-shifting digital holography,” Opt. Lett. 22(16), 1268–1270 (1997). [CrossRef] [PubMed]
].M. L. Piao, N. Kim, and J. H. Park, “Phase contrast projection display using photopolymer,” J. Opt. Soc. Korea 12(4), 319–325 (2008). [CrossRef]
- Second step: Wavelet transform is applied to each of amplitude and phase of the captured hologram with the choice of mother Wavelet function and level of transformation. In our experiment, the simplest Wavelet function, Haar, was applied.
- Third step: Hologram data after Wavelet transform is divided into sub-bands. For each sub-band, quad-tree decomposition algorithm is applied to divide the sub-band into small blocks. The size of the small block is given by × pixels. Block size is important because it affects the performance of the geometric flow detection remarkably. When the block size is small, geometric flow can be determined easily in each block due to small number of pixels, but the large number of blocks makes overall compression ratio low. When the block size is large, on the contrary, determination of the geometric flow direction in each block is not easy since the pixel to pixel variation of the hologram data is usually large. In our experiment, various block sizes from 4 × 4 pixels to 32 × 32 pixels are tested to find optimum size.
- Fourth step: The direction of the geometric flow in each block is determined. In order to find the best direction in a small block, all directions in the block are checked. The number of directions that has to be checked is 2(N-1) for a small block of N × N pixel size. The geometric flow direction is defined as a direction along which the data has minimum variation. The examination of the direction is implemented as follows. In order to examine the data variation along a test direction, the data in the block is first rearranged to 1-D series following the test direction. Haar Wavelet transform is then applied to this 1-D data series to yield 1-D Wavelet coefficient series. The resultant coefficients are compared with a threshold value T. The best direction is determined by a direction which gives minimum number of the coefficients larger than T. This procedure is illustrated in Fig. 5 and Figs. 6(a) and 6(b).
- Figure 6(a) describes the fringe of the sub-band after wavelet transform is performed with real information of a digital hologram. First step of this figure is to implement 2-D Wavelet transform with five levels transform. The examination of the fringe is implemented in each sub-band; this figure shows one small block which 8x8 pixels block size in fist sub-band with pixel value is represented in grey level. In the small block, it shows a part of the fringes in the sub-band. By using Bandelets transform, the fringe in each small block was determind. From this, we can approximate value pixel in each small block.
- Figure 6(b) describes the process which find the best direction in a small block, each small block will be checked all directions. In this figure, small block have 8x8 pixels size so it has 2*(8-1) = 14 directions are checked. The first step, the data in the small block is arranged to 1-D series following the test direction after that we applied to 1-D wavelet transform with this data (third picture of Fig. 6(b)). The resultant coefficients of wavelet transform implement absolute values and compare with a threshold value T (the red line that is described as the value T in third picture). The best direction is determined by a direction which gives minimum number of the coefficients larger than T. Cases Test (a), Test (b), Test (c) in this figure describe different cases which have different direction.
- Fifth step: The data in each block is approximated using the geometric flow direction to reduce the data amount. The original 2-D data is divided into a set of 1-D flow according to the geometric flow direction found in previous step as shown in Fig. 6(c). For the central flow, the minimum pixel value is selected as a reference and other pixels whose deviation from the reference is less than a parameter C are clipped to the reference value. For other flows, the minimum pixel value of the previous flow is selected as a reference and the same process is conducted. The result of this process is Wavelet-Bandelets coefficients and these coefficients are encoded as usual compression algorithms.
- Figure 6(c) show the example of pixel value approximation when C = 5. An small block has 4x4 pixel size with the best direction which was diagonal line from pixel position [1,1
I. Yamaguchi and T. Zhang, “Phase-shifting digital holography,” Opt. Lett. 22(16), 1268–1270 (1997). [CrossRef] [PubMed]
] to pixel position [4I. Yamaguchi and T. Zhang, “Phase-shifting digital holography,” Opt. Lett. 22(16), 1268–1270 (1997). [CrossRef] [PubMed]
,4G. A. Mills and I. Yamaguchi, “Effects of quantization in phase-shifting digital holography,” Appl. Opt. 44(7), 1216–1225 (2005). [CrossRef] [PubMed]
]. The directions is remarked (1), (2), (3), (4), (5), (6), (7), the main direction is remarked (1), respectively located above (2), (3), (4) and respectively located lower (5), (6), (7). First step, we approximate with (1) directions (it show in third picture), the minimum of value pixel in this direction is 12 so the value pixel in this direction is approximated equate 12 if the remaining pixel values in the range 12 ± 5, so we see, in this direction have three pixels which change value pixel (the red letter in third picture). Second step, we approximate with the remaining directions include (2), (3), (4), (5), (6) and (7); this process is described as follows: the minimum of value pixel in the (1) direction is approximated for (2) and (5) direction, the minimum of value pixel in the (2) direction is approximated for (3) direction; similarly, the minimum of value pixel in (5) direction is approximated for (6) direction and perform this process until the end direction in small block. In this figure, the minimum of value pixel in the (2) direction is 8 so the value pixel in (3) direction is approximated in the range 8 ± 5, similarly, the minimum of value pixel in the (5) direction is 10 so the value pixel in (6) direction is approximated in the range 10 ± 5. In this figure, the value pixel is approximated is red letter. Finally, we implement run length code with flow direction which we choose.G. A. Mills and I. Yamaguchi, “Effects of quantization in phase-shifting digital holography,” Appl. Opt. 44(7), 1216–1225 (2005). [CrossRef] [PubMed]
4. Experimental results
5. Conclusion
Acknowledgment
References and links
I. Yamaguchi and T. Zhang, “Phase-shifting digital holography,” Opt. Lett. 22(16), 1268–1270 (1997). [CrossRef] [PubMed] | |
M. L. Piao, N. Kim, and J. H. Park, “Phase contrast projection display using photopolymer,” J. Opt. Soc. Korea 12(4), 319–325 (2008). [CrossRef] | |
T. J. Naughton, Y. Frauel, B. Javidi, and E. Tajahuerce, “Compression of digital holograms for three-dimensional object reconstruction and recognition,” Appl. Opt. 41(20), 4124–4132 (2002). [CrossRef] [PubMed] | |
G. A. Mills and I. Yamaguchi, “Effects of quantization in phase-shifting digital holography,” Appl. Opt. 44(7), 1216–1225 (2005). [CrossRef] [PubMed] | |
A. Shortt, T. J. Naughton, and B. Javidi, “Compression of digital holograms of three-dimensional objects using wavelets,” Opt. Express 14(7), 2625–2630 (2006). [CrossRef] [PubMed] | |
A. E. Shortt, T. J. Naughton, and B. Javidi, “Histogram approaches for lossy compression of digital holograms of three-dimensional objects,” IEEE Trans. Image Process. 16(6), 1548–1556 (2007). [CrossRef] [PubMed] | |
H. Yoshikawa and J. Tamai, “Holographic image compression by motion picture coding,” Proc. SPIE 2652, 2–9 (1996). [CrossRef] | |
Y. H. Seo, H. J. Choi, J. S. Yoo, G. S. Lee, C. H. Kim, S. H. Lee, S. H. Lee, and D. W. Kim, “Digital hologram compression technique by eliminating spatial correlations based on MCTF,” Opt. Commun. 283(21), 4261–4270 (2010). [CrossRef] | |
E. Darakis and T. J. Naughton, “Compression of digital hologram sequences using MPEG-4,” Proc. SPIE 7358, 735811 (2009). [CrossRef] | |
Y. H. Seo, H. J. Choi, and D. W. Kim, “3D scanning-based compression technique for digital hologram video,” Signal Process. 22(2), 144–156 (2007). | |
E. Le Pennec and S. Mallat, “Sparse geometric image representation with Bandelets,” IEEE Trans. Image Process. 14(4), 423–438 (2005). [CrossRef] [PubMed] | |
E. Le Pennec and S. Mallat, “Non linear image approximation with Bandelets,” Tech. Rep. CMAP / École Polytechnique (2003). | |
E. Le Pennec and S. Mallat, “Bandelets image approximation and compression,” Multiscale Model. Simul. 4(3), 992–1039 (2005). [CrossRef] | |
C. Bernard and E. Le Pennec, “Adaptation of regular grid filtering to irregular grids,” Tech. Rep. CMAP / École Polytechnique (2003). | |
T. Bose, Digital Signal and Image Processing (Wiley, 2003), Chap. 11, pp. 623–669. | |
T. W. Ng and K. T. Ang, “Fourier-transform method of data compression and temporal fringe pattern analysis,” Appl. Opt. 44(33), 7043–7049 (2005). [CrossRef] [PubMed] |
OCIS Codes
(090.0090) Holography : Holography
(090.1760) Holography : Computer holography
(100.7410) Image processing : Wavelets
(090.1995) Holography : Digital holography
ToC Category:
Holography
History
Original Manuscript: December 7, 2010
Revised Manuscript: February 13, 2011
Manuscript Accepted: March 20, 2011
Published: April 12, 2011
Citation
Le Thanh Bang, Zulfiqar Ali, Pham Duc Quang, Jae-Hyeung Park, and Nam Kim, "Compression of digital hologram for three-dimensional object using Wavelet-Bandelets transform," Opt. Express 19, 8019-8031 (2011)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-19-9-8019
Sort: Year | Journal | Reset
References
- I. Yamaguchi and T. Zhang, “Phase-shifting digital holography,” Opt. Lett. 22(16), 1268–1270 (1997). [CrossRef] [PubMed]
- M. L. Piao, N. Kim, and J. H. Park, “Phase contrast projection display using photopolymer,” J. Opt. Soc. Korea 12(4), 319–325 (2008). [CrossRef]
- T. J. Naughton, Y. Frauel, B. Javidi, and E. Tajahuerce, “Compression of digital holograms for three-dimensional object reconstruction and recognition,” Appl. Opt. 41(20), 4124–4132 (2002). [CrossRef] [PubMed]
- G. A. Mills and I. Yamaguchi, “Effects of quantization in phase-shifting digital holography,” Appl. Opt. 44(7), 1216–1225 (2005). [CrossRef] [PubMed]
- A. Shortt, T. J. Naughton, and B. Javidi, “Compression of digital holograms of three-dimensional objects using wavelets,” Opt. Express 14(7), 2625–2630 (2006). [CrossRef] [PubMed]
- A. E. Shortt, T. J. Naughton, and B. Javidi, “Histogram approaches for lossy compression of digital holograms of three-dimensional objects,” IEEE Trans. Image Process. 16(6), 1548–1556 (2007). [CrossRef] [PubMed]
- H. Yoshikawa and J. Tamai, “Holographic image compression by motion picture coding,” Proc. SPIE 2652, 2–9 (1996). [CrossRef]
- Y. H. Seo, H. J. Choi, J. S. Yoo, G. S. Lee, C. H. Kim, S. H. Lee, S. H. Lee, and D. W. Kim, “Digital hologram compression technique by eliminating spatial correlations based on MCTF,” Opt. Commun. 283(21), 4261–4270 (2010). [CrossRef]
- E. Darakis and T. J. Naughton, “Compression of digital hologram sequences using MPEG-4,” Proc. SPIE 7358, 735811 (2009). [CrossRef]
- Y. H. Seo, H. J. Choi, and D. W. Kim, “3D scanning-based compression technique for digital hologram video,” Signal Process. 22(2), 144–156 (2007).
- E. Le Pennec and S. Mallat, “Sparse geometric image representation with Bandelets,” IEEE Trans. Image Process. 14(4), 423–438 (2005). [CrossRef] [PubMed]
- E. Le Pennec and S. Mallat, “Non linear image approximation with Bandelets,” Tech. Rep. CMAP / École Polytechnique (2003).
- E. Le Pennec and S. Mallat, “Bandelets image approximation and compression,” Multiscale Model. Simul. 4(3), 992–1039 (2005). [CrossRef]
- C. Bernard and E. Le Pennec, “Adaptation of regular grid filtering to irregular grids,” Tech. Rep. CMAP / École Polytechnique (2003).
- T. Bose, Digital Signal and Image Processing (Wiley, 2003), Chap. 11, pp. 623–669.
- T. W. Ng and K. T. Ang, “Fourier-transform method of data compression and temporal fringe pattern analysis,” Appl. Opt. 44(33), 7043–7049 (2005). [CrossRef] [PubMed]
Cited By |
OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.





OSA is a member of 