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