#3513 CleanMeshWithPatch

[creaVtk.git] / bbtk_creaVtk_PKG / src / bbcreaVtkPolyDataNormals.cxx

diff --git a/bbtk_creaVtk_PKG/src/bbcreaVtkPolyDataNormals.cxx b/bbtk_creaVtk_PKG/src/bbcreaVtkPolyDataNormals.cxx
index 92cf8834597d385f009340c36c9211d871a53d80..46bc24a56656d27a216545591e1f3005b3df1908 100644 (file)
--- a/bbtk_creaVtk_PKG/src/bbcreaVtkPolyDataNormals.cxx
+++ b/bbtk_creaVtk_PKG/src/bbcreaVtkPolyDataNormals.cxx
@@ -5,6 +5,7 @@
 #include "bbcreaVtkPackage.h"
 
 #include <vtkPolyDataNormals.h>
+//#include <vtkTriangleMeshPointNormals.h>
 
 #include <vtkPolyData.h>
 #include <vtkPointData.h>
@@ -21,7 +22,6 @@ BBTK_BLACK_BOX_IMPLEMENTATION(PolyDataNormals,bbtk::AtomicBlackBox);
 //===== 
 void PolyDataNormals::Process()
 {
-    printf("EED PolyDataNormals::Process Start\n");
 // THE MAIN PROCESSING METHOD BODY
 //   Here we simply set the input 'In' value to the output 'Out'
 //   And print out the output value
@@ -35,9 +35,17 @@ void PolyDataNormals::Process()
 //      (the one provided in the attribute 'type' of the tag 'input')
 //    bbSetOutputOut( bbGetInputIn() );
 //    std::cout << "Output value = " <<bbGetOutputOut() << std::endl;
+    
     if (bbGetInputIn()!=NULL)
     {
         vtkPolyDataNormals* normal = vtkPolyDataNormals::New();
+//      Automatically change the orientation fo normals, put it from inside to outside
+//        normal->SetAutoOrientNormals(true);
+        normal->SetConsistency(true);
+//      this force to not change the order of triangles, keep the original order        
+//        normal->SetAutoOrientNormals(false);
+//        normal->SetConsistency(false);
+
         normal->SetInputData( bbGetInputIn() );
         if (bbGetInputComputeType()==0)
         {
@@ -45,8 +53,14 @@ void PolyDataNormals::Process()
         } else {
             normal->ComputePointNormalsOn();
         }
-        //        normal->SetFeatureAngle(60.0);
+        normal->SetFeatureAngle(360);
+        normal->Update();
+    
+        /*
+        vtkTriangleMeshPointNormals* normal = vtkTriangleMeshPointNormals::New();
+        normal->SetInputData( bbGetInputIn() );
         normal->Update();
+        */
         
         // >>>>>>>>>>>>>>>>>>>> Calculate the MeanNormal
         vtkPolyData     *polydata   = normal->GetOutput();
@@ -63,12 +77,11 @@ void PolyDataNormals::Process()
         */
         dataarray   = pointdata->GetNormals();
         std::vector<double> meanNormal;
-        meanNormal.push_back(1);
+        meanNormal.push_back(0);
         meanNormal.push_back(0);
         meanNormal.push_back(0);
         if (dataarray!=NULL)
         {
-            meanNormal[0]=0;
             int i , size=dataarray->GetNumberOfTuples();
             for (i=0; i<size; i++)
             {
@@ -76,18 +89,22 @@ void PolyDataNormals::Process()
                 meanNormal[0] = meanNormal[0] + pValue[0];
                 meanNormal[1] = meanNormal[1] + pValue[1];
                 meanNormal[2] = meanNormal[2] + pValue[2];
-            }
+            } // for i
             
             if (size!=0)
             {
                 meanNormal[0] = meanNormal[0] / size;
                 meanNormal[1] = meanNormal[1] / size;
                 meanNormal[2] = meanNormal[2] / size;
+                double magnitude = sqrt( meanNormal[0]*meanNormal[0] + meanNormal[1]*meanNormal[1]  + meanNormal[2]*meanNormal[2] );
+                meanNormal[0] = meanNormal[0] / magnitude;
+                meanNormal[1] = meanNormal[1] / magnitude;
+                meanNormal[2] = meanNormal[2] / magnitude;
             } else {
                 meanNormal[0] = 1;
                 meanNormal[1] = 0;
                 meanNormal[2] = 0;
-            }
+            } // if size
         } // if dataarray
         // <<<<<<<<<<<<<<<<<<<<< Calculate the MeanNormal
 
@@ -98,37 +115,35 @@ void PolyDataNormals::Process()
             double meanNormalresult[3];
             double ang;
             double crossResult[3];
-            printf("EED PolyDataNormals::Process 1\n");
             meanN[0] = meanNormal[0];
             meanN[1] = meanNormal[1];
             meanN[2] = meanNormal[2];
-            printf("EED PolyDataNormals::Process 2\n");
+            double magnitude;
             int i , size=dataarray->GetNumberOfTuples();
             for (i=0; i<size; i++)
             {
-                printf("EED PolyDataNormals::Process 3\n");
                 pValue  = dataarray->GetTuple3( i );
                 ang     = vtkMath::AngleBetweenVectors( pValue, meanN );
-                
                 vtkMath::Cross(meanN,pValue,crossResult);
                 wxyz[0] = ang * bbGetInputParam();
                 wxyz[1] = crossResult[0];
                 wxyz[2] = crossResult[1];
                 wxyz[3] = crossResult[2];
-                printf("EED PolyDataNormals::Process 4\n");
                 vtkMath::RotateVectorByWXYZ(meanN, wxyz ,meanNormalresult);
+                magnitude = sqrt(   meanNormalresult[0]*meanNormalresult[0] +  meanNormalresult[1]*meanNormalresult[1] + meanNormalresult[2]* meanNormalresult[2]  );
+                meanNormalresult[0] = meanNormalresult[0] / magnitude;
+                meanNormalresult[1] = meanNormalresult[1] / magnitude;
+                meanNormalresult[2] = meanNormalresult[2] / magnitude;
                 dataarray->SetTuple3( i , meanNormalresult[0], meanNormalresult[1], meanNormalresult[2] );
             } // for
-
         } // if Type==1  &&  dataarray
         bbSetOutputMeanNormal( meanNormal );
         bbSetOutputOut( normal->GetOutput() );
-    } // if In
-    
-    printf("EED PolyDataNormals::Process End\n");
-
-    
+    } else {
+        bbSetOutputOut( NULL );
+    }// if In
 }
+
 //===== 
 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
 //===== 
@@ -136,34 +151,34 @@ void PolyDataNormals::bbUserSetDefaultValues()
 {
 //  SET HERE THE DEFAULT INPUT/OUTPUT VALUES OF YOUR BOX 
 //    Here we initialize the input 'In' to 0
-   bbSetInputIn(NULL);
+    bbSetInputType(0);
+    bbSetInputComputeType(0);
+    bbSetInputIn(NULL);
+    bbSetOutputOut( NULL );
 }
+
 //===== 
 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
 //===== 
 void PolyDataNormals::bbUserInitializeProcessing()
 {
-
 //  THE INITIALIZATION METHOD BODY :
 //    Here does nothing 
 //    but this is where you should allocate the internal/output pointers 
-//    if any 
-
-  
+//    if any
 }
+
 //===== 
 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
 //===== 
 void PolyDataNormals::bbUserFinalizeProcessing()
 {
-
 //  THE FINALIZATION METHOD BODY :
 //    Here does nothing 
 //    but this is where you should desallocate the internal/output pointers 
 //    if any
-  
 }
-}
-// EO namespace bbcreaVtk
+
+}// EO namespace bbcreaVtk