src/HYDROData/HYDROData_Profile.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_Profile.h"
  20
  21 #include "HYDROData_Document.h"
  22 #include "HYDROData_Iterator.h"
  23 #include "HYDROData_Tool.h"
  24 #include "HYDROData_PolylineXY.h"
  25
  26 #include <boost/math/special_functions/fpclassify.hpp>
  27
  28 #include <BRepBuilderAPI_MakeEdge.hxx>
  29 #include <BRepBuilderAPI_MakeWire.hxx>
  30 #include <BRepBuilderAPI_MakePolygon.hxx>
  31
  32 #include <BRepExtrema_ExtCC.hxx>
  33
  34 #include <BRep_Tool.hxx>
  35
  36 #include <gp_Lin.hxx>
  37 #include <gp_XY.hxx>
  38 #include <gp_XYZ.hxx>
  39 #include <gp_Pnt2d.hxx>
  40
  41 #include <TDataStd_AsciiString.hxx>
  42 #include <TDataStd_RealArray.hxx>
  43
  44 #include <TopoDS.hxx>
  45 #include <TopoDS_Edge.hxx>
  46 #include <TopoDS_Wire.hxx>
  47 #include <TopoDS_Iterator.hxx>
  48
  49 #include <OSD_File.hxx>
  50 #include <OSD_Protection.hxx>
  51
  52 #include <QColor>
  53 #include <QStringList>
  54
  55 IMPLEMENT_STANDARD_HANDLE(HYDROData_Profile, HYDROData_Object)
  56 IMPLEMENT_STANDARD_RTTIEXT(HYDROData_Profile, HYDROData_Object)
  57
  58 HYDROData_Profile::HYDROData_Profile()
  59 : HYDROData_Object()
  60 {
  61 }
  62
  63 HYDROData_Profile::~HYDROData_Profile()
  64 {
  65 }
  66
  67 QStringList HYDROData_Profile::DumpToPython( MapOfTreatedObjects& theTreatedObjects ) const
  68 {
  69   QStringList aResList = dumpObjectCreation( theTreatedObjects );
  70   QString aProfileName = GetObjPyName();
  71
  72   //TCollection_AsciiString aFilePath = GetFilePath();
  73   //if ( !aFilePath.IsEmpty() )
  74   //{
  75   //  aResList << QString( "%1.ImportFromFile( \"%2\" );" )
  76   //            .arg( aName ).arg( aFilePath.ToCString() );
  77   //}
  78
  79   bool anIsValidProfile = IsValid();
  80
  81   QStringList aPntsDefinition;
  82   QString aPntsListName = HYDROData_Tool::GenerateNameForPython( theTreatedObjects, "profile_points" );
  83
  84   QString aGap = QString().fill( ' ', aPntsListName.length() + 5 );
  85   if ( anIsValidProfile )
  86   {
  87     HYDROData_Profile::ProfilePoints aPointsList = GetProfilePoints( true );
  88     for ( int k = 1, aNbPoints = aPointsList.Size(); k <= aNbPoints; ++k )
  89     {
  90       const ProfilePoint& aPoint = aPointsList.Value( k );
  91       aPntsDefinition << QString( aGap + "gp_XYZ( %1, %2, %3 )%4" )
  92                          .arg( aPoint.X() ).arg( aPoint.Y() ).arg( aPoint.Z() )
  93                          .arg( ( k < aNbPoints ? "," : "" ) );
  94     }
  95   }
  96   else
  97   {
  98     HYDROData_IPolyline::PointsList aPointsList = GetParametricPoints();
  99     for ( int k = 1, aNbPoints = aPointsList.Size(); k <= aNbPoints; ++k )
 100     {
 101       const HYDROData_IPolyline::Point& aPoint = aPointsList.Value( k );
 102       aPntsDefinition << QString( aGap + "gp_XY( %1, %2 )%3" )
 103                          .arg( aPoint.X() ).arg( aPoint.Y() )
 104                          .arg( ( k < aNbPoints ? "," : "" ) );
 105     }
 106   }
 107
 108   if ( !aPntsDefinition.isEmpty() )
 109   {
 110     QString& aFirstStr = aPntsDefinition.first();
 111     aFirstStr = aFirstStr.trimmed();
 112     aFirstStr.prepend( QString( "%1 = [ " ).arg( aPntsListName ) );
 113
 114     aPntsDefinition.last().append( " ];" );
 115
 116     aResList << aPntsDefinition;
 117
 118     aResList << QString( "%1.%3( %2 );" )
 119                 .arg( aProfileName ).arg( aPntsListName )
 120                 .arg( anIsValidProfile ? "SetProfilePoints" : "SetParametricPoints" );
 121
 122     aResList << QString( "" );
 123   }
 124
 125   // Set a polyline type if it is not default
 126   Handle(HYDROData_ProfileUZ) aPrf = GetProfileUZ( false );
 127   if ( !aPrf.IsNull() )
 128   {
 129     HYDROData_IPolyline::SectionType aSecType = aPrf->GetSectionType( 0 );
 130     if ( aSecType != HYDROData_IPolyline::SECTION_POLYLINE )
 131     {
 132       aResList << QString( "%1.GetProfileUZ().SetSectionType( 0, %2 );" )
 133                   .arg( aProfileName ).arg( "HYDROData_IPolyline.SECTION_SPLINE" );
 134       aResList << QString( "" );
 135     }
 136   }
 137
 138   aResList << QString( "%1.Update();" ).arg( aProfileName );
 139   aResList << QString( "" );
 140
 141   return aResList;
 142 }
 143
 144 TopoDS_Shape HYDROData_Profile::GetTopShape() const
 145 {
 146   TopoDS_Wire aWire;
 147
 148   gp_XY aFirstPoint, aLastPoint;
 149   if ( !GetLeftPoint( aFirstPoint, false ) || !GetRightPoint( aLastPoint, false ) )
 150     return aWire;
 151
 152   gp_Pnt aPnt1( aFirstPoint.X(), aFirstPoint.Y(), 0 );
 153   gp_Pnt aPnt2( aLastPoint.X(),  aLastPoint.Y(),  0 );
 154
 155   BRepBuilderAPI_MakeEdge aMakeEdge( aPnt1, aPnt2 );
 156   TopoDS_Edge anEdge = aMakeEdge;
 157
 158   BRepBuilderAPI_MakeWire aMakeWire( anEdge );
 159   aWire = aMakeWire;
 160
 161   return aWire;
 162 }
 163
 164 TopoDS_Shape HYDROData_Profile::GetShape3D() const
 165 {
 166   TopoDS_Shape aShape = getShape3D();
 167   if( aShape.IsNull() )
 168     aShape = CreateProfileWire( true );
 169   return aShape;
 170 }
 171
 172 TopoDS_Shape HYDROData_Profile::CreateProfileWire( bool canUseDefaultPoints ) const
 173 {
 174   TopoDS_Wire aWire;
 175   Handle(HYDROData_ProfileUZ) aProfile = GetProfileUZ( false );
 176   if ( !aProfile.IsNull() )
 177   {
 178     ProfilePoints aProfilePoints = GetProfilePoints( false, canUseDefaultPoints );
 179     HYDROData_IPolyline::SectionType aSectionType = aProfile->GetSectionType( 0 );
 180
 181     aWire = HYDROData_PolylineXY::BuildWire( aSectionType, false, aProfilePoints );
 182   }
 183   return aWire;
 184 }
 185
 186 void HYDROData_Profile::Update()
 187 {
 188   HYDROData_Object::Update();
 189
 190   TopoDS_Shape aShape = CreateProfileWire( false );
 191   SetShape3D( aShape );
 192 }
 193
 194 QColor HYDROData_Profile::DefaultFillingColor()
 195 {
 196   return QColor( Qt::transparent );
 197 }
 198
 199 QColor HYDROData_Profile::DefaultBorderColor()
 200 {
 201   return QColor( Qt::black );
 202 }
 203
 204 QColor HYDROData_Profile::getDefaultFillingColor() const
 205 {
 206   return DefaultFillingColor();
 207 }
 208
 209 QColor HYDROData_Profile::getDefaultBorderColor() const
 210 {
 211   return DefaultBorderColor();
 212 }
 213
 214 bool HYDROData_Profile::IsValid() const
 215 {
 216   gp_XY aFirstPoint, aLastPoint;
 217   if ( !GetLeftPoint( aFirstPoint, false ) || !GetRightPoint( aLastPoint, false ) )
 218     return false;
 219
 220   int aNbPoints = NbPoints();
 221   return aNbPoints > 1;
 222 }
 223
 224 void HYDROData_Profile::SetLeftPoint( const gp_XY& theGPoint, bool IsConvertFromGlobal )
 225 {
 226   TDF_Label aLabel = myLab.FindChild( DataTag_FirstPoint );
 227   if ( aLabel.IsNull() )
 228     return;
 229
 230   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document( Label() );
 231   gp_XY aLPoint = theGPoint;
 232   if( IsConvertFromGlobal )
 233     aDoc->Transform( aLPoint, true );
 234
 235   Handle(TDataStd_RealArray) anArray;
 236   if ( !aLabel.FindAttribute( TDataStd_RealArray::GetID(), anArray ) )
 237     anArray = TDataStd_RealArray::Set( aLabel, 0, 1 );
 238
 239   anArray->SetValue( 0, aLPoint.X() );
 240   anArray->SetValue( 1, aLPoint.Y() );
 241
 242   SetToUpdate( true );
 243 }
 244
 245 bool HYDROData_Profile::GetLeftPoint( gp_XY& thePoint, bool IsConvertToGlobal,
 246                                       bool CanUseDefault ) const
 247 {
 248   HYDROData_ProfileUZ::PointsList aParametricPoints = GetParametricPoints();
 249   if ( aParametricPoints.Length() < 2 )
 250     return false;
 251
 252   thePoint = GetParametricPoints().First();
 253
 254   thePoint.SetX( 0 );
 255   thePoint.SetY( 0 ); //default left point of not-georeferenced profile
 256   TDF_Label aLabel = myLab.FindChild( DataTag_FirstPoint, false );
 257   if ( aLabel.IsNull() )
 258   {
 259     return CanUseDefault;
 260   }
 261
 262   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document( myLab );
 263   Handle(TDataStd_RealArray) anArray;
 264   if ( !aLabel.FindAttribute( TDataStd_RealArray::GetID(), anArray ) )
 265   {
 266     return CanUseDefault;
 267   }
 268
 269   thePoint.SetX( anArray->Value( 0 ) );
 270   thePoint.SetY( anArray->Value( 1 ) );
 271
 272   if( IsConvertToGlobal )
 273     aDoc->Transform( thePoint, false );
 274
 275   return true;
 276 }
 277
 278 void HYDROData_Profile::SetRightPoint( const gp_XY& theGPoint, bool IsConvertFromGlobal )
 279 {
 280   TDF_Label aLabel = myLab.FindChild( DataTag_LastPoint );
 281
 282   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document( Label() );
 283   gp_XY aLPoint = theGPoint;
 284   if( IsConvertFromGlobal )
 285     aDoc->Transform( aLPoint, true );
 286
 287   Handle(TDataStd_RealArray) anArray;
 288   if ( !aLabel.FindAttribute( TDataStd_RealArray::GetID(), anArray ) )
 289     anArray = TDataStd_RealArray::Set( aLabel, 0, 1 );
 290
 291   anArray->SetValue( 0, aLPoint.X() );
 292   anArray->SetValue( 1, aLPoint.Y() );
 293
 294   SetToUpdate( true );
 295 }
 296
 297 bool HYDROData_Profile::GetRightPoint( gp_XY& thePoint, bool IsConvertToGlobal,
 298                                        bool CanUseDefault ) const
 299 {
 300   HYDROData_ProfileUZ::PointsList aParametricPoints = GetParametricPoints();
 301   if ( aParametricPoints.Length() < 2 )
 302     return false;
 303
 304   thePoint = GetParametricPoints().Last();
 305
 306   TDF_Label aLabel = myLab.FindChild( DataTag_LastPoint, false );
 307   if ( aLabel.IsNull() )
 308   {
 309     return CanUseDefault;
 310   }
 311
 312   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document( myLab );
 313   Handle(TDataStd_RealArray) anArray;
 314   if ( !aLabel.FindAttribute( TDataStd_RealArray::GetID(), anArray ) )
 315   {
 316     return CanUseDefault;
 317   }
 318
 319   thePoint.SetX( anArray->Value( 0 ) );
 320   thePoint.SetY( anArray->Value( 1 ) );
 321
 322   if( IsConvertToGlobal )
 323     aDoc->Transform( thePoint, false );
 324
 325   return true;
 326 }
 327
 328 void HYDROData_Profile::Invalidate()
 329 {
 330   TDF_Label aFirstLabel = myLab.FindChild( DataTag_FirstPoint, false );
 331   if ( !aFirstLabel.IsNull() )
 332     aFirstLabel.ForgetAllAttributes();
 333
 334   TDF_Label aLastLabel = myLab.FindChild( DataTag_LastPoint, false );
 335   if ( !aLastLabel.IsNull() )
 336     aLastLabel.ForgetAllAttributes();
 337
 338   SetToUpdate( true );
 339 }
 340
 341 Handle(HYDROData_ProfileUZ) HYDROData_Profile::GetProfileUZ( const bool theIsCreate ) const
 342 {
 343   Handle(HYDROData_ProfileUZ) aProfileUZ;
 344
 345   TDF_Label aLabel = myLab.FindChild( DataTag_ChildProfileUZ, theIsCreate );
 346   if ( aLabel.IsNull() )
 347     return aProfileUZ;
 348
 349   aProfileUZ = Handle(HYDROData_ProfileUZ)::DownCast( HYDROData_Iterator::Object( aLabel ) );
 350   if ( aProfileUZ.IsNull() && theIsCreate )
 351   {
 352     aProfileUZ = Handle(HYDROData_ProfileUZ)::DownCast(
 353       HYDROData_Iterator::CreateObject( aLabel, KIND_PROFILEUZ ) );
 354   }
 355
 356   return aProfileUZ;
 357 }
 358
 359 int HYDROData_Profile::NbPoints() const
 360 {
 361   Handle(HYDROData_ProfileUZ) aProfileUZ = GetProfileUZ( false );
 362   return aProfileUZ.IsNull() ? 0 : aProfileUZ->NbPoints();
 363 }
 364
 365 void HYDROData_Profile::RemovePoints()
 366 {
 367   Handle(HYDROData_ProfileUZ) aProfileUZ = GetProfileUZ( false );
 368   if ( !aProfileUZ.IsNull() )
 369   {
 370     aProfileUZ->RemoveSections();
 371     SetToUpdate( true );
 372   }
 373 }
 374
 375 void HYDROData_Profile::SetParametricPoints( const HYDROData_ProfileUZ::PointsList& thePoints )
 376 {
 377   RemovePoints();
 378
 379   Handle(HYDROData_ProfileUZ) aProfileUZ = GetProfileUZ();
 380   for ( int i = 1, n = thePoints.Length(); i <= n ; ++i )
 381   {
 382     const HYDROData_ProfileUZ::Point& aPoint = thePoints.Value( i );
 383     aProfileUZ->AddPoint( 0, aPoint );
 384   }
 385
 386   SetToUpdate( true );
 387 }
 388
 389 HYDROData_ProfileUZ::PointsList HYDROData_Profile::GetParametricPoints() const
 390 {
 391   Handle(HYDROData_ProfileUZ) aProfileUZ = GetProfileUZ( false );
 392   return aProfileUZ.IsNull() ? HYDROData_ProfileUZ::PointsList() : aProfileUZ->GetPoints();
 393 }
 394
 395 void HYDROData_Profile::SetProfilePoints( const ProfilePoints& thePoints, bool IsConvertFromGlobal )
 396 {
 397   RemovePoints();
 398   if ( thePoints.Length() < 2 )
 399     return;
 400
 401   gp_XY aFirstPoint, aLastPoint;
 402
 403   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document( Label() );
 404   Handle(HYDROData_ProfileUZ) aProfileUZ = GetProfileUZ();
 405   for ( int i = 1, n = thePoints.Length(); i <= n ; ++i )
 406   {
 407     ProfilePoint aPoint = thePoints.Value( i );
 408     if( IsConvertFromGlobal )
 409       aDoc->Transform( aPoint, true );
 410
 411     gp_XY aPointXY( aPoint.X(), aPoint.Y() );
 412
 413     if ( i == 1 )
 414       aFirstPoint = aPointXY;
 415     else if ( i == n )
 416       aLastPoint = aPointXY;
 417
 418     double aDistance = gp_Pnt2d( aFirstPoint ).Distance( aPointXY );
 419
 420     HYDROData_ProfileUZ::Point aParPoint( aDistance, aPoint.Z() );
 421     aProfileUZ->AddPoint( 0, aParPoint );
 422   }
 423
 424   SetLeftPoint( aFirstPoint, false );//already converted to local CS
 425   SetRightPoint( aLastPoint, false );
 426 }
 427
 428 HYDROData_Profile::ProfilePoints HYDROData_Profile::GetProfilePoints
 429   ( bool IsConvertToGlobal, bool CanUseDefaultLeftRight ) const
 430 {
 431   ProfilePoints aResPoints;
 432
 433   gp_XY aFirstPoint, aLastPoint;
 434   if ( !GetLeftPoint( aFirstPoint, IsConvertToGlobal, CanUseDefaultLeftRight ) ||
 435        !GetRightPoint( aLastPoint, IsConvertToGlobal, CanUseDefaultLeftRight ) )
 436     return aResPoints;
 437
 438   HYDROData_ProfileUZ::PointsList aParametricPoints = GetParametricPoints();
 439   if ( aParametricPoints.Length() < 2 )
 440     return aResPoints;
 441
 442   const HYDROData_ProfileUZ::Point& aFirstParPoint = aParametricPoints.First();
 443   const HYDROData_ProfileUZ::Point& aLastParPoint = aParametricPoints.Last();
 444
 445   double aGeoDistance = gp_Pnt2d( aFirstPoint ).Distance( aLastPoint );
 446   double aParCommonDist = gp_Pnt2d( aFirstParPoint.X(), 0 ).Distance( gp_Pnt2d( aLastParPoint.X(), 0 ) );
 447
 448   // Add first point as is
 449   aResPoints.Append( ProfilePoint( aFirstPoint.X(), aFirstPoint.Y(), aFirstParPoint.Y() ) );
 450
 451   // Compute all other points
 452   for ( int i = 2, n = aParametricPoints.Length(); i < n ; ++i )
 453   {
 454     const HYDROData_ProfileUZ::Point& aParPoint = aParametricPoints.Value( i );
 455
 456     double aParPointDist = gp_Pnt2d( aFirstParPoint.X(), 0 ).Distance( gp_Pnt2d( aParPoint.X(), 0 ) );
 457
 458     double aParLen = ( aParPointDist / aParCommonDist ) * aGeoDistance;
 459
 460     double aRatio = aParLen / ( aGeoDistance - aParLen );
 461
 462     double aParX = ( aFirstPoint.X() + aRatio * aLastPoint.X() ) / ( 1 + aRatio );
 463     double aParY = ( aFirstPoint.Y() + aRatio * aLastPoint.Y() ) / ( 1 + aRatio );
 464
 465     ProfilePoint aCompPoint( aParX, aParY, aParPoint.Y() );
 466     aResPoints.Append( aCompPoint );
 467   }
 468
 469   // Add last point as is
 470   aResPoints.Append( ProfilePoint( aLastPoint.X(), aLastPoint.Y(), aLastParPoint.Y() ) );
 471
 472   return aResPoints;
 473 }
 474
 475 void HYDROData_Profile::SetFilePath( const TCollection_AsciiString& theFilePath )
 476 {
 477   TDataStd_AsciiString::Set( myLab.FindChild( DataTag_FilePath ), theFilePath );
 478 }
 479
 480 TCollection_AsciiString HYDROData_Profile::GetFilePath() const
 481 {
 482   TCollection_AsciiString aRes;
 483
 484   Handle(TDataStd_AsciiString) anAsciiStr;
 485   if ( myLab.FindChild( DataTag_FilePath ).FindAttribute( TDataStd_AsciiString::GetID(), anAsciiStr ) )
 486     aRes = anAsciiStr->Get();
 487
 488   return aRes;
 489 }
 490
 491 int HYDROData_Profile::ImportFromFile( const Handle(HYDROData_Document)& theDoc,
 492                                        const TCollection_AsciiString&    theFileName,
 493                                        NCollection_Sequence<int>&        theBadProfilesIds )
 494 {
 495   if ( theDoc.IsNull() || theFileName.IsEmpty() )
 496     return 0;
 497
 498   OSD_File aFile( theFileName );
 499   if ( !aFile.IsReadable() )
 500     return 0;
 501
 502   aFile.Open( OSD_ReadOnly, OSD_Protection() );
 503   if ( !aFile.IsOpen() )
 504     return 0;
 505
 506   NCollection_Sequence<Handle(HYDROData_Profile)> aCreatedProfiles;
 507
 508   int aProfileId = 1;
 509   Handle(HYDROData_Profile) aNewProfile;
 510   for ( ; !aFile.IsAtEnd(); ++aProfileId )
 511   {
 512     if ( aNewProfile.IsNull() )
 513       aNewProfile = Handle(HYDROData_Profile)::DownCast( theDoc->CreateObject( KIND_PROFILE ) );
 514
 515     bool anIsRead = false;
 516     if ( aNewProfile->ImportFromFile( aFile, &anIsRead ) )
 517     {
 518       aCreatedProfiles.Append( aNewProfile );
 519       aNewProfile.Nullify();
 520     }
 521     else if ( anIsRead )
 522     {
 523       theBadProfilesIds.Append( aProfileId );
 524     }
 525   }
 526
 527   if ( !aNewProfile.IsNull() )
 528     aNewProfile->Remove();
 529
 530   // Close the file
 531   aFile.Close();
 532
 533   for ( int i = 1, n = aCreatedProfiles.Length(); i <= n ; ++i )
 534   {
 535     Handle(HYDROData_Profile) aProfile = aCreatedProfiles.Value( i );
 536
 537     QString aProfileName = HYDROData_Tool::GenerateObjectName( theDoc, "Profile" );
 538     aProfile->SetName( aProfileName );
 539
 540     aProfile->SetFilePath( theFileName );
 541
 542     aProfile->SetBorderColor( HYDROData_Profile::DefaultBorderColor() );
 543   }
 544
 545   return aCreatedProfiles.Length();
 546 }
 547
 548 bool HYDROData_Profile::ImportFromFile( const TCollection_AsciiString& theFileName,
 549                                         bool*                          theIsRead )
 550 {
 551   if ( theIsRead )
 552     *theIsRead = false;
 553
 554   // Try to open the file
 555   OSD_File aFile( theFileName );
 556   if ( !aFile.IsReadable() )
 557     return false;
 558
 559   aFile.Open( OSD_ReadOnly, OSD_Protection() );
 560   if ( !aFile.IsOpen() )
 561     return false;
 562
 563   bool aRes = ImportFromFile( aFile, theIsRead );
 564
 565   // Close the file
 566   aFile.Close();
 567
 568   if ( aRes )
 569   {
 570     // Update file path
 571     SetFilePath( theFileName );
 572   }
 573
 574   return aRes;
 575 }
 576
 577 bool HYDROData_Profile::ImportFromFile( OSD_File& theFile,
 578                                         bool*     theIsRead )
 579 {
 580   if ( theIsRead )
 581     *theIsRead = false;
 582
 583   if ( !theFile.IsOpen() )
 584     return false;
 585
 586   bool aRes = true;
 587
 588   bool anIsParametric = false;
 589   bool anIsGeoref     = false;
 590
 591   HYDROData_ProfileUZ::PointsList aPointsUZ;
 592   ProfilePoints                   aPointsXYZ;
 593
 594   double aPrevVal = -DBL_MAX;
 595   while ( !theFile.IsAtEnd() )
 596   {
 597     Standard_Integer aNbRead = 0;
 598     TCollection_AsciiString aLine;
 599     theFile.ReadLine( aLine, 1024, aNbRead );
 600
 601     aLine.LeftAdjust(); aLine.RightAdjust();
 602     if ( aLine.IsEmpty() )
 603     {
 604       if ( !anIsParametric && !anIsGeoref )
 605         continue; // Definition is not started yet
 606
 607       break; // Next profile started
 608     }
 609
 610     // Set flag of read status to true
 611     if ( theIsRead )
 612       *theIsRead = true;
 613
 614     TCollection_AsciiString aValX = aLine.Token( " \t", 1 );
 615     TCollection_AsciiString aValY = aLine.Token( " \t", 2 );
 616     TCollection_AsciiString aValZ = aLine.Token( " \t", 3 );
 617
 618     if ( aValX.IsEmpty() || !aValX.IsRealValue() ||
 619          aValY.IsEmpty() || !aValY.IsRealValue() )
 620     {
 621       aRes = false;
 622       break;
 623     }
 624
 625     if ( !anIsParametric && !anIsGeoref )
 626     {
 627       anIsParametric = aValZ.IsEmpty();
 628       anIsGeoref = !aValZ.IsEmpty();
 629     }
 630
 631     double aCoordX = aValX.RealValue();
 632     double aCoordY = aValY.RealValue();
 633
 634     if ( boost::math::isnan( aCoordX ) || boost::math::isinf( aCoordX ) ||
 635          boost::math::isnan( aCoordY ) || boost::math::isinf( aCoordY ) )
 636       aRes = false;
 637
 638     if ( anIsParametric )
 639     {
 640       if ( aCoordX < aPrevVal )
 641       {
 642         // Move back readed line
 643         theFile.Seek( -( aNbRead + 1 ), OSD_FromHere );
 644         break;
 645       }
 646
 647       HYDROData_ProfileUZ::Point aPoint( aCoordX, aCoordY );
 648       aPointsUZ.Append( aPoint );
 649
 650       aPrevVal = aCoordX;
 651     }
 652     else
 653     {
 654       if ( aValZ.IsEmpty() || !aValZ.IsRealValue() )
 655       {
 656         aRes = false;
 657         break;
 658       }
 659
 660       double aCoordZ = aValZ.RealValue();
 661       if ( boost::math::isnan( aCoordZ ) || boost::math::isinf( aCoordZ ) )
 662         aRes = false;
 663
 664       ProfilePoint aPoint( aCoordX, aCoordY, aCoordZ );
 665       aPointsXYZ.Append( aPoint );
 666     }
 667   }
 668
 669   aRes = aRes && ( anIsParametric && !aPointsUZ.IsEmpty() ||
 670                    anIsGeoref && !aPointsXYZ.IsEmpty() );
 671   if ( aRes )
 672   {
 673     // Update profile points
 674     if ( anIsParametric )
 675     {
 676       SetParametricPoints( aPointsUZ );
 677     }
 678     else if ( anIsGeoref )
 679     {
 680       SetProfilePoints( aPointsXYZ, true );
 681     }
 682
 683     Update();
 684   }
 685
 686   return aRes;
 687 }
 688
 689 void HYDROData_Profile::UpdateLocalCS( double theDx, double theDy )
 690 {
 691   gp_XY aDelta( theDx, theDy );
 692   gp_XY aPnt;
 693
 694   GetLeftPoint( aPnt, false );
 695   aPnt += aDelta;
 696   SetLeftPoint( aPnt, false );
 697
 698   GetRightPoint( aPnt, false );
 699   aPnt += aDelta;
 700   SetRightPoint( aPnt, false );
 701 }
 702
 703 HYDROData_Profile::ProfilePoint HYDROData_Profile::GetBottomPoint() const
 704 {
 705   ProfilePoint aBottom;
 706
 707   // Get parametric points
 708   HYDROData_ProfileUZ::PointsList aParametricPoints = GetParametricPoints();
 709   if ( aParametricPoints.Length() < 1 ) {
 710     return aBottom;
 711   }
 712
 713   // Calculate midvalue for U parameter
 714   Standard_Real anUMidValue = aParametricPoints.First().X();
 715   Standard_Real anUMinValue = anUMidValue;
 716   Standard_Real anUMaxValue = anUMidValue;
 717
 718   for ( int i = 2, aNbPoints = aParametricPoints.Size(); i <= aNbPoints; i++ ) {
 719     const HYDROData_IPolyline::Point& aParPoint = aParametricPoints.Value( i );
 720     Standard_Real anU = aParPoint.X();
 721
 722     if ( anU < anUMinValue ) {
 723       anUMinValue = anU;
 724     } else if ( anU > anUMaxValue ) {
 725       anUMaxValue = anU;
 726     }
 727   }
 728
 729   anUMidValue = ( anUMinValue + anUMaxValue ) / 2;
 730
 731   // Find index of the parametric point with minimal Z value
 732   int aBottomIndex = 1;
 733   HYDROData_IPolyline::Point aParBottom = aParametricPoints.First();
 734
 735   for ( int i = 2, aNbPoints = aParametricPoints.Size(); i <= aNbPoints; i++ ) {
 736     const HYDROData_IPolyline::Point& aParPoint = aParametricPoints.Value( i );
 737     if ( aParPoint.Y() < aParBottom.Y() ) {
 738       aBottomIndex = i;
 739       aParBottom = aParPoint;
 740     } else if ( aParPoint.Y() == aParBottom.Y() ) {
 741       // Check which point is neares to the U = 0.5
 742       if ( fabs( aParPoint.X() - anUMidValue ) < fabs( aParBottom.X() - anUMidValue ) ) {
 743         aBottomIndex = i;
 744         aParBottom = aParPoint;
 745       }
 746     }
 747   }
 748
 749   // Find the corresponding profile point
 750   ProfilePoints aProfilePoints = GetProfilePoints( false );
 751   if ( aBottomIndex >= 1 && aBottomIndex <= aProfilePoints.Length() ) {
 752     aBottom = aProfilePoints.Value( aBottomIndex );
 753   }
 754
 755   return aBottom;
 756 }
 757
 758  HYDROData_Profile::ProfilePoint HYDROData_Profile::GetMiddlePoint() const
 759  {
 760    ProfilePoint aMiddlePoint;
 761
 762    gp_XY aLeftPnt, aRightPnt;
 763    if ( GetLeftPoint( aLeftPnt, false ) && GetRightPoint( aRightPnt, false ) ) {
 764      gp_XYZ aPnt1( aLeftPnt.X(), aLeftPnt.Y(), 0. );
 765      gp_XYZ aPnt2( aRightPnt.X(), aRightPnt.Y(), 0. );
 766      gp_Pnt aMiddlePoint2d( 0.5 * ( aPnt1 + aPnt2 ) );
 767
 768      gp_Lin aMidLin( aMiddlePoint2d, gp::DZ() );
 769      TopoDS_Edge aMidEdge = BRepLib_MakeEdge( aMidLin );
 770
 771      TopoDS_Iterator anIt( TopoDS::Wire( GetShape3D() )  );
 772      for ( ; anIt.More(); anIt.Next()) {
 773        const TopoDS_Edge& anEdge = TopoDS::Edge( anIt.Value() );
 774
 775        /*
 776        Standard_Real aStart, anEnd;
 777        Handle(Geom_Curve) aCurve = BRep_Tool::Curve( anEdge, aStart, anEnd );
 778        gp_Pnt aMiddlePointOnCurve = aCurve->Value( ( aStart + anEnd ) / 2 );
 779        */
 780
 781        BRepExtrema_ExtCC ExtremaEE( aMidEdge, anEdge);
 782        if (ExtremaEE.IsDone() && ExtremaEE.NbExt() != 0) {
 783          for ( Standard_Integer i = 1; i <= ExtremaEE.NbExt(); i++ ) {
 784            if ( ExtremaEE.SquareDistance(i) <= Precision::Confusion() ) {
 785              aMiddlePoint = ExtremaEE.PointOnE1(i).XYZ();
 786              break;
 787            }
 788          }
 789        }
 790      }
 791    }
 792
 793    return aMiddlePoint;
 794  }