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