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
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
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
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
|
/** @file sift.h
** @brief SIFT (@ref sift)
** @author Andrea Vedaldi
**/
/*
Copyright (C) 2007-12 Andrea Vedaldi and Brian Fulkerson.
All rights reserved.
This file is part of the VLFeat library and is made available under
the terms of the BSD license (see the COPYING file).
*/
#ifndef VL_SIFT_H
#define VL_SIFT_H
#include <stdio.h>
#include "generic.h"
/** @brief SIFT filter pixel type */
typedef float vl_sift_pix ;
/** ------------------------------------------------------------------
** @brief SIFT filter keypoint
**
** This structure represent a keypoint as extracted by the SIFT
** filter ::VlSiftFilt.
**/
typedef struct _VlSiftKeypoint
{
int o ; /**< o coordinate (octave). */
int ix ; /**< Integer unnormalized x coordinate. */
int iy ; /**< Integer unnormalized y coordinate. */
int is ; /**< Integer s coordinate. */
float x ; /**< x coordinate. */
float y ; /**< y coordinate. */
float s ; /**< s coordinate. */
float sigma ; /**< scale. */
} VlSiftKeypoint ;
/** ------------------------------------------------------------------
** @brief SIFT filter
**
** This filter implements the SIFT detector and descriptor.
**/
typedef struct _VlSiftFilt
{
double sigman ; /**< nominal image smoothing. */
double sigma0 ; /**< smoothing of pyramid base. */
double sigmak ; /**< k-smoothing */
double dsigma0 ; /**< delta-smoothing. */
int width ; /**< image width. */
int height ; /**< image height. */
int O ; /**< number of octaves. */
int S ; /**< number of levels per octave. */
int o_min ; /**< minimum octave index. */
int s_min ; /**< minimum level index. */
int s_max ; /**< maximum level index. */
int o_cur ; /**< current octave. */
vl_sift_pix *temp ; /**< temporary pixel buffer. */
vl_sift_pix *octave ; /**< current GSS data. */
vl_sift_pix *dog ; /**< current DoG data. */
int octave_width ; /**< current octave width. */
int octave_height ; /**< current octave height. */
vl_sift_pix *gaussFilter ; /**< current Gaussian filter */
double gaussFilterSigma ; /**< current Gaussian filter std */
vl_size gaussFilterWidth ; /**< current Gaussian filter width */
VlSiftKeypoint* keys ;/**< detected keypoints. */
int nkeys ; /**< number of detected keypoints. */
int keys_res ; /**< size of the keys buffer. */
double peak_thresh ; /**< peak threshold. */
double edge_thresh ; /**< edge threshold. */
double norm_thresh ; /**< norm threshold. */
double magnif ; /**< magnification factor. */
double windowSize ; /**< size of Gaussian window (in spatial bins) */
vl_sift_pix *grad ; /**< GSS gradient data. */
int grad_o ; /**< GSS gradient data octave. */
} VlSiftFilt ;
/** @name Create and destroy
** @{
**/
VL_EXPORT
VlSiftFilt* vl_sift_new (int width, int height,
int noctaves, int nlevels,
int o_min) ;
VL_EXPORT
void vl_sift_delete (VlSiftFilt *f) ;
/** @} */
/** @name Process data
** @{
**/
VL_EXPORT
int vl_sift_process_first_octave (VlSiftFilt *f,
vl_sift_pix const *im) ;
VL_EXPORT
int vl_sift_process_next_octave (VlSiftFilt *f) ;
VL_EXPORT
void vl_sift_detect (VlSiftFilt *f) ;
VL_EXPORT
int vl_sift_calc_keypoint_orientations (VlSiftFilt *f,
double angles [4],
VlSiftKeypoint const*k);
VL_EXPORT
void vl_sift_calc_keypoint_descriptor (VlSiftFilt *f,
vl_sift_pix *descr,
VlSiftKeypoint const* k,
double angle) ;
VL_EXPORT
void vl_sift_calc_raw_descriptor (VlSiftFilt const *f,
vl_sift_pix const* image,
vl_sift_pix *descr,
int widht, int height,
double x, double y,
double s, double angle0) ;
VL_EXPORT
void vl_sift_keypoint_init (VlSiftFilt const *f,
VlSiftKeypoint *k,
double x,
double y,
double sigma) ;
/** @} */
/** @name Retrieve data and parameters
** @{
**/
VL_INLINE int vl_sift_get_octave_index (VlSiftFilt const *f) ;
VL_INLINE int vl_sift_get_noctaves (VlSiftFilt const *f) ;
VL_INLINE int vl_sift_get_octave_first (VlSiftFilt const *f) ;
VL_INLINE int vl_sift_get_octave_width (VlSiftFilt const *f) ;
VL_INLINE int vl_sift_get_octave_height (VlSiftFilt const *f) ;
VL_INLINE int vl_sift_get_nlevels (VlSiftFilt const *f) ;
VL_INLINE int vl_sift_get_nkeypoints (VlSiftFilt const *f) ;
VL_INLINE double vl_sift_get_peak_thresh (VlSiftFilt const *f) ;
VL_INLINE double vl_sift_get_edge_thresh (VlSiftFilt const *f) ;
VL_INLINE double vl_sift_get_norm_thresh (VlSiftFilt const *f) ;
VL_INLINE double vl_sift_get_magnif (VlSiftFilt const *f) ;
VL_INLINE double vl_sift_get_window_size (VlSiftFilt const *f) ;
VL_INLINE vl_sift_pix *vl_sift_get_octave (VlSiftFilt const *f, int s) ;
VL_INLINE VlSiftKeypoint const *vl_sift_get_keypoints (VlSiftFilt const *f) ;
/** @} */
/** @name Set parameters
** @{
**/
VL_INLINE void vl_sift_set_peak_thresh (VlSiftFilt *f, double t) ;
VL_INLINE void vl_sift_set_edge_thresh (VlSiftFilt *f, double t) ;
VL_INLINE void vl_sift_set_norm_thresh (VlSiftFilt *f, double t) ;
VL_INLINE void vl_sift_set_magnif (VlSiftFilt *f, double m) ;
VL_INLINE void vl_sift_set_window_size (VlSiftFilt *f, double m) ;
/** @} */
/* -------------------------------------------------------------------
* Inline functions implementation
* ---------------------------------------------------------------- */
/** ------------------------------------------------------------------
** @brief Get current octave index.
** @param f SIFT filter.
** @return index of the current octave.
**/
VL_INLINE int
vl_sift_get_octave_index (VlSiftFilt const *f)
{
return f-> o_cur ;
}
/** ------------------------------------------------------------------
** @brief Get number of octaves.
** @param f SIFT filter.
** @return number of octaves.
**/
VL_INLINE int
vl_sift_get_noctaves (VlSiftFilt const *f)
{
return f-> O ;
}
/**-------------------------------------------------------------------
** @brief Get first octave.
** @param f SIFT filter.
** @return index of the first octave.
**/
VL_INLINE int
vl_sift_get_octave_first (VlSiftFilt const *f)
{
return f-> o_min ;
}
/** ------------------------------------------------------------------
** @brief Get current octave width
** @param f SIFT filter.
** @return current octave width.
**/
VL_INLINE int
vl_sift_get_octave_width (VlSiftFilt const *f)
{
return f-> octave_width ;
}
/** ------------------------------------------------------------------
** @brief Get current octave height
** @param f SIFT filter.
** @return current octave height.
**/
VL_INLINE int
vl_sift_get_octave_height (VlSiftFilt const *f)
{
return f-> octave_height ;
}
/** ------------------------------------------------------------------
** @brief Get current octave data
** @param f SIFT filter.
** @param s level index.
**
** The level index @a s ranges in the interval <tt>s_min = -1</tt>
** and <tt> s_max = S + 2</tt>, where @c S is the number of levels
** per octave.
**
** @return pointer to the octave data for level @a s.
**/
VL_INLINE vl_sift_pix *
vl_sift_get_octave (VlSiftFilt const *f, int s)
{
int w = vl_sift_get_octave_width (f) ;
int h = vl_sift_get_octave_height (f) ;
return f->octave + w * h * (s - f->s_min) ;
}
/** ------------------------------------------------------------------
** @brief Get number of levels per octave
** @param f SIFT filter.
** @return number of leves per octave.
**/
VL_INLINE int
vl_sift_get_nlevels (VlSiftFilt const *f)
{
return f-> S ;
}
/** ------------------------------------------------------------------
** @brief Get number of keypoints.
** @param f SIFT filter.
** @return number of keypoints.
**/
VL_INLINE int
vl_sift_get_nkeypoints (VlSiftFilt const *f)
{
return f-> nkeys ;
}
/** ------------------------------------------------------------------
** @brief Get keypoints.
** @param f SIFT filter.
** @return pointer to the keypoints list.
**/
VL_INLINE VlSiftKeypoint const *
vl_sift_get_keypoints (VlSiftFilt const *f)
{
return f-> keys ;
}
/** ------------------------------------------------------------------
** @brief Get peaks treashold
** @param f SIFT filter.
** @return threshold ;
**/
VL_INLINE double
vl_sift_get_peak_thresh (VlSiftFilt const *f)
{
return f -> peak_thresh ;
}
/** ------------------------------------------------------------------
** @brief Get edges threshold
** @param f SIFT filter.
** @return threshold.
**/
VL_INLINE double
vl_sift_get_edge_thresh (VlSiftFilt const *f)
{
return f -> edge_thresh ;
}
/** ------------------------------------------------------------------
** @brief Get norm threshold
** @param f SIFT filter.
** @return threshold.
**/
VL_INLINE double
vl_sift_get_norm_thresh (VlSiftFilt const *f)
{
return f -> norm_thresh ;
}
/** ------------------------------------------------------------------
** @brief Get the magnification factor
** @param f SIFT filter.
** @return magnification factor.
**/
VL_INLINE double
vl_sift_get_magnif (VlSiftFilt const *f)
{
return f -> magnif ;
}
/** ------------------------------------------------------------------
** @brief Get the Gaussian window size.
** @param f SIFT filter.
** @return standard deviation of the Gaussian window (in spatial bin units).
**/
VL_INLINE double
vl_sift_get_window_size (VlSiftFilt const *f)
{
return f -> windowSize ;
}
/** ------------------------------------------------------------------
** @brief Set peaks threshold
** @param f SIFT filter.
** @param t threshold.
**/
VL_INLINE void
vl_sift_set_peak_thresh (VlSiftFilt *f, double t)
{
f -> peak_thresh = t ;
}
/** ------------------------------------------------------------------
** @brief Set edges threshold
** @param f SIFT filter.
** @param t threshold.
**/
VL_INLINE void
vl_sift_set_edge_thresh (VlSiftFilt *f, double t)
{
f -> edge_thresh = t ;
}
/** ------------------------------------------------------------------
** @brief Set norm threshold
** @param f SIFT filter.
** @param t threshold.
**/
VL_INLINE void
vl_sift_set_norm_thresh (VlSiftFilt *f, double t)
{
f -> norm_thresh = t ;
}
/** ------------------------------------------------------------------
** @brief Set the magnification factor
** @param f SIFT filter.
** @param m magnification factor.
**/
VL_INLINE void
vl_sift_set_magnif (VlSiftFilt *f, double m)
{
f -> magnif = m ;
}
/** ------------------------------------------------------------------
** @brief Set the Gaussian window size
** @param f SIFT filter.
** @param x Gaussian window size (in units of spatial bin).
**
** This is the parameter @f$ \hat \sigma_{\text{win}} @f$ of
** the standard SIFT descriptor @ref sift-tech-descriptor-std.
**/
VL_INLINE void
vl_sift_set_window_size (VlSiftFilt *f, double x)
{
f -> windowSize = x ;
}
/* VL_SIFT_H */
#endif
|