Salome HOME
Lot 2: change bathy associated to natural object propagated to all cases without...
[modules/hydro.git] / src / HYDROData / HYDROData_Channel.cxx
index 09ba1c648d875557d293ebf4c927ec17431314ae..0f31fc1f1ac737f96bbd000251424402e7e918d7 100644 (file)
 #include "HYDROData_ShapesTool.h"
 #include "HYDROData_Stream.h"
 #include "HYDROData_Tool.h"
+#include "HYDROData_ChannelAltitude.h"
 
 #include <BRepBuilderAPI_MakeWire.hxx>
+#include <BRepLib_MakePolygon.hxx>
 
 #include <BRepOffsetAPI_MakePipeShell.hxx>
 #include <BRepOffsetAPI_MakePipe.hxx>
 #include <TColStd_Array1OfReal.hxx>
 
 #include <TopTools_HArray1OfShape.hxx>
+#include <TopTools_SequenceOfShape.hxx>
 
 #include <TopoDS.hxx>
 #include <TopoDS_Wire.hxx>
 #include <TopoDS_Vertex.hxx>
 
+#include <Quantity_Parameter.hxx>
+
 #define DEB_CHANNEL 1
 #ifdef DEB_CHANNEL
 #include <BRepTools.hxx>
@@ -70,7 +75,6 @@
 #define _DEVDEBUG_
 #include "HYDRO_trace.hxx"
 
-IMPLEMENT_STANDARD_HANDLE(HYDROData_Channel,HYDROData_ArtificialObject)
 IMPLEMENT_STANDARD_RTTIEXT(HYDROData_Channel,HYDROData_ArtificialObject)
 
 
@@ -84,7 +88,7 @@ HYDROData_Channel::~HYDROData_Channel()
 }
 
 QStringList HYDROData_Channel::DumpToPython( const QString& thePyScriptPath,
-                                             MapOfTreatedObjects& theTreatedObjects ) const
+  MapOfTreatedObjects& theTreatedObjects ) const
 {
   QStringList aResList = dumpObjectCreation( theTreatedObjects );
   QString aName = GetObjPyName();
@@ -92,8 +96,25 @@ QStringList HYDROData_Channel::DumpToPython( const QString& thePyScriptPath,
   Handle(HYDROData_Polyline3D) aRefGideLine = GetGuideLine();
   setPythonReferenceObject( thePyScriptPath, theTreatedObjects, aResList, aRefGideLine, "SetGuideLine" );
 
-  Handle(HYDROData_Profile) aRefProfile = GetProfile();
-  setPythonReferenceObject( thePyScriptPath, theTreatedObjects, aResList, aRefProfile, "SetProfile" );
+  bool mode = GetProfileMode();
+  QString aMode = mode ? "True" : "False" ;
+  aResList << QString( "%1.SetProfileMode( %2 )" ).arg( aName ).arg( aMode );
+
+  if (mode)
+  {
+    Handle(HYDROData_Profile) aRefProfile = GetProfile();
+    setPythonReferenceObject( thePyScriptPath, theTreatedObjects, aResList, aRefProfile, "SetProfile" );
+  }
+  else
+  {
+    QString aLC = QString::number( GetLCValue(), 'f', 3 );
+    QString aDeltaZ = QString::number( GetDeltaZValue(), 'f', 3 );
+    QString aCoteZ = QString::number( GetCoteZValue(), 'f', 3 );
+
+    aResList << QString( "%1.SetLCValue( %2 )" ).arg( aName ).arg( aLC );
+    aResList << QString( "%1.SetDeltaZValue( %2 )" ).arg( aName ).arg( aDeltaZ );
+    aResList << QString( "%1.SetCoteZValue( %2 )" ).arg( aName ).arg( aCoteZ );
+  }
 
   aResList << QString( "%1.SetEquiDistance( %2 )" ).arg( aName ).arg( GetEquiDistance() );
 
@@ -102,6 +123,7 @@ QStringList HYDROData_Channel::DumpToPython( const QString& thePyScriptPath,
   aResList << QString( "" );
 
   return aResList;
+
 }
 
 HYDROData_SequenceOfObjects HYDROData_Channel::GetAllReferenceObjects() const
@@ -122,15 +144,63 @@ HYDROData_SequenceOfObjects HYDROData_Channel::GetAllReferenceObjects() const
 bool HYDROData_Channel::CreatePresentations( const Handle(HYDROData_Polyline3D)& theGuideLine,
                                              const Handle(HYDROData_Profile)&    theProfile,
                                              PrsDefinition&                      thePrs,
-                                             double                              theEquiDistance )
+                                             double                              theEquiDistance,
+                                             bool                                ReverseXCoord)
+{
+  return internalCreatePresentations( true, theGuideLine, theProfile, TopoDS_Wire(), gp_Pnt(), thePrs, theEquiDistance, ReverseXCoord);
+}
+
+bool HYDROData_Channel::CreatePresentations( const Handle(HYDROData_Polyline3D)& theGuideLine,
+                                             double                              LC,
+                                             double                              deltaZ,
+                                             double                              coteZ,
+                                             PrsDefinition&                      thePrs,
+                                             double                              theEquiDistance,
+                                             bool                                ReverseXCoord)
+{
+  Handle(HYDROData_Profile) theProfile;
+  TopoDS_Wire W;
+  gp_Pnt MP;
+
+  gp_Pnt2d A1(-0.75*LC, coteZ-LC/2.0), A2(0.75*LC, coteZ -LC/2.0), B1(-LC/2.0, coteZ), B2(LC/2.0, coteZ), C(0, coteZ+deltaZ);
+  gp_Pnt A1_3d(A1.X(), 0, A1.Y()), A2_3d(A2.X(), 0, A2.Y()), B1_3d(B1.X(), 0, B1.Y()), B2_3d(B2.X(), 0, B2.Y()), C_3d(C.X(), 0, C.Y());
+  BRepLib_MakePolygon PM;
+  PM.Add(A1_3d);
+  PM.Add(B1_3d);
+  PM.Add(C_3d);
+  PM.Add(B2_3d);
+  PM.Add(A2_3d);
+  W = PM.Wire();
+  return internalCreatePresentations( false, theGuideLine, theProfile, W, MP, thePrs, theEquiDistance, ReverseXCoord);
+}
+
+bool HYDROData_Channel::internalCreatePresentations( bool mode,
+                                                     const Handle(HYDROData_Polyline3D)& theGuideLine,
+                                                     const Handle(HYDROData_Profile)&    theProfile,
+                                                     const TopoDS_Wire&                  theProfWire,
+                                                     const gp_Pnt&                       theMiddlePnt,
+                                                     PrsDefinition&                      thePrs,
+                                                     double                              theEquiDistance,
+                                                     bool                                ReverseXCoord)
 {
   // Check input parameters
-  if ( theGuideLine.IsNull() || theProfile.IsNull() ) {
+  if ( theGuideLine.IsNull() )
+  {
     return false;
   }
+  if (mode && theProfile.IsNull() || !mode && theProfWire.IsNull() )
+    return false;
 
   TopoDS_Wire aPathWire = TopoDS::Wire( theGuideLine->GetShape3D() );
-  TopoDS_Wire aProfileWire = TopoDS::Wire( theProfile->GetShape3D() );
+  TopoDS_Wire aProfileWire = theProfWire;
+  //ignore ReverseXCoord if mode is false
+  if (mode)
+  {
+    if (!ReverseXCoord)
+      aProfileWire = TopoDS::Wire( theProfile->GetShape3D(false, false) ); //temp force rebuild
+    else
+      aProfileWire = TopoDS::Wire( theProfile->GetShape3D(true, true));
+  }
   if ( aPathWire.IsNull() || aProfileWire.IsNull() ) {
     return false;
   }
@@ -223,7 +293,11 @@ bool HYDROData_Channel::CreatePresentations( const Handle(HYDROData_Polyline3D)&
   }
 
   // Get the profile middle point ( 3D )
-  gp_Pnt aMiddlePoint( theProfile->GetMiddlePoint( true ) );
+  gp_Pnt aMiddlePoint = theMiddlePnt;
+  if (mode)
+  {
+    aMiddlePoint = theProfile->GetMiddlePoint( true );
+  }
 
   // Translate the profile to each point on the guide line ( 3D )
   Handle(TColgp_HArray1OfPnt) anArrayOfFPnt = new TColgp_HArray1OfPnt(1, aNbPoints );
@@ -307,12 +381,30 @@ void HYDROData_Channel::Update()
   HYDROData_ArtificialObject::Update();
 
   Handle(HYDROData_Polyline3D) aGuideLine = GetGuideLine();
-  Handle(HYDROData_Profile) aProfile = GetProfile();
 
   PrsDefinition aResultPrs;
   double anEquiDistance = GetEquiDistance();
-  if ( !CreatePresentations( aGuideLine, aProfile, aResultPrs, anEquiDistance ) )
-    return;
+
+  bool invDirection = false;
+  Handle(HYDROData_IAltitudeObject) anObjAltitude = GetAltitudeObject();
+  Handle(HYDROData_ChannelAltitude) aChannelAlt = Handle(HYDROData_ChannelAltitude)::DownCast(anObjAltitude);
+  if (!aChannelAlt.IsNull())
+    invDirection = aChannelAlt->GetInvertDirection();
+
+  if (GetProfileMode())
+  {
+    Handle(HYDROData_Profile) aProfile = GetProfile();
+    if ( !CreatePresentations( aGuideLine, aProfile, aResultPrs, anEquiDistance, invDirection ) )
+      return;
+  }
+  else
+  {
+    double lc = GetLCValue();
+    double deltaz = GetDeltaZValue();
+    double cotez = GetCoteZValue();
+    if ( !CreatePresentations( aGuideLine, lc, deltaz, cotez, aResultPrs, anEquiDistance, invDirection ) )
+      return;
+  }
 
   SetShape3D( aResultPrs.myPrs3D );
   SetTopShape( aResultPrs.myPrs2D );
@@ -481,3 +573,43 @@ double HYDROData_Channel::GetEquiDistance() const
 {
   return GetDouble( DataTag_EquiDistance, 1.0 );
 }
+
+void HYDROData_Channel::SetLCValue( double val )
+{
+  SetDouble( DataTag_LC, val );
+}
+
+double HYDROData_Channel::GetLCValue() const
+{
+  return GetDouble( DataTag_LC, 1.0 );
+}
+
+void HYDROData_Channel::SetDeltaZValue( double val )
+{
+  SetDouble( DataTag_DeltaZ, val );
+}
+
+double HYDROData_Channel::GetDeltaZValue() const
+{
+  return GetDouble( DataTag_DeltaZ, 1.0 );
+}
+
+void HYDROData_Channel::SetCoteZValue( double val )
+{
+  SetDouble( DataTag_CoteZ, val );
+}
+
+double HYDROData_Channel::GetCoteZValue() const
+{
+  return GetDouble( DataTag_CoteZ, 1.0 );
+}
+
+void HYDROData_Channel::SetProfileMode( bool mode )
+{
+  SetInteger( DataTag_ProfileMode, (bool)mode );
+}
+
+bool HYDROData_Channel::GetProfileMode() const
+{
+  return (bool)GetInteger( DataTag_ProfileMode, 1 );
+}