diff --git a/_posts/2024-07-05-mesh3-improvements.md b/_posts/2024-07-05-mesh3-improvements.md new file mode 100644 index 000000000..6f24170c9 --- /dev/null +++ b/_posts/2024-07-05-mesh3-improvements.md @@ -0,0 +1,111 @@ +--- +layout: post +title: "Mesh_3 improvements" +description: "Mesh_3 improvements in CGAL-6.0" +category: "news" +tags: [""] +--- +{% include JB/setup %} + +

Jane Tournois

+

GeometryFactory

+ +
+

Introduced in CGAL 3.5, the package +3D Mesh Generation +is designed to generate a isotropic simplicial meshes discretizing 3D domains. +Since then, multiple improvements in mesh quality, performances, and versatility have been made available in CGAL. +Lately, CGAL authors have focused more on the quality of the output mesh, the ability of the +mesh generation algorithm to represent the input surface in a more accurate way, for a wider range of input data types. +

+ + + +

+

Recent releases of CGAL have introduced improvements to the 3D Tetrahedral Mesh Generation package: +detection and protection of triple lines from labeled images, +use of weighted labeled images, +dealing with self-intersecting input surfaces, +and a new approximation error criterion for polyline features.

+ +
+

Triple lines from labeled images

+ +

CGAL's tetrahedral mesh generation algorithm is able to compute tetrahedral meshes from 3D segmented images. +For all types of inputs, detection and construction of one-dimensional features must be done a priori for the +feature protection algorithm to be able to represent them properly in the generated mesh. +In 3D segmented images, such features correspond to the intersections of the bounding box of the image with the surfaces defined by the image, as well as polylines that lie at the intersection of three or more subdomains (including the outside).

+ +

The article +Hege, H. C., Seebass, M., Stalling, D., & Zöckler, M. (1997), A generalized marching cubes algorithm based on non-binary classifications, provides a list of possible voxel configurations. This description has been generalized for tetrahedral mesh generation needs, and +implemented in CGAL.

+ +

After feature lines detection, similarly to what is done for polyhedral surfaces, the polyline features +are added to the feature graph of the domain for feature protection.

+ +

Two detection functors are provided : +

+

+ +
+

Weighted Labeled Images

+ +

When a segmented image is given as input for tetrahedral mesh generation, +the generated mesh surface sometimes sticks too closely to the voxels surface, causing an aliasing effect. A solution to generate +a smooth and accurate output surface was described by Stalling et al in +Stalling, D., Zöckler, M., Sander, O., & Hege, H. C. (1998), Weighted labels for 3D image segmentation. +It consists in generating a second input image, made of integer coefficients called weights, and use those weights to define +smoother domain boundaries. The 3D image of weights can be generated using CGAL::Mesh_3::generate_label_weights().

+ +
+
+
Mesh generation from labeled image without weights (left, 88,000 vertices) + and with weights (right, 67,000 vertices). +
+ +
+

The following example shows how to use weighted labeled images: mesh_3D_weighted_image.cpp.

+ +

Utilizing tetrahedral mesh generation from labeled images, incorporating both feature detection and weighted labels, +presents challenges. A new method that combines feature detection and weighted labels has been +designed and implemented in CGAL, achieving smooth and feature preserving output meshes. +This method leverages the best of these two new functionalities.

+ + + + + +
+

Self-intersecting input surfaces

+ + +
+

Approximation Criterion on Feature Edges

+ + +

Status

+ +

All these new improvements are already integrated in CGAL's master branch on the +CGAL GitHub repository, and has just been officially released +in the latest version of CGAL, CGAL 6.0-beta1.

+ + +Documentation of the package Mesh_3 +
+ +CGAL master branch on GitHub diff --git a/images/mesh3_weighted_images.png b/images/mesh3_weighted_images.png new file mode 100644 index 000000000..69c19c4da Binary files /dev/null and b/images/mesh3_weighted_images.png differ