]> SALOME platform Git repositories - modules/hydro.git/blob - src/HYDROData/HYDROData_Profile.cxx
Salome HOME
mode normalization
[modules/hydro.git] / 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( Geom_3d )
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 = HYDROData_Object::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( true );
191   SetShape3D( aShape );
192 }
193
194 QColor HYDROData_Profile::DefaultFillingColor() const
195 {
196   return QColor( Qt::transparent );
197 }
198
199 QColor HYDROData_Profile::DefaultBorderColor() const
200 {
201   return QColor( Qt::black );
202 }
203
204 bool HYDROData_Profile::IsValid() const
205 {
206   gp_XY aFirstPoint, aLastPoint;
207   if ( !GetLeftPoint( aFirstPoint, false ) || !GetRightPoint( aLastPoint, false ) )
208     return false;
209
210   int aNbPoints = NbPoints();
211   return aNbPoints > 1;
212 }
213
214 void HYDROData_Profile::SetLeftPoint( const gp_XY& theGPoint, bool IsConvertFromGlobal )
215 {
216   TDF_Label aLabel = myLab.FindChild( DataTag_FirstPoint );
217   if ( aLabel.IsNull() )
218     return;
219
220   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document( Label() );
221   gp_XY aLPoint = theGPoint;
222   if( IsConvertFromGlobal )
223     aDoc->Transform( aLPoint, true );
224
225   Handle(TDataStd_RealArray) anArray;
226   if ( !aLabel.FindAttribute( TDataStd_RealArray::GetID(), anArray ) )
227     anArray = TDataStd_RealArray::Set( aLabel, 0, 1 );
228
229   anArray->SetValue( 0, aLPoint.X() );
230   anArray->SetValue( 1, aLPoint.Y() );
231
232   Changed( Geom_3d );
233 }
234
235 bool HYDROData_Profile::GetLeftPoint( gp_XY& thePoint, bool IsConvertToGlobal,
236                                       bool CanUseDefault ) const
237 {
238   HYDROData_ProfileUZ::PointsList aParametricPoints = GetParametricPoints();
239   if ( aParametricPoints.Length() < 2 )
240     return false;
241
242   thePoint = GetParametricPoints().First();
243
244   //thePoint.SetX( 0 );
245   thePoint.SetY( 0 ); //default left point of not-georeferenced profile
246   TDF_Label aLabel = myLab.FindChild( DataTag_FirstPoint, false );
247   if ( aLabel.IsNull() )
248   {
249     return CanUseDefault;
250   }
251
252   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document( myLab );
253   Handle(TDataStd_RealArray) anArray;
254   if ( !aLabel.FindAttribute( TDataStd_RealArray::GetID(), anArray ) )
255   {
256     return CanUseDefault;
257   }
258
259   thePoint.SetX( anArray->Value( 0 ) );
260   thePoint.SetY( anArray->Value( 1 ) );
261
262   if( IsConvertToGlobal )
263     aDoc->Transform( thePoint, false );
264
265   return true;
266 }
267
268 void HYDROData_Profile::SetRightPoint( const gp_XY& theGPoint, bool IsConvertFromGlobal )
269 {
270   TDF_Label aLabel = myLab.FindChild( DataTag_LastPoint );
271
272   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document( Label() );
273   gp_XY aLPoint = theGPoint;
274   if( IsConvertFromGlobal )
275     aDoc->Transform( aLPoint, true );
276
277   Handle(TDataStd_RealArray) anArray;
278   if ( !aLabel.FindAttribute( TDataStd_RealArray::GetID(), anArray ) )
279     anArray = TDataStd_RealArray::Set( aLabel, 0, 1 );
280
281   anArray->SetValue( 0, aLPoint.X() );
282   anArray->SetValue( 1, aLPoint.Y() );
283
284   Changed( Geom_3d );
285 }
286
287 bool HYDROData_Profile::GetRightPoint( gp_XY& thePoint, bool IsConvertToGlobal,
288                                        bool CanUseDefault ) const
289 {
290   HYDROData_ProfileUZ::PointsList aParametricPoints = GetParametricPoints();
291   if ( aParametricPoints.Length() < 2 )
292     return false;
293
294   thePoint = GetParametricPoints().Last();
295   thePoint.SetY( 0 );
296
297   TDF_Label aLabel = myLab.FindChild( DataTag_LastPoint, false );
298   if ( aLabel.IsNull() )
299   {
300     return CanUseDefault;
301   }
302
303   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document( myLab );
304   Handle(TDataStd_RealArray) anArray;
305   if ( !aLabel.FindAttribute( TDataStd_RealArray::GetID(), anArray ) )
306   {
307     return CanUseDefault;
308   }
309
310   thePoint.SetX( anArray->Value( 0 ) );
311   thePoint.SetY( anArray->Value( 1 ) );
312
313   if( IsConvertToGlobal )
314     aDoc->Transform( thePoint, false );
315
316   return true;
317 }
318
319 void HYDROData_Profile::Invalidate()
320 {
321   TDF_Label aFirstLabel = myLab.FindChild( DataTag_FirstPoint, false );
322   if ( !aFirstLabel.IsNull() )
323     aFirstLabel.ForgetAllAttributes();
324
325   TDF_Label aLastLabel = myLab.FindChild( DataTag_LastPoint, false );
326   if ( !aLastLabel.IsNull() )
327     aLastLabel.ForgetAllAttributes();
328
329   Changed( Geom_3d );
330 }
331
332 Handle(HYDROData_ProfileUZ) HYDROData_Profile::GetProfileUZ( const bool theIsCreate ) const
333 {
334   Handle(HYDROData_ProfileUZ) aProfileUZ;
335
336   TDF_Label aLabel = myLab.FindChild( DataTag_ChildProfileUZ, theIsCreate );
337   if ( aLabel.IsNull() )
338     return aProfileUZ;
339
340   aProfileUZ = Handle(HYDROData_ProfileUZ)::DownCast( HYDROData_Iterator::Object( aLabel ) );
341   if ( aProfileUZ.IsNull() && theIsCreate )
342   {
343     aProfileUZ = Handle(HYDROData_ProfileUZ)::DownCast(
344       HYDROData_Iterator::CreateObject( aLabel, KIND_PROFILEUZ ) );
345   }
346
347   return aProfileUZ;
348 }
349
350 int HYDROData_Profile::NbPoints() const
351 {
352   Handle(HYDROData_ProfileUZ) aProfileUZ = GetProfileUZ( false );
353   return aProfileUZ.IsNull() ? 0 : aProfileUZ->NbPoints();
354 }
355
356 void HYDROData_Profile::RemovePoints()
357 {
358   Handle(HYDROData_ProfileUZ) aProfileUZ = GetProfileUZ( false );
359   if ( !aProfileUZ.IsNull() )
360   {
361     aProfileUZ->RemoveSections();
362     Changed( Geom_3d );
363   }
364 }
365
366 void HYDROData_Profile::SetParametricPoints( const HYDROData_ProfileUZ::PointsList& thePoints )
367 {
368   RemovePoints();
369
370   Handle(HYDROData_ProfileUZ) aProfileUZ = GetProfileUZ();
371   for ( int i = 1, n = thePoints.Length(); i <= n ; ++i )
372   {
373     const HYDROData_ProfileUZ::Point& aPoint = thePoints.Value( i );
374     aProfileUZ->AddPoint( 0, aPoint );
375   }
376
377   Changed( Geom_3d );
378 }
379
380 HYDROData_ProfileUZ::PointsList HYDROData_Profile::GetParametricPoints() const
381 {
382   Handle(HYDROData_ProfileUZ) aProfileUZ = GetProfileUZ( false );
383   return aProfileUZ.IsNull() ? HYDROData_ProfileUZ::PointsList() : aProfileUZ->GetPoints();
384 }
385
386 void HYDROData_Profile::SetProfilePoints( const ProfilePoints& thePoints, bool IsConvertFromGlobal )
387 {
388   RemovePoints();
389   if ( thePoints.Length() < 2 )
390     return;
391
392   gp_XY aFirstPoint, aLastPoint;
393
394   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document( Label() );
395   Handle(HYDROData_ProfileUZ) aProfileUZ = GetProfileUZ();
396   for ( int i = 1, n = thePoints.Length(); i <= n ; ++i )
397   {
398     ProfilePoint aPoint = thePoints.Value( i );
399     if( IsConvertFromGlobal )
400       aDoc->Transform( aPoint, true );
401
402     gp_XY aPointXY( aPoint.X(), aPoint.Y() );
403
404     if ( i == 1 )
405       aFirstPoint = aPointXY;
406     else if ( i == n )
407       aLastPoint = aPointXY;
408
409     double aDistance = gp_Pnt2d( aFirstPoint ).Distance( aPointXY );
410     
411     HYDROData_ProfileUZ::Point aParPoint( aDistance, aPoint.Z() );
412     aProfileUZ->AddPoint( 0, aParPoint );
413   }
414
415   SetLeftPoint( aFirstPoint, false );//already converted to local CS
416   SetRightPoint( aLastPoint, false );
417 }
418
419 HYDROData_Profile::ProfilePoints HYDROData_Profile::GetProfilePoints
420   ( bool IsConvertToGlobal, bool CanUseDefaultLeftRight ) const
421 {
422   ProfilePoints aResPoints;
423
424   gp_XY aFirstPoint, aLastPoint;
425   if ( !GetLeftPoint( aFirstPoint, IsConvertToGlobal, CanUseDefaultLeftRight ) ||
426        !GetRightPoint( aLastPoint, IsConvertToGlobal, CanUseDefaultLeftRight ) )
427     return aResPoints;
428
429   HYDROData_ProfileUZ::PointsList aParametricPoints = GetParametricPoints();
430   if ( aParametricPoints.Length() < 2 )
431     return aResPoints;
432
433   const HYDROData_ProfileUZ::Point& aFirstParPoint = aParametricPoints.First();
434   const HYDROData_ProfileUZ::Point& aLastParPoint = aParametricPoints.Last();
435
436   double aFullLength = aLastPoint.X() - aFirstPoint.X();
437   double aParFullLength = aLastParPoint.X() - aFirstParPoint.X();
438
439   // Add first point as is
440   aResPoints.Append( ProfilePoint( aFirstPoint.X(), aFirstPoint.Y(), aFirstParPoint.Y() ) );
441
442   // Compute all other points
443   for ( int i = 2, n = aParametricPoints.Length(); i < n ; ++i )
444   {
445     const HYDROData_ProfileUZ::Point& aParPoint = aParametricPoints.Value( i );
446
447     double aParPointDist = aParPoint.X() - aFirstParPoint.X();
448     double aRatio = aParPointDist / aParFullLength;
449
450     double aParX = aFirstPoint.X() * (1-aRatio) + aLastPoint.X() * aRatio;
451     double aParY = aFirstPoint.Y() * (1-aRatio) + aLastPoint.Y() * aRatio;
452
453     ProfilePoint aCompPoint( aParX, aParY, aParPoint.Y() );
454     aResPoints.Append( aCompPoint );
455   }
456
457   // Add last point as is
458   aResPoints.Append( ProfilePoint( aLastPoint.X(), aLastPoint.Y(), aLastParPoint.Y() ) );
459
460   return aResPoints;
461 }
462
463 void HYDROData_Profile::SetFilePath( const TCollection_AsciiString& theFilePath )
464 {
465   TDataStd_AsciiString::Set( myLab.FindChild( DataTag_FilePath ), theFilePath );
466 }
467
468 TCollection_AsciiString HYDROData_Profile::GetFilePath() const
469 {
470   TCollection_AsciiString aRes;
471
472   Handle(TDataStd_AsciiString) anAsciiStr;
473   if ( myLab.FindChild( DataTag_FilePath ).FindAttribute( TDataStd_AsciiString::GetID(), anAsciiStr ) )
474     aRes = anAsciiStr->Get();
475
476   return aRes;
477 }
478
479 int HYDROData_Profile::ImportFromFile( const Handle(HYDROData_Document)& theDoc,
480                                        const TCollection_AsciiString&    theFileName,
481                                        NCollection_Sequence<int>&        theBadProfilesIds )
482 {
483   if ( theDoc.IsNull() || theFileName.IsEmpty() )
484     return 0;
485
486   OSD_File aFile( theFileName );
487   if ( !aFile.IsReadable() )
488     return 0;
489
490   aFile.Open( OSD_ReadOnly, OSD_Protection() );
491   if ( !aFile.IsOpen() )
492     return 0;
493
494   NCollection_Sequence<Handle(HYDROData_Profile)> aCreatedProfiles;
495
496   int aProfileId = 1;
497   Handle(HYDROData_Profile) aNewProfile;
498   for ( ; !aFile.IsAtEnd(); ++aProfileId )
499   {
500     if ( aNewProfile.IsNull() )
501       aNewProfile = Handle(HYDROData_Profile)::DownCast( theDoc->CreateObject( KIND_PROFILE ) );
502     
503     bool anIsRead = false;
504     if ( aNewProfile->ImportFromFile( aFile, &anIsRead ) )
505     {
506       aCreatedProfiles.Append( aNewProfile );
507       aNewProfile.Nullify();
508     }
509     else if ( anIsRead )
510     {
511       theBadProfilesIds.Append( aProfileId );
512     }
513   }
514
515   if ( !aNewProfile.IsNull() )
516     aNewProfile->Remove();
517
518   // Close the file
519   aFile.Close();
520
521   for ( int i = 1, n = aCreatedProfiles.Length(); i <= n ; ++i )
522   {
523     Handle(HYDROData_Profile) aProfile = aCreatedProfiles.Value( i );
524
525     QString aProfileName = HYDROData_Tool::GenerateObjectName( theDoc, "Profile" );
526     aProfile->SetName( aProfileName );
527
528     aProfile->SetFilePath( theFileName );
529
530     aProfile->SetBorderColor( aProfile->DefaultBorderColor() );
531   }
532
533   return aCreatedProfiles.Length();
534 }
535
536 bool HYDROData_Profile::ImportFromFile( const TCollection_AsciiString& theFileName,
537                                         bool*                          theIsRead )
538 {
539   if ( theIsRead )
540     *theIsRead = false;
541
542   // Try to open the file
543   OSD_File aFile( theFileName );
544   if ( !aFile.IsReadable() )
545     return false;
546
547   aFile.Open( OSD_ReadOnly, OSD_Protection() );
548   if ( !aFile.IsOpen() )
549     return false;
550
551   bool aRes = ImportFromFile( aFile, theIsRead );
552
553   // Close the file
554   aFile.Close();
555
556   if ( aRes )
557   {
558     // Update file path
559     SetFilePath( theFileName );
560   }
561
562   return aRes;
563 }
564
565 bool HYDROData_Profile::ImportFromFile( OSD_File& theFile,
566                                         bool*     theIsRead )
567 {
568   if ( theIsRead )
569     *theIsRead = false;
570
571   if ( !theFile.IsOpen() )
572     return false;
573
574   bool aRes = true;
575
576   bool anIsParametric = false;
577   bool anIsGeoref     = false;
578
579   HYDROData_ProfileUZ::PointsList aPointsUZ;
580   ProfilePoints                   aPointsXYZ;
581
582   double aPrevVal = -DBL_MAX;
583   while ( !theFile.IsAtEnd() )
584   {
585     Standard_Integer aNbRead = 0;
586     TCollection_AsciiString aLine;
587     theFile.ReadLine( aLine, 1024, aNbRead );
588
589     aLine.LeftAdjust(); aLine.RightAdjust();
590     if ( aLine.IsEmpty() )
591     {
592       if ( !anIsParametric && !anIsGeoref )
593         continue; // Definition is not started yet
594
595       break; // Next profile started
596     }
597
598     // Set flag of read status to true
599     if ( theIsRead )
600       *theIsRead = true;
601
602     TCollection_AsciiString aValX = aLine.Token( " \t", 1 );
603     TCollection_AsciiString aValY = aLine.Token( " \t", 2 );
604     TCollection_AsciiString aValZ = aLine.Token( " \t", 3 );
605
606     if ( aValX.IsEmpty() || !aValX.IsRealValue() ||
607          aValY.IsEmpty() || !aValY.IsRealValue() )
608     {
609       aRes = false;
610       break;
611     }
612
613     if ( !anIsParametric && !anIsGeoref )
614     {
615       anIsParametric = aValZ.IsEmpty();
616       anIsGeoref = !aValZ.IsEmpty();
617     }
618
619     double aCoordX = aValX.RealValue();
620     double aCoordY = aValY.RealValue();
621
622     if ( boost::math::isnan( aCoordX ) || boost::math::isinf( aCoordX ) ||
623          boost::math::isnan( aCoordY ) || boost::math::isinf( aCoordY ) )
624       aRes = false;
625
626     if ( anIsParametric )
627     {
628       if ( aCoordX < aPrevVal )
629       {
630         // Move back readed line
631         theFile.Seek( -( aNbRead + 1 ), OSD_FromHere );
632         break;
633       }
634
635       HYDROData_ProfileUZ::Point aPoint( aCoordX, aCoordY );
636       aPointsUZ.Append( aPoint );
637
638       aPrevVal = aCoordX;
639     }
640     else
641     {
642       if ( aValZ.IsEmpty() || !aValZ.IsRealValue() )
643       {
644         aRes = false;
645         break;
646       }
647
648       double aCoordZ = aValZ.RealValue();
649       if ( boost::math::isnan( aCoordZ ) || boost::math::isinf( aCoordZ ) )
650         aRes = false;
651
652       ProfilePoint aPoint( aCoordX, aCoordY, aCoordZ );
653       aPointsXYZ.Append( aPoint );
654     }
655   }
656   
657   aRes = aRes && ( anIsParametric && !aPointsUZ.IsEmpty() || 
658                    anIsGeoref && !aPointsXYZ.IsEmpty() );
659   if ( aRes )
660   {
661     // Update profile points
662     if ( anIsParametric )
663     {
664       SetParametricPoints( aPointsUZ );
665     }
666     else if ( anIsGeoref )
667     {
668       SetProfilePoints( aPointsXYZ, true );
669     }
670
671     Update();
672   }
673
674   return aRes;
675 }
676
677 void HYDROData_Profile::UpdateLocalCS( double theDx, double theDy )
678 {
679   gp_XY aDelta( theDx, theDy );
680   gp_XY aPnt;
681
682   GetLeftPoint( aPnt, false );
683   aPnt += aDelta;
684   SetLeftPoint( aPnt, false );
685
686   GetRightPoint( aPnt, false );
687   aPnt += aDelta;
688   SetRightPoint( aPnt, false );
689 }
690
691 HYDROData_Profile::ProfilePoint HYDROData_Profile::GetBottomPoint() const
692 {
693   ProfilePoint aBottom;
694
695   // Get parametric points
696   HYDROData_ProfileUZ::PointsList aParametricPoints = GetParametricPoints();
697   if ( aParametricPoints.Length() < 1 ) {
698     return aBottom;
699   }
700
701   // Calculate midvalue for U parameter
702   Standard_Real anUMidValue = aParametricPoints.First().X();
703   Standard_Real anUMinValue = anUMidValue;
704   Standard_Real anUMaxValue = anUMidValue;
705
706   for ( int i = 2, aNbPoints = aParametricPoints.Size(); i <= aNbPoints; i++ ) {
707     const HYDROData_IPolyline::Point& aParPoint = aParametricPoints.Value( i );
708     Standard_Real anU = aParPoint.X();
709
710     if ( anU < anUMinValue ) {
711       anUMinValue = anU;
712     } else if ( anU > anUMaxValue ) {
713       anUMaxValue = anU;
714     }
715   }
716
717   anUMidValue = ( anUMinValue + anUMaxValue ) / 2;
718
719   // Find index of the parametric point with minimal Z value
720   int aBottomIndex = 1;
721   HYDROData_IPolyline::Point aParBottom = aParametricPoints.First();
722
723   for ( int i = 2, aNbPoints = aParametricPoints.Size(); i <= aNbPoints; i++ ) {
724     const HYDROData_IPolyline::Point& aParPoint = aParametricPoints.Value( i );
725     if ( aParPoint.Y() < aParBottom.Y() ) {
726       aBottomIndex = i;
727       aParBottom = aParPoint;
728     } else if ( aParPoint.Y() == aParBottom.Y() ) {
729       // Check which point is neares to the U = 0.5
730       if ( fabs( aParPoint.X() - anUMidValue ) < fabs( aParBottom.X() - anUMidValue ) ) {
731         aBottomIndex = i;
732         aParBottom = aParPoint;
733       }
734     }
735   }
736
737   // Find the corresponding profile point
738   ProfilePoints aProfilePoints = GetProfilePoints( false );
739   if ( aBottomIndex >= 1 && aBottomIndex <= aProfilePoints.Length() ) {
740     aBottom = aProfilePoints.Value( aBottomIndex );
741   }
742
743   return aBottom;
744 }
745
746  HYDROData_Profile::ProfilePoint HYDROData_Profile::GetMiddlePoint( bool CanUseDefault ) const
747  {
748    ProfilePoint aMiddlePoint;
749   
750    gp_XY aLeftPnt, aRightPnt;
751    if ( GetLeftPoint( aLeftPnt, true, CanUseDefault ) && GetRightPoint( aRightPnt, true, CanUseDefault ) ) {
752      gp_XYZ aPnt1( aLeftPnt.X(), aLeftPnt.Y(), 0. );
753      gp_XYZ aPnt2( aRightPnt.X(), aRightPnt.Y(), 0. );
754      gp_Pnt aMiddlePoint2d( 0.5 * ( aPnt1 + aPnt2 ) ); 
755
756      gp_Lin aMidLin( aMiddlePoint2d, gp::DZ() );
757      TopoDS_Edge aMidEdge = BRepLib_MakeEdge( aMidLin );
758
759      TopoDS_Iterator anIt( TopoDS::Wire( GetShape3D() )  );
760      for ( ; anIt.More(); anIt.Next()) {
761        const TopoDS_Edge& anEdge = TopoDS::Edge( anIt.Value() );
762
763        /*
764        Standard_Real aStart, anEnd;
765        Handle(Geom_Curve) aCurve = BRep_Tool::Curve( anEdge, aStart, anEnd );
766        gp_Pnt aMiddlePointOnCurve = aCurve->Value( ( aStart + anEnd ) / 2 );
767        */
768
769        BRepExtrema_ExtCC ExtremaEE( aMidEdge, anEdge);
770        if (ExtremaEE.IsDone() && ExtremaEE.NbExt() != 0) {
771          for ( Standard_Integer i = 1; i <= ExtremaEE.NbExt(); i++ ) {
772            if ( ExtremaEE.SquareDistance(i) <= Precision::Confusion() ) {
773              aMiddlePoint = ExtremaEE.PointOnE1(i).XYZ();
774              break;
775            }
776          }
777        }
778      }
779    }
780
781    return aMiddlePoint;
782  }