mirror of
https://git.dev.opencascade.org/repos/occt.git
synced 2025-04-04 18:06:22 +03:00
67441d0ca3
Prs3d_Drawer, AIS_Shape, AIS_InteractiveContext - removed properties HLRDeviationCoefficient() and HLRAngle(). Prs3d_Drawer::HLRAngle() is kept as alias to Prs3d_Drawer::DeviationAngle() with deprecated flag. Prs3d_Drawer::DeviationAngle() default value is changed from 12 to 20 degrees to match Prs3d_Drawer::HLRAngle() which has been previously used in majority of cases. Removed unused property HLRBRep_PolyAlgo::Angle().
308 lines
11 KiB
C++
308 lines
11 KiB
C++
// Created on: 1993-10-27
|
|
// Created by: Jean-LOuis FRENKEL
|
|
// Copyright (c) 1993-1999 Matra Datavision
|
|
// Copyright (c) 1999-2014 OPEN CASCADE SAS
|
|
//
|
|
// This file is part of Open CASCADE Technology software library.
|
|
//
|
|
// This library is free software; you can redistribute it and/or modify it under
|
|
// the terms of the GNU Lesser General Public License version 2.1 as published
|
|
// by the Free Software Foundation, with special exception defined in the file
|
|
// OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
|
|
// distribution for complete text of the license and disclaimer of any warranty.
|
|
//
|
|
// Alternatively, this file may be used under the terms of Open CASCADE
|
|
// commercial license or contractual agreement.
|
|
|
|
#include <StdPrs_ToolTriangulatedShape.hxx>
|
|
|
|
#include <BRepMesh_DiscretFactory.hxx>
|
|
#include <BRepMesh_DiscretRoot.hxx>
|
|
#include <BRepTools.hxx>
|
|
#include <BRep_Tool.hxx>
|
|
#include <GeomAbs_SurfaceType.hxx>
|
|
#include <GeomLib.hxx>
|
|
#include <gp_XYZ.hxx>
|
|
#include <Poly.hxx>
|
|
#include <Poly_Connect.hxx>
|
|
#include <Poly_Triangulation.hxx>
|
|
#include <Precision.hxx>
|
|
#include <Prs3d.hxx>
|
|
#include <Prs3d_Drawer.hxx>
|
|
#include <TColgp_Array1OfPnt.hxx>
|
|
#include <TColgp_Array1OfPnt2d.hxx>
|
|
#include <TopAbs_Orientation.hxx>
|
|
#include <TopLoc_Location.hxx>
|
|
#include <TShort_HArray1OfShortReal.hxx>
|
|
#include <TShort_Array1OfShortReal.hxx>
|
|
#include <TopExp_Explorer.hxx>
|
|
#include <TopoDS.hxx>
|
|
#include <TopoDS_Face.hxx>
|
|
|
|
//=======================================================================
|
|
//function : IsTriangulated
|
|
//purpose :
|
|
//=======================================================================
|
|
Standard_Boolean StdPrs_ToolTriangulatedShape::IsTriangulated (const TopoDS_Shape& theShape)
|
|
{
|
|
TopLoc_Location aLocDummy;
|
|
for (TopExp_Explorer aFaceIter (theShape, TopAbs_FACE); aFaceIter.More(); aFaceIter.Next())
|
|
{
|
|
const TopoDS_Face& aFace = TopoDS::Face (aFaceIter.Current());
|
|
const Handle(Poly_Triangulation)& aTri = BRep_Tool::Triangulation (aFace, aLocDummy);
|
|
if (aTri.IsNull())
|
|
{
|
|
return Standard_False;
|
|
}
|
|
}
|
|
return Standard_True;
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : IsClosed
|
|
//purpose :
|
|
//=======================================================================
|
|
Standard_Boolean StdPrs_ToolTriangulatedShape::IsClosed (const TopoDS_Shape& theShape)
|
|
{
|
|
if (theShape.IsNull())
|
|
{
|
|
return Standard_True;
|
|
}
|
|
|
|
switch (theShape.ShapeType())
|
|
{
|
|
case TopAbs_COMPOUND:
|
|
case TopAbs_COMPSOLID:
|
|
default:
|
|
{
|
|
// check that compound consists of closed solids
|
|
for (TopoDS_Iterator anIter (theShape); anIter.More(); anIter.Next())
|
|
{
|
|
const TopoDS_Shape& aShape = anIter.Value();
|
|
if (!IsClosed (aShape))
|
|
{
|
|
return Standard_False;
|
|
}
|
|
}
|
|
return Standard_True;
|
|
}
|
|
case TopAbs_SOLID:
|
|
{
|
|
// Check for non-manifold topology first of all:
|
|
// have to use BRep_Tool::IsClosed() because it checks the face connectivity
|
|
// inside the shape
|
|
if (!BRep_Tool::IsClosed (theShape))
|
|
return Standard_False;
|
|
|
|
for (TopoDS_Iterator anIter (theShape); anIter.More(); anIter.Next())
|
|
{
|
|
const TopoDS_Shape& aShape = anIter.Value();
|
|
if (aShape.IsNull())
|
|
{
|
|
continue;
|
|
}
|
|
|
|
if (aShape.ShapeType() == TopAbs_FACE)
|
|
{
|
|
// invalid solid
|
|
return Standard_False;
|
|
}
|
|
else if (!IsTriangulated (aShape))
|
|
{
|
|
// mesh contains holes
|
|
return Standard_False;
|
|
}
|
|
}
|
|
return Standard_True;
|
|
}
|
|
case TopAbs_SHELL:
|
|
case TopAbs_FACE:
|
|
{
|
|
// free faces / shell are not allowed
|
|
return Standard_False;
|
|
}
|
|
case TopAbs_WIRE:
|
|
case TopAbs_EDGE:
|
|
case TopAbs_VERTEX:
|
|
{
|
|
// ignore
|
|
return Standard_True;
|
|
}
|
|
}
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : ComputeNormals
|
|
//purpose :
|
|
//=======================================================================
|
|
void StdPrs_ToolTriangulatedShape::ComputeNormals (const TopoDS_Face& theFace,
|
|
const Handle(Poly_Triangulation)& theTris,
|
|
Poly_Connect& thePolyConnect)
|
|
{
|
|
if (theTris.IsNull()
|
|
|| theTris->HasNormals())
|
|
{
|
|
return;
|
|
}
|
|
|
|
// take in face the surface location
|
|
const TopoDS_Face aZeroFace = TopoDS::Face (theFace.Located (TopLoc_Location()));
|
|
Handle(Geom_Surface) aSurf = BRep_Tool::Surface (aZeroFace);
|
|
const Poly_Array1OfTriangle& aTriangles = theTris->Triangles();
|
|
if (!theTris->HasUVNodes() || aSurf.IsNull())
|
|
{
|
|
// compute normals by averaging triangulation normals sharing the same vertex
|
|
Poly::ComputeNormals (theTris);
|
|
return;
|
|
}
|
|
|
|
const Standard_Real aTol = Precision::Confusion();
|
|
Handle(TShort_HArray1OfShortReal) aNormals = new TShort_HArray1OfShortReal (1, theTris->NbNodes() * 3);
|
|
const TColgp_Array1OfPnt2d& aNodesUV = theTris->UVNodes();
|
|
Standard_Integer aTri[3];
|
|
const TColgp_Array1OfPnt& aNodes = theTris->Nodes();
|
|
gp_Dir aNorm;
|
|
for (Standard_Integer aNodeIter = aNodes.Lower(); aNodeIter <= aNodes.Upper(); ++aNodeIter)
|
|
{
|
|
// try to retrieve normal from real surface first, when UV coordinates are available
|
|
if (GeomLib::NormEstim (aSurf, aNodesUV.Value (aNodeIter), aTol, aNorm) > 1)
|
|
{
|
|
if (thePolyConnect.Triangulation() != theTris)
|
|
{
|
|
thePolyConnect.Load (theTris);
|
|
}
|
|
|
|
// compute flat normals
|
|
gp_XYZ eqPlan (0.0, 0.0, 0.0);
|
|
for (thePolyConnect.Initialize (aNodeIter); thePolyConnect.More(); thePolyConnect.Next())
|
|
{
|
|
aTriangles (thePolyConnect.Value()).Get (aTri[0], aTri[1], aTri[2]);
|
|
const gp_XYZ v1 (aNodes (aTri[1]).Coord() - aNodes (aTri[0]).Coord());
|
|
const gp_XYZ v2 (aNodes (aTri[2]).Coord() - aNodes (aTri[1]).Coord());
|
|
const gp_XYZ vv = v1 ^ v2;
|
|
const Standard_Real aMod = vv.Modulus();
|
|
if (aMod >= aTol)
|
|
{
|
|
eqPlan += vv / aMod;
|
|
}
|
|
}
|
|
const Standard_Real aModMax = eqPlan.Modulus();
|
|
aNorm = (aModMax > aTol) ? gp_Dir (eqPlan) : gp::DZ();
|
|
}
|
|
|
|
const Standard_Integer anId = (aNodeIter - aNodes.Lower()) * 3;
|
|
aNormals->SetValue (anId + 1, (Standard_ShortReal )aNorm.X());
|
|
aNormals->SetValue (anId + 2, (Standard_ShortReal )aNorm.Y());
|
|
aNormals->SetValue (anId + 3, (Standard_ShortReal )aNorm.Z());
|
|
}
|
|
theTris->SetNormals (aNormals);
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : Normal
|
|
//purpose :
|
|
//=======================================================================
|
|
void StdPrs_ToolTriangulatedShape::Normal (const TopoDS_Face& theFace,
|
|
Poly_Connect& thePolyConnect,
|
|
TColgp_Array1OfDir& theNormals)
|
|
{
|
|
const Handle(Poly_Triangulation)& aPolyTri = thePolyConnect.Triangulation();
|
|
if (!aPolyTri->HasNormals())
|
|
{
|
|
ComputeNormals (theFace, aPolyTri, thePolyConnect);
|
|
}
|
|
|
|
const TColgp_Array1OfPnt& aNodes = aPolyTri->Nodes();
|
|
const TShort_Array1OfShortReal& aNormals = aPolyTri->Normals();
|
|
const Standard_ShortReal* aNormArr = &aNormals.First();
|
|
for (Standard_Integer aNodeIter = aNodes.Lower(); aNodeIter <= aNodes.Upper(); ++aNodeIter)
|
|
{
|
|
const Standard_Integer anId = 3 * (aNodeIter - aNodes.Lower());
|
|
const gp_Dir aNorm (aNormArr[anId + 0],
|
|
aNormArr[anId + 1],
|
|
aNormArr[anId + 2]);
|
|
theNormals (aNodeIter) = aNorm;
|
|
}
|
|
|
|
if (theFace.Orientation() == TopAbs_REVERSED)
|
|
{
|
|
for (Standard_Integer aNodeIter = aNodes.Lower(); aNodeIter <= aNodes.Upper(); ++aNodeIter)
|
|
{
|
|
theNormals.ChangeValue (aNodeIter).Reverse();
|
|
}
|
|
}
|
|
}
|
|
|
|
//=======================================================================
|
|
//function : IsTessellated
|
|
//purpose :
|
|
//=======================================================================
|
|
Standard_Boolean StdPrs_ToolTriangulatedShape::IsTessellated (const TopoDS_Shape& theShape,
|
|
const Handle(Prs3d_Drawer)& theDrawer)
|
|
{
|
|
return BRepTools::Triangulation (theShape, Prs3d::GetDeflection (theShape, theDrawer), true);
|
|
}
|
|
|
|
// =======================================================================
|
|
// function : Tessellate
|
|
// purpose :
|
|
// =======================================================================
|
|
Standard_Boolean StdPrs_ToolTriangulatedShape::Tessellate (const TopoDS_Shape& theShape,
|
|
const Handle(Prs3d_Drawer)& theDrawer)
|
|
{
|
|
Standard_Boolean wasRecomputed = Standard_False;
|
|
// Check if it is possible to avoid unnecessary recomputation of shape triangulation
|
|
if (IsTessellated (theShape, theDrawer))
|
|
{
|
|
return wasRecomputed;
|
|
}
|
|
|
|
Standard_Real aDeflection = Prs3d::GetDeflection (theShape, theDrawer);
|
|
|
|
// retrieve meshing tool from Factory
|
|
Handle(BRepMesh_DiscretRoot) aMeshAlgo = BRepMesh_DiscretFactory::Get().Discret (theShape,
|
|
aDeflection,
|
|
theDrawer->DeviationAngle());
|
|
if (!aMeshAlgo.IsNull())
|
|
{
|
|
aMeshAlgo->Perform();
|
|
wasRecomputed = Standard_True;
|
|
}
|
|
|
|
return wasRecomputed;
|
|
}
|
|
|
|
// =======================================================================
|
|
// function : ClearOnOwnDeflectionChange
|
|
// purpose :
|
|
// =======================================================================
|
|
void StdPrs_ToolTriangulatedShape::ClearOnOwnDeflectionChange (const TopoDS_Shape& theShape,
|
|
const Handle(Prs3d_Drawer)& theDrawer,
|
|
const Standard_Boolean theToResetCoeff)
|
|
{
|
|
if (!theDrawer->IsAutoTriangulation()
|
|
|| theShape.IsNull())
|
|
{
|
|
return;
|
|
}
|
|
|
|
const Standard_Boolean isOwnDeviationAngle = theDrawer->HasOwnDeviationAngle();
|
|
const Standard_Boolean isOwnDeviationCoefficient = theDrawer->HasOwnDeviationCoefficient();
|
|
const Standard_Real anAngleNew = theDrawer->DeviationAngle();
|
|
const Standard_Real anAnglePrev = theDrawer->PreviousDeviationAngle();
|
|
const Standard_Real aCoeffNew = theDrawer->DeviationCoefficient();
|
|
const Standard_Real aCoeffPrev = theDrawer->PreviousDeviationCoefficient();
|
|
if ((!isOwnDeviationAngle || Abs (anAngleNew - anAnglePrev) <= Precision::Angular())
|
|
&& (!isOwnDeviationCoefficient || Abs (aCoeffNew - aCoeffPrev) <= Precision::Confusion()))
|
|
{
|
|
return;
|
|
}
|
|
|
|
BRepTools::Clean (theShape);
|
|
if (theToResetCoeff)
|
|
{
|
|
theDrawer->UpdatePreviousDeviationAngle();
|
|
theDrawer->UpdatePreviousDeviationCoefficient();
|
|
}
|
|
}
|