From 48b38447493c9ca87a9d3bf7cd0db35cb82657b2 Mon Sep 17 00:00:00 2001
From: Julien Michel <julien.michel@c-s.fr>
Date: Mon, 21 Jan 2008 14:14:11 +0000
Subject: [PATCH] test
---
Code/BasicFilters/otbProlateInterpolateImageFunction.h | 8 +++++---
1 file changed, 5 insertions(+), 3 deletions(-)
diff --git a/Code/BasicFilters/otbProlateInterpolateImageFunction.h b/Code/BasicFilters/otbProlateInterpolateImageFunction.h
index 4fc23a689d..4e95649060 100644
--- a/Code/BasicFilters/otbProlateInterpolateImageFunction.h
+++ b/Code/BasicFilters/otbProlateInterpolateImageFunction.h
@@ -47,10 +47,12 @@ class ProlateFunction
inline TOutput operator()( const TInput & A ) const
{
TOutput val = itk::NumericTraits< TOutput >::Zero;
- if (A != itk::NumericTraits< TInput >::Zero)
+ if (A != itk::NumericTraits< TInput >::Zero && m_Radius!=0)
{
- unsigned int ival = static_cast<unsigned int>(m_OriginalProfilSize*static_cast<double>(vcl_abs(A))/static_cast<double>(m_Radius+1));
- val = m_OriginalProfil[ival];
+ //unsigned int ival = static_cast<unsigned int>(m_OriginalProfilSize*static_cast<double>(vcl_abs(A))/static_cast<double>(m_Radius));
+ //val = m_OriginalProfil[ival];
+ double ival = static_cast<double>(m_OriginalProfilSize)*static_cast<double>(vcl_abs(A))/static_cast<double>(m_Radius);
+ val = m_OriginalProfil[static_cast<int>(vcl_floor(ival+0.5))];
}
else
{
--
GitLab