src/HYDROData/HYDROData_Region.cxx

   1
   2 #include "HYDROData_Region.h"
   3
   4 #include "HYDROData_CalculationCase.h"
   5 #include "HYDROData_Document.h"
   6 #include "HYDROData_Iterator.h"
   7 #include "HYDROData_Zone.h"
   8
   9 #include <TopoDS.hxx>
  10 #include <TopoDS_Shape.hxx>
  11 #include <TopoDS_Shell.hxx>
  12 #include <TopoDS_Face.hxx>
  13 #include <TopTools_ListOfShape.hxx>
  14 #include <TopTools_ListIteratorOfListOfShape.hxx>
  15 #include <TopTools_IndexedMapOfShape.hxx>
  16 #include <TopExp.hxx>
  17 #include <BRep_Builder.hxx>
  18 #include <BRepAlgoAPI_Fuse.hxx>
  19 #include <ShapeUpgrade_UnifySameDomain.hxx>
  20
  21 #include <QStringList>
  22
  23 #define PYTHON_REGION_ID "KIND_REGION"
  24
  25 IMPLEMENT_STANDARD_HANDLE(HYDROData_Region, HYDROData_Entity)
  26 IMPLEMENT_STANDARD_RTTIEXT(HYDROData_Region, HYDROData_Entity)
  27
  28
  29 HYDROData_Region::HYDROData_Region()
  30  : HYDROData_Entity()
  31 {
  32 }
  33
  34 HYDROData_Region::~HYDROData_Region()
  35 {
  36 }
  37
  38 QStringList HYDROData_Region::DumpToPython( MapOfTreatedObjects& theTreatedObjects ) const
  39 {
  40   QStringList aResList;
  41
  42   Handle(HYDROData_Document) aDocument = HYDROData_Document::Document( myLab );
  43   if ( aDocument.IsNull() )
  44     return aResList;
  45
  46   QString aDocName = aDocument->GetDocPyName();
  47   QString aRegionName = GetName();
  48
  49   aResList << QString( "%1 = %2.CreateObject( %3 );" )
  50               .arg( aRegionName ).arg( aDocName ).arg( PYTHON_REGION_ID );
  51   aResList << QString( "%1.SetName( \"%2\" );" )
  52               .arg( aRegionName ).arg( aRegionName );
  53   aResList << QString( "" );
  54
  55   HYDROData_SequenceOfObjects aZones = GetZones();
  56   HYDROData_SequenceOfObjects::Iterator anIter( aZones );
  57   for ( ; anIter.More(); anIter.Next() )
  58   {
  59     Handle(HYDROData_Zone) aRefZone =
  60       Handle(HYDROData_Zone)::DownCast( anIter.Value() );
  61     if ( !aRefZone.IsNull() )
  62       setPythonReferenceObject( theTreatedObjects, aResList, aRefZone, "AddZone" );
  63   }
  64   aResList << QString( "" );
  65
  66   return aResList;
  67 }
  68
  69 void HYDROData_Region::Remove()
  70 {
  71   Handle(HYDROData_CalculationCase) aFatherCalc =
  72     Handle(HYDROData_CalculationCase)::DownCast( GetFatherObject() );
  73
  74   HYDROData_Entity::Remove();
  75
  76   if ( !aFatherCalc.IsNull() )
  77     aFatherCalc->UpdateRegionsOrder();
  78 }
  79
  80 bool HYDROData_Region::CanRemove()
  81 {
  82   return false;
  83 }
  84
  85 HYDROData_SequenceOfObjects HYDROData_Region::GetAllReferenceObjects() const
  86 {
  87   HYDROData_SequenceOfObjects aResSeq = HYDROData_Entity::GetAllReferenceObjects();
  88   aResSeq.Append( GetZones() );
  89   return aResSeq;
  90 }
  91
  92 bool HYDROData_Region::AddZone( const Handle(HYDROData_Zone)& theZone )
  93 {
  94   if ( theZone.IsNull() )
  95     return false;
  96
  97   if ( HasReference( theZone, DataTag_Zone ) )
  98     return false; // Object is already in reference list
  99
 100   // Move the zone from other region
 101   Handle(HYDROData_Region) aFatherRegion =
 102     Handle(HYDROData_Region)::DownCast( theZone->GetFatherObject() );
 103   if ( !aFatherRegion.IsNull() && aFatherRegion->Label() != myLab )
 104   {
 105     Handle(HYDROData_Zone) aNewZone = addNewZone();
 106     theZone->CopyTo( aNewZone );
 107
 108     aFatherRegion->RemoveZone( theZone );
 109
 110     theZone->SetLabel( aNewZone->Label() );
 111   }
 112   else
 113   {
 114     AddReferenceObject( theZone, DataTag_Zone );
 115   }
 116
 117   return true;
 118 }
 119
 120 HYDROData_SequenceOfObjects HYDROData_Region::GetZones() const
 121 {
 122   return GetReferenceObjects( DataTag_Zone );
 123 }
 124
 125 void HYDROData_Region::RemoveZone( const Handle(HYDROData_Zone)& theZone )
 126 {
 127   if ( theZone.IsNull() )
 128     return;
 129
 130   RemoveReferenceObject( theZone->Label(), DataTag_Zone );
 131
 132   // Remove zone from data model
 133   Handle(HYDROData_Region) aFatherRegion =
 134     Handle(HYDROData_Region)::DownCast( theZone->GetFatherObject() );
 135   if ( !aFatherRegion.IsNull() && aFatherRegion->Label() == myLab )
 136     theZone->Remove();
 137
 138   // If the last zone has been removed from region we remove this region
 139   HYDROData_SequenceOfObjects aRefZones = GetZones();
 140   if ( aRefZones.IsEmpty() )
 141     Remove();
 142 }
 143
 144 void HYDROData_Region::RemoveZones()
 145 {
 146   ClearReferenceObjects( DataTag_Zone );
 147   myLab.FindChild( DataTag_ChildZone ).ForgetAllAttributes( true );
 148 }
 149
 150 Handle(HYDROData_Zone) HYDROData_Region::addNewZone()
 151 {
 152   TDF_Label aNewLab = myLab.FindChild( DataTag_ChildZone ).NewChild();
 153
 154   Handle(HYDROData_Zone) aNewZone =
 155     Handle(HYDROData_Zone)::DownCast( HYDROData_Iterator::CreateObject( aNewLab, KIND_ZONE ) );
 156   AddZone( aNewZone );
 157
 158   return aNewZone;
 159 }
 160
 161 TopoDS_Shape HYDROData_Region::GetShape() const
 162 {
 163   TopoDS_Shape aResShape;
 164
 165   // Unite the region zones (each zone is a face) into one face (united face)
 166   // If the zones can't be united into the single face - unite them into shell
 167
 168   // Collect the list of region faces
 169   TopTools_ListOfShape aRegionFacesList;
 170
 171   HYDROData_SequenceOfObjects aZones = GetZones();
 172   HYDROData_SequenceOfObjects::Iterator aZoneIter( aZones );
 173   for ( ; aZoneIter.More(); aZoneIter.Next() ) {
 174     Handle(HYDROData_Zone) aZone =
 175       Handle(HYDROData_Zone)::DownCast( aZoneIter.Value() );
 176
 177     if ( aZone.IsNull() ) {
 178       continue;
 179     }
 180
 181     TopoDS_Face aFace = TopoDS::Face( aZone->GetShape() );
 182     if ( !aFace.IsNull() ) {
 183       aRegionFacesList.Append( aFace );
 184     }
 185   } // zones iterator
 186
 187   // Check number of faces
 188   int aNbFaces = aRegionFacesList.Extent();
 189   if ( aNbFaces > 0 ) {
 190     // The unite region face
 191     TopoDS_Face aRegionFace;
 192
 193     if ( aNbFaces == 1 ) {
 194       aRegionFace = TopoDS::Face( aRegionFacesList.First() );
 195     } else {
 196       // Fuse faces into one
 197       TopoDS_Shape aFuseShape;
 198       TopTools_ListIteratorOfListOfShape aFaceIter( aRegionFacesList );
 199       for ( ; aFaceIter.More(); aFaceIter.Next() ) {
 200         if ( aFuseShape.IsNull() ) {
 201           aFuseShape = aFaceIter.Value();
 202         } else {
 203           BRepAlgoAPI_Fuse aFuse(aFuseShape, aFaceIter.Value());
 204           if ( !aFuse.IsDone() ) {
 205             aFuseShape.Nullify();
 206             break;
 207           }
 208           aFuseShape = aFuse.Shape();
 209         }
 210       } // faces iterator
 211
 212       // Check the result of fuse operation
 213       if ( !aFuseShape.IsNull() ) {
 214         ShapeUpgrade_UnifySameDomain anUnifier( aFuseShape );
 215         anUnifier.Build();
 216         TopoDS_Shape anUnitedShape = anUnifier.Shape();
 217
 218         TopTools_IndexedMapOfShape aMapOfFaces;
 219         TopExp::MapShapes( anUnitedShape, TopAbs_FACE, aMapOfFaces );
 220         if ( aMapOfFaces.Extent() == 1 ) {
 221           aRegionFace = TopoDS::Face( aMapOfFaces(1) );
 222         }
 223       }
 224     }
 225
 226     if ( !aRegionFace.IsNull() ) { // result shape is a face
 227       aResShape = aRegionFace;
 228     } else { // result shape is a shell
 229       TopoDS_Shell aShell;
 230       BRep_Builder aBuilder;
 231       aBuilder.MakeShell( aShell );
 232
 233       TopTools_ListIteratorOfListOfShape aFaceIter( aRegionFacesList );
 234       for ( ; aFaceIter.More(); aFaceIter.Next() ) {
 235         aBuilder.Add( aShell, aFaceIter.Value() );
 236       }
 237
 238       aResShape = aShell;
 239     }
 240   }
 241
 242   return aResShape;
 243 }