An embodiment of the present invention is generally related to techniques of performing graphics processing to reduce the number of samples that need to be rendered. In particular, an embodiment of the present invention is directed to a technique to reduce the number of samples to be rendered in a tile of a current frame by reusing at least one pixel from the previous frame.
One aspect of the graphics pipeline 100 is that every region of an image is rendered at the same minimum resolution. In particular, in a conventional graphics pipeline a sampling rate (average number of samples per pixel) is typically at least one sample for every pixel of an image.
One aspect of the conventional graphics pipeline is that it is wasteful and requires more pixel shading operations then desired. In particular, there is no automation in the graphics pipeline to permit strategic choices to be automatically made to reduce a sampling rate below one sample per pixel (sub-sampling/de-sampling) in local regions of an image. In the context of mobile devices this means that the amount of power that is consumed in larger than desired.
An apparatus, system, and method are disclosed to selectively reuse pixel data from a previous frame by an advection technique. Pixel data and vertex coordinate data is buffered from a previous frame for reuse in a current frame. The speed of pixels in individual tiles is determined. If the pixels of a tile have a speed below a quasi-static speed threshold than a decision is made whether or not to reuse a fraction of the pixels from the previous frame.
One embodiment of method of performing graphics processing in a graphics pipeline via advection includes buffering pixel data and vertex coordinate data generated for a previous frame, n, into a buffer of a graphics pipeline processing a current frame n+1, where n is a positive integer. The speed of each pixel of the current frame is updated. A determination is made for a set of pixels in the current frame of corresponding positions in the previous frame. For at least one tile, color data from at least one pixel is selected from a pixel of the previous frame.
In one embodiment, an adaptive desampling (AD) sample generator stage 205 is provided to support adjusting a sampling pattern in local regions of an image, where the local region is a tile corresponding to a block of pixels (e.g., a 4×4 block of pixels, 16×16, or other size). Desampling is the reduction in the number of samples per tile that are sampled and rendered in the current frame. For example, desampling may include sampling and rendering on average less than one sample per pixel in a tile, and thus may also be described as sub-sampling. To maintain full image resolution, two different approaches may be used to obtain values of missing pixel data. A reconstruction and advection stage 210 supports two different options to reduce the number of pixels that need to be sampled and rendered in a tile while maintaining visual quality for a user. The reconstruction and advection stage 210 includes a reconstruction submodule 211 and an advection submodule 212. In one embodiment, a first option to reduce the number of pixels rendered in a tile is reconstruction via higher order polynomial interpolation and filtering in a tile to generate missing pixel data for that tile. A second option to reduce the number of pixels rendered in a tile is advection, which includes identifying locations of one or more pixels in a previous frame and reusing pixels from the previous frame for a selected fraction of pixels in the tile.
In one embodiment, pixel data of frame “n” 215 of objects from frame n 220 is saved for possible reuse of pixel data in the next frame “n+1.” Additionally, vertex coordinate data is saved for use in determining a frame-to-frame motion vector of pixels. In one embodiment, the pixel data and vertex coordinates from frame n are stored in a buffer memory for use in the next frame n+1.
Visual artifacts are less perceptible to the human eye in moving objects. Thus, one factor in whether or not a sampling rate may be reduced in a tile is whether the speed is above a threshold speed.
However, certain types of visual artifacts tend to be more noticeable in edges in color. Strictly speaking, detecting color edges in a final image is not possible without rendering the image first. However, it is possible to detect, prior to rendering, a high likelihood of edges in color. Thus, in one embodiment, an edge detection module 305 detects the likelihood of edges in color in local blocks of pixels. That is, regions are detected in which there is a high likelihood of edges in color by assuming that there is a high likelihood of color variation across objects. In one embodiment, Z values from the rasterization of the current frame are analyzed to perform edge detection. A Laplace edge detector may be defined as stencil centered on a current pixel. Any pixel in the tile is marked as having an edge if the Laplacian of the z-buffer at the pixel is greater than a threshold value multiplied by the z-value at the pixel. This defines a one bit value per tile. More generally, any type of edge detection may be used.
In one embodiment, an edge mask is generated for individual tile and an edge state bit may be generated to indicate whether or not a tile includes at least one edge. In one implementation the edge mask is generated for each 4×4 block of pixels although more generally other tile sizes may be used. This information on velocity and the presence of an edge is used by a sample generator 315 to determine a sample pattern for a tile. In one embodiment, a full sampling resolution is utilized if an edge is detected. If no edge is detected and a tile has a speed greater than a first threshold speed a first reduced sampling rate is used. If no edge is detected and a tile has a speed above a second threshold speed a second reduced sampling rate is used. Other additional optional factors could also be considered in making a sampling rate decision. In one embodiment, the sample pattern options include full sample resolution (at least one sample per pixel), one-half resolution (one-half of the pixels sampled in each tile), and one-quarter resolution (one in four pixels sampled in each tile). More generally, a plurality of sampling rates may be provided that are controlled by threshold parameters for each sample rate. Additionally, the sample rates selected may be optimized for the block/tile size selected. Thus, while an illustrative example includes three sample rates of 4, 8, and 16 samples for 4×4 blocks, the approach may be varied based on block size or other considerations to have a set of sampling rates each controlled by threshold parameters for each sample rate. Thus, the number of sampling rates, N, may be more than three, depending on implementation details, such as block/tile size and other factors.
In one embodiment, a dithering module 320 is provided to adjust the sampling pattern from a selection of sampling patterns having the same effective sampling rate. The dithering may be a repetitive sequence (e.g., sample pattern 1, sample pattern 2, sample pattern 3, sample pattern 4) or include aspects of randomization.
The dithering of the sampling pattern by the dithering module 320 reduces the visual perception of sampling artifacts by human users. The human eye and the human brain starts to blend images into a video sequence when the rate is faster than a biological threshold. That is, when images vary at a rate faster than a biological threshold, the human eye blends the images across time and perceives them as a continuously varying sequence, akin to a video. There is some contention about what the exact number of the biological threshold is. At frame rates of about 12 frames per second, the human eye and brain begins to see a sequence of images that is moving instead of individual images. However, somewhat higher frame rates of about 15 frames per second are required to experience the beginnings of relatively fluid (non-jerky) movement. However, the nature of the underlying images is also an additional factor for whether or not a human observer perceives a fluid motion at a given frame rate. Thus, the human eye will tend to average out visual artifacts that are dithered at frames rates at about 12 frames per second and higher. In one embodiment, the dithering is performed such that every pixel is rendered at least fifteen frames per second, which is faster than the human eye can discern individual images. At 60 frames per second, dithering the sample pattern in a tile every four frames corresponds to rendering each pixel at least fifteen frames per second.
In one embodiment, desampling (changing the sample pattern to reduce the sampling rate below one sample per pixel) is permitted if the speed is above a first threshold speed Kfast1. In one embodiment, the sampling rate is permitted to be reduced further if the speed exceeds a second threshold speed Kfast2. A decision whether to perform desampling may also be dependent on other conditions, such as whether or not an edge is detected.
In one embodiment, motion in the camera screen space is obtained by differencing vertex position data from the current frame and the previous frame. A speed regime of a tile is classified on a tile-by-tile basis by calculating a magnitude of a motion vector based on how much a pixel of an object has moved from one frame to another. As previously discussed, in one embodiment, splatting is used in the Z-pass to determine per-pixel motion vectors. In one embodiment, speed thresholds are defined and used as inputs to make decisions on whether adaptive desampling or advection are to be used for a current frame. One speed regime is a quasi-static regime in which an object moves slowly enough that pixels of the object are unlikely to differ significantly from their previous image counterparts. If the speed is within the quasi static speed limit, a decision may be made whether advection may be used to reuse pixels from a previous frame. In one embodiment, an upper bound on a quasi-static speed, Kstat is that a pixel in a given tile (tile m) in frame n remains in the same tile in frame n+1. In one embodiment, if the speed is below Kstat, additional checks are performed to determine if pixels from the previous frame may be used in the current frame. This may include a check that advection produced an acceptable result in a previous frame. Additionally, a check may be performed to check that the pixel values for the tile in the current frame are consistent with a small movement over the previous frame, which may be described as a discrepancy check. An advection discrepancy state bit may be associated with a tile to indicate that it has passed one or more discrepancy checks to confirm that the tile is suitable for advection of at least some of the pixel data.
In the example of
As previously discussed, in one embodiment, a maximum speed is used as condition for whether advection is permitted. In one embodiment, the criteria is that the threshold speed is low enough that local deformation of pixel positions in a small neighborhood can be classified as a rigid transformation in which the change in the positions of the pixels can be represented using one of translation and one rotation for an entire set of pixels to within a desired accuracy. For example, the maximum speed for advection can be that the magnitude of the pixel motion is below a threshold of k pixels. While it is possible that rigid transforms can occur at any speed, the likelihood decreases with increasing speed such that a speed threshold may be used as a criterion for when advection is likely to be beneficial. A discrepancy check may be performed for individual tiles to determine whether or not advection produces acceptable results. This discrepancy check can be performed in a current frame and written out as 1-bit value for each tile to let the next frame make a determination whether or not to disable advection in the neighborhood of a tile failing the discrepancy check if the check indicates that the advected results were inaccurate. That is, in this implementation advection is performed for a tile in frame n and the discrepancy check is performed in frame n and consumed by frame n+1. Frame n+1 then uses the discrepancy check (computed in frame n) to determine whether or not to perform advection in a neighborhood about a tile in frame n+1. If the discrepancy check in frame n indicates the advection result was acceptable, advection is allowed in frame n+1. If not, the advection is turned off for a selected number of frames. The discrepancy check is a check based on whether or not there is a significant change in pixel values of a tile inconsistent with the underlying assumptions of valid advection. If the pixels of an object are slowly moving then a tile is not expected to change significantly between two frames. Thus, if the state of the tile changes significantly the discrepancy check fails. A tile state discrepancy bit (e.g., 0 or 1) may be used to indicate whether the discrepancy check passes. The degree to which changes in tile state are permitted may be determined empirically or heuristically based, for example, in the tradeoffs between the computational benefits of advection and minimizing the appearance of visual artifacts.
Other ways of performing discrepancy checks may be utilized. There are computational benefits to performing advection in a tile of current frame n, performing a discrepancy check, and then utilizing the discrepancy check to determine whether or not to perform advection in frame n+1. However, it will be understood that an alternate implementation of the discrepancy check may be utilized in which the discrepancy check is performed in frame n and used to determine whether or not to utilize advection in frame n to reuse pixels from the previous frame.
The accuracy can be improved, if desired, using various enhancements. In one embodiment, back and forth error correction and compensation (BFECC) is used. BFECC utilized the position determined from Semi-Lagrangian advection and adds the velocity at that coordinate to obtain a new position in the current frame. If there is no error, this coordinate should be the same as the original position (x, y). Otherwise by subtracting half of this error from (x−vx, y−vy) the second-order accurate estimate of the position is obtained, which is accurate to half a pixel, assuming that the velocity is pixel accurate.
Traditionally, a GPU utilizes bilinear interpolation. However, bilinear interpolation has various drawbacks. In one embodiment higher order polynomials having an order of at least three, such as piece-wise cubic polynomials (also known as cubic splines), are used for efficient reconstruction of sparse samples.
Higher order polynomials, such as cubic splines, are capable of mapping a larger frequency spectrum than bilinear interpolation and provide a greater fidelity of data reconstructed from sub-sampled blocks. Additionally, when using bilinear interpolation, samples are preferred on both sides of a pixel, as one-sided linear extrapolation is likely to be inaccurate, as well as exceed the color spectrum range. In contrast, higher-order polynomials using a wider support (>1 pixel away) are more likely to accurately approximate the functional form of rendered image data. While a variety of higher order polynomials could be used, a cubic spline has continuity characteristics that are superior to those of quadratic polynomials. Owing to the edge-detection step performed before desampling, a tile undergoing reconstruction is unlikely to have sharp discontinuities, where higher-order polynomial reconstruction may perform poorly.
One aspect of performing sub-sampling is that there is a sparse sample data at runtime. In an individual block region, such as a k×k pixel region, desampling may result in a subset of pixels being rendered, such as four or eight pixels from a 4×4 block of pixels. The missing pixel data needs to be reconstructed. A consequence of having predetermined sample patterns is that there is a finite set of possible sample locations. This permits a fixed set of local stencils to be created prior to run time, stored, and used to reconstruct pixel data using cubic splines or other higher order polynomials. Conventional approaches to evaluating higher order polynomials in hardware are computationally expensive. In contrast, in embodiments of the present invention the use of a fixed set of pre-computed stencils eliminates the computational overhead, during runtime, of performing a conventional higher order polynomial evaluation. The use of a static set of samples allows the determination of possible configurations of pixels that may need to be reconstructed, and thus the stencils that are required may be pre-computed.
In one embodiment, the higher order polynomial interpolation is implemented as static stencil operations using pre-computed weights. In one embodiment, a table of stencils is stored and made available for spatial reconstruction to a reconstruction submodule 211 of the reconstruction and advection stage 210. The table of stencils provides weights based on known sample positions. In one embodiment, the table of stencils has all of the stencil weights pre-computed for each pixel location within a defined sample pattern. The pre-computed weights permit the higher order polynomial reconstruction to be performed using static stencil operations.
In one embodiment, a set of 5×5 stencils is determined for all possible pixel locations in a tile (e.g., a 4×4 tile) that may need to be interpolated during runtime. Each 5×5 stencil is computed for each pixel location and neighbor configuration. Each stencil provides a list of weight values and corresponding locations of sample points. The stencils are stored in a constant memory table available for reconstruction purposes to the reconstruction submodule 211 of reconstruction and advection stage 210. In one embodiment, at run-time, for each pixel that must be interpolated, an index is computed into this table using the pixel coordinates and sampling mask. In one implementation each stencil is addressed using: (a) the pixel's location within a tile, and (b) the sampling mask used for rendering. Thus, if dithering is employed, the stencil that is selected will depend on which sample pattern is selected for a given degree of subsampling.
In one embodiment, the higher order polynomial interpolation is performed using a multiplier/adder to accumulate the products of weights and sample color values. The accumulated value is then normalized by division, which in many cases can be performed by a bit shift for integer formats, or by subtraction for floating point formats. Thus, the use of stencils with pre-computed weights permits the higher order polynomial interpolation to be calculated at run time with comparatively little computational effort.
An example of a cubic spline functions used for reconstruction and to compute pixel color values as a weighted sum of known pixel color values is as follows:
In one embodiment, a formula expressing a weighted sum to determine a pixel color value is based weights w( ) as follows:
Where c (i, j) is the color value at pixel location (i, j), w( ) is the two dimensional spline function and “Filled” is the set of pixels that were rendered. The two-dimensional spline function is the product of two one-dimensional spline functions or w(i, j)=k(i)k(j), where the one dimensional spline function k( ) is based on the cubic filter formula described in the paper by Don P. Mitchell and Arun N. Netravali, “Reconstruction Filters in Computer Graphics,” Computer Graphics, Volume 22, Number 4, August 1988, pp. 221-228:
The distances in the Mitchell and Netravali paper are defined in a scaled pixel space:
B=C=⅓
By restricting the relative positioning of sample points, the weights and denominators can be pre-computed into stencils. Because the spline function is defined in a bounded fashion, scaling of the magnitude of x can be used to extend the functions to a desired support radius, such as a 2 pixel support radius.
For a tile of size n×n, it is possible to arrange a k×k square in (n/k)*(n/k) possible configurations. A sampling rate of 4*s requires s squares, leading to (n*n)/(k*k*s) sampling patterns.
Assuming a 5×5 stencil, there are 24 values in all to be pre-computed (the center is always 0 since the pixel itself has no color value). Of those, at most half can be rendered if 8 samples are used per 4×4 tile, leaving 12 values. In one embodiment, each stencil is implemented to include a 4-bit count of the number of non-zero weights, followed by 8-bit weights stored in one chunk, followed by 2 chunks of 3-bit coordinate offsets for the x and y coordinates relative to the center.
In one embodiment, the stencils are stored in order of the sampling patterns. In one embodiment the different sampling patterns for the same sampling rate are rotations of each other, so there are two sets of patterns. These can be stored in row major order within a 4×4, with an index list to point to the data for pixel (i, j). For rotations of the sampling mask, the coordinates are transformed appropriately.
Referring to
In one embodiment, each sampling pattern is defined as combination of sparse square patterns (e.g., four samples to be rendered in a square pattern). Selecting square patterns is useful in applications in which groups of four pixels (quads) are a default unit of processing. However, more generally other arrangements of sampling positions could be used in the sampling patterns. In one embodiment, the sample patterns are squares of size 3×3 within a 4×4 tiles. Thus, adjacent vertices are 2 pixels apart along each axis.
In one embodiment, the same sampling pattern is used for all regions of an individual frame that are sub-sampled at a given sampling rate. In this embodiment, the same sampling pattern is used in all of the tiles sub-sampled at a given sample rate because this results in consistent spacing of sample positions within every frame, simplifying reconstruction routines.
In one embodiment, the sampling patterns are based on quads to exploit Single instruction multiple data (SIMD) processing units. Consistent spacing of samples provides robust interpolation and aids in achieving full pixel resolution in a final image.
Embodiments of the present invention may be used to generate a single (non-stereoscopic) display. However, it may also be applied to stereoscopic rendering for virtual reality applications. Referring to
While the invention has been described in conjunction with specific embodiments, it will be understood that it is not intended to limit the invention to the described embodiments. On the contrary, it is intended to cover alternatives, modifications, and equivalents as may be included within the spirit and scope of the invention as defined by the appended claims. The present invention may be practiced without some or all of these specific details. In addition, well known features may not have been described in detail to avoid unnecessarily obscuring the invention. In accordance with the present invention, the components, process steps, and/or data structures may be implemented using various types of operating systems, programming languages, computing platforms, computer programs, and/or computing machines. In addition, those of ordinary skill in the art will recognize that devices, such as hardwired devices, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), or the like, may also be used without departing from the scope and spirit of the inventive concepts disclosed herein. The present invention may also be tangibly embodied as a set of computer instructions stored on a computer readable medium, such as a memory device.
The present application claims the benefit of U.S. Provisional Patent Application No. 62/018,221, filed Jun. 27, 2014; U.S. Provisional Patent Application No. 62/018,228, filed Jun. 27, 2014; U.S. Provisional Patent Application No. 62/018,254 filed Jun. 27, 2014; and U.S. Provisional Patent Application No. 62/018,274, filed Jun. 27, 2014, the contents of each of which are hereby incorporated by reference.
Number | Date | Country | |
---|---|---|---|
62018221 | Jun 2014 | US | |
62018228 | Jun 2014 | US | |
62018254 | Jun 2014 | US | |
62018274 | Jun 2014 | US |