#include "HYDROData_Iterator.h"
#include "HYDROData_NaturalObject.h"
-#include <TopoDS_Shape.hxx>
-
#include <TopTools_SequenceOfShape.hxx>
-
#include <TopExp_Explorer.hxx>
#include <QFile>
#include <QTextStream>
#include <limits>
-
+#define CLASS2D 1
+#ifdef CLASS2D
+#include <gp_Pnt.hxx>
+#include <gp_Vec.hxx>
+#include <gp_Pln.hxx>
+#include <TopExp.hxx>
+#include <Geom_Plane.hxx>
+#include <BRep_Tool.hxx>
+#include <BRepBuilderAPI_FindPlane.hxx>
+#include <TopoDS.hxx>
+#include <TopoDS_Vertex.hxx>
+#include <TopoDS_Wire.hxx>
+#include <TopoDS_Face.hxx>
+#include <TopoDS_Iterator.hxx>
+#include <TopTools_MapOfShape.hxx>
+#include <TColgp_SequenceOfVec.hxx>
+#include <TopTools_ShapeMapHasher.hxx>
+#undef _NCollection_MapHasher
+#endif
static int aMaxNameId = std::numeric_limits<int>::max();
void HYDROData_Tool::WriteStringsToFile( QFile& theFile,
return aName;
}
-Handle(HYDROData_Entity) HYDROData_Tool::FindObjectByName( const Handle(HYDROData_Document)& theDoc,
- const QString& theName,
- const ObjectKind theObjectKind )
-{
- Handle(HYDROData_Entity) anObject;
- if ( theName.isEmpty() || theDoc.IsNull() )
- return anObject;
-
- QStringList aNamesList;
- aNamesList << theName;
-
- HYDROData_SequenceOfObjects aSeqOfObjs = FindObjectsByNames( theDoc, aNamesList, theObjectKind );
- if( aSeqOfObjs.IsEmpty() )
- return anObject;
-
- anObject = aSeqOfObjs.First();
- return anObject;
-}
-
-HYDROData_SequenceOfObjects HYDROData_Tool::FindObjectsByNames( const Handle(HYDROData_Document)& theDoc,
- const QStringList& theNames,
- const ObjectKind theObjectKind )
-{
- HYDROData_SequenceOfObjects aResSeq;
- if( theDoc.IsNull() )
- return aResSeq;
-
- QStringList aNamesList = theNames;
-
- HYDROData_Iterator anIter( theDoc, theObjectKind );
- for( ; anIter.More(); anIter.Next() )
- {
- Handle(HYDROData_Entity) anObject = anIter.Current();
- if( anObject.IsNull() )
- continue;
-
- QString anObjName = anObject->GetName();
- if ( anObjName.isEmpty() || !aNamesList.contains( anObjName ) )
- continue;
-
- aResSeq.Append( anObject );
-
- aNamesList.removeAll( anObjName );
- if ( aNamesList.isEmpty() )
- break;
- }
-
- return aResSeq;
-}
-
bool HYDROData_Tool::IsGeometryObject( const Handle(HYDROData_Entity)& theObject )
{
if ( theObject.IsNull() )
}
return aName;
-}
\ No newline at end of file
+}
+//======================================================================================================
+// the Face to be: planer, whitout holes
+#ifdef CLASS2D
+TopAbs_State HYDROData_Tool::ComputePointState( const gp_XY& thePnt2d, const TopoDS_Face& theFace )
+{
+ TopAbs_State aState(TopAbs_UNKNOWN);
+ if(theFace.IsNull()) return aState;
+
+#ifdef CLASS2D
+ TopoDS_Wire aWire; int nb(0);
+ TopoDS_Iterator it(theFace);
+ for(;it.More();it.Next()) {
+ aWire = TopoDS::Wire(it.Value());
+ nb++;
+ }
+ if(nb > 1 || aWire.IsNull()) return aState;
+ gp_Pln aPlane;
+ gp_Vec aNormal;
+ BRepBuilderAPI_FindPlane fndPlane (theFace, Precision::Confusion());
+ if(fndPlane.Found())
+ aPlane = fndPlane.Plane()->Pln();
+ else
+ return aState;
+ aNormal = gp_Vec(aPlane.Axis().Direction());
+ if(theFace.Orientation() == TopAbs_REVERSED)
+ aNormal.Reverse();
+ gp_Pnt aPoint = gp_Pnt (thePnt2d.X(), thePnt2d.Y(), 0);
+ TColgp_SequenceOfVec aSeq;
+ TopTools_MapOfShape aMap;
+ it.Initialize(aWire);
+ for (;it.More(); it.Next()) {
+ const TopoDS_Vertex& aVx = TopExp::FirstVertex(TopoDS::Edge(it.Value()), Standard_True);
+ if(!aMap.Add(aVx))
+ continue;
+ const gp_Pnt& aCurPnt = BRep_Tool::Pnt(TopoDS::Vertex(aVx));
+ if(aPoint.IsEqual(aCurPnt,Precision::Confusion())) {
+ aState = TopAbs_ON;
+ return aState;
+ }
+ gp_Vec aVec (aPoint, aCurPnt);
+ aSeq.Append(aVec);
+ }
+ Standard_Real anAng(0.0);
+ for(int i = 1;i < aSeq.Length();i++) {
+ const gp_Vec& aV1 = aSeq.Value(i);
+ const gp_Vec& aV2 = aSeq.Value(i+1);
+ anAng += aV1.AngleWithRef(aV2, aNormal);
+ }
+ anAng += aSeq.Value(aSeq.Length()).AngleWithRef(aSeq.Value(1), aNormal);
+ if(abs(anAng) > Precision::Angular())
+ aState = TopAbs_IN;
+ else
+ aState = TopAbs_OUT;
+
+#else
+ BRepTopAdaptor_FClass2d aClassifier( theFace, Precision::Confusion() );
+ aState = aClassifier.Perform( gp_Pnt2d( thePoint ), Standard_False );
+#endif
+ return aState;
+}
+#endif
\ No newline at end of file