Geogram  Version 1.9.1
A programming library of geometric algorithms
LFS.h
Go to the documentation of this file.
1 /*
2  * Copyright (c) 2000-2022 Inria
3  * All rights reserved.
4  *
5  * Redistribution and use in source and binary forms, with or without
6  * modification, are permitted provided that the following conditions are met:
7  *
8  * * Redistributions of source code must retain the above copyright notice,
9  * this list of conditions and the following disclaimer.
10  * * Redistributions in binary form must reproduce the above copyright notice,
11  * this list of conditions and the following disclaimer in the documentation
12  * and/or other materials provided with the distribution.
13  * * Neither the name of the ALICE Project-Team nor the names of its
14  * contributors may be used to endorse or promote products derived from this
15  * software without specific prior written permission.
16  *
17  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
18  * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
19  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
20  * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
21  * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
22  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
23  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
24  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
25  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
26  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
27  * POSSIBILITY OF SUCH DAMAGE.
28  *
29  * Contact: Bruno Levy
30  *
31  * https://www.inria.fr/fr/bruno-levy
32  *
33  * Inria,
34  * Domaine de Voluceau,
35  * 78150 Le Chesnay - Rocquencourt
36  * FRANCE
37  *
38  */
39 
40 #ifndef GEOGRAM_DELAUNAY_LFS
41 #define GEOGRAM_DELAUNAY_LFS
42 
43 #include <geogram/basic/common.h>
45 
51 namespace GEO {
52 
61  class GEOGRAM_API LocalFeatureSize {
62  public:
69  LocalFeatureSize(index_t nb_pts, const double* pts) {
70  sliver_angle_threshold_ = 0.01;
71  init(nb_pts, pts);
72  }
73 
79  double squared_lfs(const double* p) const {
80  index_t v = spatial_search_->nearest_vertex(p);
81  const double* q = spatial_search_->vertex_ptr(v);
82  return
83  geo_sqr(p[0] - q[0]) +
84  geo_sqr(p[1] - q[1]) +
85  geo_sqr(p[2] - q[2]);
86  }
87 
92  index_t nb_poles() const {
93  return poles_.size()/3;
94  }
95 
103  const double* pole(index_t i) const {
104  geo_debug_assert(i < nb_poles());
105  return &poles_[3*i];
106  }
107 
108  protected:
116  void init(index_t nb_pts, const double* pts);
117 
118  private:
119  double sliver_angle_threshold_;
120  vector<double> poles_;
121  Delaunay_var spatial_search_;
122  };
123 }
124 
125 #endif
#define geo_debug_assert(x)
Verifies that a condition is met.
Definition: assert.h:196
Computes an approximation of lfs (local feature size).
Definition: LFS.h:61
LocalFeatureSize(index_t nb_pts, const double *pts)
Initializes lfs computation.
Definition: LFS.h:69
const double * pole(index_t i) const
Gets a reference to a pole.
Definition: LFS.h:103
double squared_lfs(const double *p) const
Computes the squared local feature size at a query point.
Definition: LFS.h:79
index_t nb_poles() const
Gets the number of poles.
Definition: LFS.h:92
void init(index_t nb_pts, const double *pts)
Constructs the internal representation used to compute the local feature size.
Abstract interface for Delaunay.
Common include file, providing basic definitions. Should be included before anything else by all head...
Global Vorpaline namespace.
Definition: basic.h:55
T geo_sqr(T x)
Gets the square value of a value.
Definition: numeric.h:301
geo_index_t index_t
The type for storing and manipulating indices.
Definition: numeric.h:329