GSoC2023 - ACVD for Simplification & Remeshing

Documentation/doc/biblio/cgal_manual.bib

@@ -2386,6 +2386,46 @@ @article{ cgal:v-tm-95
   ,update       = "98.01 kettner"
 }
 
+@article{cgal:vc-acvdupmc-04,
+ author = {Valette, Sébastien and Chassery, Jean-Marc},
+ doi = {10.1111/j.1467-8659.2004.00769.x},
+ hal_id = {hal-00534535},
+ hal_version = {v1},
+ journal = {Computer Graphics Forum},
+ note = {Eurographics 2004 proceedings},
+ number = {3},
+ pages = {381-389},
+ pdf = {https://hal.archives-ouvertes.fr/hal-00534535/file/valette.pdf},
+ publisher = {Wiley},
+ title = {Approximated Centroidal Voronoi Diagrams for Uniform Polygonal Mesh Coarsening},
+ url = {https://hal.archives-ouvertes.fr/hal-00534535},
+ volume = {23},
+ year = {2004}
+}
+
+@inproceedings{cgal:audette2011approach,
+  title={Approach-guided controlled resolution brain meshing for fe-based interactive neurosurgery simulation},
+  author={Audette, M and Rivi{\`e}re, D and Ewend, M and Enquobahrie, A and Valette, S},
+  booktitle={Workshop on Mesh Processing in Medical Image Analysis, in conjunction with MICCAI 2011},
+  pages={176--186},
+  year={2011}
+}
+
+@article{cgal:vcp-grtmmdvd-08,
+ author = {Valette, Sébastien and Chassery, Jean-Marc and Prost, Rémy},
+ doi = {10.1109/TVCG.2007.70430},
+ hal_id = {hal-00537025},
+ hal_version = {v1},
+ journal = {IEEE Transactions on Visualization and Computer Graphics},
+ number = {2},
+ pages = {369--381},
+ pdf = {https://hal.archives-ouvertes.fr/hal-00537025/file/VCP08TVCG.pdf},
+ publisher = {Institute of Electrical and Electronics Engineers},
+ title = {Generic remeshing of 3D triangular meshes with metric-dependent discrete Voronoi Diagrams},
+ url = {https://hal.archives-ouvertes.fr/hal-00537025},
+ volume = {14},
+ year = {2008}
+}
 
 @techreport{cgal:vla-lod-15,
   title={LOD Generation for Urban Scenes},
@@ -3111,14 +3151,12 @@ @article{ecvp-bhhhk-14
   bibsource = {dblp computer science bibliography, https://dblp.org/}
 }
 
-@inproceedings {dunyach2013curvRemesh,
-  booktitle = {Eurographics 2013 - Short Papers},
-  title = {{Adaptive Remeshing for Real-Time Mesh Deformation}},
-  author = {Dunyach, Marion and Vanderhaeghe, David and Barthe, Loïc and Botsch, Mario},
-  year = {2013},
-  publisher = {The Eurographics Association},
-  ISSN = {1017-4656},
-  DOI = {10.2312/conf/EG2013/short/029-032}
+@inproceedings{dunyach2013curvRemesh,
+  title={Adaptive remeshing for real-time mesh deformation},
+  author={Dunyach, Marion and Vanderhaeghe, David and Barthe, Loïc and Botsch, Mario},
+  booktitle={Eurographics 2013},
+  year={2013},
+  organization={The Eurographics Association}
 }
 
 @book{botsch2010PMP,

Installation/include/CGAL/license/Polygon_mesh_processing/acvd.h

@@ -0,0 +1,54 @@
+// Copyright (c) 2016  GeometryFactory SARL (France).
+// All rights reserved.
+//
+// This file is part of CGAL (www.cgal.org)
+//
+// $URL$
+// $Id$
+// SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-Commercial
+//
+// Author(s) : Andreas Fabri
+//
+// Warning: this file is generated, see include/CGAL/license/README.md
+
+#ifndef CGAL_LICENSE_POLYGON_MESH_PROCESSING_ACVD_H
+#define CGAL_LICENSE_POLYGON_MESH_PROCESSING_ACVD_H
+
+#include <CGAL/config.h>
+#include <CGAL/license.h>
+
+#ifdef CGAL_POLYGON_MESH_PROCESSING_ACVD_COMMERCIAL_LICENSE
+
+#  if CGAL_POLYGON_MESH_PROCESSING_ACVD_COMMERCIAL_LICENSE < CGAL_RELEASE_DATE
+
+#    if defined(CGAL_LICENSE_WARNING)
+
+       CGAL_pragma_warning("Your commercial license for CGAL does not cover "
+                           "this release of the Polygon Mesh Processing - ACVD remeshing package.")
+#    endif
+
+#    ifdef CGAL_LICENSE_ERROR
+#      error "Your commercial license for CGAL does not cover this release \
+              of the Polygon Mesh Processing - ACVD remeshing package. \
+              You get this error, as you defined CGAL_LICENSE_ERROR."
+#    endif // CGAL_LICENSE_ERROR
+
+#  endif // CGAL_POLYGON_MESH_PROCESSING_ACVD_COMMERCIAL_LICENSE < CGAL_RELEASE_DATE
+
+#else // no CGAL_POLYGON_MESH_PROCESSING_ACVD_COMMERCIAL_LICENSE
+
+#  if defined(CGAL_LICENSE_WARNING)
+     CGAL_pragma_warning("\nThe macro CGAL_POLYGON_MESH_PROCESSING_ACVD_COMMERCIAL_LICENSE is not defined."
+                          "\nYou use the CGAL Polygon Mesh Processing - ACVD remeshing package under "
+                          "the terms of the GPLv3+.")
+#  endif // CGAL_LICENSE_WARNING
+
+#  ifdef CGAL_LICENSE_ERROR
+#    error "The macro CGAL_POLYGON_MESH_PROCESSING_ACVD_COMMERCIAL_LICENSE is not defined.\
+            You use the CGAL Polygon Mesh Processing - ACVD remeshing package under the terms of \
+            the GPLv3+. You get this error, as you defined CGAL_LICENSE_ERROR."
+#  endif // CGAL_LICENSE_ERROR
+
+#endif // no CGAL_POLYGON_MESH_PROCESSING_ACVD_COMMERCIAL_LICENSE
+
+#endif // CGAL_LICENSE_POLYGON_MESH_PROCESSING_ACVD_H

Installation/include/CGAL/license/gpl_package_list.txt

@@ -67,6 +67,7 @@ Polygon_mesh_processing/miscellaneous Polygon Mesh Processing - Miscellaneous
 Polygon_mesh_processing/detect_features Polygon Mesh Processing - Feature Detection
 Polygon_mesh_processing/collision_detection Polygon Mesh Processing - Collision Detection
 Polygon_mesh_processing/autorefinement Polygon Mesh Processing - Autorefinement
+Polygon_mesh_processing/acvd Polygon Mesh Processing - ACVD remeshing
 Polyhedron 3D Polyhedral Surface
 Polyline_simplification_2 2D Polyline Simplification
 Polytope_distance_d Optimal Distances

Lab/demo/Lab/Plugins/PMP/ACVD_remeshing_dialog.ui

@@ -0,0 +1,128 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<ui version="4.0">
+ <class>ACVDDialog</class>
+ <widget class="QDialog" name="ACVDDialog">
+  <property name="geometry">
+   <rect>
+    <x>0</x>
+    <y>0</y>
+    <width>400</width>
+    <height>241</height>
+   </rect>
+  </property>
+  <property name="windowTitle">
+   <string>ACVD Remeshing</string>
+  </property>
+  <layout class="QVBoxLayout" name="verticalLayout">
+   <item>
+    <widget class="QLabel" name="label_3">
+     <property name="text">
+      <string>Target number of Vertices (Clusters):</string>
+     </property>
+    </widget>
+   </item>
+   <item>
+    <widget class="DoubleEdit" name="nb_clusters_spin_box">
+     <property name="value" stdset="0">
+      <number>100</number>
+     </property>
+    </widget>
+   </item>
+   <item>
+    <widget class="QLabel" name="label_2">
+     <property name="text">
+      <string>Cluster Metric:</string>
+     </property>
+    </widget>
+   </item>
+   <item>
+    <widget class="QRadioButton" name="IsotropicClustering">
+     <property name="text">
+      <string>Isotropic Clustering</string>
+     </property>
+     <property name="checked">
+      <bool>true</bool>
+     </property>
+     <attribute name="buttonGroup">
+      <string notr="true">buttonGroup</string>
+     </attribute>
+    </widget>
+   </item>
+   <item>
+    <widget class="QRadioButton" name="QEMClustering">
+     <property name="text">
+      <string>QEM Clustering</string>
+     </property>
+     <attribute name="buttonGroup">
+      <string notr="true">buttonGroup</string>
+     </attribute>
+    </widget>
+   </item>
+   <item>
+    <widget class="QRadioButton" name="AnisotropicClustering">
+     <property name="text">
+      <string>Anisotropic Clustering</string>
+     </property>
+     <attribute name="buttonGroup">
+      <string notr="true">buttonGroup</string>
+     </attribute>
+    </widget>
+   </item>
+   <item>
+    <widget class="QDialogButtonBox" name="buttonBox">
+     <property name="orientation">
+      <enum>Qt::Horizontal</enum>
+     </property>
+     <property name="standardButtons">
+      <set>QDialogButtonBox::Cancel|QDialogButtonBox::Ok</set>
+     </property>
+    </widget>
+   </item>
+  </layout>
+ </widget>
+ <customwidgets>
+  <customwidget>
+   <class>DoubleEdit</class>
+   <extends>QLineEdit</extends>
+   <header>CGAL_double_edit.h</header>
+  </customwidget>
+ </customwidgets>
+ <resources/>
+ <connections>
+  <connection>
+   <sender>buttonBox</sender>
+   <signal>accepted()</signal>
+   <receiver>ACVDDialog</receiver>
+   <slot>accept()</slot>
+   <hints>
+    <hint type="sourcelabel">
+     <x>248</x>
+     <y>254</y>
+    </hint>
+    <hint type="destinationlabel">
+     <x>157</x>
+     <y>274</y>
+    </hint>
+   </hints>
+  </connection>
+  <connection>
+   <sender>buttonBox</sender>
+   <signal>rejected()</signal>
+   <receiver>ACVDDialog</receiver>
+   <slot>reject()</slot>
+   <hints>
+    <hint type="sourcelabel">
+     <x>316</x>
+     <y>260</y>
+    </hint>
+    <hint type="destinationlabel">
+     <x>286</x>
+     <y>274</y>
+    </hint>
+   </hints>
+  </connection>
+ </connections>
+ <buttongroups>
+  <buttongroup name="buttonGroup"/>
+ </buttongroups>
+</ui>

Lab/demo/Lab/Plugins/PMP/ACVD_remeshing_plugin.cpp

@@ -0,0 +1,124 @@
+#include <QtCore/qglobal.h>
+
+#include "Messages_interface.h"
+#include <CGAL/Three/CGAL_Lab_plugin_helper.h>
+#include <CGAL/Three/Three.h>
+#include "Scene_polyhedron_selection_item.h"
+#include "ui_ACVD_remeshing_dialog.h"
+
+#include "SMesh_type.h"
+typedef Scene_surface_mesh_item Scene_facegraph_item;
+
+#include <CGAL/iterator.h>
+#include <CGAL/Polygon_mesh_processing/approximated_centroidal_Voronoi_diagram_remeshing.h>
+
+#include <QElapsedTimer>
+#include <QRadioButton>
+#include <QCheckBox>
+#include <QAction>
+#include <QMainWindow>
+#include <QApplication>
+#include <QDockWidget>
+#include <QEvent>
+#include <QKeyEvent>
+#include <QMouseEvent>
+#include <QFileDialog>
+
+#include <vector>
+#include <algorithm>
+#include <queue>
+
+typedef Scene_facegraph_item::Face_graph FaceGraph;
+
+using namespace CGAL::Three;
+namespace PMP = CGAL::Polygon_mesh_processing;
+
+class CGAL_Lab_acvd_remeshing_plugin :
+  public QObject,
+  public CGAL_Lab_plugin_helper
+{
+  Q_OBJECT
+  Q_INTERFACES(CGAL::Three::CGAL_Lab_plugin_interface)
+  Q_PLUGIN_METADATA(IID "com.geometryfactory.PolyhedronDemo.PluginInterface/1.0" FILE "acvd_remeshing_plugin.json")
+public:
+  bool applicable(QAction*) const {
+    return qobject_cast<Scene_facegraph_item*>(scene->item(scene->mainSelectionIndex()));
+  }
+  void print_message(QString message) { CGAL::Three::Three::information(message);}
+  QList<QAction*> actions() const { return QList<QAction*>() << actionACVD; }
+
+
+  void init(QMainWindow* mainWindow, CGAL::Three::Scene_interface* scene_interface, Messages_interface*) {
+    mw = mainWindow;
+    scene = scene_interface;
+    actionACVD = new QAction(tr(
+                                "ACVD Remeshing"
+                                ), mw);
+    //actionACVD->setProperty("subMenuName", "");
+    if (actionACVD)
+      connect(actionACVD, SIGNAL(triggered()), this, SLOT(acvd_action()));
+  }
+
+public Q_SLOTS:
+  void acvd_action() {
+
+    // Create dialog box
+    QDialog dialog(mw);
+    ui.setupUi(&dialog);
+
+    ui.nb_clusters_spin_box->setMinimum(10);
+    ui.nb_clusters_spin_box->setMaximum(1000000);
+    // get number of vertices of the mesh
+    Scene_facegraph_item* item = getSelectedItem<Scene_facegraph_item>();
+    if(!item) { return; }
+    // set default number of clusters to 5% of the number of vertices
+    ui.nb_clusters_spin_box->setValue((std::max)(item->face_graph()->number_of_faces() / 20, (unsigned int)10));
+
+    int i = dialog.exec();
+    if (i == QDialog::Rejected)
+    {
+      std::cout << "Remeshing aborted" << std::endl;
+      return;
+    }
+
+
+    Face_graph* graph = item->face_graph();
+    if(!graph) return;
+    QElapsedTimer time;
+    time.start();
+    CGAL::Three::Three::information("ACVD Remeshing...");
+
+
+    QApplication::setOverrideCursor(Qt::WaitCursor);
+    FaceGraph remeshed = *graph;
+
+    // TODO add post-processing
+
+    if (ui.IsotropicClustering->isChecked())
+      PMP::approximated_centroidal_Voronoi_diagram_remeshing(remeshed, ui.nb_clusters_spin_box->value());
+    else if (ui.QEMClustering->isChecked())
+      PMP::approximated_centroidal_Voronoi_diagram_remeshing(remeshed, ui.nb_clusters_spin_box->value(), CGAL::parameters::use_qem_based_energy(true));
+    else
+      PMP::approximated_centroidal_Voronoi_diagram_remeshing(remeshed, ui.nb_clusters_spin_box->value(), CGAL::parameters::gradation_factor(0.8));
+
+    // update the scene
+    CGAL::Three::Three::information(QString("ok (%1 ms)").arg(time.elapsed()));
+    QApplication::restoreOverrideCursor();
+    item->invalidateOpenGLBuffers();
+
+    // add the simplified mesh to the scene
+    Scene_facegraph_item* remeshed_item = new Scene_facegraph_item(remeshed);
+    remeshed_item->setName(QString(item->name() + "_simplified"));
+    scene->addItem(remeshed_item);
+    remeshed_item->invalidateOpenGLBuffers();
+    scene->itemChanged(remeshed_item);
+    item->setVisible(false);
+  }
+
+private:
+  QAction* actionACVD;
+  Ui::ACVDDialog ui;
+}; // end CGAL_Lab_acvd_remeshing_plugin
+
+
+#include "ACVD_remeshing_plugin.moc"

Lab/demo/Lab/Plugins/PMP/CMakeLists.txt

@@ -29,6 +29,11 @@ target_link_libraries(extrude_plugin PRIVATE scene_surface_mesh_item
 
 if(TARGET CGAL::Eigen3_support)
   if("${EIGEN3_VERSION}" VERSION_GREATER "3.1.90")
+
+    qt6_wrap_ui( acvd_remeshingUI_FILES ACVD_remeshing_dialog.ui)
+    cgal_lab_plugin(acvd_remeshing_plugin ACVD_remeshing_plugin ${acvd_remeshingUI_FILES} KEYWORDS PMP)
+    target_link_libraries(acvd_remeshing_plugin PUBLIC scene_surface_mesh_item scene_points_with_normal_item scene_polygon_soup_item CGAL::Eigen3_support)
+
     qt6_wrap_ui( hole_fillingUI_FILES Hole_filling_widget.ui)
     cgal_lab_plugin(hole_filling_plugin Hole_filling_plugin ${hole_fillingUI_FILES} KEYWORDS PMP)
     target_link_libraries(hole_filling_plugin PRIVATE scene_surface_mesh_item scene_polylines_item scene_selection_item CGAL::Eigen3_support)