src/BOP/BOP_WireEdgeSet.cxx

   1 // Created on: 1993-06-17
   2 // Created by: Jean Yves LEBEY
   3 // Copyright (c) 1993-1999 Matra Datavision
   4 // Copyright (c) 1999-2012 OPEN CASCADE SAS
   5 //
   6 // The content of this file is subject to the Open CASCADE Technology Public
   7 // License Version 6.5 (the "License"). You may not use the content of this file
   8 // except in compliance with the License. Please obtain a copy of the License
   9 // at http://www.opencascade.org and read it completely before using this file.
  10 //
  11 // The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
  12 // main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
  13 //
  14 // The Original Code and all software distributed under the License is
  15 // distributed on an "AS IS" basis, without warranty of any kind, and the
  16 // Initial Developer hereby disclaims all such warranties, including without
  17 // limitation, any warranties of merchantability, fitness for a particular
  18 // purpose or non-infringement. Please see the License for the specific terms
  19 // and conditions governing the rights and limitations under the License.
  20
  21 // modified by NIZNHY-PKV
  22
  23 #include <BOP_WireEdgeSet.ixx>
  24
  25 #include <Standard_ProgramError.hxx>
  26
  27 #include <Precision.hxx>
  28
  29 #include <gp_Vec.hxx>
  30 #include <gp_Vec2d.hxx>
  31 #include <gp_Pnt2d.hxx>
  32 #include <gp_Dir2d.hxx>
  33
  34 #include <Geom2d_Curve.hxx>
  35 #include <Geom2d_Line.hxx>
  36 #include <Geom_Curve.hxx>
  37 #include <Geom_Surface.hxx>
  38 #include <GeomAbs_SurfaceType.hxx>
  39 #include <GeomAbs_CurveType.hxx>
  40 #include <GeomAPI_ProjectPointOnSurf.hxx>
  41
  42 #include <BRep_Tool.hxx>
  43 #include <BRepAdaptor_Surface.hxx>
  44 #include <BRepAdaptor_Curve.hxx>
  45
  46 #include <TopoDS.hxx>
  47 #include <TopAbs.hxx>
  48 #include <TopExp_Explorer.hxx>
  49
  50 #include <BOPTools_Tools2D.hxx>
  51
  52 // modified by NIZHNY-MKK  Fri May 30 10:18:01 2003.BEGIN
  53 #include <TopTools_IndexedMapOfShape.hxx>
  54 #include <TopExp.hxx>
  55 #include <BOPTools_Tools3D.hxx>
  56
  57 static void GetOrientationVOnE(const TopoDS_Shape& V,
  58                                const TopoDS_Shape& E1,
  59                                const TopoDS_Shape& E2,
  60                                TopAbs_Orientation& o1,
  61                                TopAbs_Orientation& o2);
  62 // modified by NIZHNY-MKK  Fri May 30 10:18:03 2003.END
  63
  64 //=======================================================================
  65 //function : BOP_WireEdgeSet::BOP_WireEdgeSet
  66 //purpose  :
  67 //=======================================================================
  68   BOP_WireEdgeSet::BOP_WireEdgeSet()
  69 :
  70   BOP_ShapeSet(TopAbs_VERTEX)
  71 {
  72 }
  73 //=======================================================================
  74 //function : BOP_WireEdgeSet::BOP_WireEdgeSet
  75 //purpose  :
  76 //=======================================================================
  77   BOP_WireEdgeSet::BOP_WireEdgeSet(const TopoDS_Face& F)
  78 :
  79   BOP_ShapeSet(TopAbs_VERTEX),
  80   myFace(F)
  81 {
  82 }
  83
  84 //=======================================================================
  85 //function : Initialize
  86 //purpose  :
  87 //=======================================================================
  88   void BOP_WireEdgeSet::Initialize(const TopoDS_Face& F)
  89 {
  90   myFace=F;
  91 }
  92
  93 //=======================================================================
  94 //function : AddShape
  95 //purpose  :
  96 //=======================================================================
  97   void BOP_WireEdgeSet::AddShape(const TopoDS_Shape& S)
  98 {
  99   ProcessAddShape(S);
 100 }
 101
 102 //=======================================================================
 103 //function : AddStartElement
 104 //purpose  :
 105 //=======================================================================
 106   void BOP_WireEdgeSet::AddStartElement(const TopoDS_Shape& S)
 107 {
 108   ProcessAddStartElement(S);
 109 }
 110
 111 //=======================================================================
 112 //function : AddElement
 113 //purpose  :
 114 //=======================================================================
 115   void BOP_WireEdgeSet::AddElement(const TopoDS_Shape& S)
 116 {
 117   BOP_ShapeSet::AddElement(S);
 118 }
 119
 120 //=======================================================================
 121 //function : InitNeighbours
 122 //purpose  :
 123 //=======================================================================
 124   void  BOP_WireEdgeSet::InitNeighbours(const TopoDS_Shape& E)
 125 {
 126   mySubShapeExplorer.Init(E,mySubShapeType);
 127   myCurrentShape = E;
 128
 129   FindNeighbours();
 130 }
 131
 132
 133 //=======================================================================
 134 //function : FindNeighbours
 135 //purpose  :
 136 //=======================================================================
 137   void BOP_WireEdgeSet::FindNeighbours()
 138 {
 139   while (mySubShapeExplorer.More()) {
 140
 141     // l = list of edges neighbour of edge myCurrentShape trough
 142     // the vertex mySubShapeExplorer.Current(), which is a vertex of the
 143     // edge myCurrentShape.
 144     const TopoDS_Shape& V = mySubShapeExplorer.Current();
 145
 146     const TopTools_ListOfShape & l =
 147       MakeNeighboursList(myCurrentShape, V);
 148
 149     // myIncidentShapesIter iterates on the neighbour edges of the edge
 150     // given as InitNeighbours() argument (this edge has been stored
 151     // in the field myCurrentShape).
 152
 153     myIncidentShapesIter.Initialize(l);
 154     if (myIncidentShapesIter.More()) {
 155       break;
 156     }
 157     else{
 158       mySubShapeExplorer.Next();
 159     }
 160   }
 161 }
 162
 163 //=======================================================================
 164 //function : MakeNeighboursList
 165 //purpose  : find edges connected to Earg by Varg
 166 //=======================================================================
 167   const TopTools_ListOfShape & BOP_WireEdgeSet::MakeNeighboursList(const TopoDS_Shape& Earg,
 168                                                                    const TopoDS_Shape& Varg)
 169 {
 170   const TopoDS_Edge&   E = TopoDS::Edge(Earg);
 171   const TopoDS_Vertex& V = TopoDS::Vertex(Varg);
 172
 173   const TopTools_ListOfShape& l = mySubShapeMap.FindFromKey(Varg);
 174
 175   Standard_Integer nclosing = NbClosingShapes(l);
 176   if (nclosing) {
 177     // build myCurrentShapeNeighbours =
 178     // edge list made of connected shapes to Earg through Varg
 179     myCurrentShapeNeighbours.Clear();
 180
 181     Standard_Integer iapp = 0;
 182     for (TopTools_ListIteratorOfListOfShape it(l); it.More(); it.Next()) {
 183       iapp++;
 184       const TopoDS_Shape& curn = it.Value(); // current neighbour
 185       Standard_Boolean k = VertexConnectsEdgesClosing(Varg, Earg, curn);
 186       if (k) {
 187         myCurrentShapeNeighbours.Append(curn);
 188       }
 189     }
 190
 191     Standard_Integer newn = NbClosingShapes(myCurrentShapeNeighbours);
 192     if (newn >= 2 ) {
 193
 194       const TopoDS_Face& F = myFace;
 195       // several sewing edges connected to E by V such as :
 196       // orientation of V in E # orientation of V in its edges.
 197       // one leaves among the connected sewing edges,
 198       // only edge A that checks tg(E) ^ tg(A) > 0
 199
 200       TopAbs_Orientation Eori;
 201       gp_Vec2d d1E; gp_Pnt2d pE;
 202       Standard_Real fiE, laE, tolpc, parE;
 203
 204       parE = BRep_Tool::Parameter(V,E);
 205
 206       Handle(Geom2d_Curve) PCE;
 207       BOPTools_Tools2D::CurveOnSurface(E, F, PCE, fiE, laE, tolpc, Standard_False);
 208
 209       if (!PCE.IsNull()){
 210         PCE->D1(parE, pE, d1E);
 211       }
 212       else  {
 213         LocalD1 (F, E, V, pE, d1E);
 214       }
 215
 216       Eori = E.Orientation();
 217       if (Eori == TopAbs_REVERSED) {
 218         d1E.Reverse();
 219       }
 220
 221       TopTools_ListIteratorOfListOfShape lclo(myCurrentShapeNeighbours);
 222       Standard_Integer rang = 0;
 223       while (lclo.More()) {
 224         rang++;
 225
 226         if ( ! IsClosed(lclo.Value()) ) {
 227           lclo.Next();
 228           continue;
 229         }
 230
 231         const TopoDS_Edge& EE = TopoDS::Edge(lclo.Value());
 232         Standard_Real parEE,  fiEE, laEE, tolpcc;
 233         gp_Vec2d d1EE; gp_Pnt2d pEE;
 234
 235         parEE = BRep_Tool::Parameter(V,EE);
 236
 237         Handle(Geom2d_Curve) PCEE;
 238         BOPTools_Tools2D::CurveOnSurface (EE, F, PCEE, fiEE, laEE, tolpcc, Standard_False);//
 239         //
 240         if (!PCEE.IsNull()) {
 241           PCEE->D1(parEE,pEE,d1EE);
 242
 243         }
 244         else {
 245           LocalD1(F, EE, V, pEE, d1EE);
 246         }
 247
 248         TopAbs_Orientation EEori = EE.Orientation();
 249         if (EEori == TopAbs_REVERSED) d1EE.Reverse();
 250
 251         Standard_Real cross = d1E.Crossed(d1EE);
 252         TopAbs_Orientation oVE,oVEE;
 253
 254         // modified by NIZHNY-MKK  Fri May 30 10:17:10 2003.BEGIN
 255         GetOrientationVOnE(V,E,EE,oVE,oVEE);
 256 //      VertexConnectsEdges(V,E,EE,oVE,oVEE);
 257         // modified by NIZHNY-MKK  Fri May 30 10:17:13 2003.END
 258
 259         Standard_Boolean t2 = ( (cross > 0) && oVE == TopAbs_REVERSED ) ||
 260                               ( (cross < 0) && oVE == TopAbs_FORWARD );
 261
 262         if ( t2 ) { //-- t1
 263           // this is proper IsClosed, only it is preserved among IsClosed
 264           lclo.Next();
 265         }
 266         else {
 267           // edge IsClosed is suspended
 268           myCurrentShapeNeighbours.Remove(lclo);
 269         }
 270       }// end of while (lclo.More())
 271     }// end of if (newn >= 2 )
 272     return myCurrentShapeNeighbours;
 273   }// end of if (nclosing)
 274
 275   else {
 276     return l;
 277   }
 278 } // MakeNeighoursList
 279
 280 //=======================================================================
 281 //function : VertexConnectsEdges
 282 //purpose  :
 283 //=======================================================================
 284   Standard_Boolean BOP_WireEdgeSet::VertexConnectsEdges(const TopoDS_Shape& V,
 285                                                         const TopoDS_Shape& E1,
 286                                                         const TopoDS_Shape& E2,
 287                                                         TopAbs_Orientation& o1,
 288                                                         TopAbs_Orientation& o2) const
 289 {
 290   TopExp_Explorer ex1, ex2;
 291   ex1.Init(E1,TopAbs_VERTEX);
 292   for(; ex1.More(); ex1.Next()) {
 293
 294     if (V.IsSame(ex1.Current())) {
 295       ex2.Init(E2,TopAbs_VERTEX);
 296       for(; ex2.More(); ex2.Next()) {
 297
 298         if (V.IsSame(ex2.Current())) {
 299           o1 = ex1.Current().Orientation();
 300           o2 = ex2.Current().Orientation();
 301           if ( o1 != o2 ) {
 302             return Standard_True;
 303           }
 304         }
 305       }
 306     }
 307   }
 308   return Standard_False;
 309 }
 310
 311
 312 //=======================================================================
 313 //function : VertexConnectEdgesClosing
 314 //purpose  :
 315 //=======================================================================
 316   Standard_Boolean BOP_WireEdgeSet::VertexConnectsEdgesClosing(const TopoDS_Shape& V,
 317                                                                const TopoDS_Shape& E1,
 318                                                                const TopoDS_Shape& E2) const
 319 {
 320 //  Let S the set of incident edges on vertex V.
 321 //  S contains at least one closed edge on the periodic face to build.
 322 // (the face implied in closing test of edge is myFace)
 323 //  E1,E2 are S shapes (sharing V).
 324 //
 325 //  if E1 and E2 are not closed : edges are NOT connected
 326 //  if E1 or E2 is/are closed :
 327 //    if V changes of relative orientation between E1,E2 : edges are connected
 328 //    else : edges are NOT connected
 329 //
 330 //  example with E1 NOT closed, E2 closed :
 331 //
 332 //  E1 FORWARD, V REVERSED on E1
 333 //  E2 FORWARD, V FORWARD on E2       --> edges are connected
 334 //
 335 //  E1 FORWARD, V REVERSED on E1
 336 //  E2 REVERSED, V FORWARD on E2      --> edges are NOT connected
 337 //-----------------------------------------------------------------------
 338
 339   Standard_Boolean c1 = IsClosed(E1);
 340   Standard_Boolean c2 = IsClosed(E2);
 341
 342   Standard_Boolean testconnect = c1 || c2;
 343   Standard_Boolean resu = Standard_False;
 344   TopAbs_Orientation o1,o2;
 345
 346   // SSCONNEX = False ==> E2 is  selected to create as many
 347   // (default)           faces as there are components connected in UV.
 348   // SSCONNEX = True ==> the entire edge E2 sharing V with E1 is taken
 349   //                     so that orientation(V/E1) # orientation(V/E2)
 350   //                     ==> face of part and other sewings
 351   if ((c1 && c2)) {
 352     Standard_Boolean u1 = c1 ? IsUClosed(E1) : Standard_False;
 353     Standard_Boolean v1 = c1 ? IsVClosed(E1) : Standard_False;
 354     Standard_Boolean u2 = c2 ? IsUClosed(E2) : Standard_False;
 355     Standard_Boolean v2 = c2 ? IsVClosed(E2) : Standard_False;
 356     Standard_Boolean uvdiff = (u1 && v2) || (u2 && v1);
 357     testconnect = uvdiff;
 358   }
 359
 360   if (testconnect) {
 361     resu = VertexConnectsEdges(V,E1,E2,o1,o2);
 362   }
 363   else {
 364     // cto 012 O2 sewing edge of cylindric face
 365     // chain of components split ON and OUT of the same orientation
 366     TopAbs_Orientation oe1 = E1.Orientation();
 367     TopAbs_Orientation oe2 = E2.Orientation();
 368     Standard_Boolean iseq = E1.IsEqual(E2);
 369     if ( (c1 && c2) && (oe1 == oe2) && (!iseq) ) {
 370       resu = VertexConnectsEdges(V,E1,E2,o1,o2);
 371     }
 372   }
 373   return resu;
 374 }
 375
 376 //=======================================================================
 377 //function : NbClosingShapes
 378 //purpose  :
 379 //=======================================================================
 380   Standard_Integer BOP_WireEdgeSet::NbClosingShapes(const TopTools_ListOfShape & L) const
 381 {
 382   Standard_Integer n = 0;
 383   TopTools_ListIteratorOfListOfShape it(L);
 384   for ( ; it.More(); it.Next()) {
 385     const TopoDS_Shape& S = it.Value();
 386     if ( IsClosed(S) ) {
 387       n++;
 388     }
 389   }
 390   return n;
 391 }
 392
 393 //=======================================================================
 394 //function : LocalD1
 395 //purpose  :
 396 //=======================================================================
 397   void BOP_WireEdgeSet::LocalD1(const TopoDS_Face&  F,
 398                                 const TopoDS_Edge&  E,
 399                                 const TopoDS_Vertex& V,
 400                                 gp_Pnt2d& ,//pE,
 401                                 gp_Vec2d& d1E) const
 402 {
 403   TopLoc_Location Loc;
 404   Standard_Real parE, fiE, laE, u, v;
 405
 406   parE = BRep_Tool::Parameter(V,E);
 407
 408   Handle(Geom_Curve) CE = BRep_Tool::Curve(E, Loc, fiE, laE);
 409   CE = Handle(Geom_Curve)::DownCast(CE->Transformed(Loc.Transformation()));
 410
 411   gp_Pnt p3dE; gp_Vec d3dE;
 412   CE->D1(parE, p3dE, d3dE);
 413
 414   Handle(Geom_Surface) S = BRep_Tool::Surface(F);
 415   GeomAPI_ProjectPointOnSurf proj(p3dE,S);
 416
 417   proj.LowerDistanceParameters(u,v);
 418   gp_Pnt bid;
 419   gp_Vec d1u,d1v;
 420   S->D1(u,v,bid,d1u,d1v);
 421   u = d3dE.Dot(d1u);
 422   v = d3dE.Dot(d1v);
 423   d1E.SetCoord(u,v);
 424 }
 425
 426 //=======================================================================
 427 //function : IsClosed
 428 //purpose  :
 429 //=======================================================================
 430   Standard_Boolean BOP_WireEdgeSet::IsClosed(const TopoDS_Shape& E) const
 431 {
 432   const TopoDS_Edge& EE = TopoDS::Edge(E);
 433   Standard_Boolean closed = BRep_Tool::IsClosed(EE,myFace);
 434
 435   // modified by NIZHNY-MKK  Wed May 12 15:44:10 2004.BEGIN
 436   if(closed) {
 437     TopoDS_Edge aE1, aE2;
 438     BOPTools_Tools3D::GetSeam (myFace, aE1, aE2);
 439
 440     if(aE1.IsNull() || aE2.IsNull()) {
 441       closed = Standard_False;
 442     }
 443   }
 444   // modified by NIZHNY-MKK  Wed May 12 15:44:14 2004.END
 445
 446   return closed;
 447 }
 448
 449 //=======================================================================
 450 //function : IsUVISO
 451 //purpose  :
 452 //=======================================================================
 453   void BOP_WireEdgeSet::IsUVISO(const TopoDS_Edge& E,
 454                                 const TopoDS_Face& F,
 455                                 Standard_Boolean& uiso,
 456                                 Standard_Boolean& viso)
 457 {
 458   uiso = viso = Standard_False;
 459   Standard_Real fE, lE,tolpc;
 460
 461   Handle(Geom2d_Curve) PC;
 462   BOPTools_Tools2D::CurveOnSurface (E, F, PC, fE, lE, tolpc, Standard_False);//
 463
 464   if (PC.IsNull()) {
 465     Standard_ProgramError::Raise("BOP_WireEdgeSet::IsUVISO");
 466   }
 467
 468   Handle(Standard_Type) TheType = PC->DynamicType();
 469   if (TheType == STANDARD_TYPE(Geom2d_Line)) {
 470     const Handle(Geom2d_Line)& HL = *((Handle(Geom2d_Line)*)&PC);
 471     const gp_Dir2d&  D = HL->Direction();
 472     Standard_Real    tol = Precision::Angular();
 473     if      (D.IsParallel(gp_Dir2d(0.,1.),tol)) {
 474       uiso = Standard_True;
 475     }
 476     else if (D.IsParallel(gp_Dir2d(1.,0.),tol)) {
 477       viso = Standard_True;
 478     }
 479   }
 480 }
 481
 482 //=======================================================================
 483 //function : IsUClosed
 484 //purpose  :
 485 //=======================================================================
 486   Standard_Boolean BOP_WireEdgeSet::IsUClosed(const TopoDS_Shape& E) const
 487 {
 488   const TopoDS_Edge& EE = TopoDS::Edge(E);
 489
 490   Standard_Boolean bid, closed;
 491   IsUVISO(EE, myFace, closed, bid);
 492   return closed;
 493 }
 494
 495
 496 //=======================================================================
 497 //function : IsVClosed
 498 //purpose  :
 499 //=======================================================================
 500   Standard_Boolean BOP_WireEdgeSet::IsVClosed(const TopoDS_Shape& E) const
 501 {
 502   const TopoDS_Edge& EE = TopoDS::Edge(E);
 503
 504   Standard_Boolean bid, closed;
 505   IsUVISO (EE, myFace, bid, closed);
 506   return closed;
 507 }
 508
 509 //=======================================================================
 510 //function : Face
 511 //purpose  :
 512 //=======================================================================
 513   const TopoDS_Face& BOP_WireEdgeSet::Face() const
 514 {
 515   return myFace;
 516 }
 517 //=======================================================================
 518 //function : BOP_Tools::KeptTwice
 519 //purpose  :
 520 //=======================================================================
 521   Standard_Boolean BOP_WireEdgeSet::KeptTwice (const TopoDS_Edge& aSpEF2) const
 522
 523
 524 {
 525   Standard_Integer iCnt=0;
 526
 527   const TopTools_ListOfShape& aWESList=StartElements();
 528
 529   TopTools_ListIteratorOfListOfShape anIt(aWESList);
 530   for (; anIt.More(); anIt.Next()) {
 531     const TopoDS_Shape& aSS=anIt.Value();
 532     if (aSpEF2.IsSame(aSS)) {
 533       iCnt++;
 534       if (iCnt > 1) {
 535         return Standard_True;
 536       }
 537     }
 538   }
 539  return Standard_False;
 540 }
 541 //=======================================================================
 542 //function : BOP_Tools::RemoveEdgeFromWES
 543 //purpose  :
 544 //=======================================================================
 545   void BOP_WireEdgeSet::RemoveEdgeFromWES(const TopoDS_Edge& aSpEF2)
 546
 547 {
 548   TopTools_ListOfShape aTmpList;
 549   aTmpList=StartElements();
 550
 551   ClearContents();
 552
 553   TopTools_ListIteratorOfListOfShape anIt(aTmpList);
 554   for (; anIt.More(); anIt.Next()) {
 555     const TopoDS_Shape& aSS=anIt.Value();
 556     if (!aSpEF2.IsEqual(aSS)) {
 557       AddStartElement(aSS);
 558     }
 559   }
 560 }
 561
 562
 563 // modified by NIZHNY-MKK  Fri May 30 10:18:38 2003.BEGIN
 564 void GetOrientationVOnE(const TopoDS_Shape& V,
 565                         const TopoDS_Shape& E1,
 566                         const TopoDS_Shape& E2,
 567                         TopAbs_Orientation& o1,
 568                         TopAbs_Orientation& o2) {
 569   TopTools_IndexedMapOfShape aMap;
 570
 571   TopExp::MapShapes(E1, TopAbs_VERTEX, aMap);
 572   TopExp::MapShapes(E2, TopAbs_VERTEX, aMap);
 573
 574   if(aMap.Extent() > 1) {
 575     TopExp_Explorer ex1, ex2;
 576     ex1.Init(E1,TopAbs_VERTEX);
 577
 578     for(; ex1.More(); ex1.Next()) {
 579
 580       if (V.IsSame(ex1.Current())) {
 581         ex2.Init(E2,TopAbs_VERTEX);
 582
 583         for(; ex2.More(); ex2.Next()) {
 584
 585           if (V.IsSame(ex2.Current())) {
 586             o1 = ex1.Current().Orientation();
 587             o2 = ex2.Current().Orientation();
 588
 589             if ( o1 != o2 ) {
 590               return;
 591             }
 592           }
 593         }
 594       }
 595     }
 596   }
 597   else {
 598     TopoDS_Shape VV = V;
 599
 600     if(E1.Orientation() == TopAbs_REVERSED)
 601       VV.Reverse();
 602     o1 = VV.Orientation();
 603     VV.Reverse();
 604     o2 = VV.Orientation();
 605   }
 606   return;
 607 }
 608 // modified by NIZHNY-MKK  Fri May 30 10:18:40 2003.END