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
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
|
/*
* 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.
*/
#include "otbWrapperApplication.h"
#include "otbWrapperApplicationFactory.h"
#include "otbRPCForwardTransform.h"
#include "otbCoordinateToName.h"
namespace otb
{
namespace Wrapper
{
class ConvertSensorToGeoPoint : public Application
{
public:
/** Standard class typedefs. */
typedef ConvertSensorToGeoPoint Self;
typedef Application Superclass;
typedef itk::SmartPointer<Self> Pointer;
typedef itk::SmartPointer<const Self> ConstPointer;
/** Standard macro */
itkNewMacro(Self);
itkTypeMacro(ConvertSensorToGeoPoint, otb::Application);
/** Filters typedef */
typedef otb::RPCForwardTransform<double, 2, 2> ModelType;
typedef itk::Point<double, 2> PointType;
private:
void DoInit() override
{
SetName("ConvertSensorToGeoPoint");
SetDescription("Sensor to geographic coordinates conversion.");
// Documentation
SetDocLongDescription(
"This Application converts a sensor point of an input image to a geographic point using the Forward Sensor Model of the input image.");
SetDocLimitations("None");
SetDocAuthors("OTB-Team");
SetDocSeeAlso("ConvertCartoToGeoPoint application, otbObtainUTMZoneFromGeoPoint");
AddDocTag(Tags::Geometry);
AddParameter(ParameterType_InputImage, "in", "Sensor image");
SetParameterDescription("in", "Input sensor image.");
AddParameter(ParameterType_Group, "input", "Point Coordinates");
AddParameter(ParameterType_Float, "input.idx", "X value of desired point");
SetParameterDescription("input.idx", "X coordinate of the point to transform.");
SetDefaultParameterFloat("input.idx", 0.0);
AddParameter(ParameterType_Float, "input.idy", "Y value of desired point");
SetParameterDescription("input.idy", "Y coordinate of the point to transform.");
SetDefaultParameterFloat("input.idy", 0.0);
// Output with Output Role
AddParameter(ParameterType_Group, "output", "Geographic Coordinates");
AddParameter(ParameterType_Float, "output.idx", "Output Point Longitude");
SetParameterDescription("output.idx", "Output point longitude coordinate.");
AddParameter(ParameterType_Float, "output.idy", "Output Point Latitude");
SetParameterDescription("output.idy", "Output point latitude coordinate.");
AddParameter(ParameterType_String, "output.town", "Main town near the coordinates computed");
SetParameterDescription("output.town", "Nearest main town of the computed geographic point.");
AddParameter(ParameterType_String, "output.country", "Country of the image");
SetParameterDescription("output.country", "Country of the input image");
// Set the parameter role for the output parameters
SetParameterRole("output.idx", Role_Output);
SetParameterRole("output.idy", Role_Output);
SetParameterRole("output.town", Role_Output);
SetParameterRole("output.country", Role_Output);
// Doc example parameter settings
SetDocExampleParameterValue("in", "QB_TOULOUSE_MUL_Extract_500_500.tif");
SetDocExampleParameterValue("input.idx", "200");
SetDocExampleParameterValue("input.idy", "200");
SetOfficialDocLink();
}
void DoUpdateParameters() override
{
}
void DoExecute() override
{
// Get input Image
FloatVectorImageType::Pointer inImage = GetParameterImage("in");
// Instantiate a ForwardSensor Model
ModelType::Pointer model = ModelType::New();
model->SetMetadata(inImage->GetImageMetadata());
if (model->IsValidSensorModel() == false)
{
itkGenericExceptionMacro(<< "Unable to create a model");
}
// Convert the desired point
itk::ContinuousIndex<double, 2> inIndex;
PointType point;
inIndex[0] = GetParameterFloat("input.idx");
inIndex[1] = GetParameterFloat("input.idy");
inImage->TransformContinuousIndexToPhysicalPoint(inIndex, point);
ModelType::OutputPointType outputPoint;
outputPoint = model->TransformPoint(point);
// Set the value computed
SetParameterFloat("output.idx", outputPoint[0]);
SetParameterFloat("output.idy", outputPoint[1]);
// Set the town and the neaerest city
CoordinateToName::Pointer coord2name = CoordinateToName::New();
coord2name->SetLon(outputPoint[0]);
coord2name->SetLat(outputPoint[1]);
coord2name->Evaluate();
SetParameterString("output.town", coord2name->GetPlaceName());
SetParameterString("output.country", coord2name->GetCountryName());
}
};
}
}
OTB_APPLICATION_EXPORT(otb::Wrapper::ConvertSensorToGeoPoint)
|
