Skip to content
Snippets Groups Projects
Commit caf6f506 authored by Victor Poughon's avatar Victor Poughon
Browse files

REFAC: remove GeometriesFilter.cxx and GeometriesChangeSpatialReference.cxx

parent 62170473
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Examples/Filtering/CMakeLists.txt
+ 0
− 6
View file @ caf6f506
......@@ -26,12 +26,6 @@ target_link_libraries(CompositeFilterExample ${OTB_LIBRARIES})
add_executable(DanielssonDistanceMapImageFilter DanielssonDistanceMapImageFilter.cxx)
target_link_libraries(DanielssonDistanceMapImageFilter ${OTB_LIBRARIES})
add_executable(GeometriesChangeSpatialReference GeometriesChangeSpatialReference.cxx)
target_link_libraries(GeometriesChangeSpatialReference ${OTB_LIBRARIES})
add_executable(GeometriesFilter GeometriesFilter.cxx)
target_link_libraries(GeometriesFilter ${OTB_LIBRARIES})
add_executable(SecondDerivativeRecursiveGaussianImageFilter SecondDerivativeRecursiveGaussianImageFilter.cxx)
target_link_libraries(SecondDerivativeRecursiveGaussianImageFilter ${OTB_LIBRARIES})
......
Examples/Filtering/GeometriesChangeSpatialReference.cxx deleted 100644 → 0
+ 0
− 249
View file @ 62170473
/*
* Copyright (C) 2005-2019 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.
*/
/*===========================================================================*/
/*===============================[ Includes ]================================*/
/*===========================================================================*/
#include "otbGeometriesToGeometriesFilter.h"
#include <string>
#include <iostream>
#include "otbGeometriesSet.h"
#include "otbOGRGeometryWrapper.h"
/*===========================================================================*/
/*=========================[ Reprojection functor ]==========================*/
/*===========================================================================*/
namespace internal
{
struct Deleters
{
void operator()(OGRCoordinateTransformation *p) const { OCTDestroyCoordinateTransformation(p); }
void operator()(OGRSpatialReference *p) const
{
#if GDAL_VERSION_NUM >= 1700
OGRSpatialReference::DestroySpatialReference(p);
#else
#warning the following resource release may crash, please update your version of GDAL
delete p; // note there is no garanty
#endif
}
};
} // internal namespace
struct ReprojectTransformationFunctor
{
typedef OGRGeometry TransformedElementType;
otb::ogr::UniqueGeometryPtr operator()(OGRGeometry const* in) const;
void SetTransformation(OGRSpatialReference & in, OGRSpatialReference & out)
{
m_reprojector.reset(OGRCreateCoordinateTransformation(&in, &out));
}
private:
boost::interprocess::unique_ptr<OGRCoordinateTransformation, internal::Deleters> m_reprojector;
};
otb::ogr::UniqueGeometryPtr
ReprojectTransformationFunctor::operator()(OGRGeometry const* in) const
{
otb::ogr::UniqueGeometryPtr out(in ? in->clone() : nullptr);
if (out)
{
const OGRErr err = out->transform(m_reprojector.get());
if (err != OGRERR_NONE)
{
itkGenericExceptionMacro(<< "Cannot reproject a geometry: " << CPLGetLastErrorMsg());
}
}
return otb::move(out);
}
/*===========================================================================*/
/*==========================[ reprojection filter ]==========================*/
/*===========================================================================*/
class MyReprojectionFilter : public otb::DefaultGeometriesToGeometriesFilter<ReprojectTransformationFunctor>
{
public:
typedef MyReprojectionFilter Self;
typedef otb::DefaultGeometriesToGeometriesFilter<ReprojectTransformationFunctor> Superclass;
typedef itk::SmartPointer<Self> Pointer;
typedef itk::SmartPointer<const Self> ConstPointer;
itkNewMacro(Self);
itkTypeMacro(MyReprojectionFilter, otb::DefaultGeometriesToGeometriesFilter);
void SetOutputSpatialReference(std::string const& srs);
void UpdateTransformation(OGRSpatialReference &in)
{
(*this)->SetTransformation(in, *m_osr);
}
protected:
/** Default constructor. */
MyReprojectionFilter(){ }
/** Destructor. */
~MyReprojectionFilter() override{ }
private:
OGRSpatialReference* DoDefineNewLayerSpatialReference(otb::ogr::Layer const& itkNotUsed(source)) const override
{
return m_osr.get();
}
boost::interprocess::unique_ptr<OGRSpatialReference, internal::Deleters> m_osr;
};
void
MyReprojectionFilter::SetOutputSpatialReference(std::string const& srs)
{
OGRSpatialReference outSrs;
char *name = const_cast<char *>(srs.c_str());
OGRErr ret = outSrs.importFromWkt(&name);
if (ret)
{
ret = outSrs.importFromProj4(name);
}
if (ret)
{
ret = outSrs.importFromEPSG(atoi(name));
}
if (ret)
{
std::cout << "Error : spatial reference system not recognized ("<<srs<<")" << std::endl;
}
m_osr .reset (outSrs.Clone());
}
/*===========================================================================*/
/*=================================[ main ]==================================*/
/*===========================================================================*/
typedef MyReprojectionFilter FilterType;
struct Options
{
Options(int argc, char **argv)
{
workingInplace = true;
outputIsStdout = false;
for (int a=1; a!=argc; ++a)
{
if (!strcmp(argv[a], "-srs"))
{
if (a+1 < argc)
{
++a;
srs = argv[a];
}
}
else if (!strcmp(argv[a], "-"))
outputIsStdout = true;
else if (inputFile.empty())
inputFile = argv[a];
else if (outputFile.empty())
{
outputFile = argv[a];
workingInplace = false;
}
else
throw std::runtime_error(usage(argv[0])+"Too many (unexpected) arguments");
}
if (!outputFile.empty() && outputIsStdout)
throw std::runtime_error(usage(argv[0])+"An output file ("+outputFile+") has already been set, cannot dump to stdout");
if (inputFile.empty() || srs.empty())
throw std::runtime_error(usage(argv[0])+"Not enough parameters");
}
static std::string usage(std::string const& progname)
{
return progname + " <inputGeometriesFile> [<outputGeometriesFile>] -srs <id>\n";
}
std::string inputFile;
std::string outputFile;
std::string srs;
bool workingInplace;
bool outputIsStdout;
};
int main (int argc, char **argv)
{
if (argc < 2)
{
return EXIT_FAILURE;
}
try
{
const Options options(argc, argv);
otb::ogr::DataSource::Pointer input = otb::ogr::DataSource::New(
options.inputFile,
options.workingInplace ? otb::ogr::DataSource::Modes::Update_LayerOverwrite : otb::ogr::DataSource::Modes::Read);
otb::ogr::DataSource::Pointer output
= options.workingInplace ? input
: options.outputIsStdout ? nullptr
: otb::ogr::DataSource::New( options.outputFile, otb::ogr::DataSource::Modes::Update_LayerCreateOnly);
std::cout << "input: " << otb::ogr::GetFileListAsStringVector(&input->ogr())[0] << " should be: " << options.inputFile << "\n";
if (output)
{
std::cout << "output: " << otb::ogr::GetFileListAsStringVector(&output->ogr())[0] << " should be: " << options.outputFile << "\n";
}
// std::cout << "\n";
FilterType::Pointer filter = FilterType::New();
// TODO: this make no sense for in-place transformations
filter->SetOutputSpatialReference(options.srs);
filter->UpdateTransformation( *(const_cast<OGRSpatialReference*>(input->GetLayer(0).GetSpatialRef())) );
if (!options.workingInplace)
{
otb::GeometriesSet::Pointer in_set = otb::GeometriesSet::New(input);
filter->SetInput(in_set);
}
if (output)
{
otb::GeometriesSet::Pointer out_set = otb::GeometriesSet::New(output);
filter->SetOutput(out_set);
out_set->Update();
}
else
{
filter->UpdateOutputData(nullptr);
assert(filter->GetOutput() && "output not set");
filter->GetOutput()->Print(std::cout, 0);
}
return EXIT_SUCCESS;
}
catch (std::exception const& e)
{
std::cerr << e.what() << "\n";
}
return EXIT_FAILURE;
}
Examples/Filtering/GeometriesFilter.cxx deleted 100644 → 0
+ 0
− 112
View file @ 62170473
/*
* Copyright (C) 2005-2019 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.
*/
/*===========================================================================*/
/*===============================[ Includes ]================================*/
/*===========================================================================*/
#include "otbGeometriesToGeometriesFilter.h"
#include <string>
#include <iostream>
#include "otbGeometriesSet.h"
/*===========================================================================*/
/*=============================[ SwapXYFunctor ]=============================*/
/*===========================================================================*/
struct SwapXYFunctor
{
typedef OGRGeometry TransformedElementType;
otb::ogr::UniqueGeometryPtr operator()(OGRGeometry const* in) const
{
otb::ogr::UniqueGeometryPtr out(in ? in->clone() : nullptr);
if (out)
{
#if GDAL_VERSION_NUM >= 1900
out->swapXY();
#else
assert(!"Sorry, This example filter requires GDAL v1.9.0 or later");
#endif
}
return otb::move(out);
}
};
/*===========================================================================*/
/*=================================[ main ]==================================*/
/*===========================================================================*/
typedef otb::DefaultGeometriesToGeometriesFilter<SwapXYFunctor> FilterType;
int main (int argc, char **argv)
{
if (argc < 2)
{
std::cerr << argv[0] << " inputGeometriesFile [outputGeometriesFile]" << std::endl;
return EXIT_FAILURE;
}
try
{
const bool workingInplace = argc == 2;
const bool outputIsStdout = !workingInplace && !strcmp(argv[2], "-");
const std::string inputFile = argv[1];
const std::string outputFile = workingInplace ? argv[1] : argv[2];
otb::ogr::DataSource::Pointer input = otb::ogr::DataSource::New(
inputFile,
workingInplace ? otb::ogr::DataSource::Modes::Update_LayerOverwrite : otb::ogr::DataSource::Modes::Read);
otb::ogr::DataSource::Pointer output
= workingInplace ? input
: outputIsStdout ? nullptr
: otb::ogr::DataSource::New( outputFile, otb::ogr::DataSource::Modes::Update_LayerCreateOnly);
std::cout << "input: " << otb::ogr::GetFileListAsStringVector(&input->ogr())[0] << " should be: " << inputFile << "\n";
if (output)
{
std::cout << "output: " << otb::ogr::GetFileListAsStringVector(&output->ogr())[0] << " should be: " << outputFile << "\n";
}
// std::cout << "\n";
FilterType::Pointer filter = FilterType::New();
if (!workingInplace)
{
otb::GeometriesSet::Pointer in_set = otb::GeometriesSet::New(input);
filter->SetInput(in_set);
}
if (output)
{
otb::GeometriesSet::Pointer out_set = otb::GeometriesSet::New(output);
filter->SetOutput(out_set);
out_set->Update();
}
else
{
filter->UpdateOutputData(nullptr);
assert(filter->GetOutput() && "output not set");
filter->GetOutput()->Print(std::cout, 0);
}
return EXIT_SUCCESS;
}
catch (std::exception const& e)
{
std::cerr << e.what() << "\n";
}
return EXIT_FAILURE;
}
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment