#include <HYDROData_Stream.h>
#include <HYDROData_Profile.h>
+#include <HYDROData_Document.h>
+#include <HYDROData_ProfileUZ.h>
+#include <HYDROData_IPolyline.h>
+#include <HYDROData_PolylineXY.h>
+#include <HYDROData_InterpolatorsFactory.h>
+#include <HYDROData_IProfilesInterpolator.h>
#include <LightApp_Application.h>
#include <LightApp_SelectionMgr.h>
#include <AIS_ListOfInteractive.hxx>
#include <AIS_ListIteratorOfListOfInteractive.hxx>
+#include <BRep_Builder.hxx>
+#include <TopoDS_Compound.hxx>
+
+#include <QMessageBox>
+
HYDROGUI_ProfileInterpolateOp::HYDROGUI_ProfileInterpolateOp( HYDROGUI_Module* theModule )
: HYDROGUI_Operation( theModule )
{
aPanel->reset();
aPanel->setInterpolators( interpolators() );
}
+}
-/*
- if( myIsEdit )
- myEditedObject = Handle(HYDROData_Profile)::DownCast( HYDROGUI_Tool::GetSelectedObject( module() ) );
-
- QString aProfileName;
- if( !myEditedObject.IsNull() )
+void HYDROGUI_ProfileInterpolateOp::abortOperation()
+{
+ if ( !myPreview.IsNull() )
{
- Handle(HYDROData_ProfileUZ) aProfileUZ = myEditedObject->GetProfileUZ( false );
- if ( !aProfileUZ.IsNull() )
+ HYDROGUI_ProfileInterpolateDlg* aDlg = new HYDROGUI_ProfileInterpolateDlg( module(), getName() );
+ if ( aDlg && !aDlg->getAISContext().IsNull() )
{
- CurveCreator::Coordinates aCurveCoords;
- CurveCreator::SectionsMap aSectionsMap;
-
- HYDROData_ProfileUZ::PointsList aSectPointsList = aProfileUZ->GetPoints();
- CurveCreator::PosPointsList aPoints;
- for ( int k = 1, aNbPoints = aSectPointsList.Size(); k <= aNbPoints; ++k )
- {
- const HYDROData_ProfileUZ::Point& aSectPoint = aSectPointsList.Value( k );
- aCurveCoords.clear();
- aCurveCoords.push_back( aSectPoint.X() );
- aCurveCoords.push_back( aSectPoint.Y() );
-
- CurveCreator_PosPoint* aPosPoint = new CurveCreator_PosPoint( k, aCurveCoords );
- aPoints.push_back( aPosPoint );
- }
-
- aSectionsMap[0] = aPoints;
- myProfile->addPointsInternal( aSectionsMap );
-
- HYDROData_ProfileUZ::SectionType aSectType = aProfileUZ->GetSectionType( 0 );
-
- CurveCreator::SectionType aCurveType = CurveCreator::Polyline;
- if( aSectType == HYDROData_ProfileUZ::SECTION_SPLINE )
- aCurveType = CurveCreator::Spline;
-
- myProfile->setSectionType( 0, aCurveType );
+ aDlg->getAISContext()->Remove( myPreview, false );
+ myPreview.Nullify();
}
-
- aProfileName = myEditedObject->GetName();
- }
- else
- {
- aProfileName = HYDROGUI_Tool::GenerateObjectName( module(), tr( "DEFAULT_PROFILE_NAME" ) );
}
- aPanel->setProfileName( aProfileName );
- aPanel->setProfile( myProfile );
-*/
-// displayPreview();
-}
-
-void HYDROGUI_ProfileInterpolateOp::abortOperation()
-{
HYDROGUI_Operation::abortOperation();
}
void HYDROGUI_ProfileInterpolateOp::commitOperation()
{
+ if ( !myPreview.IsNull() )
+ {
+ HYDROGUI_ProfileInterpolateDlg* aDlg = new HYDROGUI_ProfileInterpolateDlg( module(), getName() );
+ if ( aDlg && !aDlg->getAISContext().IsNull() )
+ {
+ aDlg->getAISContext()->Remove( myPreview, false );
+ myPreview.Nullify();
+ }
+ }
+
HYDROGUI_Operation::commitOperation();
}
bool HYDROGUI_ProfileInterpolateOp::processApply( int& theUpdateFlags, QString& theErrorMsg,
QStringList& theBrowseObjectsEntries )
{
- HYDROGUI_ProfileInterpolateDlg* aPanel = ::qobject_cast<HYDROGUI_ProfileInterpolateDlg*>( inputPanel() );
- if ( !aPanel )
- return false;
-/*
- QString aProfileName = aPanel->getProfileName().simplified();
- if ( aProfileName.isEmpty() )
- {
- theErrorMsg = tr( "INCORRECT_OBJECT_NAME" );
+ HYDROGUI_ProfileInterpolateDlg* aDlg = ::qobject_cast<HYDROGUI_ProfileInterpolateDlg*>( inputPanel() );
+ if ( !aDlg )
return false;
- }
- if( !myIsEdit || ( !myEditedObject.IsNull() && myEditedObject->GetName() != aProfileName ) )
+ QString errMessage;
+ HYDROData_IProfilesInterpolator* anIterp = interpolator( aDlg->interpolator() );
+ if ( anIterp )
{
- // check that there are no other objects with the same name in the document
- Handle(HYDROData_Entity) anObject = HYDROGUI_Tool::FindObjectByName( module(), aProfileName );
- if( !anObject.IsNull() )
+ updateInterpolator( anIterp );
+ Handle(HYDROData_Stream) aRiver = Handle(HYDROData_Stream)::DownCast( HYDROGUI_Tool::FindObjectByName( module(), aDlg->river(), KIND_STREAM ) );
+ if ( !aRiver.IsNull() )
{
- theErrorMsg = tr( "OBJECT_EXISTS_IN_DOCUMENT" ).arg( aProfileName );
- return false;
+ startDocOperation();
+ //aRiver->Interpolate( anIterp );
+ if ( anIterp->GetErrorCode() == OK )
+ commitDocOperation();
+ else
+ {
+ errMessage = tr( "CALCULATE_ERROR" ), QString( anIterp->GetErrorMessage().c_str() );
+ abortDocOperation();
+ }
}
- }
+ else
+ errMessage = tr( "CANT_GET_RIVER_OBJECT" ).arg( aDlg->river() );
- Handle(HYDROData_Profile) aProfileObj;
- if( myIsEdit ){
- aProfileObj = myEditedObject;
- }
- else{
- aProfileObj = Handle(HYDROData_Profile)::DownCast( doc()->CreateObject( KIND_PROFILE ) );
+ if ( !errMessage.isEmpty() )
+ QMessageBox::critical( aDlg->topLevelWidget(), tr( "INTERPOLATION_ERROR" ), errMessage, QMessageBox::Ok );
+ else
+ commit();
}
- if( aProfileObj.IsNull() )
- return false;
+ return errMessage.isEmpty();
+}
- Handle(HYDROData_ProfileUZ) aProfileUZ = aProfileObj->GetProfileUZ();
- if ( aProfileUZ.IsNull() )
- return false;
+void HYDROGUI_ProfileInterpolateOp::updatePreview()
+{
+ HYDROGUI_ProfileInterpolateDlg* aDlg = ::qobject_cast<HYDROGUI_ProfileInterpolateDlg*>( inputPanel() );
+ if ( !aDlg )
+ return;
- aProfileObj->SetName(aProfileName);
+ Handle(AIS_InteractiveContext) aCtx = aDlg->getAISContext();
+ if ( aCtx.IsNull() )
+ return;
+
+ if ( !myPreview.IsNull() )
+ aCtx->Remove( myPreview, false );
- HYDROData_ProfileUZ::PointsList aProfileParamPoints;
+ myPreview.Nullify();
- CurveCreator::Coordinates aCurveCoords = myProfile->getPoints( 0 );
- if ( aCurveCoords.size() <= 2 )
+ HYDROData_IProfilesInterpolator* anIterp = interpolator( aDlg->interpolator() );
+ if ( anIterp )
{
- theErrorMsg = tr( "NUMBER_OF_PROFILE_POINTS_INCORRECT" );
- return false;
+ updateInterpolator( anIterp );
+ if ( anIterp->Calculate() == OK )
+ {
+ TopoDS_Shape aShape = previewShape( anIterp, aDlg->profileNumber() );
+ if ( !aShape.IsNull() )
+ {
+ myPreview = new AIS_Shape( aShape );
+ myPreview->SetColor( Quantity_NOC_RED );
+ aCtx->Display( myPreview, 0, -1, false );
+ }
+ }
}
- for ( int k = 0 ; k + 1 < aCurveCoords.size() ; k++ )
- {
- HYDROData_ProfileUZ::Point aProfileParamPoint;
+ aCtx->UpdateCurrentViewer();
+}
- aProfileParamPoint.SetX( aCurveCoords.at( k ) );
- k++;
- aProfileParamPoint.SetY( aCurveCoords.at( k ) );
+void HYDROGUI_ProfileInterpolateOp::updateInterpolator( HYDROData_IProfilesInterpolator* theInt )
+{
+ HYDROGUI_ProfileInterpolateDlg* aDlg = ::qobject_cast<HYDROGUI_ProfileInterpolateDlg*>( inputPanel() );
+ if ( !aDlg || !theInt )
+ return;
+
+ theInt->Reset();
+ theInt->SetResultProfilesNumber( aDlg->profileNumber() );
+ theInt->SetProfiles( profile( aDlg->profileStart() ), profile( aDlg->profileFinish() ) );
+ ParamsList aParams = parameters( aDlg->interpolatorParameters() );
+ for ( ParamsList::Iterator it = aParams.begin(); it != aParams.end(); ++it )
+ theInt->SetParameter( (*it).first.toStdString(), (*it).second.toStdString() );
+}
- aProfileParamPoints.Append( aProfileParamPoint );
+std::vector<double> HYDROGUI_ProfileInterpolateOp::profile( const QString& theName ) const
+{
+ std::vector<double> aPoints;
+ Handle(HYDROData_Profile) aProf = Handle(HYDROData_Profile)::DownCast( HYDROGUI_Tool::FindObjectByName( module(), theName, KIND_PROFILE ) );
+ if ( !aProf.IsNull() )
+ {
+ HYDROData_Profile::ProfilePoints aPntList = aProf->GetProfilePoints();
+ for ( int i = aPntList.Lower(); i <= aPntList.Upper(); i++ )
+ {
+ HYDROData_Profile::ProfilePoint aPnt = aPntList.Value( i );
+ aPoints.push_back( aPnt.X() );
+ aPoints.push_back( aPnt.Y() );
+ aPoints.push_back( aPnt.Z() );
+ }
}
- aProfileObj->SetParametricPoints( aProfileParamPoints );
+ return aPoints;
+}
- HYDROData_ProfileUZ::SectionType aSectType = HYDROData_ProfileUZ::SECTION_POLYLINE;
- if ( myProfile->getSectionType( 0 ) == CurveCreator::Spline )
- aSectType = HYDROData_ProfileUZ::SECTION_SPLINE;
+HYDROGUI_ProfileInterpolateOp::ParamsList HYDROGUI_ProfileInterpolateOp::parameters( const QString& theStr ) const
+{
+ // Regular expression for parsing parameters string of kind: <param_name> <param_value> <param_name> <param_value> ..."
+ // QRegExp rx( "([A-Za-z_\\d]+)\\s+([A-Za-z_\\d]+)\\s*" );
- aProfileUZ->SetSectionType( 0, aSectType );
+ // Regular expression for parsing parameters string of kind: <param_name> = <param_value> <param_name> = <param_value> ..."
+ QRegExp rx( "([A-Za-z_\\d]+)\\s*=\\s*([A-Za-z_\\d]+)\\s*" );
- if ( !myIsEdit )
+ ParamsList aParamList;
+ for ( int i = 0; i >= 0; )
{
- aProfileObj->SetBorderColor( HYDROData_Profile::DefaultBorderColor() );
+ i = rx.indexIn( theStr, i );
+ if ( i >= 0 )
+ {
+ aParamList.append( QPair<QString, QString>( rx.cap( 1 ), rx.cap( 2 ) ) );
+ i += rx.matchedLength();
+ }
}
- // At first we update the child u,z profile object
- aProfileUZ->SetToUpdate( true );
- aProfileUZ->Update();
+ return aParamList;
+}
+
+QString HYDROGUI_ProfileInterpolateOp::parameters( const ParamsList& theParamList ) const
+{
+ QStringList paramList;
+ for ( ParamsList::const_iterator it = theParamList.begin(); it != theParamList.end(); ++it )
+ paramList.append( QString( "%1 = %2" ).arg( (*it).first ).arg( (*it).second ) );
+ return paramList.join( " " );
+}
- // And now we update our edited object
- aProfileObj->Update();
- module()->setIsToUpdate( aProfileObj );
+QStringList HYDROGUI_ProfileInterpolateOp::interpolators() const
+{
+ HYDROData_InterpolatorsFactory* anIFactory = 0;
+ Handle(HYDROData_Document) aDoc = doc();
+ if ( !aDoc.IsNull() )
+ anIFactory = aDoc->GetInterpolatorsFactory();
- theUpdateFlags = UF_Model;
- if ( myIsEdit )
- theUpdateFlags |= UF_OCCViewer | UF_OCC_Forced | UF_VTKViewer;
- else
+ QStringList aNames;
+ if ( anIFactory )
{
- QString anEntry = HYDROGUI_DataObject::dataObjectEntry( aProfileObj );
- theBrowseObjectsEntries.append( anEntry );
+ std::vector<std::string> iNames = anIFactory->GetInterpolatorNames();
+ for ( std::vector<std::string>::iterator it = iNames.begin(); it != iNames.end(); ++it )
+ aNames.append( (*it).c_str() );
}
- */
- return true;
+ return aNames;
}
-void HYDROGUI_ProfileInterpolateOp::updatePreview()
+HYDROData_IProfilesInterpolator* HYDROGUI_ProfileInterpolateOp::interpolator( const QString& theName ) const
{
- HYDROGUI_ProfileInterpolateDlg* aDlg = ::qobject_cast<HYDROGUI_ProfileInterpolateDlg*>( inputPanel() );
- if ( !aDlg )
- return;
-
- Handle(AIS_InteractiveContext) aCtx = aDlg->getAISContext();
- if ( aCtx.IsNull() )
- return;
-
- aCtx->UpdateCurrentViewer();
+ HYDROData_InterpolatorsFactory* anIFactory = 0;
+ Handle(HYDROData_Document) aDoc = doc();
+ if ( !aDoc.IsNull() )
+ anIFactory = aDoc->GetInterpolatorsFactory();
+
+ HYDROData_IProfilesInterpolator* aRes = 0;
+ if ( anIFactory )
+ aRes = anIFactory->GetInterpolator( theName.toStdString() );
+ return aRes;
}
-QStringList HYDROGUI_ProfileInterpolateOp::interpolators() const
+TopoDS_Shape HYDROGUI_ProfileInterpolateOp::previewShape( HYDROData_IProfilesInterpolator* theInterp, int theNumber ) const
{
- return QStringList() << "Linear interpolator" << "Spline interpolator";
+ TopoDS_Compound aPreviewShape;
+ if ( theInterp )
+ {
+ BRep_Builder aBuilder;
+ aBuilder.MakeCompound( aPreviewShape );
+ for ( int i = 0; i < theNumber; i++ )
+ {
+ NCollection_Sequence<gp_XYZ> pointSeq;
+ std::vector<double> aPoints = theInterp->GetResultProfile( i );
+ for ( int i = 0; i < aPoints.size(); i += 3 )
+ pointSeq.Append( gp_XYZ( aPoints[i], aPoints[i+1], aPoints[i+2] ) );
+ TopoDS_Shape aWire = HYDROData_PolylineXY::BuildWire( HYDROData_IPolyline::SECTION_SPLINE, false, pointSeq );
+ if ( !aWire.IsNull() )
+ aBuilder.Add( aWire, aPreviewShape );
+ }
+ }
+
+ return aPreviewShape;
}
void HYDROGUI_ProfileInterpolateOp::onInterpolatorChanged( const QString& theInterpName )
{
HYDROGUI_ProfileInterpolateDlg* aPanel = dynamic_cast<HYDROGUI_ProfileInterpolateDlg*>( inputPanel() );
- if ( !aPanel )
+ HYDROData_IProfilesInterpolator* anInterp = interpolator( theInterpName );
+ if ( !aPanel || !anInterp )
return;
- aPanel->setInterpolatorDescription( QString( "Profile will be calculated used %1" ).arg( theInterpName ) );
+ aPanel->setInterpolatorDescription( QString( anInterp->GetDescription().c_str() ) );
updatePreview();
}