Skip to content
GitLab
Select Git revision
  • master default
Find file BlameHistoryPermalink
otbTimeSeriesStats.cxx 6.19 KiB
1 
/*=========================================================================
2
3 
  Copyright (c) Remi Cresson (IRSTEA). All rights reserved.
4
5
6 
     This software is distributed WITHOUT ANY WARRANTY; without even
7 
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
8 
     PURPOSE.  See the above copyright notices for more information.
9
10 
=========================================================================*/
11 
#include "itkFixedArray.h"
12 
#include "itkObjectFactory.h"
13 
#include "otbWrapperApplicationFactory.h"
14
15 
// Application engine
16 
#include "otbStandardFilterWatcher.h"
17 
#include "itkFixedArray.h"
18
19 
// LayerStack
20 
#include "otbImageListToVectorImageFilter.h"
21 
#include "otbMultiToMonoChannelExtractROI.h"
22 
#include "otbImageList.h"
23
24 
// Extract multichannel
25 
#include "otbMultiChannelExtractROI.h"
26
27 
// Input images list
28 
#include "otbWrapperInputImageListParameter.h"
29
30 
// Functors for NDVI operations
31 
#include "otbFunctorImageFilter.h"
32 
#include "otbTimeSeriesFunctor.h"
33
34 
// Dates
35 
#include "otbDates.h"
36
37 
namespace otb
38 
{
39
40 
namespace Wrapper
41 
{
42
43 
class TimeSeriesStats : public Application
44 
{
45 
public:
46 
  /** Standard class typedefs. */
47 
  typedef TimeSeriesStats            Self;
48 
  typedef Application                                Superclass;
49 
  typedef itk::SmartPointer<Self>                    Pointer;
50 
  typedef itk::SmartPointer<const Self>              ConstPointer;
51
52 
  /** Standard macro */
53 
  itkNewMacro(Self);
54 
  itkTypeMacro(TimeSeriesStats, Application);
55
56 
  /** Typedefs for image concatenation */
57 
  typedef otb::ImageList<FloatImageType>                                        ImageListType;
58 
  typedef ImageListToVectorImageFilter<ImageListType, FloatVectorImageType>     ListConcatenerFilterType;
59 
  typedef MultiToMonoChannelExtractROI<FloatVectorImageType::InternalPixelType,
60 
      FloatImageType::PixelType>                                                ExtractROIFilterType;
61 
  typedef ObjectList<ExtractROIFilterType>                                      ExtractROIFilterListType;
62
63 
  /** Typedefs for temporal correlation */
64 
  typedef unsigned char                                                         LabelValueType;
65 
  typedef otb::VectorImage<LabelValueType>                                      LabelVectorImageType;
66 
  typedef otb::dates::SingleDate                                                DateType;
67 
  typedef otb::Functor::MultitemporalDensityFunctor<FloatVectorImageType::PixelType,
68 
      LabelVectorImageType::PixelType>                                          DensityFunctorType;
69 
  typedef otb::FunctorImageFilter<DensityFunctorType>                           DensityFilterType;
70
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140
void DoUpdateParameters() { // Nothing to do here : all parameters are independent } void DoInit() { // Documentation SetName("TimeSeriesStats"); SetDescription("Compute time series statistics"); SetDocLongDescription("This application computes time series statistics. The first component of the output image is the number of images different than no-data."); SetDocLimitations("None"); SetDocAuthors("Remi Cresson"); // Input time series images AddParameter(ParameterType_InputImageList, "il", "Input time series"); SetParameterDescription ("il", "The user can specify multiple input images, or one big stack of images."); // Input nodata AddParameter(ParameterType_Float, "nodata", "Input no-data value"); MandatoryOff ("nodata"); SetDefaultParameterFloat ("nodata", 0.0); // Output image AddParameter(ParameterType_OutputImage, "out", "Output label image"); SetParameterDescription ("out", "Output label image"); SetDefaultOutputPixelType ("out", ImagePixelType_uint8); AddRAMParameter(); } /* * Execute the process */ void DoExecute() { // Get the input image list FloatVectorImageListType::Pointer inList = this->GetParameterImageList("il"); if( inList->Size() < 2 ) { otbAppLogFATAL("At least two images are required."); } if ( inList->Size() >= itk::NumericTraits<LabelValueType>::max()) { otbAppLogFATAL("The maximum number of input images must be below " << itk::NumericTraits<LabelValueType>::max() << ". Consider changing LabelValueType."); } // Create one stack for input NDVI images list m_Concatener = ListConcatenerFilterType::New(); m_ImageList = ImageListType::New(); m_ExtractorList = ExtractROIFilterListType::New(); // Split each input vector image into image // and generate an mono channel image list inList->GetNthElement(0)->UpdateOutputInformation(); FloatVectorImageType::SizeType size = inList->GetNthElement(0)->GetLargestPossibleRegion().GetSize(); for( unsigned int i=0; i<inList->Size(); i++ ) { FloatVectorImageType::Pointer vectIm = inList->GetNthElement(i); vectIm->UpdateOutputInformation(); if( size != vectIm->GetLargestPossibleRegion().GetSize() ) { itkExceptionMacro("Input Image size mismatch..."); } for( unsigned int j=0; j<vectIm->GetNumberOfComponentsPerPixel(); j++)
141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175
{ ExtractROIFilterType::Pointer extractor = ExtractROIFilterType::New(); extractor->SetInput( vectIm ); extractor->SetChannel( j+1 ); extractor->UpdateOutputInformation(); m_ExtractorList->PushBack( extractor ); m_ImageList->PushBack( extractor->GetOutput() ); } } // Concatenate the images m_Concatener->SetInput( m_ImageList ); m_Concatener->UpdateOutputInformation(); // Count images m_TemporalFilter = DensityFilterType::New(); m_TemporalFilter->SetInput(m_Concatener->GetOutput()); m_TemporalFilter->GetModifiableFunctor().SetInputNoDataValue(GetParameterFloat("nodata")); SetParameterOutputImage("out", m_TemporalFilter->GetOutput()); } // DoExecute() ListConcatenerFilterType::Pointer m_Concatener; ExtractROIFilterListType::Pointer m_ExtractorList; ImageListType::Pointer m_ImageList; DensityFilterType::Pointer m_TemporalFilter; }; // end of class } // namespace wrapper } // namespace otb OTB_APPLICATION_EXPORT( otb::Wrapper::TimeSeriesStats )