TECHNICAL FIELD

The present invention relates to data compression techniques, and more particularly, to methods for compactly storing and recreating script and graphic data with the aid of a digital computer.

BACKGROUND ART

In an effort to increase the ease of use and efficiency of operation of digital computers, much work has been done in perfecting the man-machine interface. Graphic input tablets have been developed which allow the operator to input data directly into the digital computer by writing on the tablet with a pen or stylus while the stylus senses the x and y coordinates of the position of the stylus and inputs this data to the computer which can then store the data in the computer's electronic memory or reproduce the graphic input on a display device. One drawback of this kind of input device is that a large number of x and y coordinates must be stored by the computer in order to reproduce the graphic input with fidelity. Prior art systems have attempted to reduce this storage requirement by storing a starting point, a slope for a line and a length of the line. For example, U.S. Pat. No. 4,181,596 issued 1 Jan. 1980 to C. E. Schwab, et al discloses such a technique. The problem with this storage technique is that the fidelity of the image reproduced by the digital computer is limited since it can only reproduce straight line segments. While the fidelity of this reproduction can be improved by using many short line segments, this defeats the concept of compressing the data by requiring many more storage points. Nowhere in the prior art is a technique for highly compressing non-linear graphic input data and reproducing high fidelity images of the data disclosed.

DISCLOSURE OF THE INVENTION

A method is provided for compactly storing and faithfully reproducing script and graphic time series data from a computer tablet or optical scanner input device. The input device sends a time sequenced stream of x and y coordinates indicating the location of the data as it moves across the tablet. The beginning character point is stored. Subsequent points are buffered until a reversal in direction is detected with regard to either coordinate axis. This point, called the retro-coordinate sample point, terminates a buffered segment and starts a new buffered segment. The ending character point terminates the last segment of the stroke. Each segment is approximated by a circular arc which is retained by storing the horizontal displacement to the mid-point of the arc along with the x and y displacements to the ending retro-coordinate sample point. The retained information is sufficient to reproduce the input script or graphic for display, print or manipulation. This is accomplished by mathematically regenerating the arc segments in the sequence they were created.

BRIEF DESCRIPTION OF THE DRAWING

FIG. 1 is a block diagram of an interactive information processing system embodying the present invention;

FIG. 2 is a functional diagram of the microprocessor shown in FIG. 1;

FIG. 3 is a functional diagram illustrating the data flow path between portions of the memory and the microprocessor and the display refresh buffer; FIG. 4 is a diagrammatic view of the display shown in FIG. 1;

FIG. 5 illustrates the retro-coordinate sample points for a script character;

FIG. 6 illustrates a graphic figure with gradually changing curves;

FIG. 7 is a graphic representation of a window surrounding a number of input points;

FIGS. 8 and 9 illustrate the bounds through which a curve will pass in accordance with the window of FIG. 7;

FIGS. 10 and 11 illustrate the trigometrical determination of the mid-point of a circular arc;

FIG. 12 illustrates the process for determining the center of a circular arc; and

FIG. 13 illustrates the determination of a circular arc utilizing polar coordinates.

BEST MODE FOR CARRYING OUT THE INVENTION

The invention will now be described as embodied in an interactive information processing system of the type shown in FIG. 1. As shown in FIG. 1, the information processing system illustrated therein comprises a graphic input tablet 10, a microprocessor 11, a display refresh buffer 12, a display device 14, a printer 15, and an auxiliary direct access storage device 16 such as a disk drive. A clock 17, for keeping the various components of the system in synchronism, is also shown in FIG. 1 and is effectively coupled to each of the units.

Tablet 10 comprises a flat writing surface on which paper, film, etc., may be placed and then traced to input data to the microprocessor 11. Tracing or writing on the tablet 10 is done with pen 9 which is electrically coupled to the tablet 10. Graphic input tablets of the type referred to herein are commercially available, for example, the demi-pad 5™ manufactured by GTCO Corporation or the TEKTRONIX® 4953/4954 Graphics Tablet manufactured by Tektronix, Inc. could be used.

The tablet 10 is connected to the microprocessor 11 by means of a bus 20. The microprocessor 11, as shown in FIG. 2, comprises an input port 21, an output port 22, a random access memory 23, and a process execution unit 24.

Functionally, memory unit 23 stores both instructions and data in specified sections which will be described in more detail later on in the specification. Data is entered into memory 23 from the tablet 10 as pairs of x-y coordinates through input port 21. As shown in FIG. 3, the section of RAM 23 which received the data from the tablet 10 is designated input queue 26. Data to be displayed is transferred by a series of instructions from queue 26 to the buffer section 27 and then to the display refresh buffer 12 through output port 22 of the microprocessor 11 executing a series of move instructions.

The microprocessor 11 may be an IBM Series 1, INTEL model 8086, or any of the recognized functionally equivalent, currently available microprocessors.

The display refresh buffer 12 is shown as a separate buffer connected between the output port 22 and the display device 14. Buffer 12, in practice, is normally a part of the display device 14 and functions to control the generation of characters on the screen of the display device 14 by exercising on-off control of the electron beam to cause selected phosphorus dots to illuminate on the display screen.

The output port 22 also supplies data stored in memory 23 to the printer 15 and disk storage unit 16, each of which may have their own internal buffers which are not shown. Commands to transfer data from the random access memory 23 to the printer 15 or storage unit 16 are sent to the microprocessor 11 by the operator from the tablet 10 or a conventional keyboard (not shown) attached to the microprocessor 11.

Printer 15 may be any suitable printer known in the art. For example, the printer may be a standard dot matrix output printer such as a wire matrix or ink jet or x-y plotter.

Dish storage 16 may also be any suitable disk storage device which is capable of storing serial by byte data supplied to it at determined sector address locations, each of which are randomly addressable by the microprocessor 11 to retrieve the data. Spatially related data supplied from disk drive 16 is stored in the display data area 28 of the memory 23 in encoded form. The other section of memory 23 shown in FIG. 3 is the display format buffer area 29 which is involved in the handling of spatially related data in decoded form.

FIG. 4 is a schematic representation of the screen of display device 14. As shown in FIG. 4, the screen has, for example, the capability of displaying a matrix of dots in rows and columns to represent information input to the display. In practice, characters consist of a string of dot positions or picture elements sometimes referred to as pels. A character is represented on device 14 as a string of six wide by ten high pels, for example, which has been designated by reference character 32 in FIG. 4. The interaction of the refresh buffer 12 and the display 14 is to convert the data stored at a location in the buffer 12 to the corresponding character on the screen; this is accomplished by addressing the appropriate dots or pels on the screen in what is commonly called an all points addressable function. Display 14 generally is provided with its own set of electronics to achieve that conversion. The microprocessor 11 need only supply the address and load the buffer 12 with the appropriate characters.

The disk storage device 16 also is generally provided with its own set of electronics for converting a byte of data supplied from the display data area 28 of memory 23 through the output port 22 to a serial by bit stream of data to be recorded at a predetermined sector of the one addressed concentric recording track on the disk. Data from the device 16 is supplied to the microprocessor 11 serial by byte from the addressed sector and storage tracks when requested.

It will be understood that all of the above described functions and interactions involving the microprocessor 11 are achieved through suitable programs which are also stored in memory 23 and which are called into operation in response to data from the tablet 10 or interrupt signals generated by the various components of the system shown in FIG. 1.

The graphic input tablet 10 produces a stream of x,y pairs which identify the location of the pen 9 as it moves on the tablet 10. The input stream of coordinates is compressed by the present invention by selecting retro-coordinate sample (RCS) transition points which segment each stroke of the pen 9 into arcs. Each such arc is then approximated by a circular arc. This approximation is captured by computing and saving the x axis displacement (called delta-d) to the arc mid-point.

              TABLE 1______________________________________RETRO COORDINATE SAMPLE DETECTION ALGORITHMCOMPUTE QUALIFYING RCS POINTSDEFINITIONS:XY =      IMMEDIATE X OR Y COORDINATE     POINT OF THE X,Y PAIRLASTXY =  LAST X OR LAST Y COORDINATE POINT     OF THE X,Y PAIRDELTXY =  LAST RCS = X OR Y COORDINATE     POINT OF THE X,Y PAIRLAST SIGN =     SIGN OF LAST X AND LAST Y COOR-     DINATE X,Y PAIRIF PEN-DOWN IS DETECTED THENSAVE LAST PEN --UP AS RCS POINTSAVE THE PEN DOWN POSITION AS AN RCS POINTCOMPUTE WINDOW LIMIT FOR EXCLUSION OF RCS POINTSGET THE HEIGHT OF THE ARCSTORE THE ARC LIFT (DX)SET LASTXY TO PEN DOWN POSITIONINITIALIZE DIRECTION = SIGN (X(2) - X(1))COMPUTE SIGN OF THE SLOPE FOR LAZY --SELSE (NOT PEN --DOWN SAMPLE)DO FROM 2ND SAMPLE THRU PEN --UPDO FOR BOTH X AND Y COORDINATESGET XY SAMPLECOMPUTE DIRECTIONCOMPUTE DELTXY (LAST RCS - XY)COMPUTE 2ND DERIVATIVE AND HOLD SIGN FOR XYSAMPLEIF (DIRECTION CHANGED FROM LAST DIRECTION)THENIF THIS RCS CANDIDATE SAMPLE IS IN WINDOWLIMIT THENDO NOT FILE THIS SAMPLE AS AN RCS POINTELSESAVE LASTXY AS AN RCS POINTCOMPUTE WINDOW LIMIT FOR EXCLUSION OFRCS POINTSLASTSIGN=SIGN SAVE DERIVATIVE SIGNGET THE HEIGHT OF THE ARCSTORE THE ARC LIFT (DX)ENDIFELSE (MOVEMENT IN SAME DIRECTION)IF SIGN = LASTSIGN OR DELTXY > 126 THENTHE DERIVATIVE SIGN HAS CHANGED WHILEDIRECTION REMAINED THE SAMESAVE LASTXY AN RCS POINTGET THE HEIGHT OF THE ARCSTORE THE ARC LIFT (DX)COMPUTE WINDOW LIMIT FOR EXCLUSION OFRCS POINTSLASTSIGN=SIGN SAVE DERIVATIVE SIGNENDIFENDIFENDDOSET LASTXY = LAST X OR LAST Y = YENDDOENDIFEND RCS______________________________________

Referring now to FIG. 5, assume that a lowercase script 1 50 is written on the tablet 10. As the pen 9 moves the sample rate is constant and each small circle exemplifies the points sampled and thereby available for consideration by the RCS algorithm shown in Table 1. When the pen touches the tablet 10, point 51 will be the first RCS point saved. As the stroke progresses up and to the right, point 52 is the second RCS point saved as the x direction changes from right to left. The y axis movement continues its direction from bottom to top until at point 53 the direction now changes to a downward movement. At point 54 the x direction changes from a left to right movement, signalling that this point is to be saved. The y movement ceases its downward movement at point 55 and the final RCS point is generated at pen up point 56.

There are four conditions which applicants have discovered that produce an RCS point: (1) pen down and pen up, (2) reversal in direction of motion with respect to the x or y coordinates, (3) the delta x or delta y from the last selected RCS point exceeds a predetermined limit, for example, 126 grid units, and (4) the second derivative of the curve changes sign while within the constraints of the other three conditions. Conditions 1 and 2 are the primary RCS criteria. The addition of conditions 3 and 4 provide special coverage conditions and may be thought of as optional processes.

The pen down and pen up points are automatically picked as RCS points so the exact end points of each stroke are retained. As the pen motion proceeds from pen down (x(1), y(1)), each new point is examined to determine the direction of motion. When a change in direction is detected with respect to either axis, the transition point (x(s), y(s)) is saved as an RCS point. The following equation defines the relationship which generates an RCS point at x(i-1), y(i-1): ##EQU1##

To handle long strokes applicants' algorithm was designed to store an RCS point when 126 units have been traversed in any direction since the last RCS point. The default value of 126 units is a designer's choice based on eight bits being used to save the displacement between RCS points which permits a value of ±126 units.

The coverage provided by the above three conditions will address most of the cases for high quality image generation using the approach defined by the present invention. However, there exists the possibility for a principal to scribe in a manner that will not capture an RCS point. To cover such situations, a higher order point of change is defined. Refer to FIG. 6a which shows the graph of a sine wave between the x values -π/2 and π/2 radians representing a segment of an input stroke. Since the first derivative of this curve as shown in FIG. 6b changes sign at the points 64 and 66, only these points would be captured by the RCS logic as described by condition 2 above. Since the midpoint 62 is on a straight line connecting 61 and 63, this curve would be represented by a straight line. Point 62 in FIG. 6a is clearly a point of change which needs to be captured. Mathematically speaking, such a point is called an "inflection point", and is defined as a point where the second derivative of a continuous curve changes signs, i.e. where the curve changes from a concave to convex or convex to concave shape. FIG. 6c illustrates the second derivative of the sine wave. Point 68 occurs at the location where the second derivative changes signs from positive to negative and so is an inflection point that can be recognized. Letting x(i), y(i), i=1,2,3 . . . represent the time sequenced coordinates of the sample points on a curve, this logic can be invoked as follows:

If, ##EQU2## then an RCS point is generated at x(i-1), y(i-1). The name "lazy-s" for these points refers to the general shape of FIG. 6a, which is but one of the orientations where the capture of inflection points is significant. However, it should be noted that the capture of these points is not crucial to the overall effectiveness of the present invention. The quality of a complex image, where many RCS points are involved but the relatively few "lazy-s" points are neglected, is usually acceptable.

Since the input tablet 10 may not be an ideal device, more than one point may be generated as the location of the pen 9 as it moves across the tablet 10. It is desirable to have the path smooth so as not to generate extraneous RCS points. To assure this, the present invention computes a window based on the input and output device height and width maximums. This window is the basis for discarding otherwise unnecessary input points and improving overall compaction. FIG. 7 is a graphic representation of the window. Point 72 is the first entry within the window, therefore no other entry is allowable in the same sequence. Points 73, 74 and 75 are discarded due to the window constraint. Point 76 is the next allowable input point to be evaluated as an RCS point.

In order to compactly store and faithfully reproduce written input from the graphic input table 10 it is necessary to emulate the naturally smooth motions of a writer while storing as little information as possible. The circle was chosen as the approximating curve in order to achieve the most fundamental design which results in a very fast execution time when implemented as a computer program. Given two RCS endpoints defined by the coordinates x(1),y(1) and x(n),y(n) and the intermediate sampling points defined by the coordinates (x(2),y(2)), (x(3),y(3)), . . . , (x(n-1),y(n-1)), the objective is to find a point x(d),y(d) that describes the circle being fit as an approximation to the curve between x(1),y(1) and x(n), y(n). Refer now to FIGS. 8 and 9. The shaded area 80 in FIG. 8 bounds the space through which a concave curve between the RCS points 81 and 82 can pass. FIG. 9 shows such a typical curve, where point 92 represents the coordinates of any of the intermediate sampling points. Referring now to FIG. 10, assume the two points x(1 ),y(1) and x(n),y(n) are on a circle 100, then the center 103 of the circle defined by the coordinates x(0),y(0) lies on the line 101 which bisects and is perpendicular to the line 102 between the two points 104 and 105. The line 102 is defined by equation III and line 101 is defined by equation IV. ##EQU3##

The circle containing x(1),y(1) and x(n),y(n) can be uniquely defined by specifying one additional point on the circle. A convenient choice is the point where the line 101 intersects the circle 100, since only one coordinate (either x or y) of this point needs to be stored. The other coordinate can be recovered through equation IV. In addition, this point is a natural choice since it represents the maximum distance from the chord between x(1),y(1) and x(n),y(n) to the concave arc of the circle between these points. Thus, the point gives a measure of the maximum "bend" of the curve between the endpoints.

The technique for finding the point 116 defined by coordinates x(d),y(d) is shown in FIG. 11. A horizontal vector 114 is drawn from point 115 to the line 101. The value of the vector 114 is the x value of the intersection of vector 114 and the line 101. The sign of the vector 114 is the critical parameter. A similar horizontal vector 113 is then drawn from each sampled point x(s),y(s), until the point is found that touches or crosses the line 101, i.e., until vector 113 equals zero or changes sign. At this point a perpendicular projection is dropped to line 101 and this becomes the point x(d),y(d). The displacement x(d)-x(1) is stored as delta-d, completing the definition of the circle. Delta-d is defined by equation V. ##EQU4##

The limiting case of the circle is a straight line. This geometry is indicated when there are no intermediate sampling points between x(1),y(1) and x(n),y(n). In addition, if dx equals zero or dy equals zero, a straight line is the only possible curve between the RCS points, as can be extrapolated from the bounds in FIG. 8. The algorithm for calculating delta-d covers this situation by setting delta-d equals dx/2. A major computational advantage to this algorithm is that only approximately half of the sampled points need be considered in the calculations to find delta-d. An example program routine for calculating delta-d is shown in Table 2 implemented in pseudo code.

              TABLE 2______________________________________BEGIN DELTA-D SELECTIONGET THE N (X,Y) PAIRS OF COORDINATES BETWEENTHE RCS POINTSIF THE RCS POINTS SPAN A STRAIGHT LINETHEN DELTA-D = (DELTA-X)/2ELSE DOT1 =      THE SIGN OF THE HORIZONTAL DISTANCE FROM THE     FIRST RCS POINT TO THE PERPENDICU-     LAR BISECTOR OF THE LINE BETWEEN     THE RCS POINTSDELTA-D = 0DO I =    2 TO N WHILE (DELTA-D = 0)T =       THE SIGN OF THE HORIZONTAL DISTANCE FROM THE     ITH INTERMEDIATE POINT TO THE     PERPENDICULAR BISECTOR OF THE     LINE BETWEEN THE RCS POINTSIF (T = 0) OR (T ≠ T1)THEN COMPUTE DELTA-D AS A FUNCTION OFT,X(I),Y(I)ENDIFENDDOENDDOENDIFOUTPUT DELTA-DEND DELTA-D SELECTION______________________________________

In a mode for implementing the invention, the input from tablet 10 consists of first quadrant values for each x and y pair. The pen down position is indicated by affixing a minus sign to the x coordinate of this x,y pair. Storage for the input pairs is 16 bits per coordinate, thereby using 32 bits per x,y pair. The technique for RCS point storage uses the full 32 bits for the pen down x,y pair. Following the pen down x,y is the delta-d (d) for the first segment. This occupies eight bits. Finally another dx and dy pair follows the delta-d, each occupying eight bits. The dx and dy pair are the displacement from the original pen down x,y. These sequence of three 8-bit values (d, dx,dy) is repeated for each segment until the pen up RCS is detected. Every dx,dy after the pen down is the displacement from the preceding x,y. It must be remembered that the dx,dy is the terminal point for one arc and an initial point for the next arc in the sequence. A negative sign is used on the first x value to indicate pen down on the RCS points.

In order to reconstruct the stored image, the arc generation routine shown in Table 3 is invoked. The image is stored as a set of triplets (two adjacent RCS points plus the intermediate point represented by delta-d, where each triplet defines an arc or segment of the image.

              TABLE 3______________________________________BEGIN ARC GENERATIONGET THE VALUES X(1),Y(1) (COORDINATES OF 1STRCS POINT)DELTA-D (DISPLACEMENT TO THE INTERMEDIATEPOINT)DELTA-X,DELTA-Y (DISPLACEMENT TO THESECOND RCS POINT)AND N (THE NUMBER OF POINTS TO GENERATE ONTHE ARC)IF DELTA-D = (DELTA-X)/2THEN GENERATE POINTS ON A STRAIGHT LINEBETWEEN THE RCS POINTSELSE DOCOMPUTE THE SECOND RCS POINT X(N),Y(N)DETERMINE THE INTERMEDIATE POINT X(D),Y(D)FROM DELTA-DDETERMINE THE RADIUS OF THE CIRCLE TO BEGENERATED (R)DETERMINE THE CENTER POINT OF THE CIRCLE(X(0),Y(0))THETA-1 =    THE ANGLE BETWEEN A HORI-        ZONTAL LINE AND THE LINE        FROM THE CENTER OF THE CIR-        CLE TO X(1),Y(1)THETA-D =    THE ANGLE BETWEEN A HORI-        ZONTAL LINE AND THE LINE        FROM THE CENTER OF THE CIR-        CLE TO X(D),Y(D)THETA-N =    THE ANGLE BETWEEN A HORI-        ZONTAL LINE AND THE LINE        FROM THE CENTER OF THE CIR-CLE TO X(N),Y(N)DELTA-THETA =        ((THETA-N) - (THETA-1))/N-1DO I = 2 TO N-1GENERATE A POINT AT THE COORDINATES:X(I) =       X(0) + R*COS(THETA-1 +        (I-1) *DELTA-THETA)Y(I) =       Y(0) + R*SIN(THETA-1 +        (I-1) *DELTA-THETA)ENDDOENDDOENDIFOUTPUT THE N GENERATED POINTS ON THE ARC OF A CIRCLEEND ARC GENERATION______________________________________

The first RCS endpoint x(1),y(1) is equivalent to the second RCS endpoint of the previous arc. The second RCS endpoint x(n),y(n) is recovered by adding the stored dx and dy values to x(1) and y(1), respectively, i.e.,

x(n)=x(1)+dx                                               VI

y(n)=y(1)+dy                                               VII

The intermediate point on the arc x(d),y(d) is found from,

x(d)=x(1)+delta-d                                          VIII ##EQU5## A circle is uniquely defined by these three points on its arc. From this information, the center and radius of the circle can be determined, thereby enabling the generation of additional points on the arc between x(1),y(1) and x(n),y(n).

The simplest procedure for finding the center 126 of the circle 120 is by drawing perpendicular bisectors of the lines joining any two pairs of points on the circle, and then solving for the intersection of these bisecting lines. FIG. 12 illustrates the process. The equation for line 121 at x(0),y(0) is ##EQU6## and the equation for line 122 at x(0),y(0) is ##EQU7## Solving for x(0), ##EQU8## y(0) is then found from either equation X or equation XI.

The radius R is the distance from the center of the circle to any point on the circle,

R=SQRT((x(1)-x(0)**2+(y(1)-y(0))**2)                       XIII

Assume now that n points are to be generated on the arc (including the endpoints x(1),y(1) and x(n),y(n)) and that these points are to define uniform sub-arcs. The use of polar coordinates is therefore indicated to facilitate this definition. FIG. 13 represents the concept. From FIG. 13, ##EQU9## The angle θ(i) of intermediate point x(i),y(i) is then, ##EQU10## And so the cartesian coordinates of the intermediate point are,

x(i)=x(0)+R cos θ(i)                                 XVII

y(i)=y(0)+R sin θ(i)                                 XVIII

By looping through equation XVI through XVIII, all intermediate points on the arc are generated.

While the invention has been particularly shown and described with respect to the preferred embodiment, it will be understood by those skilled in the art that various changes in form and detail may be made without departing from the spirit and scope of the invention.