]> SALOME platform Git repositories - modules/hydro.git/commitdiff
Salome HOME
26.12.2013. Update of SplitToZonesTool + new Transform files.
authorszy <szy@opencascade.com>
Thu, 26 Dec 2013 13:25:07 +0000 (13:25 +0000)
committerszy <szy@opencascade.com>
Thu, 26 Dec 2013 13:25:07 +0000 (13:25 +0000)
src/HYDROData/CMakeLists.txt
src/HYDROData/HYDROData_SplitToZonesTool.cxx
src/HYDROData/HYDROData_Transform.cxx [new file with mode: 0644]
src/HYDROData/HYDROData_Transform.h [new file with mode: 0644]

index 4ef25b8b9f1d86e8e9fa56078bb8b723b7c181ab..df55dbc0343d5429724968d3cab6cf29ba105088 100644 (file)
@@ -40,6 +40,7 @@ set(PROJECT_HEADERS
     HYDROData_Stream.h
     HYDROData_StreamAltitude.h
     HYDROData_Tool.h
+    HYDROData_Transform.h
     HYDROData_VisualState.h
     HYDROData_Zone.h
 )
@@ -83,6 +84,7 @@ set(PROJECT_SOURCES
     HYDROData_Stream.cxx
     HYDROData_StreamAltitude.cxx
     HYDROData_Tool.cxx
+    HYDROData_Transform.cxx
     HYDROData_VisualState.cxx
     HYDROData_Zone.cxx
 )
index 04e212a7f2deb0e7ee07a0fd4c676192a6930317..a4312ab4671110e8e16b579f9c26ffe467621b7f 100644 (file)
@@ -3,32 +3,43 @@
 
 #include "HYDROData_PolylineXY.h"
 #include "HYDROData_ShapesGroup.h"
-
-#include <BRepAlgoAPI_Common.hxx>
+#include <HYDROData_Transform.h>
 #include <BRepAlgoAPI_Cut.hxx>
+#include <BRepAlgoAPI_Common.hxx>
 #include <BRepBuilderAPI_MakeFace.hxx>
-
+#include <BRep_Builder.hxx>
 #include <TopExp_Explorer.hxx>
 #include <TopoDS.hxx>
-#include <TopoDS_Edge.hxx>
-#include <TopoDS_Iterator.hxx>
+#include <TopoDS_Shape.hxx>
+#include <TopoDS_Compound.hxx>
 #include <TopoDS_Wire.hxx>
+#include <TopoDS_Edge.hxx>
 #include <TopoDS_Iterator.hxx>
 #include <BRepCheck_Analyzer.hxx>
 #include <NCollection_IncAllocator.hxx>
 #include <BOPAlgo_Builder.hxx>
 #include <BOPAlgo_PaveFiller.hxx>
+#include <NCollection_Map.hxx>
 #include <NCollection_DataMap.hxx>
 #include <TopTools_ListOfShape.hxx>
 #include <TopTools_ListIteratorOfListOfShape.hxx>
 #include <TopTools_ShapeMapHasher.hxx>
+#include <gp_Pln.hxx>
+#include <BRepGProp.hxx>
+#include <GProp_GProps.hxx>
+#include <Geom_Plane.hxx>
+#include <BRepBuilderAPI_FindPlane.hxx>
+
+
+typedef NCollection_Map<TopoDS_Shape, TopTools_ShapeMapHasher> HYDROData_MapOfShape;
+typedef HYDROData_MapOfShape::Iterator HYDROData_MapIteratorOfMapOfShape;
 typedef NCollection_DataMap<TopoDS_Shape, TopTools_ListOfShape, TopTools_ShapeMapHasher> HYDROData_DataMapOfShapeListOfShape;
 typedef HYDROData_DataMapOfShapeListOfShape::Iterator HYDROData_DataMapIteratorOfDataMapOfShapeListOfShape;
 typedef NCollection_DataMap<TopoDS_Shape,  QStringList, TopTools_ShapeMapHasher> HYDROData_DataMapOfShapeListOfString;
 typedef HYDROData_DataMapOfShapeListOfString::Iterator HYDROData_DataMapIteratorOfDataMapOfShapeListOfString;
 #undef _NCollection_MapHasher
 
-#define DEB_SPLIT_TO_ZONES 1
+//#define DEB_SPLIT_TO_ZONES 1
 #ifdef DEB_SPLIT_TO_ZONES
 #include <BRepTools.hxx>
 #endif
@@ -54,6 +65,29 @@ TopoDS_Face HYDROData_SplitToZonesTool::SplitData::Face() const
   return aResFace;
 }
 
+Standard_Boolean buildLimFace(const TopoDS_Wire& theBndWire, TopoDS_Face& outFace) 
+{
+  GProp_GProps G;
+  BRepGProp::LinearProperties(theBndWire,G);
+  const gp_Pnt& aCPnt = G.CentreOfMass();
+  gp_Pln aPln;
+  BRepBuilderAPI_FindPlane fndPlane (theBndWire, Precision::Confusion());  
+  if(fndPlane.Found())
+    aPln = fndPlane.Plane()->Pln();
+  else
+    aPln = gp_Pln(aCPnt, gp::OZ().Direction());
+  BRepBuilderAPI_MakeFace  aMkFace(aPln, theBndWire);
+  if(aMkFace.IsDone()) {
+    outFace = aMkFace.Face();    
+    if(!outFace.IsNull()) {
+#ifdef DEB_SPLIT_TO_ZONES
+//       BRepTools::Write(limFace,"FL.brep");
+#endif
+      return Standard_True;
+       }
+  }
+  return Standard_False;
+}
 HYDROData_SplitToZonesTool::SplitDataList
   HYDROData_SplitToZonesTool::Split( const HYDROData_SequenceOfObjects&  theObjectList,
                                      const HYDROData_SequenceOfObjects&  theGroupsList,
@@ -64,7 +98,7 @@ HYDROData_SplitToZonesTool::SplitDataList
   // Preparation. 
   // Collect the object shapes to split. InputDataList will contain elements which will hold shape & name_of_shape.
   SplitDataList anInputSplitDataList;
-  //TCollection_AsciiString aNam("Faces_");
+  TCollection_AsciiString aNam("Faces_");
   for( int anIndex = 1, aLength = theObjectList.Length(); anIndex <= aLength; anIndex++ )
   {
     Handle(HYDROData_Object) aGeomObj = 
@@ -75,8 +109,8 @@ HYDROData_SplitToZonesTool::SplitDataList
     TopoDS_Shape aShape = aGeomObj->GetTopShape();
     if ( aShape.IsNull() )
       continue;  
-//TCollection_AsciiString aName = aNam + anIndex + ".brep";
-//BRepTools::Write(aShape, aName.ToCString());
+   //TCollection_AsciiString aName = aNam + anIndex + ".brep";
+   //BRepTools::Write(aShape, aName.ToCString());
 
     if ( aShape.ShapeType() == TopAbs_COMPOUND ) {
       // Create split data for each face contained in the compound
@@ -132,7 +166,6 @@ HYDROData_SplitToZonesTool::SplitDataList
          anOutputSplitDataList.append(anInputGroupList);       
     return anOutputSplitDataList;
   }
   HYDROData_DataMapOfShapeListOfString aDM3;
   if(!anInputGroupList.isEmpty()) {// Old edge ==> List_Of_Names
     QStringList aListOfNames;
@@ -146,15 +179,38 @@ HYDROData_SplitToZonesTool::SplitDataList
   BOPAlgo_PaveFiller* aPaveFiller      = new BOPAlgo_PaveFiller(pA1);
   BOPAlgo_Builder* aBuilder            = new BOPAlgo_Builder(pA2);
   BOPCol_ListOfShape aLS;  
+
   QStringList aListOfNames;
+  TopoDS_Compound aCmp;
+  BRep_Builder aBB;
+  aBB.MakeCompound(aCmp);
   for (int i=0;i < anInputSplitDataList.size() ;i++) {
     const TopoDS_Shape& aSh = anInputSplitDataList.at(i).Shape;
     aDM3.Bind(aSh, anInputSplitDataList.at(i).ObjectNames);
     aLS.Append(aSh);
+       aBB.Add(aCmp,aSh);
+    //TCollection_AsciiString aName = aNam + i + ".brep";
+    //BRepTools::Write(aSh, aName.ToCString());
   }
+  HYDROData_Transform sTool;
+  sTool.Clear();
+  sTool.SetArgument(aCmp);
+  sTool.Detect();
+  const Standard_Boolean bToTransform = sTool.ToTransform();
+  Standard_Integer anErr(0);  
+  sTool.SetToTransform(bToTransform);  
+  sTool.Perform();
+  anErr = sTool.ErrorStatus();
+  if(anErr)
+    return anOutputSplitDataList;
+  const TopoDS_Shape& aResult = sTool.Shape();
+  /* 
+  else {  
+  BOP ==
   aPaveFiller->SetArguments(aLS);
   aPaveFiller->Perform();
-  Standard_Integer anErr = aPaveFiller->ErrorStatus();
+  anErr = aPaveFiller->ErrorStatus();
   if(anErr) 
     return anOutputSplitDataList;
   BOPDS_PDS pDS = aPaveFiller->PDS();
@@ -163,16 +219,20 @@ HYDROData_SplitToZonesTool::SplitDataList
   aBuilder->Clear();
 
   // Step 2. Split faces 
-  BOPCol_ListIteratorOfListOfShape anIt(aPaveFiller->Arguments());
+  BOPCol_ListIteratorOfListOfShape anIt(aLS);
   for (; anIt.More(); anIt.Next()) {
     const TopoDS_Shape& aS = anIt.Value();
     aBuilder->AddArgument(aS);
   }
   aBuilder->PerformWithFiller(*aPaveFiller);
   anErr = aBuilder->ErrorStatus();
+
+  ////
   if(anErr) 
     return anOutputSplitDataList;
-  const TopoDS_Shape& aResult = aBuilder->Shape();
+  aResult = aBuilder->Shape();
+  } */
+
   if (aResult.IsNull()) 
     return anOutputSplitDataList;
   BRepCheck_Analyzer aCheck (aResult);
@@ -184,17 +244,17 @@ HYDROData_SplitToZonesTool::SplitDataList
     return anOutputSplitDataList;
   }
 #ifdef DEB_SPLIT_TO_ZONES
-  //BRepTools::Write(aResult, "SplitFacesV.brep");
+  BRepTools::Write(aResult, "SplitFacesV.brep");
 #endif
   
   // Step 3. Collect history  
   HYDROData_DataMapOfShapeListOfShape aDM1;
-  anIt.Init(aLS); 
+  BOPCol_ListIteratorOfListOfShape anIt(aLS); 
   //TCollection_AsciiString aNamM ("EdgM_");
   //TCollection_AsciiString aNamG ("EdgG_");
   for (int i =1;anIt.More();anIt.Next(),i++) {
        Standard_Boolean foundF(Standard_False);
-    const TopTools_ListOfShape& aListOfNew = aBuilder->Modified(anIt.Value());
+    const TopTools_ListOfShape& aListOfNew = /*aBuilder->*/sTool.Modified(anIt.Value());
        if(!aListOfNew.IsEmpty())
          foundF = Standard_True;
     TopTools_ListOfShape aList;
@@ -202,7 +262,7 @@ HYDROData_SplitToZonesTool::SplitDataList
     for(;it.More();it.Next())       
       aList.Append(it.Value());
     // Bug in History: partition should give only modified entities! => temporary solution is used
-    const TopTools_ListOfShape& aListOfGen = aBuilder->Generated(anIt.Value());
+    const TopTools_ListOfShape& aListOfGen = /*aBuilder->*/sTool.Generated(anIt.Value());
        if(!aListOfGen.IsEmpty())
          foundF = Standard_True;
        it.Initialize(aListOfGen);    
@@ -217,7 +277,7 @@ HYDROData_SplitToZonesTool::SplitDataList
          TopExp_Explorer exp (anIt.Value(), TopAbs_EDGE);
          for (int j =1;exp.More();exp.Next(),j++) {
            Standard_Boolean foundE(Standard_False);
-           const TopTools_ListOfShape& aListM = aBuilder->Modified(exp.Current());     
+           const TopTools_ListOfShape& aListM = /*aBuilder->*/sTool.Modified(exp.Current());   
                //cout << "NB_EDGE_M = " << aListM.Extent() <<endl;
                if(aListM.Extent()) foundE = Standard_True;
                it.Initialize(aListM);    
@@ -226,7 +286,7 @@ HYDROData_SplitToZonesTool::SplitDataList
                  //aName = aNamM + i + j +k +".brep";
                  //BRepTools::Write(it.Value(),aName.ToCString());
                }
-               const TopTools_ListOfShape& aListG = aBuilder->Generated(exp.Current());
+               const TopTools_ListOfShape& aListG = /*aBuilder->*/sTool.Generated(exp.Current());
                if(aListG.Extent()) foundE = Standard_True;
                it.Initialize(aListG);    
                for(int k=1;it.More();it.Next(),k++)   {  
@@ -247,25 +307,159 @@ HYDROData_SplitToZonesTool::SplitDataList
 
   // aDM2: NewShape ==> ListOfOldShapes
   HYDROData_DataMapOfShapeListOfShape aDM2;
-  HYDROData_DataMapIteratorOfDataMapOfShapeListOfShape aMIt(aDM1);
-  for(;aMIt.More();aMIt.Next()) {
-    const TopoDS_Shape& aKey = aMIt.Key();
-    TopTools_ListOfShape aList;
-    aList.Append(aKey);
-    const TopTools_ListOfShape& aListOfNew = aMIt.Value();
-    TopTools_ListIteratorOfListOfShape it(aListOfNew);
-    for(;it.More();it.Next()) {
-      if(!aDM2.IsBound(it.Value()))
-        aDM2.Bind(it.Value(), aList);
-      else {
-        TopTools_ListOfShape& aList = aDM2.ChangeFind(it.Value());
-        aList.Prepend(aKey);
-      }
-    }
+  // make limiting face
+  HYDROData_DataMapOfShapeListOfShape aDM4;
+  Standard_Boolean hasLimits(Standard_False);
+  if (! thePolyline.IsNull()) {         
+    //const TopoDS_Shape aShape = thePolyline->GetShape();     
+    const TopoDS_Wire aBndWire = TopoDS::Wire(thePolyline->GetShape());
+    if(!aBndWire.IsNull()) {
+      TopoDS_Face limFace;       
+         if(buildLimFace(aBndWire, limFace)) {
+               TopoDS_Shape aComResult;
+           BRepAlgoAPI_Common mkCom(aResult, limFace);
+               if(mkCom.IsDone()) {
+             aComResult = mkCom.Shape();
+                 BRepCheck_Analyzer aCheck (aComResult);
+          if(aCheck.IsValid()) {
+#ifdef DEB_SPLIT_TO_ZONES
+                   BRepTools::Write(aComResult,"CommonV.brep");
+#endif 
+
+                   hasLimits = Standard_True; // DM2 should filled here
+                   HYDROData_MapOfShape aView;
+                   TopExp_Explorer exp (aResult, TopAbs_FACE);
+                   for (int i =1;exp.More();exp.Next(),i++) {
+              const TopoDS_Shape& aFace = exp.Current();
+                     if(!aFace.IsNull()) {
+                       const TopTools_ListOfShape& aListOfNew = mkCom.Modified(aFace);
+                         //cout << "Modified: " << aListOfNew.Extent() <<endl;                 
+                const TopTools_ListOfShape& aListOfGen = mkCom.Generated(aFace);                       
+                           if(!aListOfNew.IsEmpty()) {
+                     aDM4.Bind(aFace, aListOfNew);
+                             //TCollection_AsciiString aName = aNam + i + ".brep";
+                                 //BRepTools::Write(aListOfNew.Last(), aName.ToCString());                             
+                               }
+                               else {
+                             if(!mkCom.IsDeleted(aFace)) {
+                                   if(!aListOfGen.IsEmpty()) {
+                         aDM4.Bind(aFace, aListOfGen);                        
+                                     //TCollection_AsciiString aName = aNam + i + "g.brep";
+                                     //BRepTools::Write(aListOfGen.Last(), aName.ToCString());
+                                       }
+                                        else {
+                                   TopTools_ListOfShape aList; 
+                                       aList.Append(aFace);
+                                       aDM4.Bind(aFace, aList); //the same face - not modified
+                                        }
+                                 }
+                        }
+                        TopExp_Explorer exp2 (aFace, TopAbs_EDGE);
+                        for (int j =1;exp2.More();exp2.Next(),j++) {
+                          const TopoDS_Shape& anEdge = exp2.Current();
+                          if(!anEdge.IsNull()) {
+                                if(aView.Contains(anEdge)) continue;
+                                aView.Add(anEdge);
+                                const TopTools_ListOfShape& aListOfNewEd = mkCom.Modified(anEdge);
+                                const TopTools_ListOfShape& aListOfGenEd = mkCom.Generated(anEdge);
+                 if(!aListOfNewEd.IsEmpty())
+                      aDM4.Bind(anEdge, aListOfNewEd);
+                                else {
+                                  if(!mkCom.IsDeleted(anEdge)) {
+                                        if(!aListOfGenEd.IsEmpty()) {
+                                          aDM4.Bind(anEdge, aListOfGenEd);
+                                        } else {
+                                          TopTools_ListOfShape aList; 
+                                      aList.Append(anEdge);
+                                      aDM4.Bind(anEdge, aList);//the same edge - not modified
+                                        }
+                                  }
+                                }
+                          }
+                        }
+                       }
+                 } //end DM4 filling (phase 1)
+                   //cout << "DM4 Ext = " <<aDM4.Extent() <<endl;
+                // phase 2 (from tool)
+         TopExp_Explorer expt (limFace, TopAbs_EDGE);
+                for(;expt.More();expt.Next()) {
+                  const TopoDS_Shape& anEdge = expt.Current();
+                  if(!anEdge.IsNull()) {
+             const TopTools_ListOfShape& aListOfNewEd = mkCom.Modified(anEdge);
+                        const TopTools_ListOfShape& aListOfGenEd = mkCom.Generated(anEdge);
+             if(!aListOfNewEd.IsEmpty())
+                  aDM4.Bind(anEdge, aListOfNewEd);
+                        else {
+                          if(!mkCom.IsDeleted(anEdge)) {
+                                if(!aListOfGenEd.IsEmpty()) {
+                                  aDM4.Bind(anEdge, aListOfGenEd);
+                                } else {
+                                  TopTools_ListOfShape aList; 
+                              aList.Append(anEdge);
+                              aDM4.Bind(anEdge, aList);//the same edge - not modified
+                                }
+                          }
+                        }
+                  }
+                }
+                //cout << "DM4 Ext = " <<aDM4.Extent() <<endl;
+         HYDROData_DataMapIteratorOfDataMapOfShapeListOfShape aMIt(aDM1); 
+                //DM4 contains Old - New after common op. DM1: old - new after Split op.
+         for(;aMIt.More();aMIt.Next()) {
+           const TopoDS_Shape& aKey = aMIt.Key();
+           TopTools_ListOfShape aList;
+           aList.Append(aKey);
+           const TopTools_ListOfShape& aListOfNew = aMIt.Value();
+           TopTools_ListIteratorOfListOfShape it(aListOfNew);
+           for(;it.More();it.Next()) {// iterate new: b1, b2, b3...
+                        if(!aDM4.IsBound(it.Value())) // bi - is deleted
+               continue; // go to the next bi
+                        else {
+                               const TopTools_ListOfShape& aListOfNew4 = aDM4.Find(it.Value());
+                           TopTools_ListIteratorOfListOfShape it4(aListOfNew4); // {c1, c2, c3,...}
+                for(;it4.More();it4.Next()) {
+                                 if(!aDM2.IsBound(it4.Value()))
+                    aDM2.Bind(it4.Value(), aList);
+                             else {
+                    TopTools_ListOfShape& aList = aDM2.ChangeFind(it4.Value());
+                    aList.Prepend(aKey);
+                                 }
+                               }
+                        }
+                  }//
+                 }
+                } else {
+                        hasLimits = Standard_False;
+#ifdef DEB_SPLIT_TO_ZONES
+                    BRepTools::Write(aComResult,"CommonNV.brep");
+#endif 
+                 }
+          }
+         }
+       }
+  }
+  if(!hasLimits) {
+    HYDROData_DataMapIteratorOfDataMapOfShapeListOfShape aMIt(aDM1);
+    for(;aMIt.More();aMIt.Next()) {
+      const TopoDS_Shape& aKey = aMIt.Key();
+      TopTools_ListOfShape aList;
+      aList.Append(aKey);
+      const TopTools_ListOfShape& aListOfNew = aMIt.Value();
+      TopTools_ListIteratorOfListOfShape it(aListOfNew);
+      for(;it.More();it.Next()) {
+        if(!aDM2.IsBound(it.Value()))
+          aDM2.Bind(it.Value(), aList);
+        else {
+          TopTools_ListOfShape& aList = aDM2.ChangeFind(it.Value());
+          aList.Prepend(aKey);
+               }
+         }
+       }
   }
-  //cout << "DM2 Ext = " <<aDM2.Extent() <<endl;
   // Step 4. Fill output structure.
-  aMIt.Initialize(aDM2);
+  //TCollection_AsciiString aNam4 ("SC_");
+  
+  HYDROData_DataMapIteratorOfDataMapOfShapeListOfShape aMIt(aDM2);
   for(int i =1;aMIt.More();aMIt.Next(),i++) {
     SplitData aDestSplitData;
     const TopoDS_Shape& aKey = aMIt.Key(); //new
@@ -274,6 +468,11 @@ HYDROData_SplitToZonesTool::SplitDataList
       aDestSplitData.Type = SplitData::Data_Zone;
     else
       aDestSplitData.Type = SplitData::Data_Edge;
+#ifdef DEB_SPLIT_TO_ZONES
+         //TCollection_AsciiString aName = aNam4 + i + ".brep";
+         //BRepTools::Write(aKey,aName.ToCString());
+#endif 
+
 
     QStringList aListOfNames; // names processing
     const TopTools_ListOfShape& aListOfOld = aMIt.Value();
@@ -289,20 +488,8 @@ HYDROData_SplitToZonesTool::SplitDataList
        aDestSplitData.ObjectNames = aListOfNames;  
        anOutputSplitDataList.append(aDestSplitData);
   }
-  //cout << "anOutputSplitDataList = " <<anOutputSplitDataList.size() <<endl;
   return anOutputSplitDataList;
-
-/* 
-  // Step 3. Extract the separate regions.
-  SplitDataList anExtractedSplitDataList;
-  SplitDataListIterator anOutputIter( anOutputSplitDataList );
-  while( anOutputIter.hasNext() )
-  {
-    const SplitData& anOutputSplitData = anOutputIter.next();
-    anExtractedSplitDataList.append( ExtractSeparateData( anOutputSplitData ) );
-  }
-*/
- // return anExtractedSplitDataList;
 }
 
 bool HYDROData_SplitToZonesTool::SplitTwoData( const SplitData& theData1,
diff --git a/src/HYDROData/HYDROData_Transform.cxx b/src/HYDROData/HYDROData_Transform.cxx
new file mode 100644 (file)
index 0000000..a4a2f12
--- /dev/null
@@ -0,0 +1,403 @@
+#include <HYDROData_Transform.h>
+
+#include <math.h>
+
+#include <Bnd_Box.hxx>
+#include <gp_Vec.hxx>
+
+#include <TopoDS_Shape.hxx>
+#include <TopoDS_Iterator.hxx>
+
+#include <BRepBndLib.hxx>
+
+#include <TopTools_ListIteratorOfListOfShape.hxx>
+
+#include <BOPCol_ListOfShape.hxx>
+#include <BOPCol_MapOfShape.hxx>
+
+#include <BOPCol_ListOfShape.hxx>
+
+#include <BOPAlgo_BuilderShape.hxx>
+#include <BOPAlgo_PaveFiller.hxx>
+
+//=======================================================================
+//function : HYDROData_Transform
+//purpose  : 
+//=======================================================================
+HYDROData_Transform::HYDROData_Transform()
+:
+  BOPAlgo_BuilderShape()
+{
+  myToTransform=Standard_False;
+  myTreshold=10000.;
+  myTransform1=NULL;
+  myTransform2=NULL;
+  myBuilder=NULL;
+}
+//=======================================================================
+//function : ~
+//purpose  : 
+//=======================================================================
+HYDROData_Transform::~HYDROData_Transform()
+{
+  Clear();
+}
+//=======================================================================
+//function : Clear
+//purpose  : 
+//=======================================================================
+void HYDROData_Transform::Clear()
+{
+  myTreshold=10000.;
+  //
+  if (myTransform1) {
+    delete myTransform1;
+    myTransform1=NULL;
+  }
+  if (myTransform2) {
+    delete myTransform2;
+    myTransform2=NULL;
+  }
+  if (myBuilder) {
+    delete myBuilder;
+    myBuilder=NULL;
+  }
+  //
+  BOPAlgo_BuilderShape::PrepareHistory();
+}
+//=======================================================================
+//function : SetArgument
+//purpose  : 
+//=======================================================================
+void HYDROData_Transform::SetArgument(const TopoDS_Shape& theShape)
+{
+  myArgument=theShape;
+}
+//=======================================================================
+//function : AddArgument
+//purpose  : 
+//=======================================================================
+const TopoDS_Shape& HYDROData_Transform::Argument()const
+{
+  return myArgument;
+}
+//=======================================================================
+//function : SetTreshold
+//purpose  : 
+//=======================================================================
+void HYDROData_Transform::SetTreshold(const Standard_Real theTreshold)
+{
+  myTreshold=theTreshold;
+}
+//=======================================================================
+//function : Treshold
+//purpose  : 
+//=======================================================================
+Standard_Real HYDROData_Transform::Treshold()const 
+{
+  return myTreshold;
+}
+//=======================================================================
+//function : SetToTransform
+//purpose  : 
+//=======================================================================
+void HYDROData_Transform::SetToTransform(const Standard_Boolean theFlag)
+{
+  myToTransform=theFlag;
+}
+//=======================================================================
+//function : ToTransform
+//purpose  : 
+//=======================================================================
+Standard_Boolean HYDROData_Transform::ToTransform()const
+{
+  return myToTransform;
+}
+//=======================================================================
+//function : CheckData
+//purpose  : 
+//=======================================================================
+void HYDROData_Transform::CheckData()
+{
+  if(myArgument.IsNull()) {
+    myErrorStatus=10;
+  }
+}
+//=======================================================================
+//function : Prepare
+//purpose  : 
+//=======================================================================
+void HYDROData_Transform::Prepare()
+{
+  TopoDS_Shape aSx;
+  //
+  myErrorStatus=0;
+  myWarningStatus=0;
+  myShape=aSx;
+}
+//=======================================================================
+//function : Perform
+//purpose  : 
+//=======================================================================
+void HYDROData_Transform::Perform()
+{
+  Prepare();
+  //
+  CheckData();
+  if(myErrorStatus) {
+    return;
+  }
+  //
+  Clear();
+  if(myErrorStatus) {
+    return;
+  }
+  //
+  ComputeTrsf();
+  if(myErrorStatus) {
+    return;
+  }
+  //
+  Transform1();
+  if(myErrorStatus) {
+    return;
+  }
+  //
+  BuildFuse();
+  if(myErrorStatus) {
+    return;
+  }
+  //
+  Transform2();
+}
+//=======================================================================
+//function : Detect
+//purpose  : 
+//=======================================================================
+void HYDROData_Transform::Detect()
+{
+  Standard_Real aXmin, aYmin, aZmin, aXmax, aYmax, aZmax;
+  Standard_Real aXc, aYc, aZc;
+  Bnd_Box aBox; 
+  gp_Vec aVec;
+  //
+  myErrorStatus=0;
+  //
+  BRepBndLib::Add(myArgument, aBox);
+  //
+  aBox.Get(aXmin, aYmin, aZmin, aXmax, aYmax, aZmax);
+  //
+  aXc=0.5*(aXmin+aXmax);
+  aYc=0.5*(aYmin+aYmax);
+  aZc=0.5*(aZmin+aZmax);
+  printf(" aVc : { %lf %lf %lf }\n", aXc, aYc, aZc);
+  //
+  myToTransform=Standard_False;
+  if (fabs(aXc)>myTreshold || 
+      fabs(aYc)>myTreshold || 
+      fabs(aZc)>myTreshold) {
+    myToTransform=!myToTransform;
+  }  
+  //
+  aVec.SetCoord(-aXc, -aYc, -aZc);
+  myTrsf1.SetTranslation(aVec);
+  //
+  aVec.SetCoord(aXc, aYc, aZc);
+  myTrsf2.SetTranslation(aVec);
+}
+//=======================================================================
+//function : ComputeTrsf
+//purpose  : 
+//=======================================================================
+void HYDROData_Transform::ComputeTrsf()
+{
+  if (!myToTransform) {
+    gp_Trsf aTrsf;
+    //
+    myTrsf1=aTrsf;
+    myTrsf2=aTrsf;
+  }
+}
+//=======================================================================
+//function : Transform1
+//purpose  : 
+//=======================================================================
+void HYDROData_Transform::Transform1()
+{
+  Standard_Boolean bIsDone;
+  //
+  myErrorStatus=0;
+  //
+  myTransform1=new BRepBuilderAPI_Transform(myTrsf1);
+  //
+  myTransform1->Perform(myArgument);
+  bIsDone=myTransform1->IsDone();
+  if (!bIsDone) {
+    myErrorStatus=11;
+    return;
+  }
+}
+//=======================================================================
+//function : BuildFuse
+//purpose  : 
+//=======================================================================
+void HYDROData_Transform::BuildFuse()
+{
+  Standard_Integer iErr;
+  TopoDS_Iterator aIt;
+  BOPCol_ListOfShape aLC;
+  BOPAlgo_PaveFiller aPF;
+  //
+  myErrorStatus=0;
+  //
+  const TopoDS_Shape& aS1=myTransform1->Shape();
+  //
+  aIt.Initialize(aS1);
+  for (; aIt.More(); aIt.Next()) {
+    const TopoDS_Shape& aS1x=aIt.Value();
+    aLC.Append(aS1x);
+  }
+  //
+  aPF.SetArguments(aLC);
+  //
+  aPF.Perform();
+  iErr=aPF.ErrorStatus();
+  if (iErr) {
+    myErrorStatus=20;
+    return;
+  }
+  //
+  myBuilder=new BOPAlgo_Builder;
+  //
+  aIt.Initialize(aS1);
+  for (; aIt.More(); aIt.Next()) {
+    const TopoDS_Shape& aS1x=aIt.Value();
+    myBuilder->AddArgument(aS1x);
+  }
+  //
+  myBuilder->PerformWithFiller(aPF);
+  iErr=myBuilder->ErrorStatus();
+  if (iErr) {
+    myErrorStatus=21;
+    return;
+  }
+}
+//=======================================================================
+//function : Transform2
+//purpose  : 
+//=======================================================================
+void HYDROData_Transform::Transform2()
+{
+  Standard_Boolean bIsDone;
+  //
+  myErrorStatus=0;
+  //
+  const TopoDS_Shape& aR1=myBuilder->Shape();
+  //
+  myTransform2=new BRepBuilderAPI_Transform(myTrsf2);
+  //
+  myTransform2->Perform(aR1);
+  bIsDone=myTransform2->IsDone();
+  if (!bIsDone) {
+    myErrorStatus=12;
+    return;
+  }
+  //
+  const TopoDS_Shape& aR2=myTransform2->Shape();
+  //
+  myShape=aR2;
+}
+//
+//=======================================================================
+//function : Modified
+//purpose  : 
+//=======================================================================
+const TopTools_ListOfShape& HYDROData_Transform::Modified
+  (const TopoDS_Shape& aS)
+{
+  TopTools_ListOfShape* pHistShapes;
+  TopTools_ListIteratorOfListOfShape aItLS;
+  //
+  pHistShapes=(TopTools_ListOfShape*)&myHistShapes;
+  pHistShapes->Clear();
+  //
+  const TopTools_ListOfShape& aLS1=myTransform1->Modified(aS);
+  const TopoDS_Shape& aS1=aLS1.First();
+  //
+  const TopTools_ListOfShape& aLS1B=myBuilder->Modified(aS1);
+  aItLS.Initialize(aLS1B);
+  for (; aItLS.More(); aItLS.Next()) {
+    const TopoDS_Shape& aS1B=aItLS.Value();
+    const TopTools_ListOfShape& aLS2B=myTransform2->Modified(aS1B);
+    const TopoDS_Shape& aS2B=aLS2B.First();
+    pHistShapes->Append(aS2B);
+  }
+  //
+  return myHistShapes;
+}
+//=======================================================================
+//function : Generated
+//purpose  : 
+//=======================================================================
+const TopTools_ListOfShape& HYDROData_Transform::Generated
+  (const TopoDS_Shape& aS)
+{
+  TopTools_ListOfShape* pHistShapes;
+  TopTools_ListIteratorOfListOfShape aItLS;
+  //
+  pHistShapes=(TopTools_ListOfShape*)&myHistShapes;
+  pHistShapes->Clear();
+  //
+  const TopTools_ListOfShape& aLS1=myTransform1->Modified(aS);
+  const TopoDS_Shape& aS1=aLS1.First();
+  //
+  const TopTools_ListOfShape& aLS1B=myBuilder->Generated(aS1);
+  aItLS.Initialize(aLS1B);
+  for (; aItLS.More(); aItLS.Next()) {
+    const TopoDS_Shape& aS1B=aItLS.Value();
+    const TopTools_ListOfShape& aLS2B=myTransform2->Modified(aS1B);
+    const TopoDS_Shape& aS2B=aLS2B.First();
+    pHistShapes->Append(aS2B);
+  }
+  //
+  return myHistShapes;
+}
+//=======================================================================
+//function : IsDeleted
+//purpose  : 
+//=======================================================================
+Standard_Boolean HYDROData_Transform::IsDeleted(const TopoDS_Shape& aS)
+{
+  Standard_Boolean bIsDeleted;
+  //
+  const TopTools_ListOfShape& aLS1=myTransform1->Modified(aS);
+  const TopoDS_Shape& aS1=aLS1.First();
+  //
+  bIsDeleted=myBuilder->IsDeleted(aS1);
+  //
+  return bIsDeleted;
+}
+//=======================================================================
+//function : HasDeleted
+//purpose  : 
+//=======================================================================
+Standard_Boolean HYDROData_Transform::HasDeleted()
+{
+  return myBuilder->HasDeleted();
+}
+//=======================================================================
+//function : HasGenerated
+//purpose  : 
+//=======================================================================
+Standard_Boolean HYDROData_Transform::HasGenerated()
+{
+  return myBuilder->HasGenerated();
+}
+//=======================================================================
+//function : HasModified
+//purpose  : 
+//=======================================================================
+Standard_Boolean HYDROData_Transform::HasModified()
+{
+  return myBuilder->HasModified();
+}
diff --git a/src/HYDROData/HYDROData_Transform.h b/src/HYDROData/HYDROData_Transform.h
new file mode 100644 (file)
index 0000000..d0d6a4e
--- /dev/null
@@ -0,0 +1,116 @@
+#ifndef HYDROData_Transform_HeaderFile
+#define HYDROData_Transform_HeaderFile
+
+#include <Standard.hxx>
+#include <Standard_Macro.hxx>
+#include <Standard_Boolean.hxx>
+#include <Standard_Integer.hxx>
+
+#include <gp_Trsf.hxx>
+
+#include <TopoDS_Shape.hxx>
+#include <TopTools_ListOfShape.hxx>
+
+#include <BRepBuilderAPI_Transform.hxx>
+
+#include <BOPCol_ListOfShape.hxx>
+#include <BOPCol_MapOfShape.hxx>
+
+#include <BOPAlgo_BuilderShape.hxx>
+#include <BOPAlgo_Builder.hxx>
+
+//=======================================================================
+//class    : HYDROData_Transform
+//purpose  :
+//=======================================================================
+class HYDROData_Transform : public BOPAlgo_BuilderShape
+{
+ public:
+  Standard_EXPORT
+    HYDROData_Transform();
+  //
+  Standard_EXPORT
+    virtual ~HYDROData_Transform();
+  //
+  Standard_EXPORT
+    void SetArgument(const TopoDS_Shape& theShape);
+  //
+  Standard_EXPORT
+    const TopoDS_Shape& Argument()const;
+  //
+  Standard_EXPORT
+    void SetToTransform(const Standard_Boolean theFlag);
+  //
+  Standard_EXPORT
+    Standard_Boolean ToTransform()const;
+  //
+  Standard_EXPORT
+    void SetTreshold(const Standard_Real theTreshold);
+  //
+  Standard_EXPORT
+    Standard_Real Treshold()const;
+  //
+  Standard_EXPORT
+    void Detect();
+  //
+  Standard_EXPORT
+    virtual void Clear();
+  //
+  Standard_EXPORT
+    virtual void Perform();
+  //
+  // History
+  Standard_EXPORT
+    virtual const TopTools_ListOfShape& Generated
+      (const TopoDS_Shape& theShape);
+  //
+  Standard_EXPORT
+    virtual const TopTools_ListOfShape& Modified
+      (const TopoDS_Shape& theShape);
+  //
+  Standard_EXPORT
+    virtual Standard_Boolean IsDeleted
+      (const TopoDS_Shape& theShape);
+  //
+  Standard_EXPORT
+    virtual Standard_Boolean HasDeleted();
+  //
+  Standard_EXPORT
+    virtual Standard_Boolean HasGenerated();
+  //
+  Standard_EXPORT
+    virtual Standard_Boolean HasModified();
+  //
+  // Inner content
+ protected:
+  Standard_EXPORT
+    virtual void CheckData();
+  //
+  Standard_EXPORT
+    void Prepare();
+  //
+  Standard_EXPORT
+    void Transform1();
+  //
+  Standard_EXPORT
+    void Transform2();
+  //
+  Standard_EXPORT
+    void BuildFuse();
+  //
+  Standard_EXPORT
+    void ComputeTrsf();
+  //
+ protected:
+  Standard_Boolean myToTransform;
+  Standard_Real myTreshold;
+  TopoDS_Shape myArgument;
+  //
+  gp_Trsf myTrsf1;
+  gp_Trsf myTrsf2;
+  BRepBuilderAPI_Transform* myTransform1;
+  BRepBuilderAPI_Transform* myTransform2;
+  BOPAlgo_Builder* myBuilder;
+};
+
+#endif