+// bbSetOutputOut( append->GetOutput() );
+
+ // Step 4. -- Define scalars ---
+ vtkCharArray *scalarsArray = vtkCharArray::New();
+ vtkPoints *points = bbGetInputMesh()->GetPoints();
+ long int iMPs,sizeMPs = points->GetNumberOfPoints();
+ scalarsArray->SetNumberOfValues(sizeMPs);
+ for (iMPs=0; iMPs<sizeMPs; iMPs++)
+ {
+ scalarsArray->SetValue(iMPs,15);
+ } // for iMPs
+ long int iIdLst,sizeIdLst = lstIdsGeodesicPerimter.size();
+ for (iIdLst=0; iIdLst<sizeIdLst; iIdLst++)
+ {
+ scalarsArray->SetValue( lstIdsGeodesicPerimter[iIdLst] ,100 );
+ }// for iIdLst
+ scalarsArray->SetName("scalarsPerimeterPatch");
+ bbGetInputMesh()->GetPointData()->AddArray( scalarsArray );
+ bbGetInputMesh()->GetPointData()->SetActiveScalars( "scalarsPerimeterPatch" );
+ bbGetInputMesh()->GetCellData()->SetActiveScalars( "scalarsPerimeterPatch" );
+
+ // Step 5. -- Connectivity Filter
+ vtkPolyDataConnectivityFilter *connectivity = vtkPolyDataConnectivityFilter::New();
+ connectivity->SetInputData( bbGetInputMesh() );
+ connectivity->SetExtractionModeToCellSeededRegions();
+ connectivity->InitializeSeedList();
+ connectivity->AddSeed(0);
+ connectivity->SetScalarConnectivity(true);
+ connectivity->SetScalarRange(10,20);
+ connectivity->SetFullScalarConnectivity(true);
+ connectivity->Update();
+// bbSetOutputOut( connectivity->GetOutput() );