From: Eduardo DAVILA Date: Tue, 20 Jun 2023 08:35:41 +0000 (+0200) Subject: #3505 bug Deformation action region X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=commitdiff_plain;h=54a9fee9effafd5d888252e6d622cdea2e2d76ed;p=bbtk.git #3505 bug Deformation action region --- diff --git a/packages/itkvtk/src/bbitkvtkGeodesicMeshDeformation.cxx b/packages/itkvtk/src/bbitkvtkGeodesicMeshDeformation.cxx index 9fb4792..1ecc52f 100644 --- a/packages/itkvtk/src/bbitkvtkGeodesicMeshDeformation.cxx +++ b/packages/itkvtk/src/bbitkvtkGeodesicMeshDeformation.cxx @@ -26,10 +26,10 @@ void GeodesicMeshDeformation::Process() // * TYPE is the C++ type of the input/output // (the one provided in the attribute 'type' of the tag 'input') std::vector lstCenter = bbGetInputCenter(); - double s = bbGetInputS(); - bool ok = false; - bool pdChanged = false; - using MeshType = itk::QuadEdgeMesh; + double s = bbGetInputS(); + bool ok = false; + bool pdChanged = false; + using MeshType = itk::QuadEdgeMesh; //Set up QuadEdge and filter every time polydata changes if ((bbGetInputIn() != polydata) && (bbGetInputActive()==true) && (bbGetInputIn() != NULL)) @@ -124,10 +124,10 @@ void GeodesicMeshDeformation::Process() double p[3]; // point double pb[3]; // point base double np[3]; // new point - double sx,sy,sz; - sx = s*2; - sy = sx; - sz = sy; +// double sx,sy,sz; +// sx = s*2; +// sy = s*2; +// sz = s*2; double displcement_x = 0; double displcement_y = 0; @@ -184,7 +184,7 @@ void GeodesicMeshDeformation::Process() using CriterionType = itk::FastMarchingThresholdStoppingCriterion; auto criterion = CriterionType::New(); sCurrent = s; - criterion->SetThreshold(s*2); + criterion->SetThreshold(s*4); fmmFilter->SetStoppingCriterion(criterion); } fmmFilter->Update();