US 20040207631 A1 Abstract A method for generating bump map data substantially in real time for use in a 3-dimensional computer graphics system. Data is received which defines an area to which a texture is to be applied. Texture data to apply to the area is also received. This data includes surface height data. A set of partially overlapping samples of texture data are then filtered and surface tangent vectors derived therefrom. A bump map surface normal is then derived from the tangent vectors.
Claims(15) 1. A method for generating bump map data substantially in real time for use in a 3-dimensional computer graphics system comprising the steps of:
receiving data defining an area to which a texture is to be applied; receiving texture data to apply to the area, the data including surface height data; filtering each of a set of partially overlapping samples of the texture data; deriving surface tangent vectors from the filtered samples; and deriving a bump map surface normal from the surface tangent vectors. 2. A method according to 3. A method according to 4. A method according to 5. A method according to 6. Apparatus for generating bump map data substantially in real time for use in a 3-dimensional computer graphics system comprising:
means for receiving data defining an area to which a texture is to be applied; means for receiving texture data to apply to the area, the data including height data; means for filtering each of a set of partially overlapping samples of the texture data; means for deriving surface tangent vectors from the filtered samples; and means for deriving a bump map surface normal from the surface tangent vectors. 7. Apparatus according to 8. Apparatus according to 9. Apparatus according to 10. Apparatus according to 11. A 3D graphics system comprising a plurality of colour data processing means for generating data for use in shading an image to be represented by the 3D graphics system;
means for supplying texture data to be applied to the image; and means for assigning the colour data processing means to the generation of bump map data for use in applying the texture data to the image. 12. (Cancelled). 13. (Cancelled). 14. A method for generating bump map data for use in a 3-dimensional computer graphics system comprising the steps of:
receiving data defining an area to which a texture is to be applied; receiving texture data to apply to the area, the data including surface height data; filtering each of a set of partially overlapping samples of the texture data; deriving surface tangent vectors from the filtered samples; and deriving a bump map surface normal from the surface tangent vectors. 15. Apparatus for generating bump map data for use in a 3-dimensional computer graphics system comprising:
means for receiving data defining an area to which a texture is to be applied; means for receiving texture data to apply to the area, the data including height data; means for filtering each of a set of partially overlapping samples of the texture data; means for deriving surface tangent vectors from the filtered samples; and means for deriving a bump map surface normal from the surface tangent vectors. Description [0001] This invention relates to a method and apparatus for generating bump map data for use in a 3 dimensional computer graphics system. [0002] In the field of 3D computer graphics, detail is often added to otherwise smooth objects though the use of Bump Mapping, which was introduced by Blinn in his paper [0003] In Blinn's technique, each perturbation is computed by first taking derivates of a bump displacement texture or ‘height map’ and subsequently applying it to the original surface normal and surface tangent vectors. The height map is a simple array of scalar values that gives the ‘vertical’ displacement or ‘height’ of a surface at regular grid points relative to that surface. Typically these are represented by monochromatic image data, e.g. a bitmap, with the brightness of any pixel being representative of the ‘height’ at that point. Standard texture mapping practices are used to access the height data. The normal perturbations and lighting calculations are done in global or model space. FIG. 8 [0004] A more ‘hardware friendly’ method was later developed by Peercy et al (“Efficient Bump Mapping Hardware”, SIGGRAPH 1997, pp 303-306, (also U.S. Pat. No. 5,949,424)). This directly stores perturbed surface normals in a texture map, often called a normal map. Unlike Blinn's method, these normals are defined in a local tangential coordinate space, which can be likened to the representation of parts of the earth's surface on a page in an atlas. In Peercy's technique, the lights used for shading are also transformed into this tangential space and thus the shading calculations are also computed locally. This process significantly reduces the number of calculations required when using bump mapping. It has become popular in recent 3D hardware systems and is sometimes known as ‘Dot3 bump mapping’. [0005] To minimize the texture memory and, more importantly, memory bandwidth required by this procedure, it is desirable to compress the normal maps. Unfortunately many of the commonly used texture compression schemes are not suitable as they cause a loss of information that, when applied to the special case of normal maps, can cause an unacceptable degradation in image quality. Two methods that are specifically tailored to normal maps, however, are described in our International patent application No. WO9909523—these typically still use 16 bits to represent each surface normal. [0006] This then leaves the task of generating the normal map. One popular method again uses an initial height map, as originally described by Blinn. From that height map, a normal map can then be pre-computed, prior to rendering, by taking the cross product of the local derivative vectors of the height function sampled at regular positions. For cases where texture filtering is required, e.g. those based on the well-known MIP mapping techniques, the height map should be repeatedly down-sampled and the associated normal map regenerated to produce the multiple MIP map levels. Problems can arise, however, when applying the texture filtering techniques, e.g. bilinear or trilinear filtering, to normal maps. [0007] It should be noted that whereas the height map could be of relatively low precision—often as low as, say 8 or even 4 bits per pixel—the normal map may require 16 to 32 bits per pixel. The pre-processing steps of generating and compressing the normal map and the process of using the compressed normal map in 3D rendering are shown in FIG. 2. In the generation phase a height map [0008] Also well known in the art is the aspect of texture filtering, primarily the application of bilinear or trilinear filtering, the latter as invented by Williams (“Pyramidal Parametrics”, Lance Williams, Computer Graphics, Vol. 7, No. 3, July 1983, pp 1-11). Bilinear filtering is briefly discussed below, since trilinear filtering is just the blending of two bilinear operations. [0009] A 2D texture can be considered to be a vector function of 2 variables (U, V). For simplicity in this discussion, we will assume that, for an N×N pixel texture, the values of U and V range from 0 to N. When bilinear filtering is applied, the pixel, or “texel”, values stored in the texture can be considered to be representative of the points in the centres of the respective texels, i.e. at coordinates (i+0.5, j+0.5), where i and j are integers and represent the texel coordinate of the particular texel. This is illustrated in FIG. 3 for texel (i,j), the centre of which is indicated by ‘20’. At this point in the texture, bilinear filtering will return the colour of that texel. Similarly, sampling at locations ‘21’, ‘22’, and ‘23’ will return the colours of texels (i+1,j),(i,j+1), and (i+l ,j+1) respectively. Now consider any sampling location within the square formed by ‘20’, ‘21’, ‘22’, and ‘23’, such as point ‘24’. Such a point has texture coordinates (u [0010] In particular, the process used in the art will be some simple variation of the following: [0011] U′:=Us−0.5;//Place stored texel value at centre of texel [0012] V′:=Vs−0.5; [0013] Ui:=floor(U′); [0014] Vi:=floor(V′); [0015] Ublend:=U′−Ui; [0016] Vblend:=V′−Vi; [0017] //Do 2 horizontal linear blends [0018] Colour0:=LinearBlend(Texel(Ui, Vi), Texel(Ui+1, Vi), Ublend); [0019] Colour1:=LinearBlend(Texel(Ui, Vi+1), Texel(Ui+1, Vi+1), Ublend); [0020] //Do 1 vertical linear blend [0021] Result:=LinearBlend(Colour0, Colour1, Vblend); [0022] The Ublend and Vblend values are thus in the range [0 . . . 1), and can be most conveniently represented by a fixed point number of, say, 8 to 16 bits precision. [0023]FIG. 4 [0024] It should be noted that the colours in 3D computer graphics are usually 4-D entities, having Red, Green, Blue, and Alpha (i.e. transparency) components. When the bilinear blending described above is performed, all four components of the various colour values are operated on in parallel. This is shown in the second stage of the bilinear operation in FIG. 4 [0025] Another known aspect of 3D computer graphics is that of fitting smooth surfaces through or near a set of control points. In particular we are interested in two types known as uniform B-spline and Bezier splines, as described in literature such as “Computer Graphics. Principles and Practice” (Foley et al) or “Curves and Surfaces for CAGD. A practical guide” (Farin). [0026] Of particular interest to this application is the case of a bi-quadratic B-spline which has C1 continuity (i.e. continuous first derivative). A bi-quadratic B-spline also has the property that, for any point on the surface, a sub-grid of 3×3 control points is needed to evaluate that point and/or derivatives at that point. A one-dimensional slice though a section of a quadratic B-spline is shown in FIG. 5. The points, ‘80’, ‘81’, and ‘82’ can be considered to be three adjacent control points in a row of the grid. The region of the curve between ‘85’, and ‘86’ depends only on these three control values (and the neighbouring 6 values in the 3×3 sub-grid in the case of a bi-quadratic surface). [0027] One popular way of evaluating such a curve is to first convert it to the equivalent Bezier representation, i.e. a different set of 3 control points, and then apply the de Casteljau algorithm which uses repeated linear interpolation (see Farin). For the simple case of quadratic curves, this amounts to using a new set of control points which are ‘88’, ‘81’ (i.e., it is re-used), and ‘89’. Points ‘88’ and ‘89’ are just the mid points of the connecting line segments and could be found by simple averaging. [0028] For the conversion of a bi-quadratic B-spline surface, the 3×3 grid of control points can be replaced by an equivalent set of 3×3 Bezier control points. An example showing the situation for a bi-quadratic surface is shown, in plan form, in FIG. 6. The original 9 B-spline control points, one example of which is shown by ‘100’, are converted into the equivalent Bezier control points, such as ‘101’. Stating this more precisely, if the grid of 3×3 B-spline points are:
[0029] then the equivalent set of Bezier points are computed from:
[0030] Referring again to FIG. 6, the region of interest is the central ‘square’, i.e. a position specified by (u,v), where 0≦u,v≦1. In the case of a bi-quadratic surface, one method based on de Casteljau would be to bi-linearly interpolate sets of 2×2 neighbouring control points, using (u,v) as weights, to produce a new set of 2×2 intermediate control points. One of the four sets of 2×2 intermediate control points is indicated by ‘102’. These four results are, in turn, bilinearly interpolated, again using the (u,v) weights, to produce the surface point. If tangents to the surface are required, a method such as given by Mann and Rose (“Computing values and derivatives of Bezier and B-spline tensor products”, CAGD, Vol 12, February 1995) can be used. For the bi-quadratic case, this can be done by performing additional linear interpolations using the 2×2 intermediate control values. Finally, taking the cross product of these tangents generates the surface normal. [0031] The height map defines ‘height’ values only at certain sample locations and so a means of computing the height surface at other points is required. In particular, bump mapping requires the surface normal which, in turn, usually implies the need for surface tangents. Blinn points out that the surface height is not actually required and proposes a function that only computes tangents. He notes that in order to avoid discontinuities in the shading, his tangent functions are continuous. Using the 3×3 grid of height samples shown in 6, Blinn's function performs 3 bilinear blends respectively of the top left, top right, and bottom left neighbours, and then computes the differences of the top left and top right result and the top left and bottom left result as part of the tangent generation. [0032] Although Blinn's function results in a continuous normal, its derivative can have discontinuities. Unfortunately, the human visual system is very sensitive to changes in the derivative of image intensity, and so ‘artefacts’ can be seen. The method also tends to emphasise the underlying grid of the height map, which can be seen in FIG. 8 [0033] Although the introduction of Peercy et al's pre-perturbed normal map method makes bump mapping more practical in real-time hardware, it still requires ‘large’ texture formats as well as the separate pre-processing step to convert a height map to normal map. The ‘large’ texture formats consume valuable bandwidth as well as memory and cache storage and, although special normal map compression techniques exist, these formats are still often larger than the original source height map. Also filtering of the normal map may also be problematic. [0034] A further limitation of Peercy et al's technique is that dynamic bump mapping, i.e. where the bump heights are computed frame-by-frame, is far more difficult to achieve. For example, the height values may be generated as the result of a separate rendering pass. The pre-processing step, including generation of the various MIP map levels, may take too much time to allow real-time rendering. [0035] Finally, it is beneficial to use a height function with C2 (or higher) continuity so that the normal interpolation is C1 (or higher). In particular, it is important to have an inexpensive means of producing this function. [0036] We have appreciated that it is possible to implement, in hardware, an additional set of functions that provides an efficient means for direct transformation of a height map into filtered perturbed surface normals that have C1 continuity. These normals can subsequently be used for various rendering purposes such as per-pixel lighting. In particular, we have devised a method which, by re-using colour texture filtering hardware that is ubiquitous in today's graphics systems in a new way with the addition of some small processing units, achieves the functions needed to compute the normal from a smooth surface controlled by a set of heights. Thus the data can be generated substantially in real time. [0037] The filtered surface normals are created ‘on demand’ and are not stored. This provides the joint benefits of reducing the amount of texture data and bandwidth needed for bump mapping, as well as overcoming some of the issues with the filtering of normal maps. This feature is also important when using dynamic height maps in real-time rendering since a pre-processing step may be prohibitive. [0038] Embodiments of the invention keep the advantages of computing bump map-based shading in local tangent space as described by Peercy et al, (although it is not restricted to doing so), with the convenience of directly using Blinn's height map but with the option of using a function with higher continuity. [0039] Preferred embodiments of the invention will now be described in detail by way of example with reference to the accompanying diagrams in which: [0040]FIG. 1 shows the process of perturbing surface Normals as described by Blinn; [0041]FIG. 2 shows a flow chart of the pipeline used for Peercy et al's method described above; [0042]FIG. 3 shows the relationship of bilinear filtering of a texture to the texels of that texture; [0043]FIG. 4 [0044]FIG. 4 [0045]FIG. 5 shows a segment of a piecewise quadratic B-spline curve, or equivalently, a section through a bi-quadratic B-spline Surface; [0046]FIG. 6 shows a plan view of a section of a height map being interpreted as a bi-quadratic B-spline surface; [0047]FIG. 7 shows an overview of a hardware system embodying the invention with modifications to support normal generation from height maps; [0048]FIG. 8 [0049]FIG. 9 shows some alternative filter patterns that could be used for computing derivatives of a bump map. [0050] The preferred embodiment will now be described. Access is provided to height map textures, which store an array of height values, using preferably 4 or 8 bits per texel. Each value will encode a fixed-point number with some number of fraction bits—preferably ¼ of the bits will be assigned to the fractional part. [0051] The embodiment fits a bi-quadratic B-spline through this set of points, thus giving the virtual height map texture C1 continuity (i.e. continuous first derivatives). In FIG. 5, the points ‘80’, ‘81’, and ‘82’ represent three adjacent height values/control points in a row of the height map. The value in the texture is allocated (preferably) to the y dimension while the other coordinate values (i.e., x and z) are implicitly defined by the texel's coordinate position. Alternative embodiments may assign these dimensions in some other permutation. The other control points needed for the surface in the texel region are shown, in plan form, in FIG. 6. [0052] The manner in which the normal is computed is now described with reference to FIG. 7. As with the standard texture filtering system described above, i.e. FIGS. 4 [0053] A modified texel fetch unit, ‘151’, which in FIG. 4 [0054] For brevity, these have be renumbered a, b, etc. [0055] It will be apparent to those skilled in the art that, with application of the address-bit interleaved texture storage format described in our British patent number GB2297886, such a height-map can be packed into the ‘equivalent’, in terms of storage, of a colour texture of ½×½ resolution of the height map. Each 2×2 group of scalar height data would occupy the space of a single four-dimensional colour. With such a format, the height map data can then be accessed using a very simple modification of exactly the same fetch mechanism used by units ‘60’ thru ‘63’ in FIG. 4 [0056] The 3×3 grid of samples is then fed into the ‘Replicate’ unit, [0057] . . . while similarly the green channel receives the top right set, the blue, the bottom left, and the alpha receives the bottom right. Clearly some values, such as b′ or e′, will be used more than once, thus the grids supplied to each unit overlap at least partially. [0058] Unit [0059] As Ublend and Vblend are typically fixed point numbers, it should be appreciated that these ‘calculations’ are completely trivial and incur no cost at all in hardware. [0060] These new blend values are distributed to bilinear units, ‘65’ thru ‘68’ as follows:
[0061] This manipulation of the blend factors eliminates the need to convert from the quadratic B-spline control points to the Bezier control points, as described previously in Equation 1. These bilinear units therefore effectively produce data which will enable surface normals with C1 continuity to subsequently be derived. [0062] The results of the 4 bilinear interpolations are fed to the tangent construction unit, [0063] Tang1[X]:=1 Texturesize; [0064] Tang1[Y]:=LinearInterpolate(VBlend, [0065] GreenResult−RedResult, [0066] AlphaResult−BlueResult); [0067] Tang1[Z]:=0; [0068] Tang2[X]:=0; [0069] Tang2[Y]:=LinearInterpolate(UBlend, [0070] BlueResult−RedResult, [0071] AlphaResult−GreenResult); [0072] Tang2[Z]:=1 Texturesize [0073] where [0074] LinearInterpolate (x, A, B):=A+x*(B−A); [0075] For reasons that will soon be apparent, unit [0076] Finally, in unit ‘156’, the cross product of these tangents is computed. It should be noted that if the preferred embodiment is chosen, the presence of zeros in the tangent components simplifies the cross product to the following calculation: [0077] N[x]:=Tang1[y]; [0078] N[y]:=1 Texturesize [0079] N[z]:=Tang2[y]; [0080] This vector is then normalised, preferably by squaring the N vector, computing the inverse of the square root of the result, and multiplying that scalar by the original components. The normalisation step may appear expensive, but it would be a requirement of any system that supported compressed normal maps, such as that described in WO9909523 or British patent application No. 0216668.4. Thus, if such texture formats were already supported, the re-normalisation hardware would be reused. An example of the output of this embodiment is shown in FIG. 8 [0081] In an alternative embodiment, trilinear filtering can be adapted to support normal generation—the only difference in procedure will be that the values fed into tangent construction unit ‘155’ will be a ‘blend’ of the values computed from two adjacent MIP map levels chosen. Other embodiments supporting improved anisotropic filtering are also feasible. [0082] In another embodiment, an interpolated scale factor may be applied to the deltas/tangents before normalisation so that a height map can be scaled differently for different models or different parts of the same model. [0083] In another embodiment, the blend factor adjust unit, [0084] In another embodiment, the actual interpolated height value would be computed by including a third linear blending operation. [0085] In another embodiment, Blinn's height interpolation function could be employed. In this embodiment, the blend factor adjust unit, [0086] In another alternative embodiment, dedicated sampling hardware could be included that takes numerous texture samples and applies an alternative derivative filter such as 4 taps, Sobel, Prewitt, Parks-McClellan derivatives filters as represented in FIG. 9. [0087] In another embodiment, colour textures are also filtered using bi-quadratic B-splines, either through the addition of bilinear filtering units, or by iterations through the colour channels, whereby the individual weights to the bilinear units are adjusted according to the previously described embodiments. Referenced by
Classifications
Legal Events
Rotate |