Skip to content
GitLab
Select Git revision
Find file BlameHistoryPermalink
otbVectorDimensionalityReduction.cxx 16.20 KiB
1 
/*
2 
 * Copyright (C) 2005-2017 Centre National d'Etudes Spatiales (CNES)
3 
 *
4 
 * This file is part of Orfeo Toolbox
5 
 *
6 
 *     https://www.orfeo-toolbox.org/
7 
 *
8 
 * Licensed under the Apache License, Version 2.0 (the "License");
9 
 * you may not use this file except in compliance with the License.
10 
 * You may obtain a copy of the License at
11 
 *
12 
 *     http://www.apache.org/licenses/LICENSE-2.0
13 
 *
14 
 * Unless required by applicable law or agreed to in writing, software
15 
 * distributed under the License is distributed on an "AS IS" BASIS,
16 
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
17 
 * See the License for the specific language governing permissions and
18 
 * limitations under the License.
19 
 */
20 
#include "otbWrapperApplication.h"
21 
#include "otbWrapperApplicationFactory.h"
22 
#include "otbOGRDataSourceWrapper.h"
23 
#include "otbOGRFeatureWrapper.h"
24 
#include "itkVariableLengthVector.h"
25 
#include "otbStatisticsXMLFileReader.h"
26 
#include "itkListSample.h"
27 
#include "otbShiftScaleSampleListFilter.h"
28 
#include "otbDimensionalityReductionModelFactory.h"
29 
#include <time.h>
30
31 
namespace otb
32 
{
33 
namespace Wrapper
34 
{
35
36 
/** Utility function to negate std::isalnum */
37 
bool IsNotAlphaNum(char c)
38 
{
39 
return !std::isalnum(c);
40 
}
41
42 
/**
43 
 * \class VectorDimensionalityReduction
44 
 *
45 
 * Apply a dimensionality reduction model on a vector file
46 
 */
47 
class VectorDimensionalityReduction : public Application
48 
{
49 
public:
50 
    /** Standard class typedefs. */
51 
  typedef VectorDimensionalityReduction Self;
52 
  typedef Application Superclass;
53 
  typedef itk::SmartPointer<Self> Pointer;
54 
  typedef itk::SmartPointer<const Self> ConstPointer;
55
56 
  /** Standard macro */
57 
  itkNewMacro(Self);
58 
  itkTypeMacro(Self, Application)
59
60 
  /** Filters typedef */
61 
  typedef float                                           ValueType;
62 
  typedef itk::VariableLengthVector<ValueType>            InputSampleType;
63 
  typedef itk::Statistics::ListSample<InputSampleType>    ListSampleType;
64 
  typedef MachineLearningModel<
65 
    itk::VariableLengthVector<ValueType>,
66 
    itk::VariableLengthVector<ValueType> >                DimensionalityReductionModelType;
67 
  typedef DimensionalityReductionModelFactory<
68 
    ValueType,ValueType>                                  DimensionalityReductionModelFactoryType;
69 
  typedef DimensionalityReductionModelType::Pointer       ModelPointerType;
70
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140
/** Statistics Filters typedef */ typedef itk::VariableLengthVector<ValueType> MeasurementType; typedef otb::StatisticsXMLFileReader<MeasurementType> StatisticsReader; typedef otb::Statistics::ShiftScaleSampleListFilter< ListSampleType, ListSampleType> ShiftScaleFilterType; protected: ~VectorDimensionalityReduction() override { DimensionalityReductionModelFactoryType::CleanFactories(); } private: void DoInit() override { SetName("VectorDimensionalityReduction"); SetDescription("Performs dimensionality reduction of the input vector data " "according to a model file."); SetDocName("Vector Dimensionality Reduction"); SetDocAuthors("OTB-Team"); SetDocLongDescription("This application performs a vector data " "dimensionality reduction based on a model file produced by the " "TrainDimensionalityReduction application."); SetDocSeeAlso("TrainDimensionalityReduction"); SetDocLimitations("None"); AddDocTag(Tags::Learning); AddParameter(ParameterType_InputVectorData, "in", "Name of the input vector data"); SetParameterDescription("in","The input vector data to reduce."); AddParameter(ParameterType_InputFilename, "instat", "Statistics file"); SetParameterDescription("instat", "A XML file containing mean and standard " "deviation to center and reduce samples before dimensionality reduction " "(produced by ComputeImagesStatistics application)."); MandatoryOff("instat"); AddParameter(ParameterType_InputFilename, "model", "Model file"); SetParameterDescription("model", "A model file (produced by the " "TrainDimensionalityReduction application,"); AddParameter(ParameterType_OutputFilename, "out", "Output vector data file " "containing the reduced vector"); SetParameterDescription("out","Output vector data file storing sample " "values (OGR format). If not given, the input vector data file is used. " "In overwrite mode, the original features will be lost."); MandatoryOff("out"); AddParameter(ParameterType_ListView, "feat", "Input features to use for reduction"); SetParameterDescription("feat", "List of field names in the input vector data used as features for reduction."); AddParameter(ParameterType_Choice, "featout", "Output feature"); SetParameterDescription("featout", "Naming of output features"); AddChoice("featout.prefix", "Prefix"); SetParameterDescription("featout.prefix", "Use a name prefix"); AddParameter(ParameterType_String, "featout.prefix.name", "Feature name prefix"); SetParameterDescription("featout.prefix.name","Name prefix for output " "features. This prefix is followed by the numeric index of each output feature."); SetParameterString("featout.prefix.name","reduced_", false); AddChoice("featout.list","List"); SetParameterDescription("featout.list", "Use a list with all names"); AddParameter(ParameterType_StringList, "featout.list.names", "Feature name list"); SetParameterDescription("featout.list.names", "List of field names for the output " "features which result from the reduction."); AddParameter(ParameterType_Int, "pcadim", "Principal component dimension");