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
|
/************************************************************************
*
* Copyright (C) 2009-2025 IRCAD France
* Copyright (C) 2012-2019 IHU Strasbourg
*
* This file is part of Sight.
*
* Sight is free software: you can redistribute it and/or modify it under
* the terms of the GNU Lesser General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* Sight is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public
* License along with Sight. If not, see <https://www.gnu.org/licenses/>.
*
***********************************************************************/
#pragma once
#include <core/progress/has_monitors.hpp>
#include <data/boolean.hpp>
#include <data/image_series.hpp>
#include <data/integer.hpp>
#include <data/model_series.hpp>
#include <data/real.hpp>
#include <service/filter.hpp>
#include <service/notifier.hpp>
#include <vtkImageData.h>
#include <vtkPolyData.h>
namespace sight::module::filter::mesh
{
/**
* @brief This service generates meshes from an image series using VTK library.
*
* This filter is composed of three steps:
* - a mesher, by default a Marching Cubes algorithm,
* - a smoother, that allows to soften edges,
* - a decimator, that reduces the number of triangles. This is the only phase that is not multithreaded and
* unfortunately it is quite slow.
*
* The service can either mesh a single value when passed as a property, or it can mesh multiple values if a <config>
* tag is provided, allowing to specify different settings for each organ.
*
* @section Signals Signals
* - \b completed(): When the mesher succeeded.
* - \b failed(): When the mesher failed.
*
* @section XML XML Configuration
* @code{.xml}
<service type="sight::module::filter::mesh::vtk_mesher" >
<in key="image_series" uid="..."/>
<inout key="model_series" uid="..." />
<config mode="add" />
<organ name="Vessels" type="Liver" value="1" color="#ff00ff" />
<organ name="Kidney" type="" value="254" color="#ff00ffe0" />
</config>
<properties pass_band="0.0001" percent_reduction="0"/>
</service>
@endcode
*
* @subsection Input Input
* - \b image_series [sight::data::image_series] : image used to generate the output mesh
* @subsection InOut InOut
* - \b model_series [sight::data::model_series]: mesh generated from ImageSeries
* @subsection Properties Properties
* - \b value: Voxel value used to determine what should considered as a mesh, unless working in config
* - \b flying_edges: Use flying edges algorithm instead of marching cubes.
* - \b pass_band: Set the passband value for the windowed sinc filter.
* - \b boundary_smoothing: Turn on/off the smoothing of points on the boundary of the mesh.
* - \b non_manifold_smoothing: Smooth non-manifold points.
* - \b feature_smoothing: Turn on/off smoothing points along sharp interior edges.
* - \b feature_angle: Specify the feature angle for sharp edge identification (requires feature_smoothing to be on)
* - \b num_iterations: Specify the number of iterations (i.e., the degree of the polynomial approximating the windowed
* sinc function).
* - \b percent_reduction: Specify the desired reduction in the total number of polygons (e.g., if
* TargetReduction is set to 90, this filter will try to reduce the data set to 10% of its original size).
* Beware, this operation can be quite slow and proportional to this setting.
* - \b preserve_topology: Turn on/off whether to preserve the topology of the original mesh when reducing the number
* - \b quadric_reduction: Improve decimation quality but increases the time by around 15%. However, it also acts as
* wonderful smoothing tool and if used with the faster flying_edges mesher, overall, this can provide faster results
* for the same quality.
* @subsection Configuration Configuration
* - \b label : This means that voxels outside the value are converted
* to black (bit value of zero), and pixels of the given value are converted to white (a bit value of one).
* - \b mode(add,replace): Specify if the reconstruction should be added or replace the existing ones.
* - \b organ: allows to mesh configure each segmentation
* - \b name: name of the organ
* - \b value: the value to use from the mask to mesh
* - \b color (optional, hexadecimal, default=#FFFFFF): color of the mesh
* - \b split (optional, boolean): Specify if the reconstructions should be split if they are not geometrically
* connected.
* - \b selected (optional, boolean): Specify if the reconstruction is marked as selected, which will can be used
* later for instance to highlight it in the visualisation.
* - \b material (optional, string): name of the Ogre material used for rendering.
* - \b uniforms (optional, string): "|" separated list of uniform keys/values, where values can be interpreted as a
* sight::data::string_serializable i.e. "u_float=1.0|u_vec3=1.0;5.0;6.0".
*/
class vtk_mesher : public service::filter,
public sight::service::notifier,
public sight::core::progress::has_monitors
{
public:
SIGHT_DECLARE_SERVICE(vtk_mesher, sight::service::filter);
struct signals
{
static inline const signal_key_t COMPLETED = "completed";
static inline const signal_key_t FAILED = "failed";
using empty_t = sight::core::com::signal<void ()>;
};
vtk_mesher() noexcept;
~vtk_mesher() noexcept override = default;
protected:
connections_t auto_connections() const override;
/// Set add or replace mode
void configuring(const config_t& _config) override;
/// Does nothing
void starting() override;
/// Does nothing
void stopping() override;
/// Performs the meshing operation
void updating() override;
private:
vtkSmartPointer<vtkPolyData> reconstruct(vtkSmartPointer<vtkImageData> _image, int _value);
enum class mode_t : std::uint8_t
{
ADD,
REPLACE
};
mode_t m_mode {mode_t::ADD};
/// Input image mask
sight::data::ptr<sight::data::image_series, sight::data::access::in> m_image {this, "image_series"};
/// Output segmentation
sight::data::ptr<sight::data::model_series, sight::data::access::inout> m_model {this, "model_series"};
/// Expose the pass-band, since it provides the biggest control over smoothing
data::property<data::boolean> m_boundary_smoothing {this, "boundary_smoothing", false};
data::property<data::boolean> m_non_manifold_smoothing {this, "non_manifold_smoothing", true};
data::property<data::boolean> m_feature_smoothing {this, "feature_smoothing", false};
data::property<data::boolean> m_preserve_topology {this, "preserve_topology", false};
data::property<data::boolean> m_use_flying_edges {this, "flying_edges", false};
data::property<data::integer> m_value {this, "value", 1};
data::property<data::integer> m_num_iterations {this, "num_iterations", 50};
data::property<data::boolean> m_quadric_reduction {this, "quadric_reduction", true};
data::property<data::real> m_reduction {this, "percent_reduction", 0.};
data::property<sight::data::real> m_pass_band {this, "pass_band", 0.01};
data::property<data::real> m_feature_angle {this, "feature_angle", 120.};
};
} // namespace sight::module::filter::mesh
|
