From 37b27427f2d349af3e9e59d0d12c349f5e2ddb3a Mon Sep 17 00:00:00 2001
From: Guillaume Pasero <guillaume.pasero@c-s.fr>
Date: Wed, 6 Aug 2014 15:55:55 +0200
Subject: [PATCH] BUG: compute gradient in index space since resolution is not
 read from image

---
 Code/FeatureExtraction/otbGenericRoadExtractionFilter.h   | 4 +++-
 Code/FeatureExtraction/otbGenericRoadExtractionFilter.txx | 1 +
 2 files changed, 4 insertions(+), 1 deletion(-)

diff --git a/Code/FeatureExtraction/otbGenericRoadExtractionFilter.h b/Code/FeatureExtraction/otbGenericRoadExtractionFilter.h
index e5778c4941..ff2eda4a91 100644
--- a/Code/FeatureExtraction/otbGenericRoadExtractionFilter.h
+++ b/Code/FeatureExtraction/otbGenericRoadExtractionFilter.h
@@ -187,7 +187,9 @@ protected:
   /** Destructor */
   ~GenericRoadExtractionFilter() {}
 
-  /** Prepare main computation method */
+  /** Prepare main computation method
+   *  Note : this function isn't called
+   */
   void BeforeGenerateData(void);
 
   /** Main computation method */
diff --git a/Code/FeatureExtraction/otbGenericRoadExtractionFilter.txx b/Code/FeatureExtraction/otbGenericRoadExtractionFilter.txx
index 9671526198..6ed3e03127 100644
--- a/Code/FeatureExtraction/otbGenericRoadExtractionFilter.txx
+++ b/Code/FeatureExtraction/otbGenericRoadExtractionFilter.txx
@@ -113,6 +113,7 @@ GenericRoadExtractionFilter<TInputImage, TOutputPath>
   m_GradientFilter->SetInput(m_SquareRootImageFilter->GetOutput());
   /** Sigma calculated with the alpha and image resolution parameters */
   m_GradientFilter->SetSigma(static_cast<SigmaType>(m_Alpha * (1.2 / m_Resolution + 1.)));
+  m_GradientFilter->SetUseImageDirection(false);
 
   m_NeighborhoodScalarProductFilter->SetInput(m_GradientFilter->GetOutput());
 
-- 
GitLab