MRG

Code/Projections/otbLambert3CartoSudProjection.h

+/*=========================================================================
+
+  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 __otbLambert3CartoSudProjection_h
+#define __otbLambert3CartoSudProjection_h
+
+
+#include "projection/ossimMapProjection.h"
+#include "projection/ossimLambertConformalConicProjection.h"
+#include "otbMapProjection.h"
+
+namespace otb
+{
+/** \class Lambert3CartoSudProjection
+* \brief This class implements the Lambert3CartoSud map projection (RGF93 french geodesic system).
+* It converts coordinates in longitude,latitude (WGS84) to Lambert 3 map coordinates.
+*
+ */
+template <InverseOrForwardTransformationEnum transform>
+class ITK_EXPORT Lambert3CartoSudProjection : public LambertConformalConicMapProjection<transform>
+{
+public :
+  /** Standard class typedefs. */
+  typedef Lambert3CartoSudProjection                  Self;
+  typedef LambertConformalConicMapProjection<transform>   Superclass;
+  typedef itk::SmartPointer<Self>                          Pointer;
+  typedef itk::SmartPointer<const Self>                    ConstPointer;
+
+  typedef typename Superclass::ScalarType                 ScalarType;
+  typedef itk::Point<ScalarType,2 >                       InputPointType;
+  typedef itk::Point<ScalarType,2 >                        OutputPointType;
+
+  /** Method for creation through the object factory. */
+  itkNewMacro( Self );
+
+  /** Run-time type information (and related methods). */
+  itkTypeMacro( Lambert3CartoSudProjection, LambertConformalConicMapProjection);
+
+
+protected:
+  Lambert3CartoSudProjection();
+  virtual ~Lambert3CartoSudProjection();
+
+private :
+  Lambert3CartoSudProjection(const Self&); //purposely not implemented
+  void operator=(const Self&);                       //purposely not implemented
+};
+
+} // namespace otb
+
+#ifndef OTB_MANUAL_INSTANTIATION
+#include "otbLambert3CartoSudProjection.txx"
+#endif
+
+#endif

Code/Projections/otbLambert3CartoSudProjection.txx

+/*=========================================================================
+
+  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 __otbLambert3CartoSudProjection_txx
+#define __otbLambert3CartoSudProjection_txx
+
+#include "otbLambert3CartoSudProjection.h"
+
+namespace otb
+{
+
+template <InverseOrForwardTransformationEnum transform>
+Lambert3CartoSudProjection<transform>
+::Lambert3CartoSudProjection()
+{
+  itk::Point<double,2> origin;
+  origin[0]=3;
+  origin[1]=46.5;
+  std::string datum = "WE"; //WGS84 datum
+
+
+  double parall1=43.1992913888888888888888888889;
+  double parall2=44.9960938888888888888888888888;
+  double falseEasting=600000;
+  double falseNorthing=3200000;
+  std::string ellipsoid = "CE";
+
+  //TODO: 29-02-2008 Emmanuel: when ossim version > 1.7.2 only
+  // SetOrigin required (remove SetEllipsoid)
+  this->SetOrigin(origin, datum);
+  this->SetEllipsoid(ellipsoid);
+  this->SetParameters(parall1, parall2, falseEasting, falseNorthing);
+}
+
+template <InverseOrForwardTransformationEnum transform>
+Lambert3CartoSudProjection<transform>
+::~Lambert3CartoSudProjection()
+{
+}
+
+
+} // namespace otb
+
+#endif

Code/Projections/otbMapProjections.h

@@ -50,6 +50,7 @@ See OTBCopyright.txt for details.
 #include "otbUtmMapProjection.h"
 #include "otbLambertConformalConicMapProjection.h"
 #include "otbLambert2EtenduProjection.h"
+#include "otbLambert3CartoSudProjection.h"
 #include "otbLambert93Projection.h"
 #include "otbEckert4MapProjection.h"
 #include "otbTransMercatorMapProjection.h"
@@ -84,6 +85,8 @@ typedef LambertConformalConicMapProjection<INVERSE>                LambertConfor
 typedef LambertConformalConicMapProjection<FORWARD>                LambertConformalConicForwardProjection;
 typedef Lambert2EtenduProjection<INVERSE>                          Lambert2EtenduInverseProjection;
 typedef Lambert2EtenduProjection<FORWARD>                          Lambert2EtenduForwardProjection;
+typedef Lambert3CartoSudProjection<INVERSE>                          Lambert3CartoSudInverseProjection;
+typedef Lambert3CartoSudProjection<FORWARD>                          Lambert3CartoSudForwardProjection;
 typedef Lambert93Projection<INVERSE>                               Lambert93InverseProjection;
 typedef Lambert93Projection<FORWARD>                               Lambert93ForwardProjection;
 typedef SVY21MapProjection<INVERSE>                                SVY21InverseProjection;