src/HYDROData/HYDROData_PolylineXY.cxx

   1 // Copyright (C) 2014-2015  EDF-R&D
   2 // This library is free software; you can redistribute it and/or
   3 // modify it under the terms of the GNU Lesser General Public
   4 // License as published by the Free Software Foundation; either
   5 // version 2.1 of the License, or (at your option) any later version.
   6 //
   7 // This library is distributed in the hope that it will be useful,
   8 // but WITHOUT ANY WARRANTY; without even the implied warranty of
   9 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  10 // Lesser General Public License for more details.
  11 //
  12 // You should have received a copy of the GNU Lesser General Public
  13 // License along with this library; if not, write to the Free Software
  14 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  15 //
  16 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  17 //
  18
  19 #include "HYDROData_PolylineXY.h"
  20
  21 #include "HYDROData_BSplineOperation.h"
  22 #include "HYDROData_Document.h"
  23 #include "HYDROData_ShapesTool.h"
  24 #include "HYDROData_Tool.h"
  25 #include "HYDROData_PolylineOperator.h"
  26 #include "HYDROData_TopoCurve.h"
  27
  28 #include <BRep_Builder.hxx>
  29 #include <BRepBuilderAPI_MakeEdge.hxx>
  30 #include <BRepBuilderAPI_MakeWire.hxx>
  31 #include <BRepBuilderAPI_MakePolygon.hxx>
  32 #include <BRepBuilderAPI_MakeFace.hxx>
  33 #include <BRepOffsetAPI_NormalProjection.hxx>
  34 #include <BRepAdaptor_Curve.hxx>
  35
  36 #ifndef LIGHT_MODE
  37 #include <GEOMBase.h>
  38 #endif
  39
  40 #include <GeomAPI_ProjectPointOnCurve.hxx>
  41 #include <GeomAdaptor_Curve.hxx>
  42 #include <Geom_Line.hxx>
  43 #include <Geom_BSplineCurve.hxx>
  44
  45 #include <GCPnts_AbscissaPoint.hxx>
  46 #include <GCPnts_QuasiUniformDeflection.hxx>
  47 #include <GCPnts_UniformDeflection.hxx>
  48
  49 #include <ImageComposer_MetaTypes.h>
  50
  51 #include <gp_Pnt.hxx>
  52 #include <gp_XY.hxx>
  53 #include <gp_Pln.hxx>
  54
  55 #include <NCollection_Map.hxx>
  56
  57 #include <TCollection_ExtendedString.hxx>
  58
  59 #include <TDataStd_ListIteratorOfListOfByte.hxx>
  60 #include <TColStd_ListIteratorOfListOfInteger.hxx>
  61 #include <TColStd_ListIteratorOfListOfReal.hxx>
  62
  63 #include <TColStd_Array1OfReal.hxx>
  64 #include <TColgp_Array1OfPnt.hxx>
  65
  66 #include <TDataStd_AsciiString.hxx>
  67 #include <TDataStd_BooleanList.hxx>
  68 #include <TDataStd_ExtStringList.hxx>
  69 #include <TDataStd_IntegerList.hxx>
  70 #include <TDataStd_ListIteratorOfListOfExtendedString.hxx>
  71 #include <TDataStd_RealList.hxx>
  72 #include <TDataStd_UAttribute.hxx>
  73
  74 #include <TopoDS_Iterator.hxx>
  75 #include <TopTools_ListIteratorOfListOfShape.hxx>
  76 #include <TopTools_HSequenceOfShape.hxx>
  77 #include <TopExp_Explorer.hxx>
  78 #include <ShapeAnalysis_FreeBounds.hxx>
  79 #include <TopoDS.hxx>
  80
  81 #include <QColor>
  82 #include <QPainterPath>
  83 #include <QVariant>
  84
  85 #define _DEVDEBUG_
  86 #include "HYDRO_trace.hxx"
  87 #include <BRepTools.hxx>
  88 #include <sstream>
  89 #include <cmath>
  90
  91 static const Standard_GUID GUID_IS_UNEDITABLE("e5799736-9030-4051-91a4-2e58321fa153");
  92
  93 const double LOCAL_SELECTION_TOLERANCE = 0.0001;
  94
  95 TCollection_AsciiString getUniqueSectionName( const NCollection_Sequence<TCollection_AsciiString>& theNamesSeq )
  96 {
  97   NCollection_Map<TCollection_AsciiString> aNamesMap;
  98
  99   for ( int i = 1, n = theNamesSeq.Size(); i <= n; ++i )
 100   {
 101     const TCollection_AsciiString& aSectName = theNamesSeq.Value( i );
 102     aNamesMap.Add( aSectName );
 103   }
 104
 105   TCollection_AsciiString aResName;
 106
 107   int aPrefIdx = 1;
 108   do
 109   {
 110     aResName = TCollection_AsciiString( "Section_" ) + aPrefIdx;
 111     ++aPrefIdx;
 112   }
 113   while ( aNamesMap.Contains( aResName ) );
 114
 115   return aResName;
 116 }
 117
 118 TCollection_AsciiString getUniqueSectionName( const Handle(TDataStd_ExtStringList)& theNamesList )
 119 {
 120   NCollection_Sequence<TCollection_AsciiString> aNamesSeq;
 121
 122   TDataStd_ListIteratorOfListOfExtendedString aNamesIter( theNamesList->List() );
 123   for ( ; aNamesIter.More(); aNamesIter.Next() )
 124     aNamesSeq.Append( aNamesIter.Value() );
 125
 126   return getUniqueSectionName( aNamesSeq );
 127 }
 128
 129 IMPLEMENT_STANDARD_HANDLE(HYDROData_PolylineXY, HYDROData_IPolyline)
 130 IMPLEMENT_STANDARD_RTTIEXT(HYDROData_PolylineXY, HYDROData_IPolyline)
 131
 132 HYDROData_PolylineXY::HYDROData_PolylineXY()
 133 : HYDROData_IPolyline(),
 134   myIsInCustomFlag( false )
 135 {
 136 }
 137
 138 HYDROData_PolylineXY::~HYDROData_PolylineXY()
 139 {
 140 }
 141
 142 QStringList HYDROData_PolylineXY::DumpToPython( const QString&       thePyScriptPath,
 143                                                 MapOfTreatedObjects& theTreatedObjects ) const
 144 {
 145   QStringList aResList = dumpObjectCreation( theTreatedObjects );
 146   QString aPolylineName = GetObjPyName();
 147
 148   // Set the wire color
 149   QStringList aWireColorDef;
 150
 151   QColor aWireColor = GetWireColor();
 152   setPythonObjectColor( aWireColorDef, aWireColor, DefaultWireColor(), "SetWireColor" );
 153
 154   if ( !aWireColorDef.isEmpty() )
 155   {
 156     aResList << aWireColorDef;
 157     aResList << QString( "" );
 158   }
 159
 160   bool anIsEditable = IsEditable();
 161   if ( !anIsEditable )
 162   {
 163     // If polyline is not editable we try to import the shape from geom
 164     TCollection_AsciiString aGeomObjectEntry = GetGeomObjectEntry();
 165     if ( !aGeomObjectEntry.IsEmpty() )
 166     {
 167       QString aSalomeObjName = HYDROData_Tool::GenerateNameForPython( theTreatedObjects, "polyline_sobj" );
 168       aResList << QString( "%1 = salome.myStudy.FindObjectID( \"%2\" )" )
 169                   .arg( aSalomeObjName ).arg( aGeomObjectEntry.ToCString() );
 170
 171       aResList << QString( "%1.ImportFromGeomIOR( %2.GetIOR() )" )
 172                   .arg( aPolylineName ).arg( aSalomeObjName );
 173
 174       aResList << QString( "%1.SetGeomObjectEntry( \"%2\" )" )
 175                   .arg( aPolylineName ).arg( aGeomObjectEntry.ToCString() );
 176     }
 177   }
 178   else
 179   {
 180     // Set polyline data
 181     NCollection_Sequence<TCollection_AsciiString>           aSectNames;
 182     NCollection_Sequence<HYDROData_PolylineXY::SectionType> aSectTypes;
 183     NCollection_Sequence<bool>                              aSectClosures;
 184     GetSections( aSectNames, aSectTypes, aSectClosures );
 185
 186     for ( int i = 1, n = aSectNames.Size(); i <= n; ++i )
 187     {
 188       const TCollection_AsciiString& aSectName = aSectNames.Value( i );
 189       const SectionType& aSectType = aSectTypes.Value( i );
 190       bool aSectClosure = aSectClosures.Value( i );
 191
 192       aResList << QString( "%1.AddSection( \"%2\", %3, %4 )" ).arg( aPolylineName )
 193                   .arg( aSectName.ToCString() ).arg( aSectType ).arg( aSectClosure );
 194
 195       HYDROData_IPolyline::PointsList aSectPointsList = GetPoints( i - 1 );
 196       for ( int k = 1, aNbPoints = aSectPointsList.Size(); k <= aNbPoints; ++k )
 197       {
 198         const Point& aSectPoint = aSectPointsList.Value( k );
 199
 200         QString anXStr = QString::number( aSectPoint.X(), 'f', 2 );
 201         QString anYStr = QString::number( aSectPoint.Y(), 'f', 2 );
 202         aResList << QString( "%1.AddPoint( %2, gp_XY( %3, %4 ) )" ).arg( aPolylineName )
 203           .arg( i - 1 ).arg( anXStr ).arg( anYStr );
 204       }
 205     }
 206   }
 207   aResList << QString( "" );
 208   aResList << QString( "%1.Update()" ).arg( aPolylineName );
 209   aResList << QString( "" );
 210
 211   return aResList;
 212 }
 213
 214 QVariant HYDROData_PolylineXY::GetDataVariant()
 215 {
 216   QPainterPath aPath = GetPainterPath();
 217
 218   QVariant aVarData;
 219   aVarData.setValue<QPainterPath>( aPath );
 220
 221   return aVarData;
 222 }
 223
 224 QColor HYDROData_PolylineXY::DefaultWireColor()
 225 {
 226   return QColor( Qt::black );
 227 }
 228
 229 bool HYDROData_PolylineXY::ImportFromGeomIOR( const TCollection_AsciiString& theIOR )
 230 {
 231 #ifdef LIGHT_MODE
 232   return false;
 233 #else
 234   if ( theIOR.IsEmpty() )
 235     return false;
 236
 237   TopoDS_Shape aShape = GEOMBase::GetShapeFromIOR( theIOR.ToCString() );
 238   if ( aShape.IsNull() )
 239     return false;
 240
 241   return ImportShape( aShape, false, NULL );
 242 #endif
 243 }
 244
 245 void HYDROData_PolylineXY::SetGeomObjectEntry( const TCollection_AsciiString& theEntry )
 246 {
 247   TDataStd_AsciiString::Set( myLab.FindChild( DataTag_GeomObjectEntry ), theEntry );
 248 }
 249
 250 TCollection_AsciiString HYDROData_PolylineXY::GetGeomObjectEntry() const
 251 {
 252   TCollection_AsciiString aRes;
 253
 254   TDF_Label aLabel = myLab.FindChild( DataTag_GeomObjectEntry, false );
 255   if ( !aLabel.IsNull() )
 256   {
 257     Handle(TDataStd_AsciiString) anAsciiStr;
 258     if ( aLabel.FindAttribute( TDataStd_AsciiString::GetID(), anAsciiStr ) )
 259       aRes = anAsciiStr->Get();
 260   }
 261
 262   return aRes;
 263 }
 264
 265 bool convertEdgeToSection( const TopoDS_Edge&                                       theEdge,
 266                            NCollection_Sequence<TCollection_AsciiString>&           theSectNames,
 267                            NCollection_Sequence<HYDROData_PolylineXY::SectionType>& theSectTypes,
 268                            NCollection_Sequence<bool>&                              theSectClosures,
 269                            NCollection_Sequence<HYDROData_PolylineXY::PointsList>&  theSectPoints,
 270                            bool                                                     IsCanBeClosed,
 271                            bool                                                     IsInterpolationAllowed,
 272                            double                                                   theDeflection )
 273 {
 274   Standard_Real aFirst = 0.0, aLast = 0.0;
 275   Handle(Geom_Curve) anEdgeGeomCurve = BRep_Tool::Curve( theEdge, aFirst, aLast );
 276   if ( anEdgeGeomCurve.IsNull() )
 277     return false;
 278
 279   TCollection_AsciiString aSectName = getUniqueSectionName( theSectNames );
 280   bool anIsEdgeClosed = anEdgeGeomCurve->IsClosed();
 281
 282   HYDROData_PolylineXY::SectionType aSectionType = HYDROData_PolylineXY::SECTION_POLYLINE;
 283   HYDROData_PolylineXY::PointsList aPointsList;
 284
 285   if( anEdgeGeomCurve->IsKind( STANDARD_TYPE(Geom_Line) ) )
 286   {
 287     Handle(Geom_Line) aGeomLine = Handle(Geom_Line)::DownCast( anEdgeGeomCurve );
 288
 289     gp_Pnt aFirstPoint, aLastPoint;
 290     aGeomLine->D0( aFirst, aFirstPoint );
 291     aGeomLine->D0( aLast, aLastPoint );
 292
 293     HYDROData_PolylineXY::Point aSectFirstPoint( aFirstPoint.X(), aFirstPoint.Y() );
 294     aPointsList.Append( aSectFirstPoint );
 295
 296     HYDROData_PolylineXY::Point aSectLastPoint( aLastPoint.X(), aLastPoint.Y() );
 297     aPointsList.Append( aSectLastPoint );
 298   }
 299   else if ( anEdgeGeomCurve->IsKind( STANDARD_TYPE(Geom_BSplineCurve) ) || IsInterpolationAllowed )
 300   {
 301     aSectionType = HYDROData_PolylineXY::SECTION_SPLINE;
 302
 303     BRepAdaptor_Curve anAdaptorCurve( theEdge );
 304     GCPnts_QuasiUniformDeflection aDiscrete( anAdaptorCurve, theDeflection );
 305     //GCPnts_UniformDeflection aDiscrete( anAdaptorCurve, theDeflection );
 306
 307     int aNbPoints = aDiscrete.NbPoints();
 308
 309     // Decrease the number of imported poles because of last one
 310     // pole is the closing point which are the start point
 311     if ( anIsEdgeClosed ) aNbPoints--;
 312
 313     for ( int i = 1; i <= aNbPoints; ++i )
 314     {
 315       const gp_Pnt& aPoint = aDiscrete.Value( i );
 316
 317       HYDROData_PolylineXY::Point aSectPoint( aPoint.X(), aPoint.Y() );
 318       aPointsList.Append( aSectPoint );
 319     }
 320   }
 321   else
 322   {
 323     // Other curve types are not supported
 324     return false;
 325   }
 326
 327   if ( aPointsList.IsEmpty() )
 328     return false;
 329
 330   theSectNames.Append( aSectName );
 331   theSectTypes.Append( aSectionType );
 332   theSectClosures.Append( anIsEdgeClosed );
 333   theSectPoints.Append( aPointsList );
 334
 335   return true;
 336 }
 337
 338
 339 bool convertEdgesToSections( const TopoDS_Edge&                                       theEdge,
 340                              NCollection_Sequence<TCollection_AsciiString>&           theSectNames,
 341                              NCollection_Sequence<HYDROData_PolylineXY::SectionType>& theSectTypes,
 342                              NCollection_Sequence<bool>&                              theSectClosures,
 343                              NCollection_Sequence<HYDROData_PolylineXY::PointsList>&  theSectPoints,
 344                              bool                                                     IsCanBeClosed,
 345                              bool                                                     IsInterpolationAllowed,
 346                              double                                                   theDeflection,
 347                              const Handle( HYDROData_PolylineXY )&                    theOldPolyline )
 348 {
 349   DEBTRACE("convertEdgesToSections")
 350   Standard_Real aFirst = 0.0, aLast = 0.0;
 351   Handle(Geom_Curve) anEdgeGeomCurve = BRep_Tool::Curve( theEdge, aFirst, aLast );
 352   if ( anEdgeGeomCurve.IsNull() )
 353     return false;
 354
 355   bool isPrevious = (theSectTypes.Size() > 0);
 356   DEBTRACE("nb sections: " << theSectTypes.Size());
 357   HYDROData_PolylineXY::SectionType prevSectType = HYDROData_PolylineXY::SECTION_SPLINE;;
 358   HYDROData_PolylineXY::PointsList prevPointList;
 359   bool isPrevClosed = true;
 360   if (isPrevious)
 361     {
 362        prevSectType = theSectTypes.Last();
 363        prevPointList = theSectPoints.Last();
 364        isPrevClosed = theSectClosures.Last();
 365     }
 366
 367   bool anIsEdgeClosed = anEdgeGeomCurve->IsClosed();
 368   anIsEdgeClosed &= IsCanBeClosed; // on split, resulting edges are normally not closed...
 369   HYDROData_PolylineXY::SectionType aSectionType ;
 370   if( anEdgeGeomCurve->IsKind( STANDARD_TYPE(Geom_Line) ) )
 371     {
 372       aSectionType = HYDROData_PolylineXY::SECTION_POLYLINE;
 373     }
 374   else if ( anEdgeGeomCurve->IsKind( STANDARD_TYPE(Geom_BSplineCurve) ) || IsInterpolationAllowed )
 375     {
 376       aSectionType = HYDROData_PolylineXY::SECTION_SPLINE;
 377     }
 378   else
 379     {
 380       // Other curve types are not supported
 381       return false;
 382     }
 383
 384   bool isNewSection = !isPrevious || isPrevClosed || anIsEdgeClosed || prevSectType != aSectionType;
 385   DEBTRACE(isNewSection <<": " << !isPrevious << " " << isPrevClosed << " " << anIsEdgeClosed << " " << (prevSectType != aSectionType));
 386
 387   HYDROData_PolylineXY::PointsList aPointsList;
 388   if (!isNewSection)
 389     {
 390       aPointsList = prevPointList;
 391     }
 392
 393   if( aSectionType == HYDROData_PolylineXY::SECTION_POLYLINE )
 394     {
 395       DEBTRACE("SECTION_POLYLINE");
 396       Handle(Geom_Line) aGeomLine = Handle(Geom_Line)::DownCast( anEdgeGeomCurve );
 397
 398       gp_Pnt aFirstPoint, aLastPoint;
 399       aGeomLine->D0( aFirst, aFirstPoint );
 400       aGeomLine->D0( aLast, aLastPoint );
 401       HYDROData_PolylineXY::Point aSectFirstPoint( aFirstPoint.X(), aFirstPoint.Y() );
 402       HYDROData_PolylineXY::Point aSectLastPoint( aLastPoint.X(), aLastPoint.Y() );
 403      if (!isNewSection)
 404         {
 405           if (aSectFirstPoint == prevPointList.Last())
 406             {
 407               DEBTRACE("points shared: a");//aPointsList.Append( aSectFirstPoint );
 408               aPointsList.Append( aSectLastPoint );
 409             }
 410           else if (aSectLastPoint == prevPointList.Last())
 411             {
 412               DEBTRACE("points shared: b");//aPointsList.Append( aSectLastPoint );
 413               aPointsList.Append( aSectFirstPoint );
 414             }
 415           else if (aSectFirstPoint == prevPointList.First())
 416             {
 417               DEBTRACE("points shared: c");//aPointsList.Prepend( aSectFirstPoint );
 418               aPointsList.Prepend( aSectLastPoint );
 419             }
 420           else if (aSectLastPoint == prevPointList.First())
 421             {
 422               DEBTRACE("points shared: d");//aPointsList.Prepend( aSectLastPoint );
 423               aPointsList.Prepend( aSectFirstPoint );
 424             }
 425           else
 426             {
 427               DEBTRACE("no point shared")
 428               isNewSection = true; // no point shared, new section
 429               aPointsList.Clear();
 430               aPointsList.Append( aSectFirstPoint );
 431               aPointsList.Append( aSectLastPoint );
 432             }
 433         }
 434      else
 435        {
 436          DEBTRACE("new section");
 437          aPointsList.Append( aSectFirstPoint );
 438          aPointsList.Append( aSectLastPoint );
 439        }
 440     }
 441   else // aSectionType == HYDROData_PolylineXY::SECTION_SPLINE
 442     {
 443       DEBTRACE("SECTION_SPLINE");
 444       isNewSection = true;
 445       aPointsList.Clear();
 446
 447       BRepAdaptor_Curve anAdaptorCurve(theEdge);
 448       if (theOldPolyline.IsNull()) // --- no previous polyline: build a set of points from scratch for the spline
 449         {
 450           GCPnts_QuasiUniformDeflection aDiscrete(anAdaptorCurve, theDeflection);
 451
 452           int aNbPoints = aDiscrete.NbPoints();
 453
 454           // Decrease the number of imported poles because of last one
 455           // pole is the closing point which are the start point
 456           if (anIsEdgeClosed)
 457             aNbPoints--;
 458
 459           for (int i = 1; i <= aNbPoints; ++i)
 460             {
 461               const gp_Pnt& aPoint = aDiscrete.Value(i);
 462               HYDROData_PolylineXY::Point aSectPoint(aPoint.X(), aPoint.Y());
 463               aPointsList.Append(aSectPoint);
 464             }
 465
 466           gp_Pnt endPts[] = {gp_Pnt(aPointsList.First().X(),aPointsList.First().Y(), 0),
 467                              gp_Pnt(aPointsList.Last().X(),aPointsList.Last().Y(), 0) };
 468           DEBTRACE("curve start: "<< endPts[0].X() << " " << endPts[0].Y());
 469           DEBTRACE("curve end: "<< endPts[1].X() << " " << endPts[1].Y());
 470         }
 471       else // --- split of a previous polyline: try to retrieve old sets of points and add intersection points
 472         {
 473           const gp_Pnt aEndPs[] = { anAdaptorCurve.Value(anAdaptorCurve.FirstParameter()).XYZ(),
 474                                     anAdaptorCurve.Value(anAdaptorCurve.LastParameter()).XYZ() };
 475           double midPar = (anAdaptorCurve.LastParameter() + anAdaptorCurve.FirstParameter())/2;
 476           gp_Pnt midPnt;
 477           anAdaptorCurve.D0(midPar, midPnt);
 478           DEBTRACE("curve first point: " << aEndPs[0].X() << " " << aEndPs[0].Y() << " " << aEndPs[0].Z());
 479           DEBTRACE("curve last point: " << aEndPs[1].X() << " " << aEndPs[1].Y() << " " << aEndPs[1].Z());
 480           DEBTRACE("curve mid point: " << midPnt.X() << " " << midPnt.Y() << " " << midPnt.Z());
 481
 482           std::vector<TopoDS_Wire> aCurves;
 483           HYDROData_PolylineOperator::GetWires(theOldPolyline, aCurves);
 484
 485           int nbSections = theOldPolyline->NbSections();
 486           DEBTRACE("nbSections: "<< nbSections << ", nbCurves: " << aCurves.size() );
 487           for (int isec = 0; isec < nbSections; isec++)
 488             {
 489               DEBTRACE("section: "<< isec);
 490               bool isOldSectionclosed = theOldPolyline->IsClosedSection(isec);
 491               TopoDS_Wire aWire = aCurves[isec]; // we suppose sections and wires are in the same order
 492               TopExp_Explorer anExp(aWire, TopAbs_EDGE);
 493               TopoDS_Edge anEdge = TopoDS::Edge(anExp.Current()); // the first is OK: only one normally with splines
 494               BRepAdaptor_Curve adaptorOldCurve(anEdge);
 495               double pfirst = adaptorOldCurve.FirstParameter();
 496               double plast = adaptorOldCurve.LastParameter();
 497               DEBTRACE("previous curve first last : "<< pfirst << " " << plast);
 498               double p[3] = {-1, -1};
 499               double d0= ProjectPointToCurve(aEndPs[0].XYZ(), adaptorOldCurve,p[0]);
 500               double d1= ProjectPointToCurve(aEndPs[1].XYZ(), adaptorOldCurve,p[1]);
 501               double d2= ProjectPointToCurve(midPnt.XYZ(), adaptorOldCurve, p[2]);
 502               DEBTRACE("d0: "<<d0<<" d1: "<<d1<<" d2: "<<d2<<" p0: "<<p[0]<<" p1: "<<p[1]<<" p2: "<<p[2]);
 503               if ((d0 < 1.e-3) and (d1 < 1.e-3) and (d2 < 1.e-3)) // we got the good old curve (and the good section)
 504                 {
 505                   double pmin = p[0];
 506                   double pmax = p[1];
 507                   bool forward = true;
 508                   if (! isOldSectionclosed) // no need to check first and last points on an open curve
 509                     {
 510                       pmin = p[1];
 511                       pmax = p[0];
 512                       forward = false;
 513                     }
 514                   else // old section closed: check if we use first or last points...
 515                     {
 516                       if((pmin < pmax) and((abs(pmin - pfirst)> 1.e-3) and (abs(pmax - plast) >1.e-3))) // internal points forward
 517                         forward = true;
 518                       else if ((pmin > pmax) and((abs(pmin - plast)> 1.e-3) and (abs(pmax - pfirst) >1.e-3))) // internal points reverse
 519                         {
 520                           pmin = p[1];
 521                           pmax = p[0];
 522                           forward = false;
 523                         }
 524                       else if ((abs(pmin - plast) <1.e-3) and (p[2] < pmax)) // forward, replace pmin par pfirst
 525                           pmin = pfirst;
 526                       else if ((abs(pmin - plast) <1.e-3) and (p[2] > pmax)) // reverse
 527                         {
 528                           pmin = p[1];
 529                           pmax = p[0];
 530                           forward = false;
 531                         }
 532                       else if ((abs(pmax - pfirst) <1.e-3) and (p[2] < pmin)) // reverse
 533                         {
 534                           pmin = p[1];
 535                           pmax = p[0];
 536                           forward = false;
 537                         }
 538                       else if ((abs(pmax - pfirst) <1.e-3) and (p[2] > pmin)) // forward, replace pmax par plast
 539                         pmax = plast;
 540                    }
 541                   HYDROData_PolylineXY::Point aFirstPoint, aLastPoint;
 542                   if (forward)
 543                     {
 544                       aFirstPoint = HYDROData_PolylineXY::Point(aEndPs[0].X(), aEndPs[0].Y());
 545                       aLastPoint = HYDROData_PolylineXY::Point(aEndPs[1].X(), aEndPs[1].Y());
 546                     }
 547                   else
 548                     {
 549                       aFirstPoint = HYDROData_PolylineXY::Point(aEndPs[1].X(), aEndPs[1].Y());
 550                       aLastPoint = HYDROData_PolylineXY::Point(aEndPs[0].X(), aEndPs[0].Y());
 551                     }
 552                   aPointsList.Append(aFirstPoint);
 553
 554                   HYDROData_PolylineXY::PointsList aSectPoints = theOldPolyline->GetPoints(isec, false);
 555                   int nbPoints = aSectPoints.Length();
 556                   if (forward)
 557                     for (int i=1; i<=nbPoints; i++)
 558                       {
 559                         HYDROData_PolylineXY::Point aPoint = aSectPoints.Value(i);
 560                         gp_XYZ p(aPoint.X(), aPoint.Y(), 0);
 561                         double param =-1;
 562                         double d = ProjectPointToCurve(p, adaptorOldCurve, param);
 563                         if ((param > pmin) and (param < pmax))
 564                           {
 565                             DEBTRACE("param: " << param);
 566                             aPointsList.Append(aPoint);
 567                           }
 568                       }
 569                   else
 570                     for (int i=nbPoints; i>0; i--)
 571                       {
 572                         HYDROData_PolylineXY::Point aPoint = aSectPoints.Value(i);
 573                         gp_XYZ p(aPoint.X(), aPoint.Y(), 0);
 574                         double param =-1;
 575                         double d = ProjectPointToCurve(p, adaptorOldCurve, param);
 576                         if ((param > pmin) and (param < pmax))
 577                           {
 578                             DEBTRACE("param: " << param);
 579                             aPointsList.Append(aPoint);
 580                           }
 581                       }
 582
 583                   aPointsList.Append(aLastPoint);
 584                 }
 585
 586             }
 587         }
 588     }
 589
 590   if ( aPointsList.IsEmpty() )
 591     return false;
 592
 593   TCollection_AsciiString aSectName = getUniqueSectionName( theSectNames );
 594   if (isNewSection)
 595     {
 596       DEBTRACE("isNewSection");
 597       theSectNames.Append( aSectName );
 598       theSectTypes.Append( aSectionType );
 599       theSectClosures.Append( anIsEdgeClosed );
 600       theSectPoints.Append( aPointsList );
 601     }
 602   else
 603     {
 604       DEBTRACE("sameSection");
 605       theSectPoints.SetValue(theSectPoints.Length(), aPointsList);
 606     }
 607
 608   return true;
 609 }
 610
 611 bool HYDROData_PolylineXY::ImportShape( const TopoDS_Shape& theShape,
 612                                         bool IsInterpolationAllowed,
 613                                         const Handle( HYDROData_PolylineXY )& theOldPolyline,
 614                                         bool IsClosureAllowed,
 615                                         double theDeviation )
 616 {
 617   DEBTRACE("ImportShape");
 618   if ( theShape.IsNull() )
 619     return false;
 620
 621   std::string brepName = this->GetName().toStdString();
 622   brepName += ".brep";
 623   BRepTools::Write( theShape, brepName.c_str() );
 624
 625   RemoveSections();
 626
 627   bool anIsCanBeImported = false;
 628
 629   NCollection_Sequence<TCollection_AsciiString> aSectNames;
 630   NCollection_Sequence<SectionType>             aSectTypes;
 631   NCollection_Sequence<bool>                    aSectClosures;
 632   NCollection_Sequence<PointsList>              aSectPoints;
 633
 634   if ( theShape.ShapeType() == TopAbs_EDGE )
 635   {
 636       DEBTRACE("TopAbs_EDGE");
 637     TopoDS_Edge anEdge = TopoDS::Edge( theShape );
 638 //    anIsCanBeImported = convertEdgeToSection( anEdge, aSectNames, aSectTypes,
 639 //      aSectClosures, aSectPoints, true, IsInterpolationAllowed, theDeviation );
 640     anIsCanBeImported = convertEdgesToSections( anEdge, aSectNames, aSectTypes, aSectClosures,
 641                                                 aSectPoints, IsClosureAllowed, IsInterpolationAllowed,
 642                                                 theDeviation, theOldPolyline );
 643   }
 644   else if ( theShape.ShapeType() == TopAbs_WIRE )
 645   {
 646       DEBTRACE("TopAbs_WIRE");
 647     TopTools_SequenceOfShape anEdges;
 648     HYDROData_ShapesTool::ExploreShapeToShapes( theShape, TopAbs_EDGE, anEdges );
 649
 650     anIsCanBeImported = !anEdges.IsEmpty();
 651     for ( int i = 1, n = anEdges.Length(); i <= n && anIsCanBeImported; ++i )
 652     {
 653       TopoDS_Edge aWireEdge = TopoDS::Edge( anEdges.Value( i ) );
 654       anIsCanBeImported = convertEdgesToSections( aWireEdge, aSectNames, aSectTypes, aSectClosures,
 655                                                   aSectPoints, IsClosureAllowed, IsInterpolationAllowed,
 656                                                   theDeviation, theOldPolyline );
 657     }
 658   }
 659
 660   if ( anIsCanBeImported )
 661   {
 662     for ( int i = 1, n = aSectNames.Length(); i <= n; ++i )
 663     {
 664       const TCollection_AsciiString& aSectName = aSectNames.Value( i );
 665       const SectionType& aSectType = aSectTypes.Value( i );
 666       bool anIsSectionClosed = aSectClosures.Value( i );
 667       const PointsList& aSectPointsList = aSectPoints( i );
 668
 669       AddSection( aSectName, aSectType, anIsSectionClosed );
 670       SetPoints( i - 1, aSectPointsList );
 671     }
 672   }
 673   else
 674   {
 675     TopoDS_Shape aShape = theShape;
 676
 677     if ( theShape.ShapeType() == TopAbs_EDGE )
 678     {
 679       // We make the wire from incoming edge because of other algorithms
 680       // are waiting at least the wire from polyline
 681       TopoDS_Edge anEdge = TopoDS::Edge( theShape );
 682       BRepBuilderAPI_MakeWire aMakeWire( anEdge );
 683       aMakeWire.Build();
 684       if ( aMakeWire.IsDone() )
 685         aShape = aMakeWire.Wire();
 686     }
 687
 688     gp_Pln aPlane( gp_Pnt( 0, 0, 0 ), gp_Dir( 0, 0, 1 ) );
 689     BRepBuilderAPI_MakeFace aMakeFace( aPlane );
 690     aMakeFace.Build();
 691     BRepOffsetAPI_NormalProjection aProj( aMakeFace.Face() );
 692     aProj.Add( aShape );
 693     aProj.Build();
 694     TopoDS_Shape aResult;
 695     if( aProj.IsDone() )
 696       aResult = aProj.Shape();
 697
 698     SetShape( aResult );
 699   }
 700
 701   setEditable( anIsCanBeImported );
 702   Update();
 703   return true;
 704 }
 705
 706 TopoDS_Wire HYDROData_PolylineXY::BuildWire( const SectionType&                  theType,
 707                                              const bool&                         theIsClosed,
 708                                              const NCollection_Sequence<gp_XYZ>& thePoints )
 709 {
 710   TopoDS_Wire aWire;
 711   if( theType == SECTION_POLYLINE )
 712   {
 713     int aNbPoints = thePoints.Length();
 714     BRepBuilderAPI_MakePolygon aMakeWire;
 715     for ( int i = 1, n = aNbPoints; i <= n ; ++i )
 716     {
 717       gp_XYZ aPoint = thePoints.Value( i );
 718       gp_Pnt aPnt( aPoint.X(), aPoint.Y(), aPoint.Z() );
 719       aMakeWire.Add( aPnt );
 720     }
 721     if( theIsClosed && ( aNbPoints > 2 ) )
 722       aMakeWire.Close();
 723
 724     if ( aMakeWire.IsDone() )
 725       aWire = aMakeWire.Wire();
 726   }
 727   else //if( theType == PolylineSection::SECTION_SPLINE )
 728   {
 729     BRepBuilderAPI_MakeWire aMakeWire;
 730
 731     if ( thePoints.Size() > 1 )
 732     {
 733       Handle(Geom_BSplineCurve) aCurve =
 734         HYDROData_BSplineOperation::ComputeCurve( thePoints, theIsClosed, LOCAL_SELECTION_TOLERANCE );
 735
 736       TopoDS_Edge anEdge = BRepBuilderAPI_MakeEdge( aCurve ).Edge();
 737       aMakeWire.Add( anEdge );
 738     }
 739     aMakeWire.Build();
 740     if ( aMakeWire.IsDone() )
 741       aWire = aMakeWire;
 742   }
 743
 744   return aWire;
 745 }
 746
 747 void HYDROData_PolylineXY::BuildPainterPath( QPainterPath&                       thePath,
 748                                              const SectionType&                  theType,
 749                                              const bool&                         theIsClosed,
 750                                              const NCollection_Sequence<gp_XYZ>& thePoints )
 751 {
 752   if ( thePoints.IsEmpty() )
 753     return;
 754
 755   if ( theType == SECTION_POLYLINE )
 756   {
 757     const gp_XYZ& aFirstPoint = thePoints.Value( 1 );
 758     thePath.moveTo( aFirstPoint.X(), aFirstPoint.Y() );
 759
 760     for( int i = 2, n = thePoints.Size(); i <= n; ++i )
 761     {
 762       const gp_XYZ& aSectPoint = thePoints.Value( i );
 763
 764       thePath.lineTo( aSectPoint.X(), aSectPoint.Y() );
 765     }
 766
 767     if( theIsClosed )
 768       thePath.closeSubpath();
 769   }
 770   else
 771   {
 772     Handle(Geom_BSplineCurve) aCurve =
 773       HYDROData_BSplineOperation::ComputeCurve( thePoints, theIsClosed, LOCAL_SELECTION_TOLERANCE );
 774     HYDROData_BSplineOperation::ComputePath( aCurve, thePath );
 775   }
 776 }
 777
 778 void HYDROData_PolylineXY::Update()
 779 {
 780   if ( !IsEditable() )
 781   {
 782     // If polyline is not editable we no need to update it wire
 783     ClearChanged();
 784     return;
 785   }
 786
 787   HYDROData_IPolyline::Update();
 788
 789   NCollection_Sequence<TCollection_AsciiString>           aSectNames;
 790   NCollection_Sequence<HYDROData_PolylineXY::SectionType> aSectTypes;
 791   NCollection_Sequence<bool>                              aSectClosures;
 792   GetSections( aSectNames, aSectTypes, aSectClosures );
 793
 794   //BRepBuilderAPI_MakeWire aMakeWire;
 795
 796   TopTools_ListOfShape aSectionWiresList;
 797
 798   for ( int aSectionId = 1, aNbSects = aSectNames.Size(); aSectionId <= aNbSects; aSectionId++ )
 799   {
 800     TCollection_AsciiString aSectName = aSectNames.Value( aSectionId );
 801     SectionType aSectionType = aSectTypes.Value( aSectionId );
 802     bool anIsSectionClosed = aSectClosures.Value( aSectionId );
 803
 804     PointsList aSectPointsList = GetPoints( aSectionId - 1 );
 805     if ( aSectPointsList.IsEmpty() )
 806       continue;
 807
 808     NCollection_Sequence<gp_XYZ> aPoints;
 809     for( int i = 1, n = aSectPointsList.Size(); i <= n; ++i )
 810     {
 811       const Point& aSectPoint = aSectPointsList.Value( i );
 812
 813       gp_XYZ aPoint( aSectPoint.X(), aSectPoint.Y(), 0.0 );
 814       aPoints.Append( aPoint );
 815     }
 816
 817     TopoDS_Wire aSectionWire = BuildWire( aSectionType, anIsSectionClosed, aPoints );
 818     if ( !aSectionWire.IsNull() ) {
 819       aSectionWiresList.Append( aSectionWire );
 820       //aMakeWire.Add( aSectionWire );
 821     }
 822   }
 823 // all input wires in the <aSectionWiresList>
 824   Handle(TopTools_HSequenceOfShape) aSeqWires = new TopTools_HSequenceOfShape;
 825   Handle(TopTools_HSequenceOfShape) aSeqEdges = new TopTools_HSequenceOfShape;
 826   TopTools_ListIteratorOfListOfShape it(aSectionWiresList);
 827   for(;it.More();it.Next())
 828   {
 829     TopExp_Explorer it2(it.Value(), TopAbs_EDGE);
 830     for(;it2.More();it2.Next())
 831       aSeqEdges->Append(it2.Current());
 832   }
 833
 834   BRep_Builder aBB;
 835   TopoDS_Compound aCmp;
 836   TopoDS_Shape aResult;
 837   aBB.MakeCompound(aCmp);
 838   if(aSeqEdges->Length() >1)
 839   {
 840     ShapeAnalysis_FreeBounds::ConnectEdgesToWires( aSeqEdges, 1E-5, Standard_True, aSeqWires );
 841
 842     if( aSeqWires->Length()==1 )
 843       aResult = aSeqWires->Value( 1 );
 844     else
 845     {
 846       for (Standard_Integer i = 1; i <= aSeqWires->Length();i++)
 847       {
 848         const TopoDS_Shape& aS1 = aSeqWires->Value(i);
 849         aBB.Add(aCmp, aS1);
 850       }
 851       aResult = aCmp;
 852     }
 853   }
 854   else if (aSeqEdges->Length() == 1)
 855   {
 856     BRepBuilderAPI_MakeWire mkWire (TopoDS::Edge(aSeqEdges->Value(1)));
 857     if (mkWire.IsDone())
 858       aResult = mkWire.Wire();
 859   }
 860
 861   SetShape( aResult );
 862 }
 863
 864 bool HYDROData_PolylineXY::IsHas2dPrs() const
 865 {
 866   return true;
 867 }
 868
 869 bool HYDROData_PolylineXY::IsEditable() const
 870 {
 871   return !myLab.IsAttribute( GUID_IS_UNEDITABLE );
 872 }
 873
 874 void HYDROData_PolylineXY::setEditable( const bool theIsEditable )
 875 {
 876   if ( !theIsEditable )
 877     TDataStd_UAttribute::Set( myLab, GUID_IS_UNEDITABLE );
 878   else
 879     myLab.ForgetAttribute( GUID_IS_UNEDITABLE );
 880 }
 881
 882 /**
 883  * Returns true if polyline is closed
 884  */
 885 bool HYDROData_PolylineXY::IsClosed( const bool theIsSimpleCheck ) const
 886 {
 887   bool anIsClosed = false;
 888
 889   TopoDS_Shape aShape = GetShape();
 890   if ( aShape.IsNull() )
 891     return anIsClosed;
 892
 893   TopTools_SequenceOfShape aWires;
 894   HYDROData_ShapesTool::ExploreShapeToShapes( aShape, TopAbs_WIRE, aWires );
 895
 896   int aNbWires = aWires.Length();
 897   if ( theIsSimpleCheck )
 898   {
 899     anIsClosed = aNbWires > 0;
 900     for ( int i = 1; i <= aNbWires && anIsClosed; ++i )
 901     {
 902       const TopoDS_Shape& aWire = aWires.Value( i );
 903       anIsClosed = BRep_Tool::IsClosed( aWire );
 904     }
 905   }
 906   else
 907   {
 908     anIsClosed = aNbWires == 1 && BRep_Tool::IsClosed( aWires.First() );
 909   }
 910
 911   return anIsClosed;
 912 }
 913
 914 double HYDROData_PolylineXY::GetDistance( const int theSectionIndex,
 915                                           const int thePointIndex ) const
 916 {
 917   double aResDistance = -1;
 918   if ( theSectionIndex < 0 || theSectionIndex >= NbSections() )
 919     return aResDistance;
 920
 921   if ( thePointIndex == 0 )
 922     return 0.0;
 923
 924   SectionType aSectionType = GetSectionType( theSectionIndex );
 925   bool anIsSectionClosed = IsClosedSection( theSectionIndex );
 926   PointsList aSectPointsList = GetPoints( theSectionIndex );
 927   if ( thePointIndex < 0 || thePointIndex >= aSectPointsList.Size()  )
 928     return aResDistance;
 929
 930   if ( aSectionType == SECTION_POLYLINE )
 931   {
 932     aResDistance = 0.0;
 933
 934     Point aPrevPoint = aSectPointsList.Value( 1 );
 935     for ( int i = 2, aNbPoints = aSectPointsList.Size(); i <= aNbPoints; ++i )
 936     {
 937       const Point& aSectPoint = aSectPointsList.Value( i );
 938       aResDistance += gp_Pnt2d( aPrevPoint ).Distance( aSectPoint );
 939       aPrevPoint = aSectPoint;
 940
 941       if ( thePointIndex == i - 1 )
 942         break;
 943     }
 944   }
 945   else
 946   {
 947     gp_XYZ aPointToTest;
 948
 949     int aSectNbPoints = aSectPointsList.Size();
 950     NCollection_Sequence<gp_XYZ> aPoints;
 951     for( int i = 1 ; i <= aSectNbPoints; ++i )
 952     {
 953       const Point& aSectPoint = aSectPointsList.Value( i );
 954
 955       gp_XYZ aPoint( aSectPoint.X(), aSectPoint.Y(), 0.0 );
 956       aPoints.Append( aPoint );
 957
 958       if ( thePointIndex == i - 1 )
 959         aPointToTest = aPoint;
 960     }
 961
 962     Handle(Geom_BSplineCurve) aCurve =
 963       HYDROData_BSplineOperation::ComputeCurve( aPoints, anIsSectionClosed, LOCAL_SELECTION_TOLERANCE );
 964
 965     Quantity_Parameter aFirstParam = aCurve->FirstParameter();
 966     Quantity_Parameter aSecondParam = aCurve->LastParameter();
 967
 968     if ( thePointIndex != aSectNbPoints - 1 )
 969     {
 970       GeomAPI_ProjectPointOnCurve aProject( aPointToTest, aCurve );
 971       aSecondParam = aProject.LowerDistanceParameter();
 972     }
 973
 974     GeomAdaptor_Curve anAdap( aCurve );
 975
 976     aResDistance = GCPnts_AbscissaPoint::Length( anAdap, aFirstParam, aSecondParam );
 977   }
 978
 979   return aResDistance;
 980 }
 981
 982 int HYDROData_PolylineXY::NbSections() const
 983 {
 984   Handle(TDataStd_ExtStringList) aNamesList;
 985   Handle(TDataStd_IntegerList)   aTypesList;
 986   Handle(TDataStd_BooleanList)   aClosuresList;
 987   getSectionsLists( aNamesList, aTypesList, aClosuresList, false );
 988
 989   return !aClosuresList.IsNull() ? aClosuresList->Extent() : 0;
 990 }
 991
 992 void HYDROData_PolylineXY::AddSection( const TCollection_AsciiString& theSectName,
 993                                        const SectionType              theSectionType,
 994                                        const bool                     theIsClosed )
 995 {
 996   Handle(TDataStd_ExtStringList) aNamesList;
 997   Handle(TDataStd_IntegerList)   aTypesList;
 998   Handle(TDataStd_BooleanList)   aClosuresList;
 999   getSectionsLists( aNamesList, aTypesList, aClosuresList );
1000
1001   TCollection_ExtendedString aSectName( theSectName );
1002   if ( aSectName.Length() <= 0 )
1003     aSectName = getUniqueSectionName( aNamesList );
1004
1005   aNamesList->Append( aSectName );
1006   aTypesList->Append( theSectionType );
1007   aClosuresList->Append( theIsClosed );
1008
1009   Changed( Geom_2d );
1010 }
1011
1012 TCollection_AsciiString HYDROData_PolylineXY::GetSectionName( const int theSectionIndex ) const
1013 {
1014   TCollection_AsciiString aResName;
1015
1016   Handle(TDataStd_ExtStringList) aNamesList;
1017   Handle(TDataStd_IntegerList) aTypesList;
1018   Handle(TDataStd_BooleanList) aClosuresList;
1019   getSectionsLists( aNamesList, aTypesList, aClosuresList, false );
1020   if ( aNamesList.IsNull() || theSectionIndex >= aNamesList->Extent() )
1021     return aResName;
1022
1023   TDataStd_ListIteratorOfListOfExtendedString aNamesIter( aNamesList->List() );
1024   for ( int i = 0; aNamesIter.More() && i != theSectionIndex; aNamesIter.Next(), ++i );
1025
1026   if ( aNamesIter.More() )
1027     aResName = aNamesIter.Value();
1028
1029   return aResName;
1030 }
1031
1032 void HYDROData_PolylineXY::SetSectionName( const int                      theSectionIndex,
1033                                            const TCollection_AsciiString& theSectionName )
1034 {
1035   Handle(TDataStd_ExtStringList) aNamesList;
1036   Handle(TDataStd_IntegerList) aTypesList;
1037   Handle(TDataStd_BooleanList) aClosuresList;
1038   getSectionsLists( aNamesList, aTypesList, aClosuresList, false );
1039   if ( aNamesList.IsNull() || theSectionIndex >= aNamesList->Extent() )
1040     return;
1041
1042   TDataStd_ListOfExtendedString anOldNamesList;
1043   anOldNamesList = aNamesList->List();
1044
1045   // Refill the existing list
1046   aNamesList->Clear();
1047
1048   TCollection_ExtendedString aNewSectName = theSectionName;
1049
1050   TDataStd_ListIteratorOfListOfExtendedString aNamesIter( anOldNamesList );
1051   for ( int i = 0; aNamesIter.More(); aNamesIter.Next(), ++i )
1052     aNamesList->Append( i == theSectionIndex ? aNewSectName : aNamesIter.Value() );
1053 }
1054
1055 HYDROData_PolylineXY::SectionType HYDROData_PolylineXY::GetSectionType( const int theSectionIndex ) const
1056 {
1057   Handle(TDataStd_ExtStringList) aNamesList;
1058   Handle(TDataStd_IntegerList) aTypesList;
1059   Handle(TDataStd_BooleanList) aClosuresList;
1060   getSectionsLists( aNamesList, aTypesList, aClosuresList, false );
1061   if ( aTypesList.IsNull() || theSectionIndex >= aTypesList->Extent() )
1062     return SECTION_POLYLINE;
1063
1064   TColStd_ListIteratorOfListOfInteger aTypesIter( aTypesList->List() );
1065   for ( int i = 0; aTypesIter.More() && i != theSectionIndex; aTypesIter.Next(), ++i );
1066
1067   return aTypesIter.More() ? (SectionType)aTypesIter.Value() : SECTION_POLYLINE;
1068 }
1069
1070 void HYDROData_PolylineXY::SetSectionType( const int         theSectionIndex,
1071                                            const SectionType theSectionType )
1072 {
1073   Handle(TDataStd_ExtStringList) aNamesList;
1074   Handle(TDataStd_IntegerList) aTypesList;
1075   Handle(TDataStd_BooleanList) aClosuresList;
1076   getSectionsLists( aNamesList, aTypesList, aClosuresList, false );
1077   if ( aTypesList.IsNull() || theSectionIndex >= aTypesList->Extent() )
1078     return;
1079
1080   TColStd_ListOfInteger anOldTypesList;
1081   anOldTypesList = aTypesList->List();
1082
1083   // Refill the existing list
1084   aTypesList->Clear();
1085
1086   TColStd_ListIteratorOfListOfInteger aTypesIter( anOldTypesList );
1087   for ( int i = 0; aTypesIter.More(); aTypesIter.Next(), ++i )
1088     aTypesList->Append( i == theSectionIndex ? theSectionType : aTypesIter.Value() );
1089
1090   Changed( Geom_2d );
1091 }
1092
1093 bool HYDROData_PolylineXY::IsClosedSection( const int theSectionIndex ) const
1094 {
1095   Handle(TDataStd_ExtStringList) aNamesList;
1096   Handle(TDataStd_IntegerList) aTypesList;
1097   Handle(TDataStd_BooleanList) aClosuresList;
1098   getSectionsLists( aNamesList, aTypesList, aClosuresList, false );
1099   if ( aClosuresList.IsNull() || theSectionIndex >= aClosuresList->Extent() )
1100     return false;
1101
1102   TDataStd_ListIteratorOfListOfByte aClosuresIter( aClosuresList->List() );
1103   for ( int i = 0; aClosuresIter.More() && i != theSectionIndex; aClosuresIter.Next(), ++i );
1104
1105   return aClosuresIter.More() ? (bool)aClosuresIter.Value() : false;
1106 }
1107
1108 void HYDROData_PolylineXY::SetSectionClosed( const int  theSectionIndex,
1109                                              const bool theIsClosed )
1110 {
1111   Handle(TDataStd_ExtStringList) aNamesList;
1112   Handle(TDataStd_IntegerList) aTypesList;
1113   Handle(TDataStd_BooleanList) aClosuresList;
1114   getSectionsLists( aNamesList, aTypesList, aClosuresList, false );
1115   if ( aClosuresList.IsNull() || theSectionIndex >= aClosuresList->Extent() )
1116     return;
1117
1118   TDataStd_ListOfByte anOldClosuresList;
1119   anOldClosuresList = aClosuresList->List();
1120
1121   // Refill the existing list
1122   aClosuresList->Clear();
1123
1124   TDataStd_ListIteratorOfListOfByte aClosuresIter( anOldClosuresList );
1125   for ( int i = 0; aClosuresIter.More(); aClosuresIter.Next(), ++i )
1126     aClosuresList->Append( i == theSectionIndex ? theIsClosed : (bool)aClosuresIter.Value() );
1127
1128   Changed( Geom_2d );
1129 }
1130
1131 void HYDROData_PolylineXY::GetSections( NCollection_Sequence<TCollection_AsciiString>& theSectNames,
1132                                         NCollection_Sequence<SectionType>&             theSectTypes,
1133                                         NCollection_Sequence<bool>&                    theSectClosures ) const
1134 {
1135   theSectNames.Clear();
1136   theSectTypes.Clear();
1137   theSectClosures.Clear();
1138
1139   if( IsCustom() )
1140   {
1141     const_cast<HYDROData_PolylineXY*>( this )->Interpolate();
1142   }
1143
1144   Handle(TDataStd_ExtStringList) aNamesList;
1145   Handle(TDataStd_IntegerList) aTypesList;
1146   Handle(TDataStd_BooleanList) aClosuresList;
1147   getSectionsLists( aNamesList, aTypesList, aClosuresList, false );
1148   if ( aNamesList.IsNull() || aTypesList.IsNull() || aClosuresList.IsNull() )
1149     return;
1150
1151   TDataStd_ListIteratorOfListOfExtendedString aNamesIter( aNamesList->List() );
1152   TColStd_ListIteratorOfListOfInteger aTypesIter( aTypesList->List() );
1153   TDataStd_ListIteratorOfListOfByte aClosuresIter( aClosuresList->List() );
1154   for ( ; aNamesIter.More() && aTypesIter.More() && aClosuresIter.More();
1155           aNamesIter.Next(), aTypesIter.Next(), aClosuresIter.Next() )
1156   {
1157     const TCollection_ExtendedString& aSectName = aNamesIter.Value();
1158     SectionType aSectType = (SectionType)aTypesIter.Value();
1159     bool aSectClosures = aClosuresIter.Value();
1160
1161     theSectNames.Append( aSectName );
1162     theSectTypes.Append( aSectType );
1163     theSectClosures.Append( aSectClosures );
1164   }
1165 }
1166
1167 void HYDROData_PolylineXY::RemoveSection( const int theSectionIndex )
1168 {
1169   Handle(TDataStd_ExtStringList) aNamesList;
1170   Handle(TDataStd_IntegerList)   aTypesList;
1171   Handle(TDataStd_BooleanList)   aClosuresList;
1172   getSectionsLists( aNamesList, aTypesList, aClosuresList, false );
1173   if ( aNamesList.IsNull() || theSectionIndex >= aNamesList->Extent() )
1174     return;
1175
1176   if ( aNamesList->Extent() == 1 )
1177   {
1178     removeSectionsLists();
1179     removePointsLists();
1180   }
1181   else
1182   {
1183     TDataStd_ListOfExtendedString anOldNamesList;
1184     anOldNamesList = aNamesList->List();
1185
1186     TColStd_ListOfInteger anOldTypesList;
1187     anOldTypesList = aTypesList->List();
1188
1189     TDataStd_ListOfByte anOldClosuresList;
1190     anOldClosuresList = aClosuresList->List();
1191
1192     // Refill the existing lists
1193     aNamesList->Clear();
1194     aTypesList->Clear();
1195     aClosuresList->Clear();
1196
1197     TDataStd_ListIteratorOfListOfExtendedString aNamesIter( anOldNamesList );
1198     TColStd_ListIteratorOfListOfInteger aTypesIter( anOldTypesList );
1199     TDataStd_ListIteratorOfListOfByte aClosuresIter( anOldClosuresList );
1200     for ( int i = 0; aNamesIter.More() && aTypesIter.More() && aClosuresIter.More();
1201                      aNamesIter.Next(), aTypesIter.Next(), aClosuresIter.Next(), ++i )
1202     {
1203       if ( i == theSectionIndex )
1204         continue; // skip index to remove
1205
1206       aNamesList->Append( aNamesIter.Value() );
1207       aTypesList->Append( aTypesIter.Value() );
1208       aClosuresList->Append( (bool)aClosuresIter.Value() );
1209     }
1210
1211     // Remove points that belongs to removed section
1212     removePointsLists( theSectionIndex );
1213   }
1214
1215   Changed( Geom_2d );
1216 }
1217
1218 void HYDROData_PolylineXY::RemoveSections()
1219 {
1220   removeSectionsLists();
1221   removePointsLists();
1222   Changed( Geom_2d );
1223 }
1224
1225 void HYDROData_PolylineXY::AddPoint( const int    theSectionIndex,
1226                                      const Point& thePoint,
1227                                      const int    thePointIndex )
1228 {
1229   Handle(TDataStd_RealList) aListX, aListY;
1230   getPointsLists( theSectionIndex, aListX, aListY );
1231
1232   if ( thePointIndex < 0 || thePointIndex >= aListX->Extent() )
1233   {
1234     aListX->Append( thePoint.X() );
1235     aListY->Append( thePoint.Y() );
1236   }
1237   else
1238   {
1239     TColStd_ListOfReal anOldListX;
1240     anOldListX = aListX->List();
1241
1242     TColStd_ListOfReal anOldListY;
1243     anOldListY = aListY->List();
1244
1245     // Refill the existing lists
1246     aListX->Clear();
1247     aListY->Clear();
1248
1249     TColStd_ListIteratorOfListOfReal anIterX( anOldListX );
1250     TColStd_ListIteratorOfListOfReal anIterY( anOldListY );
1251     for ( int i = 0; anIterX.More() && anIterY.More(); anIterX.Next(), anIterY.Next(), ++i )
1252     {
1253       double aCoordX = anIterX.Value();
1254       double aCoordY = anIterY.Value();
1255
1256       if ( i == thePointIndex )
1257       {
1258         // Insert our new point
1259         aListX->Append( thePoint.X() );
1260         aListY->Append( thePoint.Y() );
1261       }
1262
1263       aListX->Append( aCoordX );
1264       aListY->Append( aCoordY );
1265     }
1266   }
1267
1268   Changed( Geom_2d );
1269 }
1270
1271 void HYDROData_PolylineXY::SetPoint( const int    theSectionIndex,
1272                                      const Point& thePoint,
1273                                      const int    thePointIndex )
1274 {
1275   Handle(TDataStd_RealList) aListX, aListY;
1276   getPointsLists( theSectionIndex, aListX, aListY );
1277
1278   if ( thePointIndex < 0 )
1279   {
1280     aListX->Prepend( thePoint.X() );
1281     aListY->Prepend( thePoint.Y() );
1282   }
1283   else if ( thePointIndex >= aListX->Extent() )
1284   {
1285     aListX->Append( thePoint.X() );
1286     aListY->Append( thePoint.Y() );
1287   }
1288   else
1289   {
1290     TColStd_ListOfReal anOldListX;
1291     anOldListX = aListX->List();
1292
1293     TColStd_ListOfReal anOldListY;
1294     anOldListY = aListY->List();
1295
1296     // Refill the existing lists
1297     aListX->Clear();
1298     aListY->Clear();
1299
1300     TColStd_ListIteratorOfListOfReal anIterX( anOldListX );
1301     TColStd_ListIteratorOfListOfReal anIterY( anOldListY );
1302     for ( int i = 0; anIterX.More() && anIterY.More(); anIterX.Next(), anIterY.Next(), ++i )
1303     {
1304       double aCoordX = anIterX.Value();
1305       double aCoordY = anIterY.Value();
1306
1307       if ( i == thePointIndex )
1308       {
1309         // Insert our new point instead of old one
1310         aCoordX = thePoint.X();
1311         aCoordY = thePoint.Y();
1312       }
1313
1314       aListX->Append( aCoordX );
1315       aListY->Append( aCoordY );
1316     }
1317   }
1318
1319   Changed( Geom_2d );
1320 }
1321
1322 void HYDROData_PolylineXY::SetPoints( const int         theSectionIndex,
1323                                       const PointsList& thePoints )
1324 {
1325   Handle(TDataStd_RealList) aListX, aListY;
1326   getPointsLists( theSectionIndex, aListX, aListY );
1327
1328   aListX->Clear();
1329   aListY->Clear();
1330
1331   for ( int i = 1, n = thePoints.Length(); i <= n; ++i )
1332   {
1333     const Point& aPoint = thePoints.Value( i );
1334     aListX->Append( aPoint.X() );
1335     aListY->Append( aPoint.Y() );
1336   }
1337 }
1338
1339 void HYDROData_PolylineXY::RemovePoint( const int theSectionIndex,
1340                                         const int thePointIndex )
1341 {
1342   Handle(TDataStd_RealList) aListX, aListY;
1343   getPointsLists( theSectionIndex, aListX, aListY, false );
1344   if ( aListX.IsNull() || aListY.IsNull() || aListX->IsEmpty() )
1345     return;
1346
1347   if ( aListX->Extent() == 1 )
1348   {
1349     removePointsLists( theSectionIndex );
1350   }
1351   else
1352   {
1353     TColStd_ListOfReal anOldListX;
1354     anOldListX = aListX->List();
1355
1356     TColStd_ListOfReal anOldListY;
1357     anOldListY = aListY->List();
1358
1359     // Refill the existing lists
1360     aListX->Clear();
1361     aListY->Clear();
1362
1363     TColStd_ListIteratorOfListOfReal anIterX( anOldListX );
1364     TColStd_ListIteratorOfListOfReal anIterY( anOldListY );
1365     for ( int i = 0; anIterX.More() && anIterY.More(); anIterX.Next(), anIterY.Next(), ++i )
1366     {
1367       if ( i == thePointIndex )
1368         continue; // skip index to remove
1369
1370       aListX->Append( anIterX.Value() );
1371       aListY->Append( anIterY.Value() );
1372     }
1373   }
1374
1375   Changed( Geom_2d );
1376 }
1377
1378 HYDROData_PolylineXY::PointsList HYDROData_PolylineXY::GetPoints( const int theSectionIndex, bool IsConvertToGlobal ) const
1379 {
1380   PointsList aResList;
1381
1382   if( IsCustom() )
1383   {
1384     const_cast<HYDROData_PolylineXY*>( this )->Interpolate();
1385   }
1386
1387   Handle(TDataStd_RealList) aListX, aListY;
1388   getPointsLists( theSectionIndex, aListX, aListY, false );
1389   if ( aListX.IsNull() || aListY.IsNull() || aListX->IsEmpty() )
1390     return aResList;
1391
1392   TColStd_ListIteratorOfListOfReal anIterX( aListX->List() );
1393   TColStd_ListIteratorOfListOfReal anIterY( aListY->List() );
1394   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document( myLab );
1395   for ( ; anIterX.More() && anIterY.More(); anIterX.Next(), anIterY.Next() )
1396   {
1397     Point aPoint( anIterX.Value(), anIterY.Value() );
1398     if (IsConvertToGlobal)
1399       aDoc->Transform( aPoint, false );
1400     aResList.Append( aPoint );
1401   }
1402
1403   return aResList;
1404 }
1405
1406 QPainterPath HYDROData_PolylineXY::GetPainterPath() const
1407 {
1408   QPainterPath aPath;
1409
1410   NCollection_Sequence<TCollection_AsciiString>           aSectNames;
1411   NCollection_Sequence<HYDROData_PolylineXY::SectionType> aSectTypes;
1412   NCollection_Sequence<bool>                              aSectClosures;
1413   GetSections( aSectNames, aSectTypes, aSectClosures );
1414
1415   for ( int aSectionId = 1, aNbSects = aSectNames.Size(); aSectionId <= aNbSects; aSectionId++ )
1416   {
1417     TCollection_AsciiString aSectName = aSectNames.Value( aSectionId );
1418     SectionType aSectionType = aSectTypes.Value( aSectionId );
1419     bool anIsSectionClosed = aSectClosures.Value( aSectionId );
1420
1421     PointsList aSectPointsList = GetPoints( aSectionId - 1 );
1422     if ( aSectPointsList.IsEmpty() )
1423       continue;
1424
1425     NCollection_Sequence<gp_XYZ> aPoints;
1426     for( int i = 1, n = aSectPointsList.Size(); i <= n; ++i )
1427     {
1428       const Point& aSectPoint = aSectPointsList.Value( i );
1429
1430       gp_XYZ aPoint( aSectPoint.X(), aSectPoint.Y(), 0.0 );
1431       aPoints.Append( aPoint );
1432     }
1433
1434     BuildPainterPath( aPath, aSectionType, anIsSectionClosed, aPoints );
1435   }
1436
1437   return aPath;
1438 }
1439
1440 void HYDROData_PolylineXY::UpdateLocalCS( double theDx, double theDy )
1441 {
1442   NCollection_Sequence<TCollection_AsciiString>           aSectNames;
1443   NCollection_Sequence<HYDROData_PolylineXY::SectionType> aSectTypes;
1444   NCollection_Sequence<bool>                              aSectClosures;
1445   GetSections( aSectNames, aSectTypes, aSectClosures );
1446
1447   gp_XY aDelta( theDx, theDy );
1448   for ( int i = 0, aNbSects = aSectNames.Size(); i < aNbSects; i++ )
1449   {
1450     PointsList aPoints = GetPoints( i );
1451     for( int j = 1, n = aPoints.Size(); j <= n; ++j )
1452     {
1453       Point& aPoint = aPoints.ChangeValue( j );
1454       aPoint += aDelta;
1455     }
1456     SetPoints( i, aPoints );
1457   }
1458   Changed( Geom_2d );
1459 }
1460
1461 void HYDROData_PolylineXY::Transform( const QTransform& theTrsf )
1462 {
1463   NCollection_Sequence<TCollection_AsciiString>           aSectNames;
1464   NCollection_Sequence<HYDROData_PolylineXY::SectionType> aSectTypes;
1465   NCollection_Sequence<bool>                              aSectClosures;
1466   GetSections( aSectNames, aSectTypes, aSectClosures );
1467
1468   for ( int i = 0, aNbSects = aSectNames.Size(); i < aNbSects; i++ ) {
1469     PointsList aPoints = GetPoints( i );
1470     for( int j = 1, n = aPoints.Size(); j <= n; ++j ) {
1471       Point& aPoint = aPoints.ChangeValue( j );
1472
1473       QPointF aTrsfPoint = theTrsf.map( QPointF( aPoint.X(), aPoint.Y() ) );
1474
1475       aPoint.SetX( aTrsfPoint.x() );
1476       aPoint.SetY( aTrsfPoint.y() );
1477     }
1478     SetPoints( i, aPoints );
1479   }
1480
1481   Update();
1482 }
1483
1484 bool HYDROData_PolylineXY::IsCustom() const
1485 {
1486   if( myIsInCustomFlag )
1487     return false;
1488
1489   bool isNull = GetShape().IsNull();
1490   int aNbPoints = 0;
1491
1492   HYDROData_PolylineXY* aThat = const_cast<HYDROData_PolylineXY*>( this );
1493   aThat->myIsInCustomFlag = true;
1494   for( int i=0, n=NbSections(); i<n; i++ )
1495     aNbPoints += NbPoints( i );
1496   aThat->myIsInCustomFlag = false;
1497
1498   return !isNull && aNbPoints == 0;
1499 }
1500
1501 bool HYDROData_PolylineXY::GetIsInCustomFlag() const
1502 {
1503   return myIsInCustomFlag;
1504 }
1505
1506 void HYDROData_PolylineXY::SetIsInCustomFlag( bool theValue )
1507 {
1508   myIsInCustomFlag = theValue;
1509 }
1510
1511 void HYDROData_PolylineXY::Interpolate()
1512 {
1513   ImportShape( GetShape(), true, NULL );
1514 }