*Ajout d'un labelize avec une croissance de region NeighborhoodConnected

Code/BasicFilters/otbLabelizeNeighborhoodConnectedImageFilter.h

+/*=========================================================================
+
+  Program:   ORFEO Toolbox
+  Language:  C++
+  Date:      $Date$
+  Version:   $Revision$
+
+
+  Copyright (c) Centre National d'Etudes Spatiales. All rights reserved.
+  See OTBCopyright.txt for details.
+
+
+     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.
+
+=========================================================================*/
+#ifndef __otbLabelizeNeighborhoodConnectedImageFilter_h
+#define __otbLabelizeNeighborhoodConnectedImageFilter_h
+
+#include "itkNeighborhoodConnectedImageFilter.h"
+#include "otbLabelizeImageFilterBase.h"
+
+namespace otb
+{
+
+  /** \class LabelizeNeighborhoodConnectedImageFilter
+   *  \brief Labels pixels that are connected to a seed and lie
+   *         within a neighborhood
+   *  
+   *  This class implements base class otb::LabelizeImageFilterBase
+   *
+   *  First seeds are obtained with a theshold filter
+   *  For each seed, itk::NeighborhoodConnectedImageFilter labels pixel
+   *  Finally, filter produce a regions map
+   */
+  template <class TInputImage, class TOutputImage>
+    class ITK_EXPORT LabelizeNeighborhoodConnectedImageFilter
+    : public otb::LabelizeImageFilterBase<TInputImage, TOutputImage, itk::NeighborhoodConnectedImageFilter<TInputImage, TOutputImage> >
+    {
+    public:
+      /** typedef for standard classes. */
+      typedef LabelizeNeighborhoodConnectedImageFilter Self;
+      typedef otb::LabelizeImageFilterBase<TInputImage,TOutputImage, itk::NeighborhoodConnectedImageFilter<TInputImage, TOutputImage> > Superclass;
+      typedef itk::SmartPointer<Self> Pointer;
+      typedef itk::SmartPointer<const Self>  ConstPointer;
+      
+      typedef TInputImage InputImageType;
+      typedef TOutputImage OutputImageType;
+
+      typedef typename TInputImage::PixelType InputPixelType;
+      typedef typename TOutputImage::PixelType OutputPixelType;
+      typedef typename TInputImage::IndexType IndexType;
+      
+      /** "object factory" management method. */
+      itkNewMacro(Self);
+      
+      /** Return the class name. */
+      itkTypeMacro(LabelizeNeighborhoodConnectedImageFilter, LabelizeImageFilterBase);
+
+      /** Get delta lower threshold */
+      itkGetMacro(LowerThresholdDelta, InputPixelType);
+      
+      /** Set delta lower threshold */
+      itkSetMacro(LowerThresholdDelta, InputPixelType);
+      
+      /** Get delta upper threshold */
+      itkGetMacro(UpperThresholdDelta, InputPixelType);
+      
+      /** Set delta upper threshold */
+      itkSetMacro(UpperThresholdDelta, InputPixelType);
+      
+      /** Get replace value */
+      itkGetMacro(ReplaceValue, OutputPixelType);
+      
+      /** Set replace value */
+      itkSetMacro(ReplaceValue, OutputPixelType);
+      
+      /** typedef to simplify variables definition and declaration. */
+      typedef TInputImage InputImageType;
+      typedef TOutputImage OutputImageType;
+      
+    protected:
+      LabelizeNeighborhoodConnectedImageFilter();
+      virtual ~LabelizeNeighborhoodConnectedImageFilter() {};
+      virtual void PrintSelf(std::ostream& os, itk::Indent indent) const;
+      
+      /** Region growing */
+      virtual void RegionGrowing( const IndexType indexSeed );
+      
+    private:
+      LabelizeNeighborhoodConnectedImageFilter(const Self&); //purposely not implemented
+      void operator=(const Self&); //purposely not implemented
+      
+      /** Delta + threshold for growing region */
+      InputPixelType m_UpperThresholdDelta;
+      
+      /** Delta - threshold for growing region */
+      InputPixelType m_LowerThresholdDelta;
+      
+      /** Intial replace value*/
+      OutputPixelType m_ReplaceValue;
+      
+    }; // end class LabelizeconnectedThresholdImageFilter
+  
+} // end namespace otb
+
+#ifndef OTB_MANUAL_INSTANTIATION
+#include "otbLabelizeNeighborhoodConnectedImageFilter.txx"
+#endif
+
+#endif

Code/BasicFilters/otbLabelizeNeighborhoodConnectedImageFilter.txx

+/*=========================================================================
+
+Program:   ORFEO Toolbox
+Language:  C++
+Date:      $Date$
+Version:   $Revision$
+
+
+Copyright (c) Centre National d'Etudes Spatiales. All rights reserved.
+See OTBCopyright.txt for details.
+
+
+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.
+
+=========================================================================*/
+
+#include "otbLabelizeNeighborhoodConnectedImageFilter.h"
+
+namespace otb
+{
+  /**
+   * Constructor
+   */
+  template  <class TInputImage, class TOutputImage>
+  LabelizeNeighborhoodConnectedImageFilter<TInputImage, TOutputImage>
+  ::LabelizeNeighborhoodConnectedImageFilter()
+  {
+    m_LowerThresholdDelta = 10;
+    m_UpperThresholdDelta = 10;
+    
+    m_ReplaceValue = 0;
+  }
+  
+  /** Region growing
+   *
+   */
+  template <class TInputImage, class TOutputImage>
+  void
+  LabelizeNeighborhoodConnectedImageFilter<TInputImage, TOutputImage>
+  ::RegionGrowing( const IndexType indexSeed )
+  {
+    InputPixelType threshold = this->GetInput()->GetPixel(indexSeed);
+    
+    this->m_RegionGrowingFilter->SetLower(threshold-m_LowerThresholdDelta);
+    this->m_RegionGrowingFilter->SetUpper(threshold+m_UpperThresholdDelta);
+    
+    this->m_RegionGrowingFilter->SetReplaceValue(m_ReplaceValue);
+    this->m_RegionGrowingFilter->SetSeed(indexSeed);
+    m_ReplaceValue++;
+  }
+  
+  /** PrintSelf Method
+   *
+   */
+  template <class TInputImage, class TOutputImage>
+  void
+  LabelizeNeighborhoodConnectedImageFilter<TInputImage, TOutputImage>
+  ::PrintSelf(std::ostream& os, itk::Indent indent) const
+  {
+    Superclass::PrintSelf(os, indent);
+    os << indent << "LowerThresholdDelta: " << m_LowerThresholdDelta << std::endl;
+    os << indent << "UpperThresholdDelta: " << m_UpperThresholdDelta << std::endl;
+    os << indent << "ReplaceValue: " << m_ReplaceValue << std::endl;
+  }
+} // end namespace otb