In many application environments, an image sequence is encoded (or compressed) to reduce the total amount of data needed to represent the image sequence. The compressed data may then be stored or transmitted more efficiently than the original uncompressed image sequence data. The image sequence may be any sequence of images, including a sequence of video image frames and a sequence of still images. Multiple view image sequences are sequences of images corresponding to different views of a scene; the images may be captured by a single camera positioned at different viewpoints, or the images may be captured by multiple cameras positioned at different locations relative to the scene to capture the scene from different viewpoints.
Image compression methods typically fall into one or more of three main image compression classes: spectral redundancy reduction, spatial redundancy reduction, and temporal redundancy reduction. Spectral redundancy reduction methods typically reduce the amount of image data by discarding spectral data that are not strongly perceived by human eyes. Spatial redundancy reduction methods reduce higher spatial frequency components in the original image data. For example, transform coding is a common spatial redundancy compression method that involves representing an image by a set of transform coefficients. The transform coefficients are quantized individually to reduce the amount of data that is needed to represent the image. A representation of the original image is generated by applying an inverse transform to the transform coefficients. Temporal redundancy reduction methods compress a sequence of images by taking advantage of similarities between successive images. Temporal redundancy may be reduced, for example, by transmitting only those movements or changes in a given image that permit accurate reconstruction of the given image from another image (e.g., a previously received video image frame).
Various different standards of image sequence compression have been developed, often based on block-matching methods. Block-matching methods initially divide a target image (or frame in the case of video image data) to be compressed into an array of blocks (or tiles). Motion data and motion compensation difference data are generated for each block based on a set of data in a reference image (e.g., in a prior video frame) that is similar to the block. In a typical approach, the target image is completely divided into contiguous blocks and sets of pixels in the reference image that best match each block are identified. The target image is reconstructed by accessing and manipulating portions of the reference image. The motion data represents an amount of movement that repositions a suitable part of the reference image to reconstruct a given block of the target image, and the motion-compensated difference data represents intensity adjustments that are made to individual pixels within the set of data from the reference image to accurately reproduce the given block of the target image.
Various methods for computing motion vectors between blocks of a target image and corresponding blocks of a reference image have been proposed. In a typical block matching approach, a current block is compared with all the blocks of like size in a search window superimposed on the reference image. Typically, image blocks of the target image and the reference image are compared by calculating an error function value for each possible match. The motion vector with the smallest error function value is selected as the best matching motion vector for a given target image block. Exemplary block matching error functions are the sum of the absolute values of the differences of the pixels between matched blocks and the sum of the squares of the differences. Motion estimation typically requires a significant portion of the computational resources needed to implement any given image sequence compression method.
The invention features systems and methods of estimating motion for compressing multiple view images.
In one aspect, the invention features a machine-implemented method of encoding a target image of a scene captured at a first image plane. In accordance with this inventive method, a transformation is computed. The transformation maps at least three noncollinear points substantially coplanar on a scene plane in the target image to corresponding points in a references image of the scene captured at a second image plane different from the first image plane. At least one point in the target image off the scene plane and at least one corresponding point in the reference image are identified. A motion between the target image and the reference image is estimated based on the computed transformation and the identified corresponding off-scene-plane points. The target image is encoded based at least in part on the estimated motion.
The invention also features an apparatus and a machine-readable medium implementing the above-described method.
Other features and advantages of the invention will become apparent from the following description, including the drawings and the claims.
In the following description, like reference numbers are used to identify like elements. Furthermore, the drawings are intended to illustrate major features of exemplary embodiments in a diagrammatic manner. The drawings are not intended to depict every feature of actual embodiments nor relative dimensions of the depicted elements, and are not drawn to scale.
I. Overview
The image encoding embodiments described in detail below provide an efficient way to estimate motion between a target image and a reference image that may be incorporated readily into any image compression apparatus or process. These embodiments leverage geometric information computed from certain readily obtainable points in a scene to reduce the size of the search space that is searched to identify points in the reference image corresponding to points in the target image.
In the implementation of
In the implementation of
The encoder 32 and the decoder 36 may be incorporated in the same device or in difference devices. In general, encoder 32 and decoder 36 are not limited to any particular hardware or software configuration, but rather they may be implemented in any computing or processing environment, including in digital electronic circuitry or in computer hardware, firmware, device driver, a codec (e.g., an MPEG video codec), or software. The encoder 32 and the decoder 36 may be embedded in the hardware of any one of a wide variety of electronic devices, including desktop computers, laptop computers, and portable electronic devices (e.g., digital still image camera, digital video cameras, mobile phones and personal digital assistants). In some implementations, each of encoder 32 and decoder 36 is implemented, at least in part, in a computer process product tangibly embodied in a machine-readable storage device for execution by a computer processor. In some embodiments, encoder 32 and decoder 36 preferably are implemented in a high level procedural or object oriented processing language; however, the algorithms may be implemented in assembly or machine language, if desired. In any case, the processing language may be a compiled or interpreted language. Suitable processors include, for example, both general and special purpose microprocessors. Generally, a processor receives instructions and data from a read-only memory and/or a random access memory. Storage devices suitable for tangibly embodying computer process instructions include all forms of non-volatile memory, including, for example, semiconductor memory devices, such as EPROM, EEPROM, and flash memory devices; magnetic disks such as internal hard disks and removable disks; magneto-optical disks; and CD-ROM. Any of the foregoing technologies may be supplemented by or incorporated in specially designed ASICs (application-specific integrated circuits).
II. Estimating Motion Between Pairs of Multiple View Images
As mentioned above, the image encoding embodiments described herein provide an efficient way to estimate motion between a target image and a reference image that leverages geometric information computed from certain readily obtainable points in a scene to reduce the size of the search space that is searched to identify points in the reference image corresponding to points in the target image.
Referring to
Since Hπ has eight degrees of freedom up to scale, Hπ may be determined from four sets of corresponding points, x(i)1⇄x(i) 2, where i has an integer value from 1 to 4. In some implementations, more than four sets of correspondence points may be used in a singular value decomposition process to determine a least-squares solution for H.
l(x)2=x′2×x2 (3)
where x is the cross product operator. Since x′2=Hπxl,
l(x)2=(H,πx1)×x2 (4)
Similarly, the epipolar line l(y)2 is given by:
l(y)2=y′2×y2 (5)
and since y′2=Hπy1,
l(y)2=(Hπy1)×y2 (6)
The epipole e2 may be computed from l(x)2×l(y)2 which are given by equations (4) and (6), if the corresponding off-scene-plane points x1⇄x2 and y1⇄y2 are known. The epipole e2 then may be used to determine epipolar lines in the image plane 18. With this information, any given point w2 in the second image plane 18 corresponds to the homographic projection of a corresponding point w1 in the first image plane 16 plus a possible displacement along the epipolar line l(w)2. That is,
w2=(Hπw1)+Θ·e2 (7)
where Θ is a scalar.
In some implementations, the full projection model of equation (2) is approximated by an affine approximation model. The accuracy of this model improves as the distances of points in the scene to the centers of the camera(s) become much greater than the distances between points in the scene. In this model, the transformation is given by:
Since Hπ has six degrees of freedom up to scale in the affine approximation model, Hπ may be determined from three sets of corresponding points, x(i)1⇄x(i)2, where i has an integer value from 1 to 3. In the affine approximation model, the epipolar lines are assumed to be parallel, in which case only a single epipolar line is needed to determine the epipole e2. Indeed, if an epipolar line l2=(a,b,c), then the corresponding epipole e2=(−b,a,0). The epipolar line l2 may be determined from images of a single off-scene-plane point (e.g., point X).
Referring back to
Encoder 32 also identifies at least one point (e.g., point X) in the scene that does not lie on the planar surface used to compute the transformation Hπ. This point is imaged by a point (e.g., point x1) in the target image plane 16 and a corresponding point (e.g., point x2) in the reference image plane 18 (step 86). Only one pair of corresponding off-scene-plane points is used in the affine approximation model, whereas two such pairs of points are used in the full projection model of equation (2). The off-scene-plane points may be determined automatically using any automatic point correspondence matching process, or a user may identify the corresponding off-scene-plane points manually.
Based on the computed transformation and the identified corresponding off-scene-plane points, the encoder 32 estimates motion between the target image and the reference image (step 88). In some implementations, using equations (4) and (6) the encoder 32 computes the epipole e2 in the reference image plane 18 based on the identified corresponding off-scene-plane points. The epipole e2 is used to reduce the search space for blocks in the reference image that match corresponding blocks in the target image to a single dimension, as explained in detail below.
Referring to
For each block, the encoder 32 computes a respective motion vector representing motion between a target image block and a reference image block (step 92). In this process, the encoder selects the next block of the target image as the next block (step 94). For each target image block, the encoder defines a single-parameter search space relating the target image points to a block of points in the reference image (step 96). The single-parameter search space is defined with respect to the following motion model, which is derived from equation (7):
bi′=(Hπ·bi)+Θ·e2 (9)
where Hπ is the computed transformation, Θ is a scalar, e2 is the respective epipole, bi is any point within the target image block being encoded, and bi′ is the corresponding point in the reference image lying on the epipolar line li. In this way, the search space for the points within block bi′ in the reference image is parameterized by a single parameter Θ. The encoder 32 searches each single-parameter search space for a block of points in the reference image corresponding to the current target image block (step 98). This process is repeated for each block in the target image (step 100).
In some implementations, a matching score is generated for each of multiple different values for Θ and the block bi′ corresponding to the value of Θ producing the highest matching score is selected as the reference image block corresponding to the target image block. The matching score may be generated based on a correlation matching measure (e.g., the sum of the absolute values of the differences of the pixels between matched blocks, or the sum of the squares of the differences) or to some other point match (or correlation) measure. The value of Θ defines a motion vector since specification of Θ determines the reference image block matching the corresponding target image block.
For each block of the target image and for each computed motion vector, the encoder 32 also computes motion compensation difference data (step 102). The motion compensation difference data represents intensity adjustments to the points of the reference image block needed to reproduce the intensity values of the points of the target image block.
Referring back to
Other embodiments are within the scope of the claims.
Number | Name | Date | Kind |
---|---|---|---|
5535288 | Chen et al. | Jul 1996 | A |
5742710 | Hsu et al. | Apr 1998 | A |
5793985 | Natarajan et al. | Aug 1998 | A |
5953458 | Pirson et al. | Sep 1999 | A |
5963664 | Kumar et al. | Oct 1999 | A |
6144701 | Chiang et al. | Nov 2000 | A |
6205241 | Melon | Mar 2001 | B1 |
6226396 | Marugame | May 2001 | B1 |
6266158 | Hata et al. | Jul 2001 | B1 |
6339617 | Ueda | Jan 2002 | B1 |
6542547 | Wong | Apr 2003 | B1 |
6584226 | Chaddha et al. | Jun 2003 | B1 |
6628419 | So et al. | Sep 2003 | B1 |
7286689 | Damera-Venkata et al. | Oct 2007 | B2 |
20020001406 | Kochi et al. | Jan 2002 | A1 |
20020110275 | Rogina et al. | Aug 2002 | A1 |
20030044048 | Zhang et al. | Mar 2003 | A1 |
20040081239 | Patti et al. | Apr 2004 | A1 |
20040131267 | Adiletta et al. | Jul 2004 | A1 |
20040165781 | Sun | Aug 2004 | A1 |
Number | Date | Country |
---|---|---|
WO 0343342 | May 2003 | WO |
Number | Date | Country | |
---|---|---|---|
20050169543 A1 | Aug 2005 | US |