<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Finding a Strong Key Point Correspondences in Large-Scale Images and Depth Maps</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Vorontsov</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>rgii Rosokh</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>NPP "SPECPOZHTEHNIKA" Ltd</institution>
          ,
          <addr-line>Kharkiv 61017</addr-line>
          ,
          <country country="UA">Ukraine</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>National Technical University "Kharkiv Polytechnic Institute"</institution>
          ,
          <addr-line>Kharkiv 61002</addr-line>
          ,
          <country country="UA">Ukraine</country>
        </aff>
      </contrib-group>
      <abstract>
        <p>Depth map fusion is an essential task in multi-view stereo vision and structure from motion methods, but there are some issues of the process of depth maps matching of large-scale outdoor scenes, which are acquired by unmanned aerial vehicles. First, they can have low depth resolution due to long distance from a scene to a camera, in the second, as a consequence of above problem and some camera defects, they can have some noise. In our work we propose an approach of matching the depth maps based on a local shape and a color similarity correspondences. The proposed idea is to nd similar regions by calculating local key points in small patches of two images by ORB feature detector. Then we nd corresponding regions assuming spatial proximity of movement vectors of key points, this is done by the search of the main vector direction and length in feature space. We test our approach with the help of video sequences acquired through a consumer multicopter.</p>
      </abstract>
      <kwd-group>
        <kwd>Depth Map Fusion Multi-view Stereo Vision Structure From Motion Unmanned Aerial Vehicles Key Points ORB Feature Detector Spatial Similarity Nearest Neighbors Search</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>At present there is an increase in quantity of stereo vision method applications
in the production, robotics, in systems of virtual and augmented reality,
geoinformation systems, computer games and simulating. This brings about the need
for creation and development e cient methods and algorithms for solving the
problems of stereo vision. The present paper represent a process of nding the
e ective ways to solve the problem of comparing large-scale depth maps obtained
by unmanned aerial vehicles. At the same time, the obtained depth maps can
have a low depth resolution due to the large distance from the camera to the
scene and, as a result, such depth maps can have noise.</p>
      <p>
        There are several approaches for solving the problem of nding strong
correspondences of images and depth maps which aimed to increase accuracy of 3d
model reconstruction. Among them are various variations of the iterative
closest points (ICP) [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ], robust point matching [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ], methods based on probability
theory, correlation and Gaussian mixture models [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ], coherent point drift [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ],
method of sorting the correspondence space [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ]. In consequence of simpleness of
realization and low computational complexity, various modi cations of the ICP
algorithm have found wide application [
        <xref ref-type="bibr" rid="ref6">6</xref>
        ]. A detailed comparative review and
ways to increase the accuracy of algorithms for ICP matching can be found in
[
        <xref ref-type="bibr" rid="ref7">7</xref>
        ]. In [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ], the problem of sampling point's assortment for increase the accuracy
of ICP algorithms was considered. As the ICP algorithms are the greedy
optimization techniques, the initial matching has great in uence on registration of
the point sets. Therefore we need methods to nd more precise initial
approximation of surface matching to increase probability of convergence to the global
minimum. The main contribution of presented work is the algorithm of nding
of initial correspondences between two surfaces which leads to decreasing of time
complexity of the ICP algorithms and increasing of ability to converge to global
minimum.
2
      </p>
    </sec>
    <sec id="sec-2">
      <title>Method Overview</title>
      <p>Given a set of n consecutive frames in F = f1; f2; :::; fn and a set of n 1 depth
maps M = m1; m2; :::; mn 1, acquired from pairs of adjacent images in F . It is
necessary to nd such pairs of points p; q 2 R3 that for each point pj 2 mi and
the corresponding point qk 2 mi + 1; i = 1::n 1; j; k = 1::jmj, j j { number of
pixels in the depth maps will be: qk = Apj + t, where A is the 3 3 rotation
matrix, t is the 3 1 transfer vector. Alongside strict correspondences between
points of neighboring depth maps are unknown in advance. Thus, the task of
nding the transformation A and t can be divided into two sub-tasks: nding
the indices j; k for the corresponding pairs of points p; q on the depth maps and
comparison of the obtained point sets.</p>
      <p>
        In our work, we proceed from the assumption that adjacent frames have
a small di erences of the direction of movement of points, which consists of
transference and rotation combination, while camera is signi cantly distant from
the scene, we assume a rectilinear motion. Based on these assumptions, it can
be expected that the movement of points in close sections of images will have
almost similar values of the motion direction. At the same time, the initial point
correspondences are found in the original RGB image space using the robust
search of feature descriptions, such as scale-invariant feature transform (SIFT)
[
        <xref ref-type="bibr" rid="ref9">9</xref>
        ] or oriented FAST and rotated BRIEF (ORB) [
        <xref ref-type="bibr" rid="ref10">10</xref>
        ].
      </p>
      <p>
        We propose the following method for detecting strong correspondences
between points on large-scale depth maps:
1. We take set of images F as input. Each image has size W H.
2. We set the parameters for dividing images into local areas, which include
the width W and the height H of the image local part.
3. At each algorithm iteration, we take the images fi and fi+1 from F and the
corresponding to them depth maps mi from M . Thus, we form the structure
data = fi; fi+1; mi.
4. We split fi and fi+1 into t rectangular sections (patches) pj ; j = 1::t, Pi =
p1; p2; :::; pt, where t = bW=wc bH=hc, bxc - integer part of x.
5. The next step is de ning arrays of key points Ki and Ki+1 for each segment
pi+1 by the ORB algorithm, Ki = kj ; j = 1::KN - number of key points.
6. Initial mapping of points from Ki and Ki+1 is conducted. As a method of
mapping can be used any method of search approximate nearest neighbors
algorithm. We used the implementation from the FLANN library [
        <xref ref-type="bibr" rid="ref11">11</xref>
        ].
However, at this stage, there are many false-positive pairs, which have strongly
in uence on the further depth map matching process.
7. To remove incorrectly recognized correspondences of points, we propose a
novel approach based on the dividing the parameter space of the key points
into a regular grid and spatial indexing of the resulting grid for nearest points
search.
8. The obtained exact matches are transferred to the depth maps and they are
compared by the ICP method.
      </p>
      <p>Let us demonstrate our approach to determine the exact matches of key
points. For each pair of key points we set the vector vj = kpj+1 kpj , where
kpj { coordinates of kj in the patch p. Each vector vj has an angle and length
l and indicates the direction of key point movement between frames.</p>
      <p>Assumption 1. In small areas of the image, correctly mapped points will give
close values of a and l, and false-positive pairs of points will give arbitrary
corresponding values.</p>
      <p>Based on this assumption, it is possible to construct a classi er that will
nd closely spaced points in the parameter space l. To speed up the nding
process of close elements in l-space, we discretize it by splitting into a regular
grid with the number of cells A along the O axis and the number of cells L
along the Ol axis, thus the size of the grid cell is:
subject to ; l 0</p>
      <p>Then we build accumulator Acc array with size (A
cell in it with zeros. For each vj we discretize it and l:</p>
      <sec id="sec-2-1">
        <title>L) and initialize each</title>
        <p>size =
lsize =
max</p>
        <p>A
max l</p>
        <p>L
j
jdiscr = b A c
lj
ljdiscr = b L c</p>
      </sec>
      <sec id="sec-2-2">
        <title>And then:</title>
        <p>Acc[ jdiscr ][ljdiscr ] = Acc[ jdiscr ][ljdiscr ] + 1
Summary of our approach is provided on the Fig. 1.
(1)
(2)
(3)
(4)
(5)
Input images and depth maps</p>
        <p>Key points detection
(ORB or SIFT descriptors)
Initial matching
with FLANN lib
Movement vectors of
key points
90o</p>
        <p>Discretization step
(t=8)
length
cell with maximal
number of vectors
45o 8</p>
        <p>Hash-table7
building 6
0o
length
315o
5</p>
        <p>Due to Assumption 1 the cell with maximal value corresponds to the matched
key point pairs and we mark such as a "good" matches. To accelerate
computations we use a hash-table H to store indices of key points, the key value of the
elements of H we calculate as follows:
hash =
jdiscr</p>
        <p>L + ljdiscr
(6)
Therefore we can nd strong correspondences of key points and their consecutive
directions.
3</p>
      </sec>
    </sec>
    <sec id="sec-3">
      <title>Experiments</title>
      <p>We tested our algorithm on depth maps obtained from a video taken by a UAV
with a large, relatively to the size of the seen, depth.</p>
      <p>Image data and algorithm parameters:
We have frames of video with 8 Mpx resolution, so size of our all images and
depth maps is (W H) = (4096 2160) pixels each.</p>
      <p>Subdivision parameters: w = W=8, h = H=8.</p>
      <p>Number of key points for ORB detector KN = 100.</p>
      <p>Number of discretized grid cells along axis: A = L = 20.</p>
      <p>Result of depth map recti cation by our algorithm is demonstrated in Fig. 2.</p>
    </sec>
    <sec id="sec-4">
      <title>Conclusions</title>
      <p>In our paper we demonstrate an approach for estimation of strong
correspondences between the key points in two sequential video frames. The proposed
algorithm can eliminate false-positive occurrences from an array of key points.
This is done by nding of major direction of key point movement in parameter
space.</p>
      <p>The limitation of our approach is that it can delete positive correspondences
between the key points, thus, the future work is aimed to improving the quality
of correct pairs matching. Another disadvantage is the ability of our algorithm
to match key points only in the case of camera or scene movement.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Bouaziz</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Tagliasacchi</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Pauly</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          :
          <article-title>Sparse Iterative Closest Point</article-title>
          .
          <source>Computer Graphics Forum</source>
          .
          <volume>32</volume>
          ,
          <issue>113</issue>
          {
          <fpage>123</fpage>
          (
          <year>2013</year>
          ). https://doi.org/10.1111/cgf.12178
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Gold</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Rangarajan</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Lu</surname>
            ,
            <given-names>C.-P.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Suguna</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Mjolsness</surname>
          </string-name>
          , E.:
          <article-title>New algorithms for 2d and 3d point matching: pose estimation and correspondence</article-title>
          .
          <source>Pattern Recognition</source>
          .
          <volume>38</volume>
          (
          <issue>8</issue>
          ), pp.
          <volume>1019</volume>
          {
          <issue>1031</issue>
          (
          <year>1998</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>Evangelidis</surname>
            ,
            <given-names>G.D.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Horaud</surname>
          </string-name>
          , R.:
          <article-title>Joint Alignment of Multiple Point Sets with Batch and Incremental Expectation-Maximization</article-title>
          .
          <source>IEEE Transactions on Pattern Analysis and Machine Intelligence</source>
          .
          <volume>40</volume>
          ,
          <issue>1397</issue>
          {
          <fpage>1410</fpage>
          (
          <year>2018</year>
          ). https://doi.org/10.1109/TPAMI.
          <year>2017</year>
          .2717829
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Myronenko</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Song</surname>
            ,
            <given-names>X.</given-names>
          </string-name>
          :
          <string-name>
            <surname>Point-Set</surname>
            <given-names>Registration</given-names>
          </string-name>
          :
          <article-title>Coherent Point Drift</article-title>
          .
          <source>IEEE Transactions on Pattern Analysis and Machine Intelligence</source>
          .
          <volume>32</volume>
          , pp.
          <volume>2262</volume>
          {
          <issue>2275</issue>
          (
          <year>2010</year>
          ). https://doi.org/10.1109/TPAMI.
          <year>2010</year>
          .46
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5.
          <string-name>
            <surname>Assalih</surname>
          </string-name>
          , H.:
          <article-title>3D reconstruction and motion estimation using forward looking sonar</article-title>
          .
          <source>Doctoral dissertation</source>
          , Heriot-Watt University. 208 p. (
          <year>2013</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <surname>Park</surname>
          </string-name>
          , S.-Y.,
          <string-name>
            <surname>Subbarao</surname>
            ,
            <given-names>M.:</given-names>
          </string-name>
          <article-title>An accurate and fast point-to-plane registration technique</article-title>
          .
          <source>Pattern Recognition Letters</source>
          .
          <volume>24</volume>
          ,
          <issue>2967</issue>
          {
          <fpage>2976</fpage>
          (
          <year>2003</year>
          ). https://doi.org/10.1016/S0167-
          <volume>8655</volume>
          (
          <issue>03</issue>
          )
          <fpage>00157</fpage>
          -
          <lpage>0</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7.
          <string-name>
            <surname>Rusinkiewicz</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Levoy</surname>
            ,
            <given-names>M.:</given-names>
          </string-name>
          <article-title>E cient variants of the ICP algorithm</article-title>
          . In: Proceedings Third International Conference on 3-
          <string-name>
            <given-names>D</given-names>
            <surname>Digital</surname>
          </string-name>
          Imaging and Modeling. pp.
          <volume>145</volume>
          {
          <fpage>152</fpage>
          . IEEE Comput. Soc, Quebec City, Que.,
          <string-name>
            <surname>Canada</surname>
          </string-name>
          (
          <year>2001</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8.
          <string-name>
            <surname>Gelfand</surname>
            ,
            <given-names>N.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ikemoto</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Rusinkiewicz</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Levoy</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          :
          <article-title>Geometrically stable sampling for the ICP algorithm</article-title>
          . In: Fourth International Conference on 3-
          <string-name>
            <given-names>D</given-names>
            <surname>Digital</surname>
          </string-name>
          Imaging and Modeling,
          <year>2003</year>
          . 3DIM
          <year>2003</year>
          . Proceedings. pp.
          <volume>260</volume>
          {
          <fpage>267</fpage>
          . IEEE, Ban , Alberta, Canada (
          <year>2003</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          9.
          <string-name>
            <surname>Lowe</surname>
            ,
            <given-names>D.G.</given-names>
          </string-name>
          :
          <article-title>Object recognition from local scale-invariant features</article-title>
          .
          <source>In: Proceedings of the Seventh IEEE International Conference on Computer Vision</source>
          . pp.
          <volume>1150</volume>
          {1157 vol.
          <volume>2</volume>
          . IEEE, Kerkyra, Greece (
          <year>1999</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          10.
          <string-name>
            <surname>Rublee</surname>
            ,
            <given-names>E.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Rabaud</surname>
            ,
            <given-names>V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Konolige</surname>
            ,
            <given-names>K.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Bradski</surname>
          </string-name>
          , G.:
          <article-title>ORB: An e cient alternative to SIFT or SURF</article-title>
          .
          <source>In: 2011 International Conference on Computer Vision</source>
          . pp.
          <volume>2564</volume>
          {
          <fpage>2571</fpage>
          . IEEE, Barcelona, Spain (
          <year>2011</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          11.
          <string-name>
            <surname>Muja</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Lowe</surname>
            ,
            <given-names>D.G.</given-names>
          </string-name>
          :
          <article-title>Fast Approximate Nearest Neighbors with Automatic Algorithm Con guration</article-title>
          .
          <source>In: International Conference on Computer Vision Theory and Application VISSAPP'09)</source>
          . pp.
          <volume>331</volume>
          {
          <fpage>340</fpage>
          . INSTICC Press (
          <year>2009</year>
          ).
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>