# Patent

## Title
Resampling system and method for graphics data including sine-wave components

## Abstract

A method and system for calculating resample output values from input samples and their associated sample values. A resampling circuit calculates a frequency value for a sine-wave model from a sample set of the input samples and determines whether the frequency value is in a frequency range. In the case where the frequency value is in the frequency range, a sinusoidal transition model is determined based on the sample set. However, if the frequency value is outside of the frequency range, a non-sinusoidal model is determined based on the sample set. The resampling circuit then calculates resample output values from the resulting sinusoidal or non-sinusoidal model.

## Background
CROSS-REFERENCE TO RELATED APPLICATION
This application is a continuation of U.S. patent application Ser. No. 09/779,010, filed Feb. 7, 2001, now U.S. Pat. No. 7,053,896.


TECHNICAL FIELD
The present invention is related generally to the field of computer graphics, and more particularly, a system and method for resampling graphics data of a source image to produce a destination image.
BACKGROUND OF THE INVENTION
As display devices of various sizes and increased resolution have been developed and the demand for them have increased, the ability for a graphics processing system to resize and resample source images and create destination images to take advantage of the various sized and higher resolution displays is a desirable operation. In an electronic display system, color at each pixel is represented by a set of color components, and each color component is represented by a sample value. Color components such as red, green, blue (RGB) or other representations such as YCbCr are well known in the art. Whichever representation is chosen, each color component can be interpreted as a two dimensional array of samples, so three such arrays can represent images on display systems. Conceptually, resampling can be viewed as a spatial process, working on discrete input samples, represented by pixels of the source image arranged in a two-dimensional bitmap. The output samples of the destination image are spatially located at fractional sample positions within the input sample grid. Various interpolation and modeling methods are used to construct transition models between samples of the source image from which additional graphics data is produced during the resampling operation.
The additional graphics data is then used to produce larger or higher resolution destination graphics images. However, the resulting destination image must retain an acceptable image quality with respect to the source image. That is, the destination image should appear to retain at least a similar visual qualities of the source image, such as having nearly the same color balance, contrast, and brightness as the original source image. Otherwise, rather than accurately reproducing a larger or higher resolution graphics image of the source image, the resampling operation will compromise image quality by introducing image distortion. To this end, various resampling algorithms have been developed in order to create high quality destination graphics images.
With many conventional resampling algorithms, a transition model between input samples along each axis is constructed to provide output sample values. Generally good results can be obtained with separable processing along each axis for graphics images because image feature cross-sections have the same characteristics when viewed at any angle within the image plane, only at different effective sample rates. The transition models between the input samples are constructed such that the output samples interpolated from the transition model create a destination image that closely resembles the original or source image. The transition models are typically continuous so that an output sample can be generated at any position between the input samples.
Although an axis separable cubic model between two input samples can provide a model with very desirable reconstruction characteristics, algorithms for resampling and sharpening graphics data representing video often are not suitable for resizing and resampling graphics data representing test patterns containing sine-wave components. Such test patterns are called zone plates, and are characterized by a frequency component along each axis, each of which is a function of position within the pattern. The position and frequency functions are designed to change frequencies smoothly and continuously with position.
Zone plates may be embedded within patterns testing various other attributes of a video camera, storage, transmissions or display system. They are effective in testing systems with analog components (e.g., analog modulated terrestrial broadcasting), and may provide some useful tests for spectrally based compression systems (such as DCTs used in MPEG). However, these tests generally do not correspond to any attributes of the human visual system. Nevertheless, the human eye is very adept at observing large areas of inconsistency in the presentation of these patterns. Thus, to avoid viewer complaints or feelings of disappointment (whether or not they are justified), a graphics processing system having resampling and resizing capabilities should be able to accommodate these test patterns.
Therefore, there is a need for a method and system for resampling graphics data of images having sine-wave components.
SUMMARY OF THE INVENTION
The present invention relates to a method and system for calculating resample output values from input samples and their associated sample values. A resampling circuit calculates a frequency value for a sine-wave model from a sample set of the input samples and determines whether the frequency value is in a frequency range. In the case where the frequency value is in the frequency range, a sinusoidal transition model is determined based on the sample set. However, if the frequency value is outside of the frequency range, a non-sinusoidal model is determined based on the sample set. The resampling circuit then calculates resample output values from the resulting sinusoidal or non-sinusoidal model.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a block diagram of a computer system in which embodiments of the present invention are implemented.
FIG. 2 is a block diagram of a graphics processing system in the computer system of FIG. 1.
FIG. 3 is a block diagram of a resampling circuit in the graphics processing system of FIG. 2 according to an embodiment of the present invention.
FIG. 4 is a diagram representing a sample of graphics data and corresponding sample values.



DETAILED DESCRIPTION OF THE INVENTION
Embodiments of the present invention provide a method and system for calculating resampled values from a source graphics image having graphics data including sine-wave components. Certain details are set forth below to provide a sufficient understanding of the invention. However, it will be clear to one skilled in the art that the invention may be practiced without these particular details. In other instances, well-known circuits, control signals, timing protocols, and software operations have not been shown in detail in order to avoid unnecessarily obscuring the invention.
FIG. 1 illustrates a computer system 100 in which embodiments of the present invention are implemented. The computer system 100 includes a processor 104 coupled to a host memory 108 through a memory/bus interface 112. The memory/bus interface 112 is coupled to an expansion bus 116, such as an industry standard architecture (ISA) bus or a peripheral component interconnect (PCI) bus. The computer system 100 also includes one or more input devices 120, such as a keypad or a mouse, coupled to the processor 104 through the expansion bus 116 and the memory/bus interface 112. The input devices 120 allow an operator or an electronic device to input data to the computer system 100. One or more output devices 120 are coupled to the processor 104 to provide output data generated by the processor 104. The output devices 124 are coupled to the processor 104 through the expansion bus 116 and memory/bus interface 112. Examples of output devices 124 include printers and a sound card driving audio speakers. One or more data storage devices 128 are coupled to the processor 104 through the memory/bus interface 112 and the expansion bus 116 to store data in, or retrieve data from, storage media (not shown). Examples of storage devices 128 and storage media include fixed disk drives, floppy disk drives, tape cassettes and compact-disc read-only memory drives.
The computer system 100 further includes a graphics processing system 132 coupled to the processor 104 through the expansion bus 116 and memory/bus interface 112. Optionally, the graphics processing system 132 may be coupled to the processor 104 and the host memory 108 through other types of architectures. For example, the graphics processing system 132 may be coupled through the memory/bus interface 112 and a high speed bus 136, such as an accelerated graphics port (AGP), to provide the graphics processing system 132 with direct memory access (DMA) to the host memory 108. That is, the high speed bus 136 and memory bus interface 112 allow the graphics processing system 132 to read and write host memory 108 without the intervention of the processor 104. Thus, data may be transferred to, and from, the host memory 108 at transfer rates much greater than over the expansion bus 116. A display 140 is coupled to the graphics processing system 132 to display graphics images. The display 140 may be any type of display, such as a cathode ray tube (CRT), a field emission display (FED), a liquid crystal display (LCD), or the like, which are commonly used for desktop computers, portable computers, and workstation or server applications.
FIG. 2 illustrates circuitry included within the graphics processing system 132 for performing various three-dimensional (3D) graphics functions. As shown in FIG. 2, a bus interface 200 couples the graphics processing system 132 to the expansion bus 116. In the case where the graphics processing system 132 is coupled to the processor 104 and the host memory 108 through the high speed data bus 136 and the memory/bus interface 112, the bus interface 200 will include a DMA controller (not shown) to coordinate transfer of data to and from the host memory 108 and the processor 104. A graphics processor 204 is coupled to the bus interface 200 and is designed to perform various graphics and video processing functions, such as, but not limited to, generating vertex data and performing vertex transformations for polygon graphics primitives that are used to model 3D objects. The graphics processor 204 is coupled to a triangle engine 208 that includes circuitry for performing various graphics functions, such as clipping, attribute transformations, rendering of graphics primitives, and generating texture coordinates for a texture map.
A pixel engine 212 is coupled to receive the graphics data generated by the triangle engine 208. The pixel engine 212 contains circuitry for performing various graphics functions, such as, but not limited to, texture application or mapping, bilinear filtering, fog, blending, and color space conversion. A memory controller 216 coupled to the pixel engine 212 and the graphics processor 204 handles memory requests to and from an local memory 220. The local memory 220 stores graphics data, such as source pixel color values and destination pixel color values. A display controller 224 is coupled to the memory controller 216 to receive processed destination color values for pixels that are to be rendered. Coupled to the display controller 224 is a resampling circuit 228 that facilitates resizing or resampling graphics images. As will be explained below, embodiments of the resampling circuit 228 perform approximations that simplify the calculation of a model between two sample points for use during resampling. The output color values from the resampling circuit 228 are subsequently provided to a display driver 232 that includes circuitry to provide digital color signals, or convert digital color signals to red, green, and blue analog color signals, to drive the display 140 (FIG. 1).
Although the resampling circuit 228 is illustrated as being a separate circuit, it will be appreciated that the resampling circuit 228 may also be included in one of the aforementioned circuit blocks of the graphics processing system 132. For example, the resampling circuit 228 may be included in the graphics processor 204 or the display controller 224. In other embodiments, the resampling circuit 228 may be included in the display 140 (FIG. 1). Therefore, the particular location of the resampling circuit 228 is a detail that may be modified without deviating from the subject matter of the invention, and should not be used in limiting the scope of the present invention.
FIG. 3 illustrates a resampling circuit 300 that may be substituted for the resampling circuit 228 shown in FIG. 2. The resampling circuit 300 includes a sine-model resampling circuit 312 for determining if a sample of graphics data provided by the display driver 224 (FIG. 2) is likely to include sine-wave components. As mentioned previously, the resampling operations used for resampling graphics data including sine-wave components is often different than that used for resampling other graphics data representing other types of graphics images, such as video. Although one resampling algorithm may be used, the image quality of one or the other types of graphics data will be compromised. Thus, two different resampling operations are used for the different types of graphics data, one for graphics data including sine-wave components performed by the sine-model resampling circuit 312 and one for non-sine-wave graphics data performed by a non-sine-wave resampling circuit 308.
It will be appreciated that the sample values for the samples may consist of several different components. For example, the sample value may represent pixel colors which are the combination of red, green, and blue color components. Another example includes sample values representing pixel colors which are the combination of luma and chroma components. Consequently, because it is well understood in the art, although circuitry to perform graphics operation for each of the components is not expressly shown or described herein, embodiments of the present invention include circuitry, control signals, and the like necessary to perform resampling operations on each component for multi-component sample values. Moreover, it will be appreciated that embodiments of the present invention further include the circuitry, control signals, and the like necessary to perform axis separable resampling operations for graphics data represented in multiple axes. Implementation of axis separable resampling is well understood in the art, and a more detailed description of such has been omitted from herein to avoid unnecessarily obscuring the present invention.
The non-sine-wave resampling circuit 308 can perform conventional resampling operations that are well known to those of ordinary skill in the art. Alternatively, a resampling operation such as that described in co-pending application having U.S. Ser. No. 09/760,173, entitled PIXEL RESAMPLIING SYSTEM AND METHOD to Slavin, filed Jan. 12, 2001, which is incorporated herein by reference, can also be performed by the non-sine-wave resampling circuit 308. In summary, the subject matter of the aforementioned patent application includes generating a cubic model for transitions between adjacent samples from the sample values and the gradient values cosited with the two samples. The cosited gradients are approximated to facilitate generation of the transition model. The coefficients for the cubic model are determined from the known values and used by a cubic model evaluation circuit to calculate resampled values between the adjacent samples. As will be explained in more detail below, the cubic model evaluation circuit described in the aforementioned patent application may be used with the present invention to determine resampled values for graphics data including sine-wave components.
In operation, when a resampling operation is to be performed, the resampling circuit 228 (FIG. 2), becomes active and sine-model resampling circuit 312 receives sample values for the samples of graphics data of the source image to be resampled. As will be explained in more detail below, based on a sampling of the graphics data received by the resampling circuit 300, the sine-model resampling circuit 312 determines whether the graphics data includes sine-model components. If so, the sine-model resampling circuit 312 performs the resampling operation on the graphics data. All other graphics data is provided to the non-sine-model resampling circuit 308 for the resampling operation. The sine-model resampling circuit 308 performs operations of fitting a sine-model to the sample of graphics data it receives. Each of the respective resampling circuits perform various operations to resample the graphics data received from display controller 224 (FIG. 2) to produce a resampled graphics image. The resampled data generated by the resampling circuits are subsequently used in the process of rendering enlarged or higher resolution graphics images.
Although graphics data including sine-wave components may change frequency with position, such as in a zone plate test pattern, the sine-model resampling circuit 312 performs the operation with localized processing. Thus, the zone plate can be regarded as having a fixed frequency in each axis over a small region. For the small region, algorithms can be used to find the parameters for the equation:

Vp=A sin(ωp+φ)+B

where p is a local input sample position value along each axis, and Vp is an input sample value at position p. Although the previous equation has four unknowns, and consequently requires only four adjacent sample values, for reasons that will be explained later, we use five samples along each axis with a position index p of zero as the center of the samples. Initially, a set of four samples S0 . . . 3=V−2 . . . 1 is selected. The values of the selected sample set are used to solve the following equations to obtain angular frequency ω:




  
    
      d
      1
    
    =
    
      
        
          S
          0
        
        -
        
          
            S
            3
          
          .
          
            

          
          ⁢
          
            d
            2
          
        
      
      =
      
        
          
            S
            1
          
          -
          
            
              
                S
                2
              
              .
              
                

              
              ⁢
              cos
            
            ⁢
            
              
            
            ⁢
            w
          
        
        =
        
          
            (
            
              
                
                  d
                  1
                
                
                  d
                  2
                
              
              -
              1
            
            )
          
          2
        
      
    
  




  
    
      
        
          if
          ⁢
          
            
          
          ⁢
          
            (
            
              cos
              ⁢
              
                
              
              ⁢
              w
            
            )
          
        
        <
        
          
            -
            0.95
          
          ⁢
          
            
          
          ⁢
          then
          ⁢
          
            
          
          ⁢
          ω
        
      
      =
      
        arccos
        ⁢
        
          
        
        ⁢
        
          (
          
            -
            0.95
          
          )
        
      
    
    ,
    
      

    
    ⁢
    
      
        
          else
          ⁢
          
            
          
          ⁢
          if
          ⁢
          
            
          
          ⁢
          
            (
            
              cos
              ⁢
              
                
              
              ⁢
              w
            
            )
          
        
        <
        
          0.9
          ⁢
          
            
          
          ⁢
          then
          ⁢
          
            
          
          ⁢
          ω
        
      
      =
      
        arccos
        ⁢
        
          
        
        ⁢
        
          (
          
            cos
            ⁢
            
              
            
            ⁢
            w
          
          )
        
      
    
  




  
    else
    ⁢
    
      
    
    ⁢
    NOT
    ⁢
    
      -
    
    ⁢
    A
    ⁢
    
      -
    
    ⁢
    
      SINE
      .
    
  



The value of the angular frequency ω is limited to ω≧acos(−0.95) to prevent the value from going too near π=a cos(−1), the maximum angular frequency which causes ill-conditioned behavior at later stages of processing. Although the frequency limit ω≧a cos(−0.95) may introduce minor errors during the following sine-model fit operation, which will be described below, the frequency limit creates the appearance of a gradual and benign “fade-out” on zone-plate patterns near π. It will be appreciated, however, that limit values nearer to −1. are possible with low-noise, higher accuracy data.

In the case where d2 is zero, the samples are positioned symmetrically around a peak midway between samples S1 and S2. Such a situation presents an infinity of sine-wave solutions, and consequently, poorly conditioned equations. However, as shown in FIG. 4, this issue may be addressed by providing the sine-model resampling circuit 312 (FIG. 3) using five samples instead of four. That is, d2 is evaluated by the sine-model resampling circuit 312 from the middle two samples for each of the candidate sets of four samples:

S0 . . . 3=V−2 . . . 1

S0 . . . 3=V−1 . . . 2

The set {V} with the largest |d2| is selected and used to obtain a reliable estimate of the angular frequency ω.

If the maximum of d2 from both sets of four samples still results in an angular frequency ω that is near 0 (i.e., cos(ω)>0.9), then the samples are ill-conditioned, most likely the result from sine-waves components of very low amplitude or frequency. As a result, a NOT-A-SINE error is returned for the five samples. It will be appreciated that the limit of cos(ω)>0.9 may be modified for different noise and accuracy conditions. However, using the present limit will typically result in one of the sets of sample values {S} yielding a useful d2 value, and thus, provide good measurement results. Where a NOT-A-SINE error is produced, the graphics data is provided to the non-sine-model resampling circuit 308 where an alternative interpolation algorithm is performed instead. As mentioned previously, various suitable interpolation algorithms may be performed there.
Once a set of {S} values has been selected by the sine-model resampling circuit 312 and the angular frequency ω obtained, a sine fit can be obtained by finding {A, φ, B} from the sine-model equation:

Vp=A sin(ωp+φ)+B.

The values for amplitude A, phase φ, and offset B can be solved by the sine-model resampling circuit 312 using values that are already known, namely, the angular frequency ω, and the sample values of the middle three samples {V−1, V0, V1} of the five samples previously mentioned. While it would be possible to perform a least-squares fit to more than three samples, using a three-sample fit provides the benefit of simplicity, and additionally, ensures that the resulting model will go through the original three sample points. Moreover, as will be explained in further detail below, additional tests can be performed by the sine-model resampling circuit 312 on the resulting three sample fit model to confirm that it is not fitting sine-models to transitions between samples of graphics data not including sine-wave components. Solving the three-sample point equations results in:



  
    
      Acos
      ⁢
      
        
      
      ⁢
      
        (
        φ
        )
      
    
    =
    
      
        A
        ⁢
        
          
        
        ⁢
        COS
      
      =
      
        
          
            V
            1
          
          -
          
            V
            
              -
              1
            
          
        
        
          2
          ⁢
          sin
          ⁢
          
            
          
          ⁢
          
            (
            ω
            )
          
        
      
    
  




  
    
      Asin
      ⁢
      
        
      
      ⁢
      
        (
        φ
        )
      
    
    =
    
      ASIN
      =
      
        
          
            V
            1
          
          +
          
            V
            
              -
              1
            
          
          -
          
            2
            ⁢
            
              V
              0
            
          
        
        
          2
          ⁢
          
            
          
          ⁢
          
            (
            
              
                cos
                ⁢
                
                  
                
                ⁢
                
                  (
                  ω
                  )
                
              
              -
              1
            
            )
          
        
      
    
  




  
    B
    =
    
      
        V
        0
      
      -
      ASIN
    
  



which provides the offset B directly. The phase and amplitude can then be obtained directly through a rectangular to polar coordinate conversion:

φ=arctan 2(A SIN,A COS),

A=√{square root over ((A SIN)2+(A COS)2)}{square root over ((A SIN)2+(A COS)2)}

After the sine-model resampling circuit 312 resolves the {A, φ, B} values from the previous equations, the resulting sine-model can be evaluated to directly obtain resampled values from the source image. Note that the phase φ is coincident with the middle of the three samples values V0, and that the four-quadrant a tan 2(y,x) function is used. Further note that in the case where ω=0 or ω=π, a division by zero occurs. However, these values should have been excluded previously.

An alternative approach to determining resampled values according to a sine-model results from applying the A SIN and A COS values used in resolving the offset value B. Expanding the sine-model equation discussed earlier results in:

Rp=A sin(φ)cos(ωp)+A cos(φ)sin(ωp)+B

where Rp=Vp for p={−1,0,1}. As discussed previously, the values for A sin(φ) and A cos(φ), and the angular frequency ω were determined to calculate the offset value B. Thus, Rp can be evaluated at any fractional position p=Δp by substituting these values into the expanded sine-model equation to obtain a resampled result in each axis between the samples V−1 and V0.

As a means of verifying the accuracy of the sine-model generated through the three samples {V−1, V0, V1}, the model is evaluated at the positions of the first and last of the five samples (i.e., at positions V−2 and V2) using the following equation:



  
    
      diff
      A
    
    =
    
      
      
        
          R
          
            -
            2
          
        
        -
        
          V
          
            -
            2
          
        
      
      
    
  




  
    
      diff
      B
    
    =
    
      
      
        
          R
          2
        
        -
        
          V
          2
        
      
      
    
  




  
    threshold
    =
    
      A
      4
    
  




  
    if
    ⁢
    
      
    
    ⁢
    
      (
      
        
          (
          
            
              diff
              A
            
            >
            threshold
          
          )
        
        ⁢
        
          
        
        ⁢
        or
        ⁢
        
          
        
        ⁢
        
          (
          
            
              diff
              B
            
            >
            threshold
          
          )
        
      
      )
    
    ⁢
    
      
    
    ⁢
    then
    ⁢
    
      
    
    ⁢
    NOT
    ⁢
    
      -
    
    ⁢
    A
    ⁢
    
      -
    
    ⁢
    
      SINE
      .
    
  



The threshold value is set to a fraction of the amplitude of the fitted sine wave, which allows for some noise and distortions due to assumptions that the angular frequency ω is constant, or that ω may have been limited near π as previously discussed. A scaling value of 1⁄4 works quite well, and is easy to implement, but it will be appreciated that other values are possible depending upon noise levels. This test rejects fits on edges because the outlying samples will fit badly to a sine-model which was fitted to the central three samples {V−1, V0, V1}.
Note that cos(−x)=cos(x) and sin(−x)=−sin(x). Consequently, cos(−2ω) and sin(−2ω) can be calculated by sharing look-up tables when obtaining R−2 and R2. Moreover, diffA can be determined using two ROM tables to obtain the sin(2ω) and cos(2ω) values, along with two multiplies and two adders. As just discussed, only two more multipliers and adders are needed to obtain diffB.
As an alternative, rather than calculating the amplitude A precisely using the equation:

A=√{square root over ((A SIN)2+(A COS)2)}{square root over ((A SIN)2+(A COS)2)}

which involves division operations and makes the calculations more difficult and complex to solve, a usable amplitude A can be approximated for the verification operation because the value is used only as a threshold for determining the accuracy of the resulting sine-model. An economical approximation of the amplitude A to better than 5% accuracy can be obtained using:




  
    s
    =
    
      
      ASIN
      
    
  




  
    c
    =
    
      
      ACOS
      
    
  




  
    
      if
      ⁢
      
        
      
      ⁢
      
        (
        
          s
          >
          c
        
        )
      
      ⁢
      
        
      
      ⁢
      A
    
    ≈
    
      s
      +
      
        
          c
          2
        
        ⁢
        
          
        
        ⁢
        else
        ⁢
        
          
        
        ⁢
        A
      
    
    ≈
    
      c
      +
      
        
          s
          2
        
        .
      
    
  



An incrementer (for 2's complement negation), and a multiplexer can be used to obtain the absolute value of s and c. A compare, a multiplexer, and an adder are used for the remaining operations.

As mentioned previously, resampled values for a sine-model may be directly determined from the sine-model equation:

Rp=A sin(φ)cos(ωp)+A cos(φ)sin(ωp)+B.

However, the arithmetic for directly obtaining the resampled value is relatively complex, so the resulting system is expensive in hardware. As an alternative to solving the sine-model directly, a cubic model system may be used to determine resampled values. This method of determining the resampled values may be desirable where a resampling circuit is equipped with an cubic model evaluation block. The resampling operation employs a conventional cubic evaluation circuit, which is well known in the art. Although not described in greater detail herein, implementation of a cubic model evaluation block is well understood by those of ordinary skill in the art, and the description provided herein is sufficient to allow one to practice the invention without undue experimentation. Additionally, as mentioned previously, a cubic evaluation circuit suitable for implementing embodiments of the present invention is included in the system described in the aforementioned co-pending patent application, PIXEL RESAMPLING SYSTEM AND METHOD.
A cubic model may be used between two input samples p and p+1 to provide a continuous model having desirable reconstruction characteristics for graphics images. A piece-wise cubic polynomial model along an axis will be valid over a fractional input sample position Δp from 0 to 1. Consequently, the model is valid from integer sample position p to p+1:




  
    f
    ⁢
    
      
    
    ⁢
    
      (
      
        P
        +
        
          Δ
          ⁢
          
            
          
          ⁢
          p
        
      
      )
    
  
  =
  
    
      ∑
      
        i
        =
        0
      
      3
    
    ⁢
    
      
        C
        ⁡
        
          [
          
            P
            ,
            i
          
          ]
        
      
      ⁢
      
        
          
            (
            
              Δ
              ⁢
              
                
              
              ⁢
              p
            
            )
          
          i
        
        .
      
    
  




The resulting cubic model will go through the two input samples p and p+1.

As is well known, a cubic model can be solved with four constraints. Two of these constraints may be provided by the sample values fp and fp+1 at the two input samples p and p+1. These sample values are known. Two additional constraints may be provided by the gradients grp and grp+1 at, or co-sited with, the two input samples p and p+1. To solve the cosited gradients, the equation for the cubic model is differentiated with respect to Δp, resulting in:




  
    gr
    ⁢
    
      
    
    ⁢
    
      (
      
        P
        +
        
          Δ
          ⁢
          
            
          
          ⁢
          p
        
      
      )
    
  
  =
  
    
      ∑
      
        i
        =
        0
      
      3
    
    ⁢
    
      
        iC
        ⁢
        
          
        
        [
        
          P
          ,
          i
        
        ]
      
      ⁢
      
        
          
            (
            
              Δ
              ⁢
              
                
              
              ⁢
              p
            
            )
          
          
            i
            -
            1
          
        
        .
      
    
  




Evaluating the two equations at Δp={0, 1}, and solving for the four coefficients C[P, i] at the relative positions of the contributors to the cubic model are of interest results in coefficients:

k=f1−f0

C3=gr1+gr0−2k

C2=k−C3−gr0

C1=gr0

C0=f0

for the cubic equation:





  
    f
    ⁢
    
      
    
    ⁢
    
      (
      
        Δ
        ⁢
        
          
        
        ⁢
        p
      
      )
    
  
  =
  
    
      ∑
      
        i
        =
        0
      
      3
    
    ⁢
    
      
        
          
            C
            i
          
          ⁡
          
            (
            
              Δ
              ⁢
              
                
              
              ⁢
              p
            
            )
          
        
        i
      
      .
    
  




The resulting cubic equation, along with the gradients gr0 and gr1 and the sample values f0 and f1 for the two input samples p and p+1 provides a piece-wise continuous model for resampling.

Differentiating the sine-model equation with respect to the angular frequency ω to find the gradients grp results in:

grp=−A sin(φ)×ω sin(φp)+A cos(φ)×ω cos(ωp).

This model can obtain valid gradients at position p={−1,0,1}, cosited with the original fitted samples. The gradients are then passed to the cubic evaluation block to generate a resampled output point. This approach is less accurate than calculating resampled values directly through a sine-model fit because the cubic interpolation system cannot approximate the significant higher order polynomial terms in Δp that are present in sine waves at higher frequencies. This distortion along the x-axis further compounds errors along the y-axis. However, good results can be obtained up to near 0.9 of the Nyquist sampling limit. Moreover, although two output values (gradients) are evaluated instead of one for the sine model case, the values are cosited with the input samples at discrete sample times, so as p is an integer, the hardware to evaluate grp is much simpler. Note that the cubic evaluation circuit which follows should be there in any case for non-sinusoids.

As mentioned previously, embodiments of the invention have been described herein with sufficient detail to allow a person of ordinary skill in the art to practice the invention. Implementation of many of the algorithms previously described may be implemented by conventional circuitry. For example, determining the angular frequency ω can be implemented using logarithm ROMs, and the corresponding anti-logarithm and limit detection can be built into another ROM. Thus, only three ROMs and three address to obtain ω once a data {S} set has been selected. Another example is using a comparator to determine the largest |d2| calculated for the two sets of samples and a multiplexer to select the final data set of {S} to estimate the angular frequency ω. Thus, in order to prevent unnecessarily obscuring the invention, a more detailed description of the implementation of various aspects of the invention have been omitted from herein.
From the foregoing it will be appreciated that, although specific embodiments of the invention have been described herein for purposes of illustration, various modifications may be made without deviating from the spirit and scope of the invention. Accordingly, the invention is not limited except as by the appended claims.

## Claims

1. A graphics processing system, comprising:
a graphics processor operable to generate data representing graphics primitives;
a triangle engine coupled to the graphics processor and operable to render the graphics primitives;
a pixel engine coupled to the triangle engine and operable to generate data representing pixels of an image; and
a resampling circuit coupled to the pixel engine to provide resample output values, the resampling circuit operable to calculate from a sample set of the sample values an angular frequency value w for a sine-wave model and determine whether the frequency value w is in a frequency range, where the frequency value w is in the frequency range, the resampling circuit is further operable to determine from the sample set a sinusoidal model from which the resample output values are calculated, where the angular frequency value w is not in the frequency range, the resampling circuit is operable to determine from the sample set a non-sinusoidal model from which the resample output sample values are calculated and calculate resample output sample values from the resulting model.

2. The graphics processing system of claim 1 wherein the resampling circuit comprises a resampling circuit operable to determine from the sample set a cubic transition model between two of the input samples from which the resample output sample values are calculated when the angular frequency value w is not in the frequency range.

3. The graphics processing system of claim 1 wherein the resampling circuit comprises a resampling circuit operable to determine whether the frequency value w is in a frequency range between arccos(−0.95) ≦ω<arccos(0.9).

4. The graphics processing system of claim 1 wherein the resampling circuit comprises a resampling circuit operable to calculate the angular frequency value w from a sample set of sample values including first, second, third, fourth, and fifth input sample values and the angular frequency value w is calculated from:




  
    ω
    =
    
      arc
      ⁢
      
        
      
      ⁢
      
        cos
        (
        
          
            (
            
              
                
                  d
                  1
                
                
                  d
                  2
                
              
              -
              1
            
            )
          
          2
        
        )
      
    
  
  ,



where d1=(V1−V2) and d2=(V0−V1) if |V0−V1|>|V1−V0|,
otherwise d1=(V−2−V1) and d2=(V−1−V0),
where V−2, V−1, V0, V1, and V2, are the first, second, third, fourth, and fifth input sample values, respectively.

5. The graphics processing system of claim 4 wherein the resampling circuit comprises a resampling circuit operable to calculate the sine-wave model and the output sample values therefrom from the equation:


Vp=A sin(ωp+φ)+B,

where Vp is the output sample value at position p, ω is an angular frequency calculated from the input sample values,




  
    
      
        
          
            B
            =
            
            ⁢
            
              
                V
                0
              
              -
              
                A
                ⁢
                
                  
                
                ⁢
                SIN
              
            
          
          ,
        
      
    
    
      
        
          
            φ
            =
            
            ⁢
            
              arc
              ⁢
              
                
              
              ⁢
              tan
              ⁢
              
                
              
              ⁢
              2
              ⁢
              
                (
                
                  
                    A
                    ⁢
                    
                      
                    
                    ⁢
                    SIN
                  
                  ,
                  
                    A
                    ⁢
                    
                      
                    
                    ⁢
                    COS
                  
                
                )
              
            
          
          ,
          and
        
      
    
    
      
        
          
            A
            =
            
            ⁢
            
              
                
                  
                    (
                    
                      A
                      ⁢
                      
                        
                      
                      ⁢
                      SIN
                    
                    )
                  
                  2
                
                +
                
                  
                    (
                    
                      A
                      ⁢
                      
                        
                      
                      ⁢
                      COS
                    
                    )
                  
                  2
                
              
            
          
          ,
        
      
    
  




  
    
      where
      ⁢
      
        
      
      ⁢
      A
      ⁢
      
        
      
      ⁢
      COS
    
    =
    
      
        
          
            
              V
              1
            
            -
            
              V
              
                -
                1
              
            
          
          
            2
            ⁢
            
              sin
              ⁡
              
                (
                ω
                )
              
            
          
        
        ⁢
        
          
        
        ⁢
        and
        ⁢
        
          
        
        ⁢
        A
        ⁢
        
          
        
        ⁢
        SIN
      
      =
      
        
          
            
              V
              1
            
            +
            
              V
              
                -
                1
              
            
            -
            
              2
              ⁢
              
                V
                0
              
            
          
          
            2
            ⁢
            
              (
              
                
                  cos
                  ⁡
                  
                    (
                    ω
                    )
                  
                
                -
                1
              
              )
            
          
        
        .

6. The graphics processing system of claim 4 wherein the resampling circuit comprises a resampling circuit operable to calculate the sine-wave model and the output sample values therefrom from the equation:


Rp=A sin(φ)cos(ωp)+A cos(φ)sin(ωp)+B,

where Rp is the output sample value at position p, ω is the angular frequency,


B=V0−A SIN,



φ=arctan2(A SIN,A COS), and



A=√{square root over ((A SIN)2+(A COS)2)}{square root over ((A SIN)2+(A COS)2)},


where





  
    A
    ⁢
    
      
    
    ⁢
    
      cos
      ⁡
      
        (
        φ
        )
      
    
  
  =
  
    
      A
      ⁢
      
        
      
      ⁢
      COS
    
    =
    
      
        
          
            
              V
              1
            
            -
            
              V
              
                -
                1
              
            
          
          
            2
            ⁢
            
              
            
            ⁢
            
              sin
              ⁡
              
                (
                ω
                )
              
            
          
        
        ⁢
        
          
        
        ⁢
        and
        ⁢
        
          
        
        ⁢
        A
        ⁢
        
          
        
        ⁢
        
          sin
          ⁡
          
            (
            φ
            )
          
        
      
      =
      
        
          A
          ⁢
          
            
          
          ⁢
          SIN
        
        =
        
          
            
              
                V
                1
              
              +
              
                V
                
                  -
                  1
                
              
              -
              
                2
                ⁢
                
                  V
                  0
                
              
            
            
              2
              ⁢
              
                (
                
                  
                    cos
                    ⁡
                    
                      (
                      ω
                      )
                    
                  
                  -
                  1
                
                )
              
            
          
          .

7. The graphics processing system of claim 6 wherein the resampling circuit comprises a resampling circuit further operable to verify the accuracy of the sine-wave model by calculating:


diffA=|R−2−V−2| and diffB=|R2−V2|,

the resampling circuit further operable to confirm that diffA or diffB is less than a fraction of A, and if not, calculate output sample values from the non-sinusoidal model.

8. The graphics processing system of claim 7 wherein the fraction of A is one-fourth.

9. The graphics processing system of claim 7 wherein the resampling circuit comprises a resampling circuit operable to estimate A from:


A≈s+c/2if(s>c),

otherwise A≈c+s/2,
where s=|A SIN| and c=|A COS|.

10. The graphics processing system of claim 1 wherein the resampling circuit comprises a resampling circuit operable to calculate the sine-wave model and the output sample values therefrom from the equation:




  
    
      f
      ⁡
      
        (
        
          Δ
          ⁢
          
            
          
          ⁢
          p
        
        )
      
    
    =
    
      
        ∑
        
          i
          =
          0
        
        3
      
      ⁢
      
        
          
            C
            i
          
          ⁡
          
            (
            
              Δ
              ⁢
              
                
              
              ⁢
              p
            
            )
          
        
        i
      
    
  
  ,



where k=V1−V0, C3=gr1+gr0−2k, C2=k−C3−gr0, C1=gr0, C0=V0, and


grp=−A sin(φ)×ω sin(ωp)+A cos(φ)×ω cos(ωp),


where grp is the gradient value cosited at position p, ω is the angular frequency,


φ=arctan 2(A SIN ,A COS), and



A=√{square root over ((A SIN)2+(A COS)2)}{square root over ((A SIN)2+(A COS)2)},



where




  
    A
    ⁢
    
      
    
    ⁢
    
      cos
      ⁡
      
        (
        φ
        )
      
    
  
  =
  
    
      A
      ⁢
      
        
      
      ⁢
      COS
    
    =
    
      
        
          
            
              V
              1
            
            -
            
              V
              
                -
                1
              
            
          
          
            2
            ⁢
            
              
            
            ⁢
            
              sin
              ⁡
              
                (
                ω
                )
              
            
          
        
        ⁢
        
          
        
        ⁢
        and
        ⁢
        
          
        
        ⁢
        A
        ⁢
        
          
        
        ⁢
        
          sin
          ⁡
          
            (
            φ
            )
          
        
      
      =
      
        
          A
          ⁢
          
            
          
          ⁢
          SIN
        
        =
        
          
            
              
                V
                1
              
              +
              
                V
                
                  -
                  1
                
              
              -
              
                2
                ⁢
                
                  V
                  0
                
              
            
            
              2
              ⁢
              
                (
                
                  
                    cos
                    ⁡
                    
                      (
                      ω
                      )
                    
                  
                  -
                  1
                
                )
              
            
          
          .

11. A graphics processing system, comprising:
a graphics processor operable to generate data representing graphics primitives;
a triangle engine coupled to the graphics processor and operable to render the graphics primitives;
a pixel engine coupled to the triangle engine and operable to generate data representing pixels of an image; and
a resampling engine coupled to the pixel engine and operable to calculate output sample values from input sample values corresponding to graphics data of a source image, the resampling engine comprising:
a first resampling stage operable to calculate from a sample set of the input sample values an angular frequency value w for a sine-wave model and determine whether the frequency value w is in a frequency range, in response to the frequency value being in the frequency range, the first resampling stage operable to determine from the sample set a sinusoidal model from which the output sample values are calculated and calculate the output sample values from the sinusoidal model; and
a second resampling stage coupled to the first resampling stage, in response to the frequency value not being in the frequency range, the second resampling stage operable to determine from the sample set a non-sinusoidal model from which output sample values are calculated and calculate the output sample values from the non-sinusoidal model.

12. A computer system, comprising:
a processor having a processor bus;
an input device coupled to the processor through the processor bus adapted to allow data to be entered into the computer system;
an output device coupled to the processor through the processor bus adapted to allow data to be output from the computer system;
an interface circuit coupled to the processor and the input and output devices;
a memory coupled to the processor through the interface circuit and adapted to store data; and
a graphics processing system coupled to the interface circuit and adapted to generate and process graphics data, the graphics processing system comprising:
a graphics processor operable to generate data representing graphics primitives;
a triangle engine coupled to the graphics processor and operable to render the graphics primitives;
a pixel engine coupled to the triangle engine and operable to generate data representing pixels of an image; and
a resampling circuit coupled to the pixel engine to provide resample output values, the resampling circuit operable to calculate from a sample set of the sample values an angular frequency value w for a sine-wave model and determine whether the frequency value w is in a frequency range, where the frequency value w is in the frequency range, the resampling circuit is further operable to determine from the sample set a sinusoidal model from which the resample output values are calculated, where the angular frequency value w is not in the frequency range, the resampling circuit is operable to determine from the sample set a non-sinusoidal model from which the resample output sample values are calculated and calculate resample output sample values from the resulting model.

13. The computer system of claim 12 wherein the resampling circuit of the graphics processing system comprises a resampling circuit operable to determine from the sample set a cubic transition model between two of the input samples from which the resample output sample values are calculated when the angular frequency value w is not in the frequency range.

14. The computer system of claim 12 wherein the resampling circuit of the graphics processing system comprises a resampling circuit operable to determine whether the frequency value w is in a frequency range between arccos(−0.95)≦ω<arccos(0.9).

15. The computer system of claim 12 wherein the resampling circuit of the graphics processing system comprises a resampling circuit operable to calculate the angular frequency value w from a sample set of sample values including first, second, third, fourth, and fifth input sample values and the angular frequency value w is calculated from:




  
    ω
    =
    
      arc
      ⁢
      
        
      
      ⁢
      
        cos
        (
        
          
            (
            
              
                
                  d
                  1
                
                
                  d
                  2
                
              
              -
              1
            
            )
          
          2
        
        )
      
    
  
  ,



where d1=(V−1−V2) and d2=(V0−V1) if |V0−V1|>|V−1−V0|,
otherwise d1=(V−2−V1) and d2=(V−1−V0),
where V−2, V−1, V0, V1, and V2, are the first, second, third, fourth, and fifth input sample values, respectively.

16. The computer system of claim 15 wherein the resampling circuit of the graphics processing system comprises a resampling circuit operable to calculate the sine- wave model and the output sample values therefrom from the equation:


Vp=A sin(ωp+φ)+B,

where Vp is the output sample value at position p, ω is an angular frequency calculated from the input sample values,


B=V0−A SIN,



φ=arctan 2(A SIN ,A COS), and



A=√{square root over ((A SIN)2+(A COS)2)}{square root over ((A SIN)2+(A COS)2)},


where





  
    A
    ⁢
    
      
    
    ⁢
    COS
  
  =
  
    
      
        
          
            V
            1
          
          -
          
            V
            
              -
              1
            
          
        
        
          2
          ⁢
          
            
          
          ⁢
          
            sin
            ⁡
            
              (
              ω
              )
            
          
        
      
      ⁢
      
        
      
      ⁢
      and
      ⁢
      
        
      
      ⁢
      A
      ⁢
      
        
      
      ⁢
      SIN
    
    =
    
      
        
          
            V
            1
          
          +
          
            V
            
              -
              1
            
          
          -
          
            2
            ⁢
            
              V
              0
            
          
        
        
          2
          ⁢
          
            (
            
              
                cos
                ⁡
                
                  (
                  ω
                  )
                
              
              -
              1
            
            )
          
        
      
      .

17. The computer system of claim 15 wherein the resampling circuit of the graphics processing system comprises a resampling circuit operable to calculate the sine-wave model and the output sample values therefrom from the equation:


Rp=A sin(φ)cos(ωp)+A cos(φ)sin(ωp)+B,

where Rp is the output sample value at position p, ω is the angular frequency,


B=V0−A SIN,



φ=arctan 2(A SIN ,A COS), and



A=√{square root over ((A SIN)2+(A COS)2)}{square root over ((A SIN)2+(A COS)2)},


where





  
    A
    ⁢
    
      
    
    ⁢
    
      cos
      ⁡
      
        (
        φ
        )
      
    
  
  =
  
    
      A
      ⁢
      
        
      
      ⁢
      COS
    
    =
    
      
        
          
            
              V
              1
            
            -
            
              V
              
                -
                1
              
            
          
          
            2
            ⁢
            
              
            
            ⁢
            
              sin
              ⁡
              
                (
                ω
                )
              
            
          
        
        ⁢
        
          
        
        ⁢
        and
        ⁢
        
          
        
        ⁢
        A
        ⁢
        
          
        
        ⁢
        
          sin
          ⁡
          
            (
            φ
            )
          
        
      
      =
      
        
          A
          ⁢
          
            
          
          ⁢
          SIN
        
        =
        
          
            
              
                V
                1
              
              +
              
                V
                
                  -
                  1
                
              
              -
              
                2
                ⁢
                
                  V
                  0
                
              
            
            
              2
              ⁢
              
                (
                
                  
                    cos
                    ⁡
                    
                      (
                      ω
                      )
                    
                  
                  -
                  1
                
                )
              
            
          
          .

18. The computer system of claim 17 wherein the resampling circuit of the graphics processing system comprises a resampling circuit further operable to verify the accuracy of the sine-wave model by calculating:


diffA=|R−2−V−2| and diffB=|R2−V2|,

the resampling circuit further operable to confirm that diffA or diffB is less than a fraction of A, and if not, calculate output sample values from the non-sinusoidal model.

19. The computer system of claim 18 wherein the fraction of A is one-fourth.

20. The computer system of claim 18 wherein the resampling circuit of the graphics processing system comprises a resampling circuit operable to estimate A from:


A≈s+c/2if(s>c),

otherwise A≈c+s/2,
where s=|A SIN| and c=|A COS|.

21. The computer system of claim 12 wherein the resampling circuit of the graphics processing system comprises a resampling circuit operable to calculate the sine-wave model and the output sample values therefrom from the equation:




  
    
      f
      ⁡
      
        (
        
          Δ
          ⁢
          
            
          
          ⁢
          p
        
        )
      
    
    =
    
      
        ∑
        
          i
          =
          0
        
        3
      
      ⁢
      
        
          
            C
            i
          
          ⁡
          
            (
            
              Δ
              ⁢
              
                
              
              ⁢
              p
            
            )
          
        
        i
      
    
  
  ,



where k=V1−V0, C3=gr1+gr0−2k, C2=k−C3−gr0, C1=gr0, C0=V0, and


grp=−A sin(φ)×ω sin(ωp)+A cos(φ)×ω cos(ωp),


where grp is the gradient value cosited at position p, ω is the angular frequency,


φ=arctan 2(A SIN ,A COS), and



A=√{square root over ((A SIN)2+(A COS)2)}{square root over ((A SIN)2+(A COS)2)},



where




  
    A
    ⁢
    
      
    
    ⁢
    
      cos
      ⁡
      
        (
        φ
        )
      
    
  
  =
  
    
      A
      ⁢
      
        
      
      ⁢
      COS
    
    =
    
      
        
          
            
              V
              1
            
            -
            
              V
              
                -
                1
              
            
          
          
            2
            ⁢
            
              
            
            ⁢
            
              sin
              ⁡
              
                (
                ω
                )
              
            
          
        
        ⁢
        
          
        
        ⁢
        and
        ⁢
        
          
        
        ⁢
        A
        ⁢
        
          
        
        ⁢
        
          sin
          ⁡
          
            (
            φ
            )
          
        
      
      =
      
        
          A
          ⁢
          
            
          
          ⁢
          SIN
        
        =
        
          
            
              
                V
                1
              
              +
              
                V
                
                  -
                  1
                
              
              -
              
                2
                ⁢
                
                  V
                  0
                
              
            
            
              2
              ⁢
              
                (
                
                  
                    cos
                    ⁡
                    
                      (
                      ω
                      )
                    
                  
                  -
                  1
                
                )
              
            
          
          .

22. A computer system, comprising:
a processor having a processor bus;
an input device coupled to the processor through the processor bus adapted to allow data to be entered into the computer system;
an output device coupled to the processor through the processor bus adapted to allow data to be output from the computer system;
an interface circuit coupled to the processor and the input and output devices;
a memory coupled to the processor through the interface circuit and adapted to store data; and
a graphics processing system coupled to the interface circuit and adapted to generate and process graphics data, the graphics processing system comprising:
a graphics processor operable to generate data representing graphics primitives;
a triangle engine coupled to the graphics processor and operable to render the graphics primitives;
a pixel engine coupled to the triangle engine and operable to generate data representing pixels of an image; and
a resampling engine coupled to the pixel engine and operable to calculate output sample values from input sample values corresponding to graphics data of a source image, the resampling engine comprising:
a first resampling stage operable to calculate from a sample set of the input sample values an angular frequency value w for a sine-wave model and determine whether the frequency value w is in a frequency range, in response to the frequency value being in the frequency range, the first resampling stage operable to determine from the sample set a sinusoidal model from which the output sample values are calculated and calculate the output sample values from the sinusoidal model; and
a second resampling stage coupled to the first resampling stage, in response to the frequency value not being in the frequency range, the second resampling stage operable to determine from the sample set a non-sinusoidal model from which output sample values are calculated and calculate the output sample values from the non-sinusoidal model.