#include "HYDROData_Document.h"
#include "HYDROData_Iterator.h"
#include "HYDROData_Tool.h"
+#include "HYDROData_PolylineXY.h"
+
+#include <boost/math/special_functions/fpclassify.hpp>
#include <BRepBuilderAPI_MakeEdge.hxx>
#include <BRepBuilderAPI_MakeWire.hxx>
+#include <BRepBuilderAPI_MakePolygon.hxx>
#include <gp_XY.hxx>
#include <gp_XYZ.hxx>
#include <OSD_File.hxx>
#include <OSD_Protection.hxx>
+#include <QColor>
#include <QStringList>
-#define PYTHON_PROFILE_ID "KIND_PROFILE"
-
IMPLEMENT_STANDARD_HANDLE(HYDROData_Profile, HYDROData_Object)
IMPLEMENT_STANDARD_RTTIEXT(HYDROData_Profile, HYDROData_Object)
{
}
+QStringList HYDROData_Profile::DumpToPython( MapOfTreatedObjects& theTreatedObjects ) const
+{
+ QStringList aResList = dumpObjectCreation( theTreatedObjects );
+ QString aProfileName = GetObjPyName();
+
+ //TCollection_AsciiString aFilePath = GetFilePath();
+ //if ( !aFilePath.IsEmpty() )
+ //{
+ // aResList << QString( "%1.ImportFromFile( \"%2\" );" )
+ // .arg( aName ).arg( aFilePath.ToCString() );
+ //}
+
+ bool anIsValidProfile = IsValid();
+
+ QStringList aPntsDefinition;
+ QString aPntsListName = HYDROData_Tool::GenerateNameForPython( theTreatedObjects, "profile_points" );
+
+ QString aGap = QString().fill( ' ', aPntsListName.length() + 5 );
+ if ( anIsValidProfile )
+ {
+ HYDROData_Profile::ProfilePoints aPointsList = GetProfilePoints();
+ for ( int k = 1, aNbPoints = aPointsList.Size(); k <= aNbPoints; ++k )
+ {
+ const ProfilePoint& aPoint = aPointsList.Value( k );
+ aPntsDefinition << QString( aGap + "gp_XYZ( %1, %2, %3 )%4" )
+ .arg( aPoint.X() ).arg( aPoint.Y() ).arg( aPoint.Z() )
+ .arg( ( k < aNbPoints ? "," : "" ) );
+ }
+ }
+ else
+ {
+ HYDROData_IPolyline::PointsList aPointsList = GetParametricPoints();
+ for ( int k = 1, aNbPoints = aPointsList.Size(); k <= aNbPoints; ++k )
+ {
+ const HYDROData_IPolyline::Point& aPoint = aPointsList.Value( k );
+ aPntsDefinition << QString( aGap + "gp_XY( %1, %2 )%3" )
+ .arg( aPoint.X() ).arg( aPoint.Y() )
+ .arg( ( k < aNbPoints ? "," : "" ) );
+ }
+ }
+
+ if ( !aPntsDefinition.isEmpty() )
+ {
+ QString& aFirstStr = aPntsDefinition.first();
+ aFirstStr = aFirstStr.trimmed();
+ aFirstStr.prepend( QString( "%1 = [ " ).arg( aPntsListName ) );
+
+ aPntsDefinition.last().append( " ];" );
+
+ aResList << aPntsDefinition;
+
+ aResList << QString( "%1.%3( %2 );" )
+ .arg( aProfileName ).arg( aPntsListName )
+ .arg( anIsValidProfile ? "SetProfilePoints" : "SetParametricPoints" );
+
+ aResList << QString( "" );
+ }
+
+ // Set a polyline type if it is not default
+ Handle(HYDROData_ProfileUZ) aPrf = GetProfileUZ( false );
+ if ( !aPrf.IsNull() )
+ {
+ HYDROData_IPolyline::SectionType aSecType = aPrf->GetSectionType( 0 );
+ if ( aSecType != HYDROData_IPolyline::SECTION_POLYLINE )
+ {
+ aResList << QString( "%1.GetProfileUZ().SetSectionType( 0, %2 );" )
+ .arg( aProfileName ).arg( "HYDROData_IPolyline.SECTION_SPLINE" );
+ aResList << QString( "" );
+ }
+ }
+
+ aResList << QString( "%1.Update();" ).arg( aProfileName );
+ aResList << QString( "" );
+
+ return aResList;
+}
+
TopoDS_Shape HYDROData_Profile::GetTopShape() const
{
TopoDS_Wire aWire;
gp_XY aFirstPoint, aLastPoint;
- if ( !GetFirstPoint( aFirstPoint ) || !GetLastPoint( aLastPoint ) )
+ if ( !GetLeftPoint( aFirstPoint ) || !GetRightPoint( aLastPoint ) )
return aWire;
gp_Pnt aPnt1( aFirstPoint.X(), aFirstPoint.Y(), 0 );
return getShape3D();
}
-void HYDROData_Profile::UpdateShape3D()
+void HYDROData_Profile::Update()
{
- BRepBuilderAPI_MakeWire aMakeWire;
+ HYDROData_Object::Update();
- ProfilePoints aProfilePoints = GetProfilePoints();
- for ( int i = 1, n = aProfilePoints.Length(); i < n ; ++i )
+ TopoDS_Wire aWire;
+ Handle(HYDROData_ProfileUZ) aProfile = GetProfileUZ( false );
+ if ( !aProfile.IsNull() )
{
- ProfilePoint aFirstPoint = aProfilePoints.Value( i );
- ProfilePoint aSecPoint = aProfilePoints.Value( i + 1 );
-
- gp_Pnt aPnt1( aFirstPoint.X(), aFirstPoint.Y(), aFirstPoint.Z() );
- gp_Pnt aPnt2( aSecPoint.X(), aSecPoint.Y(), aSecPoint.Z() );
-
- BRepBuilderAPI_MakeEdge aMakeEdge( aPnt1, aPnt2 );
- TopoDS_Edge anEdge = aMakeEdge;
+ ProfilePoints aProfilePoints = GetProfilePoints();
+ HYDROData_IPolyline::SectionType aSectionType = aProfile->GetSectionType( 0 );
- aMakeWire.Add( anEdge );
+ aWire = HYDROData_PolylineXY::BuildWire( aSectionType, false, aProfilePoints );
}
-
- TopoDS_Wire aWire;
- if ( aMakeWire.IsDone() )
- aWire = aMakeWire;
-
SetShape3D( aWire );
}
-/**
- * Dump object to Python script representation.
- */
-QStringList HYDROData_Profile::DumpToPython( MapOfTreatedObjects& theTreatedObjects ) const
+QColor HYDROData_Profile::DefaultFillingColor()
{
- QStringList aResList;
+ return QColor( Qt::transparent );
+}
- Handle(HYDROData_Document) aDocument = HYDROData_Document::Document( myLab );
- if ( aDocument.IsNull() )
- return aResList;
-
- QString aDocName = aDocument->GetDocPyName();
- QString aProfileName = GetName();
+QColor HYDROData_Profile::DefaultBorderColor()
+{
+ return QColor( Qt::black );
+}
- aResList << QString( "%1 = %2.CreateObject( %3 );" )
- .arg( aProfileName ).arg( aDocName ).arg( PYTHON_PROFILE_ID );
- aResList << QString( "%1.SetName( \"%1\" );" ).arg( aProfileName );
+QColor HYDROData_Profile::getDefaultFillingColor() const
+{
+ return DefaultFillingColor();
+}
- return aResList;
+QColor HYDROData_Profile::getDefaultBorderColor() const
+{
+ return DefaultBorderColor();
}
bool HYDROData_Profile::IsValid() const
{
gp_XY aFirstPoint, aLastPoint;
- if ( !GetFirstPoint( aFirstPoint ) || !GetLastPoint( aLastPoint ) )
+ if ( !GetLeftPoint( aFirstPoint ) || !GetRightPoint( aLastPoint ) )
return false;
int aNbPoints = NbPoints();
return aNbPoints > 1;
}
-void HYDROData_Profile::SetFirstPoint( const gp_XY& thePoint )
+void HYDROData_Profile::SetLeftPoint( const gp_XY& thePoint )
{
TDF_Label aLabel = myLab.FindChild( DataTag_FirstPoint );
anArray->SetValue( 0, thePoint.X() );
anArray->SetValue( 1, thePoint.Y() );
+
+ SetToUpdate( true );
}
-bool HYDROData_Profile::GetFirstPoint( gp_XY& thePoint ) const
+bool HYDROData_Profile::GetLeftPoint( gp_XY& thePoint ) const
{
TDF_Label aLabel = myLab.FindChild( DataTag_FirstPoint, false );
if ( aLabel.IsNull() )
return true;
}
-void HYDROData_Profile::SetLastPoint( const gp_XY& thePoint )
+void HYDROData_Profile::SetRightPoint( const gp_XY& thePoint )
{
TDF_Label aLabel = myLab.FindChild( DataTag_LastPoint );
anArray->SetValue( 0, thePoint.X() );
anArray->SetValue( 1, thePoint.Y() );
+
+ SetToUpdate( true );
}
-bool HYDROData_Profile::GetLastPoint( gp_XY& thePoint ) const
+bool HYDROData_Profile::GetRightPoint( gp_XY& thePoint ) const
{
TDF_Label aLabel = myLab.FindChild( DataTag_LastPoint, false );
if ( aLabel.IsNull() )
return true;
}
+void HYDROData_Profile::Invalidate()
+{
+ TDF_Label aFirstLabel = myLab.FindChild( DataTag_FirstPoint, false );
+ if ( !aFirstLabel.IsNull() )
+ aFirstLabel.ForgetAllAttributes();
+
+ TDF_Label aLastLabel = myLab.FindChild( DataTag_LastPoint, false );
+ if ( !aLastLabel.IsNull() )
+ aLastLabel.ForgetAllAttributes();
+
+ SetToUpdate( true );
+}
+
Handle(HYDROData_ProfileUZ) HYDROData_Profile::GetProfileUZ( const bool theIsCreate ) const
{
Handle(HYDROData_ProfileUZ) aProfileUZ;
int HYDROData_Profile::NbPoints() const
{
Handle(HYDROData_ProfileUZ) aProfileUZ = GetProfileUZ( false );
- return aProfileUZ.IsNull() ? 0 : aProfileUZ->getNbPoints();
+ return aProfileUZ.IsNull() ? 0 : aProfileUZ->NbPoints();
}
void HYDROData_Profile::RemovePoints()
{
Handle(HYDROData_ProfileUZ) aProfileUZ = GetProfileUZ( false );
if ( !aProfileUZ.IsNull() )
- aProfileUZ->clear();
+ {
+ aProfileUZ->RemoveSections();
+ SetToUpdate( true );
+ }
}
-void HYDROData_Profile::SetParametricPoints( const CurveCreator::Coordinates& theCoords )
+void HYDROData_Profile::SetParametricPoints( const HYDROData_ProfileUZ::PointsList& thePoints )
{
RemovePoints();
Handle(HYDROData_ProfileUZ) aProfileUZ = GetProfileUZ();
-
- CurveCreator::Coordinates::const_iterator aBegIter = theCoords.begin();
- CurveCreator::Coordinates::const_iterator anEndIter = theCoords.end();
- while ( aBegIter != anEndIter )
+ for ( int i = 1, n = thePoints.Length(); i <= n ; ++i )
{
- const CurveCreator::TypeCoord& aCoordX = *aBegIter++;
- if ( aBegIter == anEndIter )
- break;
-
- const CurveCreator::TypeCoord& aCoordY = *aBegIter++;
-
- aProfileUZ->addPoint( aCoordX, aCoordY, 0 );
+ const HYDROData_ProfileUZ::Point& aPoint = thePoints.Value( i );
+ aProfileUZ->AddPoint( 0, aPoint );
}
+
+ SetToUpdate( true );
}
-CurveCreator::Coordinates HYDROData_Profile::GetParametricPoints() const
+HYDROData_ProfileUZ::PointsList HYDROData_Profile::GetParametricPoints() const
{
Handle(HYDROData_ProfileUZ) aProfileUZ = GetProfileUZ( false );
- return aProfileUZ.IsNull() ? CurveCreator::Coordinates() : aProfileUZ->getPoints();
+ return aProfileUZ.IsNull() ? HYDROData_ProfileUZ::PointsList() : aProfileUZ->GetPoints();
}
void HYDROData_Profile::SetProfilePoints( const ProfilePoints& thePoints )
double aDistance = gp_Pnt2d( aFirstPoint ).Distance( aPointXY );
- aProfileUZ->addPoint( 0, aDistance, aPoint.Z() );
+ HYDROData_ProfileUZ::Point aParPoint( aDistance, aPoint.Z() );
+ aProfileUZ->AddPoint( 0, aParPoint );
}
- SetFirstPoint( aFirstPoint );
- SetLastPoint( aLastPoint );
+ SetLeftPoint( aFirstPoint );
+ SetRightPoint( aLastPoint );
}
HYDROData_Profile::ProfilePoints HYDROData_Profile::GetProfilePoints() const
ProfilePoints aResPoints;
gp_XY aFirstPoint, aLastPoint;
- if ( !GetFirstPoint( aFirstPoint ) || !GetLastPoint( aLastPoint ) )
+ if ( !GetLeftPoint( aFirstPoint ) || !GetRightPoint( aLastPoint ) )
return aResPoints;
- CurveCreator::Coordinates aParametricPoints = GetParametricPoints();
- if ( aParametricPoints.size() < 4 )
+ HYDROData_ProfileUZ::PointsList aParametricPoints = GetParametricPoints();
+ if ( aParametricPoints.Length() < 2 )
return aResPoints;
- CurveCreator::Coordinates::const_iterator aBegIter = aParametricPoints.begin();
- CurveCreator::Coordinates::const_iterator anEndIter = aParametricPoints.end()--;
-
- gp_Pnt2d aFirstParPoint( *aBegIter++, *aBegIter++ );
- gp_Pnt2d aLastParPoint( *anEndIter--, *anEndIter );
+ const HYDROData_ProfileUZ::Point& aFirstParPoint = aParametricPoints.First();
+ const HYDROData_ProfileUZ::Point& aLastParPoint = aParametricPoints.Last();
double aGeoDistance = gp_Pnt2d( aFirstPoint ).Distance( aLastPoint );
double aParCommonDist = gp_Pnt2d( aFirstParPoint.X(), 0 ).Distance( gp_Pnt2d( aLastParPoint.X(), 0 ) );
aResPoints.Append( ProfilePoint( aFirstPoint.X(), aFirstPoint.Y(), aFirstParPoint.Y() ) );
// Compute all other points
- while ( aBegIter != anEndIter )
+ for ( int i = 2, n = aParametricPoints.Length(); i < n ; ++i )
{
- const CurveCreator::TypeCoord& aCoordU = *aBegIter++;
- if ( aBegIter == anEndIter )
- break;
+ const HYDROData_ProfileUZ::Point& aParPoint = aParametricPoints.Value( i );
- double aParPointDist = gp_Pnt2d( aFirstParPoint.X(), 0 ).Distance( gp_Pnt2d( aCoordU, 0 ) );
+ double aParPointDist = gp_Pnt2d( aFirstParPoint.X(), 0 ).Distance( gp_Pnt2d( aParPoint.X(), 0 ) );
double aParLen = ( aParPointDist / aParCommonDist ) * aGeoDistance;
double aParX = ( aFirstPoint.X() + aRatio * aLastPoint.X() ) / ( 1 + aRatio );
double aParY = ( aFirstPoint.Y() + aRatio * aLastPoint.Y() ) / ( 1 + aRatio );
-
- const CurveCreator::TypeCoord& aCoordZ = *aBegIter++;
-
- ProfilePoint aCompPoint( aParX, aParY, aCoordZ );
+ ProfilePoint aCompPoint( aParX, aParY, aParPoint.Y() );
aResPoints.Append( aCompPoint );
}
return aRes;
}
-bool HYDROData_Profile::ImportFromFile( const Handle(HYDROData_Document)& theDoc,
- const TCollection_AsciiString& theFileName )
+int HYDROData_Profile::ImportFromFile( const Handle(HYDROData_Document)& theDoc,
+ const TCollection_AsciiString& theFileName,
+ NCollection_Sequence<int>& theBadProfilesIds )
{
if ( theDoc.IsNull() || theFileName.IsEmpty() )
- return false;
+ return 0;
OSD_File aFile( theFileName );
if ( !aFile.IsReadable() )
- return false;
+ return 0;
aFile.Open( OSD_ReadOnly, OSD_Protection() );
if ( !aFile.IsOpen() )
- return false;
-
- bool aRes = true;
+ return 0;
NCollection_Sequence<Handle(HYDROData_Profile)> aCreatedProfiles;
- while ( aRes && !aFile.IsAtEnd() )
+ int aProfileId = 1;
+ Handle(HYDROData_Profile) aNewProfile;
+ for ( ; !aFile.IsAtEnd(); ++aProfileId )
{
- Handle(HYDROData_Profile) aNewProfile =
- Handle(HYDROData_Profile)::DownCast( theDoc->CreateObject( KIND_PROFILE ) );
+ if ( aNewProfile.IsNull() )
+ aNewProfile = Handle(HYDROData_Profile)::DownCast( theDoc->CreateObject( KIND_PROFILE ) );
- aRes = aNewProfile->ImportFromFile( aFile );
-
- aCreatedProfiles.Append( aNewProfile );
+ bool anIsRead = false;
+ if ( aNewProfile->ImportFromFile( aFile, &anIsRead ) )
+ {
+ aCreatedProfiles.Append( aNewProfile );
+ aNewProfile.Nullify();
+ }
+ else if ( anIsRead )
+ {
+ theBadProfilesIds.Append( aProfileId );
+ }
}
+ if ( !aNewProfile.IsNull() )
+ aNewProfile->Remove();
+
// Close the file
aFile.Close();
- for ( int i = 1, n = aCreatedProfiles.Length(); i < n ; ++i )
+ for ( int i = 1, n = aCreatedProfiles.Length(); i <= n ; ++i )
{
Handle(HYDROData_Profile) aProfile = aCreatedProfiles.Value( i );
- if ( aRes )
- {
- QString aProfileName = HYDROData_Tool::GenerateObjectName( theDoc, "Profile" );
- aProfile->SetName( aProfileName );
- aProfile->SetFilePath( theFileName );
- }
- else
- {
- aProfile->Remove();
- }
+ QString aProfileName = HYDROData_Tool::GenerateObjectName( theDoc, "Profile" );
+ aProfile->SetName( aProfileName );
+
+ aProfile->SetFilePath( theFileName );
+
+ aProfile->SetBorderColor( HYDROData_Profile::DefaultBorderColor() );
}
- return aRes;
+ return aCreatedProfiles.Length();
}
-bool HYDROData_Profile::ImportFromFile( const TCollection_AsciiString& theFileName )
+bool HYDROData_Profile::ImportFromFile( const TCollection_AsciiString& theFileName,
+ bool* theIsRead )
{
+ if ( theIsRead )
+ *theIsRead = false;
+
// Try to open the file
OSD_File aFile( theFileName );
if ( !aFile.IsReadable() )
if ( !aFile.IsOpen() )
return false;
- bool aRes = ImportFromFile( aFile );
+ bool aRes = ImportFromFile( aFile, theIsRead );
// Close the file
aFile.Close();
return aRes;
}
-bool HYDROData_Profile::ImportFromFile( OSD_File& theFile )
+bool HYDROData_Profile::ImportFromFile( OSD_File& theFile,
+ bool* theIsRead )
{
+ if ( theIsRead )
+ *theIsRead = false;
+
if ( !theFile.IsOpen() )
return false;
bool anIsParametric = false;
bool anIsGeoref = false;
- CurveCreator::Coordinates aPointsUZ;
- ProfilePoints aPointsXYZ;
+ HYDROData_ProfileUZ::PointsList aPointsUZ;
+ ProfilePoints aPointsXYZ;
double aPrevVal = -DBL_MAX;
while ( !theFile.IsAtEnd() )
aLine.LeftAdjust(); aLine.RightAdjust();
if ( aLine.IsEmpty() )
- continue;
+ {
+ if ( !anIsParametric && !anIsGeoref )
+ continue; // Definition is not started yet
+
+ break; // Next profile started
+ }
+
+ // Set flag of read status to true
+ if ( theIsRead )
+ *theIsRead = true;
TCollection_AsciiString aValX = aLine.Token( " \t", 1 );
TCollection_AsciiString aValY = aLine.Token( " \t", 2 );
double aCoordX = aValX.RealValue();
double aCoordY = aValY.RealValue();
+ if ( boost::math::isnan( aCoordX ) || boost::math::isinf( aCoordX ) ||
+ boost::math::isnan( aCoordY ) || boost::math::isinf( aCoordY ) )
+ aRes = false;
+
if ( anIsParametric )
{
if ( aCoordX < aPrevVal )
break;
}
- aPointsUZ.push_back( aCoordX );
- aPointsUZ.push_back( aCoordY );
+ HYDROData_ProfileUZ::Point aPoint( aCoordX, aCoordY );
+ aPointsUZ.Append( aPoint );
aPrevVal = aCoordX;
}
}
double aCoordZ = aValZ.RealValue();
+ if ( boost::math::isnan( aCoordZ ) || boost::math::isinf( aCoordZ ) )
+ aRes = false;
ProfilePoint aPoint( aCoordX, aCoordY, aCoordZ );
aPointsXYZ.Append( aPoint );
}
}
- aRes = aRes && ( anIsParametric && !aPointsUZ.empty() ||
+ aRes = aRes && ( anIsParametric && !aPointsUZ.IsEmpty() ||
anIsGeoref && !aPointsXYZ.IsEmpty() );
if ( aRes )
{
{
SetProfilePoints( aPointsXYZ );
}
+
+ Update();
}
return aRes;