include/vrvUtil/DtGreatCircleCurveGenerator.h

Go to the documentation of this file.

/******************************************************************************
** Copyright (c) 2012 MAK Technologies, Inc.
** All rights reserved.
******************************************************************************/


#pragma once

#include <vrvUtil/vrvUtil.h>

#include <matrix/vlVector.h>
#include <vector>

class DtReferenceEllipsoid;

namespace makVrv
{
   
   class DT_DLL_VRVUTIL DtGreatCircleCurveGenerator
   {
   public:
      typedef std::vector<DtVector> GreatCircleCurve;

      DtGreatCircleCurveGenerator( const DtVector& startGeod );

      DtGreatCircleCurveGenerator( const DtGreatCircleCurveGenerator& orig );

      DtGreatCircleCurveGenerator& operator=(const DtGreatCircleCurveGenerator&);

      virtual ~DtGreatCircleCurveGenerator();

      virtual void setRadius( double radius );

      virtual double radius() const;

      virtual const GreatCircleCurve& generateCurve( const DtVector& endGeod,
         double heightScale=0., int maxSegments=500., 
         double preferredSegmentSize=1000. );

      virtual const GreatCircleCurve& curve() const;

      virtual double arcLengthTo( const DtVector& endGeod ) const;

      virtual double initialBearingTo( const DtVector& endGeod ) const;

      virtual DtVector calcPercentAlongGreatCircle( const DtVector& endGeod, 
         double percentage, double heightScale=0. );

      virtual double altitudeForPoint( const DtVector& endGeod, 
         double percent, double heightScale ) const;
      
      virtual double radiusAtLatitude( const DtReferenceEllipsoid& ellipsoid, 
         double lat ) const;

   protected:

      virtual void generateMeridionalCurve( const DtVector& endGeod, 
         int numSegments, double heightScale );

      virtual void generateNonmeridionalCurve( const DtVector& endGeod, 
         int numSegments, double heightScale );

   private:

      DtGreatCircleCurveGenerator();

   protected:

      GreatCircleCurve myCurve;
      DtVector myStartGeod;
      double myRadius;
   };

}