Second version of RelativePositionAnalyzer

[clitk.git] / registration / clitkBSplineDeformableTransformInitializer.txx

diff --git a/registration/clitkBSplineDeformableTransformInitializer.txx b/registration/clitkBSplineDeformableTransformInitializer.txx
index 0a93ab0432ae814347454e5b47121e016e8bbedd..59af63c7d211c98d4b4e8a40c29a1b5de3a91f4e 100644 (file)
--- a/registration/clitkBSplineDeformableTransformInitializer.txx
+++ b/registration/clitkBSplineDeformableTransformInitializer.txx
@@ -1,9 +1,9 @@
-/*=========================================================================
+/*=========================================================================
   Program:   vv                     http://www.creatis.insa-lyon.fr/rio/vv
 
   Authors belong to: 
   - University of LYON              http://www.universite-lyon.fr/
-  - Léon Bérard cancer center       http://oncora1.lyon.fnclcc.fr
+  - Léon Bérard cancer center       http://www.centreleonberard.fr
   - CREATIS CNRS laboratory         http://www.creatis.insa-lyon.fr
 
   This software is distributed WITHOUT ANY WARRANTY; without even
@@ -14,10 +14,11 @@
 
   - BSD        See included LICENSE.txt file
   - CeCILL-B   http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
-======================================================================-====*/
+===========================================================================**/
 #ifndef __clitkBSplineDeformableTransformInitializer_txx
 #define __clitkBSplineDeformableTransformInitializer_txx
 #include "clitkBSplineDeformableTransformInitializer.h"
+#include "itkMath.h"
 
 namespace clitk
 {
@@ -80,7 +81,7 @@ namespace clitk
          {
            // JV
            m_ChosenSpacing[r]= m_ControlPointSpacing[r];
-           m_ControlPointSpacing[r]= ( round(m_ChosenSpacing[r]/fixedImageSpacing[r]) *fixedImageSpacing[r] ) ;
+           m_ControlPointSpacing[r]= ( itk::Math::Round<double>(m_ChosenSpacing[r]/fixedImageSpacing[r]) *fixedImageSpacing[r] ) ;
            m_NumberOfControlPointsInsideTheImage[r] = ceil( (double)fixedImageSize[r]*fixedImageSpacing[r]/ m_ControlPointSpacing[r] );
            if (  ( ceil( (double)fixedImageSize[r]*fixedImageSpacing[r]/ m_ControlPointSpacing[r] ) )
                 == ( (double)fixedImageSize[r]*fixedImageSpacing[r]/ m_ControlPointSpacing[r] ) )