Skip to content
GitLab
Select Git revision
  • zonalstats_histogram default
Find file BlameHistoryPermalink
otbVertexComponentAnalysis.cxx 4.34 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
21 
#include "otbWrapperApplication.h"
22 
#include "otbWrapperApplicationFactory.h"
23
24 
#include "otbVcaImageFilter.h"
25
26 
namespace otb
27 
{
28 
namespace Wrapper
29 
{
30
31
32 
typedef otb::VCAImageFilter<DoubleVectorImageType>                 VCAFilterType;
33 
typedef otb::VectorImageToMatrixImageFilter<DoubleVectorImageType> VectorImageToMatrixImageFilterType;
34
35 
class VertexComponentAnalysis : public Application
36 
{
37 
public:
38 
  /** Standard class typedefs. */
39 
  typedef VertexComponentAnalysis         Self;
40 
  typedef Application                   Superclass;
41 
  typedef itk::SmartPointer<Self>       Pointer;
42 
  typedef itk::SmartPointer<const Self> ConstPointer;
43
44 
  /** Standard macro */
45 
  itkNewMacro(Self);
46
47 
  itkTypeMacro(VertexComponentAnalysis, otb::Application);
48
49 
private:
50 
  void DoInit() override
51 
  {
52 
    SetName("VertexComponentAnalysis");
53 
    SetDescription("Given a set of mixed spectral vectors, estimate "
54 
    "reference substances also known as endmembers using the Vertex "
55 
    "Component Analysis algorithm.");
56
57 
    // Documentation
58 
    SetDocName("Vertex Component Analysis");
59 
    SetDocLongDescription("Apply the Vertex Component Analysis [1] to "
60 
    "an hyperspectral image to extract endmembers. Given a set of mixed "
61 
    "spectral vectors (multispectral or hyperspectral), the application "
62 
    "estimates the spectral signature of reference substances also known "
63 
    "as endmembers.");
64 
    SetDocLimitations("None");
65 
    SetDocAuthors("OTB-Team");
66 
    SetDocSeeAlso("[1] J. M. P. Nascimento and J. M. B. Dias, Vertex "
67 
    "component analysis: a fast algorithm to unmix hyperspectral data, "
68 
    "in IEEE Transactions on Geoscience and Remote Sensing, vol. 43, "
69 
    "no. 4, pp. 898-910, April 2005.J. M. P. Nascimento and "
70 
    "J. M. B. Dias, Vertex component analysis: a fast algorithm to "
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129
"unmix hyperspectral data, in IEEE Transactions on Geoscience and " "Remote Sensing, vol. 43, no. 4, pp. 898-910, April 2005."); AddDocTag(Tags::Hyperspectral); AddDocTag(Tags::DimensionReduction); AddParameter(ParameterType_InputImage, "in", "Input Image"); SetParameterDescription("in","Input hyperspectral data cube"); AddParameter(ParameterType_OutputImage, "outendm", "Output Endmembers"); SetParameterDescription("outendm","Endmembers, stored in a " "one-line multi-spectral image.Each pixel corresponds to one " "endmembers and each band values corresponds to the spectral " "signature of the corresponding endmember."); MandatoryOn("outendm"); AddParameter(ParameterType_Int, "ne", "Number of endmembers"); SetParameterDescription("ne", "The number of endmembers to extract from the hyperspectral image."); SetParameterInt("ne", 1); MandatoryOn("ne"); AddRANDParameter(); // Doc example parameter settings SetDocExampleParameterValue("in", "cupriteSubHsi.tif"); SetDocExampleParameterValue("ne", "5"); SetDocExampleParameterValue("outendm", "VertexComponentAnalysis.tif double"); SetOfficialDocLink(); } void DoUpdateParameters() override { // Nothing to do here : all parameters are independent } void DoExecute() override { DoubleVectorImageType::Pointer inputImage = GetParameterDoubleVectorImage("in"); DoubleVectorImageType::Pointer endmembersImage; const unsigned int nbEndmembers = GetParameterInt("ne"); VCAFilterType::Pointer vca = VCAFilterType::New(); vca->SetNumberOfEndmembers(nbEndmembers); vca->SetInput(inputImage); endmembersImage = vca->GetOutput(); m_Ref = vca.GetPointer(); SetParameterOutputImage<DoubleVectorImageType>("outendm", endmembersImage); } itk::LightObject::Pointer m_Ref; }; } } OTB_APPLICATION_EXPORT(otb::Wrapper::VertexComponentAnalysis)