diff --git a/Documentation/doc/biblio/cgal_manual.bib b/Documentation/doc/biblio/cgal_manual.bib
index 1c49fe6ffd2..1acdc35bcfd 100644
--- a/Documentation/doc/biblio/cgal_manual.bib
+++ b/Documentation/doc/biblio/cgal_manual.bib
@@ -3262,6 +3262,15 @@ pages = "207--221"
   year={1998}
 }
 
+
+@article{cgal:hssz-gmcabonbc-97,
+  title={A generalized marching cubes algorithm based on non-binary classifications},
+  author={H-C. Hege and M. Seebass and D. Stalling and M. Zöckler},
+  journal={Preprint},
+  number={SC 97-05},
+  year={1997}
+}
+
 % ----------------------------------------------------------------------------
 % END OF BIBFILE
 % ----------------------------------------------------------------------------
diff --git a/Mesh_3/doc/Mesh_3/Mesh_3.txt b/Mesh_3/doc/Mesh_3/Mesh_3.txt
index 29ba32bd500..232582c3a10 100644
--- a/Mesh_3/doc/Mesh_3/Mesh_3.txt
+++ b/Mesh_3/doc/Mesh_3/Mesh_3.txt
@@ -735,6 +735,24 @@ Surface of the output mesh generated with a very small `facet_distance`
 without the weights (left, 25563 vertices) and with the weights (right, 19936 vertices).
 \cgalFigureEnd
 
+\subsubsection Mesh_3DomainsFrom3DImagesWithFeatures Domains From Segmented 3D Images, with Feature Detection
+
+The example \ref Mesh_3/mesh_3D_image_with_detection_of_features.cpp shows how to use
+the feature detection functors `CGAL::Mesh_3::Detect_features_in_image` and
+`CGAL::Mesh_3::Detect_features_on_image_bbox`, in order to get a quality discretization
+of the curves lying at the intersection of three or more subdomains (including the outside).
+The original feature detection algorithm was described in \cgalCite{cgal:hssz-gmcabonbc-97},
+along with a list of possible voxel configurations. The detection implemented in \cgal
+generalizes this description and computes a graph of polyline features at the Intersection
+of three or more subdomains.
+
+\cgalExample{Mesh_3/mesh_3D_image_with_detection_of_features.cpp}
+
+\cgalFigureBegin{figure_detectFeaturesInImage, detect_features_in_image.jpg}
+Surface of the output mesh generated from a labeled image without (left)
+and with (right) feature detection and protection.
+\cgalFigureEnd
+
 
 \subsubsection Mesh_3DomainsFrom3DImagesWithCustomInitialization Domains From 3D Images, with a Custom Initialization
 
diff --git a/Mesh_3/doc/Mesh_3/fig/detect_features_in_image.jpg b/Mesh_3/doc/Mesh_3/fig/detect_features_in_image.jpg
new file mode 100644
index 00000000000..a8d96449c71
Binary files /dev/null and b/Mesh_3/doc/Mesh_3/fig/detect_features_in_image.jpg differ