|Publication number||US20020167512 A1|
|Application number||US 09/851,445|
|Publication date||Nov 14, 2002|
|Filing date||May 8, 2001|
|Priority date||May 8, 2001|
|Also published as||CN1241419C, CN1462561A, DE60209365D1, DE60209365T2, EP1393581A1, EP1393581B1, US6965379, WO2002091754A1|
|Publication number||09851445, 851445, US 2002/0167512 A1, US 2002/167512 A1, US 20020167512 A1, US 20020167512A1, US 2002167512 A1, US 2002167512A1, US-A1-20020167512, US-A1-2002167512, US2002/0167512A1, US2002/167512A1, US20020167512 A1, US20020167512A1, US2002167512 A1, US2002167512A1|
|Inventors||Mi-Suen Lee, Tomas Brodsky, Daphna Weinshall, Miroslav Trajkovic|
|Original Assignee||Koninklijke Philips Electronics N.V..|
|Export Citation||BiBTeX, EndNote, RefMan|
|Patent Citations (5), Referenced by (6), Classifications (35), Legal Events (5)|
|External Links: USPTO, USPTO Assignment, Espacenet|
 A. Field of the Invention
 The invention relates to the field of synthesizing video material, particularly with a view to enhancing a three dimensional appearance of that material.
 B. Related Art
 Synthesizing the image of a 3-D scene, as it would be captured by a camera from an arbitrary viewpoint, is a well-known research problem and several different approaches have been proposed. Given a complete 3-D model, it is possible to render the scene from any viewpoint. Rendering techniques are well known in the field of Computer Graphics (as described, for example, in J. D. Foley, A. van Damm, S. K. Feiner, and J. F. Hughes, “Computer Graphics-Principles and Practice”, Addison Wesley, Second Edition, 1992, Chapter 14). However, generation of the model is a tedious, error-prone and labor intensive task.
 On the other hand, some researchers have tried to infer 3-D information directly from images. Various methods have been developed for recovering both the shapes of the objects and the camera motion from a sequence of images (for example, O. Faugeras, “Three-Dimensional Computer Vision: a Geometric Viewpoint”, MIT press, 1993). These methods are typically computationally intensive. Furthermore, in many cases the input sequence had been captured by a pan-tilt-zoom camera, preventing recovery of the 3D model of the scene.
 Another approach is to perform image based rendering (as in L. MacMillan: “An Image Based Approach to Three-Dimensional Computer Graphics”, Ph.D. dissertation, University of North Carolina, 1997), in which case no explicit 3-D model of the object and the environment is required. These methods avoid reasoning in 3-D by using projective constraints and invariants.
 It is an object of the invention to synthesize multiple views of a scene starting from a monocular video sequence.
 This is achieved in that for each input image, two or more new views of the scene are generated. Such views are suitable for display on a 3D-TV device or any other suitable viewing apparatus. The input images are segmented into moving objects and a stationary background. The background may be further segmented into multiple regions. Applying appropriate transformations to each region generates new views. Advantageously, simple domain knowledge is used to create good approximations for the views, even when the necessary information is not available from the image sequence.
 Further objects and advantages will be apparent in the following.
 The invention will now be described by way of non-limiting example with reference to the following drawings.
FIG. 1 shows a system in which the invention can be used.
FIG. 2 shows a flowchart of the invention.
FIG. 3 shows a schematic of an image to be processed in accordance with the invention.
FIG. 4 shows segmentation of moving obj2ects.
FIG. 5 shows a flowchart relating to use of domain knowledge and further segmentation of the stationary background into multiple regions.
FIG. 6 illustrates aspects of how the new camera positions are derived.
FIG. 7 illustrates how transformations of the segmented objects are computed
FIG. 8 shows how the new views are created.
 FIGS. 9-13 show some formulae that will be used in the description of the preferred embodiments.
FIG. 1 shows a typical 3-D image viewing apparatus, sometimes also referred to as a virtual reality device. The device includes a processor 101, a memory 103, and, optionally various peripherals, such as a display 102, a network connection 104, or other types of peripherals not shown, such as a printer, keyboard, mouse, microphone, loudspeakers, camera, and/or remote control device. The connections may be wired or wireless. The processor may be of any type including embedded in a PC, television, or set top box. The device may access additional processing or memory power via the network connection 104, which may be any suitable network connection such as to the Internet, cable television, and/or a local area network (LAN). Commonly the 3-D effect is achieved using a viewing device 105 that creates a stereo image, with the difference between the images shown to each of the two eyes of the viewer giving the illusion of three dimensions.
FIG. 2 shows a flowchart of the method according to the invention. The method may be executed by the processor 101 or at a remote processor accessed via network connection 104. Though the preferred embodiment is illustrated as software, it can equally well be implemented in hardware by the skilled artisan.
 The invention will also be described with respect to the schematic image of FIG. 3. FIG. 3 shows a tennis game. Two players, 303 and 302, play on a court 304, in front of the stands 301 where the audience sits. The court 304 and stands 301 are known types of objects within a knowledge domain corresponding to tennis games.
 The input is a sequence of images I1, I2, I3, . . . , for instance, taken during a tennis match. FIG. 2 shows the steps that are applied to generate two stereoscopic views of the scene. The letter “L” is used to refer to images and transformations relating to the left view of the stereoscopic output, while the letter “R” is used to refer to images and transformations relating to the right view.
 At 201, moving objects are segmented from the stationary background. In this case, players 302 and 303 will be segmented. The details of the segmentation method will be discussed with respect to FIG. 4. Then at 202, domain knowledge is used to further segment the background into individual regions. More about box 202 will be discussed with respect to FIG. 5.
 At 203, the new views are chosen. Details will be explained with respect to FIG. 6. The operations of boxes 202 and 203 may be performed in parallel.
 At 204, for each region m, transformations TLm and TRm are computed. FIG. 7 illustrates the method.
 Then, at 205, the segmented regions are warped to create the left view, using transformations TLm. The regions are also warped to create the right view, using transformations TRm. Details are discussed with respect to FIG. 8.
 Segmentation of Moving Objects
FIG. 4 is a flowchart showing the steps performed in box 201. In boxes 402, 405, and 409, the symbol “•” is used to indicate the application of a homography to an image.
 4a. At 401, for each input image Ik, estimate the homography transformation Hk so that when Hk is applied to frame Ik, the image of the tennis court in the transformed frame is aligned with the image of the tennis court in the first frame I1.
 A homography is represented by a 3×3 matrix H, which maps a point (x, y) in the image to the point represented by formula (1) from FIG. 9, where the hij homography parameters are values from the homography matrix H. This is well known in the art, for example in O. Faugeras, Three-dimensional Computer Visions referenced above.
 The alignment of input images is done by computing spatio-temporal derivatives of the image sequence, estimating normal flow, and linearly estimating the homography parameters. One method for achieving this is found in R. Pless et al, “Detecting Independent Motion: The Statistics of Temporal Continuity”, Trans. IEEE PAMI, vol. 22, No. 8, August 2000.
 Alternatively, feature points (such as corners) and/or lines could be located and matched between the input images and the homography could be estimated from the matched features, as shown in U.S. Pat. No. 5,629,988 issued to Burt et al.
 At 402, Jk denotes the result of applying transformation Hk to the image Ik.
 4b. From domain knowledge, it is known that in most tennis sequences the camera is allowed to pan, tilt, and zoom, but otherwise is not moved around the court. For a pan-tilt-zoom camera the transformation computed in box 401 aligns not only the tennis court, but also any other stationary feature in the scene. Only the moving players are not aligned.
 The transformed images Jk from step1 can be combined to create a mosaic of the scene at 403. For each pixel in the mosaic, there are one or more pixels from images Jk. In order to filter out moving objects (the players), a median filter is applied to the multiple values at each pixel and that median value is put into the mosaic M.
 The median value is computed as follows. Suppose there are values X1, X2, . . . , XK for a certain pixel. The values are sorted and the sorted values are denoted y1, y2, . . . , yK, where y1<=y2<=. . . <=yK. The median value is then yK/2.
 An alternative way to create mosaics can be found in U.S. Pat. No. 5,629,988.
 4c. When an inverse transformation (Hk −1) is applied to mosaic M and the result is cropped appropriately, a background image Bk is obtained at 405. The background image is the same as the input image Ik, except that the moving objects are removed and replaced by pixels from the background, which were visible in other views.
 4d. Then the input image Ik is compared with the corresponding background image Bk to create a mask image Mk at 406, whose pixel values are 1 when the difference between Ik and Bk is larger than a threshold and 0 otherwise. Then at 407 from image Ik those pixels are extracted that are set to 1 in the mask image Mk and the remaining pixels are set to be black. The result is an image Fk of moving objects.
 Segmentation of Stationary Background
FIG. 5 shows an expansion of box 202. Domain identification can be obtained at 501 from electronic program guide data or by analyzing the closed captioning data in broadcast video. Even if no domain identification is available, it is possible to test whether the input video is from the correct domain by detecting the distinguishing features (box 502). If the expected features are detected, the method can proceed even without domain identification.
 Domain knowledge may also be supplied, for example, through MPEG stream meta-data, which could be available in digital TV broadcasts. If only domain identification is to be transmitted, then domain knowledge will have to be stored in advance in memory 103.
 In the example of FIG. 3, the domain knowledge would include the shape of the court 305, the fact that the center net 304 is vertical, and certain assumptions about the shape of the stands 301.
 At box 502, distinguishing features of the domain are detected. In the case of the tennis match, a distinguishing feature would be the lines of the court 305. An additional feature would be the knowledge that the color of the tennis court is uniform. The features are detected as follows:
 5a. Locate lines in the image. There are many known methods for detecting lines in images, for example.
 5b. Locate large regions of uniform color and select only those lines found in step 5a that are within such large regions.
 5c. Match detected lines with line segments in the stored model. This is done by randomly choosing two or more detected line segments and matching them with all segments in the model. Each matching determines a transformation between the image and the model (this transformation is a homography, because the court is planar).
 5d. For each homography obtained in step 5c, apply the transformation to the input image and find how many line segments in the transformed image match line segments in the model as well as how well such segments are aligned.
 5e. Steps 5c. and 5d. are repeated and a final transformation is chosen. Under the final transformation, the most segments are matched and have the best alignment between the transformed input and the model.
 At 503, the background is further segmented into the tennis court and other regions. In the case of the tennis match, these features might include the net, the umpire stand and the audience stands 301. The regions are segmented based on domain knowledge. For example, the net is in a known position with respect to the lines on the court. The umpire stand is known to be on the court itself and its color is different from the color of the court, so it can be segmented based on color. More about the type of segmentation necessary here can be found in Vishvjit S. Nalwa, A Guided Tour of Computer Vision (Addison-Wesley, 1993) section 3.3, “Image Segmentation”.
 As a result of this step, one mask image is obtained per region (analogously to the mask image obtained in 406). For each region, those mask pixels that correspond to pixels in that region are set to 1 and the other pixels are set to 0.
 Determine Transformations of the Tennis Court
 In FIG. 6, the transformations to be applied to the image of the tennis court are determined. By tennis court, only the ground is implied herein. The transformations for the other segmented background regions are determined in the next section: “Determine transformations of other regions” The transformations are homographies, since the court is planar. There are two alternative ways of determining HL and HR.
 6a. If the matching between the model and the input image is available (601, from box 502), the intrinsic and extrinsic calibration parameters (with respect to the court) can be determined 602. There are known techniques for doing so, for example in Emanuele Trucco, Alessandro Verri: “Introductory techniques for 3-D Computer Vision”, Prentice Hall, 1998.
 6b. In a fully calibrated case, any new desired camera position can be chosen 603. For example, the new cameras can be placed closer to the court to enhance 3D perception. The homographies HL and HR can be determined in closed form from the old and the new camera parameters 604. How to do this is described in the Trucco & Verri book referenced above.
 6c. If the system is not fully calibrated (calibration is not desired), heuristics 605 are used to determine HL and HR. In this particular example, suppose the net and the far baseline are approximately horizontal in the input image. For the bottom edge of the net, the Y-coordinate will be denoted by y1 and the width in the image by w1. For the far baseline, the Y-coordinate will be denoted by y2 and the width by w2.
 Because the two line segments have the same length in the 3D world, the ratio of w1/w2 is the same as the ratio of disparities for the two lines, in accordance with formula (13).
 Similarly, the Y-coordinate of the bottom scan line will be denoted in the image by yB. For each viewing device, the disparities should be smaller than a certain maximum, denoted by dMAX. The maximum dMAX might be, for example, 20 pixels.
 Utilizing formula (4), the values for SL, SR, dL and dR may be derived. Of course, this is just one of a possible set of solutions. At 606, by choosing SR=0 and dR=0, the right view can be made identical to the input view. To obtain SL and dL, the two linear equations at (14) are solved, one obtained from (13) and one from the maximum disparity constraint.
 The values y1, y2, w1, w2 and yB are measurements in the input image and dMAX is a known parameter of the viewing system, so SL and dL are the only unknowns in the equations above.
 Consequently, the homographies of formulae (2) and (3) of FIG. 9 can be used 607.
 These two homographies do not change the Y coordinate and thus create rectified stereo views, where corresponding features lie on the same scan lines in both images.
 The heuristic described above is only an example. If other known features on the court can be detected, analogous methods can be used to determine suitable homographies HL and HR.
 Determine Transformations of Other Regions
 In addition to the tennis court, the scene contains moving objects (segmented in 407) and segmented parts of the background (from 503). With reference to FIG. 7, it will now be explained how transformations TLm and TRm can be obtained for these regions. The procedure is explained for the segmented player regions, but the same (or quite similar) procedure can be applied to other regions.
 If the new view is similar to the original view (which is usually the case), the player is approximated as being in a fronto-parallel plane. This is a valid approximation considering the distance between the player and the camera. The TLm transformation is determined with respect to a bounding rectangle 701 around the region containing the foreground object 303. The lower left corner is denoted as A, the lower right comer as B, and the mid-point between A and B as C. The player 303 is assumed to be standing on the court and therefore the lower line of the bounding rectangle 701 is assumed to lie on the ground. The tennis court is transformed by homography HL. Specifically, homography HL transforms points A, B, C to points A′, B′, C′, respectively, in the new view.
 TLm must be compatible with HL: Point C should be transformed to C′ and the horizontal distance between transformed points A, B should equal the horizontal distance between points A′, B′. Finally, the aspect ratio (width divided by height) of the bounding box should remain constant. Consequently, formula (5) is derived.
 Formulae (6), (7), (8), (9), (10), and (11) of FIG. 9 are variable definitions for formulae (12) of FIG. 10. In each of the formulae (6)-(11) x and y coordinates are defined for points in FIG. 7 and their transforms.
 The transformation for the right view TRm, per formula (18) can be derived analogously, The variables appearing in the formula for TR are defined in formula (19). These in turn are derived from the assumption that homography HR transforms points A, B, C, to points A″, B″, C″, where these points are defined in formulae (15), (16) and (17) respectively.
 The objects (for example the tennis players) appear to “stand out”, because they are transformed differently from the tennis court (the ground). The transformations of the tennis court (the HL and HR homographies) are setup so that points higher in the image appear farther from the viewer. On the other hand, all the points on the tennis player appear at approximately the same distance from the viewer.
 The transformations TLm and TRm differ from the transformations HL and HR that are applied to the ground portion of tennis court. Consequently, the tennis player appears to “stand out” from the background. More specifically, the HL and HR transformations are generally constructed so that points on the far side of the court have smaller disparities and thus appear farther from the viewer. This creates a correct, slanted appearance of the court. On the other hand, the TLm and TRm transformations make all the points of the tennis player appear at approximately the same distance.
 The player's feet have approximately the same disparity as the nearby points of the court. Also, the disparity of the player's head will be approximately the same. However, points on the court that appear close to the player's head have smaller disparities and this makes the player “stand out” from the background. Other features, such as the umpire stand, are transformed analogously to the tennis players and thus also appear to “stand out” from the background.
 Creating the New Views
 Finally, the new views are created as described in FIG. 8. In previous steps, the input image was segmented into multiple regions and a transformation was computed for each region. In FIG. 8 there are two paths. The left path corresponds to generation of the left view and the right path corresponds to generating the right view. The two paths are shown as being executed in parallel; however, they might equally well be executed sequentially.
 At 801, the respective left view transformation TLm is applied to the image region, Ik, and also to each mask, Mk, to yield transformed image region ILm and transformed mask MLm. HL is also applied to the background mosaic BK, created in box 405. For the ground portion of tennis court, the transformation TLm is the homography HL.
 At 802, the transformed regions are combined to create the new left view. Whenever a pixel in the transformed mask MLm is 1, the corresponding pixel in ILm is copied to the output image.
 At 803, conflicts resulting from overlapping objects are resolved.
 If for some a, b, the objects overlap, i.e., MLa(pixel)=1 and MLb(pixel)=1, then Lk(pixel) should be set to the pixel IL(pixel) from the closest object, with object ordering defined below.
 Object Ordering:
 For each object, the bottom edge should be taken to yield its position on the reference plane. Closer objects then have the bottom edge lower in the input image. The reference plane (the tennis court) is ordered as the most distant objects.
 Also at 804, optionally, holes may be filled, where holes are pixels that are not covered by the segmented objects. If the black pixels are on the boundary, they should be left untouched, because they most likely correspond to locations never seen in the scene. On the other hand, the appropriate values from the background mosaic can fill holes inside the image, if such values are available.
 Boxes 801′, 802′, 803′, and 804′ are analogous to 801-804, respectively, with R, i.e. the right view, replacing L, i.e. the left view.
 The example of FIGS. 4-8 is worked out in terms of a device that uses stereo images to create a 3-D effect. However, the method is extendable to an implementation with an arbitrary number of views. The computation of the right view, that is, computation of HR, TRm etc. was completely independent of the computation of the left view (HL, TLm, etc.). Thus it is straightforward to generate N different views of the scene by repeating the algorithm N times. Typically, the N views would be generated according to the requirements of a particular display device.
 Alternate Embodiment for Transformation of Background Pixels
 An alternate embodiment for dealing with background areas can be derived from use of a second image of the same scene. In the following, it will be assumed that moving objects are already segmented from the image. The operations below may be used to estimate the appearance of stationary objects not on the reference plane (the tennis court), in accordance with the following
 1. First correspondences between pixels in the two images are established, manually, or automatically, see “Chapter 6 Stereo Vision”, Three Dimensional Computer Vision, Olivier Faugeras, MIT Press, 1993.
 2. The first images are registered up to an affine transform:
 a. two corresponding pairs of parallel lines on a plane are identified in accordance with formulae (20) and (21). The copairs of lines are denoted by PL1=[L11, L12] and PL2=[L21, L22], where L11, L12, L21, L22, are lines defined by endpoints Lij=[(sxij, syij), (exij, eyij)].
 b. The vanishing point (in homogeneous coordinates) for each pair of parallel lines is computed in accordance with formulae (20) and (21), where the notation of the cross product of a and b is shown in formula (22).
 c. The transform Ha that moves vp1=[x1 y1 w1] and vp2=[x2 y2 w2] to infinity is computed in accordance with formula (23).
 d. Apply Ha to the whole image, i.e., pixel according to formula (24) is moved according to formula (25).
 3. Find 4 or more points on a plane in the first image and their corresponding points in the second image. In accordance with formula (26) the transform Hb is computed. This transform moves these points in second image q1 q2 q3 q4, . . . to align with those in the first image (after applying Ha), p1′, p2′, p3′, p4′, . . . :
 Thereafter, Hb is applied to the whole image, i.e., the pixel in accordance with formula (2) is moved to the pixel in accordance with formula (28).
 4. Using 2 pairs of corresponding points pa′ and pb′ in the first image, and qa′ and qb′ in the second image, the epipoles between the transformed input images and the output image are computed. If e12 denotes the epipole between the input images, e12′ denotes the epipole between the transformed first image and the output image, e22′ denotes the epipole between the transformed second image and the output image, then—in homogeneous coordinates—formulae (29) and (30) are obtained., which corresponds to formula (31) in real coordinates: where w is a ratio.
 5. For each corresponding points pi in the transformed first image and qi in the transformed second image, the location ri′ is computed in homogenous coordinates in accordance with formula (32) for the point in the new view.
 6. All points ri′, are moved to the point ri in accordance with formula (33), where e12′=[ex12′ ey12′], and d is a constant.
 Thus, when two images are used, there is effectively a separate transformation for each pixel of the stationary background, while in the first described embodiment, the background is covered by one or more transformations, each of which relates to a number of pixels.
 The disparity in the registered images depends only on the translation between the centers of the two cameras. In generating the new view, one can change the disparity to agree with the translation with respect to the location of the center of the new camera. Such a transformation is described in: M. Irani, P. Anandan, D. Weinshall, “From Reference Frames to Reference Planes: A New Framework for 3D Scene Analysis”. Proceedings: Fifth European Conference on Computer Vision, Freiburg, June 1998.
 From reading the present disclosure, other modifications will be apparent to persons skilled in the art. Such modifications may involve other features that are already known in the design, manufacture and use of virtual reality type systems and which may be used instead of or in addition to features already described herein. Although claims have been formulated in this application to particular combinations of features, it should be understood that the scope of the disclosure of the present application also includes any novel feature or novel combination of features disclosed herein either explicitly or implicitly or any generalization thereof, whether or not it mitigates any or all of the same technical problems as does the present invention. The applicants hereby give notice that new claims may be formulated to such features during the prosecution of the present application or any further application derived therefrom. In particular, although claims are formulated herein for a method, Applicant reserves the right to add apparatus and software claims in the future.
 The word “comprising”, “comprise”, or “comprises” as used herein should not be viewed as excluding additional elements. The singular article “a” or “an” as used herein should not be viewed as excluding a plurality of elements.
|Cited Patent||Filing date||Publication date||Applicant||Title|
|US2151733||May 4, 1936||Mar 28, 1939||American Box Board Co||Container|
|CH283612A *||Title not available|
|FR1392029A *||Title not available|
|FR2166276A1 *||Title not available|
|GB533718A||Title not available|
|Citing Patent||Filing date||Publication date||Applicant||Title|
|US7742062||Jun 30, 2004||Jun 22, 2010||Thomson Licensing||Automatic resetting method intended for a geometric model of a scene over a picture of the scene, implementing device and programming medium|
|US8131062||Mar 25, 2004||Mar 6, 2012||Atelier Vision Limited||Method of processing an image|
|US8270768 *||Jan 12, 2006||Sep 18, 2012||Koninklijke Philips Electronics N. V.||Depth perception|
|US20060268009 *||Jun 30, 2004||Nov 30, 2006||Lionel Oisel||Automatic resetting method intended for a geometric model of a scene over a picture of the scene, implementing device and programming medium|
|WO2004088594A1 *||Mar 25, 2004||Oct 14, 2004||Atelier Vision Ltd||Image processing|
|WO2005006252A1 *||Jun 30, 2004||Jan 20, 2005||Kijak Ewa||Automatic resetting method intended for a geometric model of a scene over a picture of the scene, implementing device and programming medium|
|U.S. Classification||345/419, 348/E13.023, 348/E13.016, 345/634, 348/E13.065, 375/E07.086, 345/427, 348/E13.008, 375/E07.076, 382/173, 348/E13.066, 348/E13.02, 348/E13.071|
|International Classification||H04N7/26, G06T7/00, G06T1/00, H04N13/00, G06T17/40|
|Cooperative Classification||H04N19/23, H04N19/20, G06T7/002, H04N13/0289, H04N2013/0092, H04N13/0246, H04N13/0221, G06T2207/10012, H04N13/0059, H04N13/026, H04N13/0011, H04N13/0014|
|European Classification||H04N13/02C, H04N13/00P1A1, H04N7/26J, G06T7/00C1, H04N7/26J6B|
|May 8, 2001||AS||Assignment|
Owner name: KONINKLIJKE PHILIPS ELECTRONICS N.V., NETHERLANDS
Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:LEE, MI-SUEN;BRODSKY, TOMAS;WEINSHALL, DAPHNA;AND OTHERS;REEL/FRAME:011800/0969
Effective date: 20010503
|Feb 4, 2009||AS||Assignment|
|May 7, 2009||FPAY||Fee payment|
Year of fee payment: 4
|Dec 20, 2011||AS||Assignment|
Owner name: FUNAI ELECTRIC CO., LTD., JAPAN
Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:IPG ELECTRONICS 503 LIMITED;REEL/FRAME:027497/0001
Effective date: 20110824
|Mar 7, 2013||FPAY||Fee payment|
Year of fee payment: 8