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
|
// Copyright (c) 2016,2021 GeometryFactory Sarl (France).
// All rights reserved.
//
// This file is part of CGAL (www.cgal.org).
//
// $URL: https://github.com/CGAL/cgal/blob/v6.1.1/Mesh_3/include/CGAL/Mesh_3/Image_plus_weights_to_labeled_function_wrapper.h $
// $Id: include/CGAL/Mesh_3/Image_plus_weights_to_labeled_function_wrapper.h 08b27d3db14 $
// SPDX-License-Identifier: GPL-3.0-or-later OR LicenseRef-Commercial
//
//
// Author(s) : Laurent Rineau
//
#ifndef CGAL_MESH_3_IMAGE_PLUS_WEIGHTS_TO_LABELED_FUNCTION_WRAPPER_H
#define CGAL_MESH_3_IMAGE_PLUS_WEIGHTS_TO_LABELED_FUNCTION_WRAPPER_H
#include <CGAL/license/Mesh_3.h>
#include <CGAL/function_objects.h>
#include <CGAL/Image_3.h>
namespace CGAL {
namespace ImageIO {
template <typename Image_word_type,
typename Weights_type>
class Weighted_indicator_factory
{
const Image_3& image;
const Image_3& weights;
public:
Weighted_indicator_factory(const Image_3& image,
const Image_3& weights)
: image(image)
, weights(weights)
{}
class Indicator
{
const Weighted_indicator_factory& f;
const Image_word_type label;
public:
Indicator(const Weighted_indicator_factory& f,
const Image_word_type label) : f(f), label(label) {}
double operator()(const Image_word_type& x) const
{
const std::ptrdiff_t offset = &x - (Image_word_type *)f.image.data();
const Weights_type w = (std::max)(
Weights_type(128), ((Weights_type *)f.weights.data())[offset]);
return (x == label) ? w : (255 - w);
}
}; // end nested class Indicator
Indicator indicator(const Image_word_type label) const
{
return Indicator(*this, label);
}
}; // end template class ImageIO::Weighted_indicator_factory
} // end namespace CGAL::ImageIO
namespace Mesh_3 {
/**
* @class Image_plus_weights_to_labeled_function_wrapper
*
* Wraps a pair of images into a labeled function which takes his values into
* N. Uses weighted trilinear interpolation.
*/
template<typename Image_word_type = unsigned char,
typename Interpolation_type = Image_word_type,
typename Weights_type = unsigned char,
typename Return_type = int>
class Image_plus_weights_to_labeled_function_wrapper
{
public:
typedef std::function<Return_type(Interpolation_type)>
Image_values_to_labels;
// Types
typedef Return_type return_type;
typedef Image_word_type word_type;
typedef CGAL::Image_3 Image_;
/// Constructor
Image_plus_weights_to_labeled_function_wrapper
(const Image_& image,
const Image_& weights_image,
Image_values_to_labels transform = Identity<Return_type>(),
const Interpolation_type value_outside = Interpolation_type())
: r_im_(image)
, r_weights_im_(weights_image)
, transform(transform)
, value_outside(value_outside)
, indicator_factory(image, weights_image)
{
}
// Default copy constructor and assignment operator are ok
/// Destructor
~Image_plus_weights_to_labeled_function_wrapper() {}
/**
* Returns an int corresponding to the label at point `p`.
* @param p the input point
* @return the label at point `p`
*/
template <typename Point_3>
return_type operator()(const Point_3& p) const
{
return static_cast<return_type>(transform(
r_im_.template labellized_trilinear_interpolation<Image_word_type>(
CGAL::to_double(p.x() - r_im_.image()->tx),
CGAL::to_double(p.y() - r_im_.image()->ty),
CGAL::to_double(p.z() - r_im_.image()->tz),
value_outside,
indicator_factory)));
}
private:
/// Labeled image to wrap
const Image_& r_im_;
const Image_& r_weights_im_;
const Image_values_to_labels transform;
const Interpolation_type value_outside;
CGAL::ImageIO::Weighted_indicator_factory<Image_word_type,
Weights_type> indicator_factory;
}; // end class Image_plus_weights_to_labeled_function_wrapper
} // end namespace Mesh_3
} // end namespace CGAL
#endif // CGAL_MESH_3_IMAGE_PLUS_WEIGHTS_TO_LABELED_FUNCTION_WRAPPER_H
|
