WO1999005851A2 - Method of isomorphic singular manifold projection still/video imagery compression - Google Patents

Abstract

The method of still image compression uses isomorphic singular manifold projection (1) whereby surfaces of objects having singular manifold representations are represented by best match canonical polynomials to arrive at a model representation (2). The model representation (3) is compared with the original representation (4) to arrive at a difference (5). If the difference exceeds a predetermined threshold, the difference data are saved and compressed using standard lossy compression (6). The coefficients from the best match polynomial together with the difference data, if any, are then compressed using lossless compression (6).

Description

METHOD OF ISOMORPHIC SINGULAR MANIFOLD PROJECTION STILL/VIDEO IMAGERY COMPRESSION

Field of the Invention

The present invention relates to image compression systems, and in particular relates to an image compression system which provides hypercompression.

Background of the Invention

Image compression reduces the amount of data necessary to represent a digital image by eliminating spatial and/or temporal redundancies in the image information. Compression is necessary in order to efficiently store and transmit still and video image information. Without compression, most applications in which image information is stored and/or transmitted would be rendered impractical or impossible.

Generally speaking, there are two types of compression: lossless and lossy. Lossless compression reduces the amount of image data stored and transmitted without any information loss, i.e. , without any loss in the quality of the image. Lossy compression reduces the amount of image data stored and transmitted with at least some information loss, i.e., with at least some loss of quality of the image.

Lossy compression is performed with a view to meeting a given available storage and/or transmission capacity. In other words, external constraints for a given system may define a limited storage space available for storing the image information, or a limited bandwidth (data rate) available for transmitting the image information. Lossy compression sacrifices image quality in order to fit the image information within the constraints of the given available storage or transmission capacity. It follows that, in any given system, lossy compression would be unnecessary if sufficiently high compression ratios could be achieved, because a sufficiently high compression ratio would enable the image information to fit within the constraints of the given available storage or transmission capacity without information loss.

The vast majority of compression standards in existence today relate to lossy compression. These techniques typically use cosine-type transforms like DCT and wavelet compression, which are specific types of transforms, and have a tendency to lose high frequency information due to limited bandwidth. The "edges" of images typically contain very high frequency components because they have drastic gray level changes, i.e., their dynamic range is very large. Edges also have high resolution. Loss of edge information is undesirable because resolution is lost as well as high frequency information. Furthermore, human cognition of an image is primarily dependent upon edges or contours. If this information is eliminated in the compression process, human ability to recognize the image decreases.

Fractal compression, though better than most, suffers from high transmission bandwidth requirements and slow coding algorithms. Another type of motion (video) image compression technique is the ITU-recommended H.261 standard for videophone/ videoconferencing applications. It operates at integer multiples of 64 kbps and its segmentation and model based methodology splits an image into several regions of specific shapes, and then the contour and texture parameters representing the region boundaries and approximating the region pixels, respectively, are encoded. A basic difficulty with the segmentation and model-based approach is low image quality connected with the estimation of parameters in 3-D space in order to impart naturalness to the 3-D image. The shortcomings of this technique are obvious to those who have used videophone/videoconferencing applications with respect particularly to MPEG video compression.

Standard MPEG video compression is accomplished by sending an "I frame" representing motion every fifteen frames regardless of video content. The introduction of I frames asynchronously into the video bitstream in the encoder is wasteful and introduces artifacts because there is no correlation between the I frames and the B and P frames of the video. This procedure results in wasted bandwidth. Particularly, if an I frame has been inserted into B and P frames containing no motion, bandwidth is wasted because the I frame was essentially unnecessary yet, unfortunately, uses up significant bandwidth because of its full content. On the other hand, if no I frame is inserted where there is a lot of motion in the video bitstream, such overwhelming and significant errors and artifacts are created that bandwidth is exceeded. Since the bandwidth is exceeded by the creation of these errors, they will drop off and thereby create the much unwanted blocking effect in the video image. In the desired case, if an I frame is inserted where there is motion (which is where an I frame is desired and necessary) the B and P frames will already be correlated to the new motion sequence and the video image will be satisfactory. This, however, happens only a portion of the time in standard compression techniques like MPEG. Accordingly, it would be extremely beneficial to insert I frames only where warranted by video content.

The compression rates required in many applications including tactical communications are extremely high as shown in the following example making maximal compression of critical importance. Assuming 512² number of pixels, 8- bit gray level, and 30 Hz full-motion video rate, a bandwidth of 60 Mbps is required. To compress data into the required data rate of 128 kbps from such a full video uncompressed bandwidth of 60 Mbps, a 468: 1 still image compression rate is required. The situation is even more extreme for VGA full-motion video which requires 221 Mbps and thus a 1726: 1 motion video compression rate. Such compression rates, of course, greatly exceed any compression rate achievable by state of the art technology for reasonable PSNR (peak signal to noise ratio) values of approximately 30 dB. For example, the fourth public release of JPEG has only a 30: 1 compression rate and the image has many artifacts due to a PSNR of less than 20 dB, while H320 has a 300:1 compression ratio for motion and still contains many still/motion image artifacts.

The situation is even more stringent for continuity of communication when degradation of power budget or multi-path errors of wireless media further reduce the allowable data rate to far below 128 kbps. Consequently, state of the art technology is far from providing multi-media parallel channelization and continuity data rates at equal to or lower than 128 kbps. Very high compression rates, high image quality, and low transmission bandwidth are critical to modern communications, including satellite communications, which require full-motion, high resolution, and the ability to preserve high-quality fidelity of digital image transmission within a small bandwidth communication channel (e.g. TI). Unfortunately, due to the above limitations, state of the art compression techniques are not able to transmit high quality video in real-time on a band-limited communication channel. As a result, it is evident that a compression technique for both still and moving pictures that has a very high compression rate, high image quality, and low transmission bandwidth and a very fast decompression algorithm would be of great benefit. Particularly, a compression technique having the above characteristics and which preserves high frequency components as well as edge resolution would be particularly useful.

In addition to transmission or storage of compressed still or moving images, another area where the state of the art is unsatisfactory is in automatic target recognition (ATR). There are numerous applications, both civilian and military, which require the fast recognition of objects or humans amid significant background noise. Two types of ATR are used for this purpose, soft ATR and hard ATR. Soft ATR is used to recognize general categories of objects such as tanks or planes or humans whereas hard ATR is used to recognize specific types or models of objects within a particular category. Existing methods of both soft and hard ATR are Fourier transform-based. These methods are lacking in that Fourier analysis eliminates desired "soft edge" or contour information which is critical to human cognition. Improved methods are therefore needed to achieve more accurate recognition of general categories of objects by preserving critical "soft edge" information yet reducing the amount of data used to represent such objects and thereby greatly decrease processing time, increase compression rates, and preserve image quality.

Summary of the Invention The present invention is based on Isomorphic Singular Manifold Projection

(ISMP) or Catastrophe Manifold Projection (CMP). This method is based on Newtonian polynomial space and characterizes the images to be compressed with singular manifold representations called catastrophes. The singular manifold representations can be represented by polynomials which can be transformed into a few discrete numbers called "datery" (number data that represent the image) that significantly reduce information content. This leads to extremely high compression rates (CR) for both still and moving images while preserving critical information about the objects in the image.

In this method, isomorphic mapping is utilized to map between the physical boundary of a 3-D surface and its 2-D plane. A projection can be represented as a normal photometric projection by adding the physical parameters, B (luminance) to generic geometric parameters (X,Y). This projection has a unique 3-D interpretation in the form of a "canonical singular manifold" . This manifold can be described by a simple polynomial and therefore compressed into a few discrete numbers resulting in hyper compression. In essence, any image is a highly correlated sequence of data. The present invention "kills" this correlation, and image information in the form of a digital continuum of pixels almost disappears. All differences in 2-D "texture" connected with the 2-D projection of a 3-D object are "absorbed" by a contour topology, thus preserving and emphasizing the "sculpture" of the objects in the image. This allows expansion with good fidelity of a 2-D projection of a real 3-D object into an abstract (mathematical) 3-D object and is advantageous for both still and video compression and automatic target recognition.

More particularly, using catastrophe theory, surfaces of objects may be represented in the form of simple polynomials that have single-valued (isomorphic) inverse reconstructions. According to the invention, these polynomials are chosen to represent the surfaces and are then reduced to compact tabulated normal form polynomials which comprise simple numbers, i.e. , the datery, which can be represented with very few bits. This enables exceptionally high compression rates because the "sculpture" characteristics of the object are isomorphically represented in the form of simple polynomials having single-valued inverse reconstructions. Preservation of the "sculpture" and the soft edges or contours of the object is critical to human cognition of the image for both still and video image viewing and ATR. Thus, the compression technique of the present invention provides exact representation of 3-D projection edges and exact representation of all the peculiarities of moving (rotating, etc.) 3-D objects, based on a simple transition between still picture representation to moving pictures. In a preferred embodiment the following steps may be followed to compress a still image using isomorphic singular manifold projections and highly compressed datery. The first step is to subdivide the original image, I₀, into blocks of pixels, for example 16 x 16 or other sizes. These subdivisions of the image may be fixed in size or variable. The second step is to create a "canonical image" of each block by finding a match between one of fourteen canonical polynomials and the intensity distribution for each block or segment of pixels. The correct polynomial is chosen for each block by using standard merit functions. The third step is to create a model image, I_M, or "sculpture" of the entire image by finding connections between neighboring local blocks or segments of the second step to smooth out intensity (and physical structure to some degree). The fourth step is to recapture and work on the delocalized high frequency content of the image, i.e., the "texture". This is done by a subtraction of the model image, I_M, generated during the third step from the original segmented image, IQ, created during the first step. A preferred embodiment of this entire still image compression process will be discussed in detail below.

Optimal compression of video and other media containing motion may be achieved in accordance with the present invention by inserting I frames based on video content as opposed to at fixed intervals (typically every 15 frames) as in the in prior art motion estimation methods. In accordance with the motion estimation techniques of the present invention, the errors between standard "microblocks" or segments of the current frame and a predicted frame are not only sent to the decoder to reconstruct the current frame, but, in addition, are accumulated and used to determine the optimal insertion points for I frames based on video content. Where the accumulated error of all the microblocks for the current frame exceeds a predetermined threshold which itself is chosen based upon the type of video (action, documentary, nature, etc.), this indicates that the next subsequent frame after the particular frame having high accumulated error should be an I frame. Consequently, in accordance with the present invention, where the accumulated errors between the microblocks or segments of the current frame and the predicted frame exceed the threshold, the next subsequent frame is sent as an I frame which starts a new motion estimation sequence. Consequently, I frame insertion is content dependent which greatly improves the quality of the compressed video. The I frames inserted in the above compression technique may first be compressed using standard DCT based compression algorithms or the isomorphic singular manifold projection (ISMP) still image compression technique of the present invention for maximal compression. In either case, the compression techniques used are preferably MPEG compatible.

Additionally, using the motion estimation technique compression of the present invention, compression ratios can be dynamically updated from frame to frame utilizing the accumulated error information. The compression ratio may be changed based on feedback from the receiver and, for instance, where the accumulated errors in motion estimation are high, the compression ratio may be decreased, thereby increasing bandwidth of the signal to be stored. If, on the other hand, the error is low, the compression ratio can be increased, thereby decreasing bandwidth of the signal to be stored.

Because the present invention is a 3-D non-linear technique that produces high level descriptive image representation using polynomial terms that can be represented by a few discrete numbers or datery, it provides much higher image compression than MPEG (greater than 1000: 1 versus 100: 1 in MPEG), higher frame rate (up to 60 frames/ sec versus 30 frames/sec in MPEG), and higher picture quality or peak signal to noise ratio (PSNR greater than 32 dB versus PSNR greater than 23 dB in MPEG). Consequently, the compression technique of the present invention can provide more video channels than MPEG for any given channel bandwidth, video frame rate, and picture quality.

Description of the Drawings

Figure IA illustrates a singularity called a "fold"; and Figure IB illustrates a singularity called a "cusp"; Figure 2A illustrates "Newton Diagram Space" and contains "monoms" and polynomials;

Figure 2B illustrates the application of Newton diagram space in the context of ISMP Theory. Canonical and normal forms, etc.; Figure 3A depicts a fold and Figure 3B depicts tuck;

Figure 4 A illustrates a reflection from a manifold; Figure 4B illustrates the reflection from a manifold depending upon angle 0;

Figure 4C illustrates the projection on a display; Figure 5 illustrates a cylinder with constant luminance dependence;

Figure 6 illustrates an evolvement for F=X³;

Figures 7 A and 7B illustrate a reflector of a group for three mirrors R²; Figure 8 illustrates a smooth curve projection, representing movement and a physical object, and a catastrophic frame change, and positionally zoom camera changes;

Figure 9 is a abbreviated flow chart of the inventive ISMP still image compression method;

Figure 10 is a detailed flow chart of the inventive ISMP compression method in accordance with the present invention; Figure 11A illustrates an original image with an enlarged edge contour;

Figure 1 IB shows a 2-D CCD image of the enlarged edge contour; Figure 11C illustrates a model surface of the original edge contour in accordance with the present invention;

Figure 12 is an illustration of a segment of an original frame in accordance with the present invention;

Figure 13 is an illustration of connecting segments of a frame in accordance with the present invention;

Figure 14 illustrates subtracting I_M from IQ in accordance with the present invention; Figure 15 is a flow chart of the decoding process for ISMP compression in accordance with the present invention; Figure 16 is a flow chart of the motion estimation process in accordance with the present invention;

Figure 17 is a circuit schematic of the hardware uses for motion estimation in accordance with the present invention; Figure 18 is a flow chart of the error accumulation method of the present invention; and

Figure 19 is a table showing the results of data communication with varying data rates in accordance with the present invention.

Description of the Preferred Embodiments Preliminary Discussion of the Mapping of Surfaces Using Catastrophe Theory

The following is a brief introduction to Catastrophe Theory which may be helpful in understanding novel compression methodology of the present invention. Further discussion may be found in B.I. Arnold, Catastrophe Theory, Springer- Verlag 1992, which is incorporated by reference herein. Catastrophes are abrupt changes arising as a sudden response of a system to a smooth change in external conditions. In order to understand catastrophe theory, it is necessary to understand Whitney's singularity theory. A mapping of a surface onto a plane associates to each point of the surface a point of the plane. If a point on the surface is given coordinates (x_1; x₂) on the surface, and a point on the plane is given coordinates (y_l5 y₂), then the mapping is given by a pair of functions y_! = x₂) and y₂ = f₂(Xι, x₂). The mapping is said to be smooth if these functions are smooth (i.e., are differentiable a sufficient number of times, such as polynomials for example). Mappings of smooth surfaces onto a plane exist everywhere. Indeed, the majority of objects surrounding us are bounded by smooth surfaces. The visible contours of bodies are the projections of their bounding surfaces onto the retina of the eye. By examining the objects surrounding us, for instance, people's faces, the singularities of visible contours can be studied. Whitney observed that generically (for all cases bar some exceptional ones) only two kinds of singularities are encountered. All other singularities disintegrate under small movements of the body or of the direction of projection, while these two types are stable and persist after small deformations of the mapping.

An example of the first kind of singularity, which Whitney called a fold, is the singularity arising at equatorial points when a sphere is projected onto a plane such as shown in Figure IA. In suitable coordinates, this mapping is given by the formulas

i = ^χι > y₂ ^{= χ}ι

The projections of surfaces of smooth bodies onto the retina have just such a singularity at generic points, and there is nothing surprising about this. What is surprising is that besides the singularity, the fold, we encounter everywhere just one other singularity, but it is practically never noticed.

The second singularity was named the cusp by Whitney, and it arises when a surface like that in Figure IB is projected onto a plane. This surface is given by the equation

with respect to spatial coordinates (x,, x₂, y_f) and projects onto the horizontal plane (x₂, y.).

On the horizontal projection plane, one sees a semicubic parabola with a cusp (spike) at the origin. This curve divides the horizontal plane into two parts: a smaller and a larger one. The points of the smaller part have three inverse images (three points of the surface project onto them), points of the larger part only one, and points on the curve, two. On approaching the curve from the smaller part, two of the inverse images (out of three) merge together and disappear (here the singularity is a fold), and on approaching the cusp all three inverse images coalesce.

Whitney proved that the cusp is stable, i.e. , every nearby mapping has a similar singularity at an appropriate nearby point (that is, a singularity such that the deformed mapping, in suitable coordinates in a neighborhood of the point mentioned, is described by the same formulas as those describing the original mapping in a neighborhood of the original points). Whitney also proved that every singularity of a smooth mapping of a surface onto a plane, after an appropriate small perturbation, splits into folds and cusps. Thus, the visible contours of generic smooth bodies have cusps at points where the projections have cusp singularities, and they have no other singularities. These cusps can be found in the lines of every face or object. Since smooth mappings are found everywhere, their singularities must be everywhere also, and since Whitney's theory gives significant information on singularities of generic mappings, this information can be used to study large numbers of diverse phenomena and processes in all areas of science. This simple idea is the whole essence of catastrophe theory.

Technical Foundation of Catastrophic Theory

Catastrophic Manifold Projection (CMP) or Isomorphic Singular

Manifold Projection (ISMP)

The following glossary is a useful aid to understanding catastrophe theory because many of the terms used to describe it are uncommon in mathematics.

• 2-D Cartesian (Plane) Coordinates refer to standard (u, v) coordinates that describe a plane projection.

• 2-D Generalized Coordinates : (£ , v) describe a system through a minimum number of geometrical coordinates (i.e. , a number of degrees of freedom). These are usually curvilinear local coordinates, which belong to a specific surface in the vicinity of some point (i.e. , origin of coordinates).

• 3-D Cartesian Coordinates refer to (x,y,z) describing a common surface in 3-D: F (x,y,z) = 0. • 3-D Cartesian (Hyperplane) Coordinates, are: (u, v, w), where (u, v) are 2-D (plane) Cartesian coordinates; w is a third, new physical coordinate, related to luminance (B) and describing a "gray level" scale color scale.

• Arnold (Vladimir) is a Russian mathematician, who is a major contributor to catastrophe theory.

• Arnold Theorem (Local Isomorphism) A family of transformations can transform any given mapping into a set of canonical transformations by using smooth substitutions of coordinates. The Arnold theorem defines local isomorphism in a sense that defines a class of locally isomorphic functions.

1. Arnold proved that Thorn's theory can be represented in terms of group theory. 2. He also introduced an elegant theory for construction of the canonical form of singularities as they apply to wave front propagation in Lagrangian mechanics. 3. Furthermore, Arnold introduced methods based on using algebra of vector fields v . ∑

- ^R>^~ax ^• where Rj is a polynomial. and introduced a method of spectral series for reduction of arbitrary functions to normal form. 4. Finally, he introduced classification of singularities and a method that described how to determine any type of singularity within a list of singularities.

• Canonical Form is a generic mathematical term that can be defined in various ways. In the specific context of the Arnold theorem, the canonical form is the simplest polynomial, with the highest degree of monoms within the normal form area, representing a given type of catastrophe. The canonical form is represented by a segmented line in a Newton diagram.

• Canonical Transformation permits transformation of real surface form (such as F (u, v, w) = 0) into canonical form (i.e. , superposition of Morse form and singular residuum, or Thorn form) into two blocks of so-called canonical coordinates: regular and catastrophic.

• Catastrophe (a term invented by Montel) is a specific manifold mapping feature by which some points lying in the projection plane can abruptly change location in manifold. More philosophically, it "describes the emergence of discrete structures from the typical surface described in the platform of continuum."

• Catastrophes, Critical Number in 3-D (for mapping a generic surface onto a plane) is only two (2): fold and cusp (tuck). Using these two catastrophes is sufficient for static still imagery.

• Catastrophes, Total Number in 3-D (for mapping a generic surface onto a plane) is fourteen (14). Only "fold "-catastrophe does not have degenerate points; all the (13) others have. Using all 14 catastrophes is necessary in hypercompression if we consider dynamic imagery (or video).

• Catastrophic Manifold Projection is a fundamental concept of "3-D into 3-D" mapping, leading to hypercompression. This is diffeomorphic mapping, including geometrical coordinates (2-D generalized, and 2-D plane), as well as a fourth "photometric coordinate".

• Catastrophic Manifold Projection (CMP) Law is mapping:

(ξ, v, B) - (u, v, w)

Thus, the CMP is "3-D into 3-D" mapping, with two types of coordinates: "geometrical" (£, * ;(u, v), and "photometric" (B, w).

• The Critical Point is a point at which the rank of Jacobian is less than maximal (examples are maxima, minima, and bending points)

• Datery result from a novel mathematical procedure leading to a tremendous compression ratio; instead of describing some surface by a continuum, we describe this singular manifold by a few low even numbers, i.e. , datery. Therefore, during hypercompression, the surface as continuum "disappears" , leaving typical data (such as computer data).

• In a Degenerate Critical Point, the rank of Jacobian is a less than maximal rank minus one. This point can be a critical point of cusp catastrophes, for example.

• Discrete Structures are singular manifolds that can be described by a set of discrete, usually even, data (e.g., (2,5,-1,3)), leading to datery instead of description by a continuum of points (such as F(x,y,z) = 0). Such discrete structures (which are, in fact, continuums, but are still described by discrete sets), are typically referred to as singularities, bifurcation, or catastrophes. • Diffeomorphism is a stronger term than isomorphism (or homeomoφhism), and means not only isomorphism, but also smooth mapping.

• Discrete Structures are singular manifolds that can be described by a set of discrete data (e.g. , (2,5,-1 ,3)), leading to datery, instead of description by continuum of points (such as F(x,y,z) = 0). These discrete structures (which are, in fact, continuums, but are still described by discrete sets), are typically referred to as singularities, bifurcation, or catastrophes. • Field, a subset of a ring. (All non-zero field elements generate a group, by multiplication.) For example, the differential operator can be an element of a field.

• Generalized Coordinates, or, more precisely, generalized coordinates of Lagrange, are such "natural" coordinates in solid state mechanics that their number is precisely equal to the number of a body's degrees of freedom: {ζ,v,η,...)

• A Generic Surface, in the context of the CMP method, is a mathematical surface which, within infinitesimal changes, does not have the same tangent (or projection line) for more than two points along any curve lying on this surface. In other words, a surface is generic if small changes to the surface do not lead to changes in singularities (such as splitting into less complex ones.) Physical surfaces are almost always generic because of noise tolerance.

• A Group is the simplest set in mathematical models, with only a single operation.

• Hypercompression is a specific compression term which provides a datery (i.e., "stripping" a continuum surface into its discrete representation). This is possible for the surface locality in the form of catastrophe. • Ideal is another subset of a ring. A subset of the ring is an ideal if this subset is a subgroup of the ring by summation. In the context of the Arnold theorem, this summation group is a set of all monoms that lie above the canonical form segmented line.

• Isomorphic Singular Manifold Projection - see CMP.

• Jacobian is a transformation matrix whose element, Hy, can be presented in the form:

du_t

"" ^"A, where Uj = (u,v,w) is the plane projection Cartesian coordinate system, ξ_j = (ξ, v, η) is the generalized coordinate system, and

du

describes a partial derivative.

• Landau (Lev) was a Russian theoretical physicist, who won the Nobel Prize for superfluidity of the isotope helium He₃. He systematically applied the catastrophe theory approach before this theory was mathematically formulated ⁽³⁾.

• The Landau Method, applied in the second-order phase transition, applies Thorn's lemma in the form of the Taylor series, including only "important" physical terms.

• Lie Algebra is an algebra belonging to the Lie groups with a binary operation (commutator).

• A Lie Group is a group whose generator is an infinitesimal operator. • Locally Isomorphic functions have the same singular residuum

(see Thorn's lemma); thus, they can be compressed identically for "soft edges", or "object boundaries" .

• Manifold is a mathematical surface (curve or point) defined locally by a system of equations through local "canonical" coordinates, also called curvilinear (natural) coordinates, or generalized coordinates of Lagrange (known as generalized coordinates for short). • Mapping is a transformation in which u = f(ξ, v) (3) v = g(ξ, v) and vice versa. Mapping is smooth if functions f and g are smooth (i.e., differentiable a "sufficient" number of times: the highest level of

"sufficient" differentiation is equivalent to the highest power of a polynomial describing a given manifold.

• A Monom is a point in Newton diagram space, describing a given polynomial term. For example, term: x³y is equivalent to the monom (3,1), a point in a Newton diagram (Figures 2A and 2B).

• Morse was a French mathematician whose work was a precursor of catastrophe theory. At the beginning of the nineteenth century, he generalized a number of differential geometry theorems into a general class of generic surfaces. • Morse lemma: In the vicinity of a nondegenerate critical point of specific manifold mapping, a function describing specific manifold mapping in generalized coordinates can be reduced to a quadratic form.

• A Newton Diagram is discrete (Cartesian) 2-D "point" space defined in such a way that the x-axis and the y-axis describe x- polynomial and y-polynomial power, respectively. For example, the x²y -polynomial element is equivalent to point (2,1) in Newton table space. In this Newton diagram space, a given polynomial that is always normalized (i.e., with unit coefficients of x² + xy, and not x² + 3xy) is described by a segmented line. See Figures 3 A and 3B.

• Nondegenerate Critical Point For this point, only one row of Jacobian is equal to zero (this point can be maximum or minimum, as referred to in the Morse lemma).

• Normal Form is a set of monoms bounded by a canonical form segmented line (including the monoms of canonical form).

• A Ring is the second most complex set in mathematical models, with two operations. Ring sub-sets can be field and ideal. • Stable Catastrophes are always two: fold and cusp (tuck). These cannot be "easily" transferred to another catastrophe by infinitesimal transformation although others can be. See Figures 1 and 2. • A Spectral Series is a method of sequential approximation

(proposed by Arnold) that allows reduction of all catastrophe- equivalent polynomials to the canonical form, representing a given type of catastrophe.

• Thom (Rene) was a French mathematician, considered to be the "father" of catastrophe theory (1959).

• Thorn's lemma is a fundamental theorem of catastrophe theory in general and the ISMP in particular, as a generalization of the Morse lemma for degenerate critical points. It claims that, in such a case, the algebraic form describing a surface can no longer be only quadratic, but consists of a quadratic form (as in the Morse lemma) and an additive singular residuum: s

/ = ∑ ^kA ⁺ S(c _s+n, ■ ■ ^■«„), *,- = ±1 i = π

These normalized coordinates are also separated into two parts: non-generate point coordinates (NPC) (i = l,2,... ,s) and degenerate point coordinates (DPC) (i = s + l,s+2,...,n). In the residuum function g(o:_{s+ 1},... ,o;_n), the first- and second- order differentials vanish: dg = d²g = 0. Functions with the same g belong to a set of stable equivalent functions, or are locally isomorphic (Arnold).

• The Thom Statement declares that there is a finite number of catastrophes (14) in 3-D space.

• A Vector Field is a representation whose element provides a shift of polynomials in the Newton diagram (this shift does not need to be a translation). • Whitney (M.) (1955) was an American mathematician, and a major contributor to catastrophe theory. His major achievements were in studying mapping from surface to plane.

• Whitney Theorem (Two Stable Catastrophes): The local normal form of the singularities of typical stable mappings from 2-D manifolds (in 3-D) to a plane can be either fold or cusp only. (Stable Mapping): Every singularity of smooth mapping of a surface onto a plane after an appropriate small perturbation splits into stable catastrophes only (fold and cusp). This theorem is applied in CMP hypercompression into still imagery.

The following references are referred to in the text that follows and are hereby incorporated by reference.

1. M. Born, E. Wolf, Principles of Optics, Pergamon Press, 1980.

2. T. Jannson, "Radiance Transfer Function, " J. Opt. Soc. Am. Vol. 70, No. 12, 1980] pp. 1544-1549.

3. V.I. Arnold, Catastrophe Theory Springer- Verlag, NY, 1992.

4. V.I. Arnold, Singularities of Caustics and Wave Fronts, Mathematics and Its Applications (Sovien Series) Vol. 62, Kluwer Academic Publisher, 1990.

5. V.I. Arnold, The Theory of Singularity and its Applications, Academia Nazionale dei Lincei, 1993.

6. V.I. Arnold, S.M. Gusein-Zade, A.N. Varchenko, Singularities of Differential Mapping, Birkhauser, Boston-Basel-Berlin, 1988.

7. R. Gilmore, Catastrophe Theory for Scientists and Engineers, John Wiley & Sons, New York, 1981. 8. P. Grey, Psychology, Worth Publishers, New York, NY, 1991.

The following are expanded definitions, theorems, and lemmas referred to in the discussion below:

Critical Point: For a function depending on n variables ( e Rⁿ or n- dimensional magnified), a critical point is called nondegenerate if its second differential is a nondegenerate quadratic form. In other words, for this point, only one row of the Jacobian is equal to zero. Noncritical Point: In the neighborhood of regular (or noncritical) point transformations of n local coordinates ^ of a surface into coordinates u_; on a mapping plane, the transformation can be written as:

Uj = Ui (^,... ^_n), i = l ,... ,n.

In this case, the Jacobian is always nondegenerate. This means that in the vicinity of this point, it is possible to do an isomorphical transformation according to the implicit function theorem:

£i = £; (u₁, ... ,u_n), i = l,... ,n

Morse lemma: In the neighborhood of a nondegenerate critical point, a function may be reduced to its quadratic part, i.e. , it may be written into the normal form

for a certain local coordinate system (ξ,, ..., £_n).

The meaning of this lemma is as follows: Since the Jacobian of any smooth function f is nonzero in the vicinity of any nondegenerate critical point, differential replacements, such as:

U; = u_i (ξ₁, ..., £_n) (2) can transform this function into a nondegenerate quadratic form:

^«. ⁼ i Σ = 1 */?/ ^{k = ± (3)}

At a degenerate critical point, some eigenvalues of the Jacobian matrix are zero. The subspace spanned by the corresponding eigenvalues (£_s+1, ..., ξ_n) is a critical subspace and has a dimension equal to the co-rank of the critical point.

Function f can be written in the form defining Thorn's lemma, which is fundamental to the ISMP:

where g(£_s+1, ..., £„) is a function (residuum) for which dg = d²g = 0. All functions with the same g are called a differential equivalent. The term local isomorphical is used as another description of that class of function.

Thorn's lemma provides a basis for an application mapping algorithm for any surfaces that can be mapped on an image plane. Thorn's form can be used as a nondegenerate function for image approximation, but using singularities analysis allows extraction of the most important information from an image.

1st Whitney Theorem. The local normal form of the singularities of typical mappings from two-dimensional manifolds to a plane (or to another two-dimensional manifold):

cusp (" ^{= ζ 3 + ξ η} (5A) v = η

- L fold [ (5B) v = η

regular {

° v = η

Stable catastrophes (fold and cusp) are sufficient for still image compression. Fig. 3 A depicts a fold and Figure 3B depicts a cusp.

2nd Whitney Theorem: Every singularity of a smooth mapping of a surface onto a plane after an appropriate small perturbation splits into folds and cusps. Arnold Theorem: There is a family of transformations that can transform any given mapping into a set of canonical transformations by using smooth substitution.

(ξ, τ?) → (u, v) u = u (ξ, η) (6)

V = V (ξ , η) then, by using smooth diffeomorphic transformation into new "plane" coordinates (u\ v'), we obtain:

u' = a^ + a₂v + a₃u² + ... v' = bjv + b₂u + b₃v² -1- ... (7A) u" = c_tu' + c₂v' + c₃u'² + ... v" = d,v' + d₂u' + d₃v'² + ... (7B)

We can obtain

u = M, (£", V") + F_C1 (£", η") + F_C2 (ξ", v") + ... (8) v = M₂ (£", η") + _C3 (ξ η") + ... where M, and M₂ are Morse forms and F_C1, F_C2, F_C3, ... are canonical (singular) forms. In the Arnold theorem, Thorn's lemma is applied in such a sense that we represent a Thom form (as in Eq. 8) by supeφosition of Morse (smooth) forms (M), and Thom residual forms (F_c). The proof of this theorem is based on spectral series reduction to normal forms. This is a local isomorphism (or, more precisely, local diffeomorphism) , because each catastrophe is represented by a given canonical form, which, in turn, generates a normal form. Moreover, each catastrophe is represented by only one canonical form. Therefore, while general mapping is usually not isomoφhic, in this specific case, Arnold mapping is. The consequence of the Arnold theorem, proven by his students Platonova and Shcherback, - is a statement made earlier by Thom:

77ιe number of nonequivalent singularities in the projections of generic surfaces in 3- D space, defined by the families of rays issuing from different points of space outside the surface, is finite and equal to 14.

Physical Modeling by Catastrophic Manifold Projection Smooth surfaces vs. image presentation

Usually, 3-D objects, presented in the form of 2-D images, are projections of the following types of objects: 1. Smooth artificial and natural objects: This category can be described as a projection of idealized surfaces on the image plane. In accordance with the human visual system, we first try to extract objects that can be presented by smooth surfaces ("soft edges"). Edges of smooth surfaces: These soft edges, which appear during mapping, are the same as the visible contours of smooth objects. (These objects ideally fit into the proposed approach.)

2. Sharp joints of objects: One example is the corners of buildings. (These jumps can also be naturally described by the proposed approach.)

3. Textures of an object: These objects will be described by the proposed method with natural scale parameters (including fractal type textures).

Physical Model Formulation

Formation of an image can be described as light reflection from a general surface. It may be an actual radiation surface (light source, transparent surface, or semi-transparent surface) or it may be an opaque reflected surface. We have introduced a photometric projection, so each ray is reflected backward only, in accordance with the radiance projection theorem ^[2]. The reflection is the highest in the specular direction, and it is monotonically reduced, with an increase of the reflection direction separation from the specular direction, as shown in Figure 4. This photometric projection approach can be derived from Thorn's lemma and the Arnold theorem. If a reflection is identified with reflection surface luminance, B, the dependence of B on the direction will depend on the nature of the surface (whether it is smooth or rough). There is no general theory for arbitrary surfaces, although there are two limited cases: Lambert's cosine law, in which B is a constant (isotropic case), and the specular (mirror) reflection, in which incident light is reflected, without distortion, only in the specular direction. In a general case, we have intermediate distribution as seen in Figure 4 A which is a reflection from a manifold, 4B which is an explanation of reflection value depending upon angle θ, and 4C which is what would be displayed, for example, on a display. The presented photometric projection has a natural inteφretation: the reflection value decreases when the θ- value increases, and vice versa. Catastrophic Manifold Projection (CMP)

Inverse projection from a 2-D plane into the surface of a real object is analyzed. The geometry of this problem (i.e., photometric projection) has been shown in Figure 4. Now, however, this inverse problem must be formulated in precise mathematical terms, allowing design of a suitable algorithm for hypercompression. To do this, the forward problem of image formation is formulated first.

In general, image formation can be presented as differential mapping from 4-D space (x, y, z, B; where x, y, z are real 3-D space coordinates of a point and B is the luminance of that point) to the image plane (u, v, w), where u and v are coordinates of a pixel and w is a color (or gray scale level) of the pixel. The result of mapping the manifold with internal curvilinear (ξ, η) coordinates will be a 3-D surface (u, v, w); where u, v are coordinates of the point into the plane and w is the luminance of the point.

The mapping will be:

u = f j (x, y, z) (9A) v = f₂ (x, y, z) (9B) w = f ₃ (x, y, z, B). (9C) or u = F, (ξ, η) (10A) v = F₂ (ξ, η) (10B) w = F₃ (ξ, η, B) (IOC) where fj and f₂ are regular projections of a surface to a plane, f₃ is luminance projection, and F F₂, F₃ are their equivalents in the curvilinear coordinate system tt, ).

To formulate the isomoφhic singular manifold projection (ISMP) problems by applying Thom's lemma formalism (i.e., canonical form, catastrophes, etc.), one must realize first that the w (B) - dependence is a smooth monotonic one, since both w and B are various forms of luminance, in such a sense that B is the physical luminance, while w is its representation in the form of color (gray level) in the CCD plane. But, smooth dependence does not contain critical points (even nondegenerate ones). Therefore, Thom's (splitting) lemma can be applied to Function (10), in the form: w = M (B, ξ, η), or (11 A) w = M (B, + g (TJ), or (11B) w = M (B, η) + g ( , or (11C) w = M(B) + g (ξ, TJ), (11D) where the first function M represents a monotonical function of B or (function without critical points), and g (£, η) represents all singularities of projection influencing a gray scale level (color) of a given point (i.e. , g-function represents a singular Thom residuum). In order to show this, Function (IOC) is expanded into infinite Taylor series, in the vicinity of ₀, η₀, and B₀, in the form:

1 dF~ ι δ % ι ΘF, _w = _w + _ _ 2 1 (ξ - ξ_) + ± — 1 \ (η - η₀) + — — 1 \ (B - BA ⁰ 1! dξ ^ξ° °^ 1! dr\ ^η° ° 1! dB ^B° °

1 3^-F, 1 F^ + J_ 1 | (ξ - ξ ) B - BA + — | (_η - _η )(β - H ) +...

It should be noted that neither linear form of this Taylor series can be singular, by definition, and, therefore, neither is of interest. In relation to quadratic terms, coordinate substitution will be provided so that after this substitution, some free coefficients will be received that permit the zeroing of mixed quadratic terms (this approach is completely within the framework of Thom's lemma proof). On the other hand, the quadratic term (B-B₀)² is a Morse term. Therefore, it is demonstrated that there are no singular B-dependent terms. In summary, a luminance physical coordinate does not introduce new singularities, and, because g depends only on geometrical coordinates ( , TJ) (belonging to 3-D space manifold), all previous results of the Whitney-Thom-Araold theory apply in this new geometrical/physical ("geophysical") 4-D space. (See Table 1.)

Table 1

In order to explain these new results, a simple example of a homogeneous object with constant luminance B is considered.

Example 1: Arbitrary object with constant luminance

In such a case, Eq. (IOC) does not contain B-dependence; i.e., it can be written in the following form: w = F₃ (ξ, TJ) (13) Now, the first two equations (10A) and (10B) can be used without changes, to introduce (u, v) - coordinates, in the form:

w = F(u(ξ,η), v(ξ,η)) (14)

where F is some new function, and

w = F₃ (£, TJ) . (15) It is clear that w-coordinate should have the same singularities as u and v

(see Table 1).

In this case, all changes in color (gray level) will be determined only by mapping F and the contour of an object. Of course, the singularities for color w will be located at the same points as singularities of u and v. As a consequence, the singularities of color will be displaced at the contour of the object. Example 2: Cylinder with given luminance dependence Mapping of a cylinder with a given constant luminance dependence is shown in Figure 5 and described as follows:

B = f (ξ, TJ) . (16) In a cylindrical coordinate system (where axis y coincides with the axis of the cylinder), x = , where is angle Z BOA, and z is distance OB (or, radius).

Two parametric coordinates, £ = , where is angle Z BOA (A is the central point of cylinder, B is a given point); y is the axial coordinate, and z (=R, where R is const) is the radius vector (OB). That the w-parameter must be proportional to B, everywhere must be taken into account. This means that B does not create any singularities. For new coordinates on the image plane: u = R sin (£) (17 A) v = TJ (17B) w = C B + f(£, TJ); C = const ≠ 0 . (17C) On the other hand, a geometrical analysis of transformation Eq. (9A) shows that y- does not produce any singularities (since y is an axial coordinate). Therefore, it can be assumed, without loss of generality, that w depends only on x and B in the following form: w = C -B + f ( u(ξ) ) (18) or w = C B 4- f'(x) (19) where f'(x) = f (u(x)).

For critical point estimation, the Jacobian is considered, transforming coordinates (£, TJ, B) into (u, v, w) in the form:

or

The first row of the matrix is not equal to 0, except

π x = ±— 2

Therefore, it is possible to use a smooth transformation between (x, y, z) and (u, v, w).

(ξ, B) - (u, w) (21)

and there are no singularities for

π

\x \ ≠

but if π x - ±—

2

the first row becomes 0 and the determinant equals 0. This means that in the case of

π x - ±—

2

we cannot perform smooth variable substitution, and singularities exist in these points (projection of fold). Let function w be represented in the expanded Taylor series:

The significance of a nondegenerate point

π x = ±-

(fold) becomes clear if it is realized that even in the case of v, a weak dependence between w and x,

du dx grows to infinity (in the vicinity of

x = ±— ). 2

Because the singularity appears as a result of geometrical mapping (not connected to changes of color), and we assume that the color of an object is a smooth function of coordinates x_l5 x₂, it is possible to use a canonical function for representation of function f ' : f'(x_v x₃) = Bx₃ + c-F(X_j)

where F is the deformation ^[3] of a canonical polynomial (fold or x'^Λ type dependance).

The calculations presented above do not use rigorous mathematical calculations, but they are very close to Lev Landau 's approach, applied successfully to many areas of theoretical physics. In his approach, the art of throwing away "inessential" terms of the Taylor series, and preserving smaller size, yet "physically important" terms, has been rigorously proven through the course of the catastrophe theory ^[3].

Drawbacks of Fourier Analysis Describing an arbitrary function by using a standard transform, such as

Fourier or wavelet, is natural for periodic signal analysis. In image processing, however, these approaches have difficulties with describing very high redundancy regions with flat, slow-changing parts, as well as regions of abrupt change (or "soft edges"). Such classical description is unnatural for these types of objects because it creates excessively high input values in almost every coefficient of the Fourier transform as well as large coefficients in the case of the wavelet transform.

At the same time, starting from Leibniz, Huygens, and Newton, a clear geometrical (polynomial) approach was developed for an analysis of smooth curves and surfaces. As discovered recently, this approach has become strongly related to many major areas of mathematics, including group theory, regular polyhedrons, wave front propagation (caustics), and dynamic systems analysis. For a clear demonstration of the unique properties of this approach, consider the classic evolvent problem, formulated in Newton's time: For example, for f=x³, the evolvent presented in Figure 6 can be constructed. Arnold ^[6] has shown that the evolvent is directly related to the H₃ group generated by reflections of an icosahedron. (H₃ is a group of symmetry of the icosahedron.) H₃ has special properties, as described below. If complex space C³ instead of R³ is analyzed, the factor-space of C³ for this group will be isomoφhic to C³. This means there exist some basic polynomial invariants. By using these invariants, any polynomials of this group can be represented (Arnold ^[6]). To illustrate this property in 2-D, let us describe a simplified example of three mirrors on R² as seen in Figs 7 A and 7B.

The points of a plane that have an equal number of reflections (12 in Figure 4 A) belong to one (regular) orbit. Points located on the mirrors belong to another orbit. A set of all irregular orbits in a factor space is a discriminant (i.e., the manifold in a factor space). Now the plane in 3-D space can be represented, as in Figure 4, as a plane with coordinates z_l5 z₂, z₃. The plane can be determined by:

In this space, it is possible to introduce permutation of the axis, generated by reflections. z, = _Zj (24)

Orbits in this context constitute a set of numbers {z_{l t} Zj, Z₃} (with all permutations generated by reflections), with the additional condition of Eq. (23).

This unordered set will be uniquely determined by polynomials:

z³ + λ_jz² + λ₂z + λ₃ = 0 (25)

By using Eq. (25), the following is obtained:

λ_j = 0

or

The space of the orbits of this group will be naturally presented by the roots of a cubic polynomial Eq. (26). This means that in factor space, this space is just a plane with coordinating (a, b). Each point (a, b) of this space corresponds to a cubic polynomial and its roots.

If some of the roots are equal, that means we have received irregular orbits. The discriminant in this case is

4a³+27b² = 0 (27) which is of ³/₂ type curve. This curve corresponds to specific orbits (in the mirrors) in Figure 4 A.

For all other types of groups generated by reflection, analogical construction of the discriminant exists. It can be proven (Arnold; see Ref. [6]) that the surface creating the evolvent is diffeomoφhical to the discriminant of the H₃ icosahedron group. As a result, by using the group representation, the redundancy of a mathematical object that is diffeomoφhical to our mapping procedure has been greatly reduced.

Taking into account the symmetry of 3-D objects mapping (symmetry in a general sense, this means the Lie group, in this case) can optimally reduce redundancy and extract the information that describes the most important features of our object.

In summary, the polynomial representation of geometrical objects (starting from: Newton through Bernoulli, up to Thom and Arnold) seems to be more natural than the common Fourier (and wavelet) approach, because polynomials are connected to groups of symmetry that permit reduction in orbit redundancy in a most natural way.

Catastrophe Theory Applied to Still Image Compression

Because the most critical part of an object — its 3-D boundary — can be described by a 1-D contour and three or four natural digits or "coefficients" that characterize a simple catastrophic polynomial, tremendous lossless compression of object boundaries can be achieved, far exceeding state of the art compression ratios while still preserving high quality image. Since in all state of the art still image compression methods the major information loss is at the boundaries, applying ISMP compression which actually preserves the boundary or edge information, provides unparalleled fundamental compression ratio/PSNR trade off.

"Catastrophe" or alternatively isomoφhic singular manifold as used in this patent designates a mathematical object that describes the shape of 3-D object boundaries in polynomial form. The use of "catastrophe" theory for compression makes the present invention unlike all other compression methods because it helps to transmit information about 3-D object boundaries without loss, preserving the features of the object most valuable to human cognition, but with very high compression rate. By applying the present invention to still image compression, a 300: 1 still image compression ratio with practically invisible artifacts (PSNR equal 32 dB) and a 4,000: 1 full motion image compression ratio with fully developed natural motion and good image quality, is obtained.

The still image compression and related video compression technique of the present invention is extremely beneficial because, unlike other state of the art compression techniques, major information losses do not occur from compression at 3-D object boundaries (edges) that require both high dynamic range and high resolution (i.e., both high spatial and high vertical: "Lebesque" resolution). In these edges there is a vast amount of information necessary for many processing operations vital to a quality image and human cognition. The compression technique of the present invention, unlike other compression methods, preserves intact all "soft- "edge information without data loss. Hence, the present inventors have coined the term lossless-on-the-edges (LOTE) compression. LOTE compression is possible because of the fully isomoφhic projection between the 3-D object boundary vicinity and its 2-D projection on the screen. This fully isomoφhic projection between the 3-D object boundary and the 2-D projection is based on Arnold's so called "catastrophe" theory that has been adapted to still image compression here. The methodology of the present invention works especially well with objects that are closer to sculptures and objects that have mostly flat surfaces combined with edgy features, i.e., very low or high spatial frequencies. This is exactly the opposite of Fourier analysis which does not work well with very low or very high frequencies. For low frequency, Fourier analysis is unsatisfactory because the coefficients must be very well balanced and at the same time can easily be hidden in noise. For high frequency components of objects such as edges, many high frequency components exist which Fourier methods eliminate. These high frequency components are what make up all important edges of the object and eliminating them reduces human cognition. ISMP analysis on the other hand, does not have this problem because it characterizes edges and objects using manifolds and hence preserves the information that makes up those edges and that was eliminated in Fourier-based compression methods. Specific Features of the Human Perception of Visual Information and Object Recognition

An understanding of how humans recognize objects will make manifest the advantage of preserving information. The retina of the human eye contains millions of receptor cells, arranged in a mosaic-like pattern in the retinal layer. The receptor cells are cones and rods. These cones and rods provide the starting point for two separate but interacting visual systems within the human eye. Cone vision is specialized for high acuity and for perception of color. Rod vision is specialized for sensitivity and the ability to distinguish color (i.e., a person can make out the general shape of the objects, but not their colors or small details^[8]).

The main puφose of human vision is not to detect simple presence or absence of light, but rather to detect and identify objects. Objects are defined principally by their contours. The visual system registers a greater difference in brightness between adjacent visual images, faithfully recording the actual physical difference in light intensity.

David Hubel and Thornton Wiesel (Nobel prize winners in 1981) recorded the electrical activity of individual neurons in the visual cortex. They found that these cells were highly sensitive to contours, responding best not to circular spots but rather to light or dark bars or edges. They classified these cells by using a complex hierarchical system, based on their different response characteristics. In this research, the authors outlined that the perception of long and linear, bars provided maximum response in the human visual system.

Human brain zones, which decode specific properties of image recognition, are spatially organized in the brain according to their function. Thus, different localized sets of neurons in the visual cortex are specialized to carry codes for contours, color, spatial position and movement. This segregation of functions explains why a person who has had a stroke, which damaged part of the cortex, sometimes loses the ability to see contours without losing the ability to see colors.

Special mechanisms of object edge extraction in the human visual system allow extraction of important objects from a background, even if the object has bulk colors very close to the colors of the second plane. The latter feature is extremely important for registration of military targets, and makes ISMP an effective compression algorithm for ATR.

This still image compression performance can be transformed into analogous video image compression through the typical 10:1 factor for state of the art video image compression. Therefore, the inventive technique can be applied not only to high resolution digital video/still image transmission, but also to multi-media presentation, high quality video conferencing, video servers, and the storage of large amounts of video information.

Catastrophe Theory Applied to Video Compression Video compression is a four dimensional (4-D) problem where the goal is to remove spatial and temporal redundancy from the stream of video information. In video there are scenes containing an object that continuously changes without jumps and has no edges, and, on the other hand, there are also scenes where there are cuts which are big jumps in the temporal domain or big jumps in the spatial domain (such as "edges"). These abrupt changes or jumps can be described as "catastrophes. " Using catastrophe theory, these behaviors can be described by one or more elemental catastrophes. Each of these elemental catastrophes describes a particular type of abrupt change in the temporal or spatial domains. In general, categorizing catastrophes in 4-D space is even less established than catastrophe theory in general which is relatively unknown. Furthermore, 4-D space is far less understood than 3-D space but, similarities between them can be expected and can use projection-type mapping can be used, but in temporal space. One solution is to use spatial catastrophes along with temporal catastrophes.

In order to apply catastrophic theory to video imagery, a fourth "geometrical" coordinate, time leading to time-space (4-D) is preferably added. In the case of the inventive isomoφhic singular manifold projection (ISMP) methodology, five- dimensional (5-D) geometro-physical space (x, y, z, t, B), where B is brightness, or luminance, is obtained. This 4-D time-space (x, y, z, t) plus physical coordinate, B, can be split into 4-D geometro-physical space, and time (t) and treated separately except in the case of relativistic velocities. In the latter, relativistic case, the 5-D space can be analyzed by Poincare group formalism. In the common, non-relativistic case, however, temporal singularities (catastrophic) may be described in the time- luminance (t, B) domain only. The time-luminance singularities may interfere with spatial singularities (previously discussed). In such a mode of operation, each block of the image is represented by a single time- variable value. According to Figure 8, there are only two possible singularities describing any type of mapping including smooth curve projection designated (1) shown in Figure 8 where < B > is the average B-value characterizing a frame as a total structure (the smooth projection shown in Figure 8 represents movement of a physical object. Item (2) in Figure 8 representing a catastrophic frame change, and item (3) representing position/ tilt/zoom camera changes. The critical < B > -parameter may be, for example, an average block-to-block error (e.g. mean square error). In summary, temporal catastrophic formalisms can be applied to MPEG hypercompression by replacing the average error parameters by integrated luminance-changes.

Canonical Polynomials One way to represent these 4-D catastrophes is to use well known 3-D projections or mapping catastrophes which were discovered in the early 1980's. These "transformations" or "reconstructions" or "metamoφhoses" in time are 4-D problems which can be separated into two 3-D problems: 1) Spatial catastrophes may be defined in 3-D space (x, y, B) such as occurs when there is a large change in intensity, B, over a small change in x, y; 2) temporal catastrophes may also be defined temporally such as occurs where there is an abrupt change in motion over time such as is present during the rotation of an object or a cut from one scene to another. The 3-D temporal problem can be further reduced to a 2-D problem by transferring the (x, y, B) coordinates into 1-D merit space. Merit space is defined by the lack of similarity between frames in time.

Images are 3-D distributions of intensity. Abrupt changes in intensity occurring over small changes in x, y may be treated as catastrophic changes. The inventors have modified catastrophe theory to fit images and to solve the problems of image and video compression. The inventors have introduced a physical coordinate, B (luminance) into conventional geometrical coordinates to create "geometro-physical" surfaces. There exists a finite list of fourteen polynomials or "germs" which describe different edge transitions or projections in mapping in 3-D space. Typically, only about three of these polynomials or germs are necessary to describe virtually every edge effect. The others are used on occasion to describe spatial projections.

The germs of the projections are equivalent to the germs of the projections of the surfaces z = f(x,y) along the x-axis. The table below identifies the fourteen polynomials or germs.

In theory, a projection of a surface does not have any germs that are inequivalent to the fourteen germs in the above table. It should be understood that the Spectral Series for Reduction to the Normal Form (SSRNF) method is used for the unique reduction of the arbitrary polynomial to the germs presented in the above table. It is presented here only in a descriptive form:

Let e . . . , e_n - quasihomogeneous polynomial (N+p degree) that generates τ + l

for diffeomoφhism.

A → ^{τ + 1 A}P ^ΛP

Then, there is formal diffeomoφhism

and that the series f = f₀ + f , + after substitution has a form and Cj represent the numbers. y_n = ^χ _n ⁺ s_n , Ε s„— ^{e A}p-T

A^yv~>^y = fob ⁺f⁽.^{χ) +}- ⁺f_P-ι ⁺Σ ^ci^et^{(χ) +R}> ^{R e Λ}p+ι

Catastrophe theory has not been before used for studying image intensity because the number of coefficients necessary to satisfactorily describe an image using standard polynomials is simply too large and can exceed the number of pixels present in an image. Obviously, such an analysis is not worthwhile because the data that need to be handled are larger than the number of pixels, itself a very large number. The inventors have discovered that it is possible to remove many of the details or "texture" in images, leaving the important "sculpture" of the image, prior to characterizing the image with polynomials, to significantly decrease the number of coefficients in the polynomials that describe the different edge transitions in mapping and 3-D space.

Preferred Still Image Encoding Method

The following is an abbreviated description of the still image compression method as in the flow chart of Figure 9. Step 1 involves segmenting the original image into blocks of pixels, for example 16 x 16. Step 2 is to create a model surface for each segment or block corresponding to the original image so that there is isomoφhism between the original image and the polynomial surface in accordance with Arnold's Theorem. More particularly this may involve calculating the equation ^_modelled f°^{r eacn} block or segment by substituting for variables in canonical polynomials. (See Steps 3-7 of detailed flow chart which follows). This step inherently eliminates texture of the image and emphasizes the "sculpture" characteristics. Step 3 is to optimize each model segment. This is done by calculating the difference between the original and model segments and choosing coefficients for the canonical polynomial which have the lowest Q i.e., the smallest amount of difference between the original segment and the modelled segment. This is repeated on a segment by segment basis. ( See Steps 8-12 of the detailed flow chart which follows). Step 4 is to find connections between adjacent segments to create an entire image i.e. , a model image of the entire frame. ( See Steps 14-18 of the detailed flow chart). This yields an entire image that has only the "sculpture" characteristics of the original image and eliminates texture. Step 5 is to calculate the peak signal to nose ratio PSNR over the entire image and where the PSNR of the entire image is less than a threshold, the difference between the original image and the modelled image is calculated. This step recreates the texture information of the original image that was lost during the process. Thus, after this step there are two sets of data: the

"sculpture" characteristics represented by a few discrete numbers or "datery" and the texture information of the image. (See Steps 19-21) Step 6 is to use standard lossy compression on the texture portion of the data and then to combine the texture and datery and apply standard lossless compression to that combined data. (See Steps 22-

24 of the detailed flow chart).

Now the preferred still image encoding method will be described in detail in relation to the detailed flow chart. In the following description of the still image encoding process according to the present invention the following definitions are used: I₀ = original I frame I_m = modeled I frame

I_d = difference I frame (I_d = I₀ - I_m for each frame) i₀ = segment or block of original frame i_m = segment or block of modeled frame i_d = segment or block of difference frame (i_d = i₀ - i_m for each block)

Referring now to Figure 10, this figure sets forth a flow chart of the still image encoding process according to the present invention. In step 1, the next still image or I₀ frame is captured. If only still images are being compressed for still image puφoses, this image will represent one of those still images. If video is being compressed, the still image to be compressed here is one of the video's I frames which will be compressed in accordance with this method and then inserted at the appropriate location into the video bitstream. In step 2, the original image L, is segmented into blocks of pixels of any desired size such as, for example, 16 x 16 square blocks. The original image is seen in Fig. 11 A. Any segment size may be used as desired. These segments or blocks of pixels are designated i,,. This segmentation is done according to standard segmentation methods. As an example, the total number of segments or blocks for a 512 x 512 image is 512 x 512/(16 x 16) = 1024 different noninterleaving 16 x 16 segments or blocks. Step 3 is the first step involving segment by segment operation on each i₀ using matrix representation of each segment. In step 3, the Dynamic Range (R) of each segment or block is calculated using the following equation:

R = max(i₀) - min(i_o) In the above formula, the pixel having the maximum intensity is subtracted from the pixel having the minimum intensity in the segment. This difference is the Dynamic Range R.

Step 4 compares the Dynamic Range R to R_Q which is a threshold determined from trial and error. The threshold R₀ is chosen so as to eliminate unnecessary compression such as compression of background scenes. In this regard, if the value R is very small and less than R_<„ the image is most likely background and the compression technique of the present invention is not needed. In this case, the process is started over again between steps 2 and 3 and another segment or block is operated on. If R is greater than or equal to R₀ then the subsequent steps involved in choosing a canonical polynomial from the table and creating a model polynomial by solving its coefficients are then performed. This set of steps now generally described involves choosing the polynomial from the table which best matches each particular segment or block.

Turning now to step 5, a first canonical polynomial from the table is taken.

In step 6, substitutions for variables in the canonical polynomials are found. It is possible to apply (1) a nonhomogeneous linear transformation (shift of coordinates),

(2) a homogeneous linear transformation (rotation of axis) or (3) a nonhomogeneous nonlinear transformation. For example, if the canonical polynomial

J f canonical ~ ^xl ^{3 + X}\^X2

is taken from the table, variables _x and x₂ are substituted for as follows using the third example above nonhomogeneous nonlinear transformation: *ι = (y_ ⁺ øiV? + • . . l); x₂ = (y₂ + fc.y₂ ² + . . . b„y²).

From this substitution a function describing a "modeled" surface (as opposed to the original image surface) is generated as follows:

f model

= yj +

+ ay* + α³yf + 2α²y + y_χy₂ + aby_χy₂ + ay_χy₂ +

= Vi ⁺ αyf)³ + (y, + flVj)(y₂ + by₂)

At step 7, the modeled surface is created by substituting the coordinates of each pixel in the original segment or block into the equation f_modeι. A modelled surface is seen in Figure 11C. This creates a matrix containing the values f_m(u₎, ^x _m(i_,2) • • • as seen in Figure 12. Specifically, this matrix is created by substituting the coordinate of the pixel 1 , 1 from the original segment into the equation f_model to generate the element f_m(1>1) in the modeled matrix. Next, the coordinate of the pixel 1,2 from the original segment is substituted into the equation f_model to generate value _m(i_,2) which goes in the 1,2 pixel location of the modeled surface. This is done for each pixel position of the original segment to create a corresponding modeled matrix using the equation f_model.

At step 8, Q is calculated by determining the difference between the original and modeled segments, pixel by pixel, using the equation:

In other words, Q is calculated by subtracting corresponding pixels from the i₀ segment (the original segment) from the i_m (the modeled segment) and squaring this subtraction and summing up all these squares.

At step 9, Q is compared to a predetermined threshold Q₀ based on image quality desired. Q should be less than Q₀ because the point of the step 8 is to minimize the sum of the differences between the analogous pixels in the original and modeled frames so as to generate a modeled surface that is as close as possible to the original surface. If Q is greater than Q₀, the procedure loops back up to step 6 where new coefficients are tried in the same polynomial. Then steps 7, 8, and 9 are repeated, and if Q is less than Q₀ with that new set of coefficients, then the process continues into step 10 where that Q and the coefficients that produced the lowest Q for that polynomial are stored. After storage at step 10, the process loops back up to step 5 if all polynomials have not yet been tested where the next canonical polynomial from the library is chosen and tested and solved for coefficients which produce the lowest Q for that next polynomial. Hence, steps 6, 7, 8, and 9 are repeated for that next polynomial until coefficients are found which produce the lowest Q for that polynomial. At step 10, the Q and the coefficients for that next polynomial are stored. This process of steps 5, 6, 7, 8, 9, and 10 is repeated for each polynomial in the library. After each polynomial in the library is tested for the segment under test, the process moves to step 11.

At step 11, the polynomial having the lowest Q of the polynomials tested for that segment is chosen. That polynomial is transferred to step 12.

At step 12, all coefficients for the chosen polynomial (the one having the lowest Q of all the polynomials tested for that segment) are stored. These coefficients are coefficients of the equation f_modeI which describes the modeled surface.

After step 12, the next set of operations involves segment by segment operation working only with the polynomials and their coefficients whereas the above steps 5-12 worked with the matrix representation of each segment. Because only the polynomials and their coefficients are worked with in the next set of operations, a significant amount of compression has taken place because the data representing the surface is far less voluminous than when a matrix representation of the segments is worked with. The data is simply coefficients of polynomials which can be called "datery".

At step 13, the current segment is taken or captured from the above steps. At step 14, a connection is found between adjacent or neighboring segments by extending the surface of a first segment into a second segment and finding differences between the extended surface and the second segment surface. Specifically, this is done by finding the average distance "q" between the surface which extends from the first segment into the second segment and the surface of the second segment using standard methods. If the average distance "q" is smaller than a threshold value q₀, the surface of the second segment is approximated by the extended surface. In other words, if the distance q is smaller than the threshold value q₀, the second segment surface is thrown out because it can be approximated satisfactorily by substituting the extended surface in its place. If the average distance q is greater than the threshold value q₀, a connection needs to be found between the extended surface and the surface of the second segment.

Thus, at step 15, the average distance q is checked to determine whether it is less than the threshold value q₀. If it is, then the connections between the adjacent or neighboring segments which can be plotted as a graph, as seen in Figure 13, are stored on a segment by segment basis. In other words, as seen in Figure 13, the surfaces which extend from, say, a surface in segment "9" into adjacent or neighboring segments (8, 10, 14 and 15), if any, are stored in the polynomial for segment 9 (earlier calculated and then stored at step 12) which then represents that graph of connections between segment 9 and segment 8, 10, 14 and 15. In other words, the polynomial that was calculated and stored for the segment in question, here segment 9, is modified so that it now extends into adjacent segments 8, 10, 14, and 15 and represents the surfaces in those segments. The polynomials for segments

8, 9, 10, 14 will be substituted with the new bigger scale polynomial obtained from 9.

If the average distance from 9 was not less than q₀ (which indicates that the surface extended from the segment in question, segment 9 for example, into an adjacent segment, 8, 10, 14, or 15 for example, did not satisfactorily approximate the surface of the second segment), then a spline must be calculated at step 17. At step 17 splines are calculated from the segment with adjacent segments using standard spline equations which need not be detailed here.

After both steps 16 and 17, the process continues at step 18. At step 18, a model image i_m is created of the entire frame by creating a table of all segments for that frame using the information calculated for each segment in the above steps. The creation of this table representing the entire frame from its numerous segments is analogous to step 7 where a modeled segment was created by substituting the pixel coordinates from the original segment into the f_mode, polynomial to get a matrix describing the modeled surface. At step 18, however, instead of creating a modeled segment of pixels, a modeled frame is created from modeled segments. Thus, it can be seen that the smaller parts calculated above are now being combined to generate an entire modeled frame, After step 19, the peak signal to noise ratio (PSNR) is calculated over the entire image using the equation:

10 log,_n (number of gray scale levels)² PSNR = ^- -A- -

Λf-1 N-i

Q = f hv ∑ ∑ U₀(^χ,y) - ι_m(^χ,y)Ϋ ι = 0 y = 0

h and v are number of pixels in horizontal and vertical directions respectively for the entire frame image. The Q values for each of the segments were stored at step 10 above and may be retrieved for this puφose. At step 20, the PSNR of the entire frame is compared to a threshold P₀. If

PSNR is less than P₀, then no further processing according to the present invention need be accomplished and processing can continue at step 24 where a standard lossless compression such as Hoffman encoding and run-length encoding are used to further compress the frame data. The compressed data is then sent to storage or a communication link.

If PSNR is greater than P₀ at step 20, then processing continues at step 21. At step 21, the difference between the original frame I₀ and the modeled frame I_m is found by subtracting each pixel in I_m (which was created at step 18) from the corresponding pixels in I₀ and a new frame I_d is created (see Figure 14) where each pixel in that frame has as its value the difference between the corresponding pixels in the frame I₀ and the frame I_m. The frame I_d therefore corresponds to the high frequency components, such as edge information which typically is lost in conventional compression techniques. This "texture" information containing high frequency components and edge information is then compressed separately at step 22. At step 22, standard lossy texture compression of the newly created frame I_d is performed by using standard methods such as DCT, wavelet, and fractal methods. At step 22, standard additional lossless compression is also performed. The output of step 22 is I_d' which then is fed into step 23. At step 23, the I_m frame is stored and the I_d' frame is stored. This concludes the compression of the still frame or I₀ frame.

As can be seen, the polynomial surface image is highly compressed because it is stored and transmitted as a complex algorithm (polynomial) rather than as a matrix representation. Additionally, the edge contour image I_d is separated from the polynomial surface, as a by-product of characterizing the original image by a canonical polynomial and contains the high frequency and edge components and is itself compressed.

Preferred Still Image Decoding Method

The still image decoding process will now be described as seen in the flow chart of Figure 15. The input to the still image decoding process will be either just the whole frame I_m in the case where the PSNR of the whole frame at step 20 was not less than threshold P₀ or the whole frame I_m plus I_d' where the PSNR of the whole frame at step 20 was less than P₀ and the differences between the original I₀ frame and the modeled frame I_m were calculated to create new frame I_d holding the textured or high frequency and edge information. In either case, the first step in decoding is step 1 which decodes the lossless compression data from the encoder which was compressed at step 24. At step 2 of the decoding process, frame I_m is separated from the other data in the bitstream. At step 3, the first graph or segment which was stored at step 16 on a segment by segment basis is taken. At step 4, whether the segment belongs to a graph (i.e., has connections to adjacent segments) or is an isolated segment (i.e. , has no connections to neighboring segments) is tested. If the segment does belong to a graph, then at step 5, a segment i_m is constructed for each graph (analogous to the creation of the modeled matrix surface in step 7 of the encoding process) using the polynomial that was stored for that graph at step 16 of the encoding process. If the segment does not belong to a graph, then after step 4 the process skips step 5 and continues with step 6.

At step 6, the separate graphs using standard splines are connected. In other words, those segments from steps 14, 15, and 17 which were connected by splines will be reconnected here. (Recall that it was these segments for which the extended surface of another adjacent segment did not satisfactorily characterize the surface of these segments and therefore a spline equation had to be used.)

From step 5 where a segment i_m for each graph was reconstructed, and from step 6 where separate graphs were connected using standard splines, the process continues at step 7.

At step 7, the frame I_m is constructed using segments from step 5 (in the same way as the frame I_m was constructed in the encoding process at step 18 and also similar to how an individual segment or modeled surface or modeled segment was created at step 7 in the encoding process.) At step 8, the presence of a I_d' frame for or in conjunction with the frame

I_m is tested. If there is no frame I_d', then the process is finished and the still image is fully decoded for that frame. If on the other hand there is a frame I_d' in conjunction with the frame i_m, then the process continues to step 9. At step 9, the frame I_d' is decompressed. At step 10, frame I₀' is created from the combination of frame I_m from step

7 of the decoding process and frame I_d' from step 9 of the decoding process. After the frame I₀' is created, the process is finished and the still image is fully decoded.

Preferred Video Compression Method — Motion Estimation The inventive compression technique for still images can be incoφorated into standard MPEG compression to enhance video compression through spatial hypercompression of each I frame inserted into the video bitstream. Alternatively, in the preferred embodiment, a novel motion estimation technique is employed which provides significantly greater compression due to temporal compression. According to the present invention, I frames are inserted according to video content. This is done by accumulating the error or difference between all corresponding microblocks or segments of the current frame and the predicted frame and comparing that accumulated error or difference to a threshold to determine whether the next subsequent frame sent should be an I frame. If the error or difference is large (i.e., when motion error is high), the I frame is sent. If the error or difference is small, the I frame is not sent and the frame sequence is unaltered. As a consequence, full synchronization of I frame insertion with changes in scene is achieved and bandwidth is significantly reduced because I frames are inserted only where necessary, i.e. , where content requires them. Thus, the present invention, for the first time, analyzes the errors between the I frame and the B and P frames into which it will be inserted to decide whether to insert the I frame at that point or not. Consequently, the present invention significantly increases the overall image compression ratio, while offering a simultaneous benefit of increased image quality. In addition, by using the technique of the present invention for video compression, the distances between I frames are enlarged, which leads to better motion estimation and prediction.

The video compression technique of the present invention may be used with both I frames compressed using the still ISMP compression encoding process of the present invention or standard I frame compression techniques. The most significant compression will occur if both the ISMP compression encoding process of the present invention and the motion estimation process of the present invention are used. It is worth noting that in existing systems, a reasonable quality video can be produced only if I frame compression is not higher than 20: 1 to 40: 1. With the present invention I frame compression of 300: 1 is achieved. The following table illustrates the improvement over standard compression of the inventive technique of fixed separation of I frames compressed with the inventive CT algorithm used in conjunction with the inventive variable separation of I frames compressed with the CT algorithm.

Motion estimation is important to compression because many frames in full motion video are temporally correlated, e.g., a moving object on a solid background such as an image of a moving car will have high similarity from frame to frame. Efficient compression can be achieved if each component or block of the current frame to be encoded is represented by its difference with the most similar component, called the predictor, in the previous frame and by a vector expressing the relative position of the two blocks from the current frame to the predicted frame. The original block can be reconstructed from the difference, the motion vector, and the previous frame. The frame to be compensated can be partitioned into microblocks which are processed individually. In a current frame, microblocks of pixels, for example 8 x 8, are selected and the search for the closest match in the previous frame is performed. As a criterion of the best match, the mean absolute error is the most often used because of the good trade off between complexity and efficiency. The search for a match in the previous frame is performed in a, for example, 16 x 16 pixels window for an 8 x 8 reference or microblock. A total of, for example, 81 candidate blocks may be compared for the closest match. Larger search windows are possible using larger blocks 8 x 32 or l6 x l6 where the search window is 15 pixels larger in each direction leading to 256 candidate blocks and as many motion vectors to be compared for the closest match. Once the third subsequent frame is predicted, the standard methods provide that the error between a microblock in the current frame and the corresponding microblock in the predicted frame are compared and the error or difference between them is determined. This is done on a microblock by microblock basis until all microblocks in the current frame are compared to all the microblocks in the predicted frame. In the standard process these differences are sent to the decoder real time to be used by the decoder to reconstruct the original block from the difference, the motion vector, and the previous frame. The error information is not used in any other way.

In contrast, in the present invention, the error or difference calculated between microblocks in the current frame and the predicted frame are accumulated or stored and each time an error is calculated between a microblock in the current frame and the corresponding microblock in the predicted frame that error is accumulated to the existing error for that frame. Once all the errors for all the blocks in the current frame as compared to the predicted frame are generated and summed, that accumulated error is then used to determine whether a new I frame should be inserted. This methodology is MPEG compatible and yields extremely high quality video images not possible with state of the art motion estimators. The accumulated error is used to advantage by comparing it to a threshold E₀ which is preset depending upon the content or type of the video such as action, documentary, or nature. If E₀ for a particular current frame is exceeded by the accumulated error, this means that there is a significant change in the scene which warrants sending an entire new I frame. Consequently, an entire new I frame is compressed and sent, and the motion estimation sequence begins again with that new I frame. If EQ is not exceeded by the accumulated error, then the differences between the current frame and the predicted frame are sent as usual and this process continues until E₀ is exceeded and the motion estimation sequence is begun again with the sending of a new I frame.

Now turning to Figure 16, the motion estimation process is now described in detail. At step 1, the next F₀ frame is taken. This frame may be the first frame of the video in which case it is an I frame or may be a subsequent frame. At step 2, if F₀ was compressed by standard DCT methods, the left branch of the flow chart of Figure 8 is followed. If F₀ was compressed using the inventive ISMP algorithm, the right branch of the flow chart in Figure 8 is followed.

First, assuming that F₀ was compressed using standard DCT methods, step 3 involves standard segmenting of the F₀ frame into search blocks having subblocks called microblocks and defining motion vectors which are used to predict the third subsequent frame after F₀. This is accomplished using standard techniques well known in the art.

At step 4, the error or difference between each microblock in F₀ and the corresponding microblock in the predicted third subsequent frame is defined for all microblocks in F₀. At this point the inventive motion estimation processed diverges from standard techniques.

If a standard MPEG encoder-decoder scheme was being used, these microblock differences would be sent from the encoder to the decoder and used by the decoder to reconstruct F₀. By sending only the differences between F₀ and the predicted third subsequent frame, significant compression is realized because it is no longer necessary to send an entire frame of information but only the differences between them. In accordance with standard MPEG encoder-decoder techniques, however, a new I frame is necessarily transmitted every 15 frames whether an I frame is needed or not. This poses two problems. Where the I frame is not needed, bandwidth is wasted because unnecessary bits are sent from the encoder to the decoder (or stored on disc if the process is not done real time). On the other hand, where the content of the video is such that significant scene changes occur from one frame to another much more often than every 15 frames, the insertion of an I frame every 15 frames will be insufficient to ensure a high quality video image at the decoder. For these reasons, the motion estimation technique of the present invention is especially valuable because it will, dependent upon the content of the video, insert or send an I frame to the decoder when the content of the video warrants it. In this way, a high quality image is maintained.

This is accomplished in the present invention by, as seen at step 5, accumulating the error between corresponding microblocks in the F₀ and the predicted third subsequent frame as each error is defined in step 4 for each microblock of F₀. The next step, step 6, is optional and involves normalizing the total accumulated error for F₀ by defining an average error A which is the total accumulated error divided by the number of microblocks in F₀. This yields a smaller dynamic range for the errors, i.e. , smaller numbers may represent the errors. Continuing with step 7, the accumulated error (whether normalized or not) is compared to a threshold error E₀. E₀ is chosen based upon video content such as whether the video is an action film, a documentary, a nature film, or other. Action videos tend to require insertion of I frames more often because there are more drastic changes in scene from one frame to another. It is especially important when compressing such videos to use the motion estimation techmque of the present invention which can insert additional I frames based on video content where necessary to keep video image quality high. In choosing EQ, bandwidth versus quality should be considered. If E₀ is set high, a high level of errors will be tolerated and fewer I frames will need to be inserted. Quality, however, will decrease because there will be an under utilization of bandwidth. If, on the other hand, E₀ is set too low, I frames will be inserted more frequently and available bandwidth may be exceeded and frames may start to drop out as commonly happens with MPEG. So the threshold E₀ should be tuned to video content. This can be done in real time by analyzing the video off-line and varying Eg in accordance with the statistics of the video, such as the number of cuts, the amount of action, etc. This process may be enhanced by using genetic algorithms and fuzzy logic. Where the accumulated error is greater than E₀, the next frame sent will be an I frame. In accordance with standard techniques, it is preferable that the I frame be compressed prior to sending it to the decoder. This reinitiates the sequence of frames at step 8.

If the accumulated error is less than E₀, the subsequent frame is not sent as an I frame but the differences are continued to be sent at step 9 to minimize bandwidth of the signal, sent between the encoder and decoder. The process then reinitiates at step 1 where the next frame Fj is taken. That next frame may not be an I frame but may instead be a subsequent frame, and the methodology is the same in either case. The next frame, whether it is an I, B, or P frame, is compared to the predicted third subsequent frame and the method continues as described above.

In an alternative embodiment, instead of sending the I frame as the next subsequent frame, the I frame could be sent as the current frame and used to replace error data for each microblock data for each microblock stored in the decoder buffers. This could be accomplished by clearing the buffers in the decoder holding errors between each of the microblocks F₀ and the predicted third subsequent frame and replacing that data with the I frame. Although not compatible with MPEG, it may be advantageous in certain situations to clear out the buffers containing the high error frame data and replace that data with the next frame as an I frame.

The motion estimation technique of the present invention may also be used to dynamically change or update compression ratio on a frame by frame basis by providing feedback from the receiver or decoder and using that feedback to change parameters of the compression engine in the registers of the video compression chips. For example, if the accumulated error calculated in the motion estimation technique of the present invention were too frequent or extraordinarily high, this information could be used to alter the parameters of the compression engine in the video compression chips to decrease the compression ratio and thereby increase bandwidth. Conversely, if the accumulated error over time was found to be unusually low, the compression ratio could be increased and thereby the bandwidth of the signal to be stored could be decreased. This is made possible by the accumulation of errors between the corresponding microblocks of the current frame (F₀) and the predicted third subsequent frame. This is not possible in prior art techniques because, although the error between corresponding microblocks of the current frame and the predicted third subsequent frame are calculated, there is no accumulated error calculated and no use of that accumulated error anywhere in the system. In the present invention, however, the accumulated error is calculated and may, in fact, be used on a frame by frame basis to decide whether the next frame should be an entire I frame as opposed to only the difference signal.

In a bandwidth on demand system, for example, if the feedback from the receiver indicates that there is a high bit error rate (BER), the transmitter may lower the bandwidth by increasing the compression ratio. This will necessarily result in a signal having sequences of different bit rates which are not possible in prior art MPEG systems. Intelligent systems such as genetic algorithms or neural networks and fuzzy logic may be used to determine the necessary change in compression ratio and bandwidth off-line by analyzing the video frame by frame. Turning now to the right branch of Figure 16, this branch is followed if the still compression method selected was the ISMP algorithm of the present invention which compresses each frame in accordance with catastrophic theory and represents the "structure" of that image in a highly compressible form using only the coefficients of canonical polynomials. Step 3 A in the right branch would be to predict the third subsequent frame from the current frame (here F₀) using standard techniques of defining the motion vectors of microblocks within the search blocks by template matching.

Step 4A would be to define the error between microblocks in F₀ and the microblocks in the predicted third subsequent frame. This is done using standard techniques. If a particular microblock in F₀ has a match with a microblock in the predicted frame, i.e., the error is 0, then the coefficients of the polynomial that were generated for that microblock when F₀ was compressed using the ISMP algorithm are then sent to the decoder and used along with the motion vectors generated in step 3A to reconstruct F₀. The sending of just the coefficients results in much higher than normal compression because the number of bits representing those coefficients is very small. If a microblock in F₀ has no match in the predicted third subsequent frame i.e., an error exists between those corresponding microblocks, new coefficients are generated for the corresponding microblock in the predicted third subsequent frame and those coefficients are sent to the decoder and used along with the motion vectors generated in step 3 A to reconstruct F₀. As an alternative, the newly generated coefficients for the corresponding microblock in the predicted third subsequent frame could be subtracted from the coefficients of the corresponding microblock in F₀ to even further compress the data. This may be done but is not necessary because the coefficients representing each microblock constitute highly compressed data already and further compression is not necessary. At step 5 A, the errors from the above comparison of F₀ and P are accumulated.

At step 6 A, the accumulated errors are normalized by the number of microblocks.

At step 7A, the accumulated error is compared to the threshold E₀. And if the accumulated error is greater than the threshold E₀ a new I frame is sent as the new subsequent frame at step 8A. If the accumulated error is less than the threshold error E₀ the coefficients that were newly generated for a particular microblock that did not find a match are continued to be sent to the decoder at step 9A. After both steps 8A and 9A the process reinitiates at step 1. Thus, according to the present invention, the error data is used and inteφreted in a novel way which provides high compression and quality imaging.

Motion Estimation Hardware

Referring now to Figure 17, the hardware for performing motion estimation is depicted in block diagram format. All of the hardware is standard. A host computer 10 communicates with the video processor board over PCI bus 14. The host computer 10 is preferably of at least the 100 MHz pentium class. PCI bus controller 12 controls communications over the PCI bus. EPROM 14 stores the coefficients and transfers them to the PCI bus controller 12 so that all the internal registers of the PCI bus controller 12 are set upon start-up. Input video processor 16 is a standard input video processor. It is responsible for scaling and dropping pixels from a frame. It has two inputs, a standard composite NTS signal and a high resolution Y/C signal having separated luminance and chromance signals to prevent contamination. The input video processor 16 scales the normal 702 x 480 resolution of the NTS input to standard MPEG-1 resolution of 352 x 240. The input video processor 16 also contains an A/D converter which converts the input signals from analog to a digital output.

Below input video processor 16 is audio input processor 18 which has as its input left and right stereo signals. The audio input processor 18 performs A/D conversion of the input signals. The output of the audio input processor 18 is input to a digital signal processor (DSP) audio compression chip 20 which is standard. The output of the audio compression chip 20 is input into the PCI bus controller 12 which can place the compressed audio onto the PCI bus 14 for commumcation to the host computer 10. Returning to the video side, the output of the input video processor 16 is input to an ASIC 22 (Application Specific Integrated Circuit) which is one chip of a three chip video compression processor also having a DTC based compression chip 24 and a motion estimator chip 26. The ASIC 22 handles signal transport, buffering and formatting of the video data from the input video processor 16 and also controls both the DTC based compression chip 24 and motion estimator chip 26. All of these chips are standard. An output of each of the chips 22, 24, and 26 of the video compression processor 23 is input to the PCI bus controller 12 for placing the compressed video on the PCI bus for communication to the host computer 10.

The compressed video stream from the video compression processor 23 on the board 13 undergoes lossless compression in the host computer using standard lossless compression techniques such as statistical encoding and run-length coding. After that lossless compression, the audio and video are multiplexed in standard fashion into a standard video signal. In order to have synchronization of audio and video the packets containing video and audio must be interleaved into a single bit stream with proper labeling so that upon playback they can be reassembled as is well known in the art. Importantly, the errors that were calculated in the motion estimator 26 between the current frame and the predicted third subsequent frame are transmitted to the host computer 10 over the PCI bus 14 so they can be transmitted to the encoder (not shown) to recreate the current frame at the encoder using that error or difference signal and the motion vectors generated during motion estimation. This is standard in the art. In accordance with the motion estimation of the present invention, however, that error is also accumulated in the host computer in a software routine in accordance with the motion estimation techniques of the present invention.

Referring now to Figure 18 is a flow chart describing error accumulation in the motion estimation procedure. At Step 1 the error buffer in the compression processor 23 is read through the PCI bus 14. At Step 2 that error is accumulated in an error buffer created in software in the host computer 10 so that the accumulated error will equal the preexisting error plus the present error. At Step 3 the accumulated error is compared to a threshold error and if the accumulated error is larger than the threshold error then a new I frame is sent and the error buffer in the compression processor need not be read again for that particular frame. If the accumulated error is not greater than the threshold error then the process loops back up to Step 4 where the next subsequent microblock in that frame is chosen. If there is a subsequent microblock in that frame then the process repeats at Step 1 where the error buffer in the compression processor is read. That error is accumulated in the accumulated error buffer at Step 2 and that accumulated error is compared to threshold at Step 3. Note that this looping will continue from Steps 1, 2, 3, and 4 until at Step 3 the accumulated error exceeds the threshold at which point it is not longer necessary to check any more microblocks for that frame because the error became so high that the host computer determined that a new I frame should be sent to restart the motion sequence. If, on the other hand, the accumulated error for all the microblocks of an entire frame never exceeds the threshold, then after Step 4, the process will go to Step 5 and the standard MPEG compression process will continue without changes, i.e., the next B or P frame will be grabbed and compressed.

Automatic Target Recognition (ATR) The ISMP still image compression methodology of the present invention can be used to greatly enhance automatic target recognition systems because the invention emphasizes and accurately represents the natural features such as "sculpture" of the object that makes human cognition of the target easier and more accurate. Furthermore, the polynomials used to represent the sculpture of the object are stable for small variations of projection direction or changes in movement, rotation, and scale of an object. This, too, enhances automatic target recognition.

Human vision defines objects principally by their contours. The human visual system registers a greater difference in brightness between adjacent visual images that are registered, faithfully recording the actual physical difference in light intensity. Researchers have shown that individual neurons in the visual cortex are highly sensitive to contours, responding best not to circular spots but rather to light or dark bars or edges. At it turns out, the fact that ISMP compression extracts exactly these edges and emphasizes the "sculpture" characteristic of the object makes it especially advantageous for use in ATR. By preserving object edges in the compressed information, the human visual system can extract important objects from a background even if the object has bulk colors very close to the colors of other objects in the background. This feature is extremely important for registration of military targets.

In virtually all ATR applications, the structures to be identified have sculpmre. Consequently, the sculpmre portion of the image can be extracted using the inventive methodology to achieve compression ratios of at least 4: 1. Unlike prior art methods based on linear methods and Fourier transforms like JPEG and wavelet, which destroy the very information which is essential for human cognition-soft edges, the present invention preserves those soft edges that exist in sculptures in virtually all structures to be identified. In contrast, the "texture" of an object is far less critical to human cognition. The present invention takes advantage of the distinction between sculpmre and soft edges and texture by separating the sculpmre characteristics of the object from the texture characteristics and utilizing only the sculpmre information for ATR. An additional benefit of this methodology is that the sculpmre information may be transmitted using relatively little bandwidth because it can be fully represented by polynomials whereas texture information requires greater bandwidth. A preferred method of ATR involves separating the texture and sculpmre portions of the image using the ISMP compression method, using standard soft ATR on the sculpmre portion, and then using standard hard ATR methods on the entire image (both texture and sculpture). Another preferred method for ATR in accordance with the presence invention is to split the texture and sculpmre portions of the image using a portion of the ISMP compression method, using state of the art soft ATR methods on the sculpmre part, and then using state of the art hard ATR methods on the texture part. This greatly reduces the number of bits that need to be transmitted because the texture information is dropped altogether. Quality, however, remains high because the sculpmre portion of the image was derived using ISMP which retains all necessary soft edge information which is critical to human cognition. Such soft edge information would be eliminated or lost, in any event, if standard Fourier transform type compression methods are used. There are numerous applications of ATR using datery obtained from the

ISMP method. The datery can be used for autonomous object target detection, tracking, zooming, image enhancement, and almost real-time early stage recognition puφoses. The present invention provides the capability for smart network-based cooperative scene processing such as in remote intelligent consolidated operators ("RICO") where information from remote camera networks must be transmitted over a smart local area network (LAN) which interconnects a number of camera platforms for cooperative wide area surveillance and monitoring. For example, a camera platform (with the inventive ISMP method embedded therein) can extract features of the objects seen such as critical soft edge information. It can transmit those images over a smart LAN to adjacent camera platforms. This process may provide cooperative scene information transmission outside the coverage of the original or any single camera platform. Through this process, observers of a scene can perceive the "big picture" .

The images must then be transmitted from the remote camera network to a central station which may provide editing of film by computer to create the big picture. The invention will benefit such a system in two ways. First, because the ISMP compression method emphasizes the sculpmre characteristics of the object, it enhances the ability to recognize the object imaged. Second, because the "sculpmre" characteristics of the object are emphasized and represented using discrete numbers or coefficients from polynomials, the data sent is highly compressed which increases bandwidth significantly. Another application is as an autonomous movie director where standard

ATR is used and that information is compressed using the present invention for sending those images from the cameras to the central station. Because of the large volume of information that can be generated in such a system, the images must be compressed sufficiently so that they do not overwhelm the host computer. This is a real problem that is solved by the hypercompression of the present invention. These benefits apply to a wide range of systems including battlefield imaging systems and anti-terrorist recognition applications as well as full mapping capabilities.

Another significant advantage of the present invention is the ability to provide sufficiently high compression ratios for providing TV-class transmission through traditional air communication channels which are 64 kbps or less. In fact, the invention can provide such a significant compression ratio improvement of more than an order of magnitude that, generally speaking, "video through audio" is made possible. In other words, the present invention makes it possible for battlefield commanders and others to receive image information as opposed to raw data. And because the image information they receive is sent in the form of discrete numbers or coefficients of polynomials that relate to isomoφhic singular manifolds in the object, the data are highly compressed. And although highly compressed, the data preserve full information about the objects 3-D boundaries or soft edges.

An example of a real-time remote engagement (RTRE) air scenario made possible because of the present invention includes providing an aviator who is approaching his target (at sea, on the ground, or in the air) with a short TV-relay from an overflying military communications aircraft or satellites that upgrades the present target location at the last minute. This can prevent an aviator from losing track of a highly mobile target. This is made possible because the data are highly compressed and can be sent over low bandwidth air channels of 64 kbps or less and because the information that is sent preserves edge information which makes it possible for the pilot to easily recognize his target.

Because the typical air communication channels are of low bandwidth, the ability to use all that bandwidth is critical. The present invention's ability to dynamically allocate bandwidth on demand permits the use of small fractions of standard 64 kbps bandwidths for bursty compressed video/graphic image transmission. A typical air communications channel must accommodate signals of different types such as imagery, audio, sensory data, computer data and synch signals etc. The higher level protocols of the network will prioritize these different signals. Conservatively speaking, imagery is one of the lowest priority because in most cases operations can continue without it. Therefore, imagery information typically is relegated to using only the bandwidth that is available and that available bandwidth changes with time. It is extremely useful to use the ISMP method of the present invention which can be implemented with a tunable compression ratio. This is distinct from software which changes compression ratio based on the type of the object. Furthermore, intelligence systems such as genetic algorithms or fuzzy logic and neural networks can provide intelligent control of the available bandwidth and permit imagery data to be sent where otherwise it was not possible to do so. The severe constraints placed on the trade-off between the compression ratio and the PSNR by standard air channels of 64 kbps or less are highlighted by the following example. To compress data into the required data rate of 64 kbps from a fully developed synthetic aperture radar (SAR), for instance, uncompressed bandwidth of 13 Mbps, a 203:1 still image compression rate is needed. (512)² number of pixels, 10-bit grey level, and 5 Hz bursty frame rate, yields (512₂ x 10 x 5 = 13 Mbps). The situation is made even more severe for VGA full-motion video (221 Mbps) which requires a 3452: 1 motion video compression rate. The ability of the motion estimator of the present invention which inserts I frames only where the content of the video requires it can provide ten times better compression ratios than prior art systems, namely, up to 4000: 1. Thus, signals from an SAR uncompressed bandwidth of 13 Mbps may in fact, for the first time be sent through 64 kbps channels. This is made possible by the non- intuitive use of Arnold's Theorem according to which local isomoφhism (i.e., 1 : 1 direct and inverse relation) exists between the 3-D object boundary and its 2-D image. As a result, the most critical part of the object— its 3-D boundary may be described only by a 1-D contour and 3 or 4 natural digits that characterize a simple catastrophic polynomial. This creates tremendous lossless compression of object boundaries and still preserves high image quality. Experimental results show that the ISMP methodology of the present invention in contrast to state of the art compression methods can achieve a compression ratio of 160: 1 at PSNR = 38 Db with almost invisible artifacts whereas the prior art offers only CR = 60: 1 at a lower PSNR value of 26 Db. The difference in the image is significant.

Where far less than 64 kbps bandwidth channels are available, the hypercompression made possible by the present invention can permit continuity of images by "cartooning" which allows transmission of reduced real-time videoconferencing even on an 8 kbps communications channel.

Referring now to Figure 19, Figure 19 shows five categories of data reduction as a fraction of the original, five reduced data rates, and five different outcomes of those data rates. Category A represents 100% of the original which is a data rate of 64 kbps. At this data rate, the original video may be sent.

Category B represents a data reduction as a fraction of the original 75 % which is a 48 kbps data rate. The result of this transmission is that tiny details of the face or other structure are still recognizable and edges remain unchanged. Category C represents a data reduction of 50% as a fraction of the original which is a data rate of 32 kbps. The result of this transmission is that edges are hardened and there are smooth transitions for face details. Category D represents a 25% data reduction as a fraction of the original which is a reduced data rate of 16 kpbs. The result of this transmission is a heavily reduced texture and hard edges but it is still possible to recognize a human face. Category E represents a 10% reduction in data as a fraction of the original which is a reduced data rate of 12.8 kpbs. The result of this transmission is hard edges and "cartoon" type faces. While cartooning certainly does not provide optimum viewability, it may be more than adequate, for example, for soft ATR puφoses where a tank need only be distinguished from a plane or other categories of objects and the type of model of each is not required to be determined. Additionally, it was not possible to send even cartoon type images over low bandwidth communication channels using prior art methods and therefore the ability to send a cartoon type image over that communication channel where no image was possible before is a great advance. Thus, depending upon the quality of transmission required and the application, the compression techniques of the present invention may be utilized to achieve a broad array of results heretofore unobtainable using prior art compression methods. Various modes of carrying out the invention are contemplated as being within the scope of the following claims particularly pointing out and distinctly claiming the subject matter which is regarded as the invention.

Claims

CLAIMS We claim:

1. A method of compressing an image said method comprising the steps of: segmenting said image; creating a modeled surface for each segment; connecting the segments to create an entire modeled image; generating a texture image; and combining said texture image and entire modeled image.

2. The method of claim 1 wherein the step of creating a modeled surface for each segment further comprises the step of choosing a canonical polynomial to represent isomoφhic singularities in said image.

3. The method of claim 1 further comprising the steps of applying lossy compression to said texture image; and applying lossless compression to the combined said texture image and said entire modeled image.

4. A method of compressing an image, I₀ said method comprising the steps of: dividing said image I₀ into segments, each segment having a plurality of pixels; calculating the dynamic range of the pixels in each segment; selecting a best match canonical polynomial for each of said segments; finding substitutes for variables in said canonical polynomial to calculate a modeled surface equation F for each of said segments; creating a modeled surface I_M for each of said segments by substituting the coordinates of each pixel into the modeled surface equation F; storing the coefficients for the modeled surface equation F for each of said segments; finding connections between adjacent ones of said segments; creating an entire modeled image, I_M, from each of said segments; finding the difference between the image I₀ and the entire modeled image I_M to create a texture image I_d; applying standard lossy compression to said texture image I_d; storing the texture image I_d; combining the entire modeled image I_M and the texture image I_d to create a combination image, and applying lossless compression to said combination image.

5. A method of compressing a still image said method comprising the steps of: identifying at least one catastrophe in said image; representing said catastrophe with a canonical polynomial; transforming said canonical polynomial into datery.

6. A method of compressing a still image, said method comprising the step of: representing a region of abrupt changes in pixel intensity in said still image with a canonical polynomial.

7. A method of compressing an image, said method comprising the steps of: segmenting the image into blocks of pixels; creating a canonical polynomial surface for at least one catastrophe in at least one of said blocks of pixels; sending the coefficients of said at least one canonical polynomial as compressed data.

8. A method of compressing an image, said method comprising the steps of: identifying at least one isomoφhic singularity in said image by applying photometric projection to said image; characterizing said at least one isomoφhic singularity with at least one polynomial.

9. The method as defined in claim 8, further comprising the step of: creating a modeled surface of said image with said at least one polynomial, said modeled surface being isomoφhically related to said image.

10. The method of claim 8 further comprising the step of: transmitting the coefficients of said at least one polynomial as compressed data.

11. A method of compressing an image having manifolds, said method comprising the steps of: modeling the image as a photometric projection of at least one manifold in said image; mapping said at least one manifold in coordinates (x,y,B) where the coordinate B is luminance at each point (x,y); characterizing the mapping with a polynomial, said polynomial having coefficients; sending the coefficients of the polynomial as compressed data.

12. The method of claim 11, wherein the coordinate B does not introduce new singularities.

13. A method of compressing an image, said method comprising the step of describing the shape of object boundaries in the image in polynomial form.

14. A method of compressing an image, said method comprising the steps of: segmenting the image into segments; creating a modeled surface for each segment, said modeled surface for each segment being isomoφhic with respect to each segment; connecting adjacent segments to create an entire modeled image, said entire modeled image being isomoφhic with respect to said image.

15. The method as defined in claim 14, further comprising the steps of: calculating the peak signal to noise ratio over the entire modeled image; calculating the difference between said image and said entire modeled image to retrieve texture information of said image.

16. A method of compressing video, said method comprising the steps of: determining the error between a current frame and a predicted frame; inserting an I frame as the next subsequent frame after the current frame if said error exceeds a predetermined threshold.

17. A method of compressing video, said method comprising the steps of: taking a frame F₀ of the video; segmenting the frame F₀ into search blocks; predicting a subsequent frame; determining the error between the frame F₀ and said predicted frame; comparing said error to a threshold; inserting an I frame as the next subsequent frame if said error exceeds said threshold.

18. The invention as defined in claim 17, further comprising the steps of segmenting the frame F₀ into microblocks within said search blocks and normalizing the error by dividing by the number of microblocks.

19. The method as defined in claim 17, wherein said threshold is based upon video content.

20. The method as defined in claim 17, wherein said I frame is inserted in place of said frame F₀ if said error exceeds said threshold.

21. The method as defined in claim 17, further comprising the step of dynamically changing the compression ratio on a frame by frame basis based upon said error.

22. The invention as defined in claim 21, wherein one of a genetic algorithm, neural network, and fuzzy logic are used to determine the necessary change in compression ratio.

23. A method of compressed video transmission, the method the steps of: taking a first frame F₀; segmenting said frame F₀ into blocks; compressing each said block by representing it with a canonical polynomial having original coefficients; predicting a frame P subsequent to said frame F₀; determining the error between each block in said frame F₀ and said predicted frame P; accumulating the error determined between each block in said frame

F₀ and said predicted frame P; comparing said accumulated error to a threshold; if said accumulated error exceeds said threshold, inserting an I frame as the next subsequent frame to said frame F₀.

24. A method as defined in claim 23 wherein said I frame was previously compressed by representing it with canonical polynomials.

25. A method of compressed video transmission, the method comprising the steps of: taking a first frame F₀; segmenting said frame F₀ into blocks; compressing each said block by representing it with a canonical polynomial having original coefficients; predicting a frame P subsequent to said frame F₀; comparing each block in said frame F₀ with said predicted frame P to determine if there is a match; if a match is found in said frame P for a block in said frame F₀, sending said coefficients of said polynomial for that block to a decoder; if a match is not found for a block in said frame F₀, generating new coefficients of said polynomial representing said predicted frame P and sending said new coefficients to said decoder; reconstructing said frame F₀ in said decoder from said original coefficients and said new coefficients sent to said decoder; calculating an error between said frame F₀ and said predicted frame P based on said comparing each block in said frame F₀ with said predicted frame P; comparing said error to a threshold; if said error does not exceed said threshold, sending a B or P frame as the next subsequent frame to said frame F₀; if said error exceeds said threshold, inserting an I frame as the next subsequent frame to said frame F₀.

26. The method as defined in claim 25, further comprising the steps of replacing said original coefficients, and said new coefficients with an I frame.

27. A method of compressing a video image, said method comprising the steps of: taking a frame F₀; segmenting F₀ into blocks and defining motion vectors for said blocks to predict a subsequent frame P having corresponding blocks; defining errors between said blocks in said frame F₀ and said corresponding blocks in said frame P; accumulating said errors; and based on said accumulated error, sending an I frame as the next subsequent frame to said frame F₀.

28. A method of automatically recognizing a target in an image, said method comprising the steps of: segmenting said image into segments; creating a modeled surface for each segment; connecting the segments to create an entire modeled image; and comparing said entire modeled image to a library of known images to determine if there is a match with a known image in said library of images.

29. The method of claim 28 wherein the step of connecting said segments to create said entire modeled image produces a target image having the sculpmre characteristics of the image.

30. A method of automatic target recognition, said method comprising the steps of: taking an image having texture and sculpmre characteristics; processing the image so as to separate said texture characteristics from said sculpmre characteristics; comparing said sculpture characteristics to a library of known images.

31. A method of analyzing an image for automatic target recognition, said method comprising the steps of: separating the image into texture components and sculpmre components; applying soft ATR to said sculpmre components to create a soft ATR sculpmre component; combining said texture components and said soft ATR sculpmre components to form a combined image; applying hard ATR to said combined image.

32. The method of claim 1 whereby connections are found between adjacent segments and a model surface for a segment is extended to any such adjacent segment that may be satisfactorily approximated by such model surface.

33. The method of claim 32 whereby said connections are based on average distance.

34. The method of claim 1 wherein a texture image is generated only if the model image signal to noise ratio exceeds a predetermined threshold.

35. The method of claim 4 wherein said texture image is generated only if the model image signal to noise ratio exceeds a predetermined threshold.

36. The method of claim 6 wherein a model image is generated from said canonical polynomial which is a best match canonical polynomial selected based on a difference between the original image and said model image.

37. The method of claim 7 wherein a model image is generated from said canonical polynomial which is a best match canonical polynomial selected based on a difference between the original still image and said model image.

38. The method of claim 10 further comprising the step of selecting said coefficients based upon a quality determination.

39. The method of claim 10 wherein a model image is generated from said canonical polynomial, which is a best match canonical polynomial selected based on a difference between the original image and said model image.

40. The method of claim 39 wherein said difference is calculated using the equation:

Q =

41. The method of claim 13 wherein said polynomial form comprises a canonical polynomial selected based on predetermined criteria indicative of image quality.

42. The method of claim 15 further comprising compressing said texture information of said image by a lossy compression method if the difference exceeds a predetermined threshold.

43. The method of claim 13 wherein said object boundaries comprise at least one isomoφhic singularity.

44. The method of claim 43 wherein the polynomial form is selected from a table of canonical polynomials.

45. The method of claim 44 wherein said table of canonical polynomials comprises

46. The method of claim 45 wherein a transformation is applied to said selected canonical polynomial to obtain a function describing a model surface.

47. The method of claim 46 wherein said transformation is a non-homogeneous linear transformation.

48. The method of claim 47 wherein said non-homogeneous linear transformation takes the form:

f ³

J canonical ~ ^Xi ^{+ X}\^X2

49. The method of claim 48 wherein a model surface is obtained from said non-homogeneous linear transformation as follows:

= (?! + a i ^aJn)^'> ^x2 = Cy₂ ^{+ b}╧çy╬╣ bfn⁾-

50. The method of claim 49 wherein a quality determination is calculated by determining the difference between the original and model segments on a pixel-by- pixel basis using the equation:

Λf-l N-\

^Q A MN Vote?⁾ - ⁱ to^yy?

51. A method of imaging compression comprising the step of characterizing aspects of an image to be compressed with singular manifold representations.

52. The method of claim 51 wherein said aspects are surfaces of objects.

53. The method of claim 51 wherein said singular manifold representations are represented by canonical polynomials.

54. The method of claim 53 further comprising the step of reducing said polynomials to compact tabulated normal form polynomials which comprise simple numbers.

55. A method of still image encoding comprising the following steps: (a) capturing a frame;

(b) dividing said frame into segments of pixels;

(c) determining the dynamic range of a segment by subtracting the intensity of the pixel having maximum intensity from the intensity of the pixel having minimum intensity in said segment; (d) comparing said dynamic range to a threshold below which said segment is likely to represent background;

(e) selecting a canonical polynomial from a table if said threshold in (d) above is exceeded;

(f) compressing said segment using standard texture compression techniques and storing the result if the threshold of step (d) is not exceeded;

(g) performing a transformation on said canonical polynomial to obtain an equation representing a modeled surface;

(h) substituting the coordinates of each pixel from said segment into said equation representing said modeled surface to obtain a matrix of modeled surface elements of said segment;

(i) calculating the overall quality of the modeled surface of said segment compared with said original segment by (1) subtracting the difference between the pixels of said original segment and corresponding pixels of said modeled surface (2) squaring said differences (3) summing up all of said squares and (4) taking the square root of said sum to arrive at a quality determination for said modelled surface;

(j) comparing said quality determination of step (i) to a predetermined threshold;

(k) selecting new coefficients for said canonical polynomial if said quality determination is greater than said predetermined threshold of step (j) and repeating steps (i) and (j) until a best quality determination, less than said predetermined threshold of step (i) is achieved; (1) storing said best quality determination for said canonical polynomial and said coefficients that produced said best quality determination;

(m) repeating steps (f-1) for polynomials not yet tested until all canonical polynomials from said table have been tested for said segment;

(n) determining the polynomial having the overall best quality determination of the polynomials tested for said segment to arrive at a selected polynomial for said segment;

(o) storing the coefficients for said selected polynomial representing a model surface for said segment;

(p) selecting a next segment of said frame and performing steps (c) through (o) on all such next segments until all segments of the frame have been selected;

(q) calculating the average distance between said model surface of said segment and each adjacent segment of said frame to determine if connections to neighboring segments can be made; (r) comparing the average distances determined in the preceding step to a threshold average distance;

(s) extending said model surface to adjacent segments if the average distance between such segments is less than said threshold average distance;

(t) calculating a spline to approximate the surface of adjacent segments if the average distance for any such segment exceeds said threshold average distance to form a graph;

(u) constructing a model image of the entire frame by creating a table of all of the data representing the modelled segments to obtain a matrix describing the entire modeled frame surface;

(v) calculating the peak signal to noise ratio over the entire frame;

(w) comparing the peak signal to noise ratio of the entire frame to a signal to noise threshold;

(x) calculating a difference frame by subtracting the value of each pixel of the model image from each pixel of the original capmred frame if the peak signal to noise ratio exceeds said signal to noise threshold;

(y) compressing the difference frame, if any, using standard lossy compression methods;

(z) compressing the frame data comprising the coefficients for said selected polynomials, and said compressed difference frame, if any.

56. The method of claim 55 further comprising applying a nonhomogeneous linear transformation to obtain said matrix of modelled surface elements of said segment.

57. The method of claim 55 further comprising applying a homogeneous linear transformation to obtain said matrix of modelled surface elements of said segment.

58. The method of claim 55 further comprising applying a nonhomogeneous nonlinear transformation to obtain said matrix of modelled surface elements of said segment.

59. The method of claim 55 wherein said step of segmenting said frame into blocks of pixels comprises segmenting into blocks that are fixed in size.

60. The method of claim 55 wherein said step of segmenting said frame into blocks of pixels comprises segmenting into blocks that are variable in size.

61. A method of motion estimation for enhancing video compression comprising the steps of: accumulating the error between segments of a current frame and segments of a predicted frame; comparing the accumulated error to a threshold to determine whether the next frame sent should be an I frame.

62. A method of video compression utilizing motion estimation comprising the steps of:

(a) comparing a segment of the current frame to be encoded to a corresponding segment of the previous frame to arrive at an error difference; (b) repeating the calculation of Step (a) until error differences have been determined for all segments of the frame;

(c) summing said error differences to arrive at a cumulative error difference for that frame;

(d) comparing said cumulative error difference to a threshold to determine whether a new I frame should be inserted in the compressed video data stream.

63. The method of claim 62 further comprising adjusting the compression ratio on a frame by frame basis based on feedback and further using such feedback to change parameters of the compression engine.

64. The method of claim 62 further comprising lowering bandwidth by increasing the compression ratio based on feedback indicating that there is a high bit error rate.

65. The method of claim 63 wherein the change in compression ratio and bandwidth are determined by one of, a genetic algorithm, a neural network and fuzzy logic.

66. The method of claim 62 further comprising the step of presetting said threshold based upon the program content or type of video images to compressed such as action, documentary, or nature.

67. The method of claim 62 further comprising applying spatial hypercompression utilizing isomoφhic singular manifold representations to each I frame inserted into the compressed video data stream.

68. A method of image compression comprising the following steps:

(a) subdividing an original image I₀ into blocks of pixels;

(b) creating a canonical image of each block by finding a best match between one of a predetermined set of canonical polynomials and the intensity distribution for each block of pixels by using standard merit functions;

(c) creating a model image I_ra comprised of said canonical images and by finding connections between neighboring blocks of pixels thereby smoothing out intensity physical structure of said modelled image; (d) recapturing lost high frequency content of said image, if any, by subtracting said model Image I_m from said original image I₀.

69. A method of automatic target recognition comprising the following:

(a) separating the texture and sculpmre portions of an image using isomoφhic singular manifold projection; (b) utilizing standard hard automatic target recognition methods on the entire image including both of said texture and sculpture portions.

70. A method of automatic target recognition comprising the following steps:

(a) providing the texture and sculpture portions of an image using a portion of the isomoφhic singular manifold projection compression method; (b) using a soft automatic target recognition method on said sculpture portion of said image;

(c) using a hard automatic target recognition method on said texture portions of said image.

71. A method of automatic target recognition comprising using datery obtained from an isomoφhic singular manifold projection.

72. A method of object target detection comprising using datery obtained from an isomoφhic singular manifold projection.

73. A method of tracking a target comprising using datery obtained from an isomoφhic singular manifold projection.

74. A method of zooming comprising using datery obtained from an isomoφhic singular manifold projection.

75. A method of image enhancement comprising using datery obtained from an isomoφhic singular manifold projection method.

76. A method of automatic target recognition utilizing isomoφhic singular mamfold projection whereby critical soft edge information may be extracted and transmitted over a smart local area network between adjacent camera platforms to provide cooperative scene information for an observer.

77. A method of still image decoding using isomoφhic singular mamfold projection in which surfaces of objects have been represented by canonical polynomials comprising the following steps:

(a) decoding lossless compression data;

(b) separating the data representing model image frame I_m from any other data in the bit stream;

(c) extracting a segment of data;

(d) testing said data to determine whether it belongs to a graph;

(e) if said data does belong to a graph, constructing a modeled segment i_m for the segment represented by said data using the canonical polynomial that was stored for said graph;

(f) connecting separate graphs utilizing spline functions;

(g) constructing a frame I_m from such graphs and segments; (h) testing to determine if an I_d' frame representing high frequency components of the image to be decoded is present; (i) decompressing frame I_d' if present;

(j) deriving output frame I₀' from the combination of frame I_m and any such frame I_d' found to be present.

78. A method of video compression, said method comprising the steps of deriving an error between microblocks in a current frame of video and a predicted frame; accumulating said errors; comparing said accumulated error to a threshold; sending an I frame when said threshold is exceeded.

79. The method of claim 78 further comprising sending the I frame as the current frame so as to replace error data.

80. The method of claim 78 wherein the threshold is varied in accordance with statistical information of the video based on program content such as the number of cuts, the amount of action, etc.

81. The method of claim 78 further comprising varying said threshold by using genetic algorithms.

82. The method of claim 78 further comprising varying said threshold by using fuzzy logic.

83. The method of claim 78 further comprising applying spatial hypercompression utilizing isomoφhic singular manifold representations to each I frame inserted into the video bitstream.

84. An apparams for performing motion estimation comprising: a host computer; a video processor board linked to said host computer by a PCI bus for providing a stream of compressed video data to said PCI bus; an input video processor for converting a standard NTS signal and a high resolution Y/C signal to digital video data; an audio input processor for converting analog audio signals to digital audio data; a digital signal processor for processing digital audio data output from said audio input processor; an application specific integrated circuit for receiving output from said input video processor and for buffering and formatting said data; a DTC based compression chip for receiving video data from the application specific integrated circuit; a motion estimator chip for forming motion estimation calculations on data received from said application specific integrated circuit; a buffer for accumulating errors between a current frame and a predicted frame.

85. The apparatus of claim 84 wherein said motion estimation is performed in a video teleconferencing system.

86. The apparams of claim 84 wherein said motion estimation is performed as part of an MPEG editor.

87. The apparams of claim 84 further comprising a comparator for comparing said accumulated error to a predetermined threshold.

88. The apparatus of claim 87 further comprising means for sending an I frame whenever said accumulated error exceeds said predetermined threshold.

89. The apparams of claim 87 further comprising means for adjusting said predetermined threshold based on the program content of said video data.

90. The apparams of claim 87 further comprising means for adjusting said predetermined threshold based on fuzzy logic.

91. The apparams of claim 87 further comprising means for adjusting said predetermined threshold based on a neural network.

92. The apparatus of claim 87 further comprising means for adjusting said predetermined threshold based on a genetic algorithm.

93. The apparams of claim 87 further comprising program means for applying spatial hypercompression utilizing isomoφhic singular manifold representations to each I frame inserted into the video bitstream.

94. A method of video compression, said method comprising the step of applying temporal catastrophic formalisms to MPEG hypercompression by replacing the average error parameters by integrated luminance-changes.