CURVE GENERATION IN A DISPLAY SYSTEM

Abstract: A curve generator for a display system comprises arc generation means for generating an arc 10 of a circle 14 from data defining the locations of two end points P₁ and P₃ and an intermediate point P₂ on the arc. The arc generation means comprises initialisation means 40 for calculating the angle subtended between a first vector 21, from a first of the end points P₁ to the intermediate point P₂, and a second vector 32, from the second of the end points P₃ to the intermediate point P₂, and arc plotting means 44 for defining a succession of further vectors n1 from said first end point P₁ and for calculating, for each further vector, its point of intersection P n with a counterpart vector 3n, from said second end point P₃, with which it subtends said angle, whereby a succession of further points P n are plotted on the circular arc. The plotting logic thus plots the points of the arc with respect to a given point on the arc itself by generating vectors from that given point and enables the computation of the arc to be performed substantially within the system co-ordinate space in which the arc exists, which reduces the number of places of accuracy needed in order to accurately compute the arc.

Description: The present invention relates to a curve generator for a display system, to a display system including such a generator and to methods of generating curves in a display system. To generate curve such as an arc of a circle from data defining the arc is a non-trivial task in a computer due to the significant calculations involved. A circle centred at the origin of a coordinate system can be defined by the following, well known, equation: x² + y² = R² (1) where x and y are the variable horizontal and vertical coordinates of the coordinate system, and R is the radius of the circle. Solving this equation for y gives: y = ± (R² - x²) ½ (2) A simple approach to drawing one quadrant of a circle using this equation is to increment x in unit steps from O to R solving for + y at each of the steps. The other three quadrants of the circle can then be determined about the origin using the symmetry of the circle. Although this technique works, it is inefficient because of the multiply and square root operations. A similarly inefficient method is simply to plot R cos ϑ or R sin ϑ by stepping ϑ from 0 to 90°, and then generating the other three quadrants using the symmetry of the circle. An improved method of generating a circular arc is described in J E Bresenham's article "A Linear Algorithm for Incremental Digital Display of Circular Arcs" published on pages 100-106 of "Communications of the ACM", Vol 20, No 2, in February 1977. This method, which is based on equation (1) above, was conceived for use with pen plotters, although it is applicable generally to pixel based display systems. In accordance with the method, points on a circle centred at the origin are generated by stepping round the circle. At each step the pixel point which is closest to the true circle is selected for display by employing an error term: D (Pt i ) = (x² i + y² i ) - R² (3) where D(pt i ) is the difference between the true circle for the ith point (Pt i ) and x i and y i are the x and y values calculated for the ith point. References to these and to many other techniques can be found in Foley and Van Dam's book "Fundamentals of Interactive Computer Graphics" (published in 1982 by Addison-Wesley Publishing Company), on pages 442-446. In theory the prior techniques work well, although with varying degrees of efficiency. They are designed essentially to draw complete circles, although they can of course be used to draw an arc of a circle. In practice, however, when implemented in a graphics processing system, severe limitations as to their applicability to the drawing of arcs are encountered. In particular, the prior techniques run into difficulty in accurately plotting an arc which forms part of a very large circle where only part of the circle is within the system coordinate space. This results primarily from the need to compute, either explicitly or implicitly, the actual centre of the circle of which the arc forms part. If an arc is to be drawn which is almost a straight line, it will be appreciated that a very large number will be needed to define the radius of circle of which it forms part, as the centre of the circle may be well outside the bounds of the coordinate space within which the arcs are to be generated. In accordance with a first aspect of the present invention there is provided, a curve generator for a display system, the curve generator comprising arc generation means for generating an arc of a circle from data defining the locations of two end points and an intermediate point on the arc, wherein the arc generation means comprises initialisation means for calculating the angle subtended between a first vector, from a first of the end points to the intermediate point, and a second vector, from the second of the end points to the intermediate point, and arc plotting means for defining a succession of further vectors from said first end point and for calculating, for each further vector, its point of intersection with a counterpart vector, from said second end point, with which it subtends said angle, whereby a succession of further points of intersection are plotted on the circular arc. In accordance with a second aspect of the present invention there is provided, a method of generating an arc of a circle in a display system comprising processing and storage means from data defining the locations of two end points and an intermediate point on the arc comprising the steps of: (a) calculating the angle subtended between a first vector, from a first of the end points to the intermediate point, and a second vector, from the second of the end points to the intermediate point; (b) storing the calculated angle in said storage means; (c) defining a further vector from said first end point; (d) calculating, for said further vector, its point of intersection with a counterpart vector, from the second end point, with which it subtends said calculated angle; and (e) repeating steps (c) and (d) for yet further vectors from said first end point whereby a succession of further points of intersection are plotted on the circular arc. The present invention is based on the well known theorem in geometry which states that for a triangle with vertices P₁, P₂ P₃ inscribed on a circle radius R: where a, b and c are the interior angles subtended at the vertices P₁, P₂ and P₃ respectively and A, B, and C are the lengths of the sides opposite to the vertices P₁, P₂ and P₃ respectively. Although this theorem is featured in good school textbooks in geometry, its application to solving the problem of arc generation in computers has not heretofore been recognised despite the extensive efforts in the art to improve the performance of arc generators. The advantages of the present invention result primarily from the fact that the computation of the points of the arc is not performed with respect to the centre of the circle of which the arc forms part. The points on the arc are plotted instead with respect to a given point on the arc itself by generating vectors from that given point. The present invention thus enables the computation of the arc to be performed substantially within the system coordinate space in which the arc exists, which reduces the number of places of accuracy needed in order to accurately compute the arc. In addition, the present invention has the advantage that the computation of the arc can be performed in integer arithmetic, which enhances the performance of the arc generator. In order that the present invention may be more fully appreciated, there follows a description of the principles behind the present invention and of particular embodiments of the invention with reference to the accompanying drawings in which: Figure 1 is an illustration, used to explain the principles of operation of the present invention, of an arc to be drawn within a rectangular coordinate space; Figure 2 is a further illustration, used to explain the principles of operation of the present invention; Figure 3 is a schematic block diagram showing the logical structure of a particular embodiment of the present invention; Figure 4 is an illustration to explain how a generalised curve may be drawn using the principles behind the present invention; Figure 5 is a schematic block diagram showing the logical structure of another particular embodiment of the present invention; and Figure 6 is a schematic diagram of a workstation which can incorporate an embodiment of the present invention such as shown in Figure 3 or 5. Before dealing with the embodiments of the present invention, the principles on which the present invention is based will be explained in the following. Let us assume that an arc is defined in three dimensional (x, y, z) coordinate space in terms of the two end points of the arc and a third point which lies on the arc intermediate the end points. Figure 1 illustrates an example of such an arc 10 within a coordinate space 12. For reasons of ease of illustration only, a two dimensional (x and y) coordinate space is shown, although it will be apparent from the following description that the treatment of the two-dimensional case can easily be expanded to three-dimensional (x, y, z) space. The arc 10 forms part of a circle 14 of radius R, the centre 16 of which lies outside the coordinate space 12 in the illustrated example. The arc 10 is defined in terms of the coordinates (x₁, y₁) of a first end point P₁ of the arc, the coordinates (x₃, y₃) of a second end point P₃ of the arc and the coordinates (x₂, y₂) of a third point P₂ which lies on the arc. The three points P₁, P₂ and P₃ can be considered to form respectively first, second and third vertices of a triangle 18. Figure 1 shows the interior angles subtended at the three vertices P₁, P₂ and P₃ to have the values a, b and c respectively and the sides 32, 31, 21 of the triangle to have the lengths A, B and C respectively. Equation (4) above applies to any triangle inscribed on a circle. Thus the following equation applies to triangle 20, shown in Figure 2, which comprises vertices P₁, P₃ and P n inscribed on the circle 14, the interior angles subtended at the vertices P₁, P₃ and P n being respectively a n , b n and c n and the lengths of the opposite sides 3n, 31, n1 being respectively A n , B n and C n : As the points P₁ and P₃ of the triangle 20 are common to the triangle 18, it follows that the length of the side 31 joining those points will have the length B for both triangles. Thus it also follows that for any triangle formed from points P₁ and P₃ and a third point P n on the arc, the interior angle b n is constant ie: b n = b = arcsin(B/2R) (6) Given that the sum of the interior angles of a triangle is Π radians, and that a straight line through a point subtends Π radians about that point, it can be shown by the following that the angle δ n subtended between a tangent 24 to the circle 14 at the point P₁ and the line n1 joining points P₁ and P n of the triangle 20 will be equal to the interior angle C n subtended at the vertex P₃ of that triangle 20. That is: Π = t + δ n + a n ; and (7) Π= b + a n + c n . (8) δ n = Π - t - (Π - b - c n ) = + b + c n - t (9) As P n → P₁ δ n → 0 and c n → 0. Thus t = b, both being constants, and δ n = c n (10) From equations (5) and (10) it can be seen that the length of the side n1 of the triangle 20 opposite the vertex P₃ is: If, then, an angle dir is defined as the angle subtended by the line n1 and the x co-ordinate axis (ie. the gradient of that line), the following equation can be used to compute dir: where g₃₁ is the angle of the line 31 to the x axis (ie. g₃₁ is the gradient of that line). Using equations (11) and (12), it is then possible to calculate the x n and y n coordinates of the point P n by evaluating the following expressions: Figure 3 is a schematic block diagram showing the logical structure of an arc generator forming a particular embodiment of the present invention. Only those parts of the arc generator which are necessary for explaining the present invention are shown in Figure 3. The arc generator would normally be incorporated in a display system (eg. a graphics workstation) of an otherwise conventional construction (see for example Figure 6). Initialisation logic 40 computes a number of initialisation values from the co-ordinate positions of the points P₁, P₂ and P₃ stored in input storage 38. The input storage can be part of the general purpose memory of a display system, or an input buffer or input registers, or, actually part of the arc generator. The co-ordinate positions can have been generated in a display system in response to positions indicated on the screen of a display system by mouse movements or have been generated in any other suitable way. As the initialisation values are calculated by the initialisation logic they are stored in the intermediate storage 42. The intermediate storage can be formed from dedicated registers, or can be configured in general purpose storage. The gradient or slope g₃₁ of the line 31 joining P₃ and P₁ is computed from the co-ordinate values for those points: g₃₁ = arctan ((y₃-y₁)/(x₃-x₁)) (15) The gradients (g₃₂ and g₂₁) for the lines 32 and 21 are similarly computed. From the gradients of the sides 32 and 21, stored in the intermediate storage 42, the internal angle b subtended at the point P₃ and the absolute value of the sine of that angle are then computed as: The length (B) of the line 31 joining points P₁ and P₃ is also computed as: The choice of sin or cos for the calculation of B depends on whether the angle g₃₁ is greater than, or less than π/4 in order to minimise errors. A stepping angle δ is also computed in order to determine the angle to be swept between successive points on the arc. This can be computed on the basis of δ = (π + b)/N where the number N is chosen to give a desired density of points for the arc. The number N can also be provided as an input from the input storage 38, or elsewhere as desired. It could also be calculated as a function of a desired resolution. Once the initialisation values have been computed and stored in the intermediate storage 42, the individual points on the arc can be plotted by plotting logic 44. The plotting logic defines a succession of further vectors from the first point in terms of a vector angle δ n = nδ where δ ₁ = δ for the first vector, δ n = nδ for the nth vector until δ N-1 =(N-1)δ for the N-1th vector. As is illustrated in Figure 2, the angle δ n is the angle between the tangent t to the circle and the vector in question. By incrementing n from 1 to N-1 the plotting logic sweeps out the arc from one of the points of the arc (eg. the first) rather than, as is done by the prior methods, from the centre of the circle. The x and y values of the end points of the arc P₁ and P₃ are given for n=o and n=N respectively. The plotting logic determines the intersection point of each of the vectors with a counterpart vector from the third point, with which is intersects at the angle b, stored in the intermediate storage. The plotting logic does this, not by defining a set of vectors from the third point, but by evaluating equations (13) and (14) for each of the vectors from the first point using the values stored in the intermediate storage. The end points of the arc are defined by the points P₁ and P₃. The equations (13) and (14) can be simplified for certain special cases as is explained below and the plotting logic contains specific logic for detecting these special cases. When the three points are nearly in a straight line, the angle B will approximate π. Thus, if on testing the value of B it is found to be within 2⁻¹² of π (the precise number depends on the required precision of the calculations), then the relationships b = Nδ, b = sin(b) and δ n = sinδ n can be used to simplify equation (13) to x n = x₁ + (nB)/(Ncos(dir)) (19) and equation (14) to y n = y₁ +(nB)/(Nsin(dir)). (20) If the angle B is found to be less than 2⁻¹⁶ radians (the precise number depends on the required precision of the calculations), then the two end points are nearly identical and the relationships 2R=B/sin(b)=m , where m is a constant, and m = (x₂ - x₁)/cos(g₂₁) or m = (y₂ - y₁)/sin(g₂₁) can be used to simplify equation (13) to x n = x₁ + (m|sin(δ n )|cos(dir)) (21) and equation (14) to y n = y₁ + (m|sinδ n |sin(dir)). (22) It should be noted that if B is found to be zero, then the two end points of the arc are identical and an arc cannot then be drawn as there is no orientation information. The calculation of the values above can be performed using integer arithmetic. Not only does this give a significant performance advantage over floating point arithmetic, but it also has the advantage that on overflow during binary integer arithmetic, the sign changes which occur mirror the sign changes for angle. This greatly simplifies the sign (+/-) management. The various trigonometric functions which are required can conveniently be provided using look-up tables in the plotting logic. The trig tables could, alternatively be provided in the control logic, which is responsible for the overall control of the arc generation, or elsewhere. Another advantage of the use of binary arithmetic is that the values resulting from evaluating the equations 13/14 etc. can be mapped onto individual display positions (pixels) by the plotting logic clipping the binary number at an appropriate accuracy. In the same way, intermediate points on the arc which fall outside the co-ordinate space may simply be discarded by clipping. The results of the calculations, (ie. the plots for the arc) are stored in results storage 46. The results storage shown in Figure 3 is the display buffer of the display system, and as such does not form part of the arc generator. The output of the plotting logic (ie. the clipped values) is (are) used to address the display buffer for storing data indicative of pixels to be displayed on a display screen. This need not be the case however, and the results storage could instead form part of general storage, or may even form part of the arc generator itself, in which case the output of the display plotting logic would be merely stored therein for further processing Although this particular embodiment of the present invention uses the equations 13/14 etc. to evaluate the appropriate intersection points, it will be understood that an alternative embodiment of the present invention could in fact find a point on the arc by defining a vector from the first point which is rotated from the vector 21 by a given amount, defining a vector from the third point which is rotated by the same amount from the vector 32, and then determining the intersection point of those vectors. A set of points on the arc could then be plotted by repeating this process for different degrees of rotation. This works because the interior angle at the point of intersection remain constant due once more to the fact that the interior angle of a triangle add up radians. The techniques described above can be used, not only to generate the arc of a circle, but also to generate any other curve for which a transform operator is known for mapping the curve onto a circular arc. For example, in the Graphics Object Content Architecture (GOCA) which forms a particular part of the Systems Network Architecture (SNA), for defining graphics objects, part of an ellipse in real, or reference space, is defined in elliptical space as the corresponding arc of a circle in terms of the two end points and an intermediate point on the arc. Figure 4 illustrates how an ellipse can be represented in this way. GOCA assumes real space to be a cube 60 with 2¹⁶ locations in each of the x, y, and z directions. In other words, 3 times 16 (ie. 48) bits are needed to identify a location in real space. In order to transform points between real and mapped space a transformation matrix (termed a PQRS matrix) is employed. Elliptical space is also defined as a cube, but as a result of the matrix calculation it has 2³³ locations on each side (2¹⁶ * 2¹⁶ + carry) . For reasons of ease of representation, only two-dimensional (x and y) real 60 and mapped 62 space is shown in Figure 4. Within the real space 60, a curve 64, which forms part of an ellipse, 65 is defined in terms of the two end points C₁ and C₃ and an intermediate point C₂ on the curve. In order to map the curve onto an arc of a circle, a PQRS matrix 68 is provided. The points C₁, C₂ and C₃ become the mapped points P₁, P₂ and P₃ in the mapped space after the inverse of the PQRS transformation. The points P₁ and P₃ form the two end points and the point P₂ forms the intermediate point on the on an arc 67 of a circle 68. If the techniques described above for plotting an arc of a circle are then used to plot the arc 68, and the resulting point plots in mapped space are transformed back into real space using the PQRS transformation, the individual points on the curve 64 can be plotted in real space. Clearly this technique can be employed for any general curve for which a transform is known which will allow the curve to be mapped onto an arc of a circle and vice versa. Figure 5 is a schematic block diagram showing the logical structure of a curve generator forming a second particular embodiment of the present invention. This curve generator is capable of generating a curve such as the curve 64. The curve generator comprises PQRS matrix storage for storing the transform operator information for transforming data between real and mapped space. In addition it comprises initial input storage 72 for the storage of the three points C₁, C₂ and C₃ for the ellipse to be drawn in real space. The coordinates of the points on the ellipse are transformed into mapped space by the inverse transform logic 74 and are stored in the input storage 38 of the arc generator as the mapped points P₁, P₂ and P₃. These mapped points are then processed by the initialisation logic 40 of the arc generator and are stored in the intermediate storage 42 as described above for the curve 10 with respect to Figure 3. The intermediate values stored in the intermediate storage are then processed as explained for the plotting logic 44 of the arc generator using the equations 13/14 etc. As shown in Figure 5, the results of the plotting logic (ie. the plots of the circular arc), are not passed directly to the results storage 46, but are instead passed to transform logic 76 which uses the inverse of the transform used be the transform logic 74 to transform the plots of the circular arc in mapped space into corresponding plots on the curve 64 for storage in the results storage 46. It will be appreciated that the clipping function described for the plotting logic in Figure 3 could be performed instead by the transform logic in the embodiment of Figure 5. It will also be appreciated that the plotting and transform logic could be incorporated into one unit as the output of one feeds directly into the other. Alternatively, the transform and inverse transform logic could be combined as their input parameters are similar, and appropriate control provided by the control logic 48. The advantages provided by the present invention are perhaps most apparent when the number of binary places of accuracy needed in the present invention is compared to that needed by prior approaches. In the example of the GOCA 33 bit co-ordinate space, it can be shown that the prior techniques can require up to 81 binary places of accuracy in order to accurately calculate all the arcs that need to be drawn. The present invention on the other hand allows curves to be plotted using at most 32 binary places of accuracy for the computation of the points on a circular arc, although 48 places are needed for the reverse transform into drawing, or real, space. Much of the computation can be performed in 16 places of accuracy. Mapping the resulting curve, or arc, points onto locations in real space (eg. onto pixel positions) can be simply done by truncating or clipping the calculated values at the 16th most significant bit. Figure 6 shows an overview of a workstation which can incorporate an embodiment of the present invention such as is shown in Figure 3 or Figure 5. The workstation comprises a number of different systems units connected via a system bus 82. The system bus comprises a data bus 84, an address bus 86 and a control bus 88. Connected to the system bus are a microprocessor 80, random access memory 90, a keyboard adapter 98, a display adapter 102, an I/O adapter 92, and a communications adapter 96. The keyboard adapter is used to connect a keyboard 100 to the system bus The display adapter connects the system bus to a display device 104. The I/O adapter likewise connects the system bus to other I/O devices such as disk units, and the communications adapter allows the workstation to be connected to and to communicate with an external processor or processors such as a host processor. A curve generator incorporating an arc generator in accordance with the present invention is implemented in software in the workstation shown in Figure 1. Control code for implementing the logic shown in Figures 3/5 is provided in the workstation storage 90 and the storage elements shown in Figures 3/5 are provided by configuring the workstation RAM. In the illustrated embodiment the display buffer is configured in RAM as well. A detailed listing of the code for implementing the logic and the storage elements shown in Figures 3/5 is not supplied with this description as the implementation of that logic is merely a matter of routine for the skilled person given the above description of the functions to be performed. Although particular embodiments and a particular implementation of the present invention have been described herein, it will be appreciated that many modifications and additions are within the scope of the appended claims. Instead of implementing the present invention in software, it will be appreciated, for example, that the present invention could equally be implemented with special purpose hardware logic, with or without the provision of special registers for the intermediate storage of variables. The logic units shown in Figures 3 and 5 could, for example, be implemented using programmable logic arrays. Also, if the display buffer were instead included in the display adapter, the arc and/or curve generators could be incorporated in the display adapter as well in order to relieve the system processor of the task of plotting individual display points.

Claims: A curve generator for a display system, the curve generator comprising arc generation means for generating an arc of a circle from data defining the locations of two end points and an intermediate point on the arc, wherein the arc generation means comprises initialisation means (40) for calculating the angle subtended between a first vector, from a first of the end points to the intermediate point, and a second vector, from the second of the end points to the intermediate point, and arc plotting means (44) for defining a succession of further vectors from said first end point and for calculating, for each further vector, its point of intersection with a counterpart vector, from said second end point, with which it subtends said angle, whereby a succession of further points of intersection are plotted on the circular arc. A curve generator as claimed in claim 1 wherein said arc plotting means (40) defines each said further vector in terms of a vector angle δ = nδ, where n is stepped from 1 to N for each of a succession of N further vectors each separated from the next by a stepping angle δ and wherein said arc plotting means calculates said point of intersection for each said further vector from the gradient of that further vector with respect to a reference direction and from a distance from said first end point along that further vector. A curve generator as claimed in claim 2 wherein said initialisation means (40) also comprises means for computing the length of a third vector, from said first to said second end points, and wherein said plotting means calculates the gradient of a said further vector as: π - δ n - b g₃₁ , and calculates said distance from said first end point along that further vector as: where b is said angle subtended between said first and second vectors, g₃₁ is the gradient of the third vector, and B is said distance between said first and second end points. A curve generator as claimed in any one of the preceding claims, wherein the curve generator is for generating a curve from data defining the locations of two end points and an intermediate point on the curve in reference space and a transform operator for mapping the curve onto the arc of a circle in mapped space, wherein the curve generator additionally comprises transform logic (74) for transforming data representative of the curve between reference and mapped space using the transform operator or its inverse as appropriate, and wherein the arc generation means is for generating an arc of a circle in mapped space from the locations of the two end points and the intermediate point on the curve when mapped into mapped space. A method of generating an arc of a circle in a display system comprising processing (80) and storage (90) means from data defining the locations of two end points and an intermediate point on the arc comprising the steps of: (a) calculating the angle subtended between a first vector, from a first of the end points to the intermediate point, and a second vector, from the second of the end points to the intermediate point; (b) storing the calculated angle in said storage means; (c) defining a further vector from said first end point; (d) calculating, for said further vector, its point of intersection with a counterpart vector, from the second end point, with which it subtends said calculated angle; and (e) repeating steps (c) and (d) for yet further vectors from said first end point whereby a succession of further points of intersection are plotted on the circular arc. A method of generating a curve in a display system comprising processing (80) and storage (90) means from data defining the locations of two end points and an intermediate point on the curve in reference space and a transform operator for mapping the curve onto the arc of a circle in mapped space, said method comprising the steps of: (a) mapping the three points from reference space into mapped space, (b) calculating the angle subtended between a first vector, from a first of the mapped end points to the mapped intermediate point, and a second vector, from the second of the mapped end points to the mapped intermediate point; (c) storing the calculated angle in said storage means; (d) defining a further vector from said first mapped end point; (e) calculating, for said further vector, its point of intersection with a counterpart vector, from the second mapped end point, with which it subtends said calculated angle; (f) transforming the point of intersection from mapped space into reference space using the transform operator; and (g) repeating steps (d) (e) and (f) for yet further vectors whereby a succession of further points of intersection are plotted on the curve in reference space.

Inventors: CLARKE, DAVID ANDREW, FARR, ROBERT WILLIAM ERIC
Applicants: INTERNATIONAL BUSINESS MACHINES CORPORATION
Application year: 1988
Main classification: G09G1
IPC codes: G09G1, G09G5, G06T11