Developer Documentation
OVMPropertyVisualizerVectorFieldDifferenceT_impl.hh
1/*===========================================================================*\
2* *
3* OpenFlipper *
4 * Copyright (c) 2001-2015, RWTH-Aachen University *
5 * Department of Computer Graphics and Multimedia *
6 * All rights reserved. *
7 * www.openflipper.org *
8 * *
9 *---------------------------------------------------------------------------*
10 * This file is part of OpenFlipper. *
11 *---------------------------------------------------------------------------*
12 * *
13 * Redistribution and use in source and binary forms, with or without *
14 * modification, are permitted provided that the following conditions *
15 * are met: *
16 * *
17 * 1. Redistributions of source code must retain the above copyright notice, *
18 * this list of conditions and the following disclaimer. *
19 * *
20 * 2. Redistributions in binary form must reproduce the above copyright *
21 * notice, this list of conditions and the following disclaimer in the *
22 * documentation and/or other materials provided with the distribution. *
23 * *
24 * 3. Neither the name of the copyright holder nor the names of its *
25 * contributors may be used to endorse or promote products derived from *
26 * this software without specific prior written permission. *
27 * *
28 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS *
29 * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED *
30 * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A *
31 * PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER *
32 * OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, *
33 * EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, *
34 * PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR *
35 * PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF *
36 * LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING *
37 * NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS *
38 * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. *
39* *
40\*===========================================================================*/
41
42
43
44
45#define OVM_PROPERTY_VISZUALIZTER_VECTOR_FIELD_DIFFERENCE_CC
46
47#include "OVMPropertyVisualizerVectorFieldDifference.hh"
48
49template <typename MeshT>
50OVMPropertyVisualizerVectorFieldDifference<MeshT>::OVMPropertyVisualizerVectorFieldDifference(MeshT* _mesh, int objectID, PropertyInfo _propertyInfo1, PropertyInfo _propertyInfo2)
51 : OVMPropertyVisualizer<MeshT>(_mesh, objectID, _propertyInfo1),
52 propertyInfo2(_propertyInfo2)
53{
54 if (PropertyVisualizer::widget) delete PropertyVisualizer::widget;
55 VectorFieldDifferenceWidget* w = new VectorFieldDifferenceWidget();
56 w->paramVectorFieldDifference->setTitle(QString("3D Vector Field Difference Parameters of ").append(PropertyVisualizer::propertyInfo.propName().c_str()));
57 PropertyVisualizer::widget = w;
58}
59
60namespace OVMPVVFD{ //OVMPVVFD = OpenMeshProeprtyVisualizerVectorFieldDifference
61
62template <typename MeshT>
63float scalarFn_norm_of_diff(const ACG::Vec3d &a, const ACG::Vec3d &b) {
64 return (a - b).norm();
65}
66
67template <typename MeshT>
68float scalarFn_diff_of_norms(const ACG::Vec3d &a, const ACG::Vec3d &b) {
69 return std::fabs(a.norm() - b.norm());
70}
71
72template <typename MeshT>
73float scalarFn_4_symm_diff(const ACG::Vec3d &a, const ACG::Vec3d &b) {
74 double alpha = std::acos((a|b) / a.norm() / b.norm());
75 alpha -= std::floor((alpha + M_PI_4) / M_PI_2) * M_PI_2;
76 return std::fabs(alpha);
77}
78
79
80template<typename Prop1, typename Prop2, float (*ScalarFn)(const typename Prop1::value_type &, const typename Prop2::value_type &)>
81class ScalarAssigner {
82 public:
83 ScalarAssigner(const Prop1 &prop1, const Prop2 &prop2) :
84 prop1(prop1), prop2(prop2) { }
85
86 template<typename Handle>
87 float operator() (const Handle &handle) const {
88 return ScalarFn(prop1[handle], prop2[handle]);
89 }
90
91 protected:
92 const Prop1 &prop1;
93 const Prop2 &prop2;
94};
95
96template<typename MeshT, typename IteratorT, typename PropHandle, float (*ScalarFn)(const ACG::Vec3d &, const ACG::Vec3d &)>
97void colorElements(int objectID,
98 PropHandle prop1,
99 PropHandle prop2,
100 IteratorT primitivesBegin,
101 IteratorT primitivesEnd)
102{
103
104 std::vector<float> scalars;
105 std::transform(primitivesBegin, primitivesEnd, std::back_inserter(scalars),
106 ScalarAssigner<PropHandle, PropHandle, ScalarFn >(prop1, prop2));
107
108 const float min = *std::min_element(scalars.begin(), scalars.end());
109 const float max = *std::max_element(scalars.begin(), scalars.end());
110
111 VolumeMeshObject<MeshT>* object;
112 PluginFunctions::getObject(objectID, object);
113
114 ACG::ColorCoder colCod(min,max,false);
115 for (std::vector<float>::iterator i = scalars.begin(); i != scalars.end(); ++i)
116 object->colors()[*(primitivesBegin++)] = colCod(*i);
117}
118
119}
120
121
122template <typename MeshT>
123void OVMPropertyVisualizerVectorFieldDifference<MeshT>::visualizeCellProp(bool _setDrawMode)
124{
125 MeshT* mesh = OVMPropertyVisualizer<MeshT>::mesh;
126
127 OpenVolumeMesh::CellPropertyT<ACG::Vec3d> prop1 = mesh->template request_cell_property<ACG::Vec3d>(PropertyVisualizer::propertyInfo.propName());
128 OpenVolumeMesh::CellPropertyT<ACG::Vec3d> prop2 = mesh->template request_cell_property<ACG::Vec3d>(propertyInfo2.propName());
129
130 VectorFieldDifferenceWidget* w = static_cast<VectorFieldDifferenceWidget*>(PropertyVisualizer::widget);
131
132 if (w->vecFieldDiff_4symm_rb->isChecked())
133 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::CellIter, OpenVolumeMesh::CellPropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_4_symm_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->cells_begin(), mesh->cells_end());
134 if (w->vecFieldDiff_diff_norm_rb->isChecked())
135 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::CellIter, OpenVolumeMesh::CellPropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_diff_of_norms<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->cells_begin(), mesh->cells_end());
136 if (w->vecFieldDiff_norm_diff_rb->isChecked())
137 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::CellIter, OpenVolumeMesh::CellPropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_norm_of_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->cells_begin(), mesh->cells_end());
138 if (_setDrawMode)
139 {
140 VolumeMeshObject<MeshT>* object;
141 PluginFunctions::getObject(OVMPropertyVisualizer<MeshT>::mObjectID, object);
142 object->setObjectDrawMode(OVMPropertyVisualizer<MeshT>::drawModes.cellsColoredPerCell);
143 }
144}
145
146template <typename MeshT>
147void OVMPropertyVisualizerVectorFieldDifference<MeshT>::visualizeFaceProp(bool _setDrawMode)
148{
149 MeshT* mesh = OVMPropertyVisualizer<MeshT>::mesh;
150
151 OpenVolumeMesh::FacePropertyT<ACG::Vec3d> prop1 = mesh->template request_face_property<ACG::Vec3d>(PropertyVisualizer::propertyInfo.propName());
152 OpenVolumeMesh::FacePropertyT<ACG::Vec3d> prop2 = mesh->template request_face_property<ACG::Vec3d>(propertyInfo2.propName());
153
154 VectorFieldDifferenceWidget* w = static_cast<VectorFieldDifferenceWidget*>(PropertyVisualizer::widget);
155
156 if (w->vecFieldDiff_4symm_rb->isChecked())
157 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::FaceIter, OpenVolumeMesh::FacePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_4_symm_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->faces_begin(), mesh->faces_end());
158 if (w->vecFieldDiff_diff_norm_rb->isChecked())
159 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::FaceIter, OpenVolumeMesh::FacePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_diff_of_norms<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->faces_begin(), mesh->faces_end());
160 if (w->vecFieldDiff_norm_diff_rb->isChecked())
161 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::FaceIter, OpenVolumeMesh::FacePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_norm_of_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->faces_begin(), mesh->faces_end());
162 if (_setDrawMode)
163 {
164 VolumeMeshObject<MeshT>* object;
165 PluginFunctions::getObject(OVMPropertyVisualizer<MeshT>::mObjectID, object);
166 object->setObjectDrawMode(OVMPropertyVisualizer<MeshT>::drawModes.facesColoredPerFace);
167 }
168}
169
170template <typename MeshT>
171void OVMPropertyVisualizerVectorFieldDifference<MeshT>::visualizeHalffaceProp(bool _setDrawMode)
172{
173 MeshT* mesh = OVMPropertyVisualizer<MeshT>::mesh;
174
175 OpenVolumeMesh::HalfFacePropertyT<ACG::Vec3d> prop1 = mesh->template request_halfface_property<ACG::Vec3d>(PropertyVisualizer::propertyInfo.propName());
176 OpenVolumeMesh::HalfFacePropertyT<ACG::Vec3d> prop2 = mesh->template request_halfface_property<ACG::Vec3d>(propertyInfo2.propName());
177
178 VectorFieldDifferenceWidget* w = static_cast<VectorFieldDifferenceWidget*>(PropertyVisualizer::widget);
179
180 if (w->vecFieldDiff_4symm_rb->isChecked())
181 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::HalfFaceIter, OpenVolumeMesh::HalfFacePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_4_symm_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->halffaces_begin(), mesh->halffaces_end());
182 if (w->vecFieldDiff_diff_norm_rb->isChecked())
183 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::HalfFaceIter, OpenVolumeMesh::HalfFacePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_diff_of_norms<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->halffaces_begin(), mesh->halffaces_end());
184 if (w->vecFieldDiff_norm_diff_rb->isChecked())
185 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::HalfFaceIter, OpenVolumeMesh::HalfFacePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_norm_of_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->halffaces_begin(), mesh->halffaces_end());
186 if (_setDrawMode)
187 {
188 VolumeMeshObject<MeshT>* object;
189 PluginFunctions::getObject(OVMPropertyVisualizer<MeshT>::mObjectID, object);
190 object->setObjectDrawMode(OVMPropertyVisualizer<MeshT>::drawModes.halffacesColoredPerHalfface);
191 }
192}
193
194template <typename MeshT>
195void OVMPropertyVisualizerVectorFieldDifference<MeshT>::visualizeEdgeProp(bool _setDrawMode)
196{
197 MeshT* mesh = OVMPropertyVisualizer<MeshT>::mesh;
198
199 OpenVolumeMesh::EdgePropertyT<ACG::Vec3d> prop1 = mesh->template request_edge_property<ACG::Vec3d>(PropertyVisualizer::propertyInfo.propName());
200 OpenVolumeMesh::EdgePropertyT<ACG::Vec3d> prop2 = mesh->template request_edge_property<ACG::Vec3d>(propertyInfo2.propName());
201
202 VectorFieldDifferenceWidget* w = static_cast<VectorFieldDifferenceWidget*>(PropertyVisualizer::widget);
203
204 if (w->vecFieldDiff_4symm_rb->isChecked())
205 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::EdgeIter, OpenVolumeMesh::EdgePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_4_symm_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->edges_begin(), mesh->edges_end());
206 if (w->vecFieldDiff_diff_norm_rb->isChecked())
207 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::EdgeIter, OpenVolumeMesh::EdgePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_diff_of_norms<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->edges_begin(), mesh->edges_end());
208 if (w->vecFieldDiff_norm_diff_rb->isChecked())
209 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::EdgeIter, OpenVolumeMesh::EdgePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_norm_of_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->edges_begin(), mesh->edges_end());
210 if (_setDrawMode)
211 {
212 VolumeMeshObject<MeshT>* object;
213 PluginFunctions::getObject(OVMPropertyVisualizer<MeshT>::mObjectID, object);
214 object->setObjectDrawMode(OVMPropertyVisualizer<MeshT>::drawModes.edgesColoredPerEdge);
215 }
216}
217
218template <typename MeshT>
219void OVMPropertyVisualizerVectorFieldDifference<MeshT>::visualizeHalfedgeProp(bool _setDrawMode)
220{
221 MeshT* mesh = OVMPropertyVisualizer<MeshT>::mesh;
222
223 OpenVolumeMesh::HalfEdgePropertyT<ACG::Vec3d> prop1 = mesh->template request_halfedge_property<ACG::Vec3d>(PropertyVisualizer::propertyInfo.propName());
224 OpenVolumeMesh::HalfEdgePropertyT<ACG::Vec3d> prop2 = mesh->template request_halfedge_property<ACG::Vec3d>(propertyInfo2.propName());
225
226 VectorFieldDifferenceWidget* w = static_cast<VectorFieldDifferenceWidget*>(PropertyVisualizer::widget);
227
228 if (w->vecFieldDiff_4symm_rb->isChecked())
229 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::HalfEdgeIter, OpenVolumeMesh::HalfEdgePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_4_symm_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->halfedges_begin(), mesh->halfedges_end());
230 if (w->vecFieldDiff_diff_norm_rb->isChecked())
231 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::HalfEdgeIter, OpenVolumeMesh::HalfEdgePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_diff_of_norms<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->halfedges_begin(), mesh->halfedges_end());
232 if (w->vecFieldDiff_norm_diff_rb->isChecked())
233 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::HalfEdgeIter, OpenVolumeMesh::HalfEdgePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_norm_of_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->halfedges_begin(), mesh->halfedges_end());
234 if (_setDrawMode)
235 {
236 VolumeMeshObject<MeshT>* object;
237 PluginFunctions::getObject(OVMPropertyVisualizer<MeshT>::mObjectID, object);
238 object->setObjectDrawMode(OVMPropertyVisualizer<MeshT>::drawModes.halfedgesColoredPerHalfedge);
239 }
240}
241
242template <typename MeshT>
243void OVMPropertyVisualizerVectorFieldDifference<MeshT>::visualizeVertexProp(bool _setDrawMode)
244{
245 MeshT* mesh = OVMPropertyVisualizer<MeshT>::mesh;
246
247 OpenVolumeMesh::VertexPropertyT<ACG::Vec3d> prop1 = mesh->template request_vertex_property<ACG::Vec3d>(PropertyVisualizer::propertyInfo.propName());
248 OpenVolumeMesh::VertexPropertyT<ACG::Vec3d> prop2 = mesh->template request_vertex_property<ACG::Vec3d>(propertyInfo2.propName());
249
250 VectorFieldDifferenceWidget* w = static_cast<VectorFieldDifferenceWidget*>(PropertyVisualizer::widget);
251
252 if (w->vecFieldDiff_4symm_rb->isChecked())
253 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::VertexIter, OpenVolumeMesh::VertexPropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_4_symm_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->vertices_begin(), mesh->vertices_end());
254 if (w->vecFieldDiff_diff_norm_rb->isChecked())
255 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::VertexIter, OpenVolumeMesh::VertexPropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_diff_of_norms<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->vertices_begin(), mesh->vertices_end());
256 if (w->vecFieldDiff_norm_diff_rb->isChecked())
257 OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::VertexIter, OpenVolumeMesh::VertexPropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_norm_of_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->vertices_begin(), mesh->vertices_end());
258 if (_setDrawMode)
259 {
260 VolumeMeshObject<MeshT>* object;
261 PluginFunctions::getObject(OVMPropertyVisualizer<MeshT>::mObjectID, object);
262 object->setObjectDrawMode(OVMPropertyVisualizer<MeshT>::drawModes.verticesColored);
263 }
264}
265
ACG::ColorCoder
Class for generating nice colors for doubles.
Definition: ColorCoder.hh:68
OpenMesh::VectorT::norm
auto norm() const -> decltype(std::sqrt(std::declval< VectorT< S, DIM > >().sqrnorm()))
compute euclidean norm
Definition: Vector11T.hh:434
PropertyInfo
Cellection of information about a property.
Definition: Utils.hh:109
PluginFunctions::getObject
bool getObject(const int _identifier, BaseObject *&_object)
Get the object which has the given identifier.
Definition: PluginFunctions.cc:298