X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=itk%2FitkBSplineInterpolateImageFunctionWithLUT.txx;h=b9354a5e13ce2d13762de9a7f5eefef4f7d410e8;hb=573d80d0f7a17607d2ee883c21c940c0ba020282;hp=c776afae94d3f338ef752443c0c06bf575250fd5;hpb=a87fe1bfb8c6b85456952f0b584d44b1da6cb3fa;p=clitk.git diff --git a/itk/itkBSplineInterpolateImageFunctionWithLUT.txx b/itk/itkBSplineInterpolateImageFunctionWithLUT.txx index c776afa..b9354a5 100644 --- a/itk/itkBSplineInterpolateImageFunctionWithLUT.txx +++ b/itk/itkBSplineInterpolateImageFunctionWithLUT.txx @@ -1,23 +1,24 @@ +/*========================================================================= + 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://www.centreleonberard.fr + - CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr + + This software is distributed WITHOUT ANY WARRANTY; without even + the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR + PURPOSE. See the copyright notices for more information. + + It is distributed under dual licence + + - BSD See included LICENSE.txt file + - CeCILL-B http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html +===========================================================================**/ #ifndef _ITKINTERPOLATEIMAGEFUNCTIONWITHLUT_TXX #define _ITKINTERPOLATEIMAGEFUNCTIONWITHLUT_TXX +#include -/* ========================================================================= - -@file itkBSplineInterpolateImageFunctionWithLUT.txx -@author David Sarrut - -Copyright (c) -* CREATIS (Centre de Recherche et d'Applications en Traitement de l'Image). -All rights reserved. See Doc/License.txt or -http://www.creatis.insa-lyon.fr/Public/Gdcm/License.html for details. -* Léon Bérard cancer center, 28 rue Laënnec, 69373 Lyon cedex 08, France -* http://www.creatis.insa-lyon.fr/rio - -This software is distributed WITHOUT ANY WARRANTY; without even the -implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR -PURPOSE. See the above copyright notices for more information. - -========================================================================= */ namespace itk { //==================================================================== @@ -190,7 +191,7 @@ namespace itk // Compute index in precomputed weights table TCoefficientType t2 = mSamplingFactors[l]*t1; - index[l] = (IndexValueType)lrint(t2); + index[l] = (IndexValueType)Math::Round(t2); // For even order : test if too close to 0.5 (but lower). In this // case : take the next coefficient @@ -240,7 +241,7 @@ namespace itk evaluateIndex[l] = indx; } else { // Use this index calculation for even splineOrder - if (mSplineOrders[l] == 0) evaluateIndex[l] = (long)rint(x[l]); + if (mSplineOrders[l] == 0) evaluateIndex[l] = Math::Round(x[l]); else { indx = (long)vcl_floor((x[l]+ 0.5)) - mSplineOrders[l] / 2; //this->m_SplineOrder / 2; evaluateIndex[l] = indx;