From c1d3a15e7ba65b8b94ca07dfebb0e0be9e119ca6 Mon Sep 17 00:00:00 2001
From: Jonathan Guinet <jonathan.guinet@c-s.fr>
Date: Fri, 29 Jun 2012 15:42:19 +0200
Subject: [PATCH] ENH: VectorDataRepojection application.
---
Applications/Projections/CMakeLists.txt | 6 +
.../Projections/otbVectorDataReprojection.cxx | 194 ++++++++++++++++++
2 files changed, 200 insertions(+)
create mode 100644 Applications/Projections/otbVectorDataReprojection.cxx
diff --git a/Applications/Projections/CMakeLists.txt b/Applications/Projections/CMakeLists.txt
index eea4d18662..a66fa88484 100644
--- a/Applications/Projections/CMakeLists.txt
+++ b/Applications/Projections/CMakeLists.txt
@@ -30,6 +30,12 @@ OTB_CREATE_APPLICATION(NAME ObtainUTMZoneFromGeoPoint
OTB_CREATE_APPLICATION(NAME ConvertSensorToGeoPoint
SOURCES otbConvertSensorToGeoPoint.cxx
LINK_LIBRARIES OTBBasicFilters)
+
OTB_CREATE_APPLICATION(NAME GridBasedImageResampling
SOURCES otbGridBasedImageResampling.cxx
LINK_LIBRARIES OTBBasicFilters)
+
+OTB_CREATE_APPLICATION(NAME VectorDataReprojection
+ SOURCES otbVectorDataReprojection.cxx
+ LINK_LIBRARIES OTBBasicFilters;OTBIO;OTBOGRAdapters;OTBProjections)
+
\ No newline at end of file
diff --git a/Applications/Projections/otbVectorDataReprojection.cxx b/Applications/Projections/otbVectorDataReprojection.cxx
new file mode 100644
index 0000000000..3d5aa3459a
--- /dev/null
+++ b/Applications/Projections/otbVectorDataReprojection.cxx
@@ -0,0 +1,194 @@
+/*=========================================================================
+
+ 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 "otbWrapperNumericalParameter.h"
+#include "otbGeometriesProjectionFilter.h"
+#include "otbGeometriesSet.h"
+
+
+// MapProjection handler
+#include "otbWrapperMapProjectionParametersHandler.h"
+
+
+namespace otb
+{
+
+
+namespace Wrapper
+{
+
+class VectorDataReprojection : public Application
+{
+public:
+ /** Standard class typedefs. */
+ typedef VectorDataReprojection Self;
+ typedef Application Superclass;
+ typedef itk::SmartPointer<Self> Pointer;
+ typedef itk::SmartPointer<const Self> ConstPointer;
+
+ /** filter typedef **/
+ typedef otb::GeometriesProjectionFilter ProjectionFilterType;
+
+ typedef otb::GeometriesSet InputGeometriesType;
+ typedef otb::GeometriesSet OutputGeometriesType;
+
+ /** Standard macro */
+ itkNewMacro(Self)
+ ;
+
+ itkTypeMacro(VectorDataReprojection, otb::Application)
+ ;
+
+private:
+ void DoInit()
+ {
+ SetName("VectorDataReprojection");
+ std::ostringstream oss;
+ oss << "This application allows to reproject a vector data using support image projection reference"
+ ", or a user specified map projection" << std::endl;
+ SetDescription(oss.str());
+ // Documentation
+ SetDocName("Vector Data reprojection");
+ oss.str("");
+ oss <<" This application allows to reproject a vector data using support image projection reference"
+ ", or a user given map projection." << std::endl;
+ oss <<" If given, image keywordlist can be added to reprojected vectordata.";
+ SetDocLongDescription(oss.str());
+ SetDocLimitations(" ");
+ SetDocAuthors("OTB-Team");
+ SetDocSeeAlso(" ");
+
+ AddDocTag(Tags::Geometry);
+ AddDocTag(Tags::Vector);
+ AddDocTag(Tags::Coordinates);
+
+ // Set the parameters
+
+ AddParameter(ParameterType_Group, "in", "Input data");
+ AddParameter(ParameterType_InputFilename, "in.vd", "Input vector data");
+ SetParameterDescription("in.vd", "The input vector data to reproject");
+ AddParameter(ParameterType_InputImage, "in.kwl", "Use image keywords list");
+ SetParameterDescription("in.kwl", "Optional input image to fill vector data with image kwl.");
+ MandatoryOff("in.kwl");
+ DisableParameter("in.kwl");
+
+ AddParameter(ParameterType_Group, "out", "Output data");
+ AddParameter(ParameterType_OutputFilename, "out.vd", "Output vector data");
+ SetParameterDescription("out.vd", "The reprojected vector data");
+
+ // output projection choice
+ AddParameter(ParameterType_Choice, "out.proj", "Output Projection choice");
+ AddChoice("out.proj.image", "Use image projection ref");
+ SetParameterDescription("out.proj.image", "Vector data will be reprojected in image projection ref.");
+ AddParameter(ParameterType_InputImage, "out.proj.image.in", "Image used to get projection map");
+ SetParameterDescription("out.proj.image.in", "Projection map will be found using image metadata");
+ AddChoice("out.proj.user", "User defined projection");
+ MapProjectionParametersHandler::AddMapProjectionParameters(this, "out.proj.user.map");
+
+ // Doc example parameter settings
+ SetDocExampleParameterValue("in.vd", "VectorData_QB1.shp");
+ SetDocExampleParameterValue("out.proj", "image");
+ SetDocExampleParameterValue("out.proj.image", "ROI_QB_MUL_1.tif");
+ SetDocExampleParameterValue("out.vd","reprojected_vd.shp");
+ }
+
+ void DoUpdateParameters()
+ {
+
+ }
+
+ void DoExecute()
+ {
+ GetLogger()->Debug("Entering DoExecute");
+
+ // Get the input image
+
+ otb::ogr::DataSource::Pointer OGRDSin = otb::ogr::DataSource::New(GetParameterString("in.vd"),
+ otb::ogr::DataSource::Modes::Read);
+
+ m_InputGeomSet = InputGeometriesType::New(OGRDSin);
+
+ // Filter instanciation
+ m_GeometriesProjFilter = ProjectionFilterType::New();
+ m_GeometriesProjFilter->SetInput(m_InputGeomSet);
+
+ if (HasValue("in.kwl"))
+ {
+ FloatVectorImageType::Pointer inImage = GetParameterFloatVectorImage("in.kwl");
+ m_GeometriesProjFilter->SetInputOrigin(inImage->GetOrigin()); // nec qd capteur
+ m_GeometriesProjFilter->SetInputSpacing(inImage->GetSpacing()); // nec qd capteur
+ m_GeometriesProjFilter->SetInputKeywordList(inImage->GetImageKeywordlist());
+ //otbAppLogINFO(<<"kwl."<<std::endl);
+ }
+
+ if (GetParameterInt("out.proj") == 0)
+ {
+ FloatVectorImageType::Pointer outImage = GetParameterFloatVectorImage("out.proj.image.in");
+
+ if (outImage)
+ {
+ m_GeometriesProjFilter->SetOutputOrigin(outImage->GetOrigin()); // nec qd capteur
+ m_GeometriesProjFilter->SetOutputSpacing(outImage->GetSpacing()); // nec qd capteur
+ m_OutputProjectionRef = outImage->GetProjectionRef(); // ~ wkt
+ if (m_OutputProjectionRef.empty())
+ {
+ m_GeometriesProjFilter->SetOutputKeywordList(outImage->GetImageKeywordlist()); // nec qd capteur
+ }
+ }
+ else
+ {
+ // LogFATAl
+ }
+ }
+ else
+ {
+
+ // Get the output projection Ref
+ m_OutputProjectionRef = MapProjectionParametersHandler::GetProjectionRefFromChoice(this, "out.proj.user.map");
+
+ }
+ otbAppLogINFO(<<"Geometries reprojection with following projection reference :"
+ <<std::endl<<m_OutputProjectionRef<<std::endl);
+
+ m_GeometriesProjFilter->SetOutputProjectionRef(m_OutputProjectionRef);
+
+ OGRDSout = otb::ogr::DataSource::New(GetParameterString("out.vd"),
+ otb::ogr::DataSource::Modes::Update_LayerOverwrite);
+
+ m_OutputGeomSet = OutputGeometriesType::New(OGRDSout);
+
+ m_GeometriesProjFilter->SetOutput(m_OutputGeomSet);
+ m_GeometriesProjFilter->Update();
+ OGRDSout->SyncToDisk();
+ }
+
+ otb::ogr::DataSource::Pointer OGRDSout;
+ std::string m_OutputProjectionRef;
+ ImageKeywordlist m_KeywordList;
+ InputGeometriesType::Pointer m_InputGeomSet;
+ OutputGeometriesType::Pointer m_OutputGeomSet;
+ //FloatVectorImageType::Pointer m_inImage;
+ // FloatVectorImageType::Pointer m_outImage;
+ ProjectionFilterType::Pointer m_GeometriesProjFilter;
+};
+
+ } // namespace Wrapper
+} // namespace otb
+
+OTB_APPLICATION_EXPORT(otb::Wrapper::VectorDataReprojection)
--
GitLab