1 // Created on: 1993-09-23
2 // Created by: Jean-Louis FRENKEL
3 // Copyright (c) 1993-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
6 // This file is part of Open CASCADE Technology software library.
8 // This library is free software; you can redistribute it and/or modify it under
9 // the terms of the GNU Lesser General Public License version 2.1 as published
10 // by the Free Software Foundation, with special exception defined in the file
11 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
12 // distribution for complete text of the license and disclaimer of any warranty.
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
17 #include <StdPrs_ShadedShape.hxx>
19 #include <Bnd_Box.hxx>
20 #include <BRepTools.hxx>
21 #include <BRepBndLib.hxx>
22 #include <BRep_Builder.hxx>
23 #include <BRep_Tool.hxx>
24 #include <Graphic3d_ArrayOfSegments.hxx>
25 #include <Graphic3d_ArrayOfTriangles.hxx>
26 #include <Graphic3d_AspectFillArea3d.hxx>
27 #include <Graphic3d_Group.hxx>
31 #include <NCollection_List.hxx>
32 #include <Precision.hxx>
33 #include <Prs3d_Drawer.hxx>
34 #include <Prs3d_IsoAspect.hxx>
35 #include <Prs3d_LineAspect.hxx>
36 #include <Prs3d_Presentation.hxx>
37 #include <Prs3d_ShadingAspect.hxx>
38 #include <Poly_Connect.hxx>
39 #include <Poly_PolygonOnTriangulation.hxx>
40 #include <Poly_Triangulation.hxx>
41 #include <StdPrs_ToolTriangulatedShape.hxx>
42 #include <StdPrs_WFShape.hxx>
44 #include <TopExp_Explorer.hxx>
46 #include <TopoDS_Compound.hxx>
47 #include <TopoDS_Face.hxx>
48 #include <TopoDS_Shape.hxx>
49 #include <TColgp_Array1OfDir.hxx>
50 #include <TColgp_Array1OfPnt2d.hxx>
51 #include <TColgp_HArray1OfPnt.hxx>
52 #include <TopTools_ListOfShape.hxx>
53 #include <TopTools_IndexedDataMapOfShapeListOfShape.hxx>
58 //! Computes wireframe presentation for free wires and vertices
59 void wireframeFromShape (const Handle (Prs3d_Presentation)& thePrs,
60 const TopoDS_Shape& theShape,
61 const Handle (Prs3d_Drawer)& theDrawer)
63 TopExp_Explorer aShapeIter (theShape, TopAbs_FACE);
64 if (!aShapeIter.More())
66 StdPrs_WFShape::Add (thePrs, theShape, theDrawer);
70 const Standard_Boolean aDrawAllVerticesFlag = (theDrawer->VertexDrawMode() == Prs3d_VDM_All);
71 if (!aDrawAllVerticesFlag && theShape.ShapeType() != TopAbs_COMPOUND)
76 // We have to create a compound and collect all subshapes not drawn by the shading algo.
79 // - isolated vertices, if aDrawAllVerticesFlag == Standard_False
80 // - all shape's vertices, if aDrawAllVerticesFlag == Standard_True
81 TopoDS_Compound aCompoundWF;
82 BRep_Builder aBuilder;
83 aBuilder.MakeCompound (aCompoundWF);
84 Standard_Boolean hasElement = Standard_False;
87 for (aShapeIter.Init (theShape, TopAbs_EDGE, TopAbs_FACE); aShapeIter.More(); aShapeIter.Next())
89 hasElement = Standard_True;
90 aBuilder.Add (aCompoundWF, aShapeIter.Current());
92 // isolated or all vertices
93 aShapeIter.Init (theShape, TopAbs_VERTEX, aDrawAllVerticesFlag ? TopAbs_SHAPE : TopAbs_EDGE);
94 for (; aShapeIter.More(); aShapeIter.Next())
96 hasElement = Standard_True;
97 aBuilder.Add (aCompoundWF, aShapeIter.Current());
101 StdPrs_WFShape::Add (thePrs, aCompoundWF, theDrawer);
105 //! Computes special wireframe presentation for faces without triangulation.
106 void wireframeNoTriangFacesFromShape (const Handle(Prs3d_Presentation)& thePrs,
107 const TopoDS_Shape& theShape,
108 const Handle(Prs3d_Drawer)& theDrawer)
110 TopoDS_Compound aCompoundWF;
111 BRep_Builder aBuilder;
112 aBuilder.MakeCompound (aCompoundWF);
113 TopLoc_Location aLoc;
114 Standard_Boolean hasElement = Standard_False;
116 for (TopExp_Explorer aShapeIter(theShape, TopAbs_FACE); aShapeIter.More(); aShapeIter.Next())
118 const TopoDS_Face& aFace = TopoDS::Face (aShapeIter.Current());
119 const Handle(Poly_Triangulation) aTriang = BRep_Tool::Triangulation (aFace, aLoc);
120 if (aTriang.IsNull())
122 hasElement = Standard_True;
123 aBuilder.Add (aCompoundWF, aFace);
129 Standard_Integer aPrevUIsoNb = theDrawer->UIsoAspect()->Number();
130 Standard_Integer aPrevVIsoNb = theDrawer->VIsoAspect()->Number();
131 theDrawer->UIsoAspect()->SetNumber (5);
132 theDrawer->VIsoAspect()->SetNumber (5);
134 StdPrs_WFShape::Add (thePrs, aCompoundWF, theDrawer);
136 theDrawer->UIsoAspect()->SetNumber (aPrevUIsoNb);
137 theDrawer->VIsoAspect()->SetNumber (aPrevVIsoNb);
141 //! Gets triangulation of every face of shape and fills output array of triangles
142 static Handle(Graphic3d_ArrayOfTriangles) fillTriangles (const TopoDS_Shape& theShape,
143 const Standard_Boolean theHasTexels,
144 const gp_Pnt2d& theUVOrigin,
145 const gp_Pnt2d& theUVRepeat,
146 const gp_Pnt2d& theUVScale)
148 Handle(Poly_Triangulation) aT;
149 TopLoc_Location aLoc;
151 Standard_Integer aNbTriangles = 0;
152 Standard_Integer aNbVertices = 0;
154 // Precision for compare square distances
155 const Standard_Real aPreci = Precision::SquareConfusion();
157 TopExp_Explorer aFaceIt(theShape, TopAbs_FACE);
158 for (; aFaceIt.More(); aFaceIt.Next())
160 const TopoDS_Face& aFace = TopoDS::Face(aFaceIt.Current());
161 aT = BRep_Tool::Triangulation (aFace, aLoc);
164 aNbTriangles += aT->NbTriangles();
165 aNbVertices += aT->NbNodes();
168 if (aNbVertices < 3 || aNbTriangles <= 0)
170 return Handle(Graphic3d_ArrayOfTriangles)();
173 Handle(Graphic3d_ArrayOfTriangles) anArray = new Graphic3d_ArrayOfTriangles (aNbVertices, 3 * aNbTriangles,
174 Standard_True, Standard_False, theHasTexels);
175 Standard_Real aUmin (0.0), aUmax (0.0), aVmin (0.0), aVmax (0.0), dUmax (0.0), dVmax (0.0);
176 for (aFaceIt.Init (theShape, TopAbs_FACE); aFaceIt.More(); aFaceIt.Next())
178 const TopoDS_Face& aFace = TopoDS::Face(aFaceIt.Current());
179 aT = BRep_Tool::Triangulation (aFace, aLoc);
184 const gp_Trsf& aTrsf = aLoc.Transformation();
186 // Determinant of transform matrix less then 0 means that mirror transform applied.
187 Standard_Boolean isMirrored = aTrsf.VectorialPart().Determinant() < 0;
189 // Extracts vertices & normals from nodes
190 const TColgp_Array1OfPnt& aNodes = aT->Nodes();
191 const TColgp_Array1OfPnt2d& aUVNodes = aT->UVNodes();
192 StdPrs_ToolTriangulatedShape::ComputeNormals (aFace, aT);
193 const TShort_Array1OfShortReal& aNormals = aT->Normals();
194 const Standard_ShortReal* aNormArr = &aNormals.First();
198 BRepTools::UVBounds (aFace, aUmin, aUmax, aVmin, aVmax);
199 dUmax = (aUmax - aUmin);
200 dVmax = (aVmax - aVmin);
203 const Standard_Integer aDecal = anArray->VertexNumber();
204 for (Standard_Integer aNodeIter = aNodes.Lower(); aNodeIter <= aNodes.Upper(); ++aNodeIter)
206 aPoint = aNodes (aNodeIter);
207 const Standard_Integer anId = 3 * (aNodeIter - aNodes.Lower());
208 gp_Dir aNorm (aNormArr[anId + 0], aNormArr[anId + 1], aNormArr[anId + 2]);
209 if (aFace.Orientation() == TopAbs_REVERSED)
213 if (!aLoc.IsIdentity())
215 aPoint.Transform (aTrsf);
216 aNorm.Transform (aTrsf);
219 if (theHasTexels && aUVNodes.Upper() == aNodes.Upper())
221 const gp_Pnt2d aTexel = (dUmax == 0.0 || dVmax == 0.0)
222 ? aUVNodes (aNodeIter)
223 : gp_Pnt2d ((-theUVOrigin.X() + (theUVRepeat.X() * (aUVNodes(aNodeIter).X() - aUmin)) / dUmax) / theUVScale.X(),
224 (-theUVOrigin.Y() + (theUVRepeat.Y() * (aUVNodes(aNodeIter).Y() - aVmin)) / dVmax) / theUVScale.Y());
225 anArray->AddVertex (aPoint, aNorm, aTexel);
229 anArray->AddVertex (aPoint, aNorm);
233 // Fill array with vertex and edge visibility info
234 const Poly_Array1OfTriangle& aTriangles = aT->Triangles();
235 Standard_Integer anIndex[3];
236 for (Standard_Integer aTriIter = 1; aTriIter <= aT->NbTriangles(); ++aTriIter)
238 if ((aFace.Orientation() == TopAbs_REVERSED) ^ isMirrored)
240 aTriangles (aTriIter).Get (anIndex[0], anIndex[2], anIndex[1]);
244 aTriangles (aTriIter).Get (anIndex[0], anIndex[1], anIndex[2]);
247 gp_Pnt aP1 = aNodes (anIndex[0]);
248 gp_Pnt aP2 = aNodes (anIndex[1]);
249 gp_Pnt aP3 = aNodes (anIndex[2]);
251 gp_Vec aV1 (aP1, aP2);
252 if (aV1.SquareMagnitude() <= aPreci)
256 gp_Vec aV2 (aP2, aP3);
257 if (aV2.SquareMagnitude() <= aPreci)
261 gp_Vec aV3 (aP3, aP1);
262 if (aV3.SquareMagnitude() <= aPreci)
269 if (aV1.SquareMagnitude() > aPreci)
271 anArray->AddEdges (anIndex[0] + aDecal,
273 anIndex[2] + aDecal);
280 //! Prepare shaded presentation for specified shape
281 static Standard_Boolean shadeFromShape (const TopoDS_Shape& theShape,
282 const Handle(Prs3d_Presentation)& thePrs,
283 const Handle(Prs3d_Drawer)& theDrawer,
284 const Standard_Boolean theHasTexels,
285 const gp_Pnt2d& theUVOrigin,
286 const gp_Pnt2d& theUVRepeat,
287 const gp_Pnt2d& theUVScale,
288 const bool theIsClosed)
290 Handle(Graphic3d_ArrayOfTriangles) aPArray = fillTriangles (theShape, theHasTexels, theUVOrigin, theUVRepeat, theUVScale);
291 if (aPArray.IsNull())
293 return Standard_False;
296 Handle(Graphic3d_Group) aGroup = Prs3d_Root::NewGroup (thePrs);
297 aGroup->SetClosed (theIsClosed);
298 aGroup->SetGroupPrimitivesAspect (theDrawer->ShadingAspect()->Aspect());
299 aGroup->AddPrimitiveArray (aPArray);
300 return Standard_True;
303 //! Compute boundary presentation for faces of the shape.
304 static Handle(Graphic3d_ArrayOfSegments) fillFaceBoundaries (const TopoDS_Shape& theShape)
306 // collection of all triangulation nodes on edges
307 // for computing boundaries presentation
308 Standard_Integer aNodeNumber = 0;
309 Standard_Integer aNbPolylines = 0;
311 TopLoc_Location aTrsf;
312 TColgp_SequenceOfPnt aSeqPntsExtra;
313 for (TopExp_Explorer aFaceIter (theShape, TopAbs_FACE); aFaceIter.More(); aFaceIter.Next())
315 const TopoDS_Face& aFace = TopoDS::Face (aFaceIter.Current());
316 if (aFace.NbChildren() == 0)
318 // handle specifically faces without boundary definition (triangulation-only)
319 StdPrs_WFShape::AddEdgesOnTriangulation (aSeqPntsExtra, aFace, Standard_False);
323 // explore all boundary edges
324 TopTools_IndexedDataMapOfShapeListOfShape anEdgesMap;
325 TopExp::MapShapesAndAncestors (theShape, TopAbs_EDGE, TopAbs_FACE, anEdgesMap);
326 for (TopTools_IndexedDataMapOfShapeListOfShape::Iterator anEdgeIter (anEdgesMap); anEdgeIter.More(); anEdgeIter.Next())
329 if (anEdgeIter.Value().Extent() == 0)
334 // take one of the shared edges and get edge triangulation
335 const TopoDS_Face& aFace = TopoDS::Face (anEdgeIter.Value().First());
336 Handle(Poly_Triangulation) aTriangulation = BRep_Tool::Triangulation (aFace, aTrsf);
337 if (aTriangulation.IsNull())
342 const TopoDS_Edge& anEdge = TopoDS::Edge (anEdgeIter.Key());
343 Handle(Poly_PolygonOnTriangulation) anEdgePoly = BRep_Tool::PolygonOnTriangulation (anEdge, aTriangulation, aTrsf);
344 if (!anEdgePoly.IsNull()
345 && anEdgePoly->Nodes().Length() >= 2)
347 aNodeNumber += anEdgePoly->Nodes().Length();
351 if (aNodeNumber == 0)
353 if (aSeqPntsExtra.Size() < 2)
355 return Handle(Graphic3d_ArrayOfSegments)();
358 Standard_Integer aNbVertices = aSeqPntsExtra.Size();
359 Handle(Graphic3d_ArrayOfSegments) aSegments = new Graphic3d_ArrayOfSegments (aNbVertices);
360 for (Standard_Integer aPntIter = 1; aPntIter <= aNbVertices; aPntIter += 2)
362 aSegments->AddVertex (aSeqPntsExtra.Value (aPntIter));
363 aSegments->AddVertex (aSeqPntsExtra.Value (aPntIter + 1));
368 // create indexed segments array to pack polylines from different edges into single array
369 const Standard_Integer aSegmentEdgeNb = (aNodeNumber - aNbPolylines) * 2;
370 Handle(Graphic3d_ArrayOfSegments) aSegments = new Graphic3d_ArrayOfSegments (aNodeNumber + aSeqPntsExtra.Size(), aSegmentEdgeNb + aSeqPntsExtra.Size());
371 for (TopTools_IndexedDataMapOfShapeListOfShape::Iterator anEdgeIter (anEdgesMap); anEdgeIter.More(); anEdgeIter.Next())
373 if (anEdgeIter.Value().Extent() == 0)
378 const TopoDS_Face& aFace = TopoDS::Face (anEdgeIter.Value().First());
379 Handle(Poly_Triangulation) aTriangulation = BRep_Tool::Triangulation (aFace, aTrsf);
380 if (aTriangulation.IsNull())
385 const TopoDS_Edge& anEdge = TopoDS::Edge (anEdgeIter.Key());
386 Handle(Poly_PolygonOnTriangulation) anEdgePoly = BRep_Tool::PolygonOnTriangulation (anEdge, aTriangulation, aTrsf);
387 if (anEdgePoly.IsNull()
388 || anEdgePoly->Nodes().Length () < 2)
393 // get edge nodes indexes from face triangulation
394 const TColgp_Array1OfPnt& aTriNodes = aTriangulation->Nodes();
395 const TColStd_Array1OfInteger& anEdgeNodes = anEdgePoly->Nodes();
397 // collect the edge nodes
398 Standard_Integer aSegmentEdge = aSegments->VertexNumber() + 1;
399 for (Standard_Integer aNodeIdx = anEdgeNodes.Lower(); aNodeIdx <= anEdgeNodes.Upper(); ++aNodeIdx)
401 // node index in face triangulation
402 // get node and apply location transformation to the node
403 const Standard_Integer aTriIndex = anEdgeNodes.Value (aNodeIdx);
404 gp_Pnt aTriNode = aTriNodes.Value (aTriIndex);
405 if (!aTrsf.IsIdentity())
407 aTriNode.Transform (aTrsf);
410 aSegments->AddVertex (aTriNode);
411 if (aNodeIdx != anEdgeNodes.Lower())
413 aSegments->AddEdge ( aSegmentEdge);
414 aSegments->AddEdge (++aSegmentEdge);
420 Standard_Integer aSegmentEdge = aSegments->VertexNumber();
421 const Standard_Integer aNbVertices = aSeqPntsExtra.Size();
422 for (Standard_Integer aPntIter = 1; aPntIter <= aNbVertices; aPntIter += 2)
424 aSegments->AddVertex (aSeqPntsExtra.Value (aPntIter));
425 aSegments->AddEdge (++aSegmentEdge);
426 aSegments->AddVertex (aSeqPntsExtra.Value (aPntIter + 1));
427 aSegments->AddEdge (++aSegmentEdge);
434 } // anonymous namespace
436 // =======================================================================
437 // function : ExploreSolids
439 // =======================================================================
440 void StdPrs_ShadedShape::ExploreSolids (const TopoDS_Shape& theShape,
441 const BRep_Builder& theBuilder,
442 TopoDS_Compound& theClosed,
443 TopoDS_Compound& theOpened,
444 const Standard_Boolean theIgnore1DSubShape)
446 if (theShape.IsNull())
451 switch (theShape.ShapeType())
453 case TopAbs_COMPOUND:
454 case TopAbs_COMPSOLID:
456 for (TopoDS_Iterator anIter (theShape); anIter.More(); anIter.Next())
458 ExploreSolids (anIter.Value(), theBuilder, theClosed, theOpened, theIgnore1DSubShape);
464 for (TopoDS_Iterator anIter (theShape); anIter.More(); anIter.Next())
466 const TopoDS_Shape& aSubShape = anIter.Value();
467 const Standard_Boolean isClosed = aSubShape.ShapeType() == TopAbs_SHELL &&
468 BRep_Tool::IsClosed (aSubShape) &&
469 StdPrs_ToolTriangulatedShape::IsTriangulated (aSubShape);
470 theBuilder.Add (isClosed ? theClosed : theOpened, aSubShape);
477 theBuilder.Add (theOpened, theShape);
484 if (!theIgnore1DSubShape)
486 theBuilder.Add (theOpened, theShape);
496 // =======================================================================
499 // =======================================================================
500 void StdPrs_ShadedShape::Add (const Handle(Prs3d_Presentation)& thePrs,
501 const TopoDS_Shape& theShape,
502 const Handle(Prs3d_Drawer)& theDrawer,
503 const StdPrs_Volume theVolume)
506 StdPrs_ShadedShape::Add (thePrs, theShape, theDrawer,
507 Standard_False, aDummy, aDummy, aDummy, theVolume);
510 // =======================================================================
513 // =======================================================================
514 void StdPrs_ShadedShape::Add (const Handle (Prs3d_Presentation)& thePrs,
515 const TopoDS_Shape& theShape,
516 const Handle (Prs3d_Drawer)& theDrawer,
517 const Standard_Boolean theHasTexels,
518 const gp_Pnt2d& theUVOrigin,
519 const gp_Pnt2d& theUVRepeat,
520 const gp_Pnt2d& theUVScale,
521 const StdPrs_Volume theVolume)
523 if (theShape.IsNull())
528 // add wireframe presentation for isolated edges and vertices
529 wireframeFromShape (thePrs, theShape, theDrawer);
531 // Use automatic re-triangulation with deflection-check logic only if this feature is enable
532 if (theDrawer->IsAutoTriangulation())
534 // Triangulation completeness is important for "open-closed" analysis - perform tessellation beforehand
535 StdPrs_ToolTriangulatedShape::Tessellate (theShape, theDrawer);
538 // add special wireframe presentation for faces without triangulation
539 wireframeNoTriangFacesFromShape (thePrs, theShape, theDrawer);
541 // The shape types listed below need advanced analysis as potentially containing
542 // both closed and open parts. Solids are also included, because they might
543 // contain non-manifold parts inside (internal open shells)
544 if ((theShape.ShapeType() == TopAbs_COMPOUND
545 || theShape.ShapeType() == TopAbs_COMPSOLID
546 || theShape.ShapeType() == TopAbs_SOLID)
547 && theVolume == StdPrs_Volume_Autodetection)
549 // collect two compounds: for opened and closed (solid) sub-shapes
550 TopoDS_Compound anOpened, aClosed;
551 BRep_Builder aBuilder;
552 aBuilder.MakeCompound (aClosed);
553 aBuilder.MakeCompound (anOpened);
554 ExploreSolids (theShape, aBuilder, aClosed, anOpened, Standard_True);
556 if (aClosed.NbChildren() > 0)
558 shadeFromShape (aClosed, thePrs, theDrawer,
559 theHasTexels, theUVOrigin, theUVRepeat, theUVScale, true);
562 if (anOpened.NbChildren() > 0)
564 shadeFromShape (anOpened, thePrs, theDrawer,
565 theHasTexels, theUVOrigin, theUVRepeat, theUVScale, false);
570 // if the shape type is not compound, composolid or solid, use autodetection back-facing filled
571 shadeFromShape (theShape, thePrs, theDrawer,
572 theHasTexels, theUVOrigin, theUVRepeat, theUVScale,
573 theVolume == StdPrs_Volume_Closed);
576 if (theDrawer->FaceBoundaryDraw())
578 Handle(Graphic3d_ArrayOfSegments) aBndSegments = fillFaceBoundaries (theShape);
579 if (!aBndSegments.IsNull())
581 Handle(Graphic3d_AspectLine3d) aBoundaryAspect = theDrawer->FaceBoundaryAspect()->Aspect();
582 Handle(Graphic3d_Group) aPrsGrp = Prs3d_Root::CurrentGroup (thePrs);
583 aPrsGrp->SetGroupPrimitivesAspect (aBoundaryAspect);
584 aPrsGrp->AddPrimitiveArray (aBndSegments);
589 // =======================================================================
590 // function : FillTriangles
592 // =======================================================================
593 Handle(Graphic3d_ArrayOfTriangles) StdPrs_ShadedShape::FillTriangles (const TopoDS_Shape& theShape,
594 const Standard_Boolean theHasTexels,
595 const gp_Pnt2d& theUVOrigin,
596 const gp_Pnt2d& theUVRepeat,
597 const gp_Pnt2d& theUVScale)
599 return fillTriangles (theShape, theHasTexels, theUVOrigin, theUVRepeat, theUVScale);
602 // =======================================================================
603 // function : FillFaceBoundaries
605 // =======================================================================
606 Handle(Graphic3d_ArrayOfSegments) StdPrs_ShadedShape::FillFaceBoundaries (const TopoDS_Shape& theShape)
608 return fillFaceBoundaries (theShape);
611 // =======================================================================
612 // function : AddWireframeForFreeElements
614 // =======================================================================
615 void StdPrs_ShadedShape::AddWireframeForFreeElements (const Handle (Prs3d_Presentation)& thePrs,
616 const TopoDS_Shape& theShape,
617 const Handle (Prs3d_Drawer)& theDrawer)
619 wireframeFromShape (thePrs, theShape, theDrawer);
622 // =======================================================================
623 // function : AddWireframeForFacesWithoutTriangles
625 // =======================================================================
626 void StdPrs_ShadedShape::AddWireframeForFacesWithoutTriangles (const Handle(Prs3d_Presentation)& thePrs,
627 const TopoDS_Shape& theShape,
628 const Handle(Prs3d_Drawer)& theDrawer)
630 wireframeNoTriangFacesFromShape (thePrs, theShape, theDrawer);