Depth Analysis of Greyscale Integral Images Using Continuous Multiview Wavelet Transform.
There are many 3D technologies [1-3] in 3D imaging. Particularly, the wavelets are used in 3D imaging, for instance, in multiview video compression , in image coding , in image fusion , as a quality metric , etc. The disparity of stereoscopic images can be estimated , the shape of the photographed/synthesized objects can be analysed, and the depth can be extracted using the wavelet analysis of the integral images [9-11]. In this paper, we propose the technique to eliminate the texture effect. Based on the similarity between the 3D images, where a 3D content is represented in a single image plane consisting of the logical image cells , our results can be applied to integral [13-15], multiview [16-18], plenoptic images, and light-field displays [19-21].
The result of the wavelet analysis of 3D images, however, depends on the colour/texture of the surface of the object. This effect was noticed before, and the most images in [9-11] were binary BW images and most results were presented in a qualitative visual form. This unexpected undesirable effect must be reduced, but a solution has been unknown yet.
The intensity in any point of the image in the image plane is proportional to the brightness of the point of the object in that point. The intensity of all separated parts of voxel pattern  is equal to the brightness of the corresponding point of the object. The voxel patterns precede the wavelets; thus, this important property is kept in wavelets. Therefore, the wavelet coefficients depend on the colour; and the result of the wavelet analysis of a multiview image is proportional to the brightness of voxels (or pieces of surface in a texture model).
Define the central view (CV) of an integral image as an image, where the central pixels of all image cells of the integral image are assembled in accordance with the location of the cells; i.e., the centre pixel of the left top cell goes to the left top corner of the CV, the centre pixel of the right top cell goes to the right top corner of the CV, etc.; such image can be seen by a hypothetic (nonexistent) camera located at the centre of the lens array.
The CV can be calculated by applying the known interlacing technique (see, e.g., ) along both dimensions of the integral image. The CV can be calculated for binary, grey-scale, or colour images; see Figure 3. In our paper, the CV is solely used as a simple graphical picture of all objects of a 3D-image, because in the original integral image, the objects may look unclear (nonsharp, or out-of-focus, or blurred, etc.); compare Figures 1(a) and 1(b).
What is important is that it is unnecessary to calculate the CV for the wavelet transform, and therefore not a one operation with the wavelets needs the CV. In this context, the central view is nothing but a descriptive illustration for the journal paper.
Later in this paper we will name an integral/plenoptic image by its CV (as "in Figure 1(b)"), but actually such a referencing will mean the integral image (i.e., "Figure 1(a)") which corresponds to the named CV.
Previously, the colour dependence remained behind the scene, and the depth was mixed with colours. In this paper, as soon as we know an answer, we tried to open a door, slightly. To restore the depth more correctly, we propose to use the normalized image. Firstly, we would like to illustrate the dependence on colours in the wavelet analysis of the integral images.
2. Dependence of Wavelet Coefficients on Colours
Consider colours of three digits in the 123-image. If all colours are identical (as in the B image in Figure 1), then the wavelet coefficients are the same for every digit. Note that it is a repainted colour image, originally provided by Prof. B. Lee.
The wavelet analysis of the binary BW image Figure 1 shows almost uniform shape (a flat-top curve) for the depth planes between -6 and +6, where the digits of this 3D image are presumably located in space; see Figure 2(a). The wavelet coefficients for each digit processed separately are shown in Figure 2(b), where all three maxima are close to each other.
Based on Figure 2(b), one may conclude that the digit 1 is located between the +1st and +2nd planes, digit 2 in 0th, digit 3 between -2nd and -3rd planes. Outside of this region ([absolute value of (depth)] > 3), the wavelet coefficients monotonously decay. Therefore later in the related Figures 4 and 5, we will show the wavelet coefficients within the depth region [-2, +2] only, where the expected result is a flat-top horizontal line.
Then, if the shades of grey of digits or their colours are not the same as in Figure 3, the expected results of the wavelet analysis would be different for every digit.
The results of the wavelet analysis (wavelet coefficients) of greyscale images, Figures 3(a) and 3(b), are shown in Figure 4(a). The wavelet coefficients of the colour digits, Figures 3(c) and 3(d), are shown in Figure 4(b).
There is an essential difference between the wavelet coefficients of digits of various colours. In all cases, instead of a flat top and decay (as in Figure 2(a) for the binary image), the wavelet coefficients of grey images may rise up, fall down, or have a maximum or minimum in the middle; see Figure 4. These graphs can be treated ambiguously: either the digits in the same plane have different colours, or the digits of the same colour are located at different distances. A variety of intermediate interpretations is also possible.
This confirms a strong dependence of the results on the colours of voxels (texture); it is an undesirable side effect in the depth (shape) extraction. Below, we describe how to eliminate it.
3. Materials and Methods
In our examples, the source images are integral/multiview/ plenoptic images with the square grid of the image cells (a cell is an area under a lenticular lens). The images were taken from different independent sources, either photographs or synthesized (computer-generated) images.
The multiview wavelets and the algorithm of the continuous multiview wavelet transform are exactly the same as presented in our previous paper about the multiview wavelets [10, 11].
In order to avoid the undesirable effect of texture and to restore the spatial structure without regard to the colours, we propose to process so-called normalized image . A normalized image is typically used in order to reduce a nonuniform illumination in a local neighbourhood. The normalized image can be built by the algorithm . The normalized image corresponding to the original image of digits Figure 1 is shown in Figure 5(a). The result of the wavelet analysis of this normalized image is shown in Figure 5(b).
Despite insignificant imperfectness (such as non-100% horizontal flat top), the graphs in Figures 2(a) and 5(b) are similar to each other more than any graphs in Figure 4. This means that the normalized image predominantly contains the information about the 3D structure (depth), rather than about colours (texture). For instance, the RMS difference of the wavelet coefficients between the processed RGB and BW images is 1.2 A.U. (normalized RMS error 60%); in grey and BW images it is 0.7 A.U. (48%); however, in BW and normalized images it is 2.5-4.2 times less (0.3 A.U, 19%). This is not a complete elimination of the undesirable effect, but its essential reduction.
Before processing, the colour images were transformed into the grey-scaled images, but it is not necessary. The full-colour images can be processed as well, by processing the R, G, B colour components separately.
The dimensions of the images are as follows: the image of digits is 800 x 800 pixels (CV and normalized image are 40 x 40 pixels), the rabbit 1350 x 1350 (CV and normalized 45 x 45), the books 4031 x 4031 pixels (CV and normalized 139 x 139), and the house 4575 x 4575 pixels (CV and normalized 75 x 90).
To illustrate the proposed technique, we applied it to plenoptic/ integral images from various sources. The normalization is followed by the wavelet transform. In the examples, we will compare the results (i.e., the wavelet coefficients for the source and normalized images) by the depth planes and along the rows (horizontal lines) in the array of the coefficients.
The original and normalized images of the house are shown in Figure 6 (recall, this is a CV).
Consider, for example, the depth plane 0 and the row 70 of this image. The boundary line between the path and the lawn is at the same time the line of the change of colour of the texture. Because of that, the depth of this line is concealed (hidden) in the original image; see the wavelet coefficients and the graph along the horizontal row in Figures 7(a) and 7(b).What is important is that in the normalized image, this line can be clearly seen as a separate pulse in Figures 7(c) and 7(d).
N.B. The first image of each pair in Figure 7 displays the modulus of the wavelet coefficients (the black colour means maximum, and the white colour means zero), while the second graph is the full profile along the selected row. The same layout will be used later in Figures 9, 11, and 12.
Also, note the difference in the average (mean) level on the body of the car between the columns 45 and 65 as indicated by the dashed line in Figures 7(b) and 7(d). The influence of the texture colour of the body of the car is clearly reduced in Figure 7(d) down to the average level of that row.
The source and normalized images of books are shown in Figure 8.
In this example, we consider the depth plane -1, the row 95. Figures 9(c) and 9(d) clearly show a recognized 3D edge (backbone of the book) as a separate pulse, but not a texture-induced effect (a step pulse) as in Figures 9(a) and 9(b).
The influence of the texture is reduced in this image too. The average level at the cover of the book is almost the same along the row in the processed normalized image (shown in Figures 9(b) and 9(d) by the dashed line).
N.B. The letters "EG" on the cover of the book in Figure 9(a) are not a restored 3D structure, rather a texture of the surface.
The source and normalized images of the rabbit are shown in Figure 10. In the 3D analysis, two planes will be considered.
(1) Plane 2, row 13. The eye of the rabbit becomes clearly recognized in the normalized image, while it is completely invisible in the unnormalized one; see Figure 11.
(2) Plane 5, row 17. The same is valid for the nose. The nose is clearly recognized in the processed normalized image, while it is hidden in the source (unnormalized) image; see Figure 12.
These two examples also demonstrate that the influence of the texture is reduced in the normalized images, so as some previously unrecognized features appear.
We processed the greyscale images. Before processing, the colour images (if any) were transformed into the grey-scaled images, but it is not necessary. The full-colour images can be processed as well, by processing the RGB colour components. Also, we used the BW normalized image; however in general, each colour component can be normalized individually, and a colourful normalized image could be obtained.
The numeration of the rows and columns in the array of the wavelet coefficients is slightly different from that of rows and columns in the CV, because of the different sizes. (The size of the CV is fixed, but the array size at each depth plane varies, because we did not make any assumptions about the behaviour of the image beyond the sides of the original image.) The difference is less than the number of the current depth plane, but this small difference is out of the scope of the current paper.
The 3D images used in this paper were obtained by the independent authors either from a plenoptic camera (Figures 6 and 8) or by means of the computer simulation of the integral imaging (Figures 1 and 10). In either case, a lens array with a square grid of microlenses or its computational equivalent was used. The high quality lens arrays are known for their very uniform structure (a small deviation of the lens pitch across the array).The proposed image processing procedure can be applied to the 3D images built on the hexagonal grid of the lenses (with properly redesigned wavelets).
Instead of a lens array, the 3D mages can be also obtained from a camera array, as described in  and in the references in ; on the other hand, the layout of cameras in the camera array might be less uniform than that of the lens array. Moreover, the 3D imaging with the sensors at random locations is demonstrated . Therefore, we hope that a wavelet processing can be in principle applied to 3D images from the camera arrays; however the wavelets have to be radically modified in this case.
In the normalization of the most images (only a few of many processed images are presented in this paper), we used the default values [[sigma].sub.1] = 18, [[sigma].sub.2] = 24; however in a few cases, these values were reduced to [[sigma].sub.1] = 4.5, [[sigma].sub.2] = 6.
The influence of colours (texture) is essentially reduced by means of the normalized image instead of the source full-colour or greyscale image. The technique is confirmed by processing integral/plenoptic images from independent sources. The numerical comparison became possible between the planes and within each plane (because the wavelet coefficients are normalized by definition, and because the integral image itself is used in the normalized form). Integral, multiview, plenoptic (light-field) images with the square grid of cells can be processed. Colour images can be processed as well. The proposed technique can be efficiently used in 3D imaging for the depth extraction and the shape reconstruction without regard to the colour/texture.
The plenoptic/integral image data used to support the findings of this study were obtained from two sources. Namely, two of four images we used in the examples are available for free on the website "http://www.tgeorgiev.net/" by Dr. T. Georgiev ; the other two images were kindly supplied by Prof. B. Lee upon our personal request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
This work was partially supported for one author (Prof. I. Palchikova) by the Russian Foundation for Basic Research and by the Ministry of Education, Science, and Innovative Policy of the Novosibirsk Region within the framework of the research [project No. 17-47-540269]. We greatly appreciate Prof. B. Lee for the images provided to us personally (the digits and the rabbit). Also, we are thankful to Dr. T. Georgiev for his wonderful and very useful website (research until 2014), particularly, for the images available at that site for free (the input and the house crop).
 J. Hong, Y. Kim, H. Choi et al., "Three-dimensional display technologies of recent interest: principles, status, and issues," Applied Optics, vol. 50, no. 34, pp. H87-H115, 2011.
 J.-Y. Son, B. Javidi, and K.-D. Kwack, "Methods for displaying three-dimensional images," Proceedings of the IEEE, vol. 94, no. 3, pp. 502-522, 2006.
 M. Martinez-Corral and B. Javidi, "Fundamentals of 3D imaging and displays: a tutorial on integral imaging, light-field, and plenoptic systems," Advances in Optics and Photonics, vol. 10, no. 3, p. 512, 2018.
 M. Flierl and B. Girod, "Multiview video compression," IEEE Signal Processing Magazine, vol. 24, no. 6, pp. 66-76, 2007.
 W. Yang, Y. Lu, F. Wu, J. Cai, K. N. Ngan, and S. Li, "4-D-wavelet-based multiview video coding," IEEE Transactions on Circuits and Systems for Video Technology, vol. 16, no. 11, pp. 1385-1395, 2006.
 J. L. Rubio-guivernau, V. Gurchenkov, M. A. Luengo-oroz et al., "Wavelet-based image fusion in multi-view three-dimensional microscopy," Bioinformatics, vol. 28, no. 2, pp. 238-245, 2012.
 P. B. Zadehand C. V. Serdean, "Stereo video disparity estimation using multi-wavelets," in Proceedings of the Seventh International Conference on Digital Telecommunications (ICDT), pp. 50-54, 2012.
 E. Bosc, F. Battisti, M. Carli, and P. Le Callet, "A wavelet-based image quality metric for the assessment of 3D synthesized views," in Proceedings of the 24th IS and T/SPIE Stereoscopic Displays and Applications Conference, SD and A 2013, vol. 8648, USA, February 2013.
 V. Saveljev, "Wavelets and continuous wavelet transform for autostereoscopic multiview images," Journal of Electrical Engineering, vol. 4, no. 1, pp. 19-23, 2016.
 V. Saveljev and I. Palchikova, "Analysis of autostereoscopic three-dimensional images using multiview wavelets," Applied Optics, vol. 55, no. 23, pp. 6275-6284, 2016.
 V. Saveljev and I. Palchikova, "Analytical model of multiview autostereoscopic 3D display with a barrier or a lenticular plate," Journal of Information Display, vol. 19, no. 2, pp. 99-110, 2018.
 V. V. Saveljev and S.-J. Shin, "Layouts and cells in integral photography and point light source model," Journal of the Optical Society of Korea, vol. 13, no. 1, pp. 131-138, 2009.
 S.-G. Park, J. Yeom, Y. Jeong, N. Chen, J.-Y. Hong, and B. Lee, "Recent issues on integral imaging and its applications," Journal of Information Display, vol. 15, no. 1, pp. 37-46, 2014.
 X. Xiao, B. Javidi, M. Martinez-Corral, and A. Stern, "Advances in three-dimensional integral imaging: Sensing, display, and applications," Applied Optics, vol. 52, no. 4, pp. 546-560, 2013.
 J.-H. Park, K. Hong, and B. Lee, "Recent progress in three-dimensional information processing based on integral imaging," Applied Optics, vol. 48, no. 34, pp. H77-H94, 2009.
 N. A. Dodgson, "Analysis of the viewing zone of multi-view autostereoscopic displays," in Proceedings of the Stereoscopic Displays and Virtual Reality Systems IX, vol. 4660, pp. 254-265, USA, January 2002.
 B.-R. Lee, J.-J. Hwang, and J.-Y. Son, "Characteristics of composite images in multiview imaging and integral photography," Applied Optics, vol. 51, no. 21, pp. 5236-5243, 2012.
 Y. Takaki, K. Tanaka, and J. Nakamura, "Super multi-view display with a lower resolution flat-panel display," Optics Express, vol. 19, no. 5, pp. 4129-4139, 2011.
 R. Ng, M. Levoy, M. Bredif, G. Duval, M. Horowitz, and P. Hanrahan, "Light field photography with a hand-held plenoptic camera," Stanford [Tech Report], CTSR 2005-02, April 2005.
 T. Georgiev and A. Lumsdaine, "Focused plenoptic camera and rendering," Journal of Electronic Imaging, vol. 19, no. 2, 2010.
 J. Kim, S. Moon, Y. Jeong, C. Jang, Y. Kim, and B. Lee, "Dual-dimensional microscopy: real-time in vivo three-dimensional observation method using high-resolution light-field microscopy and light-field display," Journal of Biomedical Optics, vol. 23, no. 6, pp. 1-11, 2018.
 M.-C. Park, S. Ju Park, V. V. Saveljev, and S. Hwan Kim, "Synthesizing 3-D images with voxels," in Three-Dimensional Imaging, Visualization, and Display, pp. 207-225, Springer, 2009.
 P. V. Johnson, J. Kim, and M. S. Banks, "Stereoscopic 3D-display technique using spatiotemporal interlacing has improved spatial and temporal properties," Optics Express, vol. 23, no. 7, pp. 9252-9275, 2015.
 D. Sage and M. Unser, "Teaching image-processing programming in Java," IEEE Signal Processing Magazine, vol. 20, no. 6, pp. 43-52, 2003.
 D. Sage, "Local Normalization," Web-Page, 07 Sep. 2018, http://bigwww.epfl.ch/sage/soft/localnormalization/.
 T. Georgiev, "Research in the area of light fields until 2014," Web-Page, 13 Sep. 2018, http://tgeorgiev.net/.
 Y. Xing, Z.-L. Xiong, M. Zhao, and Q.-H. Wang, "Real-time integral imaging pickup system using camera array," in Proceedings of the Advances in Display Technologies VIII, vol. 10556, p. 12, San Francisco, USA, January 2018.
 M. Danesh Panah, B. Javidi, and E. A. Watson, "Three dimensional imaging with randomly distributed sensors," Optics Express, vol. 16, no. 9, pp. 6368-6377, 2008.
Vladimir Saveljev (1) and Irina Palchikova (2,3)
(1) NEMO Lab., Department of Physics, Myongji University, Yongin-si, Gyeonggi-do, 17058, Republic of Korea
(2) Technological Design Institute of Scientific Instrument Engineering, Siberian Branch of the Russian Academy of Sciences (TDISIE SB RAS), ul. Russkaya, 41, Novosibirsk, 630058, Russia
(3) Department of Physics, Novosibirsk State University (NSU), ul. Pirogova, 2, Novosibirsk, 630090, Russia
Correspondence should be addressed to Vladimir Saveljev; firstname.lastname@example.org
Received 10 October 2018; Accepted 28 October 2018; Published 2 December 2018
Guest Editor: Xiaowei Li
Caption: Figure 1: (a) Binary BW integral image of digits 800 x 800 pixels (credits for the original colour image to Prof. B. Lee, SNU); (b) the central view (CV) of this image (as seen through the camera located behind the central microlens of the lens array) 40 x 40 pixels.
Caption: Figure 2: Wavelet coefficients of binary BW image Figure 1: (a) whole image, (b) each digit separately. (The depth plane is an integer number.)
Caption: Figure 3: Various grey levels and colours of digits (CVs).
Caption: Figure 4: Wavelet coefficients of images from Figure 2: (a) grey-scale, (b) colour.
Caption: Figure 5: (a) Normalized image (CV). (b) Wavelet coefficients of normalized image.
Caption: Figure 6: Source and normalized images of the house (CVs). The original colour image "house crop" is available "for free" at the website by T. Georgiev .
Caption: Figure 7: Wavelet coefficients of source image (picture and graph on the top) and normalized image (on the bottom): modulus of the wavelet coefficients in (a), (c); the wavelet coefficients along the row 70 in (b), (d).
Caption: Figure 8: Source and normalized images of books (CVs). The original colour image "input" is available "for free" at the website of T. Georgiev .
Caption: Figure 9: Wavelet coefficients (modules on the left, graph along the row 95 on the right) of source (top) and normalized (bottom) images. The layout is the same as in Figure 7.
Caption: Figure 10: Source and normalized images of rabbit (CVs). The original colour image was provided by Prof. Lee by a personal request.
Caption: Figure 11: Wavelet coefficients in 2nd plane for source and normalized images, row 13 on the right. The layout is the same as in Figure 6.
Caption: Figure 12: Wavelet coefficients in 5th plane for source and normalized images, row 17 on the right. The layout is the same as in Figure 6.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Saveljev, Vladimir; Palchikova, Irina|
|Publication:||International Journal of Optics|
|Date:||Jan 1, 2018|
|Previous Article:||Solvent Effect on the Solid-Surface Fluorescence of Pyrene on Cellulose Diacetate Matrices.|
|Next Article:||Continuous-Wave and Passively Q-Switched Er:[Y.sub.2][O.sub.3] Ceramic Laser at 2.7 [micro]m.|