AnisotropicInterpolator.h

Go to the documentation of this file.

#pragma once
//------------------------------------------------------------------------------
//------------------------------------------------------------------------------

//----- Included files ---------------------------------------------------------
#include <xmscore/misc/base_macros.h>
#include <xmscore/misc/boost_defines.h>
#include <xmscore/points/ptsfwd.h>
#include <xmscore/stl/vector.h>

//----- Forward declarations ---------------------------------------------------

//----- Namespace declaration --------------------------------------------------
namespace xms
{
//----- Constants / Enumerations -----------------------------------------------

//----- Forward declarations ---------------------------------------------------

//----- Structs / Classes ------------------------------------------------------

struct SNResult
{
  SNResult(size_t a_index, double a_param, double a_cross, const Pt3d& a_position)
  : m_position(a_position)
  , m_param(a_param)
  , m_cross(a_cross)
  , m_index(a_index)
  {
  }

  Pt3d m_position; 
  double m_param;  
  double m_cross;  
  size_t m_index;  
};

struct LineParameters
{
  LineParameters(double a_slope, double a_intercept);
  LineParameters NormalThrough(const Pt3d& a_p) const;
  bool Intersection(const LineParameters& a_other, Pt3d& a_p) const;

  double m_slope;     
  double m_intercept; 
};

typedef std::vector<SNResult> VecSNResult;

class AnisotropicInterpolator
{
public:
  AnisotropicInterpolator();

  void SetPoints(const VecPt3d& a_centerlinePts,
                 const VecPt3d& a_interpolationPts,
                 bool a_pickClosest = false);
  const VecPt3d& GetInterpolationPoints() const;

  double InterpolatePoint(const Pt3d& a_target,
                          const VecPt3d& a_interpolationPoints,
                          double a_xScale,
                          double a_IDWPower = 2,
                          double a_tol = 1.0e-5);

  // Implementation details - for testing only

  VecPt3d ComputeInterpolationPoints(const VecPt3d& a_interpolationPts, bool a_pickClosest);

  double Distance(const Pt3d& a_p0, const Pt3d& a_p1) const;

  double AnisotropicDistanceSquared(const Pt3d& a_p0, const Pt3d& a_p1, double a_xScale);

  double CalculateIDW(double a_xScale,
                      const VecPt3d& a_points,
                      const Pt3d& a_target,
                      double a_power,
                      double a_tol = 1.0e-5);

  double CrossProduct(size_t a_segmentIndex, const Pt3d& a_p) const;

  LineParameters GetLineParameters(const Pt3d& a_p0, const Pt3d& a_p1) const;

  double GetParameterInSegment(size_t a_segmentIndex, const Pt3d& a_p) const;

  bool GetIntersectionOfSegmentWithPoint(size_t a_segmentIndex,
                                         const Pt3d& a_p,
                                         Pt3d& a_intersection,
                                         double& a_param) const;

  void GetNormalPointParameter(size_t a_segmentIndex,
                               const Pt3d& a_p,
                               double& a_lastParam,
                               VecSNResult& a_results) const;

  void GetAllNormalPointParameters(const Pt3d& a_p, VecSNResult& a_results);

  VecPt3d TransformPoint(const Pt3d& a_p, bool a_onlyClosest = false);

private:
  XM_DISALLOW_COPY_AND_ASSIGN(AnisotropicInterpolator);

  // data members:
  VecPt3d m_centerLinePts; 
  VecDbl m_accLengths;     
  std::vector<LineParameters> m_lineParams; 
  VecPt3d m_SNPoints;                       
};

//----- Function prototypes ----------------------------------------------------

} // namespace xms