# Patent

## Title
Multi-perspective stereoscopy from light fields

## Abstract

Methods and systems for generating stereoscopic content with granular control over binocular disparity based on multi-perspective imaging from representations of light fields are provided. The stereoscopic content is computed as piecewise continuous cuts through a representation of a light field, minimizing an energy reflecting prescribed parameters such as depth budget, maximum binocular disparity gradient, desired stereoscopic baseline. The methods and systems may be used for efficient and flexible stereoscopic post-processing, such as reducing excessive binocular disparity while preserving perceived depth or retargeting of already captured scenes to various view settings. Moreover, such methods and systems are highly useful for content creation in the context of multi-view autostereoscopic displays and provide a novel conceptual approach to stereoscopic image processing and post-production.

## Background
BACKGROUND
The present invention relates to image generation, and in particular to methods and apparatus for generating stereoscopic image pairs based on multi-perspective imaging from a light field.
Three-dimensional (“3D”) television, movies, and games and displays for other purposes have been gaining more and more popularity both within the entertainment industry and among consumers. An ever increasing amount of content is being created, distribution channels including live-broadcast are being developed, and 3D monitors and TV sets are being sold in all major electronic stores.
One approach to 3D imagery is to generate a stereoscopic image pair where one image of the pair is provided to the viewer's left eye and the other image of the pair is provided to the viewer's right eye with something to eliminate or lessen cross-over. Where the images in the image pair are related (e.g., captured from the same scene) but different, those differences are interpreted by the viewer's brain as depth information, thus creating 3D effects. The image pair might be of a physical scene, where light from objects in the scene are captured by a camera or other optical sensor, and the stereoscopic effect is generated by capturing the scene imagery using two cameras offset by some baseline amount. The image pair might be of a virtual scene, such as a scene generated entirely using computer processing and/or geometric models.
Binocular parallax (i.e., binocular disparity) is one cue that stereoscopic image generation systems use for generating stereoscopic scene perception. In stereography, one common method for controlling the amount of binocular parallax is based on setting the baseline, or the inter-axial distance, of two cameras prior to image acquisition. However, the range of admissible baselines is quite limited, since most scenes exhibit more disparity than humans can tolerate when viewing the content on a stereoscopic display. Reducing the baseline of cameras decreases the amount of binocular disparity, but it also causes scene elements to appear overly flat.
Another, more sophisticated, approach to disparity control requires remapping image disparities (or remapping the depth of scene elements) and re-synthesizing new images. However, this approach typically requires accurate disparity computation and hole filling (filling in gaps that appear in the image because scene elements are moved in the re-synthesized views). For computer-generated images, depth remapping of scene elements implies severe changes of lighting, shading, and the scene composition in general.
The computer graphics and computer vision community has studied the geometry and applications of multi-perspective imaging, for example [Wood et al. 1997] describe a computer-assisted method to compute multi-perspective panoramas from a collection of perspective images, and employed multi-perspective imaging in movie production in order to provide a richer and more complete visualization of stereoscopic contents, for example for drawing backgrounds for two-dimensional (“2D”) cell animation [Thomas and Johnston 1995]. In the recent years, many types of multi-perspective cameras and corresponding images have been introduced. Examples include push-broom cameras [Hartley and Gupta 1997] and related Multiple-Center-of-Projection Images [Rademacher and Bishop 1998], cross slit cameras [Pajdla 2002; Zomet et al. 2003], or general linear cameras [Yu and McMillan 2004]. However, creating a richer and perceptually pleasing multi-perspective stereoscopic content without any specific camera model remains a difficult problem.
It has been shown that remapping binocular disparities may be used to refine and optimize stereoscopic content for display on different output devices or according to user preferences. A number of technical approaches for disparity remapping for stereoscopic content have been proposed. [Jones et al. 2001] analyze the scene depth range and adjust the stereoscopic camera baseline to a given disparity budget. Feldman et al. [2003] present a system that uses a nonlinear depth-scaling for transmitting a three-dimensional (3D) scene to be rendered from multiple views. [Holliman 2004] describes a system that compresses the scene depth for stereoscopic displays by identifying a region of interest and compressing it differently compared to the rest of the scene. Koppal et al. [2011] discusses optimal stereo and describes basic post-processing tools with their main focus on shot planning during capture. Ward et al. [2011] proposed a system for 2D-to-3D conversion that relies on image warping and requires manual interaction. [Kim et al. 2008] discuss how to perform non-linear depth remapping for multi-view autostereoscopic displays. [Zwicker et al. 2006] present a remapping and a pre-filtering framework for automultiscopic displays that adapts an in-put light field to the display capabilities.
All these works, however, are restricted in the type of disparity remapping operators they support. In particular, they do not provide a solution for detailed control of disparity in real world images. Although a nonlinear and local disparity remapping to control and retarget the depth of a stereoscopic content has been presented by [Lang et al. 2010], this method is limited in amount of remapping that it applied without producing noticeable distortions of the image content. In particular, this method causes salient scene structures to bend (e.g., straight lines) and fails to allow for per-pixel disparity control.
Accordingly, it is desirable to develop methods and systems to overcome the aforementioned deficiencies to have per-pixel disparity control by selecting actual light rays from an input light field instead of using image deformations or inpainting.
REFERENCES


    [Adelson and Bergen 1991] Adelson, E., and Bergen, J., The plenoptic function and the elements of early vision, Computational Models of Usual Processing, 1991
    [Adelson and Wang 1992] Adelson, E. H., and Wang, J., Single lens stereo with a plenoptic camera, IEEE Transactions on Pattern Analysis and Machine Intelligence, 14, 2, 99-106, 1992
    [Agrawal and Raskar 2007] Agrawal, A., and Raskar, R., Gradient domain manipulation techniques in vision and graphics, In ICCV Courses, 2007
    [Boykov and Kolmogorov 2004] Boykov, Y., and Kolmogorov, V., An experimental comparison of min-cut/max-flow algorithms for energy minimization in vision,
    IEEE Trans. Pattern Anal. Mach Intell., 26, 9, 1124-1137, 2004
    [Boykov et al. 2001] Boykov, Y., Veksler, 0., and Zabih, R., Fast approximate energy minimization via graph cuts, IEEE Trans. Pattern Anal. Mach. Intell., 23, 11, 1222-1239, 2001
    [Chai et al. 2000] Chai, J., Chan, S.-C., Shum, H.-Y., and Tong, X., Plenoptic sampling, In SIGGRAPH, 307-318, 2000
    [Feldmann et al. 2003] Feldmann, L, Schreer, 0., and Kauff, P., Nonlinear depth scaling for immersive video applications, In Proceedings of WIAMIS, 2003
    [Georgiev et al. 2006] Georgiev, T., Zheng, C., Nayar, S., Curless, B., Salesin, D.,
    and Intwala, C., Spatio-angular resolution trade-offs in integral photography, In Rendering Techniques: 17th Eurographics Workshop on Rendering, 2006
    [Gortler et al. 1996] Gortler, S. J., Grzeszczuk, R., Szeliski, R., and Cohen, M. F., The lumigraph, In SIGGRAPH, 43-54, 1996
    [Halle 1998] Halle, M., Multiple viewpoint rendering, In Proceedings of SIGGRAPH 98, Computer Graphics Proceedings, Annual Conference Series, 243 254, 1998
    [Hartley and Gupta 1997] Hartley, R. L, and Gupta, R., Linear pushbroom cameras, IEEE Transactions on Pattern Analysis and Machine Intelligence, 19, 9, 963-975, 1997
    [Holliman 2004] Holliman, N., Mapping perceived depth to regions of interest in stereoscopic images, In Stereoscopic Displays and Applications XV, Proceedings of SPIE, Vol. 5291, 117-128, 2004
    [Jones et al. 2001] Jones, G., Lee, D., Holliman, N., and Ezra, D., Controlling perceived depth in stereoscopic images, In Stereoscopic Displays and Virtual Systems VIII, Proceedings of SPIE, Vol. 4297, 42-53, 2001
    [Kang and Szeliski 2004] Kang, S. B., and Szeliski, R., Extracting view-dependent depth maps from a collection of images, IJCV 58, 2, 139-163, 2004
    [Kim et al. 2008] Kim, M., Lee, S., Choi, C., Urn, G.-M., Hur, N., and Kim, J., Depth scaling of multiview images for automultiscopic 3d monitors, in 3DTV Conference, 181-184, 2008
    [Kohli and Torr 2010] Kohli, P., and Torr, P. H. S., Dynamic graph cuts and their applications in computer vision, In Computer Vision: Detection, Recognition and Reconstruction, 51-108, 2010
    [Lang et al. 2010] Lang, M., Hornung, A., Wang, 0., Poulakos, S., Smolic, A., and Gross, M., Nonlinear disparity mapping for stereoscopic 3d, ACM Transactions on Graphics, 29, 4 (July), 75:1 75:10, 2010
    [Levoy and Hanrahan 1996] Levoy, M., and Hanrahan, P., Light field rendering, In SIGGRAPH, 31-42, 1996
    [Mendiburu 2009] Mendiburu, B., 3D Movie Making: Stereoscopic Digital Cinema from Script to Screen, Focal Press, 2009
    [Ng et al. 2005] Ng, R., Levoy, M., Bredif, M., Duval, G., Horowitz, M., and Hanrahan, P., Light field photography with a hand-held plenoptic camera, Stanford University Computer Science Technical Report CSTR, 2005-02, 2005
    [Pajdla 2002] Pajdla, T., Geometry of two-slit camera, Research Report CTU-CMP-2002-02, Czech Technical University, Prague, 2002
    [Peleg et al. 2001] Peleg, S., Ben-Ezra, M., and Pritch, Y., Omnistereo: Panoramic stereo imaging, IEEE Transactions on Pattern Analysis and Machine Intelligence, 23, 3, 279-290, 2001
    [Rademacher 1998] Rademacher, P., and Bishop, G., Multiple-center-of-projection images, In Proceedings of SIGGRAPH, 98, Computer Graphics Proceedings, Annual Conference Series, 199-206, 1998
    [Rubinstein 2008] Rubinstein, M., Shamir, A., and Avidan, S., Improved seam carving for video retargeting, ACM Trans. Graph, 27, 3, 2008
    [Seitz 2001] Seitz, S., The space of all stereo images, In IEEE International Conference on Computer Vision, vol. 1, 26-33, 2001
    [Shum and He 1999] Shum, H.-Y., and He, L.-W., Rendering with concentric mosaics, In Proceedings of SIGGRAPH, 99, Computer Graphics Proceedings, Annual Conference Series, 299-306, 1999
    [Thomas and Johnston 1995] Thomas, F., and Johnston, 0., Disney Animation: The Illusion Of Life, Hyperion, Los Angeles, 1995
    [Veeraraghavan 2007] Veeraraghavan, A., Raskar, R., Agrawal, A. K., Mohan, A., and Tumblin, J., Dappled photography: mask enhanced cameras for heterodyned light fields and coded aperture refocusing, ACM Trans. Graph., 26, 3. 2007
    [Ward et al. 2011] WARD, B., KANG, S. B., AND BENNETT, E. P., Depth director: A system for adding depth to movies. IEEE CG&A 31, 1, 36-48.
    [Wilburn et al. 2005] Wilburn, B., Joshi, N., Vaish, V., Talvala, E.-V., Antunez, E., Barth, A., Adams, A., Horowitz, M., and Levoy, M., High performance imaging using large camera arrays, ACM Transactions on Graphics, 24, 3 (August), 765-776, 2005
    [Wood 1997] Wood, D. N., Finkelstein, A., Hughes, J. F., Thayer, C. E., and Salesin, D. H., Multiperspective panoramas for cel animation, In Proceedings of SIGGRAPH, 97, Computer Graphics Proceedings, Annual Conference Series, 243-250, 1997
    [Yang et al. 2002] Yang, J. C., Everett, M., Buehler, C., and McMillan, L., A real-time distributed light field camera, In Rendering Techniques 2002: 13th Eurographics Workshop on Rendering, 77-86, 2002
    [Yu and McMillan 2004] Yu, J., and McMillan, L., General linear cameras, In 8th European Conference on Computer Vision (ECCV 2004), 14-27, 2004
    [Yu et al. 2009] Yu, X., Yu, J., and McMillan, L., Towards multi-perspective rasterization, The Visual Computer, 25, 5-7 (May), 549-557, 2009
    [Yu et al. 2010] Yu, J., McMillan, L., and Sturm, P., Multi-perspective modelling, rendering and imaging, Computer Graphics Forum, 29, 1, 1, 227 246, 2010
    [Zomet et al. 2003] Zomet, A., Feldman, D., Peleg, S., and Weinshall, D., Mosaicing new views: the crossed-slits projection, IEEE Transactions on Pattern Analysis and Machine Intelligence, 25, 61741-754, 2003
    [Zwicker 2006] Zwicker, M., Matusik, W., Durand, F., and Pfister, H., Antialiasing for automultiscopic 3d displays, In Rendering Techniques 2006: 17th Eurographics Workshop on Rendering, 73-82, 2006


BRIEF SUMMARY OF EXEMPLARY EMBODIMENTS
Methods and systems for stereoscopic content generation with granular control over binocular disparity based on multi-perspective imaging from representations of light fields are provided. To create the stereoscopic content with desired binocular disparity properties, an optimal set of light rays from a representation of a light fields is selected. The optimal set of light rays are selected by determining non-planar, piecewise continuous surface cuts through the representation of the light field. In some embodiments, the result is optimized by minimizing an energy reflecting prescribed parameters such as depth budget, maximum disparity gradient, and desired stereoscopic baseline. The methods and systems may be used for efficient and flexible stereoscopic post-processing, such as reducing excessive disparity while preserving perceived depth or retargeting of already captured scenes to various view settings.
In some embodiments, a method of generating stereoscopic content from a 3D representation of a light field is provided. The method typically includes capturing a 3D representation of a light field representing light rays from a scene using a set of standard perspective RGB (or other color space) input images, determining a plurality of 2D subsets of the 3D representation of the light field, wherein at least one 2D subset corresponds to a non-planar surface in the 3D representation of the light field, and generating, by a computer system, a stereoscopic image pair using the plurality of 2D subsets of the 3D representation of the light field.
In certain embodiments, the representation of the light field includes a 3D epipolar plane image (“EPI”) volume, the EPI volume having a plurality of horizontally parallel EPIs, each EPI having a plurality of vertically parallel scanlines of the set of input images, the plurality of scanlines corresponding to a particular scanline of the set of input images. In some embodiments, determining the plurality of 2D subsets of the 3D representation of the light field includes determining a first image of the stereoscopic image pair by selecting a reference image from the set of input images, the reference image corresponding to a 2D planar surface of the EPI volume, the 2D planar surface being perpendicular to the plurality of EPIs of the EPI volume, defining a set of binocular disparity goals for each pixel of a second image of the stereoscopic image pair with respect to the corresponding pixel of the first image of the stereoscopic image pair, and determining the second image of the stereoscopic image pair that corresponds to a non-linear surface of the 3D representation of the light field based on the set of binocular disparity goals, the non-linear surface corresponding to a 2D subset of the EPI volume perpendicular to each EPI of the EPI volume. In some embodiments, for light rays of each EPI, a set of pixel positions from the set of input images corresponding to each light ray, a light ray of an EPI for each pixel of the first image based at least upon the set of pixel positions, and the corresponding pixel of the second image based at least upon the set of binocular disparity goals and the light ray of the EPI for each pixel of the first image is determined.
In many embodiments, one or more additional features may be included to provide for additional functionality. For example, the set of input images are acquired using at least 2D subsets of one or more representations of a light field with three or higher dimensions, a camera mounted to a linear stage, a linear camera array, microlens array, or corresponding rendering of virtual scenes. In some other embodiments, a set of binocular disparity goals may be defined for each sub-pixel or group of pixels (e.g. a patch of 64×64 pixels).
In certain embodiments, two or more 2D subsets of the 3D representation of the light field correspond to parallel non-planar surfaces in the 3D representation of the light field, where individual sets of binocular disparity goals is defined for each pixel of the two or more 2D subset of the of the 3D representation of the light field with respect to the corresponding pixel of the reference image of the stereoscopic image pair.
In other embodiments, the stereoscopic image pair is generated using an energy minimization problem formulation that includes a content adaptive smoothness term to minimize visual artifacts of the generated stereoscopic image pair based on a smoothness constraint, wherein the smoothness constraint is adjusted for each region of the generated stereoscopic image pair based on visual saliency of the regions.
In other embodiments, a solution for the formulated energy minimization problem is determined using a graph cut optimization with binary s-t-cuts. In certain embodiments, the generated stereoscopic image pair is refined using at least one of linear, nonlinear, or gradient-based binocular disparity remappings. In some embodiments, the generated stereoscopic image pair is refined based on a manual control of binocular disparities.
According to some embodiments, a method of deferred rendering of computer-generated content is provided. Such method includes generating a depth map for each input image of a set of computer-generated input images of a scene, generating a depth volume from the generated depth maps of the set of input images, the depth volume corresponding to a 3D representation of a light field representing light rays from the scene, determining a plurality of 2D subsets of the depth volume corresponding to a plurality of 2D subsets of the 3D representation of the light field, where at least one 2D subset corresponds to a non-planar, surface of the depth volume, determining a plurality of input images from the set of input images included in the 2D subsets of the depth volume, rendering, by a computer system, the plurality of input images; and generating the plurality of 2D subsets of the 3D representation of the light field.
In certain aspects, a non-transitory computer-readable medium containing program instructions that, when executed by a computer, generates stereoscopic content from a representation of a light field. The program code for capturing a 3D representation of a light field representing light rays from a scene using a set of input images includes program code for determining a plurality of 2D subsets of the 3D representation of the light field, wherein at least one 2D subset corresponds to a non-planar surface in the 3D representation of the light field, and program code for generating, by a computer system, a stereoscopic image pair using the plurality of 2D subsets of the 3D representation of the light field. In some embodiments, the representation of the light field includes a 3D EPI volume, the EPI volume having a plurality of horizontally parallel EPIs, each EPI having a plurality of vertically parallel scanlines of the set of input images, the plurality of scanlines corresponding to a particular scanline of the set of input images.
In some embodiments, the non-transitory computer-readable medium containing contains program code for determining a first image of the stereoscopic image pair by selecting a reference image from the set of input images, the reference image corresponding to a 2D planar surface of the EPI volume, the 2D planar surface being perpendicular to the plurality of EPIs of the EPI volume, program code for defining a set of binocular disparity goals for each pixel of a second image of the stereoscopic image pair with respect to the corresponding pixel of the first image of the stereoscopic image pair, and program code for determining the second image of the stereoscopic image pair that corresponds to a non-linear surface of the 3D representation of the light field based on the set of binocular disparity goals, the non-linear surface corresponding to a 2D subset of the EPI volume perpendicular to each EPI of the EPI volume.
The following detailed description together with the accompanying drawings will provide a better understanding of the nature and advantages of the present invention.



BRIEF DESCRIPTION OF THE DRAWINGS
The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee.
In order to more fully understand the present invention, reference is made to the accompanying drawings. Understanding that these drawings are not to be considered limitations in the scope of the invention, the presently described embodiments and the presently understood best mode of the invention are described with additional detail through use of the accompanying drawings.
FIG. 1 illustrates an operation of a system for generating multi-perspective stereoscopic output images from a 3D light field; FIG. 1 comprises FIGS. 1(a)-1(d).
FIG. 2 is a simple example that illustrates a light field parameterization; FIG. 2 comprises FIGS. 2(a) and 2(b).
FIG. 3 illustrates examples of a planar single perspective surface cut of the light field and a non-planar multi-perspective surface cut of the light field; FIG. 3 comprises FIGS. 3(a)-3(b).
FIG. 4 illustrates an example 2D EPI of a light field with two planar u-t-cuts and the corresponding stereoscopic image pair generated from Is and Is′.
FIG. 5 illustrates an example 2D u-s-slice of a normalized disparity volume D and an example 2D u-s-slice of a true image disparity volume with respect to a reference view Is; FIG. 5 comprises FIGS. 5(a) and 5(b).
FIG. 6 is an illustration of multi-perspective light field surface cuts for changing stereoscopic disparity; FIG. 6 comprises FIGS. 6(a)-6(d).
FIG. 7 illustrates an example 1D slice of a 2D goal disparity map G and the corresponding difference volume of Ts−G illustrated as an unsigned function.
FIG. 8 illustrates the effects of disocclusions and occlusions with respect to goal disparities; FIG. 8 comprises FIGS. 8(a) and 8(b).
FIG. 9 illustrates examples of stereo pair with large baseline, stereo pair with small baseline, and nonlinear disparity remapping; FIG. 9 comprises FIGS. 9(a)-9(h).
FIG. 10 illustrates examples of gradient-based disparity remapping; FIG. 10 comprises FIGS. 10(a)-10(c).
FIG. 11 illustrates an example gradient domain compression; FIG. 11 comprises FIGS. 11(a)-11(e).
FIG. 12 illustrates an example of artistic control over disparity constraints; FIG. 12 comprises FIGS. 12(a) and 12(b).
FIG. 13 illustrates multiple views (i.e. slices of the light field), along with examples of optimized and unoptimized stereo images for multi-view autostereoscopic displays; FIG. 13 comprises FIGS. 13(a)-13(c).
FIG. 14 illustrates an example hardware setup for performing techniques described herein.



DETAILED DESCRIPTION OF EXEMPLARY EMBODIMENTS
Methods and systems for generating stereoscopic view from a light field are described herein. Such stereoscopic image generation is useful, as part of graphics rendering, for creating convincing yet perceptually pleasing stereoscopic content using novel stereoscopic processing and post-processing techniques. The generated stereoscopic content using these novel techniques can be particularly useful in context of multi-view autostereoscopic displays.
Many of embodiments described herein allow for generating image pairs for stereoscopic viewing with accurate control over binocular disparity, such as range and gradient, based on multi-perspective imaging from a light field, which is an improvement over prior approaches. The accurate control over binocular disparity allows for per-pixel control of binocular disparity of the generated stereoscopic content. For example, the binocular disparity for each neighboring pixel can be aligned based on a predefined range, in accordance with some embodiments. The per-pixel control over binocular disparity provides flexible and efficient stereoscopic post-processing, such as rendering excessive disparity while preserving while preserving the depth perception.
Lack of proper stereoscopic tools, particularly post-processing tools, has led to usage of tradition monoscopic tools and workflows for stereo-specific tasks [Mendiburu 2009], which are not well-suited for stereoscopic domain. As the result, a variety of prototypes have been developed for light field acquisition [Adelson and Wang 1992; Yang et al. 2002; Ng et al. 2005; Wilburn et al. 2005; Georgiev et al 2006; Veeraraghavan et al. 2007], which enable post-acquisition content modification such as depth-of-field (“DOP”), focus, or viewpoint changes. Although some of these developed prototypes offer post-acquisition content modification, the concept of post-acquisition control and editing is missing in stereoscopic processing. The methods and systems described herein can be used in the context of post-processing pipeline for stereoscopic content creation and editing.
Image Generation from a Light Field
As described above, binocular disparity is one of the most important visual cues that human vision uses to determine the depth perception. The binocular disparity is defined as the difference in location of an object in a scene, due to the eyes' horizontal separation, as seen by each left and right eye. As the result, the generated stereoscopic images of a scene have to only feature horizontal parallax, to simulate the eye's horizontal separation, without any vertical displacement of scene points between the generated stereoscopic images.
[Seitz 2001] showed that in order to only provide horizontal parallax, the generated stereoscopic images have to be constructed from a very specific three-parameter family of light rays. This allows for processing and generating stereoscopic images from a 3D light field instead of four-dimensional (“4D”) of higher dimensions as described by [Adelson and Bergen 1991; Levoy and Hanrahan 1996; Gortler et al. 1996].
FIG. 1 illustrates an operation of a system for generating multi-perspective stereoscopic content from a light field. FIG. 1 comprises FIGS. 1(a)-1(d). FIG. 1(a) shows multiple standard perspective RGB images that are used to create a 3D light field. FIG. 1(b) depicts a 3D light field that is created from the set of standard perspective RGB images. FIG. 1(c) illustrates multi-perspective surface cuts from the 3D light field in order to generate stereoscopic content from the light field. FIG. 1(d) illustrated an example stereoscopic content output generated as the result of the multi-perspective cuts, where the generated stereoscopic content shows an emphasis on the depth of local objects of a scene while suppressing the depth perception in the rest of the scene.
In some embodiments, each light ray of the 3D light field L: R3→R3 is parameterized by three parameters, L(u, v, s). Specifically, parameter s may correspond to the one-dimensional (“1D”) positional degree of freedom of the ray origin, whereas parameters (u, v) represent the ray direction. For example, FIG. 1 shows an actual 3D light field in the form of an epipolar plane image (“EPI”) volume, as described in [Gortler et al. 1996], which can be intuitively interpreted as a stack of 2D input images. Since the capture process naturally results in a discrete set of rays, each of the parameters u, v, and s may be implicitly treated as integers. Therefore, s may correspond to an index to one of the input images, while (u, v) indexes a pixel in image Is, i.e., L(u, v, s)=Is(u, v).
The parameterization of the light field is illustrated in FIG. 2 with a simple example. FIG. 2 comprises FIGS. 2(a) and 2(b). FIG. 2(a) depicts a 2D image of a scene and the corresponding imaging setup to generate a light field. Assuming a uniform sampling of the ray space with respect to these parameters, FIG. 2(b) illustrates a 2D light field corresponding or EPI to FIG. 2(a), where each point in ray space corresponds to a ray in the light field. Scene points seen in multiple images become EPI lines in ray space (see FIG. 1 or 3). The slope of each EPI line is proportional to the distance of the corresponding scene point. As illustrated in FIG. 2, a vertical EPI line may be vertical for the black point at infinity.
Stereoscopic Light Field Surface Cuts
According to some embodiments, a 2D standard (single) perspective or multi-perspective view can be generated by selecting a 2D subset of rays of the light field, L. As illustrated in FIG. 3(a), a planar u-v-slice or 2D cut at a particular parameter position s extracts the original input image, Is, with standard perspective. Surface cuts made by varying parameter s can yield images with varying centers of projection, in accordance with some embodiments. A planar v-s-cut with constant parameter u results in a so-called “pushbroom” panorama, which corresponds to a sensor with a single pixel column and a linearly varying position of the camera center [Yu et al. 2010]. As further illustrated in FIG. 3, a u-s-cut represents a single EPI, i.e., a 2D stack of the same scan line across all input images.
Any 2D subset of rays, other than planar cuts, with a certain ray coherence can be used to generate “meaningful” images. As described by [Seitz 2001] and [Peleg et al. 2001], multi-perspective images that feature only horizontal parallax can be fused stereoscopically. This would allow for generation of multi-perspective stereoscopic image pairs with controlled disparity by computing corresponding surface cuts through a light field. For example, FIG. 1(c) illustrates piecewise continuous surface cuts through an EPI volume shown in red that may be generated, which enable per-pixel disparity control.
In order to produce an output image from a light field, light rays lying on the cut surface may be sampled through parameterization of a cut through the light field. In some embodiments, multiple simultaneous surface cuts of the light field may be sampled through parameterization for stereoscopic image generation. As illustrated in FIG. 2(a), different parameters s represent input images captured at different, linearly translated camera positions. Correspondingly, the difference Δ(s′, s)=s′−s is proportional to the camera baseline between two images Is′ and Is, i.e., b=Δ(s′, s)/γ, where γ depends on the distance between the cameras, the field-of-view, and the horizontal image resolution. Hence, a stereoscopic image pair with baseline b can be generated by picking a reference view Is and selecting the second view at s′=s+γb corresponding to two parallel u-v-cuts through L. The convergence for such a stereoscopic image pair can be modified by shifting Is horizontally with respect to Is, as illustrated in FIG. 4.
In some embodiments, a 3D disparity volume D: R3→R+ is defined to provide maximum control of the perceived depth when viewing a stereo image pair generated from L.
As shown in FIG. 5(a), the 3D disparity volume stores the scaled reciprocal of the distance to the corresponding scene point for each ray in L. The 3D disparity volume may be created from standard depth or disparity maps and can be interpreted as a normalized disparity, such that the actual image disparity for a pixel p in Is′, to a reference image Is can be formulated as shown in Equation 1, where p is a shorthand notation for the coordinate pair (u, v).

Ts(p,s′)=Δ(s′,s)D(p,s′)  (Eqn. 1)

An example of the actual image disparity for a pixel p in Is′, to a reference image Is, i.e. the true disparity volume Ts for a particular view Is, is illustrated in FIG. 5b. 

Goal-Based Multi-Perspective Cuts

Having a reference view Is and the true disparities Ts, a second view Is′ such that Ts (*, *, s′) may be found that does not exceed a certain disparity range. In some embodiments, the distance Δ(s′, s) between the cuts, for example Is and Is′, may be adjusted to ensure the proper disparity range. In some other embodiments, the proper disparity range may be adjusted for only a particular point of a scene rather than changing the depth impression of the whole scene through determining non-planar and multi-perspective views that satisfy more general, content-dependent disparity constraints.
FIG. 6(a) shows a normalized disparity volume d and two images Is and Is′′. According to some embodiments, the horizontal parallax or image space disparity d of a pixel in Is to the corresponding pixel in Is′ can be computed as Ts(p′, s′)=Δ(s, s′)D(p, s). FIG. 6(b) illustrates a modified stereo pair that features a different depth impression only for the particular scene point seen at Is(p) by keeping s and s′ fixed and updating the actual scene depth, D(p, s). The actual scene depth may by updated through deformation of the actual geometry of the scene, which modifies the disparity d to d′ and changes the slope of the corresponding lines in the EPI volume. However, updating the depth of a scene implies changes to the complete underlying light field since changing the depth of a scene point changes the slope of the corresponding line in a ray space. As can be seen in FIG. 6(c), changing depth of the orange region results in different disocclusion patterns, with missing information in the light field depicted as red and blue regions. Since the corresponding light rays for these regions have not been captured in the original light field, completing these regions would require complex resampling and hole-filling operations on the light field.
According to certain embodiments, a non-planar surface cut sC: R2→R through the light field is determined, as illustrated in FIG. 6(d), which maps rays p to parameters s in order to meet a certain set of goal disparity constraints. The corresponding multi-perspective image IsC generated from the non-planar cut sC effectively results in the same change of disparity from d to d′ as other approaches previously described without deforming the actual geometry of the scene. A second view satisfying the disparity constraint for pixel p with respect to reference image Is can be generated from a cut sC that intersects the EPI line corresponding to Is(p) at parameter position u+d′. The cut sC, corresponding to the second view, selects a pixel from some input image for each pixel Is(p) such that the desired disparity constraints are fulfilled. As the result, the cut sC produces a multi-perspective output image IC that, together with the reference view Is, forms a stereoscopic image pair effectively controlling the camera baseline per-pixel.
A 2D map, G: R2→R, is defined as a set of goal disparities that, for each pixel of the output view IC, defines the desired disparity with respect to the reference view Is as follows. For example, if the disparity of pixel u in Is to the multi-perspective image IC should be changed to d′, the value of the goal disparity map at position u+d′ has to be set to G(u+d′, v=d′. In order to create a corresponding stereoscopic image pair, the goal disparity map has to be defined as shown by Equation 2, which can be constructed by iterating over all pixels u in the reference image Is.

G(u+φ(D(u,v,s)),v)=φ(D(u,v,s))  (Eqn. 2)

In Equation 2, φ: R→R may be a disparity mapping function that defines how the normalized disparity D may be mapped to a new disparity range.

Considering that the true disparity volume Ts(u, v, s′) represents the actual disparity of a point (u, v, s′) with respect to Is, the difference volume Ts(u, v, s′)−G (u, v) then represents the deviation of a pixel's disparity from the desired goal disparity. In some embodiments, this deviation representation is used to determine a cut sC that passes close to the zero set of this difference volume in order to satisfy the disparity goal for each point of the scene. For example, a difference volume Ts(u, v, s′)−G (u, v) and a vertically stretched 1D slice of the 2D goal disparity map G are illustrated in FIG. 7. As can be seen, the difference volume Ts−G is shown as an unsigned function, which represent deviation of each point in the light field from the desired disparity. The cut, sC, can be determined such that it passes close to the zero set of the difference volume. This would lead to generation of the resulting image IC and Is that may form a multi-perspective stereoscopic image pair with the desired goal disparity.
The construction of G is neither surjective nor injective due to occlusions and disocclusions in the scene. In other words, disparity constraints cannot define for regions that are occluded in Is but most likely disoccluded in an output view IC. Since a reference image Is contains depth information for scene elements visible in Is, the construction of G by forward mapping of remapped disparities φ(D(*, *, s)) for G construction is not surjective, which leads to undefined segments in G, such as left blue segment as shown in FIG. 8. In order to show an undistorted standard perspective view of all undefined disoccluded regions, a certain smoothness may be imposed on the cut sc, as the undefined regions generally only cover a small number of pixels, in accordance to some embodiments. Similarly, the construction of G is not infective due to visible regions in Is that may be occluded in other views. Since differently remapped disparities of close and distant objects that may compete for the same slot in G, such as overlapping orange and pink region as shown in FIG. 8(b), the disparity constraints for the closer object may be selected and stored.
Disparity Remapping
Linear remapping of the disparity range corresponds to changing the camera baseline between two standard perspective views. As shown in FIG. 9(b), two cuts are selected closer than the cuts shown in FIG. 9(a) in order to produce a decreasing baseline. This way, the quite abstract and unintuitive concept of the camera baseline is eliminated, and one can directly specify the desired goal disparity range of the output images as the image disparities resulting from a particular camera baseline are dependent on the scene content.
When nonlinear changes of the disparity space are required, arbitrary remapping functions φ can be applied to construct the desired goal disparity volume and even constant disparities are possible. For example, φ could be any of the nonlinear disparity mapping operators introduced by [Lang et al. 2010] for display adaptation, stereoscopic error correction, or artistic effects. These functions can act globally on the complete domain as well as locally by remapping disparity gradients.
FIG. 9(a) shows a standard stereo pair with a large baseline where the foreground provides a good impression of depth. Typically, the large baseline causes quite large disparities in the background region of the image that can lead to ghosting artifacts or even the inability to stereoscopic fuse, when viewed on a larger screen. Although decreasing the baseline (i.e. reducing the distance of cuts) would reduce the ghosting problems with respect to the background, as shown in FIG. 9(b), but also considerably flattens the foreground of the scene. In such cases, nonlinear disparity mapping function can enhance the depth impression of the foreground, while keeping the maximum disparities in the background bounded as shown in FIG. 9(c). As illustrated by the histograms and graphs of FIG. 9(d)-(h), the depth between the foreground and the background is preserved using the non-linear disparity remapping while the depth is lost in case of linear disparity remapping.
In some embodiments, the final goal disparity may be set using the gradient based remapping. For gradient based remapping, gradients in x and y directions are computed and processed non-uniformly using the gradient magnitudes, for example to suppress big disparity jumps. Then, the height field is reconstructed by integration of the gradients using a Poisson Solver. As shown in by examples of FIG. 10, non-linear disparity gradient remapping may reduce the overall disparity range while preserving the perception of depth discontinuities and local depth variations. The top portion of each FIG. 10(a)-(c) shows the stereo image pair with two perspective images and a fixed baseline while the bottom portion of each FIG. 10(a)-(c) depicts the end result after non-linear disparity remapping. As expected, elements closer to camera (e.g. car or plane) would keep the depth perception while the background ghosting effect has been reduced.
In some embodiments, scene elements with strong negative disparity may be compressed using the gradient remapping function. For example as shown in FIG. 11, the partially cropped couple feature strong negative disparity, which results in a window violation, as described by [Mendiburu 2009]. Since changing the convergence would increase the background disparities, the scene elements with strong negative disparity may be compressed in order to resolve the disparity problem. This way, the scene element with strong negative disparity would be pushed closer to the screen while keeping the background disparities unchanged. By comparing the histograms of the FIG. 11, it is evident that the background and foreground disparities are compressed after gradient domain remapping.
In some embodiments, a concise manual control of disparities is provided, which is an important requirement in any stereoscopic production environment. An artist may construct a desired goal disparity map by modifying the depth map of the reference view using existing 2D segmentation and brush tools, which allows for interesting artistic effects as well as fine-scale correction of the stereoscopic impression of a scene. FIG. 12 shows examples of manual control over disparity. In some other embodiments, an artist may manually “mask out” certain undesirable regions of a light field, e.g., for removing certain visually distracting effects from a stereo pair, such as monocular specular reflections. For example, FIG. 12(b) illustrates disparity map and histogram for a resulting stereo pair after an artist's modifications. As used herein, it should be understood that “artist” refers to a user who is generating an image and may encompass someone just creating an image for their own use, a professional creating images for some sort of image production, or other types of users that desire to create images.
Formulation as an Energy Minimization Problem
A wide variety of problems may be solved in terms of energy minimization problem. Such problems may be formulated as an energy minimization problem, where the solution typically corresponds to a maximum posteriori estimate of a solution. The value of the formulated function represents the energy of a system that is being modeled.
Accordingly, the problem of determining the proper cut sC may be formulated as an energy minimization problem. The energy measuring the deviation of a 2D cut sC may be expressed as shown in Equation 3.




  
    
      
        
          
            E
            d
          
          ⁡
          
            (
            
              s
              C
            
            )
          
        
        =
        
          
            ∑
            p
          
          ⁢
          
            
              
              
                
                  
                    T
                    s
                  
                  ⁡
                  
                    (
                    
                      p
                      ,
                      
                        
                          s
                          C
                        
                        ⁡
                        
                          (
                          p
                          )
                        
                      
                    
                    )
                  
                
                -
                
                  G
                  ⁡
                  
                    (
                    p
                    )
                  
                
              
              
            
            .
          
        
      
    
    
      
        (
        
          Eqn
          .
          
            
          
          ⁢
          3
        
        )
      
    
  




Although a cut sC computed from this data term closely follows the prescribed goal disparities, it does not enforce any coherence between neighboring output rays or output pixels. The lack of coherence between output image points may cause visual artifacts in noisy or ambiguous estimates of Ts. Such artifacts are particularly noticeable in highly textured regions or at depth discontinuities. Therefore, an additional content-adaptive smoothness term may be introduced to the formulated energy minimization problem, in accordance to some embodiments.
The smoothness constraint may be applied differently to each region of an image. In some embodiments, a higher smoothness is enforced to increase the coherence of the rays selected by sC in the proximity of visually salient parts of an image, such as depth discontinuities and highly textured regions. In particular, a higher saliency may be assigned to scene elements close to the camera and cut through more distant regions. Conversely, the smoothness constraint may be more relaxed in order to increase the flexibility of the cut to perform multi-perspective view transitions in the light field. The energy formulation reflecting the described properties may be expressed as shown in Equation 4, subject to the equalities of Equations 5 and 6, where Nu and Nv are the sets of all neighboring pixels along the u-axis and v-axis, respectively.




  
    
      
        
          
            E
            s
          
          ⁡
          
            (
            
              s
              C
            
            )
          
        
        =
        
          
            
              ∑
              
                
                  (
                  
                    p
                    ,
                    q
                  
                  )
                
                ∈
                
                  N
                  u
                
              
            
            ⁢
            
              
                
                
                  
                    
                      s
                      C
                    
                    ⁡
                    
                      (
                      p
                      )
                    
                  
                  -
                  
                    
                      s
                      C
                    
                    ⁡
                    
                      (
                      q
                      )
                    
                  
                
                
              
              ⁢
              
                
                  p
                  u
                
                ⁢
                
                  (
                  *
                
                )
              
            
          
          +
          
            
              ∑
              
                
                  (
                  
                    p
                    ,
                    q
                  
                  )
                
                ∈
                
                  N
                  u
                
              
            
            ⁢
            
              
                
                
                  
                    
                      s
                      C
                    
                    ⁡
                    
                      (
                      p
                      )
                    
                  
                  -
                  
                    
                      s
                      C
                    
                    ⁡
                    
                      (
                      q
                      )
                    
                  
                
                
              
              ⁢
              
                
                  p
                  v
                
                ⁢
                
                  (
                  *
                
                )
              
            
          
        
      
    
    
      
        (
        
          Eqn
          .
          
            
          
          ⁢
          4
        
        )
      
    
  




pu(*)=min(Pmax,|∂sD(*)|+λD(*)+κ|∂sL(*)|)  (Eqn. 5)

pv(*)=min(pmax,|∂sD(*)|+λD(*)+κ|∂uL(*)|)  (Eqn. 6)

In Equations 4, 5 and 6, (*) stands for (p, sC(p)). The term |sC (p)−sC (q)| penalizes variation of the cut sC along the s-axis, i.e., view transitions. This penalty is weighted by the content-adaptive terms pu(*) and pv(*), respectively.

For both axes, the weighting depends on depth discontinuities as ∂sD and the absolute normalized disparity D. Depending on the sampling of the light field, noticeable disparity jumps may be introduced by view transitions to an adjacent view for scene elements very close to the viewer. In some embodiments, increasing smoothness for those regions and over depth discontinuities effectively moves view transitions to less noticeable regions in the background. In some other embodiments, the smoothness for other regions of the image with different saliency may also be increased in order to encourage view transitions in less salient regions.
The smoothness of the cut sC may be controlled effectively using the depth-based terms introduced in equation 4. For example, change of radiance between different input images Is may be taken into account for the u-axis, while jumps of the cut in the proximity of vertical image edges of v-axis may be penalized. In some embodiment, similar to the concept of robust nonlinear error functions, the maximum penalty pmax ensures that the cut sC can be discontinuous when required by the defined goal disparity.
According to some embodiments, constants are introduced in equation 4 to bring all terms to a similar scale. For example, defining λ=0.05, κ=0.1, and pmax=0.3 would cause the final energy to be defined as shown in Equation 7, where k=250.

E(sC)=Ed(sC)+kEs(sC)  (Eqn. 7)

According to some embodiments, increasing the smoothness constraint through increasing the value of k may produce “flattened” cuts, i.e. output images closer to a standard perspective image. Such increase simply makes the output image to simply fall back to the original input images without compromising the image quality.

Optimization via Graph Min-Cuts

The energy minimization problem as described by Equation 7 may be solved using the graph cut optimization using, for example, standard procedure for binary s-t cuts. The graph used for optimization may comprise a set of nodes that correspond to pixels or some other feature and a set of edges connection the nodes. Such graph may also include some additional special nodes of the graph called terminals corresponding to a set of labels that can be assigned to pixels. The terminals are typically called the source s and the sink t. An s-t cut on a graph with two terminals would lead to partitioning of the nodes of the graph into two disjoint subsets of S and T, respectively containing the source s and the sink t.
According to some embodiments, employing the standard procedure for binary s-t cuts as follows:
1. For n input images of dimension w×h, construct a 3D regular graph of size w×h×(n+1).
2. Associate a ray at position (u, v, s′) with a directional graph edge between the corresponding two nodes along the s-axis, and choose the edge weight as |Ts(u, v, s′)−G (u, v)|.
3. Weight bi-directional edges between neighboring nodes along the w axis and v-axis with the corresponding smoothness values kpu and kpv respectively.
4. Connect boundary nodes corresponding to parameters s=0 and s=n to the source and sink of the graph, respectively, with infinite weights.
As the result, the minimum s-t cut of this graph may yield the desired cut surface sC that minimizes Eq. (5). In some embodiments, the cut sC may be determined by a modification of the s-t cut algorithm and/or energy minimization formulation. For example, additional penalty edges for enforcing the cut) (C0) continuity may be used as described by [Rubinstein et al. 2008]. As another example, piecewise continuous cuts may be used to provide more flexibility due to the support for sudden view transitions.
According to some other embodiments, different methods for solving the formulated minimization energy problem may be used, such as an alternative formulation may be used on multi-labeling via a-expansion, where each label is associated with a particular u-v-slice along the s-axis of the EPI volume.
Extensions and Applications
Above, details have been described for a basic methods and systems for generating stereoscopic view from the light field, which includes formulation of energy minimization problem and optimizing it via graph min-cuts. These embodiments are useful in many situations as-is. However, stereoscopic image generation from a light field often has unique, specialized requirements, and the disclosed embodiments might not cater to all situations. However, the additional features described herein can provide a variety of embodiments that may satisfy the needs of some particular situations.
N-View Stereo from Multiple Cuts
In some embodiments, two or more multi-perspective cuts may be created instead of creating a stereo pair consisting of a standard perspective image Is and a multi-perspective image IC. For example, where the goal disparities are evenly distributed to both views and the reference view is centered between the two multi-perspective output cuts, two goal disparity maps GL and GR may be defined and expressed as equations 8 and 9.

GL(u−1/2φ(D(u,v,s)),v)=−1/2p(D(u,v,s))  (Eqn. 8)

GR(u+1/2φ(D(u,v,s)),v)=+1/2p(D(u,v,s))  (Eqn. 9)

It should be understood that more than two views may be created and handled in an identical manner.

While separately defining a goal disparity map for each view provides high flexibility, many application scenarios, such as multi-view autostereoscopic displays, often require simply a linear change of disparity between views. This can be exploited for an efficient interpolation based algorithm to generate multiple views with just the reference view s and one multi-perspective cut sC. For example, sC may be determined using a mapping function φ(D(u, v, s)) while the two views s and sC may be converted into n views with linearly interpolated disparities. Using Equation 2, a computer process operate knowing that the goal disparities of view kε[0, n] may be expressed as in Equation 10.

G(u+−k/nφ)(D(u,v,s)),v)=−k/nφ)(D(u,v,s))  (Eqn. 10)

In this example, a cut sC may contain the interpolated points of all EPI lines connecting corresponding points of s and sC.

Multi-View Autostereoscopic Displays
A particularly useful application is to compute image sets for multi-view autostereoscopic displays, which generally support only a very restricted disparity range without leading to excessive ghosting artifacts. Similar to stereoscopic displays, the multi-view autostereoscopic displays have a limited depth budget. In order to avoid inter-perspective aliasing, an input light field may be pre-filtered and remapped to the available spatio-angular display bandwidth as described by [Zwicker et al. 2006] to optimized the output image displayed. Properly remapped data may be obtained to drive a particular automultiscopic display by determining multiple surface cuts through a light field volume in order to create n-view stereoscopic datasets. For example, an 8-view autostereoscopic display from Alioscopy is illustrated in FIG. 13.
FIG. 13(a) illustrates multi-perspective 8-view stereo that is depth-enhanced and optimized with respect to the disparity requirements of an actual multi-view autostereoscopic display. FIG. 13(b) shows an unoptimized content with considerable ghosting artifact. In contrast, FIG. 13(c) illustrates an optimized n-view stereo images that are designed to meet the disparity requirements of the output device while enhancing perceived depth our method can automatically compute.
Stereoscopic Video Processing
It is generally desirable to enforce a certain continuity between two cuts at consecutive time steps in order to process video. Currently, a temporal dimension is added to the graph structure to enforce temporal smoothness. Each time step has its own 3D subgraph, and corresponding nodes of subgraphs from consecutive time steps are connected via an additional edge. However, this leads to processing the whole 4D spatio-temporal light field volume at once.
In some embodiments, the concept of decaying influence of previous time steps on the data and smoothness terms of the current time step is introduced to avoid processing the whole 4D spatio-temporal light field volume at once. For example, et may denote the edge weight for a given time step t according to Equation (3) and Equation (4). During the update of the graph structure from time t−1 to t, for any edge, the temporally averaged edge weight may be set as shown by Equation 11. In some embodiments, α=0.9 may be used since the temporal evolution of a light field is quite coherent.

e′t=αet+(1+α)et-1  (Eqn. 11)

Deferred Rendering for Computer-Generated Content
Another particular useful application is in computer-generated content processing, such as 3D animation movies. Implementing multi-perspective camera models into the computer-generated rendering pipeline to meet the expectations of a director regarding control and flexibility is often a difficult problem. Also, changing (e.g. warping) the 3D geometry does not allow for arbitrary complex disparity constraints without compromising the scene composition, lighting calculations, or visual effects.
In some embodiments, deferred rendering method may be used to provide a practical solution for processing the computer-generated content to reduce rendering time for generating a light field for computer-generated content. This way, only the depth maps of the input views may be rendered and used to generate a disparity volume D, which is typically several orders of magnitude faster than rendering fully shaded color images and creating a light field from these images. Once the required set of input views is known from the surface cut sC of the disparity volume, those fully shaded color images corresponding to the surface cut sC or just the optimal set of light rays determined from the surface cut sC can be rendered and combined for computer-generated content processing. In some embodiments, even lower resolution proxy geometry could be used instead of the highly tessellated subdivision surfaces often used in rendering.
Different Light Field Parameterizations
In other application scenarios, it may be preferred to produce other forms of stereoscopic images, such as Omnistereo panoramas as discussed by [Peleg et al. 2001], or stereo Pushbroom panoramas and Cyclographs as discussed by [Seitz 2001]. For these types of images, the light field parameterization, light field cut, and sampling scheme may have to be reconsidered. For example, having a standard perspective reference image or sampling along the u-v-dimension of the light field may not be desirable in all scenarios.
As mentioned above, in our light field representation, a stereoscopic Pushbroom panorama simply corresponds to a v-s-cut instead of a u-v-cut. In other words, stereoscopic Pushbroom images may be created by mere swapping the dimensions u and s in our disclosed formulation. For Omnistereo panoramas and Cyclographs, the 3D light fields are essentially constructed with a rotating camera with a certain offset orthogonal to the rotation axis, yielding a u-v-a volume. In these scenarios, planar v-α slices can be used to produce stereoscopic panoramas.
Required Number of Input Images
The number and spacing of input images needed generally depend on the desired level of control over disparity, which is application dependent. For instance, a typical problem scenario in stereoscopic movie production is that parts of a sequence shot with a standard stereo camera violate certain disparity constraints, such as window violation, maximal disparity, etc. In order to fix such issues, a second view with the remaining scene may be manually composed, so that a few additional views of a scene would be able to resolve certain problems. As an example, as few as 4 input images may be used for the final output view for correcting excessive disparities and stereoscopic window violation. Generally, the number of images used for image generation rarely exceeded 10 images. One of the main reasons for such performance improvement is the normalized disparity volume D, which can in practice often be computed from a much sparser set of depth maps or images. Accordingly, the proposed methods and systems not only perform well with densely sampled light fields, but also scale well and are applicable to scenarios with a considerably smaller number of inputs views.
Variations
In addition to basic systems described herein, more complex systems might be implemented. As described above, per-pixel goal disparity constraints are mainly used to derive the energy minimization problem to determine a surface cut through the light field. In addition to the goal-based cut optimization mainly focused on disparity constraints, there might be more emphasis placed on other constraints, such as stereoscopic perception, visual saliency, or temporal coherence.
In order to provide a better performance, potential resampling issues that might arise may be addressed explicitly. To address the resampling issues, different method and techniques may be used. For example, gradient-based image reconstruction, gradient-domain cuts, sub-pixel resolution techniques, and methods for upsampling of the light field may be used.
Determining a cut through the light field is largely dependent on per-pixel goal disparity constraints. The disparity constraints are defined with respect to a reference image view. In another system variation, disparity constraints may be defined using pairwise constraints between neighboring views for regions occluded in the reference view.
Hardware Example
FIG. 14 illustrates an example system 100 for animation management according to an embodiment of the present invention. In the presently described embodiment, computer system 100 typically includes a monitor 110, computer 120, a keyboard 130, a user input device 140, computer interfaces 150, and the like. Images and frames can be generated, stored in memory or hard disk storage, and processed. The interfaces and/or memory might also be used to provide the metadata about an image, objects in the virtual space and other considerations.
In various embodiments, display/monitor 110 may be embodied as a CRT display, an LCD display, a plasma display, a direct-projection or rear-projection DLP, a microdisplay, or the like. In various embodiments, monitor 110 may be used to visually display images, or the like, as well as being part of an interactive environment.
In the present embodiment, user input device 140 is typically embodied as a computer mouse, a trackball, a track pad, a joystick, wireless remote, drawing tablet, voice command system, eye tracking system, and the like. Embodiments of network interface 150 typically include an Ethernet card, a modem (telephone, satellite, cable, ISDN), (asynchronous) digital subscriber line (DSL) unit, FireWire interface, USB interface, and the like. For example, network interface 150 may be coupled to a computer network, to a FireWire bus, or the like. In other embodiments, network interface 150 may be physically integrated on the motherboard of computer 120 and/or include software drivers, or the like.
In various embodiments, computer 120 typically includes familiar computer components such as a processor 160, and memory storage devices, such as a random access memory (RAM) 170, disk drives 180, and system bus 190 interconnecting the above components. RAM 170 or other memory might hold computer instructions to be executed by one or more processors as a mechanism for effecting some functionality described herein that is implemented in software. In one embodiment, computer 120 includes one or more CoreTM microprocessors from Intel. Further, in the present embodiment, computer 120 typically includes a UNIX-based operating system.
RAM 170 and disk drive 180 are examples of computer-readable tangible media configured to store embodiments of the present invention including computer-executable code implementing techniques described herein, data such as image files, object/scene models including geometric descriptions of virtual objects, images, camera paths, procedural descriptions, a rendering engine, executable computer code, and/or the like. Other types of tangible media may include magnetic storage media such as floppy disks, networked hard disks, or removable hard disks, optical storage media such as CD-ROMS, DVDs, holographic memories, and/or bar codes, semiconductor memories such as flash memories, read-only-memories (ROMS), battery-backed volatile memories, networked storage devices, and the like.
In various embodiments, computer system 100 may also include software that enables communications over a network such as the HTTP, TCP/IP, RTP/RTSP protocols, and the like. In alternative embodiments of the present invention, other communications software and transfer protocols may also be used, for example IPX, UDP or the like.
In some embodiments of the present invention, a graphical processor unit or “GPU”, may be used to accelerate various operations. Such operations may include determining image rendering, camera movements, view alterations, camera paths or other processes.
FIG. 14 is representative of a computer system capable of embodying the present invention. It will be readily apparent to one of ordinary skill in the art that many other hardware and software configurations are suitable for use with the present invention. For example, the computer may be a desktop, portable, rack-mounted or tablet configuration. Additionally, the computer may be a series of networked computers. Further, the use of other microprocessors are contemplated, such as XeonTM, PentiumTM or ItaniumTM microprocessors from Intel; TurionTM 64 or OpteronTM microprocessors from Advanced Micro Devices, Inc; and the like. Further, other types of operating systems are contemplated, such as VistaTM or Windows XPTM or the like from Microsoft Corporation, SolarisTM from Sun Microsystems, Linux, Unix, or the like.
In still other embodiments, the platform is a gaming system, such as Microsoft's Xbox 360TM game console, Sony's Playstation 3TM console, or Nintendo's WiiTM console.
In still other embodiments, the techniques described above may be implemented upon a chip or an auxiliary processing board. Many types of configurations for computational devices can be used to implement various methods described herein. Further, processing components having different levels of computational power, e.g., microprocessors, graphics processors, RISC processors, embedded processors, or the like can also be used to implement various embodiments.
As has now been described, methods and systems for generating stereoscopic content with per-pixel control over disparity based on multi-perspective imaging from light fields are provided. This might find uses in many fields. Merely by way of example, stereoscopic post-processing of a captured live action movie, scene composition modification of computer generated image content without requiring a change of scene geometry, for single view rendering of moving all or part of scene objects with occlusion and disocclusion handled in a proper way, for rendering views for multi-view autostereoscopic displays, and/or for advanced 2D movie editing effects. For instance, a cut/transition between two scenes, where one first replaces the background and then the foreground or vice versa (e.g., the actors from the previous scene remain visible, while the whole rest of the scene is gradually replaced, might be a situation where some of the teaching herein might find a use. With light field cut computation, one could selectively mix light rays of two or more different scenes and compose them together.
As described herein, per-pixel disparity control can be obtained, at least by synthesizing the entire output images from the light rays sampled in a light field. A mathematically valid framework is provided to locate each light ray needed to synthesize output images. The image generation can be configured to be represented by a statement of an optimization problem, which can then be solved by standard solvers.
Unlike other approaches, all operations might be done in the output images space, so that information loss and sampling issues are minimized, spatial distortion caused by insufficient input data is eliminated or reduced. Handling can then fully exploit a light field without a limitation on the number of output views. In some cases, an additional degree of freedom for the selection of rays is provided, which allows a vertical image column to be synthesized from multiple views. A more sophisticated ray selection strategy to cope with the visual smoothness through the optimization is then possible. In some embodiments, image generation can employ arbitrary disparity constraints.
Further embodiments can be envisioned to one of ordinary skill in the art after reading this disclosure. In other embodiments, combinations or sub combinations of the above disclosed invention can be advantageously made. The example arrangements of components are shown for purposes of illustration and it should be understood that combinations, additions, rearrangements, and the like are contemplated in alternative embodiments of the present invention. Thus, while the invention has been described with respect to exemplary embodiments, one skilled in the art will recognize that numerous modifications are possible.
For example, the processes described herein may be implemented using hardware components, software components, and/or any combination thereof. The specification and drawings are, accordingly, to be regarded in an illustrative rather than a restrictive sense. It will, however, be evident that various modifications and changes may be made thereunto without departing from the broader spirit and scope of the invention as set forth in the claims and that the invention is intended to cover all modifications and equivalents within the scope of the following claims.

## Claims

1. A method, using a computer, for generating stereoscopic content from a three-dimensional representation of a light field representing virtual or actual light rays from a scene, wherein the representation is readable in electronic form, the method comprising:
capturing the three-dimensional representation of the light field, wherein the light field includes a plurality of input images;
receiving a user-defined two-dimensional map including a set of goal binocular disparities for each pixel of a first image of a stereoscopic image pair with respect to a corresponding pixel of a second image of the stereoscopic image pair, wherein a first goal binocular disparity defined in the two-dimensional map for a first group of pixels is greater than a second goal binocular disparity defined in the two-dimensional map for a second group of pixels;
determining a non-planar surface cut through the three-dimensional representation of the light field, wherein the non-planar surface cut is determined using the received two-dimensional map;
determining a first two-dimensional subset and a second two-dimensional subset of the three-dimensional representation of the light field, wherein at least one of the first two-dimensional subset and the second two-dimensional subset corresponds to the non-planar surface cut; and
generating, using the computer, the first image of the stereoscopic image pair using the first two-dimensional subset and the second image of the stereoscopic image pair using the second two-dimensional subset.

2. The method of claim 1, wherein the plurality of input images includes standard perspective images.

3. The method of claim 1, wherein the plurality of input images are acquired using at least two-dimensional subsets of one or more representations of a light field with three or higher dimensions, a camera mounted onto a linear stage, a linear camera array, microlens array, or corresponding rendering of virtual scenes.

4. The method of claim 1, wherein some of the light rays represented by the representation of the light field are parameterized using at least three parameters, the at least three parameters defining the three-dimensional representation of the light field.

5. The method of claim 4, wherein the at least three parameters indicate a direction and one-dimensional positional degree of freedom of ray origin for each light ray.

6. The method of claim 1, wherein the three-dimensional representation of the light field is in form of a three-dimensional epipolar plane image volume, the epipolar plane image volume having a plurality of epipoloar plane images, each epipolar place image having a plurality of scanlines of a set of images from the scene.

7. The method of claim 6, wherein a height of the three-dimensional epipolar plane image volume corresponds to a horizontally parallel stack of the plurality of epipolar plane images, and wherein a depth of the three-dimensional epipolar plane image volume corresponds to a vertically parallel stack of the set of images from the scene such that each epipolar plane image includes a vertically parallel stack of a particular scanline of each image of the set of images from the scene.

8. The method of claim 1, wherein
at least one two-dimensional subset of the three-dimensional representation of the light field is a two-dimensional reference subset that corresponds to an input image.

9. The method of claim 1, wherein determining the two-dimensional subset that corresponds to the non-planar surface cut comprises:
determining a set of locations that corresponds to each light ray;
determining a light ray for each location of the set of locations within the non-planar surface cut of the three-dimensional representation of the light field based at least upon the set of locations that corresponds to each light ray; and
for each location within the two-dimensional reference subset, determining the corresponding location within the two-dimensional subset that corresponds to the non-planar surface cut based at least upon a set of binocular disparity goals and the determined light ray.

10. The method of claim 9, wherein the set of disparity goals are defined using at least one of linear, nonlinear, or gradient-based binocular disparity remappings.

11. The method of claim 9, wherein the set of disparity goals are defined based upon a manual modification.

12. The method of claim 9, wherein the stereoscopic image pair is generated using an energy minimization problem formulation.

13. The method of claim 12, wherein the energy minimization problem formulation includes a data term indicating the set of binocular disparity goals and the set of locations that corresponds to each light ray.

14. The method of claim 12, wherein the energy minimization problem formulation includes a content adaptive smoothness term to minimize visual artifacts of the generated stereoscopic image pair based upon a smoothness constraint.

15. The method of claim 14, wherein the smoothness constraint is adjusted for each region of the generated stereoscopic image pair based upon visual saliency of the regions.

16. The method of claim 15, wherein a solution for the formulated energy minimization problem is determined using a graph cut optimization with binary s-t-cuts.

17. The method of claim 1, wherein multiple subsets of three-dimensional representation of the light field with linear change in binocular disparity are determined using an interpolation based algorithm based upon at least the two-dimensional reference subset and at least one two-dimensional subset that corresponds to the non-planar surface cut in the three-dimensional representation of the light field.

18. A non-transitory computer-readable medium containing program instructions that, when executed by a computer, generates stereoscopic content from a three-dimensional representation of a light field, comprising:
program code for capturing a three-dimensional representation of a light field representing virtual or actual light rays from a scene;
program code for receiving a user-defined two-dimensional map including a set of goal binocular disparities for each pixel of a first image of a stereoscopic image pair with respect to a corresponding pixel of a second image of the stereoscopic image pair, wherein a first goal binocular disparity defined in the two-dimensional map for a first group of pixels is greater than a second goal binocular disparity defined in the two-dimensional map for a second group of pixels;
program code for determining a non-planar surface cut through the three-dimensional representation of the light field, wherein the non-planar surface cut is determined using the received two-dimensional map;
program code for determining a first two-dimensional subset and a second two-dimensional subset of the three-dimensional representation of the light field, wherein at least one of the first two-dimensional subset and the second two-dimensional subset corresponds to the non-planar surface cut; and
program code for generating, by a computer system, the first image of the stereoscopic image pair using the first two-dimensional subset and the second image of the stereoscopic image pair using the second two-dimensional subset.

19. The non-transitory computer-readable medium of claim 18, wherein some of the light rays represented by the representation of the light field are parameterized using at least three parameters, the at least three parameters defining the three-dimensional representation of the light field.

20. The non-transitory computer-readable medium of claim 18, wherein the at least three parameters indicate a direction and one-dimensional positional degree of freedom of ray origin for each light ray.

21. The non-transitory computer-readable medium of claim 18, wherein the three-dimensional representation of the light field is in form of a three-dimensional epipolar plane image volume, the epipolar plane image volume having a plurality of epiploar place images, each epipolar plane image having a plurality of scanlines of a set of images from the scene.

22. The non-transitory computer-readable medium of claim 21, wherein a height of the three-dimensional epipolar plane image volume corresponds to a horizontally parallel stack of the plurality of epipolar plane images, and wherein a depth of the three-dimensional epipolar plane image volume corresponds to a vertically parallel stack of the set of images from the scene such that each epipolar plane image includes a vertically parallel stack of a particular scanline of each image of the set of images from the scene.

23. The non-transitory computer-readable medium of claim 18, wherein the program code for determining the first two-dimensional subset and the second two-dimensional subset of the three-dimensional representation of the light field comprises:
program code for defining a set of binocular disparity goals for each location within at least one two-dimensional subset of the three-dimensional representation of the light field with respect to a two-dimensional reference subset of the three-dimensional representation of the light filed; and
program code for determining a non-planar surface of the three-dimensional representation of the light field that corresponds to the at least one two-dimensional subset of the three-dimensional representation of the light field based upon the set of binocular disparity goals.

24. The non-transitory computer-readable medium of claim 23, wherein the program code for determining the plurality of two-dimensional subsets of the three-dimensional representation of the light field further comprises:
program code for determining the two-dimensional reference subset of the three-dimensional representation of the light field, the two-dimensional reference subset corresponds to an input image of the three-dimensional representation of the light field.

25. The non-transitory computer-readable medium of claim 23, wherein the program code for determining the non-planar surface of the three-dimensional representation of the light field comprises:
program code for determining a set of locations that corresponds to each light ray;
program code for determining a light ray for each location of the set of locations within the non-planar surface cut of the three-dimensional representation of the light field based at least upon the set of locations that corresponds to each light ray; and
program code for determining the corresponding location within the two-dimensional subset of the three-dimensional representation of the light field for each location within two-dimensional reference subset based at least upon the set of binocular disparity goals and the determined light ray.

26. The non-transitory computer-readable medium of claim 23, wherein the set of disparity goals are defined using at least one of linear, nonlinear, gradient-based binocular disparity remappings, or a manual modification.

27. The non-transitory computer-readable medium of claim 18, wherein the stereoscopic image pair is generated using an energy minimization problem formulation.

28. The non-transitory computer-readable medium of claim 27, wherein the energy minimization problem formulation includes a content adaptive smoothness term to minimize visual artifacts of the generated stereoscopic image pair based upon a smoothness constraint.