EP2923331A1 - Verfahren zur erstellung panoramischer 3d-mosaike einer szene - Google Patents

Verfahren zur erstellung panoramischer 3d-mosaike einer szene

Info

Publication number
EP2923331A1
EP2923331A1 EP13795213.1A EP13795213A EP2923331A1 EP 2923331 A1 EP2923331 A1 EP 2923331A1 EP 13795213 A EP13795213 A EP 13795213A EP 2923331 A1 EP2923331 A1 EP 2923331A1
Authority
EP
European Patent Office
Prior art keywords
reconstruction
scene
mosaicing
panoramic
image
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
EP13795213.1A
Other languages
English (en)
French (fr)
Inventor
Raphaël HORAK
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Thales SA
Original Assignee
Thales SA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Thales SA filed Critical Thales SA
Publication of EP2923331A1 publication Critical patent/EP2923331A1/de
Withdrawn legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/60Control of cameras or camera modules
    • H04N23/698Control of cameras or camera modules for achieving an enlarged field of view, e.g. panoramic image capture
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/04Context-preserving transformations, e.g. by using an importance map
    • G06T3/047Fisheye or wide-angle transformations
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4038Image mosaicing, e.g. composing plane images from plane sub-images
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/38Registration of image sequences
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/55Depth or shape recovery from multiple images
    • G06T7/579Depth or shape recovery from multiple images from motion
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/55Depth or shape recovery from multiple images
    • G06T7/593Depth or shape recovery from multiple images from stereo images
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/40Extracting pixel data from image sensors by controlling scanning circuits, e.g. by modifying the number of pixels sampled or to be sampled
    • H04N25/41Extracting pixel data from a plurality of image sensors simultaneously picking up an image, e.g. for increasing the field of view by combining the outputs of a plurality of sensors
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00Indexing scheme for image data processing or generation, in general
    • G06T2200/32Indexing scheme for image data processing or generation, in general involving image mosaicing
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10032Satellite or aerial image; Remote sensing
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20221Image fusion; Image merging
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30181Earth observation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N13/00Stereoscopic video systems; Multi-view video systems; Details thereof
    • H04N2013/0074Stereoscopic image analysis
    • H04N2013/0088Synthesising a monoscopic image signal from stereoscopic images, e.g. synthesising a panoramic or high resolution monoscopic image

Definitions

  • the field of the invention is that of the 3D mosaicing of a scene from successive panoramic images of this scene, one or more 3D reconstructions of this scene having previously been performed.
  • the 3D reconstruction of a scene consists in obtaining, from successive 2D images of this scene taken from different points of view, a so-called 3D reconstructed image such that at each pixel of the reconstructed image, that is to say at any point where the reconstruction declares that there is a scene element, are associated the coordinates of the point of the corresponding scene, defined in a reference X, Y, Z linked to this scene.
  • the classic mosaicization called 2D mosaicing consists of successive images of a scene to project successively them on a main plane of the scene and to assemble them to make a mosaic.
  • Reference works include:
  • the finalized 3D reconstruction is not obvious, since it consists of the assembly of local 3D reconstructions (resulting from the process of stereoscopic restitution of 2 images often small field) which can be very noisy because of the limited number images that made it possible to build it, the limited field of cameras and the fact that the reconstruction plans dependent on the respective attitudes of the cameras have a geometry that is difficult to measure with precision (the relative position and geometry of the cameras used for the reconstruction 3D is often imprecise in practice because they are cameras that are 1 or 2 meters apart from each other and can vibrate relative to each other: it is even more obvious when these cameras are motorized).
  • the stereoscopic system renders poorly the plans that are almost perpendicular to one of the two cameras (this is the problem of the return of sloping roofs in stereoscopic aerial or satellite imagery).
  • a low or medium field camera moving but the 3D reconstruction is limited by the course and orientation of the camera and is not omnidirectional; moreover, the reconstruction may have holes due to uncontrolled movements of the camera or non-overlapping of the camera during its movement.
  • the algorithms used for the 3D reconstruction impose a reconstruction in a reference linked or close to the plane focal point of the camera which limits the possibilities of reconstruction (only one main plane of reconstruction and very limited reconstruction when the camera changes orientation).
  • the result of the reconstruction is also very noisy and can present many errors due to the low overlap between images, a constant reconstruction plan of the reconstructed scene (and a camera that can deviate from this plane) and the use of algorithms that exploit for 3D reconstruction only two images separated by a relatively small distance.
  • the mosaicization obtained by the plating on the ground of the successive images does not work and is not compliant when the scene is not flat and / or comprises elements 3D.
  • the 3D reconstruction is not omnidirectional and is not necessarily segmented, the measurements being obtained in the form of point clouds difficult to use automatically.
  • the mesh obtained by these active sensors has the disadvantage of being non-angularly dense (typically less than 4 points per m 2 for airborne applications at 1 km height). The technique is not suitable for the moment to be able to produce a textured image of the scene and must almost always be corrected manually.
  • All of the above solutions are unsuitable for 3D mosaicing for a large 3D scene, for example covering the scene over more than 120 ° of angular extent on all sides and can be built continuously along the entire path.
  • the 3D instantaneous mosaics obtained have deformations and are limited in angular (typically ⁇ 30 °) or spatial extent.
  • the assembly of mosaics is complex when the terrain is 3D and the final result does not conform to the geometry of the scene.
  • 3D mosaicing according to the invention is defined as a generalization of 2D mosaicing by operating this mosaic on any non-planar scene comprising a set of curved surfaces and 3D objects in relief.
  • the object of the invention is to be able to perform a 3D mosaic on any 3D scene from a panoramic sensor performing successive acquisitions of the scene along any trajectory and following different points of view, which ensure both a significant spatial and angular extent of mosaicing and conformity thereof.
  • This 3D mosaic design assumes that one or more 3D reconstructions of the scene have been established on the same scene as the panoramic sensor.
  • the proposed solution is based on the use:
  • mosaicing exploiting this information to produce a mosaic with a very large spatial and angular extent representing this scene according to all the points of view of the panoramic system during its movement along a path called mosaicing.
  • the panoramic system comprises one or more sensors whose images do not necessarily have overlapping between them, and allow to cover the entire scene to reconstruct instantly (with holes if the sensors are not overlapping) or during the movement .
  • the subject of the invention is a method of mosaicing a scene into a 3D mosaic. It is mainly characterized in that at least one 3D reconstruction of the scene has been obtained during the following steps of: acquisition of successive images by a sensor moving along a 3D reconstruction trajectory,
  • the textures recovered are those of the rectified image (s).
  • step C) includes a step of determining sectors of the current panoramic image that correspond to the selected projection surfaces. This characteristic is preferably applied when the previous case (with rectified images) can not be used. These steps are preferably repeated with each new 2D panoramic image acquisition.
  • a 3D mosaicing is thus carried out which is a generalization of the 2D mosaicing in the sense that the projection can be done on any 3D surface, which itself can be made of several flat or non-planar surfaces with discontinuities.
  • This 3D mosaicking consists of successive 2D images of a scene (taken from different points of view) and the 3D reconstruction of the scene in the previous sense, to project and assemble the different 2D images on the geometric modeling of the 3D reconstruction, allowing to restore the whole scene in the form of a textured mosaic plated on the various 3D elements of this scene. It makes it possible to reproduce in an appropriate manner an assembly of images on any scene presenting relief or 3D elements.
  • the reconstituted 3D mosaic is a 3D textured reconstruction of the scene.
  • each texture having a resolution it comprises a step for determining the resolution of the textures and the fusion of the textures of the step E) is carried out according to these resolutions.
  • the reconstruction and mosaicization trajectories can be the same trajectory.
  • the one or more 3D reconstructions have been obtained by using several reconstruction planes to benefit from the possibility of projecting the textures of the mosaic on all projection planes, in particular in very different planes of directions, and thus extend to maximum the possibility of angular and spatial extent of 3D mosaicking as indicated in the first example of 3D reconstruction method described later.
  • the invention also relates to a 3D mosaicization equipment of a scene that comprises: a panoramic system capable of forming 2D images of the scene, called 2D panoramic images, provided with localization means and,
  • a calculator comprising:
  • automatic processing means for complementary images possibly associated with or replaced by a man-machine interface for complementary images possibly associated with or replaced by a man-machine interface.
  • FIG. 1 diagrammatically represents an example of equipment for implementing the 3D reconstruction and mosaicing method according to the invention
  • FIG. 2 diagrammatically represents different steps of an example of a multi-plane 3D reconstruction method
  • FIG. 3 schematically represents various steps of the mosaicing method according to the invention
  • FIG. 4 illustrates measurement ambiguities produced by a concave object when there is only one reconstruction plane
  • FIG. 5 represents an example of sectoral decomposition of a panoramic image resulting from a panoramic system
  • FIG. 6 represents examples of rectified images of sectors of the panoramic image of FIG. 5 projected on different grinding planes
  • FIG. 7 schematically represents, for an exemplary trajectory, examples of rectification planes, lines of sight of the panoramic system Ldv1 and Ldv2 being independent of these planes,
  • FIG. 8 schematically represents an example of time evolution of grinding planes and 3D reconstruction planes according to the invention, for a given trajectory.
  • FIG. 9 represents the result of a mosaic of successive panoramic images in the case where each projected image takes into account the 3D reconstruction according to the invention (FIG. 9b) and in the case where it does not take it into account (FIG. 9a).
  • the exploitation of the transverse field is done by reconstructing the relief and the texture according to all the lateral points of view seen by the system. panoramic view that can be presented to the operator according to different reconstruction plans.
  • temporal stereoscopy with panoramic system which differs from conventional stereoscopy using two low field cameras
  • a trusted map related to the hierarchy of the quality of the information extracted from the 2D images targeting an object of the scene on very different points of view and which is directly related to the temporal exploitation of the 2D images of a moving panoramic system. Exploitation of a 3D reconstruction to treat only the visible parts of the textures in a 3D mosaic applied on any 3D surface, which can include complex objects, including concave
  • FIG. 1 which comprises:
  • a panoramic system 1 capable of forming 2D panoramic images of the scene, comprising a sensor 14 associated with an optical 11 and provided with locating means such as a GPS 12 and an inertial unit 13, and
  • a computer 2 comprising: means 21 for implementing the mosaicking method
  • the 2D images come from the panoramic system 1 traveling along a known reconstruction trajectory, which can be measured in relative image-to-image as the displacement progresses, thanks to the means location and calculator 2.
  • the system is panoramic in that it provides a 2D panoramic image.
  • it may include a lens 1 1 large fish-eye type field, or any conventional optical medium or catadioptric large field capable of providing a 2D panoramic image, or from a weaker field optics but animated more or less ample movements to capture the different portions of scenes that we want to rebuild in their entirety.
  • a 2D image covering a large field greater than 60 ° is for example obtained from a 45 ° field system 1 with a movement allowing it to cover this total field of 60 °.
  • the choice of the panoramic system 1 technology is not limited: it can be passive but it can be generalized to an active system as long as it allows to implement the multi-plane fusion step presented here.
  • This panoramic system may also comprise a set of optical sensors that are not independent of each other, covering together a determined panoramic angular coverage, for example identical from one image to the other, or maximum. All these optical sensors may not be overlapping, that is to say that the overall image obtained at a time by this set is not continuous (may include holes), the "holes" being filled when moving this set.
  • An example of a 2D panoramic image obtained with fish-eye type optics, and sectors (5 in this example) is shown in FIG.
  • This trajectory can be arbitrary and / or determined as the 3D reconstruction process unfolds.
  • the trajectory can be calculated as the panning system is moved by locating means measuring the relative displacements of position and attitude of the panoramic system in the scene such as GPS 12, inertial unit 13 or other.
  • This movement can be controlled by an operator via a human-machine interface 22 or be autonomous.
  • the images thus obtained are such that the image of at least one point of the scene is in at least 3 panoramic images respectively obtained according to different panoramic-point system directions of the scene.
  • the processing step of these 2D panoramic images respectively obtained at successive instants, by the processing unit 21 comprises the following substeps described in connection with FIG. a 3D reconstruction process preceding the actual mosaicization process.
  • Step a) Determine reconstruction plans in the scene.
  • Different reconstruction plans Cj can be chosen to establish the 3D reconstructions by highlighting different aspects of the scene, for example to cover the scene over a large spatial and angular extent, or which will allow to have a better representation of it .
  • Each reconstruction plane in the scene is determined experimentally by the operator or automatically according to the scene already restored; it can also be determined according to the trajectory of the panoramic system, typically around the average of this trajectory calculated between two shots, and according to the complexity of the scene and is independent of the line of sight of the panoramic system.
  • the selected reconstruction planes may be, for example, 3 or four planes that are tangent to a cylinder which surround the average trajectory of the system, so as to ensure a reconstruction in the different directions visible by the panoramic system.
  • the plane of the ground a plane perpendicular to the ground and tangent on one side to the cylinder surrounding the trajectory, a plane perpendicular to the ground and tangent on the other side of the cylinder, a plane parallel to the ground located at a height greater than 100m.
  • these previously defined reconstruction planes can be updated to approach or merge with the planar surfaces of the current reconstruction, extractable automatically or experimentally by an operator.
  • these previously defined reconstruction planes can be updated to approach or merge with the planar surfaces of the current reconstruction, extractable automatically or experimentally by an operator.
  • several parallel or perpendicular planes are used to restore the uniqueness and completeness of the 3D representation. This is the case for example when the scene comprises a concave object, or in the case where a single reconstruction plane provides different measurements of 3D dimensions depending on the angle at which the measurement is made, and is therefore unable to provide a single measure, as shown in Figure 4.
  • This figure illustrates the Z reconstruction ambiguity of the point (X, Y): the acquisitions at positions 1 and 2 of the trajectory reconstruct z1 on the reconstruction plane P1, but the acquisitions at the 2 and 3 of the trajectory reconstruct z2 on the same plane of projection P1.
  • a plane P3 is also chosen to find the lateral limits of the concave object.
  • This rectification consists of calculating at least one projection plane best adapted to the rectification and applying the transformation transforming any sector of each of the two 2D panoramic images on each plane.
  • Each projection plane used for the rectification can be chosen freely by the operator according to the trajectory of the panoramic system, from an infinite choice of positions and orientations all parallel to the trajectory of the panoramic system; the plane or each of them is independent of the evolution of the line of sight of the panoramic system (which can rotate on itself during its movement along its trajectory), contrary to the classic stereoscopy where the plane of chosen correction depends on the evolution of the line of sight and the choice of grinding plans is very limited.
  • FIG. 7 An example of rectification planes referenced R1, R2, R3 is shown in FIG. 7; they are parallel to the trajectory. Is also indicated direction of the LdV (LdV1, LdV2) of the panoramic two-point trajectory sensor, which illustrates that the choice of these plans is independent of the LdV.
  • FIG. 8 is a top view of a scene comprising 3D objects.
  • position pairs (1, 2, 3) of the panoramic sensor corresponding to 3 pairs of panoramic images acquired during this step b); each pair of positions is associated with two grinding planes (R1 1, R21 for the pair 1, R12, R22 for the pair 2 and R13, R23 for the pair 3).
  • Three reconstruction planes P1, P2, P3 were chosen in the example of this figure.
  • the rectification plane which is linked to the reconstruction planes chosen in a), is chosen so as to be the closest in the geometric sense of the reconstruction plane determined in step a). It is also chosen so as to guarantee a minimum of pixellic resolution and must not be too far in the angular direction of the associated reconstruction plane (typical angular difference of less than 30 °).
  • This transformation is an example of transformation in the case of a rectilinear panoramic optics (fisheye type); it does not understand the distortion parameters of the panoramic system that can be calculated and offset elsewhere.
  • the transformation can easily be generalized and adapted to any panoramic system having its own optical formula.
  • a sector of the panoramic image corresponds to an equivalent portion projected on a rectification plane.
  • the sectoral decomposition of the panoramic image depends on the rectification plans chosen and the influence of the projection on these plans.
  • FIG. 6 Examples of rectified images are shown in FIG. 6.
  • the first results from the projection of the sector 1 of the image of FIG. 5 on a vertical rectification plane
  • the second results from the projection of the sector 2 of the image of the Figure 5 on another vertical rectification plane
  • the third results from the projection of the sector 3 on a vertical rectification plane different from the first two
  • the fourth results from the projection of the sector 5 on a horizontal rectification plane.
  • This projection is reiterated in the same rectification plane P, for a sector of another 2D panoramic image obtained at time t + ⁇ to obtain another corrected image, ⁇ being predetermined experimentally or determined so that the displacement From the system between t and t + ⁇ realizes a stereo base large enough to be compatible with the desired accuracy for 3D reconstruction.
  • De must at least be equal to 5m, which corresponds to a minimum angle difference of 6 ° between 2 acquisitions by the panoramic system.
  • the use of a panoramic system makes it possible to increase the reconstruction accuracy by increasing the distance De and the angular separation between two acquisitions, beyond what a weak or medium-field sensor can do for the same spatial coverage of reconstruction.
  • 3D The 3D stereoscopic base for 3D reconstruction may be larger than that of a conventional stereoscopic process due to the use of a panoramic field (and the longer presence of objects in this field), and this allows the This process has a greater ultimate reconstruction accuracy, which is also increased by the fusion of the measurements offered by the process.
  • the actual acquisition of the panoramic system can be faster while keeping the displacement between the pairs of imagery rectified 2D used to reconstruct the 3D of the scene.
  • the method then consists in taking a first pair of 2D panoramic images from a displacement De, doing an intermediate 3D reconstruction with this pair, then taking another pair of 2D images always from a displacement of to the next acquisition to remake an intermediate 3D reconstruction as long as the points of the scene concerned by these different pairs of images remain in the field of the panoramic system.
  • the intermediate 3D reconstruction in a 3D reference linked to the P is obtained by matching point-to-point the two rectified images in P, and with the knowledge of the movement of the panoramic system.
  • This mapping is a dense process that matches, as far as possible, each of the points in a 2D image of the stereoscopic pair with a point in the other image.
  • Step d): transform this intermediate 3D reconstruction into a fixed 3D reference ( absolute) including the reconstruction plane determined in step a), called 3D reconstruction mark. A transformed intermediate 3D reconstruction is thus obtained.
  • steps b) to d) are successively linked in this order; these iterations can also be performed in parallel (several steps b) are performed in parallel with several rectification planes P, determined in parallel, etc.).
  • these steps b) to d) are repeated as long as at least one point of the reconstructed scene remains in the field of view of the panoramic system.
  • This 3D reconstruction method makes it possible to find the most appropriate dense 3D mesh for representing the scene, such that at each point of this mesh are associated the coordinates of the corresponding point in a frame ⁇ , ⁇ , ⁇ linked to the scene.
  • Step g) possible: repeat steps b) to f) for each reconstruction plane chosen in a), with the same panoramic images but with different sectors, so as to obtain as many 3D reconstructions as selected reconstruction plans.
  • These 3D reconstructions or the intermediate 3D reconstructions obtained during these iterations are advantageously spatially merged to update the final 3D reconstruction (s), and thus increase the accuracy and robustness of these reconstructions.
  • the spatial fusion of the 3D reconstructions constructed according to different planes takes into account the reconstruction precision of the different elements of each reconstruction which is not the same according to the different planes and which can be predicted mathematically. This spatial fusion is obtained by exploiting several rectification planes corresponding to the different sectors of each image used.
  • the set of steps a) to g) are also preferably repeated at least once with new pairs of panoramic images, for example with temporally offset intermediate images of the previous ones, or with other sectors of the panoramic images already considered. This allows for a continuous process of updating the final 3D reconstructions. These new pairs of panoramic images can come from each panoramic image acquisition but not necessarily.
  • the exploitation of the redundancies and the quality of the 2D rectified images allows the process to produce a confidence map reflecting the quality of the final reconstruction.
  • This confidence map is constructed pixel by pixel for each 3D reconstruction, considering the number of times each pixel has been built and if the conditions of this construction were good, these being for example defined according to a threshold of match quality determined experimentally or mathematically.
  • the rectification plans of the various pairs of images chosen to carry out the reconstruction are those which are the most Geometrically close to the construction plan and only one sector is chosen to perform the correction on these different rectification plans.
  • the final 3D reconstruction is performed only on a single reconstruction plane (the one initially chosen). It follows that only the faces and surfaces that are not too far angularly ( ⁇ 30 °) from the construction plan will be restored, the others may be either poorly restored or may have gaps.
  • the sectors that will be extracted for mosaicing will be those corresponding to the rectification plans used with this panoramic image in the reconstruction. They will be well adapted to the mosaic of 3D reconstruction, but will be less suitable and sometimes poorly adapted to mosaicing badly reconstructed parts, for example perpendicular or near parts to be perpendicular to the reconstruction plane.
  • a third example of 3D reconstruction that can be used is performed by repeated stereoscopy: the rectification plane depends entirely on the movement of the line of sight of the panoramic sensor and is made in a plane that is as close as possible to the focal planes of the two images. used for successive reconstructions.
  • the reconstruction plan is often confused or close to the first rectification plan used for reconstruction (most often we take the focal plane of the first image or a parallel plan as a reference).
  • the projection used for the rectification is not direct but passes through a first projection along the focal plane, followed by a second transformation linked to that which connects the focal plane of each image to the plane of rectification.
  • This second transformation may present a significant rotation because of the movements of the line of sight of the sensor, thus inducing a ground image filling very little the plane of rectification and may even induce the impossibility of matching with the second ground image, or worst matching errors.
  • the second problem is the significant distance possible in the angular direction of the rectification plane with respect to the reconstruction plane inducing in the final reconstruction, if it is made in a frame linked to the scene, significant noise reconstruction. It follows path portions without 3D reconstruction or with 3D reconstruction erroneous. At best, the method decreases compared to the two preceding methods the number of images that can be used to correctly reconstruct a point and therefore increases the reconstruction noise.
  • the 3D reconstruction presents in this case the following characteristics: correctly restored but strongly noisy for the geometrically closest surfaces of the reconstruction plane, which may present significant errors in places or gaps, and restored in a very partial way or poorly rendered for faces or surfaces almost perpendicular to the reconstruction plane.
  • a fourth example of a 3D reconstruction is one in which the 3D reconstruction has been carried out according to another trajectory and from other points of view than those used for the 3D mosaicization envisaged. This is referred to as external 3D reconstruction as opposed to previous examples of internal 3D reconstruction.
  • the sensor used for the 3D reconstruction is not necessarily panoramic.
  • the 3D reconstruction and 3D mosaicking paths may be spatially different and / or temporally different.
  • the rectification plans used by this 3D reconstruction can no longer be used to match the texture of the panoramic image with the 3D reconstruction, and it becomes necessary to use in this case only the point of view to delimit the image. texture to project.
  • the projection on the surfaces assuming that the choice of surfaces is adapted to the texture to be projected during the mosaicing, will be tainted by the relative point of view estimation error between the sensor and the 3D model.
  • mosaic 2D images of the scene acquired by a panoramic sensor moving along a mosaic path
  • This mosaic generally has several 2D textured planes present in the 3D scene or that approximate it, but can also be on a 3D surface.
  • projection surfaces are the surfaces on which the mosaic is built; they can be freely chosen by the operator or can be determined automatically from the 3D reconstruction. As indicated above, some of these surfaces may be left (curves) or even be a 3D surface of which modeling is known.
  • mosaicization planes In the case of a panoramic system viewing a highly 3D scene having different faces, several mosaicization planes (or surfaces) may (and have interest to) be chosen.
  • highly 3D scene we mean a scene containing many 3D elements producing significant disparities between two successive acquisitions, as is the case for example of a drone flying over a low-flying urban environment.
  • the mosaicing method exploits the fact that the surfaces or projection planes have different orientations for projecting the textures of the images on each of the projection surfaces or planes.
  • the texture is a set of intensities of pixels on an image region.
  • the exploitation of the 3D reconstruction makes it possible to keep only the visible parts of the projected images. This makes it possible to avoid projecting portions of images that belong to other portions of the scene and that would be hidden in the projection onto a mosaicization plane.
  • the multi-plane (or multi-surface) mosaicing process is reiterated preferably with each new 2D image acquisition performed by the panoramic system, and the new mosaic is merged with the old one (obtained at t-1). to update this one.
  • the 3D reconstruction of the scene and the projections of the textures on it are calculated at each acquisition of the so-called initial 2D image, an acquisition being separated from the previous one by a time interval ⁇ defined above.
  • the 3D reconstruction of the scene and the projections of the textures on it are calculated at each image acquisition of the panoramic system and at a high frequency from previous images previously stored.
  • the reconstruction and mosaicization trajectories are then the same.
  • the intermediate images comprised between two successive images of ⁇ successive for the 3D reconstruction are stored so that they can also be used for 3D reconstruction in the manner of a FIFO acronym for the English expression "First In First Out "(each new image acquired is compared with the first image stored to establish a new 3D reconstruction instance, this first image is then deleted from the list and the last added to the updated list).
  • the intermediate images can also be used to facilitate the correspondence between the first and last image, or to fill in "holes" in the 3D model.
  • a mosaic is then obtained at the end of the following steps A) to E) described in connection with Figure 3, for each new 2D image acquired by the panoramic system.
  • 3D reconstruction and mosaicing are carried out successively after each acquisition; this supposes that a new 3D reconstruction of reference has just been made following (one or) 3D reconstructions already made.
  • 3D reconstruction and mosaicing are performed in parallel after each acquisition; this assumes that the mosaicization is performed while a new 3D reconstruction is still in progress, in which case the reference 3D reconstruction is the one performed to one of the previous acquisitions of 2D images, or a 3D reconstruction performed previously.
  • This step consists of choosing the 3D projection planes or surfaces on which the mosaic is built.
  • These 3D projection planes or surfaces can be freely chosen by the operator at a given moment of the mosaicing or computed automatically from the current 3D reconstruction (or reference) of the scene according to predetermined criteria (for example, planes parallel to the reconstituted soil surface or main planes extracted from the 3D reconstruction).
  • Non-planar 3D projection surfaces may also be used if the scene is suitable and if the operator sees an interest therein; this allows for example to represent objects of the scene or a backdrop that have particular geometric shapes, but this does not detract from the conformity that can be obtained by multiple projections that would be exclusively flat.
  • the choice of projection surfaces can be made for each panoramic image in all possible directions of the panoramic image, these surfaces being chosen preferentially from those of the reconstructed 3D model and as close as possible to the geometric sense of the available correction plans containing rectified images from the current panoramic image.
  • the possibilities in the direction of the projection surfaces that can be used by mosaicing are directly related to the number of reconstruction planes used by the reconstruction method.
  • the projection is not performed immediately but the projection parameters are stored in memory for use by step E).
  • step E Select with the help of the (or) reconstruction (s) 3D, among the previous textures the parts of these which are visible on each projection surface and preferably from the latter those which are exploitable (that is to say having a sufficient resolution) for the projection of step E, from:
  • step C the corrected 2D image from the current image if the first case of step C applies
  • step C the second case of step C applies.
  • the newly calculated 3D reconstruction automatically calculates the hidden or weakly resolved parts (and conversely the exploitable visible parts) in the projection plane that would result from masks present in the scene. This is like selecting the textures to keep in the mosaic dataset. This calculation is accurate because the 3D reconstruction was first performed in the reference linked to the panoramic system.
  • One of the peculiarities of 3D mosaicing according to the invention is to take advantage of the calculation of the hidden parts to eliminate the masks generated by the scene on the different projection planes (surfaces) and to consider only the visible parts on these planes (surfaces).
  • This makes it possible to temporally mosaic only portions of scenes that are always visible and thus avoid deformations due to the projection of parts of the scene not belonging to the projection plane (default present in a conventional mosaic which projects from all the image on the projection plane without being able to take into account the masking by the elements of the scene evolving as the sensor moves
  • Figure 9 shows the case of a mosaic with taking into account that the parts visible and that resulting from the projection of the textures without delimitation of the always visible parts: it results in the second case of a progressive deformation of the textures due to the projection of textures of which a part should have been masked).
  • the method of limiting the textures projected to the visible parts in the model depends on the accuracy of the 3D model and the precision of rendering point of view between the current image and this 3D model.
  • the imperfections of restitution of surfaces perpendicular or close to being perpendicular to the reconstruction plane in the case of a single-plane process prevent correct calculation of the projection of textures on these surfaces and to delimit correctly the masking of these imperfections on d other surfaces.
  • a direct projection of the rectified image used or the selected sector in the panoramic panoramic image is also preferred. current, in order to maximize the surface of the projected textures on the projection plane (s).
  • the textures selected in the previous step are projected onto the projection planes (or surfaces), and temporally merged with the current mosaic to form a new mosaic.
  • the mosaic presents the mosaic to the operator according to different planes or more generally according to different 3D surfaces of presentation, by projection of the 3D reconstruction (s) textured (s) on these plans of presentation.
  • These presentation plans are freely chosen by the operator and serve only to present the results of the mosaic according to different perspectives chosen by it.
  • the mosaic can be presented to the presentation plan operator by presentation plan, or according to the plans representing the unfolding of the curved surfaces on which the textures have been projected (in the case of a projection on a cylinder for example).
  • One can of course also present the textured 3D result directly in the form of virtual 3D using a suitable software.
  • the projection result provides a still consistent image, which is not necessarily the case, as has been explained, with a conventional mosaic method.
  • FIG. 9a shows the mosaicization results obtained without applying the method according to the invention, in particular without applying a selection criterion on the textures: deformations can be seen.
  • FIG. 9b shows the mosaicization results obtained by using, according to the invention, a 3D reconstruction for projecting the unmasked textures onto the projection surface bound to the ground: there is no deformation.
  • This method of 3D reconstruction and simultaneous omnidirectional mosaicking is not limited to an optical panoramic system.
  • the textures measured over a large directional field can be very well exploited by means other than optics, for example by an active means of lidar or sonar type; the process could then also exploit the distances given by the instruments.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Signal Processing (AREA)
  • Multimedia (AREA)
  • Stereoscopic And Panoramic Photography (AREA)
  • Studio Devices (AREA)
  • Processing Or Creating Images (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)
  • Image Generation (AREA)
EP13795213.1A 2012-11-21 2013-11-20 Verfahren zur erstellung panoramischer 3d-mosaike einer szene Withdrawn EP2923331A1 (de)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
FR1203126A FR2998400B1 (fr) 2012-11-21 2012-11-21 Procede de reconstruction 3d et de mosaiquage 3d panoramique d'une scene
FR1302583A FR2998401B1 (fr) 2012-11-21 2013-11-08 Procede de reconstruction 3d et de mosaiquage 3d panoramique d'une scene
PCT/EP2013/074278 WO2014079880A1 (fr) 2012-11-21 2013-11-20 Procede de mosaiquage 3d panoramique d'une scene

Publications (1)

Publication Number Publication Date
EP2923331A1 true EP2923331A1 (de) 2015-09-30

Family

ID=48224824

Family Applications (3)

Application Number Title Priority Date Filing Date
EP16197081.9A Active EP3144881B1 (de) 2012-11-21 2013-11-20 Verfahren zur 3d-panorama-mosaikbildung einer szene
EP13792927.9A Active EP2923330B1 (de) 2012-11-21 2013-11-20 Verfahren zur 3d-rekonstruktion und panoramischen 3d-mosaikerstellung einer szene
EP13795213.1A Withdrawn EP2923331A1 (de) 2012-11-21 2013-11-20 Verfahren zur erstellung panoramischer 3d-mosaike einer szene

Family Applications Before (2)

Application Number Title Priority Date Filing Date
EP16197081.9A Active EP3144881B1 (de) 2012-11-21 2013-11-20 Verfahren zur 3d-panorama-mosaikbildung einer szene
EP13792927.9A Active EP2923330B1 (de) 2012-11-21 2013-11-20 Verfahren zur 3d-rekonstruktion und panoramischen 3d-mosaikerstellung einer szene

Country Status (7)

Country Link
US (3) US9729789B2 (de)
EP (3) EP3144881B1 (de)
JP (2) JP6407878B2 (de)
DK (1) DK3144881T3 (de)
FR (2) FR2998400B1 (de)
IL (3) IL238939B (de)
WO (2) WO2014079880A1 (de)

Families Citing this family (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105096251A (zh) * 2015-07-18 2015-11-25 哈尔滨工业大学 一种利用超分辨率重建技术提高拼接图像分辨率的方法
US10291845B2 (en) * 2015-08-17 2019-05-14 Nokia Technologies Oy Method, apparatus, and computer program product for personalized depth of field omnidirectional video
TWI555378B (zh) * 2015-10-28 2016-10-21 輿圖行動股份有限公司 一種全景魚眼相機影像校正、合成與景深重建方法與其系統
US9823664B2 (en) 2016-02-25 2017-11-21 A.M.T.S., Llc Unmanned aircraft for positioning an instrument for inspection purposes and methods of inspecting a target surface
JP6835536B2 (ja) * 2016-03-09 2021-02-24 株式会社リコー 画像処理方法、表示装置および点検システム
CN107203965B (zh) * 2016-03-18 2020-03-10 中国科学院宁波材料技术与工程研究所 一种基于多通道图像融合的全景图像拼接方法
US10313575B1 (en) * 2016-11-14 2019-06-04 Talon Aerolytics, Inc. Drone-based inspection of terrestrial assets and corresponding methods, systems, and apparatuses
KR101837403B1 (ko) 2016-12-13 2018-04-19 국방과학연구소 고속 영상 모자이킹 방법 및 장치
US10192111B2 (en) 2017-03-10 2019-01-29 At&T Intellectual Property I, L.P. Structure from motion for drone videos
US10861359B2 (en) * 2017-05-16 2020-12-08 Texas Instruments Incorporated Surround-view with seamless transition to 3D view system and method
NL2019516B1 (en) * 2017-09-08 2019-03-19 Fugro N V System and method for determination of a spatial property of a submerged object in a 3D-space
CN108053409B (zh) * 2017-12-11 2022-05-13 中南大学 一种遥感影像分割基准库自动构建方法及其系统
TWI643498B (zh) * 2018-01-02 2018-12-01 元智大學 以單一影像計算鏡頭視角的方法與攝像裝置
US11430145B2 (en) 2018-06-17 2022-08-30 Foresight Ai Inc. Identification of local motions in point cloud data
CN110968715B (zh) * 2018-09-28 2022-08-05 哈尔滨理工大学 基于非均匀投影的三维模型检索
CN109064410B (zh) * 2018-10-24 2023-03-14 清华大学深圳研究生院 一种基于超像素的光场图像拼接方法
CN110503621B (zh) * 2019-08-23 2022-07-01 上海圭目机器人有限公司 基于定位数据的图像全局拼接方法
CN111445544A (zh) * 2020-02-24 2020-07-24 东莞职业技术学院 一种基于路线实况检测的地图绘制方法
CN111882936A (zh) * 2020-07-18 2020-11-03 纽智医疗科技(苏州)有限公司 一种基于力反馈与虚拟现实的手术预演及教学系统
CN113971628B (zh) * 2020-07-24 2025-05-06 株式会社理光 图像匹配方法、装置和计算机可读存储介质
WO2022081717A1 (en) * 2020-10-13 2022-04-21 Flyreel, Inc. Generating measurements of physical structures and environments through automated analysis of sensor data
CN113344795B (zh) * 2021-08-05 2021-10-29 常州铭赛机器人科技股份有限公司 一种基于先验信息的快速图像拼接方法
CN114630037B (zh) * 2022-03-22 2024-02-02 北京有竹居网络技术有限公司 拍摄装置及其控制方法、装置、可读存储介质和电子设备
CN115291225B (zh) * 2022-07-27 2025-05-16 四川中烟工业有限责任公司 一种基于超声的场景全息成像的方法
CN117132716B (zh) * 2023-10-24 2024-01-26 北京瑞祺皓迪技术股份有限公司 一种动力环境监控方法、装置、电子设备及存储介质
CN118840486A (zh) * 2024-07-05 2024-10-25 三星电子(中国)研发中心 场景的三维重建方法、装置、电子设备和存储介质

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6858826B2 (en) * 1996-10-25 2005-02-22 Waveworx Inc. Method and apparatus for scanning three-dimensional objects
ATE280978T1 (de) * 1997-01-30 2004-11-15 Yissum Res Dev Co Mosaikbildverarbeitungssystem
JP2003141562A (ja) * 2001-10-29 2003-05-16 Sony Corp 非平面画像の画像処理装置及び画像処理方法、記憶媒体、並びにコンピュータ・プログラム
JP4255272B2 (ja) * 2002-12-06 2009-04-15 学校法人早稲田大学 画像生成装置
JP4355535B2 (ja) * 2003-08-07 2009-11-04 株式会社岩根研究所 360度画像変換処理装置
JP4511147B2 (ja) * 2003-10-02 2010-07-28 株式会社岩根研究所 三次元形状生成装置
JP2006350852A (ja) * 2005-06-17 2006-12-28 Nabura:Kk 画像生成システム
US8818076B2 (en) * 2005-09-01 2014-08-26 Victor Shenkar System and method for cost-effective, high-fidelity 3D-modeling of large-scale urban environments
JP4889351B2 (ja) * 2006-04-06 2012-03-07 株式会社トプコン 画像処理装置及びその処理方法
KR101188588B1 (ko) * 2008-03-27 2012-10-08 주식회사 만도 모노큘러 모션 스테레오 기반의 주차 공간 검출 장치 및방법
EP2686827A1 (de) * 2011-03-18 2014-01-22 C3 Technologies AB 3d-strassen

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
None *
See also references of WO2014079880A1 *

Also Published As

Publication number Publication date
IL259954B (en) 2019-05-30
WO2014079880A1 (fr) 2014-05-30
WO2014079879A3 (fr) 2014-08-07
JP2017097871A (ja) 2017-06-01
EP2923330B1 (de) 2016-12-21
IL259954A (en) 2018-07-31
IL238939A0 (en) 2015-07-30
JP6407878B2 (ja) 2018-10-17
IL238939B (en) 2018-07-31
FR2998400B1 (fr) 2016-01-15
US9998660B2 (en) 2018-06-12
EP2923330A2 (de) 2015-09-30
US20150304576A1 (en) 2015-10-22
JP6483075B2 (ja) 2019-03-13
FR2998400A1 (fr) 2014-05-23
FR2998401A1 (fr) 2014-05-23
FR2998401B1 (fr) 2017-03-10
IL238938A0 (en) 2015-07-30
US9729789B2 (en) 2017-08-08
EP3144881B1 (de) 2018-08-15
WO2014079879A2 (fr) 2014-05-30
EP3144881A1 (de) 2017-03-22
JP2016502720A (ja) 2016-01-28
US20170310892A1 (en) 2017-10-26
DK3144881T3 (en) 2018-12-03
US20150304558A1 (en) 2015-10-22

Similar Documents

Publication Publication Date Title
EP2923330B1 (de) Verfahren zur 3d-rekonstruktion und panoramischen 3d-mosaikerstellung einer szene
Remondino et al. State of the art in high density image matching
EP2715662B1 (de) Verfahren zum lokalisieren einer kamera und 3d-rekonstruktion in einer teilweise bekannten umgebung
FR2853121A1 (fr) Dispositif de surveillance des alentours d'un vehicule
EP2984625B1 (de) Verfahren zur messung der winkelabweichung eines mobilen elements bezüglich einer referenzrichtung
WO2006064051A1 (fr) Procede de traitement d'images mettant en oeuvre le georeferencement automatique d'images issues d'un couple d'images pris dans le meme plan focal
FR2953313A1 (fr) Systeme optronique et procede d'elaboration d'images en trois dimensions dedies a l'identification
WO2012175888A1 (fr) Génération de données de carte
EP3377854B1 (de) Verfahren zur erfassung von bildern einer szene, von einem sensor an bord eines beweglichen trägers, mit servosteuerung von dessen sichtlinie
EP2428934A1 (de) Verfahren zur Schätzung der Bewegung eines Trägers im Verhältnis zu einer Umgebung, und Berechnungsvorrichtung für ein Navigationssystem
EP1756771B1 (de) Fein-stereoskopischer bildvergleich und dediziertes instrument mit einem geringen stereoskopischen koeffizient
FR2952743A3 (fr) Procede d'estimation du mouvement d'un instrument d'observation a defilement survolant un corps celeste
EP2432660A1 (de) Verfahren und vorrichtung zur erweiterung eines sichtbereichs
KR102225321B1 (ko) 복수 영상 센서로부터 취득한 영상 정보와 위치 정보 간 연계를 통한 도로 공간 정보 구축을 위한 시스템 및 방법
EP4526843B1 (de) Verfahren und vorrichtung zur bestimmung der absoluten pose einer kamera an bord eines luft- oder raumfahrzeuges
Barazzetti et al. 3D Modeling and Rendering with a Tesla Model 3 Highland
WO2015075085A1 (fr) Procédé d'estimation du mouvement d'un objet
EP1958157A1 (de) Verfahren zur einfügung stereoskopischer bilder in eine korrespondenz

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20150522

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

DAX Request for extension of the european patent (deleted)
17Q First examination report despatched

Effective date: 20170127

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN

18D Application deemed to be withdrawn

Effective date: 20170808