src/EntityGUI/EntityGUI_FeatureDetectorDlg.cxx

   1 // Copyright (C) 2007-2013  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.
  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 // GEOM GEOMGUI : GUI for Geometry component
  24 // File   : EntityGUI_SketcherDlg.cxx
  25 // Author : Renaud NEDELEC, Open CASCADE S.A.S.
  26
  27 // SALOME includes
  28 #include "EntityGUI_FeatureDetectorDlg.h"
  29 #include <ShapeRec_FeatureDetector.hxx>
  30
  31 #include <OCCViewer_ViewWindow.h>
  32 #include <OCCViewer_ViewManager.h>
  33
  34 #include <SOCC_ViewModel.h>
  35
  36 #include <DlgRef.h>
  37 #include <GeometryGUI.h>
  38 #include <EntityGUI.h>
  39 #include <GEOMBase.h>
  40
  41 #include <SUIT_Desktop.h>
  42 #include <SUIT_ResourceMgr.h>
  43 #include <SUIT_Session.h>
  44 #include <SUIT_ViewWindow.h>
  45 #include <SalomeApp_Application.h>
  46 #include <LightApp_Application.h>
  47 #include <LightApp_SelectionMgr.h>
  48
  49 #include <SALOME_ListIteratorOfListIO.hxx>
  50
  51 #include <SalomeApp_Study.h>
  52
  53 #include <utilities.h>
  54 #include <Precision.hxx>
  55
  56 // OCCT includes
  57 #include <gp_Pnt.hxx>
  58 #include <TopoDS_Shape.hxx>
  59 #include <TopoDS_Wire.hxx>
  60
  61 #include <BRepBuilderAPI_MakeVertex.hxx>
  62 #include <BRepBuilderAPI_MakeFace.hxx>
  63 #include <BRepBuilderAPI_MakePolygon.hxx>
  64 #include <BRepBuilderAPI_Transform.hxx>
  65
  66 #include <AIS_TexturedShape.hxx>
  67 #include <TCollection_AsciiString.hxx>
  68 #include <StdSelect_DisplayMode.hxx>
  69
  70 #include <Graphic3d_MaterialAspect.hxx>
  71
  72 // C++ includes
  73 #include <set>
  74 #include <utility>
  75
  76 // boost includes
  77 #include <boost/utility.hpp>
  78
  79 // Constructors
  80 enum{
  81   CONTOURS,
  82   CORNERS,
  83   LINES
  84 };
  85
  86 enum {
  87   RADIO_BUTTONS,
  88   MSG,
  89   PUSH_BUTTON,
  90 };
  91
  92 enum {
  93   KERNEL_SIZE,
  94   FIND_CONTOURS_METHOD,
  95   QUALITY_LEVEL,
  96   MIN_DISTANCE,
  97   TYPE_CRITERIA,
  98   MAX_ITER,
  99   EPSILON,
 100   LOW_THRESHOLD,
 101   RATIO,
 102   L2GRADIENT,
 103   SMOOTH_SIZE,
 104   HBINS,
 105   SBINS,
 106   HIST_TYPE,
 107   THRESHOLD_VALUE,
 108   MAX_THRESHOLD,
 109 };
 110
 111 // // // View
 112 // // enum {
 113 // //   XY,
 114 // //   YZ,
 115 // //   XZ
 116 // // };
 117
 118
 119 //=================================================================================
 120 // class    : EntityGUI_FeatureDetectorDlg()
 121 // purpose  : Constructs a EntityGUI_FeatureDetectorDlg which is a child of 'parent', with the
 122 //            name 'name' and widget flags set to 'f'.
 123 //            The dialog will by default be modeless, unless you set 'modal' to
 124 //            TRUE to construct a modal dialog.
 125 //=================================================================================
 126 EntityGUI_FeatureDetectorDlg::EntityGUI_FeatureDetectorDlg( GeometryGUI* theGeometryGUI, QWidget* parent,
 127                                               bool modal, Qt::WindowFlags fl)
 128   : GEOMBase_Skeleton(theGeometryGUI, parent, modal, fl)
 129 {
 130   QPixmap image0(SUIT_Session::session()->resourceMgr()->loadPixmap("GEOM", tr("ICON_SELECT")));
 131   QPixmap image1(SUIT_Session::session()->resourceMgr()->loadPixmap("OCCViewer", tr("ICON_OCCVIEWER_VIEW_DUMP")));
 132
 133   setWindowTitle(tr("GEOM_DETECT_TITLE"));
 134
 135   /***************************************************************/
 136
 137   mainFrame()->GroupConstructors->setTitle(tr("GEOM_FEATURES"));
 138   mainFrame()->RadioButton1->setText(tr("GEOM_CONTOURS"));
 139   mainFrame()->RadioButton2->setText(tr("GEOM_CORNERS"));
 140   mainFrame()->RadioButton3->setAttribute(Qt::WA_DeleteOnClose);
 141   mainFrame()->RadioButton3->close();
 142
 143 //   myViewGroup = new DlgRef_3Radio(centralWidget());
 144 //   myViewGroup->GroupBox1->setTitle(tr("GEOM_VIEW"));
 145 //   myViewGroup->RadioButton1->setText(tr( "GEOM_TOP"  ));
 146 //   myViewGroup->RadioButton2->setText(tr( "GEOM_FRONT"));
 147 //   myViewGroup->RadioButton3->setText(tr( "GEOM_LEFT" ));
 148 //   myViewButtonGroup = new QButtonGroup( this );
 149 //   myViewButtonGroup->addButton( myViewGroup->RadioButton1, XY ); // Top view
 150 //   myViewButtonGroup->addButton( myViewGroup->RadioButton2, YZ ); // Front View
 151 //   myViewButtonGroup->addButton( myViewGroup->RadioButton3, XZ ); // Left View
 152 //
 153 //   myViewGroup->hide();
 154
 155   // Widgets for the selection of the picture and the Region Of Interest
 156
 157
 158   mySelectionGroup = new DlgRef_1Sel1Frame(centralWidget());
 159
 160   mySelectionGroup->PushButton1->setIcon(image0);
 161   mySelectionGroup->PushButton1->setCheckable(true);
 162   mySelectionGroup->PushButton1->setAutoExclusive(true);
 163
 164   mySelectionGroup->PushButton2->setIcon(image1);
 165   mySelectionGroup->PushButton2->setCheckable(true);
 166   mySelectionGroup->PushButton2->setAutoExclusive(true);
 167
 168   mySelectionGroup->TextLabel1->setText(tr( "GEOM_PICTURE" ));
 169   mySelectionGroup->FrameLabel->setText("");
 170
 171   myOutputGroup = new DlgRef_3Radio(centralWidget());
 172   myOutputGroup->GroupBox1->setTitle(tr("GEOM_DETECT_OUTPUT"));
 173   myOutputGroup->RadioButton2->setText(tr( "GEOM_POLYLINE"  ));
 174   myOutputGroup->RadioButton1->setText(tr( "GEOM_SPLINE"));
 175   myOutputGroup->RadioButton3->hide();
 176
 177 //   myOutputGroup->hide(); //caché pour la demo
 178
 179   // Parameters
 180   QGridLayout* parametersLayout = new QGridLayout();
 181   parametersLayout->setMargin(9);
 182   QComboBox* kernelSize = new QComboBox();
 183   kernelSize->addItems( QStringList() << "3" << "5" << "7" );
 184   myWidgets.insert( KERNEL_SIZE, kernelSize );
 185   myUseROI = new QCheckBox( tr("USE_ROI") );
 186   myUseROI->setChecked( true );
 187
 188   parametersLayout->addWidget( new QLabel( tr("KERNEL_SIZE") ), 0, 0 );
 189   parametersLayout->addWidget( kernelSize, 0, 1 );
 190   parametersLayout->addWidget( myUseROI, 0, 2 );
 191
 192   // Corners Parameters
 193   myCornersParameters = new QFrame();
 194   myCornersParameters->setFrameStyle(QFrame::NoFrame);
 195   QDoubleSpinBox* qualityLevel = new QDoubleSpinBox();
 196   qualityLevel->setValue( 0.2 );
 197   qualityLevel->setSingleStep( 0.1 );
 198   qualityLevel->setRange( 0.01, 0.99 );
 199   myWidgets.insert( QUALITY_LEVEL, qualityLevel );
 200   QDoubleSpinBox* minDistance = new QDoubleSpinBox();
 201   minDistance->setValue( 1 );
 202   myWidgets.insert( MIN_DISTANCE, minDistance );
 203   QComboBox* typeCriteria = new QComboBox();
 204   typeCriteria->addItems( QStringList() << tr("CV_TERMCRIT_ITER") << tr("CV_TERMCRIT_EPS") << tr("CV_TERMCRIT_ITER | CV_TERMCRIT_EPS") );
 205   typeCriteria->setCurrentIndex( 2 );
 206   myWidgets.insert( TYPE_CRITERIA, typeCriteria );
 207   QSpinBox* maxIter = new QSpinBox();
 208   maxIter->setValue( 20 );
 209   maxIter->setRange( 1, 100 );
 210   myWidgets.insert( MAX_ITER, maxIter );
 211   QDoubleSpinBox* epsilon = new QDoubleSpinBox();
 212   epsilon->setValue( 0.03 );
 213   epsilon->setSingleStep( 0.01 );
 214   myWidgets.insert( EPSILON, epsilon );
 215
 216   QGridLayout* cornersLayout = new QGridLayout();
 217   cornersLayout->setMargin( 0 );
 218   cornersLayout->addWidget( new QLabel( tr("QUALITY_LEVEL") ), 0, 0 );
 219   cornersLayout->addWidget( qualityLevel, 0, 1 );
 220   cornersLayout->addWidget( new QLabel( tr("MIN_DISTANCE") ), 0, 2 );
 221   cornersLayout->addWidget( minDistance, 0 ,3 );
 222   cornersLayout->addWidget( new QLabel( tr("TYPE_CRITERIA") ), 1, 0 );
 223   cornersLayout->addWidget( typeCriteria, 1, 1, 1, 3 );
 224   cornersLayout->addWidget( new QLabel( tr("MAX_ITER") ), 2, 0 );
 225   cornersLayout->addWidget( maxIter, 2, 1 );
 226   cornersLayout->addWidget( new QLabel( tr("EPSILON") ), 2, 2 );
 227   cornersLayout->addWidget( epsilon, 2 , 3 );
 228
 229   myCornersParameters->setLayout( cornersLayout );
 230
 231   // Contours Parameters
 232   myContoursParameters = new QFrame();
 233   myContoursParameters->setFrameStyle(QFrame::NoFrame);
 234
 235   QComboBox* findContoursMethod = new QComboBox();
 236   findContoursMethod->addItems( QStringList() << tr("CV_CHAIN_APPROX_NONE") << tr("CV_CHAIN_APPROX_SIMPLE") << tr("CV_CHAIN_APPROX_TC89_L1") << tr("CV_CHAIN_APPROX_TC89_KCOS") );
 237   myWidgets.insert( FIND_CONTOURS_METHOD, findContoursMethod );
 238
 239   myCannyParameters = new QFrame();
 240   myCannyParameters->setFrameStyle(QFrame::NoFrame);
 241   QSpinBox* threshold = new QSpinBox();
 242   threshold->setRange( 0, 255 );
 243   threshold->setValue( 100 );
 244   myWidgets.insert( LOW_THRESHOLD, threshold );
 245   QSpinBox* ratio = new QSpinBox();
 246   ratio->setValue( 3 );
 247   ratio->setRange( 0, 255 );
 248   myWidgets.insert( RATIO, ratio );
 249   QCheckBox* L2gradient = new QCheckBox(tr("L2GRADIENT"));
 250   L2gradient->setChecked(true);
 251   myWidgets.insert( L2GRADIENT, L2gradient );
 252
 253   QGridLayout* cannyLayout = new QGridLayout();
 254   cannyLayout->setMargin( 0 );
 255   cannyLayout->addWidget( new QLabel( tr("LOWTHRESHOLD") ), 0, 0 );
 256   cannyLayout->addWidget( threshold, 0, 1 );
 257   cannyLayout->addWidget( new QLabel( tr("RATIO") ), 0, 2 );
 258   cannyLayout->addWidget( ratio, 0, 3 );
 259   cannyLayout->addWidget( L2gradient, 1, 0 );
 260
 261   myCannyParameters->setLayout( cannyLayout );
 262   myCannyParameters->setHidden( true );
 263
 264   myColorFilterParameters = new QFrame();
 265   myColorFilterParameters->setFrameStyle(QFrame::NoFrame);
 266
 267   QSpinBox* smoothSize = new QSpinBox();
 268   smoothSize->setValue( 3 );
 269   smoothSize->setSingleStep( 2 );
 270   smoothSize->setRange( 1, 1000 );
 271   myWidgets.insert( SMOOTH_SIZE, smoothSize );
 272   QSpinBox* hbins = new QSpinBox();
 273   hbins->setValue( 30 );
 274   hbins->setRange( 0, 179 );
 275   myWidgets.insert( HBINS, hbins );
 276   QSpinBox* sbins = new QSpinBox();
 277   sbins->setValue( 32 );
 278   sbins->setRange( 0, 255 );
 279   myWidgets.insert( SBINS, sbins );
 280   QComboBox* histType = new QComboBox();
 281   histType->addItems( QStringList() << tr("CV_HIST_ARRAY") << tr("CV_HIST_SPARSE") );
 282   myWidgets.insert( HIST_TYPE, histType );
 283   QDoubleSpinBox* thresholdValue = new QDoubleSpinBox();
 284   thresholdValue->setRange( 0, 254 );
 285   thresholdValue->setValue( 128 );
 286   myWidgets.insert( THRESHOLD_VALUE, thresholdValue );
 287   QDoubleSpinBox* maxThreshold = new QDoubleSpinBox();
 288   maxThreshold->setRange( 1, 255 );
 289   maxThreshold->setValue( 255 );
 290   myWidgets.insert( MAX_THRESHOLD, maxThreshold );
 291
 292   QGridLayout* colorFilterLayout = new QGridLayout();
 293   colorFilterLayout->setMargin( 0 );
 294   colorFilterLayout->addWidget( new QLabel( tr("SMOOTH_SIZE") ), 0, 0 );
 295   colorFilterLayout->addWidget( smoothSize, 0, 1 );
 296   colorFilterLayout->addWidget( new QLabel( tr("HBINS") ), 1, 0 );
 297   colorFilterLayout->addWidget( hbins, 1, 1 );
 298   colorFilterLayout->addWidget( new QLabel( tr("SBINS") ), 1, 2 );
 299   colorFilterLayout->addWidget( sbins, 1, 3 );
 300   colorFilterLayout->addWidget( new QLabel( tr("HIST_TYPE") ), 2, 0 );
 301   colorFilterLayout->addWidget( histType, 2, 1, 1, 3 );
 302   colorFilterLayout->addWidget( new QLabel( tr("THRESHOLD_VALUE") ), 3, 0 );
 303   colorFilterLayout->addWidget( thresholdValue, 3, 1 );
 304   colorFilterLayout->addWidget( new QLabel( tr("MAX_THRESHOLD") ), 3, 2 );
 305   colorFilterLayout->addWidget( maxThreshold, 3, 3 );
 306
 307   myColorFilterParameters->setLayout( colorFilterLayout );
 308   QGridLayout* contoursLayout = new QGridLayout();
 309   contoursLayout->setMargin( 0 );
 310   contoursLayout->addWidget( new QLabel( tr("FIND_CONTOURS_METHOD") ), 0, 0 );
 311   contoursLayout->addWidget( findContoursMethod, 0, 1 );
 312   contoursLayout->addWidget( myCannyParameters, 1, 0, 1, 3 );
 313   contoursLayout->addWidget( myColorFilterParameters, 1, 0, 1, 3  );
 314
 315   myContoursParameters->setLayout( contoursLayout );
 316
 317   parametersLayout->addWidget( myCornersParameters, 1, 0, 1, 3 );
 318   parametersLayout->addWidget( myContoursParameters, 1, 0, 1, 3 );
 319
 320   QVBoxLayout* layout = new QVBoxLayout(centralWidget());
 321   layout->setMargin(0); layout->setSpacing(6);
 322 //   layout->addWidget( myViewGroup);
 323   layout->addLayout( parametersLayout );
 324   layout->addWidget( mySelectionGroup);
 325   layout->addWidget( myOutputGroup);
 326
 327 //   mainFrame()->GroupBoxName->hide();
 328
 329   // Build an instance of detection used to perform image processing operations
 330   aDetector = new ShapeRec_FeatureDetector();
 331
 332   setHelpFileName( "shape_recognition_page.html" );
 333
 334   Init();
 335
 336 }
 337
 338 //=================================================================================
 339 // function : ~EntityGUI_FeatureDetectorDlg()
 340 // purpose  : Destroys the object and frees any allocated resources
 341 //=================================================================================
 342 EntityGUI_FeatureDetectorDlg::~EntityGUI_FeatureDetectorDlg()
 343 {
 344
 345 }
 346
 347 //=================================================================================
 348 // function : Init()
 349 // purpose  :
 350 //=================================================================================
 351 void EntityGUI_FeatureDetectorDlg::Init()
 352 {
 353   connect( myGeomGUI,         SIGNAL( SignalCloseAllDialogs() ), this, SLOT( ClickOnCancel() ) );
 354   connect( buttonOk(),        SIGNAL( clicked() ),               this, SLOT( ClickOnOk() ) );
 355   connect( buttonApply(),     SIGNAL( clicked() ),               this, SLOT( ClickOnApply() ) );
 356   connect( myUseROI,          SIGNAL( toggled(bool) ),           this, SLOT( onCheckBoxClicked(bool) ) );
 357   connect( this,              SIGNAL(constructorsClicked(int)),  this, SLOT(ConstructorsClicked(int)));
 358   connect( mySelectionGroup->PushButton2,      SIGNAL( clicked() ),               this, SLOT( onButtonClicked() ) );
 359   connect( mySelectionGroup->PushButton1,       SIGNAL( clicked() ),               this, SLOT( onButtonClicked() ) );
 360 //   connect( myViewButtonGroup, SIGNAL( buttonClicked( int ) ),    this, SLOT( onViewClicked( int ) ) );
 361   connect( myGeomGUI->getApp()->selectionMgr(), SIGNAL( currentSelectionChanged() ),this, SLOT( SelectionIntoArgument() ) );
 362
 363   myConstructorId = 0;
 364   mySelectionGroup->PushButton1->click();
 365
 366 //   SetEditCurrentArgument();
 367   SelectionIntoArgument();
 368
 369   initName(tr("GEOM_CONTOURS"));
 370   resize(100,100);
 371
 372 //   myViewGroup->RadioButton1->setChecked(true);
 373   myOutputGroup->RadioButton1->setChecked(true);
 374
 375   gp_Pnt aOrigin = gp_Pnt(0, 0, 0);
 376   gp_Dir aDirZ = gp_Dir(0, 0, 1);
 377   gp_Dir aDirX = gp_Dir(1, 0, 0);
 378
 379   aGlobalCS = gp_Ax3(aOrigin, aDirZ, aDirX);
 380
 381   myStartPnt = gp_Pnt(0,0,0);
 382   myEndPnt = myStartPnt;
 383
 384   myGeomGUI->SetWorkingPlane( aGlobalCS );
 385   myGeomGUI->ActiveWorkingPlane();
 386
 387   ConstructorsClicked(myConstructorId);
 388 }
 389
 390 //=================================================================================
 391 // function : SelectionIntoArgument()
 392 // purpose  : Called when selection as changed or other case
 393 //=================================================================================
 394 void EntityGUI_FeatureDetectorDlg::SelectionIntoArgument()
 395 {
 396
 397   // TODO supprimer les lignes qui ne servent à rien le cas échéant
 398   SUIT_ViewWindow*       theViewWindow  = getDesktop()->activeWindow();
 399   std::map< std::string , std::vector<Handle(AIS_InteractiveObject)> >::iterator AISit;
 400   SOCC_Viewer* soccViewer = (SOCC_Viewer*)(theViewWindow->getViewManager()->getViewModel());
 401
 402   if (!myEditCurrentArgument->isEnabled())
 403     return;
 404
 405   myEditCurrentArgument->setText( "" );
 406
 407   LightApp_SelectionMgr* aSelMgr = myGeomGUI->getApp()->selectionMgr();
 408   SALOME_ListIO aSelList;
 409   aSelMgr->selectedObjects(aSelList);
 410   SALOME_ListIteratorOfListIO anIt( aSelList );
 411
 412   for( ; anIt.More(); anIt.Next() )
 413     if( !anIt.Value().IsNull() )
 414     {
 415       myFaceEntry = anIt.Value()->getEntry();
 416     }
 417
 418   if (aSelList.Extent() != 1) {
 419     if (myEditCurrentArgument == mySelectionGroup->LineEdit1)
 420       myFace.nullify();
 421     return;
 422   }
 423
 424   TopAbs_ShapeEnum aNeedType = TopAbs_FACE ;
 425   GEOM::GeomObjPtr aSelectedObject = getSelected( aNeedType );
 426   TopoDS_Shape aShape;
 427   if ( aSelectedObject && GEOMBase::GetShape( aSelectedObject.get(), aShape ) && !aShape.IsNull() ) {
 428     QString aName = GEOMBase::GetName( aSelectedObject.get() );
 429     myEditCurrentArgument->setText( aName );
 430
 431     if ( myEditCurrentArgument == mySelectionGroup->LineEdit1 ) {
 432       myFace = aSelectedObject;
 433       AISit = soccViewer->entry2aisobjects.find(myFaceEntry.toStdString());
 434       if (AISit == soccViewer->entry2aisobjects.end())
 435         return;
 436
 437       Handle(AIS_InteractiveObject) myAIS = (*AISit).second[0];
 438       Handle(GEOM_AISShape) myAISShape;
 439       if( myAIS->IsInstance( STANDARD_TYPE(GEOM_AISShape) ) ) {
 440         myAISShape = Handle(GEOM_AISShape)::DownCast( myAIS );
 441       }
 442       else
 443         return ;
 444
 445       std::string theImgFileName = myAISShape->TextureFile();
 446       if ( theImgFileName == "" )
 447         return ;
 448
 449       // Setting the image caracteristics
 450       aDetector->SetPath( theImgFileName );
 451       height            =  aDetector->GetImgHeight();
 452       width             =  aDetector->GetImgWidth();
 453       pictureLeft       = -0.5 * width;              // X coordinate of the top left  corner of the background image in the view
 454       pictureTop        =  0.5 * height;             // Y coordinate of both top corners
 455
 456     }
 457   }
 458
 459 }
 460
 461 //=================================================================================
 462 // function : acceptMouseEvent()
 463 // purpose  :
 464 //=================================================================================
 465 bool EntityGUI_FeatureDetectorDlg::acceptMouseEvent() const
 466 {
 467   return mySelectionGroup->PushButton2->isChecked();
 468 }
 469
 470 //=======================================================================
 471 // function : ClickOnOk()
 472 // purpose  :
 473 //=======================================================================
 474 void EntityGUI_FeatureDetectorDlg::ClickOnOk()
 475 {
 476   setIsApplyAndClose( true );
 477   if (onAccept())
 478     ClickOnCancel();
 479 }
 480
 481 //=================================================================================
 482 // function : ClickOnApply()
 483 // purpose  :
 484 //=================================================================================
 485 bool EntityGUI_FeatureDetectorDlg::ClickOnApply()
 486 {
 487   if (!onAccept())
 488     return false;
 489
 490 //   initName();
 491   ConstructorsClicked(getConstructorId());
 492   return true;
 493 }
 494
 495 //=================================================================================
 496 // function : onConstructorClicked(int)
 497 // purpose  :
 498 //=================================================================================
 499 void EntityGUI_FeatureDetectorDlg::ConstructorsClicked(int id)
 500 {
 501   MESSAGE("Constructor id ="<<id)
 502   myConstructorId = id;
 503   switch (id)
 504   {
 505     case CORNERS:
 506       myCornersParameters->show();
 507       myContoursParameters->hide();
 508 //       myViewGroup->show();
 509 //       mySelectionGroup->show();
 510       myOutputGroup->hide();
 511 //       mySelectionGroup->TextLabel2->setText(tr("GEOM_DETECT_ZONE"));
 512       initName(tr("GEOM_CORNERS"));
 513       break;
 514     case CONTOURS:
 515       myCornersParameters->hide();
 516       myContoursParameters->show();
 517 //       myViewGroup->hide();
 518 //       mySelectionGroup->hide();
 519 //       mySelectionGroup->show();
 520       myOutputGroup->show();
 521       mySelectionGroup->TextLabel2->setText(tr("GEOM_COLOR_FILTER"));
 522       initName(tr("GEOM_CONTOURS"));
 523       break;
 524     case LINES:
 525 //       myViewGroup->hide();
 526 //       mySelectionGroup->hide();
 527 //       mySelectionGroup->show();
 528       myOutputGroup->hide();
 529       mySelectionGroup->TextLabel2->setText(tr(""));
 530       initName(tr("GEOM_LINES"));
 531       break;
 532   }
 533 }
 534
 535 // //=================================================================================
 536 // // function : onViewClicked()
 537 // // purpose  :
 538 // //=================================================================================
 539 // void EntityGUI_FeatureDetectorDlg::onViewClicked(int id)
 540 // {
 541 //   gp_Pnt aOrigin = gp_Pnt(0, 0, 0);
 542 //   gp_Dir aDirZ;
 543 //   gp_Dir aDirX;
 544 //
 545 //   switch(id)
 546 //   {
 547 //     case XY:
 548 //         aDirZ = gp_Dir(0, 0, 1);
 549 //         aDirX = gp_Dir(1, 0, 0);
 550 //       break;
 551 //     case YZ:
 552 //         aDirZ = gp_Dir(1, 0, 0);
 553 //         aDirX = gp_Dir(0, 1, 0);
 554 //       break;
 555 //     case XZ:
 556 //         aDirZ = gp_Dir(0, -1, 0);
 557 //         aDirX = gp_Dir(1, 0 , 0);
 558 //       break;
 559 //   }
 560 //
 561 //   myWPlane = gp_Ax3(aOrigin, aDirZ, aDirX);
 562 //   myGeomGUI->SetWorkingPlane( myWPlane );
 563 //   myGeomGUI->ActiveWorkingPlane();
 564 //
 565 // }
 566
 567 //=================================================================================
 568 // function : onButtonClicked()
 569 // purpose  :
 570 //=================================================================================
 571 void EntityGUI_FeatureDetectorDlg::onButtonClicked()
 572 {
 573   QPushButton* send = (QPushButton*)sender();
 574   if (send == mySelectionGroup->PushButton2)
 575   {
 576     mySelectionGroup->LineEdit1->setEnabled(false);
 577   }
 578   else if (send == mySelectionGroup->PushButton1)
 579   {
 580     myStartPnt = gp_Pnt(0,0,0);
 581     myEndPnt   = myStartPnt;
 582     myEditCurrentArgument = mySelectionGroup->LineEdit1;
 583     mySelectionGroup->LineEdit1->setEnabled(true);
 584   }
 585 }
 586
 587 //=================================================================================
 588 // function : onCheckBoxClicked(bool)
 589 // purpose  :
 590 //=================================================================================
 591 void EntityGUI_FeatureDetectorDlg::onCheckBoxClicked( bool isChecked )
 592 {
 593   mySelectionGroup->TextLabel2->setVisible( isChecked );
 594   mySelectionGroup->Frame->setVisible( isChecked );
 595   mySelectionGroup->PushButton2->setVisible( isChecked );
 596   myCannyParameters->setHidden( isChecked );
 597   myColorFilterParameters->setHidden( !isChecked );
 598 }
 599
 600 //=================================================================================
 601 // function : parametersChanged()
 602 // purpose  :
 603 //=================================================================================
 604 ShapeRec_Parameters* EntityGUI_FeatureDetectorDlg::parametersChanged()
 605 {
 606   ShapeRec_Parameters* aParameters = new ShapeRec_Parameters();
 607
 608   if ( myConstructorId == CORNERS ) {
 609     ShapeRec_CornersParameters* aCornersParameters = dynamic_cast<ShapeRec_CornersParameters*>( aParameters );
 610     if ( !aCornersParameters ) aCornersParameters = new ShapeRec_CornersParameters();
 611     aCornersParameters->qualityLevel = (dynamic_cast<QDoubleSpinBox*>(myWidgets[QUALITY_LEVEL]))->value();
 612     aCornersParameters->minDistance  = (dynamic_cast<QDoubleSpinBox*>(myWidgets[MIN_DISTANCE]))->value();
 613     switch ( (dynamic_cast<QComboBox*>(myWidgets[TYPE_CRITERIA]))->currentIndex() ) {
 614     case 0: aCornersParameters->typeCriteria = CV_TERMCRIT_ITER;
 615     case 1: aCornersParameters->typeCriteria = CV_TERMCRIT_EPS;
 616     case 2: aCornersParameters->typeCriteria = CV_TERMCRIT_ITER | CV_TERMCRIT_EPS;
 617     }
 618     aCornersParameters->maxIter = (dynamic_cast<QSpinBox*>(myWidgets[MAX_ITER]))->value();
 619     aCornersParameters->epsilon = (dynamic_cast<QDoubleSpinBox*>(myWidgets[EPSILON]))->value();
 620     aParameters = aCornersParameters;
 621   }
 622   else if ( myConstructorId == CONTOURS ) {
 623     if ( !myUseROI->isChecked() ) {
 624       ShapeRec_CannyParameters* aCannyParameters = dynamic_cast<ShapeRec_CannyParameters*>( aParameters );
 625       if ( !aCannyParameters ) aCannyParameters = new ShapeRec_CannyParameters();
 626       aCannyParameters->lowThreshold = (dynamic_cast<QSpinBox*>(myWidgets[LOW_THRESHOLD]))->value();
 627       aCannyParameters->ratio        = (dynamic_cast<QSpinBox*>(myWidgets[RATIO]))->value();
 628       aCannyParameters->L2gradient   = (dynamic_cast<QCheckBox*>(myWidgets[L2GRADIENT]))->isChecked();
 629       aParameters = aCannyParameters;
 630     }
 631     else {
 632       ShapeRec_ColorFilterParameters* aColorFilterParameters = dynamic_cast<ShapeRec_ColorFilterParameters*>( aParameters );
 633       if ( !aColorFilterParameters ) aColorFilterParameters = new ShapeRec_ColorFilterParameters();
 634       int aSmoothSize = (dynamic_cast<QSpinBox*>(myWidgets[SMOOTH_SIZE]))->value();
 635       aColorFilterParameters->smoothSize  = aSmoothSize % 2 == 0 ? aSmoothSize - 1 : aSmoothSize;
 636       aColorFilterParameters->histSize[0] = (dynamic_cast<QSpinBox*>(myWidgets[HBINS]))->value();
 637       aColorFilterParameters->histSize[1] = (dynamic_cast<QSpinBox*>(myWidgets[SBINS]))->value();
 638       switch ( (dynamic_cast<QComboBox*>(myWidgets[HIST_TYPE]))->currentIndex() ) {
 639       case 0: aColorFilterParameters->histType = CV_HIST_ARRAY;  break;
 640       case 1: aColorFilterParameters->histType = CV_HIST_SPARSE; break;
 641       }
 642       aColorFilterParameters->threshold    = (dynamic_cast<QDoubleSpinBox*>(myWidgets[THRESHOLD_VALUE]))->value();
 643       aColorFilterParameters->maxThreshold = (dynamic_cast<QDoubleSpinBox*>(myWidgets[MAX_THRESHOLD]))->value();
 644       aParameters = aColorFilterParameters;
 645     }
 646   }
 647
 648   aParameters->kernelSize = ( (dynamic_cast<QComboBox*>(myWidgets[KERNEL_SIZE]))->currentText() ).toInt();
 649   switch ( (dynamic_cast<QComboBox*>(myWidgets[FIND_CONTOURS_METHOD]))->currentIndex() ) {
 650   case 0: aParameters->findContoursMethod = CV_CHAIN_APPROX_NONE;      break;
 651   case 1: aParameters->findContoursMethod = CV_CHAIN_APPROX_SIMPLE;    break;
 652   case 2: aParameters->findContoursMethod = CV_CHAIN_APPROX_TC89_L1;   break;
 653   case 3: aParameters->findContoursMethod = CV_CHAIN_APPROX_TC89_KCOS; break;
 654   }
 655
 656   return aParameters;
 657 }
 658
 659 //=================================================================================
 660 // function : setStartPnt( const QPoint& )
 661 // purpose  :
 662 //=================================================================================
 663 void EntityGUI_FeatureDetectorDlg::setStartPnt(const gp_Pnt& theStartPnt)
 664 {
 665   myStartPnt = theStartPnt;
 666   MESSAGE("myStartPnt = ("<<theStartPnt.X()<<", "<<theStartPnt.Y()<<")")
 667 }
 668
 669 //=================================================================================
 670 // function : setEndPnt( const QPoint& )
 671 // purpose  :
 672 //=================================================================================
 673 void EntityGUI_FeatureDetectorDlg::setEndPnt(const gp_Pnt& theEndPnt)
 674 {
 675   myEndPnt = theEndPnt;
 676   MESSAGE("myEndPnt = ("<<theEndPnt.X()<<", "<<theEndPnt.Y()<<")")
 677   if (setSelectionRect() && aDetector->GetImgHeight() > 0)
 678     showImageSample();
 679 }
 680
 681 //=================================================================================
 682 // function : setSelectionRect()
 683 // purpose  :
 684 //=================================================================================
 685 bool EntityGUI_FeatureDetectorDlg::setSelectionRect()
 686 {
 687   // Set detection rectangle in the background image coordinates system
 688   double left    = std::min( myStartPnt.X(), myEndPnt.X() );
 689   double top     = std::max( myStartPnt.Y(), myEndPnt.Y() );
 690   double right   = std::max( myStartPnt.X(), myEndPnt.X() );
 691   double bottom  = std::min( myStartPnt.Y(), myEndPnt.Y() );
 692
 693   QPoint topLeft     = QPoint(left  - pictureLeft, pictureTop - top   );
 694   QPoint bottomRight = QPoint(right - pictureLeft, pictureTop - bottom);
 695
 696   myRect = QRect(topLeft, bottomRight);
 697
 698   return (!myRect.isEmpty() && myRect.width() > 1);
 699 }
 700
 701 //=================================================================================
 702 // function : showImageSample()
 703 // purpose  : Display a preview of the image sample selected by the user
 704 //=================================================================================
 705 void EntityGUI_FeatureDetectorDlg::showImageSample()
 706 {
 707   // Cropp the image to the selection rectangle given by the user
 708   aDetector->SetROI( myRect );
 709   std::string samplePicturePath = aDetector->CroppImage();
 710
 711   // Display the result
 712   QPixmap pixmap(QString(samplePicturePath.c_str()));
 713   mySelectionGroup->FrameLabel->setPixmap(pixmap);
 714   mySelectionGroup->FrameLabel->setMask(pixmap.mask());
 715 }
 716
 717 //=================================================================================
 718 // function : createOperation
 719 // purpose  :
 720 //=================================================================================
 721 GEOM::GEOM_IOperations_ptr EntityGUI_FeatureDetectorDlg::createOperation()
 722 {
 723   return myGeomGUI->GetGeomGen()->GetIShapesOperations( getStudyId() );
 724 }
 725
 726 //=================================================================================
 727 // function : execute()
 728 // purpose  :
 729 //=================================================================================
 730 bool EntityGUI_FeatureDetectorDlg::execute( ObjectList& objects )
 731 {
 732   bool res = false;
 733
 734   GEOM::GEOM_IBasicOperations_var  aBasicOperations  = myGeomGUI->GetGeomGen()->GetIBasicOperations( getStudyId() );
 735   GEOM::GEOM_IShapesOperations_var aShapesOperations = GEOM::GEOM_IShapesOperations::_narrow( getOperation() );
 736
 737   ShapeRec_Parameters* parameters = parametersChanged();
 738   if ( !parameters ) parameters = new ShapeRec_Parameters();
 739
 740   bool useROI = myUseROI->isChecked() && !myRect.isEmpty();
 741   try {
 742     if (myConstructorId == CORNERS)
 743     {
 744       double subPictureLeft;
 745       double subPictureTop;
 746       if( useROI )
 747       {
 748         subPictureLeft    = myStartPnt.X();
 749         subPictureTop     = myStartPnt.Y();
 750       }
 751       else
 752       {
 753         subPictureLeft    = pictureLeft;
 754         subPictureTop     = pictureTop;
 755       }
 756       aDetector->ComputeCorners( useROI, parameters );
 757       CvPoint2D32f* corners     = aDetector->GetCorners();
 758       int cornerCount           = aDetector->GetCornerCount();
 759       int i;
 760
 761       // Build the geom objects associated to the detected corners and returned by execute
 762       if( !aBasicOperations->_is_nil() && !aShapesOperations->_is_nil() )
 763       {
 764         GEOM::GEOM_Object_var  aGeomCorner;
 765         GEOM::ListOfGO_var     geomCorners = new GEOM::ListOfGO();
 766         geomCorners->length( cornerCount );
 767         for (i = 0; i < cornerCount; i++)
 768         {
 769           double x = subPictureLeft + corners[i].x;
 770           double y = subPictureTop  - corners[i].y;
 771           double z =  0;
 772
 773           aGeomCorner = aBasicOperations->MakePointXYZ( x,y,z );
 774
 775           geomCorners[i] = aGeomCorner;
 776         }
 777         GEOM::GEOM_Object_var aCompound = aShapesOperations->MakeCompound(geomCorners);
 778         if ( !aCompound->_is_nil() )
 779         {
 780           // Dark blue color
 781           SALOMEDS::Color aColor;
 782           aColor.R = 0;
 783           aColor.G = 0;
 784           aColor.B = 0.8;
 785
 786           aCompound->SetColor(aColor);
 787           aCompound->SetMarkerStd(GEOM::MT_POINT,GEOM::MS_30);
 788           objects.push_back( aCompound._retn() );
 789           res = true;
 790         }
 791       }
 792     }
 793     else if (myConstructorId == CONTOURS)
 794     {
 795       GEOM::GEOM_ICurvesOperations_var aCurveOperations = myGeomGUI->GetGeomGen()->GetICurvesOperations( getStudyId() );
 796
 797       aDetector->ComputeContours( useROI, parameters );
 798       std::vector< std::vector<cv::Point> >   contours  = aDetector->GetContours();
 799       std::vector<cv::Vec4i>                  hierarchy = aDetector->GetContoursHierarchy();
 800
 801       std::vector< cv::Point >                contour;
 802       int idx = 0;
 803
 804       GEOM::ListOfGO_var                      geomContours = new GEOM::ListOfGO();
 805       int contourCount = 0;
 806
 807       bool insert;
 808
 809       MESSAGE("hierarchy.size() =" << hierarchy.size())
 810       if ( hierarchy.size() < 1 ) {
 811         getOperation()->SetErrorCode( "Impossible detected contours" );
 812         return false;
 813       }
 814
 815       for( ; idx >= 0; idx = hierarchy[idx][0])
 816       {
 817 //       for(int count=0, child=idx; child>=0, count<1; child=hierarchy[idx][2], count++)
 818 //       {
 819 //         contour = contours[child];
 820         contour = contours[idx];
 821         std::vector< cv::Point >::iterator it;
 822         std::vector< cv::Point >::iterator it_previous;
 823         std::vector< cv::Point >::iterator it_next;
 824         GEOM::GEOM_Object_var  aGeomContourPnt;
 825         GEOM::ListOfGO_var     geomContourPnts = new GEOM::ListOfGO();
 826
 827         geomContourPnts->length( contour.size() );
 828
 829         int j = 0;
 830         std::set< std::vector<int> > existing_points;
 831         std::pair< std::set< std::vector<int> >::iterator,bool > pnt_it;
 832         for ( it=contour.begin() ; it < contour.end()-1; it++ )
 833         {
 834           int pnt_array[] = {it->x,it->y};
 835           std::vector<int> pnt (pnt_array, pnt_array + sizeof(pnt_array) / sizeof(int) );
 836
 837           pnt_it=existing_points.insert(pnt);
 838           if (pnt_it.second == true)         // To avoid double points in the contours
 839           {
 840             insert = true;
 841             if (it!=contour.begin())         // From the second point on perform some checking to avoid loops in the contours we build
 842             {
 843               it_previous = boost::prior(it);
 844               it_next = boost::next(it);
 845
 846               double u_v_scalar_product = (it->x - it_previous->x) * (it_next->x - it->x) +
 847                                           (it->y - it_previous->y) * (it_next->y - it->y);
 848               if (u_v_scalar_product < 0)
 849               {
 850                 double u_v_det = (it->x - it_previous->x) * (it_next->y - it->y) -
 851                                  (it->y - it_previous->y) * (it_next->x - it->x);
 852
 853                 double norme_u = sqrt ( (it->x - it_previous->x)*(it->x - it_previous->x) +
 854                                         (it->y - it_previous->y)*(it->y - it_previous->y) );
 855
 856                 double norme_v = sqrt ( (it->x - it_next->x)*(it->x - it_next->x) +
 857                                         (it->y - it_next->y)*(it->y - it_next->y) );
 858
 859                 double u_v_sinus = u_v_det / (norme_u * norme_v);
 860
 861                 if (fabs(u_v_sinus) < Precision::Confusion())
 862                 {
 863                   insert = false;
 864                 }
 865               }
 866             }
 867             double x = -0.5 *width  + it->x;
 868             double y =  0.5 *height - it->y;
 869             double z =  0;
 870             if (insert)
 871             {
 872               aGeomContourPnt    = aBasicOperations->MakePointXYZ( x,y,z );
 873               geomContourPnts->length( j+1 );
 874               geomContourPnts[j] = aGeomContourPnt;
 875               j++;
 876             }
 877           }
 878         }
 879
 880         GEOM::GEOM_Object_var aWire;
 881         if(myOutputGroup->RadioButton2->isChecked())
 882         {
 883           aWire = aCurveOperations->MakePolyline(geomContourPnts.in(), false);
 884         }
 885         else if(myOutputGroup->RadioButton1->isChecked())
 886         {
 887           aWire = aCurveOperations->MakeSplineInterpolation(geomContourPnts.in(), /*closed =*/ false, /*reordering =*/ false);
 888         }
 889         else
 890           return res;
 891
 892         if ( !aWire->_is_nil() )
 893         {
 894           geomContours->length(contourCount + 1);
 895           geomContours[contourCount] = aWire;
 896           contourCount++;
 897         }
 898 //       }
 899       }
 900       GEOM::GEOM_Object_var aContoursCompound = aShapesOperations->MakeCompound(geomContours);
 901       if ( !aContoursCompound->_is_nil() )
 902       {
 903         objects.push_back( aContoursCompound._retn() );
 904       }
 905       res=true;
 906     }
 907   }
 908   catch( ... ) {
 909     getOperation()->SetErrorCode( "Impossible detected contours/corners" );
 910     return false;
 911   }
 912
 913   // TEST not very conclusive
 914
 915 //   else if(myConstructorId ==LINES)
 916 //   {
 917 //     aDetector->ComputeLines();
 918 //     std::vector<cv::Vec4i>  lines = aDetector->GetLines();
 919 //     GEOM::GEOM_Object_var  Pnt1;
 920 //     GEOM::GEOM_Object_var  Pnt2;
 921 //     GEOM::GEOM_Object_var  aLine;
 922 //
 923 //     GEOM::ListOfGO_var     geomLines = new GEOM::ListOfGO();
 924 //     int linesCount=0;
 925 //     for( int i = 0; i < lines.size(); i++ )
 926 //     {
 927 //       Pnt1 = aBasicOperations->MakePointXYZ( -0.5 *width + lines[i][0], 0.5 *height - lines[i][1], 0 );
 928 //       Pnt2 = aBasicOperations->MakePointXYZ( -0.5 *width + lines[i][2], 0.5 *height - lines[i][3], 0 );
 929 //       aLine = aBasicOperations->MakeLineTwoPnt( Pnt1, Pnt2 );
 930 //       if ( !aLine->_is_nil() )
 931 //       {
 932 //         geomLines->length(linesCount + 1);
 933 //         geomLines[linesCount] = aLine;
 934 //         linesCount++;
 935 //       }
 936 //     }
 937 //     GEOM::GEOM_Object_var aLinesCompound = aShapesOperations->MakeCompound(geomLines);
 938 //     if ( !aLinesCompound->_is_nil() )
 939 //     {
 940 //       objects.push_back( aLinesCompound._retn() );
 941 //     }
 942 //
 943 //     res=true;
 944 //   }
 945
 946   return res;
 947 }