src/HYDROGUI/HYDROGUI_ProfileInterpolateOp.cxx

   1 // Copyright (C) 2007-2015  CEA/DEN, EDF R&D, OPEN CASCADE
   2 //
   3 // Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   5 //
   6 // This library is free software; you can redistribute it and/or
   7 // modify it under the terms of the GNU Lesser General Public
   8 // License as published by the Free Software Foundation; either
   9 // version 2.1 of the License, or (at your option) any later version.
  10 //
  11 // This library is distributed in the hope that it will be useful,
  12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14 // Lesser General Public License for more details.
  15 //
  16 // You should have received a copy of the GNU Lesser General Public
  17 // License along with this library; if not, write to the Free Software
  18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  19 //
  20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  21 //
  22
  23 #include <HYDROGUI_ProfileInterpolateOp.h>
  24
  25 #include "HYDROGUI_Tool.h"
  26 #include "HYDROGUI_Module.h"
  27 #include "HYDROGUI_DataObject.h"
  28 #include "HYDROGUI_UpdateFlags.h"
  29 #include "HYDROGUI_ProfileInterpolateDlg.h"
  30
  31 #include <HYDROData_Stream.h>
  32 #include <HYDROData_Profile.h>
  33 #include <HYDROData_Document.h>
  34 #include <HYDROData_ProfileUZ.h>
  35 #include <HYDROData_IPolyline.h>
  36 #include <HYDROData_PolylineXY.h>
  37 #include <HYDROData_InterpolatorsFactory.h>
  38 #include <HYDROData_IProfilesInterpolator.h>
  39
  40 #include <LightApp_Application.h>
  41 #include <LightApp_SelectionMgr.h>
  42 #include <LightApp_UpdateFlags.h>
  43
  44 #include <OCCViewer_ViewModel.h>
  45 #include <OCCViewer_ViewPort3d.h>
  46 #include <OCCViewer_ViewWindow.h>
  47 #include <OCCViewer_ViewManager.h>
  48 #include <OCCViewer_AISSelector.h>
  49
  50 #include <AIS_ListOfInteractive.hxx>
  51 #include <AIS_ListIteratorOfListOfInteractive.hxx>
  52
  53 #include <BRep_Builder.hxx>
  54 #include <TopoDS_Compound.hxx>
  55
  56 #include <QMessageBox>
  57
  58 HYDROGUI_ProfileInterpolateOp::HYDROGUI_ProfileInterpolateOp( HYDROGUI_Module* theModule )
  59     : HYDROGUI_Operation( theModule )
  60 {
  61     setName( tr( "PROFILE_INTERPOLATION" ) );
  62 }
  63
  64 HYDROGUI_ProfileInterpolateOp::~HYDROGUI_ProfileInterpolateOp()
  65 {
  66 }
  67
  68 void HYDROGUI_ProfileInterpolateOp::startOperation()
  69 {
  70     HYDROGUI_Operation::startOperation();
  71
  72     HYDROGUI_ProfileInterpolateDlg* aPanel = ::qobject_cast<HYDROGUI_ProfileInterpolateDlg*>( inputPanel() );
  73     if ( aPanel )
  74     {
  75         aPanel->reset();
  76         aPanel->setInterpolators( interpolators() );
  77     }
  78 }
  79
  80 void HYDROGUI_ProfileInterpolateOp::abortOperation()
  81 {
  82     if ( !myPreview.IsNull() )
  83     {
  84         HYDROGUI_ProfileInterpolateDlg* aDlg = new HYDROGUI_ProfileInterpolateDlg( module(), getName() );
  85         if ( aDlg && !aDlg->getAISContext().IsNull() )
  86         {
  87             aDlg->getAISContext()->Remove( myPreview, false );
  88             myPreview.Nullify();
  89         }
  90     }
  91
  92     HYDROGUI_Operation::abortOperation();
  93 }
  94
  95 void HYDROGUI_ProfileInterpolateOp::commitOperation()
  96 {
  97     if ( !myPreview.IsNull() )
  98     {
  99         HYDROGUI_ProfileInterpolateDlg* aDlg = new HYDROGUI_ProfileInterpolateDlg( module(), getName() );
 100         if ( aDlg && !aDlg->getAISContext().IsNull() )
 101         {
 102             aDlg->getAISContext()->Remove( myPreview, false );
 103             myPreview.Nullify();
 104         }
 105     }
 106
 107     HYDROGUI_Operation::commitOperation();
 108 }
 109
 110 HYDROGUI_InputPanel* HYDROGUI_ProfileInterpolateOp::createInputPanel() const
 111 {
 112     HYDROGUI_ProfileInterpolateDlg* aDlg = new HYDROGUI_ProfileInterpolateDlg( module(), getName() );
 113
 114     connect( aDlg, SIGNAL( riverChanged( const QString& ) ), this, SLOT( onRiverChanged( const QString& ) ) );
 115     connect( aDlg, SIGNAL( interpolatorChanged( const QString& ) ), this, SLOT( onInterpolatorChanged( const QString& ) ) );
 116     connect( aDlg, SIGNAL( interpolatorParametersChanged( const QString& ) ), this, SLOT( updatePreview() ) );
 117     connect( aDlg, SIGNAL( profileNumberChanged( int ) ), this, SLOT( updatePreview() ) );
 118     connect( aDlg, SIGNAL( profileStartChanged( const QString& ) ), this, SLOT( updatePreview() ) );
 119     connect( aDlg, SIGNAL( profileFinishChanged( const QString& ) ), this, SLOT( updatePreview() ) );
 120
 121     return aDlg;
 122 }
 123
 124 bool HYDROGUI_ProfileInterpolateOp::processApply( int& theUpdateFlags, QString& theErrorMsg,
 125                                                   QStringList& theBrowseObjectsEntries )
 126 {
 127     HYDROGUI_ProfileInterpolateDlg* aDlg = ::qobject_cast<HYDROGUI_ProfileInterpolateDlg*>( inputPanel() );
 128     if ( !aDlg )
 129         return false;
 130
 131     QString errMessage;
 132     HYDROData_IProfilesInterpolator* anIterp = interpolator( aDlg->interpolator() );
 133     if ( anIterp )
 134     {
 135         updateInterpolator( anIterp );
 136         Handle(HYDROData_Stream) aRiver = Handle(HYDROData_Stream)::DownCast( HYDROGUI_Tool::FindObjectByName( module(), aDlg->river(), KIND_STREAM ) );
 137         if ( !aRiver.IsNull() )
 138         {
 139             startDocOperation();
 140             aRiver->Interpolate( anIterp );
 141             if ( anIterp->GetErrorCode() == OK )
 142                 commitDocOperation();
 143             else
 144             {
 145                 errMessage = tr( "CALCULATE_ERROR" ), QString( anIterp->GetErrorMessage().c_str() );
 146                 abortDocOperation();
 147             }
 148         }
 149         else
 150             errMessage = tr( "CANT_GET_RIVER_OBJECT" ).arg( aDlg->river() );
 151
 152         if ( !errMessage.isEmpty() )
 153             QMessageBox::critical( aDlg->topLevelWidget(), tr( "INTERPOLATION_ERROR" ), errMessage, QMessageBox::Ok );
 154         else
 155             commit();
 156     }
 157
 158     bool res = errMessage.isEmpty();
 159
 160     if ( res )
 161         theUpdateFlags = UF_Model | UF_OCCViewer | UF_OCC_Forced | UF_VTKViewer;
 162
 163     return res;
 164 }
 165
 166 void HYDROGUI_ProfileInterpolateOp::updatePreview()
 167 {
 168     HYDROGUI_ProfileInterpolateDlg* aDlg = ::qobject_cast<HYDROGUI_ProfileInterpolateDlg*>( inputPanel() );
 169     if ( !aDlg )
 170         return;
 171
 172     Handle(AIS_InteractiveContext) aCtx = aDlg->getAISContext();
 173     if ( aCtx.IsNull() )
 174         return;
 175
 176     if ( !myPreview.IsNull() )
 177         aCtx->Remove( myPreview, false );
 178
 179     myPreview.Nullify();
 180
 181     HYDROData_IProfilesInterpolator* anIterp = interpolator( aDlg->interpolator() );
 182     if ( anIterp )
 183     {
 184         updateInterpolator( anIterp );
 185         anIterp->Calculate();
 186         if ( anIterp->GetErrorCode() == OK )
 187         {
 188             TopoDS_Shape aShape = previewShape( anIterp );
 189             if ( !aShape.IsNull() )
 190             {
 191                 myPreview = new AIS_Shape( aShape );
 192                 myPreview->SetColor( Quantity_NOC_RED );
 193                 aCtx->Display( myPreview, 1, -1, false );
 194             }
 195         }
 196     }
 197
 198     aCtx->UpdateCurrentViewer();
 199 }
 200
 201 void HYDROGUI_ProfileInterpolateOp::updateInterpolator( HYDROData_IProfilesInterpolator* theInt )
 202 {
 203     HYDROGUI_ProfileInterpolateDlg* aDlg = ::qobject_cast<HYDROGUI_ProfileInterpolateDlg*>( inputPanel() );
 204     if ( !aDlg || !theInt )
 205         return;
 206
 207     theInt->Reset();
 208     theInt->SetResultProfilesNumber( aDlg->profileNumber() );
 209     theInt->SetProfiles( profile( aDlg->profileStart() ), profile( aDlg->profileFinish() ) );
 210     ParamsList aParams = parameters( aDlg->interpolatorParameters() );
 211     for ( ParamsList::Iterator it = aParams.begin(); it != aParams.end(); ++it )
 212         theInt->SetParameter( (*it).first.toStdString(), (*it).second.toStdString() );
 213 }
 214
 215 std::vector<double> HYDROGUI_ProfileInterpolateOp::profile( const QString& theName ) const
 216 {
 217     std::vector<double> aPoints;
 218     Handle(HYDROData_Profile) aProf = Handle(HYDROData_Profile)::DownCast( HYDROGUI_Tool::FindObjectByName( module(), theName, KIND_PROFILE ) );
 219     if ( !aProf.IsNull() )
 220     {
 221         HYDROData_Profile::ProfilePoints aPntList = aProf->GetProfilePoints();
 222         for ( int i = aPntList.Lower(); i <= aPntList.Upper(); i++ )
 223         {
 224             HYDROData_Profile::ProfilePoint aPnt = aPntList.Value( i );
 225             aPoints.push_back( aPnt.X() );
 226             aPoints.push_back( aPnt.Y() );
 227             aPoints.push_back( aPnt.Z() );
 228         }
 229     }
 230     return aPoints;
 231 }
 232
 233 HYDROGUI_ProfileInterpolateOp::ParamsList HYDROGUI_ProfileInterpolateOp::parameters( const QString& theStr ) const
 234 {
 235     // Regular expression for parsing parameters string of kind: <param_name> <param_value> <param_name> <param_value> ..."
 236     // QRegExp rx( "([A-Za-z_\\d]+)\\s+([A-Za-z_\\d]+)\\s*" );
 237
 238     // Regular expression for parsing parameters string of kind: <param_name> = <param_value> <param_name> = <param_value> ..."
 239     QRegExp rx( "([A-Za-z_\\d]+)\\s*=\\s*([A-Za-z_\\d]+)\\s*" );
 240
 241     ParamsList aParamList;
 242     for ( int i = 0; i >= 0; )
 243     {
 244         i = rx.indexIn( theStr, i );
 245         if ( i >= 0  )
 246         {
 247             aParamList.append( QPair<QString, QString>( rx.cap( 1 ), rx.cap( 2 ) ) );
 248             i += rx.matchedLength();
 249         }
 250     }
 251
 252     return aParamList;
 253 }
 254
 255 QString HYDROGUI_ProfileInterpolateOp::parameters( const ParamsList& theParamList ) const
 256 {
 257     QStringList paramList;
 258     for ( ParamsList::const_iterator it = theParamList.begin(); it != theParamList.end(); ++it )
 259         paramList.append( QString( "%1 = %2" ).arg( (*it).first ).arg( (*it).second ) );
 260     return paramList.join( " " );
 261 }
 262
 263 QStringList HYDROGUI_ProfileInterpolateOp::interpolators() const
 264 {
 265     HYDROData_InterpolatorsFactory* anIFactory = 0;
 266     Handle(HYDROData_Document) aDoc = doc();
 267     if ( !aDoc.IsNull() )
 268         anIFactory = aDoc->GetInterpolatorsFactory();
 269
 270     QStringList aNames;
 271     if ( anIFactory )
 272     {
 273         std::vector<std::string> iNames = anIFactory->GetInterpolatorNames();
 274         for ( std::vector<std::string>::iterator it = iNames.begin(); it != iNames.end(); ++it )
 275             aNames.append( (*it).c_str() );
 276     }
 277
 278     return aNames;
 279 }
 280
 281 HYDROData_IProfilesInterpolator* HYDROGUI_ProfileInterpolateOp::interpolator( const QString& theName ) const
 282 {
 283     HYDROData_InterpolatorsFactory* anIFactory = 0;
 284     Handle(HYDROData_Document) aDoc = doc();
 285     if ( !aDoc.IsNull() )
 286         anIFactory = aDoc->GetInterpolatorsFactory();
 287
 288     HYDROData_IProfilesInterpolator* aRes = 0;
 289     if ( anIFactory )
 290         aRes = anIFactory->GetInterpolator( theName.toStdString() );
 291     return aRes;
 292 }
 293
 294 TopoDS_Shape HYDROGUI_ProfileInterpolateOp::previewShape( HYDROData_IProfilesInterpolator* theInterp ) const
 295 {
 296     TopoDS_Compound aPreviewShape;
 297     if ( theInterp )
 298     {
 299         BRep_Builder aBuilder;
 300         aBuilder.MakeCompound( aPreviewShape );
 301         for ( int i = 0; i < theInterp->GetCalculatedProfilesNumber(); i++ )
 302         {
 303             NCollection_Sequence<gp_XYZ> pointSeq;
 304             std::vector<double> aPoints = theInterp->GetResultProfile( i );
 305             for ( int i = 0; i < (int) aPoints.size(); i += 3 )
 306                 pointSeq.Append( gp_XYZ( aPoints[i], aPoints[i+1], aPoints[i+2] ) );
 307             TopoDS_Shape aWire = HYDROData_PolylineXY::BuildWire( HYDROData_IPolyline::SECTION_SPLINE, false, pointSeq );
 308             if ( !aWire.IsNull() )
 309                 aBuilder.Add( aPreviewShape, aWire );
 310         }
 311     }
 312
 313     return aPreviewShape;
 314 }
 315
 316 void HYDROGUI_ProfileInterpolateOp::onInterpolatorChanged( const QString& theInterpName )
 317 {
 318     HYDROGUI_ProfileInterpolateDlg* aPanel = dynamic_cast<HYDROGUI_ProfileInterpolateDlg*>( inputPanel() );
 319     HYDROData_IProfilesInterpolator* anInterp = interpolator( theInterpName );
 320     if ( !aPanel || !anInterp )
 321         return;
 322
 323     aPanel->setInterpolatorDescription( QString( anInterp->GetDescription().c_str() ) );
 324
 325     updatePreview();
 326 }
 327
 328 void HYDROGUI_ProfileInterpolateOp::onRiverChanged( const QString& theRiver )
 329 {
 330     HYDROGUI_ProfileInterpolateDlg* aDlg = ::qobject_cast<HYDROGUI_ProfileInterpolateDlg*>( inputPanel() );
 331     if ( !aDlg )
 332         return;
 333
 334     Handle(AIS_InteractiveContext) aCtx = aDlg->getAISContext();
 335     if ( !aCtx.IsNull() )
 336     {
 337         AIS_ListOfInteractive aList;
 338         aCtx->DisplayedObjects( aList );
 339         for ( AIS_ListIteratorOfListOfInteractive it( aList ); it.More(); it.Next() )
 340         {
 341             Handle(AIS_Shape) anObj = Handle(AIS_Shape)::DownCast( it.Value() );
 342             if ( !anObj.IsNull() )
 343                 aCtx->Remove( anObj, false );
 344         }
 345
 346         Handle(HYDROData_Stream) aStream = Handle(HYDROData_Stream)::DownCast( HYDROGUI_Tool::FindObjectByName( module(), theRiver, KIND_STREAM ) );
 347         if ( !aStream.IsNull() )
 348         {
 349             HYDROData_SequenceOfObjects aSeq = aStream->GetProfiles();
 350             for ( int i = aSeq.Lower(); i <= aSeq.Upper(); i++ )
 351             {
 352                 Handle(HYDROData_Profile) aProfile = Handle(HYDROData_Profile)::DownCast( aSeq.Value( i ) );
 353                 if ( !aProfile.IsNull() )
 354                 {
 355                     Handle(AIS_Shape) aPrs = new AIS_Shape( aProfile->GetShape3D() );
 356                     aPrs->SetOwner( aProfile );
 357                     aPrs->SetColor( Quantity_NOC_BLACK );
 358                     aCtx->Display( aPrs, 1, 0, false );
 359                 }
 360             }
 361             OCCViewer_ViewManager* vm = aDlg->viewManager();
 362             if ( vm )
 363             {
 364                 QVector<SUIT_ViewWindow*> winList = vm->getViews();
 365                 for ( QVector<SUIT_ViewWindow*>::iterator it = winList.begin(); it != winList.end(); ++it )
 366                 {
 367                     OCCViewer_ViewWindow* occWin = ::qobject_cast<OCCViewer_ViewWindow*>( *it );
 368                     if ( occWin )
 369                     {
 370                         OCCViewer_ViewPort3d* vp = occWin->getViewPort();
 371                         if ( vp )
 372                             vp->fitAll();
 373                     }
 374                 }
 375             }
 376         }
 377     }
 378     updatePreview();
 379 }