# Patent

## Title
Method and apparatus for managing artifacts in frequency domain processing of light-field images

## Abstract

Various methods and apparatus for removing artifacts in frequency domain processing of light-field images are described. Methods for the reduction or removal of the artifacts are described that include methods that may be applied during frequency domain processing and a method that may be applied during post-processing of resultant angular views. The methods may be implemented in software as or in a light-field frequency domain processing module. The described methods include an oversampling method to determine the correct centers of slices, a phase multiplication method to determine the correct centers of slices, a method to exclude low-energy slices, and a cosmetic correction method.

## Background
BACKGROUND
Description of the Related Art
Conventional cameras fail to capture a large amount of optical information. In particular, a conventional camera does not capture information about the location on the aperture where different light rays enter the camera. During operation, a conventional digital camera captures a two-dimensional (2-D) image representing a total amount of light that strikes each point on a photosensor within the camera. However, this 2-D image contains no information about the directional distribution of the light that strikes the photosensor. Directional information at the pixels corresponds to locational information at the aperture.
In contrast, light-field cameras sample the four-dimensional (4-D) optical phase space or light-field and in doing so capture information about the directional distribution of the light rays. This information captured by light-field cameras may be referred to as the light-field, the plenoptic function, or radiance. In computational photography, a light-field is a 4-D record of all light rays in 3-D. Radiance describes both spatial and angular information, and is defined as density of energy per unit of area per unit of stereo angle (in radians). A light-field camera captures radiance; therefore, light-field images originally taken out-of-focus may be refocused, noise may be reduced, viewpoints may be changed, and other light-field effects may be achieved.
Light-fields may be captured with a conventional camera. In one conventional method, M×N images of a scene are captured from different positions with a conventional camera. If, for example, 8×8 images are captured from 64 different positions, 64 images are produced. The pixel from each position (i,j) in each image are taken and placed into blocks, to generate 64 blocks.
FIG. 1A illustrates an exemplary prior art light-field camera, or camera array, which employs an array of two or more objective lenses 110. Each objective lens focuses on a particular region of photosensor 108, or alternatively on a separate photosensor 108. This light-field camera 100 may be viewed as a combination of two or more conventional cameras that each simultaneously records an image of a subject on a particular region of photosensor 108 or alternatively on a particular photosensor 108. The captured images may then be combined to form one image.
FIG. 1B illustrates an exemplary prior art plenoptic camera, another type of light-field camera, that employs a single objective lens and a microlens or lenslet array 106 that includes, for example, about 100,000 lenslets. Lenslet array 106 is typically placed a small distance ( ̃0.5 mm) from a photosensor 108, e.g. a charge-coupled device (CCD). The raw image captured with a plenoptic camera 102 is made up of an array of small images, typically circular, of the main camera lens 108. These small images may be referred to as microimages. The lenslet array 106 enables the plenoptic camera 102 to capture the light-field, i.e. to record not only image intensity, but also the distribution of intensity in different directions at each point. Each lenslet splits a beam coming to it from the main lens 104 into rays coming from different “pinhole” locations on the aperture of the main lens 108. Each of these rays is recorded as a pixel on photosensor 108, and the pixels under each lenslet collectively form an n-pixel image. This n-pixel area under each lenslet may be referred to as a macropixel, and the camera 102 generates a microimage at each macropixel. The plenoptic photograph captured by a camera 102 with, for example, 100,000 lenslets will contain 100,000 macropixels, and thus generate 100,000 microimages of a subject. Each macropixel contains different angular samples of the light rays coming to a given microlens. Each macropixel contributes to only one pixel in the different angular views of the scene. As a result, each angular view contains 100,000 pixels.
Another type of integral or light-field camera is similar to the plenoptic camera of FIG. 1B, except that an array of pinholes are used between the main lens and the photosensor instead of an array of lenslets.
Yet another type of integral or light-field camera is similar to the plenoptic camera of FIG. 1B, except that a non-refractive mask is used between the main lens and the photosensor instead of an array of lenslets. The mask is a non-refractive element, and attenuates the incoming rays but does not bend them. The captured image is the convolution of the incoming light field with the mask light field. This camera design captures the 4-D light field directly in the Fourier domain. Thus, a 2-D sensor pixel represents a coded linear combination of several rays. The linear combination can be decoded by software to obtain the 4-D light field.
Frequency Domain Processing of Light-Field Images
Light-field capture may be analyzed in spatio-angular representation. A light-field camera may sample the optical signal within a single photograph by multiplexing the 4-D radiance onto the physical 2-D surface of the sensor. In addition to sampling the light field spatially, methods have been developed for multiplexing the radiance in the frequency domain by optically mixing different spatial and angular frequency components. The idea of multiplexing the 4-D light field onto a 2-D sensor is extremely note-worthy and powerful. More specifically, the idea that the use of a simple mask can turn a conventional camera into a light field camera is very practical. This approach requires processing the captured light-fields in the frequency domain. The general concept of processing data in the frequency domain with traditional tools like Fourier transform is fairly straightforward, but there are implementation details and caveats that need to be considered in order to obtain visually pleasing results.
Various techniques for analyzing radiance in the frequency domain have been developed. Along with techniques for analyzing radiance in the frequency domain, a number of good results have been derived, such as application of Poisson summation formula to depth representation, general light fields and displays, optical transforms, Fourier slice theorem applied to refocusing, and others. However, one area that has not achieved full attention is the quality of the resultant images. This is because achieving high quality results by doing frequency domain processing as a primary method of decoding the 4-D light field, is nontrivial. Although the overall method of decoding the light field is relatively straightforward, there are various limitations in getting things just right. Thus, results are prone to artifacts due to these limitations that are inherent to the frequency-domain processing of images and in image processing techniques in general.
SUMMARY
Various embodiments of a method and apparatus for managing artifacts in frequency domain processing of light-field images are described. Various limitations inherent to conventional methods of frequency-domain processing of light-field images may result in artifacts in the resulting output images. Various methods for the reduction or removal of the artifacts are described that include methods applied during frequency domain processing of light-field images and methods applied as or during post-processing of the results of frequency domain processing. The described methods for the removal of artifacts may be implemented in software as or in one or more light-field frequency domain processing modules. The module(s) may, for example, be implemented in a light-field image processing application or library. In some embodiments, one or more of the described methods for artifact removal may be implemented as components in or sub-modules of a light-field frequency domain processing module.
The process of light-field capture through mask-based or lens-based light-field cameras multiplexes the angular information within a 4-D light field visible in the 2-D Fourier transform of the acquired image. This 2-D Fourier transform reveals slices which correspond to each angular view. These modulated slices may be demodulated (or simply cropped) and rearranged into a 2-D stack. The center of a slice is generally determined as the pixel with the highest absolute value or frequency. However, a center determined as the pixel with the highest absolute value may not be the true center, i.e. the DC point, of the underlying analog signal. Incorrect centers of the slices may result in wave artifacts in the resultant angular views. In one embodiment, an oversampling method may be implemented to determine the correct analog centers of the slices. In one embodiment, a phase multiplication method may be implemented to determine the correct analog centers of the slices. Angular views may be generated from the slices using the determined correct centers. Determining and using the correct centers (the DC points) of the underlying analog signals may help to reduce or eliminate wave artifacts from the resultant angular views.
In one embodiment, two or more low-energy slices may be excluded from the demultiplexing process. The excluded slices may be the slices closest to the central slice or DC point, and are the slices containing the lowest energy. Excluding the low-energy slices may help to reduce or eliminate wave artifacts from the resultant angular views.
Angular views may be generated from the slices. The angular views may still include some artifacts. In one embodiment, a cosmetic correction of the intermediate angular views may be performed to generate corrected output angular views. After frequency domain processing is performed to generate a set of angular views, a base angular view may be selected. Cosmetic correction of the other angular views, or of at least a portion of the other angular views, may then be performed to adjust the other angular views according to the base angular view. In one embodiment, the pixel values of an angular view may be adjusted according to a Gaussian filter comparison of the angular view to the base angular view to generate a corrected angular view.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1A illustrates an exemplary prior art light-field camera, or camera array.
FIG. 1B illustrates an exemplary prior art plenoptic camera.
FIG. 2 illustrates an exemplary method for the frequency domain processing of a captured light-field image.
FIGS. 3A and 4A illustrate exemplary angular views in which intensity waves appear as an artifact of conventional frequency domain processing of respective raw light-field images.
FIGS. 3B and 4B illustrate the angular views of FIGS. 3B and 4B with artifacts removed according to embodiments of the methods described herein.
FIG. 5 illustrates a light-field camera in which the F-number of the optical systems consisting of the mask elements and the photosensor is equal to the F-number of the optical system consisting of the photosensor and the main lens.
FIG. 6 illustrates a light-field camera in which the F-number of the optical systems consisting of the mask elements and the photosensor is not equal to the F-number of the optical system consisting of the photosensor and the main lens.
FIGS. 7A through 7C graphically illustrate an exaggerated effect of a transfer function on the signal in the frequency domain.
FIG. 8A shows an angular view without applying a phase multiplication method to detect the correct centers of slices.
FIG. 8B shows the angular view of FIG. 8A for which an embodiment of a phase multiplication method was applied to detect the correct centers of slices.
FIGS. 9A and 9B illustrate an angular view before and after application of a cosmetic correction method, according to one embodiment.
FIG. 10A illustrates a light-field frequency domain processing module, according to one embodiment.
FIG. 10B illustrates a light-field frequency domain processing module in more detail, according to one embodiment.
FIG. 11 is a flowchart of a method for removing artifacts in frequency domain processing of light-field images according to one embodiment.
FIG. 12 is a flowchart of a method for removing artifacts in frequency domain processing of light-field images by finding the correct center of slices according to one embodiment.
FIG. 13 is a flowchart of a method for removing artifacts in frequency domain processing of light-field images by excluding low-energy slices according to one embodiment.
FIG. 14 is a flowchart of a method for performing cosmetic correction of intermediate angular views to generate corrected angular views according to one embodiment.
FIG. 15 is a high-level flowchart of a method for removing artifacts in frequency domain processing of light-field images by finding the correct center of slices using phase multiplication according to one embodiment.
FIG. 16 is a flowchart of a method for finding the correct center of a slice in frequency domain processing of light-field images using phase multiplication according to one embodiment.
FIG. 17 illustrates an exemplary computer system that may be used in embodiments.



While the invention is described herein by way of example for several embodiments and illustrative drawings, those skilled in the art will recognize that the invention is not limited to the embodiments or drawings described. It should be understood, that the drawings and detailed description thereto are not intended to limit the invention to the particular form disclosed, but on the contrary, the intention is to cover all modifications, equivalents and alternatives falling within the spirit and scope of the present invention as defined by the appended claims. The headings used herein are for organizational purposes only and are not meant to be used to limit the scope of the description or the claims. As used throughout this application, the word “may” is used in a permissive sense (i.e., meaning having the potential to), rather than the mandatory sense (i.e., meaning must). Similarly, the words “include”, “including”, and “includes” mean including, but not limited to.
DETAILED DESCRIPTION OF EMBODIMENTS
Various embodiments of a method and apparatus for managing artifacts in frequency domain processing of light-field images are described. Various limitations inherent to conventional methods of frequency-domain processing of light-field images may result in artifacts in the resulting output images. A discussion of the characteristics, sources and effects of artifacts that affect the quality of the results of frequency domain processing of light-field images is provided. Various methods for the reduction or removal of the artifacts are described. The methods include methods applied during frequency domain processing of light-field images and methods applied as or during post-processing of the results of frequency domain processing. The described methods for the removal of artifacts may be implemented in software as or in one or more light-field frequency domain processing modules. The module(s) may, for example, be implemented in a light-field image processing application or library. In some embodiments, one or more of the described methods for artifact removal may be implemented as components in or sub-modules of a light-field frequency domain processing module.
A process of light-field capture through mask-based or lens-based light-field cameras may multiplex the angular information within a 4-D light field visible in the 2-D Fourier transform of the acquired image. This 2-D Fourier transform reveals slices which correspond to each angular view. These modulated slices may be demodulated (or simply cropped) and rearranged into a 2-D stack. An inverse 2-D Fourier transform may then be performed on slices in the stack to generate angular views of the scene in the acquired image. Each angular view may represent a scene captured in the light-field image at a particular focal plane.
FIG. 2 illustrates an exemplary method for the frequency domain processing of a captured light-field image. An exemplary raw light-field image 200 is captured with a light-field camera, for example a plenoptic camera 102 as illustrated in FIG. 1B. Region 202 shows an expanded region of exemplary light-field image 200 including multiple microimages. Note that raw light-field images captured with other types of light-field cameras may look visually different than exemplary images 200 and 202. A two-dimensional (2-D) Fourier transform is performed on light-field image 200 to generate a frequency domain image 204. Two-dimensional (2-D) slice extraction is performed to generate multiple slices 206. Centers of the slices may be determined, and one or more stacks of slices may be generated. A 4-D inverse Fourier transform may be performed on the stacks of slices to generate results, in this example multiple images 208 focused at different depths. The multiple images 208 each represent a distinct angular view of the captured scene. The frequency domain processing may be performed on each color channel of an image separately (e.g., Red, Green, and Blue, for RGB images), and the color channels may then be aligned and merged after the frequency-based processing is complete.
Frequency Domain Processing Artifacts
The artifacts that result from frequency domain processing of light-field images are best visible as an intensity wave across the various angular views (e.g., images 208 of FIG. 2). These waves may appear as random low frequency waves that are present spatially within a single angular view, but that appear to travel across multiple angular views. Despite the fact that the waves appear in individual angular views, the waves may partially or completely disappear when the views are mixed for an application such as refocusing. FIGS. 3A and 4A illustrate exemplary angular views in which intensity waves appear as an artifact of conventional frequency domain processing of respective raw light-field images.
In addition to the wave-like artifact, a set of angular views generated by conventional frequency domain processing of a light-field image may appear to be warped in a peculiar way such that when the angular views are viewed in a movie, a wind-like movement or artifact may be apparent.
Another type of artifact that may appear is chromatic aberration, which may result from an incorrect alignment of the separately-processed color channels (e.g., Red, Green, and Blue, for RGB images) after frequency-based processing.
Sources of Wave Artifacts
The process of light-field capture through mask-based or lens-based light-field cameras multiplexes the angular information within a 4-D light field visible in the 2-D Fourier transform of the acquired image. as the modulated slices (e.g., slices 206 of FIG. 2) may be demodulated (or simply cropped) and rearranged into a 2-D stack. As a result, a crucial step involved in processing is to detect the correct centers of these slices. The center of a slice is generally determined as the pixel with the highest absolute value. However, a center determined as the pixel with the highest absolute value may not be at the true center or DC point of the underlying analog signal. The true center or DC point of the analog signal may be assumed to be within one pixel of the pixel with the highest absolute value. Assuming the center of the slice is determined to be (x, y), the actual location of the center (the center or DC point of the underlying signal) may be expressed by (x+δx, y+δy) where |δx, δy|<1.
Incorrect detection of the centers of these slices directly corresponds to modulation of these cropped slices by a very low frequency. Since a discrete Fourier transform (DFT) samples the underlying discrete-time Fourier transform (DTFT) of the image, detection of the discrete location (a pixel) of the peak (the center) can be correct only within a 1-pixel neighborhood. In other words, the DC point of the underlying analog signal, the correct center of the slice, is within a 1-pixel neighborhood of the detected pixel location. Also, there may be cases where the detected center can be completely incorrect, due to the fact that the original signal contains large low-frequency side-lobes. Incorrect centers for slices is responsible for one of the most prominent artifacts that arise in the results obtained by light field cameras—the traveling intensity waves, as illustrated in FIGS. 3A and 4A.
Sources of Wind Artifacts
Another type of artifact may be seen as a wind-like movement across a set of angular views generated by conventional frequency domain processing of a light-field image. One possible explanation for this wind-like artifact may be related to the circular aperture of cameras and the F-numbers of the main lens and the elements in the lens array or mask situated between the main lens and the photosensor. Referring to FIG. 5, in a light-field camera, each of the elements 310 (lenslet, pinhole, etc.) in the mask 306 captures a different part or view (aperture slices 314) of the main aperture or lens 304 on a different region 312 of the photosensor 308. The photosensor 308 thus captures light from different slices 314 of the main lens 304 as microimages at different locations or regions 312 on the photosensor 308.
FIG. 5 illustrates the case where the F-number of the optical systems consisting of the mask elements 310 and the photosensor 308 is equal to the F-number of the optical system consisting of the photosensor 308 and the main lens 304. Each slice 314 of the aperture, because the lens is circular, captures a different amount of light; the central slice 314 captures the largest amount of light. As a result, each microimage captured on the photosensor 308 contains a different amount of light. When these microimages are put together to form one coherent angular view, the views contain different intensity. This may result in the described wind-like artifact.
The wind-like artifact may be worsened when considering the case of FIG. 6. In the light-field camera of FIG. 6, the F-number of the optical systems consisting of the mask elements 310 and the photosensor 308 is different than the F-number of the optical system consisting of the photosensor 308 and the main lens 304. As a result, the microimages in regions 312 overlap each other on the photosensor 308. This results in the edges of the microimages containing overlapping information from neighboring microimages. However, this spatially neighboring information is not well-neighboring from an angular information point of view. This means that there may be redundant information in the farthest angular views.
The above assumes a best-case scenario where the physical mask 306 is optically high quality and is placed equidistant from the photosensor 308. Distance of the mask 306 from the photosensor 308 directly effects amount of overlapping angular information. This may also be viewed as spatial aliasing for which the frequency counterpart can be understood in the angular dimension. If the mask 306 were physically spaced unevenly in relation to the photosensor 308, this would correspond to the diffusion of angular information across microimages. This diffusion can be confined to different regions of the image as it is simply controlled by the distance of the mask 306 from the photosensor 308.
In addition, when all of the angular views are extracted, each pixel does not necessarily correspond to the same view because of the above-described phenomena. This observation leads to the conclusion that the true angular view is not a simple frame in the space-time volume of angular views, but rather a surface within the volume.
Another source of artifacts may due to the optical transfer function of the lens or the pre-mask optical system. FIGS. 7A through 7C graphically illustrate an exaggerated effect of this transfer function on the signal in the frequency domain. FIG. 7A illustrates the Fourier transform of the original light field before it passes through the camera optical system. FIG. 7B illustrates the hypothetical transfer function of the camera's optical system. FIG. 7C illustrates the resulting light field as captured by the photosensor. It is evident in FIG. 7C that the angular slices are not symmetric, causing the overall result of the entire process to be complex instead of real. The effect of this attenuation function is more difficult to correct when approaching high frequencies, as it becomes harder to invert. This problem is similar to a traditional signal processing problem of amplitude demodulation in the presence of a memory-less channel.
Methods of Artifact Removal
This section presents various methods of artifact removal that may be implemented in various embodiments of a light-field frequency domain processing module. The methods include methods applied during frequency domain processing of light-field images and methods applied as or during post-processing of the results of frequency domain processing.
The first two described methods of artifact removal in frequency domain processing of light-field images are directed at finding the correct center of the Fourier slices (e.g., slices 206 in FIG. 2). As described above, determining the correct center of the Fourier slices (e.g., slices 206 in FIG. 2) is important in removing artifacts, specifically wave artifacts. In an extremely well-calibrated light-field camera, the correct location of the centers may be predetermined. However, for most light-field cameras, this is difficult or impossible to achieve. Another way to view this problem is to find the location of the carrier frequencies (of the mask or microlens array) in the observed image. Referring to traditional signal processing, putting a mask in front of a photosensor corresponds to modulating the signal with a carrier frequency. Modulation simply replicates the FFT of the signal at the location of the carrier frequencies in the FFT. Therefore, the centers of the slices are the same as the locations of the carrier frequencies contained in the optical mask or the lens array.
Oversampling
One method to more correctly detect the center of a Fourier slice (e.g., a slice 206 in FIG. 2) that may be implemented in embodiments of a light-field frequency domain processing module is to oversample the Discrete Fourier transform (DFT) of the light-field image, for example by two or three times, extract the slices from the oversampled DFT of the image, and find the centers of the slices. In one embodiment, oversampling is performed by padding the acquired light-field image with the correct amount of zeros according to the oversampling rate and taking the Fourier transform of the padded image. Oversampling the DFT of the light-field image may provide better clarity as to where the underlying peak or center (the DC point) of the analog signal actually lies.
However, due to practical limitations of memory and other resources on various operating systems and computing devices, this may be difficult to achieve with straightforward FFT implementations.
Phase Multiplication
An alternative method to oversampling for detecting the exact centers of the Fourier slices (e.g., slices 206 of FIG. 2) that may be implemented in embodiments of a light-field frequency domain processing module employs phase multiplication as described below. For most practical examples, the neighborhood of certainty for the center of the slices lies within an area of 1×1 pixel, and subpixel-accurate location of the center of the slices is being sought. Assume the center is determined to be (x, y); however, the actual location of the center is (x+δx, y+δy) where |δx, δy|<1. Instead of shifting the entire slice by (δx, δy), embodiments may use the following property of Fourier-transform:

F{x(x+δt)}=F( ω)e−iδt 

A slice Si may be multiplied by the phase factor:




  e
  
    -
    
      i
      ⁡
      
        (
        
          
            
              
                2
                ⁢
                
                  
                
                ⁢
                πδ
                ⁢
                
                  
                
                ⁢
                x
              
              N
            
            ⁢
            x
          
          +
          
            
              
                2
                ⁢
                
                  
                
                ⁢
                πδ
                ⁢
                
                  
                
                ⁢
                x
              
              M
            
            ⁢
            y
          
        
        )
      
    
  




After multiplying by the above phase factor, the inverse Fourier transform may be computed to get Si,mid. Si,mid is then rearranged in 2-D, and another 2-D-inverse Fourier-transform is performed in order to obtain the angular views (e.g., views 208 of FIG. 2). As a result, the overall process, e.g. as illustrated in FIG. 2:

4-DIFT(arrangein4-D(R(ωx,ωθ)))

is replaced with:

2-DIFT(arrangein2-D(correct phase(arrangein2-D(R(ωx,ωθ)))).

Various automated or manual methods may be used in various embodiments to detect the motion (δx, δy). In on embodiment, a “brute force” or manual method may be used in order to detect the shift in which a user visually determines a best set of output angular views. FIG. 8A shows a resultant angular view without applying the above phase multiplication method, and FIG. 8B shows a resultant angular view from the same light-field image to which an embodiment of the phase multiplication method was applied.
Skipping or Excluding the Lowest Slices
One embodiment of a light-field frequency domain processing module may exclude the “lowest slices” from the demultiplexing process. The lowest slices are the slices with the lowest energy, and are the slices closest to the center slice. The center slice contains the DC point; the lowest slices are the slices closest to the DC point that are not the center slice. Slices higher (or further) from the center of the transform contain higher and higher frequencies. Slices closer to the center, or the DC point, contain more low energy or low frequencies.
Observation and experiments with various light-field camera designs indicates that there are cases when the effect of low intensity waves is severe. That is, major regions of all angular views are nearly or completely black. This indicates that the intensity wave is fairly stationary across all the angular views. This may correspond to a very low-frequency component that corrupts the results in the angular dimension. Such corruption appears to be introduced through the carrier (or modulator). Excluding the lowest modulated slices, the slices closest to the central slice or DC point, which are the slices containing the lowest energy, from the entire demultiplexing process may result in the wave artifacts almost completely disappearing. FIG. 3B illustrates an exemplary angular view for which the lowest slices were excluded during processing.
Cosmetic Correction
The above methods may be applied by embodiments of a light-field frequency domain processing module during frequency domain processing of light-field images. However, these methods may remove some, or most, artifacts, but may not remove all artifacts; even with the application of the above-described methods, vignetting-type artifacts, varying color-shifts, or other wavy artifacts in image intensity may still be visible in some output angular views. Thus, as an alternative to, or in addition to, the above-described methods, some embodiments may provide a method for artifact removal during a post-processing stage. If the previously described methods are not applied during frequency domain processing, or if one or more of the methods are applied but the resulting angular views still display some artifacts, a post-processing cosmetic correction technique may be applied to remove some or all remaining artifacts to improve the look of the resultant angular view(s) generated by the frequency domain processing. A supposition is made that any artifacts are present in the low-frequency content of the image. In accordance with this supposition, the cosmetic corrections are targeted at the low-end of the image spectrum.
The cosmetic correction method is additive in nature; however, the cosmetic correction method may greatly increase the visual quality of the extract angular views. As previously mentioned, frequency-domain processing of a light-field image generates multiple angular views, e.g. views 206 of FIG. 2. In the cosmetic correction method, one of the angular views is selected as a base view; the base view is the best view according to color/intensity as determined by a human observer or operator or by an automated technique, and then one or more of the other angular views are adjusted to match (color/intensity) the selected, or base, view. The cosmetic correction method is essentially a color/intensity matching process. A base view is selected, and the pixel values in one or more of the other angular views are adjusted according to the pixel values in the base view to generate cosmetically corrected angular views from the one or more of the other angular views. Referring to the base view as Sbase, one or more other angular views are adjusted according to the following formula:

Si,new=Si+Filter(Sbase−Si)

Here, Filter( ) is a low-pass filter. For the sake of speed, this may be accomplished in one embodiment by:

Filter(x)=gaussian(σ)*downsample(x)

The exemplary angular view illustrated in FIG. 4B was obtained using the above-described cosmetic correction method. In addition, FIG. 9A shows an original angular view, and FIG. 9B shows the same view after application of the cosmetic correction method. Note that, while the results are chromatically correct, darker regions may be replaced with blurry content by the cosmetic correction method.
Implementations
Various methods for the reduction or removal of the artifacts have been described. The methods include methods applied during frequency domain processing of light-field images and methods applied as or during post-processing of the results of frequency domain processing. The described methods for the removal of artifacts may be implemented in software as or in one or more light-field frequency domain processing modules. The module(s) may, for example, be implemented in a light-field image processing application or library. In some embodiments, one or more of the above-described methods for artifact removal may be implemented as component(s) in a light-field frequency domain processing module. An exemplary hardware system on which various embodiments may be instantiated is illustrated in FIG. 17.
FIG. 10A illustrates a light-field frequency domain processing module, according to one embodiment. Light-field frequency domain processing module 400 receives an input light-field image 410. Light-field frequency domain processing module 400 performs frequency domain processing on the input image 410 to generate one or more output angular views 440. During frequency domain processing, one or more of the above-described methods for removing artifacts may be applied. After frequency domain processing, the above-described method for cosmetic correction may be applied. Output angular views 440 may, for example, be stored to a storage medium 450, such as system memory, a disk drive, DVD, CD, etc.
Light-field frequency domain processing module 400 may provide a user interface that provides one or more textual and/or graphical user interface elements, modes or techniques via which a user may view or control various aspects of frequency domain processing including, but not limited to, implementations of the artifact removal methods as described herein. For example, the user interface may include user interface elements that allow a user to select input and output files, to initiate the cosmetic correction method, to select a base image Sbase for the cosmetic correction process, to specify a method to be used to detect the correct center of slices, and so on.
FIG. 10B illustrates a light-field frequency domain processing module in more detail, according to one embodiment. This Figure shows light-field frequency domain processing module 400 as including frequency domain processing 420 and post-processing 430. During frequency domain processing 420, a slice center detection component 422 may implement one or both of the above-described methods for determining the correct center of the slices extracted from the 2-D Fourier image. If both methods are implemented, the user interface may provide a user interface mechanism whereby the user may select one, both, or even neither of the methods to be performed on a particular input image 410. Also during frequency domain processing 420, a slice exclusion component 424 may implement the slice exclusion method as described above. In the slice exclusion method, two or more slices with the lowest energy (e.g., slices nearest to the center slice) may be determined and excluded from further processing to generate output angular views 440. Note that, while FIG. 10B shows slice exclusion component 424 as logically after slice center detection component 422, slice center exclusion may be performed before or after slice center detection in different embodiments.
Frequency domain processing 420 may generate intermediate angular views 440A. During post-processing 430, a cosmetic correction component 432 may implement the cosmetic correction method described above on intermediate angular views 440A to generate output angular views 440B with artifacts removed. In one embodiment, cosmetic correction may be optional; thus, the user interface may provide one or more mechanisms whereby the user may select to perform or not perform cosmetic correction. In one embodiment, the user may be asked, via the user interface, to select and specify one of intermediate angular views 440A to be used as the base view (Sbase) in the cosmetic correction process. In one embodiment, the user may be asked, via the user interface, to select and specify, or to exclude, particular ones of the other angular views to which the cosmetic correction process is to be applied. Note that the base view itself may not be modified in the cosmetic correction process.
FIG. 11 is a flowchart of a method for removing artifacts in frequency domain processing of light-field images according to one embodiment. The method may be implemented in a light-field frequency domain processing module, such as module 400 in FIGS. 10A and 10B. Various embodiments may implement some or all of the elements illustrated in FIG. 11. In some embodiments, various ones of the elements may be optionally performed or skipped according to user preference. In some embodiments, the elements may be performed in a different order than shown.
As indicated at 500, a light-field image may be obtained. As indicated at 502, a 2-D Fourier transform of the light-field image may be generated. As indicated at 504, multiple 2-D slices may be extracted from the 2-D Fourier transform.
Incorrect centers of the slices may result in wave artifacts in the resultant angular views. In one embodiment, as indicated at 506, the correct centers of the 2-D slices may be determined. In one embodiment, an oversampling method may be implemented to determine the correct centers of the slices, as previously described and as illustrated in FIG. 12. In the oversampling method, the 2-D Fourier transform of the light-field image may be oversampled prior to slice extraction, e.g. prior to element 504. In one embodiment, oversampling is done by padding the acquired light-field image with the correct amount of zeros according to the oversampling rate and taking the Fourier transform of the padded image. In one embodiment, a phase multiplication method may be implemented to determine the correct centers of the slices after the slices are extracted, as previously described and as illustrated in FIGS. 15 and 16.
As indicated at 508, two or more low-energy slices may be excluded from the demultiplexing process. The excluded slices may be the slices closest to the central slice or DC point, and are the slices containing the lowest energy. Excluding the low-energy slices may help to reduce or eliminate wave artifacts from the resultant angular views. FIG. 13 further illustrates a method of excluding low-energy slices.
As indicated at 510, intermediate angular views may be generated from the slices. In one embodiment, to generate an angular view, multiple ones of the slices may be stacked according to the determined centers of the slices, and an inverse 4-D Fourier transform may be performed on the stack to generate the angular view.
The intermediate angular views may still include some artifacts. As indicated at 512, a cosmetic correction of the intermediate angular views may be performed to generate corrected output angular views. In one embodiment, the user may be given the option to perform the cosmetic correction, if desired. FIG. 14 further illustrates a method of performing cosmetic correction.
FIG. 12 is a flowchart of a method for removing artifacts in frequency domain processing of light-field images by finding the correct center of slices according to one embodiment. As indicated at 520, the 2-D Fourier transform of a light-field image may be oversampled. As indicated at 522, 2-D slices may be extracted from the oversampled 2-D Fourier transform of the light-field image. As indicated at 524, the correct centers of the 2-D slices may be determined. Angular views may then be generated from the slices using the correct centers, as indicated at 526.
FIG. 13 is a flowchart of a method for removing artifacts in frequency domain processing of light-field images by excluding low-energy slices according to one embodiment. As indicated at 530, 2-D slices may be extracted from a 2-D Fourier transform of a light-field image. As indicated at 532, low-energy slices may be determined. In one embodiment, the low-energy slices are determined as the slices closest to the central slice or DC point. As indicated at 534, the low-energy slices are excluded from the demultiplexing process. Angular views may then be generated from the slices using the remaining, as indicated at 536.
FIG. 14 is a flowchart of a method for performing cosmetic correction of intermediate angular views to generate corrected angular views according to one embodiment. After frequency domain processing is performed to generate a set of angular views, a base angular view may be selected, as indicated at 540. In one embodiment, a user may select a visually preferred view as the base angular view. Cosmetic correction of the other angular views, or of at least a portion of the other angular views, may then be performed to adjust the other angular views according to the base angular view. As indicated at 542, another, or next, angular view may be obtained. As indicated at 544, the pixel values of the other angular view may be adjusted according to a Gaussian filter comparison of the angular view to the base angular view to generate a corrected angular view. Referring to the base view as Sbase, the other angular view may be adjusted according to the following formula:

Si,new=Si+Filter(Sbase−Si)

Here, Filter( ) is a low-pass filter. For the sake of speed, this may be accomplished in one embodiment by:

Filter(x)=gaussian(σ)*downsample(x)

At 546, if there are more angular views to be corrected, the method returns to 542. Otherwise, the cosmetic correction process is done.
FIG. 15 is a high-level flowchart of a method for removing artifacts in frequency domain processing of light-field images by finding the correct center of slices using phase multiplication according to one embodiment. As indicated at 550, 2-D slices may be extracted from a 2-D Fourier transform of a light-field image. As indicated at 552, a phase multiplication method may be performed to determine the correct centers of the 2-D slices. Assume the center is determined to be (x, y); however, the actual location of the center is (x+δx, y+δy) where |δx, δy|<1. Instead of shifting the entire slice by (δx, δy), a slice Si may be multiplied by the phase factor:




  e
  
    -
    
      i
      ⁡
      
        (
        
          
            
              
                2
                ⁢
                
                  
                
                ⁢
                πδ
                ⁢
                
                  
                
                ⁢
                x
              
              N
            
            ⁢
            x
          
          +
          
            
              
                2
                ⁢
                
                  
                
                ⁢
                πδ
                ⁢
                
                  
                
                ⁢
                x
              
              M
            
            ⁢
            y
          
        
        )
      
    
  




After multiplying by the above phase factor, the inverse Fourier transform may be computed to get Si,mid. Si,mid is then rearranged in 2-D, and another 2-D-inverse Fourier-transform is performed in order to obtain the angular views (e.g., views 208 of FIG. 2). As a result, the overall process may be represented by:

2-DIFT(arrangein2-D(correct phase(arrangein2-D(R(ωx,ωθ)))).

Angular views may then be generated from the slices using the correct centers, as indicated at 554.
FIG. 16 is a flowchart of a method for finding the correct center of a slice in frequency domain processing of light-field images using phase multiplication according to one embodiment. 2-D slices may be extracted from a 2-D Fourier transform of a light-field image. The true center of the analog signal is assumed to be within one pixel of the initial center, generally determined as the pixel with the highest absolute value. Thus, one method to determine the true center is to iteratively perform phase multiplication at different intervals, or shifts, within one pixel of the center, generate results, and determine or select the best results. As indicated at 560, the shift is initialized. The shift is a value of less than one pixel, since the true center of the analog signal is assumed to be within one pixel of the initial center, generally determined as the pixel with the highest absolute value. As indicated at 562, a slice Si is multiplied by the phase factor:




  e
  
    -
    
      i
      ⁡
      
        (
        
          
            
              
                2
                ⁢
                
                  
                
                ⁢
                πδ
                ⁢
                
                  
                
                ⁢
                x
              
              N
            
            ⁢
            x
          
          +
          
            
              
                2
                ⁢
                
                  
                
                ⁢
                πδ
                ⁢
                
                  
                
                ⁢
                x
              
              M
            
            ⁢
            y
          
        
        )
      
    
  




After multiplying by the above phase factor, the inverse Fourier transform may be computed to get Si,mid. Si,mid is then rearranged in 2-D, and another 2-D-inverse Fourier-transform is performed in order to obtain the angular view, as indicated at 568.
At 570, if there are more shifts to be tried, the next shift is set at 568, and the method returns to 562. Otherwise, all iterations have been performed, and a best or preferred result may be selected, as indicated at 572.
Exemplary System
Various components of embodiments of a method and apparatus for managing artifacts in frequency domain processing of light-field images may be executed on one or more computer systems, which may interact with various other devices. One such computer system is illustrated by FIG. 17. In the illustrated embodiment, computer system 700 includes one or more processors 710 coupled to a system memory 720 via an input/output (I/O) interface 730. Computer system 700 further includes a network interface 740 coupled to I/O interface 730, and one or more input/output devices 750, such as cursor control device 760, keyboard 770, audio device 790, and display(s) 780. In some embodiments, it is contemplated that embodiments may be implemented using a single instance of computer system 700, while in other embodiments multiple such systems, or multiple nodes making up computer system 700, may be configured to host different portions or instances of embodiments. For example, in one embodiment some elements may be implemented via one or more nodes of computer system 700 that are distinct from those nodes implementing other elements.
In various embodiments, computer system 700 may be a uniprocessor system including one processor 710, or a multiprocessor system including several processors 710 (e.g., two, four, eight, or another suitable number). Processors 710 may be any suitable processor capable of executing instructions. For example, in various embodiments, processors 710 may be general-purpose or embedded processors implementing any of a variety of instruction set architectures (ISAs), such as the x86, PowerPC, SPARC, or MIPS ISAs, or any other suitable ISA. In multiprocessor systems, each of processors 710 may commonly, but not necessarily, implement the same ISA.
System memory 720 may be configured to store program instructions and/or data accessible by processor 710. In various embodiments, system memory 720 may be implemented using any suitable memory technology, such as static random access memory (SRAM), synchronous dynamic RAM (SDRAM), nonvolatile/Flash-type memory, or any other type of memory. In the illustrated embodiment, program instructions and data implementing desired functions, such as those described above for a method and apparatus for managing artifacts in frequency domain processing of light-field images, are shown stored within system memory 720 as program instructions 725 and data storage 735, respectively. In other embodiments, program instructions and/or data may be received, sent or stored upon different types of computer-accessible media or on similar media separate from system memory 720 or computer system 700. Generally speaking, a computer-accessible medium may include storage media or memory media such as magnetic or optical media, e.g., disk or CD/DVD-ROM coupled to computer system 700 via I/O interface 730. Program instructions and data stored via a computer-accessible medium may be transmitted by transmission media or signals such as electrical, electromagnetic, or digital signals, which may be conveyed via a communication medium such as a network and/or a wireless link, such as may be implemented via network interface 740.
In one embodiment, I/O interface 730 may be configured to coordinate I/O traffic between processor 710, system memory 720, and any peripheral devices in the device, including network interface 740 or other peripheral interfaces, such as input/output devices 750. In some embodiments, I/O interface 730 may perform any necessary protocol, timing or other data transformations to convert data signals from one component (e.g., system memory 720) into a format suitable for use by another component (e.g., processor 710). In some embodiments, I/O interface 730 may include support for devices attached through various types of peripheral buses, such as a variant of the Peripheral Component Interconnect (PCI) bus standard or the Universal Serial Bus (USB) standard, for example. In some embodiments, the function of I/O interface 730 may be split into two or more separate components, such as a north bridge and a south bridge, for example. In addition, in some embodiments some or all of the functionality of I/O interface 730, such as an interface to system memory 720, may be incorporated directly into processor 710.
Network interface 740 may be configured to allow data to be exchanged between computer system 700 and other devices attached to a network, such as other computer systems, or between nodes of computer system 700. In various embodiments, network interface 740 may support communication via wired or wireless general data networks, such as any suitable type of Ethernet network, for example; via telecommunications/telephony networks such as analog voice networks or digital fiber communications networks; via storage area networks such as Fibre Channel SANs, or via any other suitable type of network and/or protocol.
Input/output devices 750 may, in some embodiments, include one or more display terminals, keyboards, keypads, touchpads, scanning devices, voice or optical recognition devices, or any other devices suitable for entering or retrieving data by one or more computer system 700. Multiple input/output devices 750 may be present in computer system 700 or may be distributed on various nodes of computer system 700. In some embodiments, similar input/output devices may be separate from computer system 700 and may interact with one or more nodes of computer system 700 through a wired or wireless connection, such as over network interface 740.
As shown in FIG. 17, memory 720 may include program instructions 725, configured to implement embodiments of a method and apparatus for managing artifacts in frequency domain processing of light-field images as described herein, and data storage 735, comprising various data accessible by program instructions 725. In one embodiment, program instructions 725 may include software elements of a method and apparatus for managing artifacts in frequency domain processing of light-field images as illustrated in the above Figures. Data storage 735 may include data that may be used in embodiments. In other embodiments, other or different software elements and data may be included.
Those skilled in the art will appreciate that computer system 700 is merely illustrative and is not intended to limit the scope of a method and apparatus for managing artifacts in frequency domain processing of light-field images as described herein. In particular, the computer system and devices may include any combination of hardware or software that can perform the indicated functions, including computers, network devices, internet appliances, PDAs, wireless phones, pagers, etc. Computer system 700 may also be connected to other devices that are not illustrated, or instead may operate as a stand-alone system. In addition, the functionality provided by the illustrated components may in some embodiments be combined in fewer components or distributed in additional components. Similarly, in some embodiments, the functionality of some of the illustrated components may not be provided and/or other additional functionality may be available.
Those skilled in the art will also appreciate that, while various items are illustrated as being stored in memory or on storage while being used, these items or portions of them may be transferred between memory and other storage devices for purposes of memory management and data integrity. Alternatively, in other embodiments some or all of the software components may execute in memory on another device and communicate with the illustrated computer system via inter-computer communication. Some or all of the system components or data structures may also be stored (e.g., as instructions or structured data) on a computer-accessible medium or a portable article to be read by an appropriate drive, various examples of which are described above. In some embodiments, instructions stored on a computer-accessible medium separate from computer system 700 may be transmitted to computer system 700 via transmission media or signals such as electrical, electromagnetic, or digital signals, conveyed via a communication medium such as a network and/or a wireless link. Various embodiments may further include receiving, sending or storing instructions and/or data implemented in accordance with the foregoing description upon a computer-accessible medium. Accordingly, the present invention may be practiced with other computer system configurations.
CONCLUSION
Various embodiments may further include receiving, sending or storing instructions and/or data implemented in accordance with the foregoing description upon a computer-accessible medium. Generally speaking, a computer-accessible medium may include storage media or memory media such as magnetic or optical media, e.g., disk or DVD/CD-ROM, volatile or non-volatile media such as RAM (e.g. SDRAM, DDR, RDRAM, SRAM, etc.), ROM, etc. As well as transmission media or signals such as electrical, electromagnetic, or digital signals, conveyed via a communication medium such as network and/or a wireless link.
The various methods as illustrated in the Figures and described herein represent exemplary embodiments of methods. The methods may be implemented in software, hardware, or a combination thereof. The order of method may be changed, and various elements may be added, reordered, combined, omitted, modified, etc.
Various modifications and changes may be made as would be obvious to a person skilled in the art having the benefit of this disclosure. It is intended that the invention embrace all such modifications and changes and, accordingly, the above description to be regarded in an illustrative rather than a restrictive sense.

## Claims

1. A computer-implemented method, comprising:
processing a light-field image in the frequency domain to render multiple angular views of a scene, wherein, said processing comprises:
extracting multiple slices from a Fourier transform of the light-field image; and
rendering the multiple angular views from the slices in accordance with subpixel centers of the slices, wherein each subpixel center corresponds to the DC point of the analog signal underlying the respective slice, wherein rendering the multiple angular views from the slices in accordance with the subpixel centers of the slices prevents or lessens artifacts in the angular views;

wherein said processing further comprises applying a technique to determine the subpixel centers for the slices, wherein the technique is one of an oversampling technique or a phase multiplication technique;
wherein said oversampling technique comprises generating the Fourier transform as an oversampled Fourier transform of the light-field image and determining the subpixel centers from the slices extracted from the oversampled Fourier transform; and
wherein said phase multiplication technique comprises iteratively applying phase multiplication to each slice to locate the DC point of the analog signal underlying the respective slice.

2. The method as recited in claim 1, wherein the angular views each represent a scene captured in the light-field image at a particular focal plane.

3. The method as recited in claim 1, wherein said generating the Fourier transform as an oversampled Fourier transform comprises padding the light-field image according to an oversampling rate and generating the Fourier transform from the padded light-field image.

4. The method as recited in claim 1, wherein, in said phase multiplication technique, (x, y) are pixel coordinates of an initial center for a respective slice and (δx, δy) are subpixel shifts for the pixel coordinates determined during said iterative phase multiplication, where |δx, δy|<1, and where (x+δx, y+δy) are coordinates of the subpixel center corresponding to the DC point of the analog signal underlying the respective slice.

5. The method as recited in claim 1, wherein said applying phase multiplication to each slice comprises multiplying the respective slice by the phase factor




  e
  
    -
    
      i
      ⁡
      
        (
        
          
            
              
                2
                ⁢
                
                  
                
                ⁢
                πδ
                ⁢
                
                  
                
                ⁢
                x
              
              N
            
            ⁢
            x
          
          +
          
            
              
                2
                ⁢
                
                  
                
                ⁢
                πδ
                ⁢
                
                  
                
                ⁢
                x
              
              M
            
            ⁢
            y
          
        
        )
      
    
  




to locate the DC point of the analog signal underlying the respective slice, where (x, y) are pixel coordinates of an initial center for the respective slice and (δx, δy) are subpixel shifts for the pixel coordinates.

6. The method as recited in claim 1, wherein said rendering the multiple angular views from the slices comprises:
determining a subset of the slices as low-energy slices; and
excluding the low-energy slices when generating the angular views from the slices.

7. The method as recited in claim 6, wherein the low-energy slices are slices that are closest to a center slice, wherein the center slice includes the DC point of the Fourier transform of the light-field image.

8. The method as recited in claim 1, further comprising:
determining one of the angular views as a base view; and
adjusting the pixel values in one or more of the other angular views according to the pixel values in the base view to generate cosmetically corrected angular views from the one or more of the other angular views.

9. The method as recited in claim 1, wherein said rendering the multiple angular views from the slices in accordance with the subpixel centers of the slices comprises, for each angular view to be generated:
generating a stack including at least a portion of the slices, wherein the slices are aligned in the stack in accordance with the subpixel centers of the slices; and
performing an inverse Fourier transform on the stack to generate the angular view.

10. A system, comprising:
at least one processor; and
a memory comprising program instructions, wherein the program instructions are executable by the at least one processor to process a light-field image in the frequency domain to render multiple angular views of a scene, wherein, to process the light-field image, the program instructions are executable by the at least one processor to:
extract multiple slices from a Fourier transform of the light-field image; and
render the multiple angular views from the slices in accordance with subpixel centers of the slices, wherein each subpixel center corresponds to the DC point of the analog signal underlying the respective slice, wherein rendering the multiple angular views from the slices in accordance with the subpixel centers of the slices prevents or lessens artifacts in the angular views;
wherein, to process the light-field image, the program instructions are further executable by the at least one processor to apply a technique to determine the subpixel centers for the slices, wherein the technique is one of an oversampling technique or a phase multiplication technique;
wherein said oversampling technique generates the Fourier transform as an oversampled Fourier transform of the light-field image and determines the subpixel centers from the slices extracted from the oversampled Fourier transform; and
wherein said phase multiplication technique iteratively applies phase multiplication to each slice to locate the DC point of the analog signal underlying the respective slice.

11. The system as recited in claim 10, wherein the angular views each represent a scene captured in the light-field image at a particular focal plane.

12. The system as recited in claim 10, wherein, to generate the Fourier transform as an oversampled Fourier transform, the program instructions are executable by the at least one processor to pad the light-field image according to an oversampling rate and generate the Fourier transform from the padded light-field image.

13. The system as recited in claim 10, wherein, in said phase multiplication technique, (x, y) are pixel coordinates of an initial center for a respective slice and (δx, δy) are subpixel shifts for the pixel coordinates determined during said iterative phase multiplication, where |δx, δy|<1, and where (x+δx, y+δy) are coordinates of the subpixel center corresponding to the DC point of the analog signal underlying the respective slice.

14. The system as recited in claim 10, wherein, to apply phase multiplication to each slice, the program instructions are executable by the at least one processor to multiply the respective slice by the phase factor




  e
  
    -
    
      i
      ⁡
      
        (
        
          
            
              
                2
                ⁢
                
                  
                
                ⁢
                πδ
                ⁢
                
                  
                
                ⁢
                x
              
              N
            
            ⁢
            x
          
          +
          
            
              
                2
                ⁢
                
                  
                
                ⁢
                πδ
                ⁢
                
                  
                
                ⁢
                x
              
              M
            
            ⁢
            y
          
        
        )
      
    
  




to locate the DC point of the analog signal underlying the respective slice, where (x, y) are pixel coordinates of an initial center for the respective slice and (δx, δy) are subpixel shifts for the pixel coordinates.

15. The system as recited in claim 10, wherein, to render the multiple angular views from the slices, the program instructions are executable by the at least one processor to:
determine a subset of the slices as low-energy slices; and
exclude the low-energy slices when generating the angular views from the slices.

16. The system as recited in claim 15, wherein the low-energy slices are slices that are closest to a center slice, wherein the center slice includes the DC point of the Fourier transform of the light-field image.

17. The system as recited in claim 10, wherein the program instructions are executable by the at least one processor to:
determine one of the angular views as a base view; and
adjust the pixel values in one or more of the other angular views according to the pixel values in the base view to generate cosmetically corrected angular views from the one or more of the other angular views.

18. The system as recited in claim 10, wherein, to render the multiple angular views from the slices in accordance with the subpixel centers of the slices, the program instructions are executable by the at least one processor to, for each angular view to be generated:
generate a 2-dimensional (2-D) stack including at least a portion of the slices, wherein the slices are aligned in the stack in accordance with the subpixel centers of the slices; and
perform an inverse four-dimensional (4-D) Fourier transform on the stack to generate the angular view.

19. A non-transitory computer-readable storage medium storing program instructions, wherein the program instructions are computer-executable to implement:
processing a light-field image in the frequency domain to render multiple angular views of a scene, wherein, in said processing, the program instructions are computer-executable to implement:
extracting multiple slices from a Fourier transform of the light-field image; and
rendering the multiple angular views from the slices in accordance with subpixel centers of the slices, wherein each subpixel center corresponds to the DC point of the analog signal underlying the respective slice, wherein rendering the multiple angular views from the slices in accordance with the subpixel centers of the slices prevents or lessens artifacts in the angular views;

wherein, in said processing, the program instructions are further computer-executable to implement applying a technique to determine the subpixel centers for the slices, wherein the technique is one of an oversampling technique or a phase multiplication technique;
wherein said oversampling technique generates the Fourier transform as an oversampled Fourier transform of the light-field image and determines the subpixel centers from the slices extracted from the oversampled Fourier transform; and
wherein said phase multiplication technique iteratively applies phase multiplication to each slice to locate the DC point of the analog signal underlying the respective slice.

20. The non-transitory computer-readable storage medium as recited in claim 18, wherein the angular views each represent a scene captured in the light-field image at a particular focal plane.

21. The non-transitory computer-readable storage medium as recited in claim 18, wherein, in said generating the Fourier transform as an oversampled Fourier transform, the program instructions are computer-executable to implement padding the light-field image according to an oversampling rate and generating the Fourier transform from the padded light-field image.

22. The non-transitory computer-readable storage medium as recited in claim 18, wherein, in said phase multiplication technique, (x, y) are pixel coordinates of an initial center for a respective slice and (δx, δy) are subpixel shifts for the pixel coordinates determined during said iterative phase multiplication, where |δx, δy|<1, and where (x+δx, y+δy) are coordinates of the subpixel center corresponding to the DC point of the analog signal underlying the respective slice.

23. The non-transitory computer-readable storage medium as recited in claim 18, wherein, in said applying phase multiplication to each slice, the program instructions are computer-executable to implement multiplying the respective slice by the phase factor




  e
  
    -
    
      i
      ⁡
      
        (
        
          
            
              
                2
                ⁢
                
                  
                
                ⁢
                πδ
                ⁢
                
                  
                
                ⁢
                x
              
              N
            
            ⁢
            x
          
          +
          
            
              
                2
                ⁢
                
                  
                
                ⁢
                πδ
                ⁢
                
                  
                
                ⁢
                x
              
              M
            
            ⁢
            y
          
        
        )
      
    
  




to determine the DC point of the analog signal underlying the respective slice, where (x, y) are pixel coordinates of an initial center for the respective slice and (δx, δy) are subpixel shifts for the pixel coordinates.

24. The non-transitory computer-readable storage medium as recited in claim 18, wherein, in said rendering the multiple angular views from the slices, the program instructions are computer-executable to implement:
determining a subset of the slices as low-energy slices; and
excluding the low-energy slices when generating the angular views from the slices.

25. The non-transitory computer-readable storage medium as recited in claim 24, wherein the low-energy slices are slices that are closest to a center slice, wherein the center slice includes the DC point of the Fourier transform of the light-field image.

26. The non-transitory computer-readable storage medium as recited in claim 18, wherein the program instructions are computer-executable to implement:
determining one of the angular views as a base view;
adjusting the pixel values in one or more of the other angular views according to the pixel values in the base view to generate cosmetically corrected angular views from the one or more of the other angular views.

27. The non-transitory computer-readable storage medium as recited in claim 18, wherein, in said rendering the multiple angular views of the scene from the slices in accordance with the subpixel centers of the slices, the program instructions are computer-executable to implement, for each angular view to be generated:
generating a stack including at least a portion of the slices, wherein the slices are aligned in the stack in accordance with the determined subpixel centers of the slices; and
performing an inverse Fourier transform on the stack to generate the angular view.