File: otbWrapperImageInterface.cxx

package info (click to toggle)
otb 8.1.1%2Bdfsg-1
  • links: PTS, VCS
  • area: main
  • in suites: bookworm
  • size: 1,030,436 kB
  • sloc: xml: 231,007; cpp: 224,490; ansic: 4,592; sh: 1,790; python: 1,131; perl: 92; makefile: 72
file content (92 lines) | stat: -rw-r--r-- 2,900 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
/*
 * Copyright (C) 2005-2022 Centre National d'Etudes Spatiales (CNES)
 *
 * This file is part of Orfeo Toolbox
 *
 *     https://www.orfeo-toolbox.org/
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *     http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */

#if defined(_MSC_VER)
#pragma warning(disable : 4786)
#endif

#include "otbWrapperApplicationRegistry.h"


int otbWrapperImageInterface(int argc, char* argv[])
{
  if (argc < 4)
  {
    std::cerr << "Usage: " << argv[0] << " application_path  input_image  output_txt" << std::endl;
    return EXIT_FAILURE;
  }

  std::string path   = argv[1];
  std::string input  = argv[2];
  std::string output = argv[3];

  otb::Wrapper::ApplicationRegistry::SetApplicationPath(path);

  otb::Wrapper::Application::Pointer app1 = otb::Wrapper::ApplicationRegistry::CreateApplication("Smoothing");

  if (app1.IsNull())
  {
    std::cerr << "Failed to create application" << std::endl;
    return EXIT_FAILURE;
  }

  std::ofstream ofs(output);
  if (!ofs.is_open())
  {
    fprintf(stderr, "Error, can't open file");
    return EXIT_FAILURE;
  }

  app1->SetParameterString("in", input);
  app1->Execute();

  ofs << "Size: " << app1->GetImageSize("out") << std::endl;
  ofs << "Origin: " << app1->GetImageOrigin("out") << std::endl;
  ofs << "Spacing: " << app1->GetImageSpacing("out") << std::endl;

  ofs << "ProjectionRef:" << std::endl;
  ofs << app1->GetImageProjection("out") << std::endl;

  itk::MetaDataDictionary dict = app1->GetMetadataDictionary("out");
  ofs << "Dictionary keys:" << std::endl;
  for (auto& key : dict.GetKeys())
  {
    ofs << "  - " << key << std::endl;
  }

  otb::Wrapper::ImageBaseType::RegionType region;
  region.SetIndex(0, 10);
  region.SetIndex(1, 10);
  region.SetSize(0, 7);
  region.SetSize(1, 7);

  ofs << "RAM usage (in Bytes): " << app1->PropagateRequestedRegion("out", region) << std::endl;
  ofs << "Input requested: " << app1->GetImageRequestedRegion("in") << std::endl;

  otb::Wrapper::Application::Pointer app2 = otb::Wrapper::ApplicationRegistry::CreateApplication("ConcatenateImages");
  app2->AddParameterStringList("il", input);
  app2->AddParameterStringList("il", input);
  app2->Execute();

  ofs << "Number of bands [il]: " << app2->GetImageNbBands("il", 0) << "+" << app2->GetImageNbBands("il", 1) << std::endl;
  ofs << "Number of bands [out]: " << app2->GetImageNbBands("out") << std::endl;
  ofs.close();
  return EXIT_SUCCESS;
}