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
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
|
/*=========================================================================
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 "otbGenericRSResampleImageFilter.h"
#include "otbImageFileReader.h"
#include "otbImageFileWriter.h"
#include <ogr_spatialref.h>
// Extract ROI
#include "otbMultiChannelExtractROI.h"
// Images definition
const unsigned int Dimension = 2;
typedef double PixelType;
typedef otb::VectorImage<PixelType, Dimension> ImageType;
typedef ImageType::SizeType SizeType;
typedef otb::GenericRSResampleImageFilter<ImageType,
ImageType> ImageResamplerType;
typedef ImageResamplerType::OriginType OriginType;
typedef ImageResamplerType::SpacingType SpacingType;
typedef otb::ImageFileReader<ImageType> ReaderType;
typedef otb::ImageFileWriter<ImageType> WriterType;
int otbGenericRSResampleImageFilterNew(int itkNotUsed(argc), char * itkNotUsed(argv)[])
{
// SmartPointer instantiation
ImageResamplerType::Pointer resampler = ImageResamplerType::New();
std::cout << resampler << std::endl;
return EXIT_SUCCESS;
}
int otbGenericRSResampleImageFilter(int itkNotUsed(argc), char* argv[])
{
// SmartPointer instantiation
ImageResamplerType::Pointer resampler = ImageResamplerType::New();
const char * infname = argv[1];
const char * outfname = argv[6];
unsigned int isize = atoi(argv[2]);
double iGridSpacing = atof(argv[3]);
int useInRpc = atoi(argv[4]);
int useOutRpc = atoi(argv[5]);
ReaderType::Pointer reader = ReaderType::New();
// Read the input image
reader->SetFileName(infname);
reader->UpdateOutputInformation();
// Fill the output size with the user selection
SizeType size;
size.Fill(isize);
// Set the origin & the spacing of the output
OriginType origin;
origin[0] = 367340;
origin[1] = 4.83467e+06;
SpacingType spacing;
spacing[0] = 0.6;
spacing[1] = -0.6;
// Build the output projection ref : UTM ref
OGRSpatialReference oSRS;
oSRS.SetProjCS("UTM");
oSRS.SetUTM(31, true);
char * utmRef = ITK_NULLPTR;
oSRS.exportToWkt(&utmRef);
// Displacement Field spacing
SpacingType gridSpacing;
gridSpacing[0] = iGridSpacing;
gridSpacing[1] = -iGridSpacing;
// Default value builder
ImageType::PixelType defaultValue;
itk::NumericTraits<ImageType::PixelType>::SetLength(defaultValue, reader->GetOutput()->GetNumberOfComponentsPerPixel());
// Set the Resampler Parameters
resampler->SetInput(reader->GetOutput());
resampler->SetDisplacementFieldSpacing(gridSpacing);
resampler->SetOutputOrigin(origin);
resampler->SetOutputSize(size);
resampler->SetOutputSpacing(spacing);
resampler->SetOutputProjectionRef(utmRef);
resampler->SetEdgePaddingValue(defaultValue);
if (useInRpc)
{
resampler->SetInputRpcGridSize(20);
resampler->EstimateInputRpcModelOn();
}
if (useOutRpc)
{
resampler->SetOutputRpcGridSize(20);
resampler->EstimateOutputRpcModelOn();
}
// Write the resampled image
WriterType::Pointer writer= WriterType::New();
writer->SetNumberOfDivisionsTiledStreaming(4);
writer->SetFileName(outfname);
writer->SetInput(resampler->GetOutput());
writer->Update();
std::cout << resampler << std::endl;
return EXIT_SUCCESS;
}
int otbGenericRSResampleImageFilterFromMap(int itkNotUsed(argc), char* argv[])
{
typedef otb::MultiChannelExtractROI<PixelType, PixelType> ExtractROIType;
// SmartPointer instantiation
ExtractROIType::Pointer extractor = ExtractROIType::New();
ImageResamplerType::Pointer resampler = ImageResamplerType::New();
const char * infname = argv[1];
const char * outfname = argv[4];
double iGridSpacing = atof(argv[2]);
int useInRpc = atoi(argv[3]);
// Reader Instantiation
ReaderType::Pointer reader = ReaderType::New();
reader->SetFileName(infname);
reader->UpdateOutputInformation();
SpacingType spacing;
spacing[0] = 2.5;
spacing[1] = -2.5;
// Displacement Field spacing
SpacingType gridSpacing;
gridSpacing[0] = iGridSpacing;
gridSpacing[1] = -iGridSpacing;
// Default value builder
ImageType::PixelType defaultValue;
itk::NumericTraits<ImageType::PixelType>::SetLength(defaultValue, reader->GetOutput()->GetNumberOfComponentsPerPixel());
// Extract a roi centered on the input center
ImageType::RegionType roi;
ImageType::IndexType roiIndex;
SizeType roiSize;
// Fill the size
roiSize.Fill(250);
// Fill the start index
roiIndex[0] = (unsigned int)((reader->GetOutput()->GetLargestPossibleRegion().GetSize()[0] - roiSize[0]) /2);
roiIndex[1] = (unsigned int)((reader->GetOutput()->GetLargestPossibleRegion().GetSize()[1] - roiSize[1]) /2);
roi.SetIndex(roiIndex);
roi.SetSize(roiSize);
extractor->SetExtractionRegion(roi);
extractor->SetInput(reader->GetOutput());
extractor->UpdateOutputInformation();
// Set the Resampler Parameters
resampler->SetInput(extractor->GetOutput());
resampler->SetDisplacementFieldSpacing(gridSpacing);
resampler->SetOutputParametersFromMap("UTM", spacing);
if (useInRpc)
{
resampler->SetInputRpcGridSize(20);
resampler->EstimateInputRpcModelOn();
}
// Write the resampled image
WriterType::Pointer writer= WriterType::New();
writer->SetAutomaticTiledStreaming();
writer->SetFileName(outfname);
writer->SetInput(resampler->GetOutput());
writer->Update();
std::cout << resampler << std::endl;
return EXIT_SUCCESS;
}
|
