Julien Osman · db1fe22b · fb4d22f4 · fb4d22f4 · fb4d22f4 · fb4d22f4
--- a/Modules/Core/Transform/include/otbSarSensorModel.h 0 → 100644
+++ b/Modules/Core/Transform/include/otbSarSensorModel.h 0 → 100644
+ * See the License for the specific language governing permissions and
+ * limitations under the License.
+ */
+
+#ifndef otbSarSensorModel_h
+#define otbSarSensorModel_h
+
+#include "otbImageMetadata.h"
+
+
+#include "itkPoint.h"
+
+namespace otb
+{
+
+class SarSensorModel
--- a/Modules/Core/Metadata/include/otbImageMetadataInterfaceBase.h
+++ b/Modules/Core/Metadata/include/otbImageMetadataInterfaceBase.h

  virtual bool CanRead() const {return false;}

-  virtual void Parse(const MetadataSupplierInterface &)
-    {
+  void SupplierCheck() const
+  {
+    if (m_MetadataSupplierInterface == nullptr)
--- a/Modules/Core/Metadata/src/otbCosmoImageMetadataInterface.cxx
+++ b/Modules/Core/Metadata/src/otbCosmoImageMetadataInterface.cxx

 }

-std::vector<Orbit> CosmoImageMetadataInterface::getOrbits(const MetadataSupplierInterface & mds, const std::string & reference_UTC)
+std::vector<Orbit> CosmoImageMetadataInterface::getOrbits(const std::string & reference_UTC)
--- a/Modules/Core/Metadata/src/otbCosmoImageMetadataInterface.cxx
+++ b/Modules/Core/Metadata/src/otbCosmoImageMetadataInterface.cxx
  return orbitVector ;
 }

-void CosmoImageMetadataInterface::ParseGdal(const MetadataSupplierInterface & mds)
+void CosmoImageMetadataInterface::ParseGdal(ImageMetadata & imd)
 {
  // Check acquisition mode and product type
-  Fetch(MDStr::Mode, mds, "Acquisition_Mode");
-  if((m_Imd[MDStr::Mode] != "HIMAGE") &&
-        (m_Imd[MDStr::Mode] != "SPOTLIGHT") &&
-          (m_Imd[MDStr::Mode] != "ENHANCED SPOTLIGHT"))
+  Fetch(MDStr::Mode, imd, "Acquisition_Mode");
+  if((imd[MDStr::Mode] != "HIMAGE") &&
--- a/Modules/Core/Metadata/include/otbImageMetadataInterfaceBase.h
+++ b/Modules/Core/Metadata/include/otbImageMetadataInterfaceBase.h

  virtual bool CanRead() const {return false;}

-  virtual void Parse(const MetadataSupplierInterface &)
-    {
+  void SupplierCheck() const
+  {
+    if (m_MetadataSupplierInterface == nullptr)
+      throw std::runtime_error("In ImageMetadataInterface, the MetadataSupplier needs to be set before calling the Parse function.");
+  }
+
+  virtual void Parse(ImageMetadata &)
--- a/Modules/Core/Metadata/include/otbRadarsat2ImageMetadataInterface.h
+++ b/Modules/Core/Metadata/include/otbRadarsat2ImageMetadataInterface.h
  double GetCenterIncidenceAngle() const override;

  /*get lookup data for calculating backscatter */
-  void CreateCalibrationLookupData(const short type) override;
+  const LookupDataPointerType CreateCalibrationLookupData(const short type) const override;
--- a/Modules/Core/Transform/src/otbSarSensorModel.cxx
+++ b/Modules/Core/Transform/src/otbSarSensorModel.cxx
  line = (timeSinceStart/m_SarParam.azimuthTimeInterval) + currentBurst->startLine;
 }

+void SarSensorModel::SlantRangeToGroundRange(const double & slantRange, const TimeType & azimuthTime, double & groundRange) const
+{
+  ApplyCoordinateConversion(slantRange, azimuthTime, m_SarParam.slantRangeToGroundRangeRecords , groundRange);
+}
+
+void SarSensorModel::ApplyCoordinateConversion(const double & in,
+                                 const TimeType& azimuthTime,
+                                 const std::vector<CoordinateConversionRecord> & records,
+                                 double & out) const
+{
--- a/Modules/Core/Transform/include/otbSarSensorModel.h
+++ b/Modules/Core/Transform/include/otbSarSensorModel.h
                                          Point3DType& sensorPos, 
                                          Vector3DType& sensorVel) const;

+  void LineSampleHeightToWorld(const Point2DType& imPt,
+                               const double & heightAboveEllipsoid,
--- a/Modules/Core/Transform/src/otbSarSensorModel.cxx
+++ b/Modules/Core/Transform/src/otbSarSensorModel.cxx
+  double minDistance = std::numeric_limits<double>::max();
+
+  for (const auto & elem : gcpParam.GCPs)
+  {
+    auto currentDistance = squaredDistance(imPt, elem);
+    if (currentDistance < minDistance)
+    {
+      minDistance = currentDistance;
+      closest = &elem;
+    }
+  }
+
+  return *closest;
+}
+
+void SarSensorModel::projToSurface(const GCP & gcp, const Point2DType& imPt, double heightAboveEllipsoid, Point3DType & ecefPoint) const
--- a/Modules/Core/Transform/src/otbSarSensorModel.cxx
+++ b/Modules/Core/Transform/src/otbSarSensorModel.cxx

 }

+void SarSensorModel::LineSampleHeightToWorld(const Point2DType& imPt,
+                                             const double & heightAboveEllipsoid,
+                                             Point3DType& worldPt) const
+{
+  assert(m_Imd.Has(MDGeom::GCP));
+
+  const auto gcp = findClosestGCP(imPt, m_Imd.GetGCPParam());