From: Pablo Garzon <gapablo2001@gmail.com>
Date: Fri, 28 Jul 2023 09:18:14 +0000 (+0200)
Subject: #3501 Geodesic Deformation
X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=commitdiff_plain;h=296a644550aec58bcab0a7c55e486da19eef78f1;p=bbtk.git

#3501 Geodesic Deformation
---

diff --git a/packages/itkvtk/src/bbitkvtkGeodesicMeshDeformation.cxx b/packages/itkvtk/src/bbitkvtkGeodesicMeshDeformation.cxx
index 9bcf5fe..9855950 100644
--- a/packages/itkvtk/src/bbitkvtkGeodesicMeshDeformation.cxx
+++ b/packages/itkvtk/src/bbitkvtkGeodesicMeshDeformation.cxx
@@ -32,6 +32,7 @@ void GeodesicMeshDeformation::Process()
     using MeshType      = itk::QuadEdgeMesh<double, 3>;
     std::vector<double> deformInfo;
 	bbSetOutputOut(deformInfo);
+	//std::vector<double> displacementVector;
     
     //Set up QuadEdge and filter every time polydata changes
     if ((bbGetInputIn() != polydata) && (bbGetInputActive()==true) && (bbGetInputIn() != NULL))
@@ -215,8 +216,10 @@ void GeodesicMeshDeformation::Process()
 						points->SetPoint(BegProcessedIt.Value().GetNode(), np);
 					}
 					++BegProcessedIt;
-				}				
-				std::vector<double> info{pb[0],pb[1],pb[2], (double) bbGetInputEdgeId(), s};
+				}
+				//double directionMoved[3] = {lstCenter[0]-displcement_x, lstCenter[1]-displcement_y, lstCenter[2]-displcement_z};
+				//vtkMath::Normalize(directionMoved);
+				std::vector<double> info{lstCenter[0],lstCenter[1],lstCenter[2], bbGetInputDirection()[0] , bbGetInputDirection()[1], bbGetInputDirection()[2],(double) bbGetInputEdgeId(), s};
 				bbSetOutputOut(info);
 				cout << "info updated" << endl;
 				points->Modified();