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
|
/*=========================================================================
Program: ORFEO Toolbox
Language: C++
Date: $Date$
Version: $Revision$
Copyright (c) Centre National d'Etudes Spatiales. All rights reserved.
See OTBCopyright.txt for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notices for more information.
=========================================================================*/
#include "otbWrapperApplication.h"
#include "otbWrapperApplicationFactory.h"
#include "otbForwardSensorModel.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::ForwardSensorModel<double> ModelType;
typedef itk::Point<double, 2> PointType;
private:
void DoInit() ITK_OVERRIDE
{
SetName("ConvertSensorToGeoPoint");
SetDescription("Sensor to geographic coordinates conversion.");
// Documentation
SetDocName("Convert Sensor Point To Geographic Point");
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");
}
void DoUpdateParameters() ITK_OVERRIDE
{
}
void DoExecute() ITK_OVERRIDE
{
// Get input Image
FloatVectorImageType::Pointer inImage = GetParameterImage("in");
// Instantiate a ForwardSensor Model
ModelType::Pointer model = ModelType::New();
model->SetImageGeometry(inImage->GetImageKeywordlist());
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)
|
