knn_matching.cpp 6.31 KB
Newer Older
biagio montesano's avatar
biagio montesano committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41
/*M///////////////////////////////////////////////////////////////////////////////////////
 //
 //  IMPORTANT: READ BEFORE DOWNLOADING, COPYING, INSTALLING OR USING.
 //
 //  By downloading, copying, installing or using the software you agree to this license.
 //  If you do not agree to this license, do not download, install,
 //  copy or use the software.
 //
 //
 //                           License Agreement
 //                For Open Source Computer Vision Library
 //
 // Copyright (C) 2014, Biagio Montesano, all rights reserved.
 // Third party copyrights are property of their respective owners.
 //
 // Redistribution and use in source and binary forms, with or without modification,
 // are permitted provided that the following conditions are met:
 //
 //   * Redistribution's of source code must retain the above copyright notice,
 //     this list of conditions and the following disclaimer.
 //
 //   * Redistribution's in binary form must reproduce the above copyright notice,
 //     this list of conditions and the following disclaimer in the documentation
 //     and/or other materials provided with the distribution.
 //
 //   * The name of the copyright holders may not be used to endorse or promote products
 //     derived from this software without specific prior written permission.
 //
 // This software is provided by the copyright holders and contributors "as is" and
 // any express or implied warranties, including, but not limited to, the implied
 // warranties of merchantability and fitness for a particular purpose are disclaimed.
 // In no event shall the Intel Corporation or contributors be liable for any direct,
 // indirect, incidental, special, exemplary, or consequential damages
 // (including, but not limited to, procurement of substitute goods or services;
 // loss of use, data, or profits; or business interruption) however caused
 // and on any theory of liability, whether in contract, strict liability,
 // or tort (including negligence or otherwise) arising in any way out of
 // the use of this software, even if advised of the possibility of such damage.
 //
 //M*/

42
#include <iostream>
43
#include <opencv2/opencv_modules.hpp>
44

45 46 47
#ifdef HAVE_OPENCV_FEATURES2D

#include <opencv2/line_descriptor.hpp>
48
#include <opencv2/core/utility.hpp>
49 50 51 52 53
#include <opencv2/imgproc.hpp>
#include <opencv2/features2d.hpp>
#include <opencv2/highgui.hpp>

#include <vector>
54
#include <time.h>
55 56

using namespace cv;
57
using namespace cv::line_descriptor;
58 59

static const char* keys =
60 61
{ "{@image_path1 | | Image path 1 }"
    "{@image_path2 | | Image path 2 }" };
62 63 64

static void help()
{
65 66
  std::cout << "\nThis example shows the functionalities of descriptors matching\n" << "Please, run this sample using a command in the form\n"
            << "./example_line_descriptor_matching <path_to_input_image 1>" << "<path_to_input_image 2>" << std::endl;
67 68 69

}

70 71
uchar invertSingleBits( uchar dividend_char, int numBits );

72
/* invert numBits bits in input char */
73
uchar invertSingleBits( uchar dividend_char, int numBits )
74
{
75 76 77
  std::vector<int> bin_vector;
  long dividend;
  long bin_num;
78

79 80
  /* convert input char to a long */
  dividend = (long) dividend_char;
81

82 83 84
  /*if a 0 has been obtained, just generate a 8-bit long vector of zeros */
  if( dividend == 0 )
    bin_vector = std::vector<int>( 8, 0 );
85

86 87 88 89
  /* else, apply classic decimal to binary conversion */
  else
  {
    while ( dividend >= 1 )
90
    {
91 92 93
      bin_num = dividend % 2;
      dividend /= 2;
      bin_vector.push_back( bin_num );
94
    }
95
  }
96

97 98 99 100 101 102
  /* ensure that binary vector always has length 8 */
  if( bin_vector.size() < 8 )
  {
    std::vector<int> zeros( 8 - bin_vector.size(), 0 );
    bin_vector.insert( bin_vector.end(), zeros.begin(), zeros.end() );
  }
103

104 105 106 107 108
  /* invert numBits bits */
  for ( int index = 0; index < numBits; index++ )
  {
    if( bin_vector[index] == 0 )
      bin_vector[index] = 1;
109

110 111 112
    else
      bin_vector[index] = 0;
  }
113

114
  /* reconvert to decimal */
115
  uchar result = 0;
116
  for ( int i = (int) bin_vector.size() - 1; i >= 0; i-- )
Bleach666's avatar
Bleach666 committed
117
    result += (uchar) ( bin_vector[i] * (1 << i) );
118

119
  return result;
120 121 122 123
}

int main( int argc, char** argv )
{
124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167
  /* get parameters from comand line */
  CommandLineParser parser( argc, argv, keys );
  String image_path1 = parser.get<String>( 0 );
  String image_path2 = parser.get<String>( 1 );

  if( image_path1.empty() || image_path2.empty() )
  {
    help();
    return -1;
  }

  /* load image */
  cv::Mat imageMat1 = imread( image_path1, 1 );
  cv::Mat imageMat2 = imread( image_path2, 1 );

  if( imageMat1.data == NULL || imageMat2.data == NULL )
  {
    std::cout << "Error, images could not be loaded. Please, check their paths" << std::endl;
  }

  /* create binary masks */
  cv::Mat mask1 = Mat::ones( imageMat1.size(), CV_8UC1 );
  cv::Mat mask2 = Mat::ones( imageMat2.size(), CV_8UC1 );

  /* create a pointer to a BinaryDescriptor object with default parameters */
  Ptr<BinaryDescriptor> bd = BinaryDescriptor::createBinaryDescriptor();

  /* compute lines */
  std::vector<KeyLine> keylines1, keylines2;
  bd->detect( imageMat1, keylines1, mask1 );
  bd->detect( imageMat2, keylines2, mask2 );

  /* compute descriptors */
  cv::Mat descr1, descr2;
  bd->compute( imageMat1, keylines1, descr1 );
  bd->compute( imageMat2, keylines2, descr2 );

  /* create a BinaryDescriptorMatcher object */
  Ptr<BinaryDescriptorMatcher> bdm = BinaryDescriptorMatcher::createBinaryDescriptorMatcher();

  /* make a copy of descr2 mat */
  Mat descr2Copy = descr1.clone();

  /* randomly change some bits in original descriptors */
168
  srand( (unsigned int) time( NULL ) );
169 170 171 172 173 174 175 176 177 178 179

  for ( int j = 0; j < descr1.rows; j++ )
  {
    /* select a random column */
    int randCol = rand() % 32;

    /* get correspondent data */
    uchar u = descr1.at<uchar>( j, randCol );

    /* change bits */
    for ( int k = 1; k <= 5; k++ )
180
    {
181 182
      /* copy current row to train matrix */
      descr2Copy.push_back( descr1.row( j ) );
183

184 185
      /* invert k bits */
      uchar uc = invertSingleBits( u, k );
186

187 188
      /* update current row in train matrix */
      descr2Copy.at<uchar>( descr2Copy.rows - 1, randCol ) = uc;
189
    }
190
  }
191

192 193
  /* prepare a structure to host matches */
  std::vector<std::vector<DMatch> > matches;
194

195 196
  /* require knn match */
  bdm->knnMatch( descr1, descr2, matches, 6 );
197 198 199

}

200 201 202 203 204 205 206
#else

int main()
{
    std::cerr << "OpenCV was built without features2d module" << std::endl;
    return 0;
}
207

208
#endif // HAVE_OPENCV_FEATURES2D