From ceebee0da4f9b890c5c00f6f3b7961983743f87f Mon Sep 17 00:00:00 2001
From: Otmane Lahlou <otmane.lahlou@c-s.fr>
Date: Thu, 8 Dec 2011 17:37:33 +0100
Subject: [PATCH] ENH: add handler class for the MapProjection parameters
---
...bWrapperMapProjectionParametersHandler.cxx | 238 ++++++++++++++++++
...otbWrapperMapProjectionParametersHandler.h | 83 ++++++
2 files changed, 321 insertions(+)
create mode 100644 Code/ApplicationEngine/otbWrapperMapProjectionParametersHandler.cxx
create mode 100644 Code/ApplicationEngine/otbWrapperMapProjectionParametersHandler.h
diff --git a/Code/ApplicationEngine/otbWrapperMapProjectionParametersHandler.cxx b/Code/ApplicationEngine/otbWrapperMapProjectionParametersHandler.cxx
new file mode 100644
index 0000000000..755af2a5f2
--- /dev/null
+++ b/Code/ApplicationEngine/otbWrapperMapProjectionParametersHandler.cxx
@@ -0,0 +1,238 @@
+/*=========================================================================
+
+ 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 "otbWrapperMapProjectionParametersHandler.h"
+
+// Needed for methods relative to projections
+#include "otbMapProjections.h"
+#include "otbImageToGenericRSOutputParameters.h"
+
+namespace otb
+{
+namespace Wrapper
+{
+
+void MapProjectionParametersHandler::AddMapProjectionParameters( Application::Pointer app, const std::string & key)
+{
+ app->AddParameter(ParameterType_Choice, key, "Output Map Projection");
+ app->SetParameterDescription(key,"Parameters of the ouptut map projection.");
+
+ // utm
+ std::ostringstream oss;
+ oss << key<<".utm";
+ app->AddChoice(oss.str(), "Universal Trans-Mercator (UTM)");
+ app->SetParameterDescription(oss.str(),"A system of transverse mercator projections dividing the surface of Earth between 80S and 84N latitude.");
+
+ oss << ".zone";
+ std::cout <<"zone key "<< oss.str() << std::endl;
+ app->AddParameter(ParameterType_Int, oss.str(), "Zone number");
+ app->SetParameterDescription(oss.str(),"The zone number ranges from 1 to 60 and allows to define the transverse mercator projection (along with the hemisphere)");
+ app->SetMinimumParameterIntValue(oss.str(), 1);
+
+ oss.str("");
+ oss <<key<<".utm" <<".hem";
+ std::cout <<"hem key "<< oss.str() << std::endl;
+ app->AddParameter(ParameterType_Empty, oss.str(), "Northern Hemisphere");
+ app->SetParameterDescription(oss.str(),"The transverse mercator projections are defined by their zone number as well as the hemisphere. Activate this parameter if your image is in the northern hemisphere.");
+
+ // lambert2
+ oss.str("");
+ oss << key<<".lambert2";
+ app->AddChoice(oss.str(), "Lambert II Etendu");
+ app->SetParameterDescription(oss.str(),"This is a Lambert Conformal Conic projection mainly used in France.");
+
+ // lambert93
+ oss.str("");
+ oss << key<<".lambert93";
+ app->AddChoice(oss.str(), "Lambert93");
+ app->SetParameterDescription(oss.str(), "This is a Lambert 93 projection mainly used in France.");
+
+ // Transmercator
+ oss.str("");
+ oss << key<<".transmercator";
+ app->AddChoice(oss.str(), "Transmercator");
+ app->SetParameterDescription(oss.str(),
+ "This is a Transverse Mercator, a cylindrical projection widely used. The required values are set by default.");
+
+ std::ostringstream subParamKey;
+ subParamKey<<oss.str()<<".falseeasting";
+ app->AddParameter(ParameterType_Float, subParamKey.str(), "False easting");
+ app->SetParameterDescription( subParamKey.str(), " Transmercator false easting value.");
+
+ subParamKey.str("");
+ subParamKey<<oss.str()<<".falsenorthing";
+ app->AddParameter(ParameterType_Float, subParamKey.str(), "False northing");
+ app->SetParameterDescription(subParamKey.str(), " Transmercator false northing value.");
+
+ subParamKey.str("");
+ subParamKey<<oss.str()<<".scale";
+ app->AddParameter(ParameterType_Float, subParamKey.str(), "Scale factor");
+ app->SetParameterDescription(subParamKey.str(), " Transmercator scale factor value.");
+
+ // wgs84
+ oss.str("");
+ oss << key<<".wgs";
+ app->AddChoice(oss.str(), "WGS 84");
+ app->SetParameterDescription(oss.str(),"This is a Geographical projection");
+
+ // epsg code
+ oss.str("");
+ oss<<key<<".epsg";
+ std::cout <<"epsg key "<< oss.str() << std::endl;
+ app->AddChoice(oss.str(),"EPSG Code");
+ app->SetParameterDescription(oss.str(),
+ "This code is a generic way of identifying map projections, and allows to specify a large amount of them. See www.spatialreference.org to find which EPSG code is associated to your projection;");
+
+ oss <<".code";
+ std::cout <<"code key "<< oss.str() << std::endl;
+ app->AddParameter(ParameterType_Int, oss.str(), "EPSG Code");
+ app->SetParameterDescription(oss.str(),"See www.spatialreference.org to find which EPSG code is associated to your projection");
+ app->SetDefaultParameterInt(oss.str(), 32631);
+
+ app->SetParameterString(key, "epsg");
+}
+
+/**
+ * Helper method : Compute the ProjectionRef knowing the map
+ * projection picked up by the user
+ *
+ */
+const std::string MapProjectionParametersHandler::GetProjectionRefFromChoice(Application::Pointer app,
+ const std::string & key)
+{
+ std::ostringstream zoneKey;
+ zoneKey << key<<".utm.zone";
+
+ std::ostringstream hemKey;
+ hemKey << key<<".utm.hem";
+
+ std::ostringstream epsgKey;
+ epsgKey << key <<".epsg.code";
+
+ // Get the user choice
+ switch ( app->GetParameterInt(key) )
+ {
+ case Map_Utm:
+ {
+ typedef UtmInverseProjection UtmProjectionType;
+ UtmProjectionType::Pointer utmProjection = UtmProjectionType::New();
+
+ // Set the zone
+ utmProjection->SetZone(app->GetParameterInt(zoneKey.str()));
+
+ // Set the hem
+ std::string hem = "N";
+ if (!app->IsParameterEnabled(hemKey.str()))
+ hem = "S";
+ utmProjection->SetHemisphere(hem[0]);
+
+ // Get the projection ref
+ return utmProjection->GetWkt();
+ }
+ break;
+ case Map_Lambert2:
+ {
+ typedef Lambert2EtenduForwardProjection Lambert2ProjectionType;
+ Lambert2ProjectionType::Pointer lambert2Projection = Lambert2ProjectionType::New();
+ return lambert2Projection->GetWkt();
+ }
+ break;
+ case Map_Lambert93:
+ {
+ typedef otb::Lambert93InverseProjection Lambert93ProjectionType;
+ Lambert93ProjectionType::Pointer lambert93Projection = Lambert93ProjectionType::New();
+ return lambert93Projection->GetWkt();
+ }
+ break;
+ case Map_Transmercator:
+ {
+ typedef otb::TransMercatorInverseProjection TransMercatorProjectionType;
+ TransMercatorProjectionType::Pointer transMercatorProjection = TransMercatorProjectionType::New();
+
+ // Get the Transmercator parameters
+ std::ostringstream falseeasting;
+ std::ostringstream falsenorthing;
+ std::ostringstream scale;
+
+ falseeasting <<key <<".transmercator.falseeasting";
+ falsenorthing <<key <<".transmercator.falsenorthing";
+ scale <<key <<".transmercator.scale";
+
+ transMercatorProjection->SetParameters(app->GetParameterFloat( falseeasting.str() ),
+ app->GetParameterFloat( falsenorthing.str() ),
+ app->GetParameterFloat( scale.str()) );
+ return transMercatorProjection->GetWkt();
+ }
+ break;
+ case Map_WGS84:
+ {
+ return otb::GeoInformationConversion::ToWKT(4326);
+ }
+ break;
+ case Map_Epsg:
+ {
+ return otb::GeoInformationConversion::ToWKT(app->GetParameterInt(epsgKey.str()));
+ }
+ break;
+ }
+}
+
+/**
+ * Helper method : Compute the UTM paramaters relative to an image Origin
+ * Note: The key of the image must be set to be able to get the image.
+ * The key must be totally if the InputImageParameter belongs
+ * to a ParamaterGroup, ie set io.in
+ */
+void MapProjectionParametersHandler::InitializeUTMParameters(Application::Pointer app,
+ const std::string & imageKey,
+ const std::string & mapKey )
+{
+ // Get the UTM params keys
+ std::ostringstream zoneKey;
+ zoneKey << mapKey<<".utm.zone";
+
+ std::ostringstream hemKey;
+ hemKey << mapKey<<".utm.hem";
+
+ // Compute the zone and the hemisphere if not UserValue defined
+ if(!app->HasUserValue(zoneKey.str())
+ && app->HasValue(imageKey)
+ && !app->HasAutomaticValue(zoneKey.str()))
+ {
+ // Compute the Origin lat/long coordinate
+ typedef otb::ImageToGenericRSOutputParameters<FloatVectorImageType> OutputParametersEstimatorType;
+ OutputParametersEstimatorType::Pointer genericRSEstimator = OutputParametersEstimatorType::New();
+ genericRSEstimator->SetInput(app->GetParameterImage(imageKey));
+ genericRSEstimator->SetOutputProjectionRef(otb::GeoInformationConversion::ToWKT(4326));
+ genericRSEstimator->Compute();
+
+ int zone = otb::Utils::GetZoneFromGeoPoint(genericRSEstimator->GetOutputOrigin()[0],
+ genericRSEstimator->GetOutputOrigin()[1]);
+ // Update the UTM Gui fields
+ app->SetParameterInt(zoneKey.str(), zone);
+ if(genericRSEstimator->GetOutputOrigin()[0] > 0.)
+ {
+ app->EnableParameter(hemKey.str());
+ }
+
+ app->AutomaticValueOn(zoneKey.str());
+ app->AutomaticValueOn(hemKey.str());
+ }
+}
+
+}// End namespace Wrapper
+}// End namespace otb
diff --git a/Code/ApplicationEngine/otbWrapperMapProjectionParametersHandler.h b/Code/ApplicationEngine/otbWrapperMapProjectionParametersHandler.h
new file mode 100644
index 0000000000..4e51a9cddf
--- /dev/null
+++ b/Code/ApplicationEngine/otbWrapperMapProjectionParametersHandler.h
@@ -0,0 +1,83 @@
+/*=========================================================================
+
+ 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.
+
+ =========================================================================*/
+#ifndef __otbWrapperMapProjectionParametersHandler_h
+#define __otbWrapperMapProjectionParametersHandler_h
+
+#include "otbWrapperApplication.h"
+
+namespace otb
+{
+namespace Wrapper
+{
+
+// enum used for the map projection choices
+enum
+{
+ Map_Utm,
+ Map_Lambert2,
+ Map_Lambert93,
+ Map_Transmercator,
+ Map_WGS84,
+ Map_Epsg
+};
+
+
+/** \class MapProjectionParametersHandler
+ * \brief This class represent a helper class.
+ *
+ * This class allow creating a Group with several mapProjection
+ * paramaters
+ */
+class MapProjectionParametersHandler
+{
+public:
+ /**
+ * Add a Group containing several choices for map projections
+ *
+ */
+ static void AddMapProjectionParameters(Application::Pointer app, const std::string & key);
+
+ /**
+ * Helper method : Compute the ProjectionRef knowing the map
+ * projection picked up by the user
+ *
+ */
+ static const std::string GetProjectionRefFromChoice(Application::Pointer app,
+ const std::string & key);
+
+
+ /**
+ * Helper method : Compute the UTM paramaters relative an image
+ * Note: The key of the image must be set to be able to get the image.
+ * The key must be totally if the InputImageParameter belongs
+ * to a ParamaterGroup, ie set io.in
+ */
+ static void InitializeUTMParameters(Application::Pointer app,
+ const std::string & imageKey,
+ const std::string & mapKey );
+
+protected:
+ MapProjectionParametersHandler(); // not implemented
+ virtual ~MapProjectionParametersHandler(); // not implemented
+};
+
+}
+}
+
+
+#endif // __otbWrapperMapProjectionParametersHandler_h_
--
GitLab