src/PIPELINE/VISU_MaskPointsFilter.cxx

   1 // Copyright (C) 2007-2011  CEA/DEN, EDF R&D, OPEN CASCADE
   2 //
   3 // This library is free software; you can redistribute it and/or
   4 // modify it under the terms of the GNU Lesser General Public
   5 // License as published by the Free Software Foundation; either
   6 // version 2.1 of the License.
   7 //
   8 // This library is distributed in the hope that it will be useful,
   9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 // Lesser General Public License for more details.
  12 //
  13 // You should have received a copy of the GNU Lesser General Public
  14 // License along with this library; if not, write to the Free Software
  15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 //
  17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 //
  19
  20 //  VISU OBJECT : interactive object for VISU entities implementation
  21 // File:    VISU_StreamLinesPL.cxx
  22 // Author:  Alexey PETROV
  23 // Module : VISU
  24 //
  25 #include "VISU_MaskPointsFilter.hxx"
  26
  27 #include <vtkObjectFactory.h>
  28 #include <vtkPointSet.h>
  29 #include <vtkPointData.h>
  30 #include <vtkCellData.h>
  31 #include <vtkPoints.h>
  32 #include <vtkIdList.h>
  33
  34 vtkStandardNewMacro(VISU_MaskPointsFilter);
  35
  36 VISU_MaskPointsFilter::VISU_MaskPointsFilter(){
  37   PercentsOfUsedPoints = 1.0;
  38 }
  39
  40 VISU_MaskPointsFilter::~VISU_MaskPointsFilter(){}
  41
  42 void VISU_MaskPointsFilter::Execute(){
  43   vtkPointSet *anInput = this->GetInput(), *anOutput = this->GetOutput();
  44   anOutput->GetPointData()->CopyAllOff();
  45   anOutput->GetCellData()->CopyAllOff();
  46   anOutput->CopyStructure(anInput);
  47
  48   vtkPoints* aPoints = vtkPoints::New();
  49   vtkIdList *anIdList = vtkIdList::New();
  50   vtkIdType iEnd = anInput->GetNumberOfPoints();
  51   for(vtkIdType i = 0; i < iEnd; i++){
  52     anInput->GetPointCells(i,anIdList);
  53     if(anIdList->GetNumberOfIds() > 0)
  54       aPoints->InsertNextPoint(anInput->GetPoint(i));
  55   }
  56   vtkPoints* aNewPoints = vtkPoints::New();
  57   iEnd = aPoints->GetNumberOfPoints();
  58   if (PercentsOfUsedPoints > 0){
  59     vtkIdType anOffset = vtkIdType(1.0/PercentsOfUsedPoints);
  60     if(anOffset < 1) anOffset = 1;
  61     for(vtkIdType i = 0; i < iEnd; i += anOffset)
  62       aNewPoints->InsertNextPoint(aPoints->GetPoint(i));
  63   }
  64   anOutput->SetPoints(aNewPoints);
  65   aNewPoints->Delete();
  66   aPoints->Delete();
  67 }