A moving video image estimation system

Abstract: A moving video image estimation system is provided, in which, when a video image occurring at an arbitrary time t between time n and time n + 1 or at an arbitrary time t after time n + 1 is produced for time series video images given at discrete times 1, 2, ... n, n + 1, ..., even if the position and shape of an object in the video image are considerably changed between the frames, the video image of the object in the original frame can be prevented from being reproduced doubly in the interpolated video image. Based on an original video image of time n and time n + 1 , the centroid, the principal axis of inertia, the moment about the principal axis of inertia and the moment about the axis perpendicular to the principal axis of inertia are obtained. By using them, an affine transformation for transforming the original video image at time n to the original video image at time n + 1. Further, based on the infinitesimal transformation (A), {exp (At), and exp(A (t - 1))} obtained by making the affine transformation continuous with regard to time is executed on the original video image at time n and time n + 1. The results are synthesized to perform an interpolation between the frames. {exp (A (t - i))} is applied to the original video image at time n + 1, the video image after time n + 1 is predicted.

Description: Field of the Invention The present invention generally relates to interframe estimation systems for interpolating moving pictures between frames and, more particularly, to a system for generating a video image occurring at an arbitrary time t between video images provided at times n and n + 1 and at an arbitrary time t after time n + 1 for a time series of video images provided at discrete times 1, 2, ... n, n + 1, ... Description of the Related Art In currently available television systems for mass media broadcasts, video images are output every one-thirtieth of a second to accurately express a moving picture. If, for example, only ten video images are output per second, a moving picture becomes hard to watch. In the field of so-called animation, a number of slightly different video images are formed and output at the rate of 30 images per second, thereby obtaining a smooth moving picture. This causes the problem that a vast number of processes are required to form these images manually. However, if some intermediate images are omitted and these omitted images are then interpolated to form a moving picture, the number of processes can be considerably reduced. Further, in the medical field, a user of an ultrasonic diagnostic machine may want to obtain continuous ultrasonic video images of, for example, a cardiac blood stream flow. However, in the present technology, there is the problem that only 30 continuous video images can be obtained per second. Under such circumstances, an interframe interpolation system for effectively obtaining a video image at an arbitrary time t between times n and n + 1 and at an arbitrary time t after time n + 1 on the basis of video images given at discrete times n and n + 1 is needed. Figures 1A and 1B are schematic diagrams used to explain a conventional moving video image estimation system. Figure 1A shows examples of time series video images at time n and time n + 1, and Figure 1B shows problems encountered with existing interframe interpolation systems. When a video image at an arbitrary time t between time n and time n + 1 is produced for time series video images given at discrete times 1, 2, ... n, n + 1, ... {see Figure 1A}. According to the prior-art real moving picture interpolation system, concentrations at every pixel of the video image at time n and the video image at time n + 1 are interpolated linearly. More specifically, when the video image at time n is expressed as I(x,y;n) and the video image at time n + 1 is expressed as I(x,y;n+1), a video image J(x,y;t) to be produced at an intermediate time t between times n and n + 1 is calculated from the following equation. J(x,y;t) = (n+1-t)I(x,y;n) + (t-n)I(x,y;n+1) In accordance with this conventional system, when the position and shape of an object in a video image are considerably changed between the frames, there is then the disadvantage that the video image of the object in the original frame will be doubly reproduced in the interpolated video image. Summary of the Invention Accordingly, it is an object of the present invention to provide an improved moving video image estimation system which can eliminate the aforenoted shortcomings and disadvantages encountered with the conventional moving video image estimation system. More specifically, it is an object of the present invention to provide a moving video image estimation system which performs an estimation of a real moving picture. It is another object of the present invention to provide a moving video image estimation system in which, when a video image of an arbitrary time t between time n and time n + 1 is produced for time series video images given at discrete times 1, 2, ... n, n + 1, ..., even if the position and shape of an object in the video image are considerably changed between the frames, the video image of the object in the original frame is not doubly reproduced in the interpolated video image. According to a first aspect of the present invention, a moving video image estimation system is comprised of moment calculating means for calculating primary and secondary moments of a video image; affine transformation calculating means utilizing an output of the moment calculating means to determine affine transformation (G) by which a first-frame video image (I 1 ) is geometrically transformed into a second-frame video image (I 2 ); first affine transformation continuing means for multiplying by an arbitrary time constant (t) an infinitesimal transformation (A) given by the relation G = exp (A) with respect to said affine transformation (G), and obtaining an exponential transformation {exp(At)} as a continued affine transformation; and affine transformation executing means for performing said continued affine transformation {exp(At)} to thereby obtain a video image at a time after an arbitrary time (t) from a time corresponding to the first-frame video image. Further, said moving video image estimation system may be comprised of second affine transformation continuing means for multiplying said infinitesimal transformation (A) by another arbitrary time constant (t - 1) and obtaining an exponential transformation {exp(A (t - 1))}; whereby said affine transformation executing means being further adapted for obtaining a video image at the arbitrary time (t) between the first and second frames as a linear sum of a result which is provided by executing a transformation {exp(At)}, provided as an output of the first affine transformation continuing means, onto the first frame video image (I 1 ) and of a result provided by executing a transformation {exp(A (t - 1))}, provided as an output of the second affine transformation continuing means, onto the second frame video image (I 2 ). Further preferred embodiments of the present invention are defined in the subclaims. The above, and other objects, features and advantages of the present invention will become apparent in the following detailed description of illustrative embodiments to be read in conjunction with the accompanying drawings, in which like reference numerals are used to identify the same or similar parts in the several views. Brief Description of the Drawings Figures 1A and 1B are schematic diagrams used to explain a conventional interframe interpolation system, Figures 2A and 2B are function block diagrams, used to explain the principle of the present invention, Figure 3A is a schematic diagram used to explain the affine transformation of the present invention and Figure 3B is a flowchart used to explain how to estimate a video image, Figure 4 is a schematic block diagram showing a first embodiment of a system which utilizes the moving video image estimation system of the present invention, Figure 5 is a schematic block diagram showing an arrangement of a hardware of the first embodiment shown in Figure 4, Figure 6 is a flowchart used to explain a video image moment characteristic extracting process, Figure 7 is a schematic representation used to explain how to calculate a projection onto an oblique axis, Figure 8 is a flowchart used to explain a processing for determining affine transformation, Figure 9 is a flowchart used to explain a first affine transformation continuing process, Figure 10 is a flowchart used to explain a second affine transformation continuing process, Figure 11 is a flowchart used to explain how to assume a video image at an arbitrary time, Figures 12A to 12F are views showing examples of video images provided by the moving video image estimation system of the present invention, Figure 13 shows a block diagram of an arrangement of the second embodiment of the moving video image data estimation system, and Figure 14 shows a block diagram of an arrangement of the third embodiment of the moving video image data estimation system. Detailed Description of the Preferred Embodiments The principle of the present invention is described below with reference initially to Figures 2A and 2B. Figure 2A shows in a block diagram form the principle of the first principle of the present invention. As shown in Figure 2A, a moment calculating means 1 calculates primary and secondary moments of a video image and utilizes these primary and secondary moments to calculate, for example, coordinates of centroid of video image and central moments. An affine transformation calculating means 2 utilizes the coordinates of centroid of video images and central moments, calculated by the moment calculating means 1 from first-frame video image I 1 and second-frame video image I 2 , to determine an affine transformation (G) which geometrically transforms the first frame video image into the second frame video image. Coefficients calculated for the affine transformation herein are, for example, direction  of the principal axis of inertia of each video image, moment λ about the principal axis of inertia, moment κ about the axis perpendicular to the principal axis of inertia and matrix elements and vector elements for transforming coordinates indicating the first frame picture into coordinates indicating the second frame picture. Then, an affine transformation continuing means 3 multiplies by a constant t infinitesimal transformation A given by the relation G = expA, as the affine transformation G transforms the first frame video image calculated by the affine transformation calculating means 2 into the second frame video image, and then obtains exponential transformation exp(At) as a continued affine transformation. An affine transformation executing means 4 executes the continued affine transformation exp(At) onto the first frame video image I 1 to thereby obtain a video image at a time after an arbitrary time (t seconds) from a time corresponding to the video image of the first frame. Figure 2B shows in a function block form the second principle of the present invention. Referring to Figure 2B, the actions of the moment calculating means 1 and the affine transformation calculating means 2 are exactly the same as those of the first embodiment shown in Figure 2A. As shown in Figure 2B, in exactly the same manner as the affine transformation continuing means 3 in the first principle, a first affine transformation continuing means 5 multiplies by the constant t the infinitesimal transformation A, given by the relation G = expA. The affine transformation G geometrically transforms the first video image I 1 calculated by the affine transformation calculating means 2 into the second frame video image I 2 , and obtains the exponential transformation exp(At) as the first continued affine transformation. A second affine transformation continuing means 6 multiplies the infinitesimal transformation A by a different constant (t - 1) and obtains exponential transformation exp(A (t - 1)) as a second continued affine transformation. An affine transformation executing means 7 obtains a video image at an arbitrary time t after the first frame and before the second frame as a linear sum of a result provided by executing the transformation exp(At) provided as the output of the first affine transformation continuing means 5 to the first frame video image I 1 and a result provided by executing the transformation exp(A (t - 1)) provided as the output of the second affine transformation continuing means 6 to the second frame video image I 2 . As described above, the first principle is equivalent to the case in which the affine transformation executing means 7 obtains the video image at an arbitrary time using only the output of the first affine transformation continuing means 5 and not the output of the second affine transformation continuing means 6 used in the second principle. Therefore, it can be considered that the first principle forms one portion of the second principle, and the action of the second principle will be mainly explained hereinafter. For example, in a system for interpolating a real moving picture between frames, or in a system in which a video image at an arbitrary time t between time n and time n + 1 is produced for time-series video images given at discrete times 1, 2, ... n, n + 1,..., most specific feature of the present invention is that an approximately elliptic video image can be obtained from primary and secondary moments of an arbitrary video image and that the arbitrary video image can be transformed linearly by "rotation", "translation" and "contraction", according to the known affine transformation. Therefore, in the present invention, approximate ellipses of original video images at time n and time n + 1 are obtained, the affine transformation between the two video images is calculated from these two approximate elliptic video images, and a continued affine transformation {exp(At), exp(A (t - 1))} obtained by making the affine transformation continuous with regard to time based on the infinitesimal transformation (A) of the affine transformation is executed on the original video images, that is, the original video images at time n and time n + 1. Thus, the video image at time t is obtained. The principle of the present invention will be described hereinafter in detail. To follow the change in position and shape of the object, let us consider an affine transformation by which a first frame video image (video image at time n) is transformed into a second frame video image (video image at time n + 1). In the present invention, as shown in Figure 3A, centroid, principal axes of inertia, and moments about principal axes of inertia and moments about axes perpendicular to the principal axes of inertia are initially calculated from the first and second frame video images. Then, the affine transformation through which the characteristics of the first frame video image are transferred into those of the second frame video image is obtained. Further, this affine transformation is continued with respect to time and this continued affine transformation is executed, for example, onto the first frame video image to thereby estimate a video image at the arbitrary time t. The affine transformation for connecting characteristics is explained below. Assume that (X 1 , Y 1 ) and (X 2 , Y 2 ) are coordinates of centroid of the first video image I 1 (x, y) and the second video image I 2 (x, y),  1 and  2 are angles formed by the principal axis of inertia forms x axis, and EXP( λ 1 ), EXP( λ 2 ), EXP( κ 1 ) and EXP( κ 2 ) are moments about axes perpendicular to the principal axis of inertia. Herein, the following are established: X 1 = ∫∫ x I 1 (x, y) d x d y/∫∫ I 1 (x, y) d x d y Y 1 = ∫∫ y I 1 (x, y) d x d y/∫∫ I 1 (x, y) d x d y tan (2  1 ) - 2M 1xy /(M 1xx - M 1yy ) EXP (2κ 1 ) = (M 1xx + M 1yy ) + √ ((M 1xx - M 1yy ) 2 + 4 M 1xy ) 2 EXP (2λ 1 ) = (M 1xx + M 1yy ) - √ ((M 1xx - M 1xy ) 2 + 4 M 1xy ) 2 M 1XX = ∫∫ (x-X 1 ) 2 I 1 (x, y) d x d y/∫∫ I 1 (x, y) d x d y M 1YY = ∫∫ (y-Y 1 ) 2 I 1 (x, y) d x d y/∫∫ I 1 (x, y) d x d y M 1XY = ∫∫ (x-X 1 ) (y-Y 1 ) I 1 (x, y) d x d y/∫∫ I 1 (x, y) d x d y X 2 = ∫∫ x I 2 (x, y) d x d y/∫∫ I 2 (x, y) d x d y Y 2 = ∫∫ y I 2 (x, y) d x d y/∫∫ I 2 (x, y) d x d y t a n (2  2 ) = 2 M 2xy /(M 2xx - M 2yy ) EXP (2κ 2 ) = (M 2xx + M 2yy ) + √ ((M 2xx - M 2yy ) 2 + 4 M 2xy ) 2 EXP (2λ 2 ) = (M 2xx + M 2yy ) - √ ((M 2xx - M 2xy ) 2 + 4 M 2xy ) 2 M 2XY = ∫∫ (x-X 2 ) I 2 (x, y) d x d y/∫∫I 2 (x, y) d x d y M 2YY = ∫∫ (y-Y 2 ) I 2 (x, y) d x d y/∫∫I 2 (x, y) d x d y M 2XY = ∫∫ (x-X 2 ) (y-Y 2 ) I 2 (x, y) d x d y/∫∫I 2 (x, y) d x d y M 1XX and M 1YY are equations of central moments and are normalized by the division of ∫∫I 1 (x, y)dxdy, i.e., the zero-order moment m 0 . Thus, when the second video image I 2 (x, y), for example, is α times the first video image I 1 (x, y ), that is, the increased α times by the ambient illumination, the change of shape of the video image can be extracted regardless of the change of the brightness of the entirety. Further, assume that T(X, Y) is a transformation in which the video image is translated by X in the x axis direction and by Y in the y axis direction, that R(  ) is a transformation in which the image is rotated about an origin by  and that M( κ, λ ) is a transformation in which the video image is expanded and contracted exp(κ ) times in the x axis direction and exp(λ ) times in the y axis direction. These transformations are operations (affine transformations) which act not only on the video image space but also on the characteristic space. At that time, characteristic F 1 = (X 1 , Y 1 ,  1 , κ 1 λ 1 ) of the first video image and characteristic F 2 = (X 2 , Y 2 ,  2 κ 2 , λ 2 ) of the second video image are obtained by executing, as shown in Figure 3A, the following transformations on origin F 0 = (0, 0, 0, 0, 0) of characteristic space: F 1 = G 1 *F 0 G 1 = T(X 1 ,Y 1 ) R( 1 )M(κ 1, λ 1 ) F 2 = G 2 *F 0 G 2 = T(X 2 ,Y 2 ) R( 2 )M(κ 2 ,λ 2 ) From equations (1) and (2), it is clear that the characteristic F 2 of the second video image is equal to the result provided by executing the following transformation on the characteristic F 1 of the first video image: F 2 = G 2 G 1 -1 * F 1 The process for continuing the affine transformation is explained below. Initially, a continued affine transformation G(t) for transforming the moment characteristic F 1 of the first video image into characteristic F(t) of an arbitrary time from a time standpoint is obtained from the equation (3). Requesting that F(t) coincides with F 2 after the unit time yields the following equation: F(t) = G(t) * F 1 , F(0) = F 1 , F(1) = F 2 Although an infinite number of transformations satisfy equation (4), the one that satisfies the following condition is obtained: F(t) = exp(At)F 1 where A represents a transformation independent of time, Equation (5) means that the characteristic is successively changed by the infinitesimal transformation A. That is, where h is the step-size of time ≒ 0. Accordingly, G(t) is obtained by using infinitesimal transformation A which satisfies G 2 * G 1 -1 = exp(A) as follows: G(t) = exp(At) Equations (4) and (6) express the position and shape of an object at an arbitrary time t. From these equations, it is possible to understand the moving state of a physical object, such as a change in its position and shape, from a time standpoint. A process for estimating a video image at an arbitrary time is explained below. Assume that I 1 (x, y) is the first video image and I 2 (x, y) is the second video image. A video image at an arbitrary time is obtained by executing the affine transformation determined by equation (6) onto the first video image I 1 (x, y). More precisely, I (x,y;t) = G(t) * I 1 (x, y) However, G(1) * I 1 (x, y) and I 2 (x, y) do not always coincide. Accordingly, let us consider a method for estimating the video image I (x,y;t) at time t by overlapping G(t) * I 1 (x, y) and G -1 (1 - t) * I 2 (x, y) as: I (x,y;t) = (1 - t)G(t) * I 1 (x, y) + tG -1 (1 - t)* I 2 (x, y) In accordance with equation (8), I (x,y;t) coincides with the first video image at time t = 0 and with the second video image at time t = 1. Figure 3B is a flowchart used to explain the principle of the video image estimating system of the present invention. Referring to Figure 3B, in step S10, the moment characteristic F 1 of the video image is extracted from the first video image and in step S11, the moment characteristic F 2 is extracted from the second video image. In step S12, the affine transformation g(t) which connects these two characteristics is determined, and in the next step S13, the affine transformation is continued using equation (6). In the next step S14, the video image at an arbitrary time t is estimated from equation (8). As described above, according to the present invention, a video image at an arbitrary time t between time n and time n + 1 in a time sequence as one portion of the sequential change from the video image at time n to the video image at time n + 1 of video images at discrete times 1, 2, . . . n, n + 1, . . . , is estimated based on the infinitesimal transformation (A) of the affine transformation which connects the characteristics of two video images. Therefore, a video image at time t can be produced selectively. The present invention will be described more fully with reference to the following drawings. Figure 4 is a schematic block diagram showing the first embodiment of a video image data interpolating system which utilizes the moving video image estimation system of the present invention. Referring to Figure 4, an input unit 20 converts video images input in a time series into a digital video image. The digitized video image is input to a second video image memory 21, and the video image stored in the second video image memory 21 is then input to a first video image memory 22. Accordingly, the first video image memory 22 stores the video image at time n and the second video image memory 21 stores the video image at time n + 1. A moment calculating unit 23 calculates from the aforementioned equations primary and secondary moments of the video image at time n + 1 stored in the second video image memory 21. The calculated result is input to a second moment memory 24, and the moment characteristic stored in the second moment memory 24 is then input to a first moment memory 25. Accordingly, the first moment memory 25 stores the moment of the video image at time n and the second moment memory 24 stores the moment of the video image at time n + 1. From these moments, an affine transformation calculating unit 26 calculates from equations (1) to (3) an affine transformation which transforms the moment characteristic of the video image at time n into the moment characteristic of the video image at time n + 1. First and second affine transformation continuing units 27 and 28 calculate infinitesimal transformation A of the affine transformation within equation (5) and then calculate continuing affine transformation from a time standpoint from equation (6) or the following equation: G(t) = exp(A (t - 1)) where the variable t, indicative of time, uses time n as its starting point. A first affine transformation unit 30 and a second affine transformation unit 31 within an affine transformation executing unit 29 execute continued affine transformations on the video image at time n and the video image at time n + 1, respectively, and results thereof are synthesized by an image synthesizing unit 32 in accordance with equation (8). Figure 5 shows a hardware arrangement of the first embodiment of the system shown in Figure 4. The hardware of Figure 5 is different from that in the block diagram of Figure 4 in that in Figure 5 the input unit 20 is composed of a television camera 35 and an analog-to-digital converter 36; the moment calculating unit 23 is composed of a video image projection circuit 37, three projection memories 38 to 40 for storing projections to the X axis, the Y axis and the oblique axis, and a moment calculating circuit 41; the affine transformation calculating unit 26, the first affine transformation continuing unit 27 and the second affine transformation continuing unit 28 form one portion of a microprocessor 42; and the first affine transformation unit 30 is composed of a first affine transformation coefficient memory 43 and a first affine circuit 45; the second affine transformation unit 31 is composed of a second affine transformation coefficient memory 44 and a second affine circuit 46: The video image synthesizing unit 32 in Figure 4 is equivalent to the video image synthesizing circuit 47. Prior to the detailed explanation of the processing of the respective units of the first embodiment of the system shown in Figures 4 and 5, let us explain how to obtain the affine transformation G(t) given by equation (6). The affine transformation G(t) is given by equation (6) and A in this equation is determined as: G 2 * G 1 -1 = exp(A) However, it is not easy to directly calculate the above equation, and in actual practice, G(t) is calculated from the following mathematical relation. Affine transformation G given by G 2 * G 1 -1 = G is expressed by the following 3 x 3 matrix: In equation (9), L is a 2 x 2 matrix and B is a vector, and by this transformation, point (x, y) is moved according to: L 11 and B 1 are given by the following equations in accordance with the equations (1) to (3): The matrix L is expressed by matrix u of an infinitesimal transformation as: L = expu The shape of u is obtained later, and A which is expressed by the relation G = expA is expressed by the following equation: This equation can be confirmed in actual practice by calculating the following equation: expA = 1 + A + A 2 /2 + A 3 /3! + ... Incidentally, E in equation (14) is the unit matrix. Then, G(t) = expAt can be expressed as where ℓ (t) = exput. In equation (16), B and L are given already and hence, ℓ (t) = exput may be calculated. In this case, if eigenvalues of L are exp α and exp β , then L is expressed by utilizing regular matrix P as follows: Multiplying by P from the left of this equation and P -1 from the left thereof yields Since L = expu, Thus, Accordingly, In other words, respective matrix elements of ℓ(t) can be expressed by linear expression of exp αt and exp βt, i.e. by substituting ℓ(t) into the following equation: From the condition of ℓ(0) = identical transformation and ℓ(1) = L, r ij and s ij in the above equations are obtained. That is, From these relations, r ij and s ij are obtained, and equation (22) yields From equation (11), it is to be understood that the eigenvalues of L are exp( κ 2 - κ 1 ) and exp( λ 2 - λ 1 ), and these eigenvalues yield where h 1 1 =exp λ t (exp κ cos  sin t-sin (t- 1) ) /sin  h 1 2 =exp λ t (-exp (-κ) sin sin t/sin  h 2 1 =exp λ t (exp κsin sin t) /sin  h 2 2 =exp λ t (-exp (-κ) cos sin t-sin (t-1) ) /sin   = 2 - 1 , λ = (κ 2 - κ 1 + λ 2 - λ 1 )/2 κ = (κ 2 - κ 1 - λ 2 + λ 1 )/2 cos =cos h κ cos  In actual calculation, if ℓ(t) is calculated from equations (26) and (27), then all elements of G(t) are determined by equation (16). As described above, in actual practice, G(t) is calculated directly without calculating specific values of A. However, the obtained value becomes coincident with that of G(t) = exp At in which A as given by G = exp A is employed. Figure 6 is a flowchart used to explain the process in which the moment characteristic of the video image is extracted by the moment calculating unit. Referring to Figure 6, from the digital video image I(i, j), projection to the x axis is calculated in step S50, projection to the y axis is obtained in step S51, and projection to the oblique axis is obtained in step S52. In the equations for calculating these projections, as in the diagram used to explain the projection to the oblique axis of Figure 7, N represents the size of the video image, that is, the length of one side. When the video image is expressed as I (i, j), i and j take integer values from 1 to N. However, in the equation P X+Y(κ) ,κ takes integer values ranging from 1 to 2N. If κ is less than N, the projection to the oblique axis is obtained from the upper equation, while if κ is larger than N, the projection to the oblique axis is calculated from the lower equation. The projection to the oblique axis is equivalent to the process in which the sum of pixel data is calculated in the oblique direction as shown in Figure 7. Referring to Figure 6, zero-order moment m * is calculated in step S53, primary moment m X and m Y are calculated in steps S54 and S55 and centroid coordinates X and Y are calculated at steps S56 and S57. Further, secondary moments m XX , m YY and m XY are calculated in steps S58, S59 and S60, and central moments M XX , M YY and M XY are calculated in steps S61, S62 and S63. Then, the processing is ended. Incidentally, the equation of central moments M XX , M YY and M XY is substantially the same that earlier noted and, by way of example, the following equation is established: M 1 XX =∫∫ (x-X 1 ) 2 I 1 (x, y) d x d y/∫∫ I 1 (x, y) d x d y =∫∫ (x 2 -2 xX 1 +X 2 ) I 1 (x, y) d x d y/∫∫ I 1 (x, y) d x d y =∫∫ X 2 I 1 (x, y) d x d y/∫∫ I 1 (x, y) d x d y -2X 1 ∫∫ x I 1 (x, y) d x d y/∫∫ I 1 (x, y) d x d y+X 1 2 ∫∫ I 1 (x, y) d x d y/∫∫ I 1 (x, y) d x d y =m xx /m 0 -2X 1 X 1 +X 1 2 =m xx /m 0 -X 1 2 Figure 8 is a flowchart used to explain the process wherein affine transformation is determined by the affine transformation calculating unit. Referring to Figure 8, in steps S65, S66 and S67, by utilizing the centroid and central moment of the first video image, direction  1 of the principal axis of inertia, moment κ 1 of the axis perpendicular to the principal axis of inertia and moment λ 1 about the principal axis of inertia are calculated. Then, in steps S68, S69 and S70, by utilizing the centroid and central moment of the second video image, direction  2 of the principal axis of inertia, moment κ 2 of the axis perpendicular to the principal axis of inertia and moment λ 2 about the principal axis of inertia are calculated. Then, the affine transformation is calculated by using these elements at step S71, and the processing is ended. Figure 9 is a flowchart used to explain the process in which the first affine transformation is continued by the first affine transformation continuing unit. This continuing process is executed by using affine transformation coefficients , λ and κ, matrix element Lij and vector element Bi calculated by the affine transformation calculating unit in Figure 8. Referring to Figure 8, time t is set in step S73, matrix element hij is calculated by using the equation (27) in step S74 and matrix element lij(t) of the affine transformation given by the following equation is calculated in step S75: In step S76, element Pij of the matrix for calculating the vector element bi(t) in the above equation is calculated, and in step S77, the vector element bi(t) is calculated and then the processing is ended. Incidentally, the thus obtained b 1 (t) and b 2 (t) are the calculated results of (1 (t) - E) (L - E) -1 B in equation (16). Figure 10 is a flowchart used to explain the process wherein the second affine transformation is continued by the second affine transformation unit. Figure 10 is similar to Figure 9. Referring to Figure 10, time t is set in step S80. In step S81, (t - 1) is substituted for t, and then processes similar to those in steps S74 to S77 are performed in steps S82 to S85, and matrix element lij' (t) and vector element bi'(t) are calculated. Then, the processing is ended. Figure 11 is a flowchart used to explain the process wherein a video image at an arbitrary time is estimated by the video image synthesizing unit. Referring to Figure 11, at step S87, the affine transformation calculated in Figure 9 is executed on the first video image I 1 , and in step S88 the affine transformation obtained in Figure 10 is executed on the second video image I 2 . In step S89, the results of these executions are synthesized to obtain an interpolated video image at time t. Figures 12A to 12F show examples of the results of interpolated video images. Figure 12A shows the first video image, Figure 12F shows the second video image, and Figures 12B to 12E show interpolated video images at times between the first and second video images. The embodiment described above is a moving video image estimation system for interpolating an image at an arbitrary time between time n for a first video frame image and a time n + 1 for a second frame video image. A moving video image estimation system for predicting an image after a time n + 1, namely, an arbitrary time after the second frame video image, is explained hereinafter. A video image prediction can be performed by making time t after the first frame video image larger than "1", i.e. providing t > 1 in the above embodiment. In this case, first affine transformation unit 30 shown in Figure 4 executes transformation {exp (At)} as an output of the first affine transformation continuing unit on the first frame video image data I 1 stored in the first video image memory 22 and only the result of this transformation is output from video image synthesizing unit 32. However, it is preferable for a prediction of a video image after time t + 1 to be performed based on the second frame video image I 2 (x,y). Figure 13 is a block diagram of the structure of a second embodiment of the moving video image data estimation system according to the present invention. The present invention predicts a video image at an arbitrary time after time n + 1 based on the second frame video image data I 2 (x,y) at time n + 1. As in the previous embodiment, shown in Figure 4, this embodiment shown in Figure 13, comprises an input unit 90, a video image memory b91 corresponding to the second video image memory 21, moment calculating unit 92, moment memory a 94 and moment memory b93 corresponding to the first and second moment memories 25 and 24, affine transformation determining unit a95 corresponding to affine transformation calculating unit 26 and affine transformation continuing unit 97 corresponding to the first affine transformation continuing unit 27. The embodiment of Figure 13 executes the same operation as the embodiment of Figure 4. However, the embodiment of Figure 13 does not have a first video image memory 22 or a second affine transformation continuing unit 28, but instead has an affine transformation determining unit b96, which is not provided in the first embodiment shown in Figure 4. Because of this difference, an affine transformation executing unit 98 executes an operation different from the first embodiment shown in Figure 4. The first frame video image at time n is input to input unit 90 and is converted to digital data I 1 (x,y). The digital data I 1 (x,y)is stored in the image video memory b91 and is transferred to moment calculating unit 92. Moment calculating unit 92 calculates a moment characteristic F 1 (X 1 , Y 1 ,  1 , κ 1 , λ 1 ) from digital data I 1 (x,y) and transmits the calculated F 1 data to moment memory b93. Then, moment characteristic F 1 of the first frame video image at time n is stored in moment memory b93. When first frame video image data I 1 (x,y) stored in video image memory b91 is output to moment calculating unit 92, it becomes possible to input the second frame video image at time n + 1. The second frame video image at time n + 1 is first input to input unit 90 and converted to digital data I 2 (x,y), and then stored in video image memory b90. Therefore, the second frame video image memory data I 2 (x,y) at time n + 1 is stored in video image b91 and the first frame video image data I 1 (x,y) at time n is disposed of. The second frame video image data I 2 (x,y) stored in image memory b91 is then transferred to moment calculating unit 92. Moment calculating unit 92 calculates moment characteristic F 2 (X 2 , Y 2 ,  2 , κ 2 , κ 2 ). The calculated moment characteristic F 3 of the second frame video memory is transferred to moment memory b93 and at this point, the moment characteristic F 1 of the first frame video image stored in the moment memory b93 is transferred to moment memory a94 and stored therein. Therefore, the moment characteristic F 1 of the first frame video image is stored in moment memory a 94 and moment characteristic F 2 of the second frame video image is stored in moment memory b93. When respective moment characteristics F 1 and F 2 are stored in moment memories a and b (94 and 93), affine transformation determining units a and b (95 and 96), which receive moment characteristics F 1 and F 2 , are activated. Affine transformation determining unit a95 obtains affine transformation G which converts the moment characteristic F 1 of the first frame video image to the moment characteristic F 2 of the second frame video image. Affine transformation determining unit b96 obtains affine transformation G -1 , which converts the moment characteristic F 2 of the second frame video image to moment characteristic F 1 of the first frame video image. G -1 is an inverse transformation of the affine transformation G. Next, affine transformation G calculated by affine transformation determining unit a95 is transferred to affine transformation continuing unit 97, which obtains affine transformation G(t). The affine transformation G(t) is continuous with regard to time. It is obtained by equation (6) in the same manner as in the first embodiment. Finally, affine transformation executing unit 98 predicts the video image at an arbitrary time t after time n + 1 (t > 1). Affine transformation executing unit 98 receives affine transformation G(t) which is obtained by affine transformation continuing unit 97 and is continuous with time, affine transformation G -1 for converting the moment characteristic of the second frame video image to the moment characteristic of the first frame video image, and the second frame video image data I 2 (x,y). Affine transformation executing unit 98 utilizes the fact that I 1 (x,y) is approximately equal to G -1 *I 2 (x,y), and instead of using equation (7) uses the equation I (x,y;t) = G(t)G -1 (1) * I 2 (x,y) where the prediction is performed by using t > 1. G (t) at an arbitrary time t (t > 1) after time n + 1 is applied to video image obtained by executing the inverse transformation G -1 on the second frame video image at time n + 1. Figure 14 shows the system structure of a third embodiment of the present invention. The embodiment shown in Figure 14 simplifies the embodiment shown in Figure 13. Equation (28) can be transformed to the following equation (28)' by using the following equations (29) and (30). G(t 1 )G(t 2 ) = G(t 1 + t 2 ) G -1 (t) = G(-t) I (x,y;t) = G(t) G -1 * I 2 (x,y) Equation (31) for predicting a video image at an arbitrary time t of t > 1 can be obtained. Figure 14 shows the structure of a system for predicting a video image by using the equation (31). The embodiment shown in Figure 14 comprises input unit 90, video image memory b91, moment calculating unit 92, moment memory a94, moment memory b93, affine transformation determining unit a95 and affine transformation continuing unit 97, as in the embodiment shown in Figure 13. Thus, it performs an operation similar to that of the embodiment of Figure 13. Affine transformation executing unit 98 shown in Figure 14 is different from that in the embodiment of Figure 13 in that it executes the image prediction process in accordance with equation (31). The embodiment shown in Figure 14 performs the same operation as in the embodiment shown in Figure 13 until respective moment characteristics F 1 and F 2 of the first frame video image and the second frame video image are calculated. The first frame video image is first received by input unit 90 and converted to video image data I 1 (x,y). It is then stored in video image memory b91. Video image data I 1 (x,y) stored in video image memory b91 is then transmitted to moment calculating unit 92, which calculates the moment characteristic F 1 (X 1 , Y 1 ,  1 , κ 1 , λ 1 )of the first frame video image and stores it in the moment memory b93. Next, input unit 90 receives the second frame video image, and converts it into video image data I 2 (x,y), and stored in video image memory b91. At this time, the first frame video data I 1 (x,y) stored in video image memory b91 is disposed of. Video image data I 2 (x,y) stored in video image memory b91 is then transmitted to moment calculating unit 92 and moment calculating unit 92 calculates the moment characteristic F 2 (X 2 , Y 2 ,  2 , κ 2 , λ 2 ,) of the second frame video image and stores it in moment memory b93. At this time, the moment characteristic F 1 of the first frame video image stored in moment memory b93 at the time is transmitted to moment memory a94 and is stored therein. When respective moment characteristics F 1 and F 2 are stored in moment memories a, b (94 and 93), then affine transformation determining unit a95 is activated. Affine transformation determining unit a95 receives moment characteristics F 1 and F 2 and obtains affine transformation G to convert moment characteristic F1 to moment characteristic F 2 . Affine transformation G calculated by affine transformation determining unit a95 is transmitted to affine transformation continuing unit 97 and affine transformation continuing unit 97 obtains affine transformation G(t), which is continuous with time. This process can be achieved by using equation (6) as in the first embodiment. Finally, affine transformation executing unit 98 predicts the video image at an arbitrary time t after time n (t > 1). Affine transformation executing unit 98 receives affine transformation G(t), which is obtained by affine transformation continuing unit 97 and is continuous with time, and the second frame video image data I 2 (x,y). It thereby predicts a video image I(x,y;t) at an arbitrary time t(t > 1) after the second frame of the video image, in accordance with equation (31). At described in detail above, according to the present invention, a video image at an arbitrary time t between a time n and a time n + 1 is estimated as one portion of the sequential change of the image at time n to the video image at time n + 1 so that the video image at time t can be alternatively produced. Further, it becomes possible to estimate a video image at an arbitrary time t after time n + 1 by using a video image at time n and a video image at time n + 1. Therefore, the present invention can be applied to a wide variety of systems utilizing real moving pictures, such as television, animation, medical treatment industries and so on. Having described preferred embodiments of the invention with reference to the accompanying drawings, it is to be understood that the invention is not limited to those precise embodiments and that various changes and modifications thereof could be effected by one skilled in the art without departing from the scope of the invention as defined in the appended claims.

Claims: A moving video image estimation system comprising: moment calculating means (1) for calculating primary and secondary moments of a video image; affine transformation calculating means (2) utilizing an output of said moment calculating means to determine affine transformation (G) by which a first frame video image (I 1 ) is geometrically transformed into a second frame video image (I 2 ); first affine transformation continuing means (3, 5) for multiplying by an arbitrary time constant an infinitesimal transformation (A), given by the relation G = exp(A) with respect to said affine transformation (G), and obtaining an exponential transformation {exp(At)}; and affine transformation executing means (4, 7) for performing said continued affine transformation {exp(At)} to thereby obtain a video image of a time after an arbitrary time (t) from a time corresponding to said first frame video image. The moving video image estimation system according to claim 1, wherein said affine transformation executing means (4, 7) executes said continued affine transformation {exp(At)} on the first frame video image (I 1 ) to thereby obtain a video image occurring after said arbitrary time (t) from a time corresponding to said first frame video image. The moving video image estimation system according to claim 1, wherein moments of said first and second frame video images (I 1 ; I 2 ) are calculated to extract characteristics of respective video images, the affine transformation for coupling the characteristics of said first and second frame video image is determined, the exponential transformation {exp(At)} is obtained as the continued affine transformation after the infinitesimal transformation (A) of said affine transformation (G) is multiplied by the arbitrary time constant and said continued affine transformation {exp(At)} is executed on said first frame video image (I 1 ) to thereby assume a video image at the arbitrary time (t). The moving video image estimation system according to claim 1, wherein said affine transformation calculating means (2) calculates an inverse transformation (G -1 ) of said affine transformation (G) and said affine transformation executing means (4, 7) executes a transformation of the product {exp(At)}G -1 of said continued affine transformation and said inverse transformation on said second frame video image (I 2 ) to thereby obtain a video image after said second frame video image (I 2 ). The moving video image estimation system according to claim 4, wherein said moment calculating means (1) comprises an input portion (20) for receiving video image data, a video image memory (21, 22) for first storing said first frame video image and finally storing said second frame video image, a first moment memory (25) for storing a moment of said first frame video image and a second moment memory (24) for storing a moment of said second frame video memory; and said affine transformation calculating means (2) comprises a first affine transformation determining unit (27) for determining said affine transformation (G) and a second affine transformation determining unit (28) for calculating the inverse transformation (G -1 ) of said affine transformation. The moving video image estimation system according to claim 1, wherein said affine transformation executing means (4,7) executes a transformation {exp(A (t-1))} on said second frame video image (I 2 ) by using said continued affine transformation {exp(At)} to thereby obtain a video image after said second frame video image. The moving video image estimation system according to claim 6, wherein said moment calculating means (1) comprises an input portion (20) for receiving video image data, a video image memory (21, 22) for first storing said first frame video image (I 1 ) and finally storing said second frame video image (I 2 ), a first moment memory (25) for storing a moment of said first frame video image and a second moment memory (24) for storing a moment of said second frame video memory; and said affine transformation calculating means (2) comprises an affine transformation determining unit (26) for determining said affine transformation (G). A moving video image estimation system according to claim 1, further comprising: second affine transformation continuing means (6) for multiplying said infinitesimal transformation (A) by another arbitrary time constant and obtaining an exponential transformation {exp(A(t-1))}; whereby said affine transformation executing means (4, 7) being further adapted for obtaining a video image at the arbitrary time (t) between said first and second frames as a linear sum of results which is provided by executing a transformation {exp(At)}, provided as an output of said first affine transformation continuing means, onto said first frame video image (I 1 ) and of a result which is provided by executing a transformation {exp(A(t-1))}, provided as an output of said second affine transformation continuing means, onto said second frame video image (I 2 ). The moving video image estimation system according to claim 8, wherein an input unit (20) is provided for inputting video image data, a first video image memory (22) is provided for storing video image data of said first frame, a second video image memory (21) is provided for storing video image data of said second frame, a first moment memory (25) is provided for storing a moment of said first frame video image, a second moment memory (24) is provided for storing a moment of said second frame video image, a first affine transformation unit (27) is provided for performing a continued affine transformation output from said first affine transformation continuing unit onto the first frame video image stored in said first video image memory, a second affine transformation unit (28) is provided for performing a continued affine transformation output from said second affine transformation continuing unit onto the second frame video image stored in said second video image memory, and a video image synthesizing unit (32) is provided for calculating a linear sum of outputs of said first and second affine transformation units to thereby obtain a video image at an arbitrary time, wherein said system includes an affine transformation executing unit equivalent to said affine transformation executing means. The moving video image estimation system according to claim 9, wherein said input unit (20) comprises a television camera (35) for picking up a video image and an analog-to-digital converter (36) for converting an output of said television camera into a digital signal, said moment calculating unit (23) comprises of a video image projection circuit (37) for obtaining projections to the x axis, the y axis and the oblique axis of a video image; three projection memories for storing projections of the x axis, the y axis and the oblique axis output from said video image projection circuit; and a moment calculating circuit (41) for calculating a moment of a video image by utilizing outputs of said three projection memories, said affine transformation calculating unit (2) and said first and second affine transformation continuing units (27, 28) are formed of one portion of a microprocessor (42), said first affine transformation unit (27) is composed of a first affine transformation coefficient memory (43) and a first affine circuit (45), said second affine transformation unit (28) is composed of a second affine transformation coefficient memory (44) and a second affine transformation circuit (46), and said video image synthesizing unit (32) is a video image synthesizing circuit (47). The moving video image estimation system according to claim 9, wherein said moment calculating unit (23) calculates projections of the x axis, the y axis and the oblique axis; zero-order moment m 0 ; primary moments m X and m y ; coordinates X and Y of gravity of a video image; secondary moments m XX , m YY and m XY ; and central moments M XX , M YY and M XY . The moving video image estimation system according to claim 9, wherein said affine transformation calculating unit (26) utilizes coordinates of centroid X and Y and central moments M XX , M YY and M XY of said first and second frame video images output from said moment calculating unit to calculate a direction Θ of the principal axis of inertia of each video image, a moment κ about the axis perpendicular to said principal axis of inertia and a moment λ about said principal axis of inertia, thereby obtaining matrix element Lij and vector element Bi of matrix which transforms a coordinate (x, y) into a coordinate (x', y') as an affine transformation coefficient. The moving video image estimation system according to claim 9, wherein said first affine transformation continuing unit (27) determines the arbitrary time t, the matrix element lij(t) of affine transformation is calculated, and the vector element bi(t) of affine transformation is calculated. The moving video image estimation system according to claim 9, wherein said second affine transformation continuing unit (28) determines the arbitrary time t, t-1 is substituted for t, the matrix element lij'(t) of affine transformation is calculated, and the vector element bi'(t) of affine transformation is calculated. The moving video image estimation system according to claim 9, wherein said affine transformation executing unit executes the affine transformation on said first frame video image by using the elements lij(t) and bi(t) calculated by said first affine transformation continuing unit (27), said affine transformation executing unit executes the affine transformation on said second frame video image by using the elements lij'(t) and bi'(t) calculated by said second affine transformation continuing unit (28), and executed results of said two affine transformations are synthesized to obtain an interpolated video image.

Inventors: SATOSHI NAOI, TAKASHI TORIU, SATOSHI, NAOI, TAKASHI, TORIU
Applicants: FUJITSU LTD, FUJITSU LIMITED
Application year: 1991
Main classification: G06T7
IPC codes: G06T3, G06F17, H04N5, G06T7, H04N7