From dea96df77ca574d96ad1abe25f7ffc6b04e776d4 Mon Sep 17 00:00:00 2001
From: Gregoire Mercier <gregoire.mercier@telecom-bretagne.eu>
Date: Fri, 25 Mar 2011 16:28:17 +0100
Subject: [PATCH] ADD Sparse Unmixing filter
---
.../otbSparseUnmixingImageFilter.h | 124 ++++++++++++++++++
1 file changed, 124 insertions(+)
create mode 100644 Code/Hyperspectral/otbSparseUnmixingImageFilter.h
diff --git a/Code/Hyperspectral/otbSparseUnmixingImageFilter.h b/Code/Hyperspectral/otbSparseUnmixingImageFilter.h
new file mode 100644
index 0000000000..80ec6b3860
--- /dev/null
+++ b/Code/Hyperspectral/otbSparseUnmixingImageFilter.h
@@ -0,0 +1,124 @@
+/*=========================================================================
+
+ 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 __otbSparseUnmixingImageFilter_h
+#define __otbSparseUnmixingImageFilter_h
+
+#include "otbMacro.h"
+#include "itkImageToImageFilter.h"
+#include "otbWaveletTransformImageFilter.h"
+#include "otbSparseWvltToAngleMapperListFilter.h"
+#include "otbImageList.h"
+#include "itkListSample.h"
+#include "itkFixedArray.h"
+#include "itkHistogram.h"
+#include "otbMatrixMultiplyImageFilter.h"
+
+namespace otb {
+
+/** \class SparseUnmixingImageFilter
+ * \brief This class detects linear dependencies from 2 wavelet decompositions
+ *
+ * This class perform a linear unmixing from a sparse representation of 2
+ * signals (Through multiscale wavelet transforms) by means of line detection
+ * in the scatterplot.
+ *
+ * It is implemeted as a mutli-filter that performs wavelet decomposition of the 2
+ * images, Sample list construction of the scatter plot sample, histogram estimation
+ * of the phase and source detection...
+ *
+ */
+template < class TInputImage, class TOutputImage,
+ class TPixelPrecision = double,
+ Wavelet::Wavelet TMotherWaveletOperator = Wavelet::SYMLET8 >
+class ITK_EXPORT SparseUnmixingImageFilter
+ : public itk::ImageToImageFilter<TInputImage, TOutputImage>
+{
+public:
+ /** Standard typedefs */
+ typedef SparseUnmixingImageFilter Self;
+ typedef itk::ImageToImageFilter<TInputImage, TOutputImage> Superclass;
+ typedef itk::SmartPointer<Self> Pointer;
+ typedef itk::SmartPointer<const Self> ConstPointer;
+
+ /** Type macro */
+ itkNewMacro(Self);
+
+ /** Creation through object factory macro */
+ itkTypeMacro(SparseUnmixingImageFilter, ImageToImageFilter);
+
+ /** Dimension */
+ itkStaticConstMacro(InputImageDimension, unsigned int, TInputImage::ImageDimension);
+ itkStaticConstMacro(OutputImageDimension, unsigned int, TOutputImage::ImageDimension);
+
+ /** Template parameters typedefs */
+ typedef TInputImage InputImageType;
+ typedef TOutputImage OutputImageType;
+ typedef TPrecision PrecisionType;
+ typedef TMotherWaveletOperator MotherWaveletOperatorID;
+
+ /** Filter types and related */
+ typedef Image< PrecisionType, InputImageDimension > InternalImageType;
+
+ typedef WaveletOperator< MotherWaveletOperatorID, Wavelet::FORWARD, Precision, InputImageDimension > WaveletOperator;
+ typedef WaveletFilterBank< InputImageType, InternalImageType, WaveletOperator, Wavelet::FORWARD > FilterBankType;
+ typedef WaveletTransform< InputImageType, InternalImageType, FilterBankType, Wavelet::FORWARD > WvltFilterType;
+ typedef typename WvltFilterType::Pointer WvltFilterPointerType;
+ typedef typename WvltFilterType::OutputImageListType InternalImageListType;
+
+ typedef itk::Statistics::ListSample< itk::FixedArray< PrecisionType, 1 > > SampleType;
+ typedef SparseWvltToAngleMapperListFilter< InternalImageListType, SampleType > ListFilterType;
+ typedef typename ListFilterType::Pointer ListFilterPointerType;
+
+ typedef typename itk::Statistics::Histogram< PrecisionType > HistogramType;
+ typedef typename HistogramType::Pointer HistogramPointerType;
+
+ typedef MatrixMultiplyImageFilter< InternalImageType, OutputImageType, PrecisionType > TransformFilterType;
+ typedef typename TransformFilterType::Pointer TransformFilterPointerType;
+
+ /** Sets/Gets */
+ void SetNumberOfDecomposition ( unsigned int nb )
+ {
+ m_WvltFilter1->SetNumberOfDecomposition( nb );
+ m_WvltFilter2->SetNumberOfDecomposition( nb );
+ this->Modified();
+ }
+ unsigned int GetNumberOfDecomposition () const
+ {
+ return m_WvltFilter1->GetNumberOfDecomposition();
+ }
+
+protected:
+ SparseUnmixingImageFilter();
+ virtual ~SparseUnmixingImageFilter() { }
+
+private:
+ SparseUnmixingImageFilter(const Self &); // not implemented
+ void operator=(const Self &);
+
+
+ WvltFilterPointerType m_WvltFilter1;
+ WvltFilterPointerType m_WvltFilter2;
+ ListFilterPointerType m_ListFilter;
+ HistogramPointerType m_Histogram;
+ TransformFilterPointerType m_Transformer;
+}; // end of class
+
+} // end of namespace otb
+#endif
+
+
--
GitLab