GDAL
gdal_simplesurf.h
Go to the documentation of this file.
1/******************************************************************************
2 * $Id$
3 * Project: GDAL
4 * Purpose: Correlator
5 * Author: Andrew Migal, migal.drew@gmail.com
6 *
7 ******************************************************************************
8 * Copyright (c) 2012, Andrew Migal
9 *
10 * SPDX-License-Identifier: MIT
11 ****************************************************************************/
12
19#ifndef GDALSIMPLESURF_H_
20#define GDALSIMPLESURF_H_
21
22#include "gdal_priv.h"
23#include "cpl_conv.h"
24#include <list>
25
36{
37 public:
43
49
62 GDALFeaturePoint(int nX, int nY, int nScale, int nRadius, int nSign);
63 virtual ~GDALFeaturePoint();
64
67
78 double &operator[](int nIndex);
79
90 double operator[](int nIndex) const;
91
93 static const int DESC_SIZE = 64;
94
100 int GetX() const;
101
107 void SetX(int nX);
108
114 int GetY() const;
115
121 void SetY(int nY);
122
128 int GetScale() const;
129
135 void SetScale(int nScale);
136
142 int GetRadius() const;
143
149 void SetRadius(int nRadius);
150
156 int GetSign() const;
157
163 void SetSign(int nSign);
164
165 private:
166 // Coordinates of point in image
167 int nX;
168 int nY;
169 // --------------------
170 int nScale;
171 int nRadius;
172 int nSign;
173 // Descriptor array
174 double *padfDescriptor;
175};
176
187{
189
190 public:
192 virtual ~GDALIntegralImage();
193
201 void Initialize(const double **padfImg, int nHeight, int nWidth);
202
212 double GetValue(int nRow, int nCol);
213
225 double GetRectangleSum(int nRow, int nCol, int nWidth, int nHeight);
226
236 double HaarWavelet_X(int nRow, int nCol, int nSize);
237
247 double HaarWavelet_Y(int nRow, int nCol, int nSize);
248
254 int GetHeight();
255
261 int GetWidth();
262
263 private:
264 double **pMatrix = nullptr;
265 int nWidth = 0;
266 int nHeight = 0;
267};
268
279{
281
282 public:
284
293 GDALOctaveLayer(int nOctave, int nInterval);
294 virtual ~GDALOctaveLayer();
295
305 void ComputeLayer(GDALIntegralImage *poImg);
306
322 int scale;
326 int width;
334 double **detHessians;
338 int **signs;
339};
340
349{
351
352 public:
361 GDALOctaveMap(int nOctaveStart, int nOctaveEnd);
362 virtual ~GDALOctaveMap();
363
370 void ComputeMap(GDALIntegralImage *poImg);
371
391 static bool PointIsExtremum(int row, int col, GDALOctaveLayer *bot,
393 double threshold);
394
399
403 static const int INTERVALS = 4;
404
409
414};
415
429{
430 private:
435 class MatchedPointPairInfo
436 {
437 public:
438 MatchedPointPairInfo(int nInd_1, int nInd_2, double dfDist)
439 : ind_1(nInd_1), ind_2(nInd_2), euclideanDist(dfDist)
440 {
441 }
442
443 int ind_1;
444 int ind_2;
445 double euclideanDist;
446 };
447
449
450 public:
475 GDALSimpleSURF(int nOctaveStart, int nOctaveEnd);
476 virtual ~GDALSimpleSURF();
477
495 GDALRasterBand *green,
496 GDALRasterBand *blue, int nXSize,
497 int nYSize, double **padfImg,
498 int nHeight, int nWidth);
499
515 std::vector<GDALFeaturePoint> *
516 ExtractFeaturePoints(GDALIntegralImage *poImg, double dfThreshold);
517
538 static CPLErr
539 MatchFeaturePoints(std::vector<GDALFeaturePoint *> *poMatchPairs,
540 std::vector<GDALFeaturePoint> *poFirstCollect,
541 std::vector<GDALFeaturePoint> *poSecondCollect,
542 double dfThreshold);
543
544 private:
554 static double GetEuclideanDistance(const GDALFeaturePoint &firstPoint,
555 const GDALFeaturePoint &secondPoint);
556
562 static void NormalizeDistances(std::list<MatchedPointPairInfo> *poList);
563
570 static void SetDescriptor(GDALFeaturePoint *poPoint,
571 GDALIntegralImage *poImg);
572
573 private:
574 int octaveStart;
575 int octaveEnd;
576 GDALOctaveMap *poOctMap;
577};
578
579#endif /* GDALSIMPLESURF_H_ */
Class of "feature point" in raster.
Definition: gdal_simplesurf.h:36
GDALFeaturePoint()
Standard constructor.
Definition: gdal_simplesurf.cpp:22
void SetY(int nY)
Set Y coordinate of point.
Definition: gdal_simplesurf.cpp:80
GDALFeaturePoint & operator=(const GDALFeaturePoint &point)
Assignment operator.
Definition: gdal_simplesurf.cpp:43
int GetSign() const
Fetch sign of Hessian determinant of point.
Definition: gdal_simplesurf.cpp:105
void SetX(int nX)
Set X coordinate of point.
Definition: gdal_simplesurf.cpp:70
static const int DESC_SIZE
Descriptor length.
Definition: gdal_simplesurf.h:93
int GetY() const
Fetch Y-coordinate (line) of point.
Definition: gdal_simplesurf.cpp:75
int GetRadius() const
Fetch radius of point.
Definition: gdal_simplesurf.cpp:95
int GetScale() const
Fetch scale of point.
Definition: gdal_simplesurf.cpp:85
int GetX() const
Fetch X-coordinate (pixel) of point.
Definition: gdal_simplesurf.cpp:65
void SetScale(int nScale)
Set scale of point.
Definition: gdal_simplesurf.cpp:90
void SetRadius(int nRadius)
Set radius of point.
Definition: gdal_simplesurf.cpp:100
void SetSign(int nSign)
Set sign of point.
Definition: gdal_simplesurf.cpp:110
double & operator[](int nIndex)
Provide access to point's descriptor.
Definition: gdal_simplesurf.cpp:115
Integral image class (summed area table).
Definition: gdal_simplesurf.h:187
int GetHeight()
Fetch height of integral image.
Definition: gdal_octave.cpp:22
double HaarWavelet_Y(int nRow, int nCol, int nSize)
Get value of vertical Haar wavelet in specified square grid.
Definition: gdal_octave.cpp:127
double GetValue(int nRow, int nCol)
Fetch value of specified position in integral image.
Definition: gdal_octave.cpp:74
double HaarWavelet_X(int nRow, int nCol, int nSize)
Get value of horizontal Haar wavelet in specified square grid.
Definition: gdal_octave.cpp:121
void Initialize(const double **padfImg, int nHeight, int nWidth)
Compute integral image for specified array.
Definition: gdal_octave.cpp:32
double GetRectangleSum(int nRow, int nCol, int nWidth, int nHeight)
Get sum of values in specified rectangular grid.
Definition: gdal_octave.cpp:82
int GetWidth()
Fetch width of integral image.
Definition: gdal_octave.cpp:27
Class for computation and storage of Hessian values in SURF-based algorithm.
Definition: gdal_simplesurf.h:279
int height
Image height in pixels.
Definition: gdal_simplesurf.h:330
int ** signs
Hessian signs for speeded matching.
Definition: gdal_simplesurf.h:338
double ** detHessians
Hessian values for image pixels.
Definition: gdal_simplesurf.h:334
int filterSize
Length of the side of filter.
Definition: gdal_simplesurf.h:314
int radius
Length of the border.
Definition: gdal_simplesurf.h:318
int octaveNum
Octave which contains this layer (1,2,3...)
Definition: gdal_simplesurf.h:310
int width
Image width in pixels.
Definition: gdal_simplesurf.h:326
int scale
Scale for this layer.
Definition: gdal_simplesurf.h:322
void ComputeLayer(GDALIntegralImage *poImg)
Perform calculation of Hessian determinants and their signs for specified integral image.
Definition: gdal_octave.cpp:156
Class for handling octave layers in SURF-based algorithm.
Definition: gdal_simplesurf.h:349
static const int INTERVALS
Value for constructing internal octave space.
Definition: gdal_simplesurf.h:403
int octaveStart
Number of bottom octave.
Definition: gdal_simplesurf.h:408
GDALOctaveLayer *** pMap
2-dimensional array of octave layers
Definition: gdal_simplesurf.h:398
GDALOctaveMap(int nOctaveStart, int nOctaveEnd)
Create octave space.
Definition: gdal_octave.cpp:228
void ComputeMap(GDALIntegralImage *poImg)
Calculate Hessian values for octave space (for all stored octave layers) using specified integral ima...
Definition: gdal_octave.cpp:240
int octaveEnd
Number of top octave.
Definition: gdal_simplesurf.h:413
static bool PointIsExtremum(int row, int col, GDALOctaveLayer *bot, GDALOctaveLayer *mid, GDALOctaveLayer *top, double threshold)
Method makes decision that specified point in middle octave layer is maximum among all points from 3x...
Definition: gdal_octave.cpp:247
A single raster band (or channel).
Definition: gdal_priv.h:1519
Class for searching corresponding points on images.
Definition: gdal_simplesurf.h:429
std::vector< GDALFeaturePoint > * ExtractFeaturePoints(GDALIntegralImage *poImg, double dfThreshold)
Find feature points using specified integral image.
Definition: gdal_simplesurf.cpp:238
GDALSimpleSURF(int nOctaveStart, int nOctaveEnd)
Prepare class according to specified parameters.
Definition: gdal_simplesurf.cpp:148
static CPLErr MatchFeaturePoints(std::vector< GDALFeaturePoint * > *poMatchPairs, std::vector< GDALFeaturePoint > *poFirstCollect, std::vector< GDALFeaturePoint > *poSecondCollect, double dfThreshold)
Find corresponding points (equal points in two collections).
Definition: gdal_simplesurf.cpp:369
static CPLErr ConvertRGBToLuminosity(GDALRasterBand *red, GDALRasterBand *green, GDALRasterBand *blue, int nXSize, int nYSize, double **padfImg, int nHeight, int nWidth)
Convert image with RGB channels to grayscale using "luminosity" method.
Definition: gdal_simplesurf.cpp:156
Various convenience functions for CPL.
CPLErr
Error category.
Definition: cpl_error.h:37
#define CPL_DISALLOW_COPY_ASSIGN(ClassName)
Helper to remove the copy and assignment constructors so that the compiler will not generate the defa...
Definition: cpl_port.h:1030
C++ GDAL entry points.