mode T1&T2 avec la grille

diff --git a/src/RosetteCIH/plugin/RosetteCIHFilters/vtkRosetteCIH.cxx b/src/RosetteCIH/plugin/RosetteCIHFilters/vtkRosetteCIH.cxx
index 8f38c11a4a7c7d8cf110d604e52c3e9ab839f32a..6a2b19144db7da1bc24addf2848a2ece352acac1 100644 (file)
--- a/src/RosetteCIH/plugin/RosetteCIHFilters/vtkRosetteCIH.cxx
+++ b/src/RosetteCIH/plugin/RosetteCIHFilters/vtkRosetteCIH.cxx
@@ -176,8 +176,23 @@ void vtkRosetteCIH::PostTraitementT1etT2(
   vtkSmartPointer<vtkDataSet> gl2 =
     this->GenerateGlyphLinesFor(usgIn, "T2", COMPRESS_TRACTION);
   //
+  vtkNew<vtkDataSetSurfaceFilter> surface;
+  surface->SetNonlinearSubdivisionLevel(0);
+  surface->SetInputData(usgIn);
+  surface->Update();
+  vtkNew<vtkPolyData> surfaceCpy;
+  surfaceCpy->ShallowCopy(surface->GetOutput());
+  vtkNew<vtkDoubleArray> compressionOrTraction;
+  auto nbOfTuples(surface->GetOutput()->GetNumberOfPoints());
+  compressionOrTraction->SetNumberOfComponents(1);
+  compressionOrTraction->SetNumberOfTuples(nbOfTuples);
+  compressionOrTraction->SetName(COMPRESS_TRACTION);
+  compressionOrTraction->Fill(NAN);
+  surfaceCpy->GetPointData()->AddArray(compressionOrTraction);
+  //
   vtkNew<vtkMultiBlockDataGroupFilter> mb;
   vtkNew<vtkCompositeDataToUnstructuredGridFilter> cd;
+  mb->AddInputData(surfaceCpy);
   mb->AddInputData(gl1);
   mb->AddInputData(gl2);
   cd->SetInputConnection(mb->GetOutputPort());